]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PDF/spdf/abfkw2.F
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / PDF / spdf / abfkw2.F
CommitLineData
fe4da5cc 1*
2* $Id$
3*
4* $Log$
5* Revision 1.1.1.2 1996/10/30 08:27:27 cernlib
6* Version 7.04
7*
8* Revision 1.1.1.1 1996/04/12 15:28:53 plothow
9* Version 7.01
10*
11*
12#include "pdf/pilot.h"
13 SUBROUTINE ABFKW2(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DGL)
14#include "pdf/expdp.inc"
15 + PARPI(40),CALCPI(8,20,25),CALCPIO(8,20,25),ZEROD,
16 + DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DGL
17 REAL X, Q, UPV, DNV, USEA, DSEA, STR, CHM, GL
18 COMMON/W5051I5/CALCPIO
19 REAL XPDF(7)
20 DATA ZEROD/0.D0/
21C----------------------------------------------------------------------
22 DATA (PARPI(I),I=1,40) /
23 + 0.18099999D+00, 0.60926290D+00, 0.10000000D+01, 0.00000000D+00,
24 + 0.00000000D+00, 0.60926290D+00, 0.10000000D+01, 0.00000000D+00,
25 + 0.20329999D+01, 0.96999998D+01, 0.83280414D+00, 0.00000000D+00,
26 + 0.14280005D+01, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
27 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
28 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
29 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
30 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
31 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
32 + 0.00000000D+00, 0.10000002D+00, 0.20000000D+01, 0.10000000D+09/
33 DATA ((CALCPI(I,K,1 ),I=1,8), K=1,9) /
34 + 0.35012248D+00, 0.35012248D+00, 0.33556090D+00, 0.33556090D+00,
35 + 0.33556090D+00, 0.33556090D+00, 0.83161515D+00,-0.52041704D-17,
36 + 0.37201120D+00, 0.37201120D+00, 0.35644488D+00, 0.35643017D+00,
37 + 0.35643017D+00, 0.35644488D+00, 0.16943822D+01, 0.21627137D-01,
38 + 0.39940177D+00, 0.39940177D+00, 0.38279767D+00, 0.38276892D+00,
39 + 0.38276892D+00, 0.38279767D+00, 0.27090403D+01, 0.48760048D-01,
40 + 0.43261897D+00, 0.43261897D+00, 0.41494360D+00, 0.41490152D+00,
41 + 0.41490152D+00, 0.41494360D+00, 0.38788270D+01, 0.81724886D-01,
42 + 0.47194124D+00, 0.47194124D+00, 0.45316071D+00, 0.45310598D+00,
43 + 0.45310598D+00, 0.45316071D+00, 0.52042871D+01, 0.12080131D+00,
44 + 0.51763555D+00, 0.51763555D+00, 0.49771556D+00, 0.49764887D+00,
45 + 0.49764887D+00, 0.49771556D+00, 0.66837094D+01, 0.16625757D+00,
46 + 0.56994488D+00, 0.56994488D+00, 0.54885080D+00, 0.54877279D+00,
47 + 0.54877279D+00, 0.54885080D+00, 0.83134438D+01, 0.21833799D+00,
48 + 0.62908077D+00, 0.62908077D+00, 0.60677763D+00, 0.60668894D+00,
49 + 0.60668894D+00, 0.60677763D+00, 0.10087784D+02, 0.27725545D+00,
50 + 0.69521725D+00, 0.69521725D+00, 0.67166978D+00, 0.67157102D+00,
51 + 0.67157102D+00, 0.67166978D+00, 0.11999100D+02, 0.34318535D+00/
52 DATA ((CALCPI(I,K,1 ),I=1,8), K=10,18) /
53 + 0.76848370D+00, 0.76848370D+00, 0.74365640D+00, 0.74354814D+00,
54 + 0.74354814D+00, 0.74365640D+00, 0.14037940D+02, 0.41625848D+00,
55 + 0.84895905D+00, 0.84895905D+00, 0.82281626D+00, 0.82269906D+00,
56 + 0.82269906D+00, 0.82281626D+00, 0.16193174D+02, 0.49655522D+00,
57 + 0.93666656D+00, 0.93666656D+00, 0.90917250D+00, 0.90904687D+00,
58 + 0.90904687D+00, 0.90917250D+00, 0.18452166D+02, 0.58410027D+00,
59 + 0.10315695D+01, 0.10315695D+01, 0.10026884D+01, 0.10025548D+01,
60 + 0.10025548D+01, 0.10026884D+01, 0.20800966D+02, 0.67885837D+00,
61 + 0.11335678D+01, 0.11335678D+01, 0.11032639D+01, 0.11031228D+01,
62 + 0.11031228D+01, 0.11032639D+01, 0.23224522D+02, 0.78073106D+00,
63 + 0.12424964D+01, 0.12424964D+01, 0.12107340D+01, 0.12105858D+01,
64 + 0.12105858D+01, 0.12107340D+01, 0.25706912D+02, 0.88955454D+00,
65 + 0.13581235D+01, 0.13581235D+01, 0.13248675D+01, 0.13247125D+01,
66 + 0.13247125D+01, 0.13248675D+01, 0.28231576D+02, 0.10050988D+01,
67 + 0.14801520D+01, 0.14801520D+01, 0.14453671D+01, 0.14452057D+01,
68 + 0.14452057D+01, 0.14453671D+01, 0.30781557D+02, 0.11270681D+01,
69 + 0.16082199D+01, 0.16082199D+01, 0.15718716D+01, 0.15717042D+01,
70 + 0.15717042D+01, 0.15718716D+01, 0.33339740D+02, 0.12551020D+01/
71 DATA ((CALCPI(I,K,1 ),I=1,8), K=19,20) /
72 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
73 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
74 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
75 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
76 DATA ((CALCPI(I,K,2 ),I=1,8), K=1,9) /
77 + 0.35450801D+00, 0.35450801D+00, 0.33231685D+00, 0.33231685D+00,
78 + 0.33231685D+00, 0.33231685D+00, 0.83042667D+00,-0.39898640D-16,
79 + 0.37318654D+00, 0.37318654D+00, 0.34955323D+00, 0.34953492D+00,
80 + 0.34953492D+00, 0.34955323D+00, 0.15666705D+01, 0.18428098D-01,
81 + 0.39580234D+00, 0.39580234D+00, 0.37069078D+00, 0.37065508D+00,
82 + 0.37065508D+00, 0.37069078D+00, 0.23994112D+01, 0.40810506D-01,
83 + 0.42247525D+00, 0.42247525D+00, 0.39584907D+00, 0.39579685D+00,
84 + 0.39579685D+00, 0.39584907D+00, 0.33260917D+01, 0.67268661D-01,
85 + 0.45329787D+00, 0.45329787D+00, 0.42512047D+00, 0.42505257D+00,
86 + 0.42505257D+00, 0.42512047D+00, 0.43429635D+01, 0.97896597D-01,
87 + 0.48834142D+00, 0.48834142D+00, 0.45857604D+00, 0.45849322D+00,
88 + 0.45849322D+00, 0.45857604D+00, 0.54448091D+01, 0.13276701D+00,
89 + 0.52765497D+00, 0.52765497D+00, 0.49626478D+00, 0.49616776D+00,
90 + 0.49616776D+00, 0.49626478D+00, 0.66253462D+01, 0.17193050D+00,
91 + 0.57125877D+00, 0.57125877D+00, 0.53820690D+00, 0.53809633D+00,
92 + 0.53809633D+00, 0.53820690D+00, 0.78771848D+01, 0.21540883D+00,
93 + 0.61914185D+00, 0.61914185D+00, 0.58439154D+00, 0.58426801D+00,
94 + 0.58426801D+00, 0.58439154D+00, 0.91919815D+01, 0.26319261D+00/
95 DATA ((CALCPI(I,K,2 ),I=1,8), K=10,18) /
96 + 0.67125862D+00, 0.67125862D+00, 0.63477328D+00, 0.63463732D+00,
97 + 0.63463732D+00, 0.63477328D+00, 0.10560562D+02, 0.31523779D+00,
98 + 0.72752673D+00, 0.72752673D+00, 0.68927008D+00, 0.68912216D+00,
99 + 0.68912216D+00, 0.68927008D+00, 0.11973067D+02, 0.37146361D+00,
100 + 0.78782554D+00, 0.78782554D+00, 0.74776173D+00, 0.74760225D+00,
101 + 0.74760225D+00, 0.74776173D+00, 0.13419105D+02, 0.43175105D+00,
102 + 0.85199542D+00, 0.85199542D+00, 0.81008913D+00, 0.80991844D+00,
103 + 0.80991844D+00, 0.81008913D+00, 0.14887917D+02, 0.49594204D+00,
104 + 0.91983784D+00, 0.91983784D+00, 0.87605444D+00, 0.87587283D+00,
105 + 0.87587283D+00, 0.87605444D+00, 0.16368538D+02, 0.56383966D+00,
106 + 0.99111626D+00, 0.99111626D+00, 0.94542193D+00, 0.94522965D+00,
107 + 0.94522965D+00, 0.94542193D+00, 0.17849959D+02, 0.63520897D+00,
108 + 0.10655579D+01, 0.10655579D+01, 0.10179197D+01, 0.10177170D+01,
109 + 0.10177170D+01, 0.10179197D+01, 0.19321289D+02, 0.70977869D+00,
110 + 0.11428559D+01, 0.11428559D+01, 0.10932421D+01, 0.10930291D+01,
111 + 0.10930291D+01, 0.10932421D+01, 0.20771906D+02, 0.78724363D+00,
112 + 0.12226727D+01, 0.12226727D+01, 0.11710528D+01, 0.11708296D+01,
113 + 0.11708296D+01, 0.11710528D+01, 0.22191594D+02, 0.86726773D+00/
114 DATA ((CALCPI(I,K,2 ),I=1,8), K=19,20) /
115 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
116 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
117 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
118 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
119 DATA ((CALCPI(I,K,3 ),I=1,8), K=1,9) /
120 + 0.35969733D+00, 0.35969733D+00, 0.32591301D+00, 0.32591301D+00,
121 + 0.32591301D+00, 0.32591301D+00, 0.82805124D+00,-0.43368087D-16,
122 + 0.37500840D+00, 0.37500840D+00, 0.33921282D+00, 0.33918988D+00,
123 + 0.33918988D+00, 0.33921282D+00, 0.14320118D+01, 0.15200902D-01,
124 + 0.39298618D+00, 0.39298618D+00, 0.35514484D+00, 0.35509987D+00,
125 + 0.35509987D+00, 0.35514484D+00, 0.20896087D+01, 0.33103685D-01,
126 + 0.41363247D+00, 0.41363247D+00, 0.37371096D+00, 0.37364473D+00,
127 + 0.37364473D+00, 0.37371096D+00, 0.27961368D+01, 0.53711536D-01,
128 + 0.43693480D+00, 0.43693480D+00, 0.39489890D+00, 0.39481209D+00,
129 + 0.39481209D+00, 0.39489890D+00, 0.35466994D+01, 0.77013287D-01,
130 + 0.46285814D+00, 0.46285814D+00, 0.41867397D+00, 0.41856714D+00,
131 + 0.41856714D+00, 0.41867397D+00, 0.43355598D+01, 0.10297524D+00,
132 + 0.49135036D+00, 0.49135036D+00, 0.44498453D+00, 0.44485811D+00,
133 + 0.44485811D+00, 0.44498453D+00, 0.51565328D+01, 0.13154665D+00,
134 + 0.52233727D+00, 0.52233727D+00, 0.47375703D+00, 0.47361138D+00,
135 + 0.47361138D+00, 0.47375703D+00, 0.60029852D+01, 0.16265474D+00,
136 + 0.55572286D+00, 0.55572286D+00, 0.50489627D+00, 0.50473167D+00,
137 + 0.50473167D+00, 0.50489627D+00, 0.68679706D+01, 0.19620494D+00/
138 DATA ((CALCPI(I,K,3 ),I=1,8), K=10,18) /
139 + 0.59138844D+00, 0.59138844D+00, 0.53828453D+00, 0.53810115D+00,
140 + 0.53810115D+00, 0.53828453D+00, 0.77443392D+01, 0.23207998D+00,
141 + 0.62919277D+00, 0.62919277D+00, 0.57378174D+00, 0.57357971D+00,
142 + 0.57357971D+00, 0.57378174D+00, 0.86248540D+01, 0.27014009D+00,
143 + 0.66897271D+00, 0.66897271D+00, 0.61122609D+00, 0.61100547D+00,
144 + 0.61100547D+00, 0.61122609D+00, 0.95023062D+01, 0.31022354D+00,
145 + 0.71054429D+00, 0.71054429D+00, 0.65043515D+00, 0.65019595D+00,
146 + 0.65019595D+00, 0.65043515D+00, 0.10369625D+02, 0.35214781D+00,
147 + 0.75370431D+00, 0.75370431D+00, 0.69120741D+00, 0.69094963D+00,
148 + 0.69094963D+00, 0.69120741D+00, 0.11219984D+02, 0.39571110D+00,
149 + 0.79823235D+00, 0.79823235D+00, 0.73332435D+00, 0.73304796D+00,
150 + 0.73304796D+00, 0.73332435D+00, 0.12046891D+02, 0.44069438D+00,
151 + 0.84389311D+00, 0.84389311D+00, 0.77655277D+00, 0.77625773D+00,
152 + 0.77625773D+00, 0.77655277D+00, 0.12844281D+02, 0.48686372D+00,
153 + 0.89043914D+00, 0.89043914D+00, 0.82064746D+00, 0.82033373D+00,
154 + 0.82033373D+00, 0.82064746D+00, 0.13606581D+02, 0.53397296D+00,
155 + 0.93761367D+00, 0.93761367D+00, 0.86535407D+00, 0.86502166D+00,
156 + 0.86502166D+00, 0.86535407D+00, 0.14328774D+02, 0.58176661D+00/
157 DATA ((CALCPI(I,K,3 ),I=1,8), K=19,20) /
158 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
159 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
160 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
161 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
162 DATA ((CALCPI(I,K,4 ),I=1,8), K=1,9) /
163 + 0.36476665D+00, 0.36476665D+00, 0.31343633D+00, 0.31343633D+00,
164 + 0.31343633D+00, 0.31343633D+00, 0.82330650D+00,-0.38163916D-16,
165 + 0.37664555D+00, 0.37664555D+00, 0.32260136D+00, 0.32256942D+00,
166 + 0.32256942D+00, 0.32260136D+00, 0.12945379D+01, 0.12128303D-01,
167 + 0.39022327D+00, 0.39022327D+00, 0.33344256D+00, 0.33337871D+00,
168 + 0.33337871D+00, 0.33344256D+00, 0.17895988D+01, 0.26008755D-01,
169 + 0.40543373D+00, 0.40543373D+00, 0.34589471D+00, 0.34579887D+00,
170 + 0.34579887D+00, 0.34589471D+00, 0.23035577D+01, 0.41576230D-01,
171 + 0.42220829D+00, 0.42220829D+00, 0.35989022D+00, 0.35976223D+00,
172 + 0.35976223D+00, 0.35989022D+00, 0.28318563D+01, 0.58763054D-01,
173 + 0.44046046D+00, 0.44046046D+00, 0.37534387D+00, 0.37518352D+00,
174 + 0.37518352D+00, 0.37534387D+00, 0.33696021D+01, 0.77483727D-01,
175 + 0.46009403D+00, 0.46009403D+00, 0.39216094D+00, 0.39196800D+00,
176 + 0.39196800D+00, 0.39216094D+00, 0.39118767D+01, 0.97643063D-01,
177 + 0.48099992D+00, 0.48099992D+00, 0.41023409D+00, 0.41000832D+00,
178 + 0.41000832D+00, 0.41023409D+00, 0.44537567D+01, 0.11913303D+00,
179 + 0.50305752D+00, 0.50305752D+00, 0.42944471D+00, 0.42918588D+00,
180 + 0.42918588D+00, 0.42944471D+00, 0.49904042D+01, 0.14183407D+00/
181 DATA ((CALCPI(I,K,4 ),I=1,8), K=10,18) /
182 + 0.52613523D+00, 0.52613523D+00, 0.44966342D+00, 0.44937137D+00,
183 + 0.44937137D+00, 0.44966342D+00, 0.55171440D+01, 0.16561567D+00,
184 + 0.55009155D+00, 0.55009155D+00, 0.47075121D+00, 0.47042580D+00,
185 + 0.47042580D+00, 0.47075121D+00, 0.60295332D+01, 0.19033738D+00,
186 + 0.57477639D+00, 0.57477639D+00, 0.49256067D+00, 0.49220187D+00,
187 + 0.49220187D+00, 0.49256067D+00, 0.65234244D+01, 0.21585013D+00,
188 + 0.60003260D+00, 0.60003260D+00, 0.51493758D+00, 0.51454544D+00,
189 + 0.51454544D+00, 0.51493758D+00, 0.69950194D+01, 0.24199778D+00,
190 + 0.62569766D+00, 0.62569766D+00, 0.53772257D+00, 0.53729726D+00,
191 + 0.53729726D+00, 0.53772257D+00, 0.74409132D+01, 0.26861880D+00,
192 + 0.65160554D+00, 0.65160554D+00, 0.56075296D+00, 0.56029479D+00,
193 + 0.56029479D+00, 0.56075296D+00, 0.78581281D+01, 0.29554809D+00,
194 + 0.67758862D+00, 0.67758862D+00, 0.58386467D+00, 0.58337408D+00,
195 + 0.58337408D+00, 0.58386467D+00, 0.82441367D+01, 0.32261891D+00,
196 + 0.70347962D+00, 0.70347962D+00, 0.60689412D+00, 0.60637171D+00,
197 + 0.60637171D+00, 0.60689412D+00, 0.85968749D+01, 0.34966480D+00,
198 + 0.72911352D+00, 0.72911352D+00, 0.62968019D+00, 0.62912673D+00,
199 + 0.62912673D+00, 0.62968019D+00, 0.89147454D+01, 0.37652154D+00/
200 DATA ((CALCPI(I,K,4 ),I=1,8), K=19,20) /
201 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
202 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
203 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
204 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
205 DATA ((CALCPI(I,K,5 ),I=1,8), K=1,9) /
206 + 0.36743196D+00, 0.36743196D+00, 0.28976010D+00, 0.28976010D+00,
207 + 0.28976010D+00, 0.28976010D+00, 0.81384158D+00,-0.38163916D-16,
208 + 0.37601016D+00, 0.37601016D+00, 0.29484821D+00, 0.29478982D+00,
209 + 0.29478982D+00, 0.29484821D+00, 0.11623402D+01, 0.94648758D-02,
210 + 0.38554691D+00, 0.38554691D+00, 0.30090458D+00, 0.30078741D+00,
211 + 0.30078741D+00, 0.30090458D+00, 0.15148449D+01, 0.19952699D-01,
212 + 0.39594768D+00, 0.39594768D+00, 0.30783659D+00, 0.30766050D+00,
213 + 0.30766050D+00, 0.30783659D+00, 0.18673407D+01, 0.31369239D-01,
214 + 0.40712527D+00, 0.40712527D+00, 0.31555913D+00, 0.31532421D+00,
215 + 0.31532421D+00, 0.31555913D+00, 0.22164485D+01, 0.43627684D-01,
216 + 0.41898255D+00, 0.41898255D+00, 0.32397756D+00, 0.32368417D+00,
217 + 0.32368417D+00, 0.32397756D+00, 0.25588348D+01, 0.56631374D-01,
218 + 0.43141994D+00, 0.43141994D+00, 0.33299510D+00, 0.33264383D+00,
219 + 0.33264383D+00, 0.33299510D+00, 0.28913769D+01, 0.70281217D-01,
220 + 0.44433366D+00, 0.44433366D+00, 0.34251111D+00, 0.34210286D+00,
221 + 0.34210286D+00, 0.34251111D+00, 0.32111869D+01, 0.84473997D-01,
222 + 0.45761704D+00, 0.45761704D+00, 0.35242237D+00, 0.35195829D+00,
223 + 0.35195829D+00, 0.35242237D+00, 0.35156443D+01, 0.99103611D-01/
224 DATA ((CALCPI(I,K,5 ),I=1,8), K=10,18) /
225 + 0.47116134D+00, 0.47116134D+00, 0.36262385D+00, 0.36210540D+00,
226 + 0.36210540D+00, 0.36262385D+00, 0.38024294D+01, 0.11406192D+00,
227 + 0.48485678D+00, 0.48485678D+00, 0.37300979D+00, 0.37243870D+00,
228 + 0.37243870D+00, 0.37300979D+00, 0.40695430D+01, 0.12923975D+00,
229 + 0.49859366D+00, 0.49859366D+00, 0.38347470D+00, 0.38285301D+00,
230 + 0.38285301D+00, 0.38347470D+00, 0.43153217D+01, 0.14452797D+00,
231 + 0.51226343D+00, 0.51226343D+00, 0.39391450D+00, 0.39324452D+00,
232 + 0.39324452D+00, 0.39391450D+00, 0.45384440D+01, 0.15981863D+00,
233 + 0.52575992D+00, 0.52575992D+00, 0.40422763D+00, 0.40351196D+00,
234 + 0.40351196D+00, 0.40422763D+00, 0.47379298D+01, 0.17500612D+00,
235 + 0.53898036D+00, 0.53898036D+00, 0.41431616D+00, 0.41355766D+00,
236 + 0.41355766D+00, 0.41431616D+00, 0.49131333D+01, 0.18998820D+00,
237 + 0.55182652D+00, 0.55182652D+00, 0.42408677D+00, 0.42328855D+00,
238 + 0.42328855D+00, 0.42408677D+00, 0.50637299D+01, 0.20466714D+00,
239 + 0.56420565D+00, 0.56420565D+00, 0.43345175D+00, 0.43261714D+00,
240 + 0.43261714D+00, 0.43345175D+00, 0.51896975D+01, 0.21895064D+00,
241 + 0.57603135D+00, 0.57603135D+00, 0.44232978D+00, 0.44146236D+00,
242 + 0.44146236D+00, 0.44232978D+00, 0.52912945D+01, 0.23275273D+00/
243 DATA ((CALCPI(I,K,5 ),I=1,8), K=19,20) /
244 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
245 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
246 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
247 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
248 DATA ((CALCPI(I,K,6 ),I=1,8), K=1,9) /
249 + 0.36371961D+00, 0.36371961D+00, 0.24715937D+00, 0.24715937D+00,
250 + 0.24715937D+00, 0.24715937D+00, 0.79501053D+00,-0.41633363D-16,
251 + 0.36918406D+00, 0.36918406D+00, 0.24848601D+00, 0.24838839D+00,
252 + 0.24838839D+00, 0.24848601D+00, 0.10329147D+01, 0.71774690D-02,
253 + 0.37506236D+00, 0.37506236D+00, 0.25030064D+00, 0.25011013D+00,
254 + 0.25011013D+00, 0.25030064D+00, 0.12630821D+01, 0.14825257D-01,
255 + 0.38126493D+00, 0.38126493D+00, 0.25251687D+00, 0.25223842D+00,
256 + 0.25223842D+00, 0.25251687D+00, 0.14831325D+01, 0.22854718D-01,
257 + 0.38771416D+00, 0.38771416D+00, 0.25506046D+00, 0.25469924D+00,
258 + 0.25469924D+00, 0.25506046D+00, 0.16912725D+01, 0.31189161D-01,
259 + 0.39432999D+00, 0.39432999D+00, 0.25785517D+00, 0.25741654D+00,
260 + 0.25741654D+00, 0.25785517D+00, 0.18859428D+01, 0.39749592D-01,
261 + 0.40103450D+00, 0.40103450D+00, 0.26082726D+00, 0.26031679D+00,
262 + 0.26031679D+00, 0.26082726D+00, 0.20658360D+01, 0.48459182D-01,
263 + 0.40775133D+00, 0.40775133D+00, 0.26390484D+00, 0.26332830D+00,
264 + 0.26332830D+00, 0.26390484D+00, 0.22299159D+01, 0.57242734D-01,
265 + 0.41440630D+00, 0.41440630D+00, 0.26701852D+00, 0.26638184D+00,
266 + 0.26638184D+00, 0.26701852D+00, 0.23774071D+01, 0.66027293D-01/
267 DATA ((CALCPI(I,K,6 ),I=1,8), K=10,18) /
268 + 0.42092808D+00, 0.42092808D+00, 0.27010197D+00, 0.26941122D+00,
269 + 0.26941122D+00, 0.27010197D+00, 0.25077938D+01, 0.74742768D-01,
270 + 0.42724871D+00, 0.42724871D+00, 0.27309241D+00, 0.27235378D+00,
271 + 0.27235378D+00, 0.27309241D+00, 0.26208077D+01, 0.83322457D-01,
272 + 0.43330415D+00, 0.43330415D+00, 0.27593111D+00, 0.27515086D+00,
273 + 0.27515086D+00, 0.27593111D+00, 0.27164144D+01, 0.91703603D-01,
274 + 0.43903486D+00, 0.43903486D+00, 0.27856389D+00, 0.27774834D+00,
275 + 0.27774834D+00, 0.27856389D+00, 0.27947967D+01, 0.99827898D-01,
276 + 0.44438619D+00, 0.44438619D+00, 0.28094152D+00, 0.28009698D+00,
277 + 0.28009698D+00, 0.28094152D+00, 0.28563334D+01, 0.10764193D+00,
278 + 0.44930884D+00, 0.44930884D+00, 0.28302004D+00, 0.28215280D+00,
279 + 0.28215280D+00, 0.28302004D+00, 0.29015781D+01, 0.11509756D+00,
280 + 0.45375911D+00, 0.45375911D+00, 0.28476105D+00, 0.28387731D+00,
281 + 0.28387731D+00, 0.28476105D+00, 0.29312349D+01, 0.12215221D+00,
282 + 0.45769919D+00, 0.45769919D+00, 0.28613187D+00, 0.28523773D+00,
283 + 0.28523773D+00, 0.28613187D+00, 0.29461347D+01, 0.12876909D+00,
284 + 0.46109722D+00, 0.46109722D+00, 0.28710561D+00, 0.28620701D+00,
285 + 0.28620701D+00, 0.28710561D+00, 0.29472105D+01, 0.13491727D+00/
286 DATA ((CALCPI(I,K,6 ),I=1,8), K=19,20) /
287 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
288 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
289 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
290 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
291 DATA ((CALCPI(I,K,7 ),I=1,8), K=1,9) /
292 + 0.35031868D+00, 0.35031868D+00, 0.17838643D+00, 0.17838643D+00,
293 + 0.17838643D+00, 0.17838643D+00, 0.75774795D+00,-0.17347235D-17,
294 + 0.35319490D+00, 0.35319490D+00, 0.17707442D+00, 0.17700394D+00,
295 + 0.17700394D+00, 0.17707442D+00, 0.89825513D+00, 0.51655107D-02,
296 + 0.35610678D+00, 0.35610678D+00, 0.17599078D+00, 0.17585755D+00,
297 + 0.17585755D+00, 0.17599078D+00, 0.10264673D+01, 0.10419464D-01,
298 + 0.35898919D+00, 0.35898919D+00, 0.17507296D+00, 0.17488466D+00,
299 + 0.17488466D+00, 0.17507296D+00, 0.11417077D+01, 0.15700814D-01,
300 + 0.36178843D+00, 0.36178843D+00, 0.17427026D+00, 0.17403450D+00,
301 + 0.17403450D+00, 0.17427026D+00, 0.12436244D+01, 0.20959492D-01,
302 + 0.36445318D+00, 0.36445318D+00, 0.17353467D+00, 0.17325888D+00,
303 + 0.17325888D+00, 0.17353467D+00, 0.13320989D+01, 0.26147937D-01,
304 + 0.36693594D+00, 0.36693594D+00, 0.17282226D+00, 0.17251366D+00,
305 + 0.17251366D+00, 0.17282226D+00, 0.14071820D+01, 0.31222316D-01,
306 + 0.36919317D+00, 0.36919317D+00, 0.17209324D+00, 0.17175880D+00,
307 + 0.17175880D+00, 0.17209324D+00, 0.14691039D+01, 0.36142724D-01,
308 + 0.37118536D+00, 0.37118536D+00, 0.17131184D+00, 0.17095825D+00,
309 + 0.17095825D+00, 0.17131184D+00, 0.15182466D+01, 0.40873144D-01/
310 DATA ((CALCPI(I,K,7 ),I=1,8), K=10,18) /
311 + 0.37287724D+00, 0.37287724D+00, 0.17044653D+00, 0.17008011D+00,
312 + 0.17008011D+00, 0.17044653D+00, 0.15551267D+01, 0.45381665D-01,
313 + 0.37423788D+00, 0.37423788D+00, 0.16946996D+00, 0.16909670D+00,
314 + 0.16909670D+00, 0.16946996D+00, 0.15803743D+01, 0.49640539D-01,
315 + 0.37524080D+00, 0.37524080D+00, 0.16835903D+00, 0.16798451D+00,
316 + 0.16798451D+00, 0.16835903D+00, 0.15947127D+01, 0.53626232D-01,
317 + 0.37586403D+00, 0.37586403D+00, 0.16709480D+00, 0.16672417D+00,
318 + 0.16672417D+00, 0.16709480D+00, 0.15989375D+01, 0.57319427D-01,
319 + 0.37609000D+00, 0.37609000D+00, 0.16566240D+00, 0.16530040D+00,
320 + 0.16530040D+00, 0.16566240D+00, 0.15938966D+01, 0.60704969D-01,
321 + 0.37590560D+00, 0.37590560D+00, 0.16405095D+00, 0.16370185D+00,
322 + 0.16370185D+00, 0.16405095D+00, 0.15804711D+01, 0.63771761D-01,
323 + 0.37530195D+00, 0.37530195D+00, 0.16225331D+00, 0.16192094D+00,
324 + 0.16192094D+00, 0.16225331D+00, 0.15595582D+01, 0.66512618D-01,
325 + 0.37427429D+00, 0.37427429D+00, 0.16026589D+00, 0.15995363D+00,
326 + 0.15995363D+00, 0.16026589D+00, 0.15320546D+01, 0.68924063D-01,
327 + 0.37282174D+00, 0.37282174D+00, 0.15808841D+00, 0.15779918D+00,
328 + 0.15779918D+00, 0.15808841D+00, 0.14988422D+01, 0.71006089D-01/
329 DATA ((CALCPI(I,K,7 ),I=1,8), K=19,20) /
330 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
331 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
332 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
333 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
334 DATA ((CALCPI(I,K,8 ),I=1,8), K=1,9) /
335 + 0.33891365D+00, 0.33891365D+00, 0.12193780D+00, 0.12193780D+00,
336 + 0.12193780D+00, 0.12193780D+00, 0.71647517D+00,-0.22551405D-16,
337 + 0.34048312D+00, 0.34048312D+00, 0.12014300D+00, 0.12010115D+00,
338 + 0.12010115D+00, 0.12014300D+00, 0.80298019D+00, 0.39753828D-02,
339 + 0.34189642D+00, 0.34189642D+00, 0.11848590D+00, 0.11841177D+00,
340 + 0.11841177D+00, 0.11848590D+00, 0.87734632D+00, 0.78824510D-02,
341 + 0.34310846D+00, 0.34310846D+00, 0.11692131D+00, 0.11682346D+00,
342 + 0.11682346D+00, 0.11692131D+00, 0.93976006D+00, 0.11682687D-01,
343 + 0.34408357D+00, 0.34408357D+00, 0.11541412D+00, 0.11530029D+00,
344 + 0.11530029D+00, 0.11541412D+00, 0.99054686D+00, 0.15345970D-01,
345 + 0.34479032D+00, 0.34479032D+00, 0.11393328D+00, 0.11381038D+00,
346 + 0.11381038D+00, 0.11393328D+00, 0.10302042D+01, 0.18846087D-01,
347 + 0.34520108D+00, 0.34520108D+00, 0.11245157D+00, 0.11232570D+00,
348 + 0.11232570D+00, 0.11245157D+00, 0.10593167D+01, 0.22160222D-01,
349 + 0.34529249D+00, 0.34529249D+00, 0.11094584D+00, 0.11082237D+00,
350 + 0.11082237D+00, 0.11094584D+00, 0.10785614D+01, 0.25269373D-01,
351 + 0.34504522D+00, 0.34504522D+00, 0.10939667D+00, 0.10928023D+00,
352 + 0.10928023D+00, 0.10939667D+00, 0.10886816D+01, 0.28158041D-01/
353 DATA ((CALCPI(I,K,8 ),I=1,8), K=10,18) /
354 + 0.34444397D+00, 0.34444397D+00, 0.10778826D+00, 0.10768282D+00,
355 + 0.10768282D+00, 0.10778826D+00, 0.10904685D+01, 0.30814209D-01,
356 + 0.34347740D+00, 0.34347740D+00, 0.10610827D+00, 0.10601714D+00,
357 + 0.10601714D+00, 0.10610827D+00, 0.10847430D+01, 0.33229207D-01,
358 + 0.34213798D+00, 0.34213798D+00, 0.10434759D+00, 0.10427349D+00,
359 + 0.10427349D+00, 0.10434759D+00, 0.10723392D+01, 0.35397545D-01,
360 + 0.34042185D+00, 0.34042185D+00, 0.10250016D+00, 0.10244523D+00,
361 + 0.10244523D+00, 0.10250016D+00, 0.10540880D+01, 0.37316742D-01,
362 + 0.33832865D+00, 0.33832865D+00, 0.10056269D+00, 0.10052859D+00,
363 + 0.10052859D+00, 0.10056269D+00, 0.10308048D+01, 0.38987108D-01,
364 + 0.33586130D+00, 0.33586130D+00, 0.98534461D-01, 0.98522331D-01,
365 + 0.98522331D-01, 0.98534461D-01, 0.10032768D+01, 0.40411516D-01,
366 + 0.33302574D+00, 0.33302574D+00, 0.96416984D-01, 0.96427558D-01,
367 + 0.96427558D-01, 0.96416984D-01, 0.97225332D+00, 0.41595142D-01,
368 + 0.32983067D+00, 0.32983067D+00, 0.94213779D-01, 0.94247393D-01,
369 + 0.94247393D-01, 0.94213779D-01, 0.93843835D+00, 0.42545203D-01,
370 + 0.32628730D+00, 0.32628730D+00, 0.91930065D-01, 0.91986707D-01,
371 + 0.91986707D-01, 0.91930065D-01, 0.90248427D+00, 0.43270681D-01/
372 DATA ((CALCPI(I,K,8 ),I=1,8), K=19,20) /
373 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
374 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
375 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
376 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
377 DATA ((CALCPI(I,K,9 ),I=1,8), K=1,9) /
378 + 0.33238608D+00, 0.33238608D+00, 0.70040595D-01, 0.70040595D-01,
379 + 0.70040595D-01, 0.70040595D-01, 0.66031793D+00, 0.17347235D-17,
380 + 0.33271345D+00, 0.33271345D+00, 0.68816053D-01, 0.68818244D-01,
381 + 0.68818244D-01, 0.68816053D-01, 0.70473315D+00, 0.29681705D-02,
382 + 0.33270733D+00, 0.33270733D+00, 0.67648993D-01, 0.67656660D-01,
383 + 0.67656660D-01, 0.67648993D-01, 0.73875315D+00, 0.57814180D-02,
384 + 0.33234675D+00, 0.33234675D+00, 0.66510921D-01, 0.66526622D-01,
385 + 0.66526622D-01, 0.66510921D-01, 0.76312692D+00, 0.84221450D-02,
386 + 0.33161675D+00, 0.33161675D+00, 0.65380695D-01, 0.65406431D-01,
387 + 0.65406431D-01, 0.65380695D-01, 0.77857886D+00, 0.10877773D-01,
388 + 0.33050725D+00, 0.33050725D+00, 0.64241407D-01, 0.64278674D-01,
389 + 0.64278674D-01, 0.64241407D-01, 0.78592436D+00, 0.13139682D-01,
390 + 0.32901128D+00, 0.32901128D+00, 0.63079069D-01, 0.63128910D-01,
391 + 0.63128910D-01, 0.63079069D-01, 0.78598987D+00, 0.15201618D-01,
392 + 0.32712559D+00, 0.32712559D+00, 0.61882946D-01, 0.61946014D-01,
393 + 0.61946014D-01, 0.61882946D-01, 0.77961309D+00, 0.17060166D-01,
394 + 0.32485031D+00, 0.32485031D+00, 0.60645114D-01, 0.60721713D-01,
395 + 0.60721713D-01, 0.60645114D-01, 0.76762637D+00, 0.18714363D-01/
396 DATA ((CALCPI(I,K,9 ),I=1,8), K=10,18) /
397 + 0.32218876D+00, 0.32218876D+00, 0.59360240D-01, 0.59450376D-01,
398 + 0.59450376D-01, 0.59360240D-01, 0.75084187D+00, 0.20165536D-01,
399 + 0.31914731D+00, 0.31914731D+00, 0.58025321D-01, 0.58128746D-01,
400 + 0.58128746D-01, 0.58025321D-01, 0.73004005D+00, 0.21417078D-01,
401 + 0.31573511D+00, 0.31573511D+00, 0.56639423D-01, 0.56755673D-01,
402 + 0.56755673D-01, 0.56639423D-01, 0.70596015D+00, 0.22474230D-01,
403 + 0.31196386D+00, 0.31196386D+00, 0.55203419D-01, 0.55331851D-01,
404 + 0.55331851D-01, 0.55203419D-01, 0.67929239D+00, 0.23343841D-01,
405 + 0.30784758D+00, 0.30784758D+00, 0.53719730D-01, 0.53859557D-01,
406 + 0.53859557D-01, 0.53719730D-01, 0.65067186D+00, 0.24034136D-01,
407 + 0.30340232D+00, 0.30340232D+00, 0.52192070D-01, 0.52342393D-01,
408 + 0.52342393D-01, 0.52192070D-01, 0.62067497D+00, 0.24554476D-01,
409 + 0.29864590D+00, 0.29864590D+00, 0.50625206D-01, 0.50785040D-01,
410 + 0.50785040D-01, 0.50625206D-01, 0.58981742D+00, 0.24915126D-01,
411 + 0.29359760D+00, 0.29359760D+00, 0.49024713D-01, 0.49193014D-01,
412 + 0.49193014D-01, 0.49024713D-01, 0.55855334D+00, 0.25127030D-01,
413 + 0.28827795D+00, 0.28827795D+00, 0.47396764D-01, 0.47572451D-01,
414 + 0.47572451D-01, 0.47396764D-01, 0.52727628D+00, 0.25201593D-01/
415 DATA ((CALCPI(I,K,9 ),I=1,8), K=19,20) /
416 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
417 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
418 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
419 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
420 DATA ((CALCPI(I,K,10),I=1,8), K=1,9) /
421 + 0.33311612D+00, 0.33311612D+00, 0.38900848D-01, 0.38900848D-01,
422 + 0.38900848D-01, 0.38900848D-01, 0.60555744D+00, 0.00000000D+00,
423 + 0.33199868D+00, 0.33199868D+00, 0.38440375D-01, 0.38454832D-01,
424 + 0.38454832D-01, 0.38440375D-01, 0.62458520D+00, 0.22866471D-02,
425 + 0.33043676D+00, 0.33043676D+00, 0.37969307D-01, 0.37997962D-01,
426 + 0.37997962D-01, 0.37969307D-01, 0.63520631D+00, 0.43898213D-02,
427 + 0.32843120D+00, 0.32843120D+00, 0.37472250D-01, 0.37514590D-01,
428 + 0.37514590D-01, 0.37472250D-01, 0.63841151D+00, 0.63055117D-02,
429 + 0.32598518D+00, 0.32598518D+00, 0.36938515D-01, 0.36993857D-01,
430 + 0.36993857D-01, 0.36938515D-01, 0.63506941D+00, 0.80320343D-02,
431 + 0.32310653D+00, 0.32310653D+00, 0.36361096D-01, 0.36428622D-01,
432 + 0.36428622D-01, 0.36361096D-01, 0.62607358D+00, 0.95709924D-02,
433 + 0.31980523D+00, 0.31980523D+00, 0.35734938D-01, 0.35813729D-01,
434 + 0.35813729D-01, 0.35734938D-01, 0.61228013D+00, 0.10925353D-01,
435 + 0.31609407D+00, 0.31609407D+00, 0.35057287D-01, 0.35146350D-01,
436 + 0.35146350D-01, 0.35057287D-01, 0.59450584D+00, 0.12099894D-01,
437 + 0.31198829D+00, 0.31198829D+00, 0.34327272D-01, 0.34425566D-01,
438 + 0.34425566D-01, 0.34327272D-01, 0.57351938D+00, 0.13100839D-01/
439 DATA ((CALCPI(I,K,10),I=1,8), K=10,18) /
440 + 0.30750533D+00, 0.30750533D+00, 0.33545668D-01, 0.33652126D-01,
441 + 0.33652126D-01, 0.33545668D-01, 0.55003264D+00, 0.13935675D-01,
442 + 0.30266466D+00, 0.30266466D+00, 0.32714650D-01, 0.32828198D-01,
443 + 0.32828198D-01, 0.32714650D-01, 0.52469578D+00, 0.14612938D-01,
444 + 0.29748745D+00, 0.29748745D+00, 0.31837544D-01, 0.31957117D-01,
445 + 0.31957117D-01, 0.31837544D-01, 0.49809326D+00, 0.15142004D-01,
446 + 0.29199638D+00, 0.29199638D+00, 0.30918602D-01, 0.31043158D-01,
447 + 0.31043158D-01, 0.30918602D-01, 0.47074269D+00, 0.15532884D-01,
448 + 0.28621535D+00, 0.28621535D+00, 0.29962788D-01, 0.30091319D-01,
449 + 0.30091319D-01, 0.29962788D-01, 0.44309512D+00, 0.15796036D-01,
450 + 0.28016920D+00, 0.28016920D+00, 0.28975574D-01, 0.29107116D-01,
451 + 0.29107116D-01, 0.28975574D-01, 0.41553632D+00, 0.15942176D-01,
452 + 0.27388348D+00, 0.27388348D+00, 0.27962754D-01, 0.28096394D-01,
453 + 0.28096394D-01, 0.27962754D-01, 0.38838946D+00, 0.15982109D-01,
454 + 0.26738416D+00, 0.26738416D+00, 0.26930282D-01, 0.27065163D-01,
455 + 0.27065163D-01, 0.26930282D-01, 0.36191917D+00, 0.15926575D-01,
456 + 0.26069742D+00, 0.26069742D+00, 0.25884124D-01, 0.26019451D-01,
457 + 0.26019451D-01, 0.25884124D-01, 0.33633600D+00, 0.15786113D-01/
458 DATA ((CALCPI(I,K,10),I=1,8), K=19,20) /
459 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
460 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
461 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
462 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
463 DATA ((CALCPI(I,K,11),I=1,8), K=1,9) /
464 + 0.33679579D+00, 0.33679579D+00, 0.20800817D-01, 0.20800817D-01,
465 + 0.20800817D-01, 0.20800817D-01, 0.55224311D+00, 0.11926224D-17,
466 + 0.33387951D+00, 0.33387951D+00, 0.20911204D-01, 0.20925082D-01,
467 + 0.20925082D-01, 0.20911204D-01, 0.55443659D+00, 0.17809761D-02,
468 + 0.33046999D+00, 0.33046999D+00, 0.20958273D-01, 0.20984584D-01,
469 + 0.20984584D-01, 0.20958273D-01, 0.55013117D+00, 0.33754846D-02,
470 + 0.32658729D+00, 0.32658729D+00, 0.20937303D-01, 0.20974641D-01,
471 + 0.20974641D-01, 0.20937303D-01, 0.54040051D+00, 0.47884497D-02,
472 + 0.32225035D+00, 0.32225035D+00, 0.20845861D-01, 0.20892872D-01,
473 + 0.20892872D-01, 0.20845861D-01, 0.52613566D+00, 0.60250537D-02,
474 + 0.31748196D+00, 0.31748196D+00, 0.20684374D-01, 0.20739765D-01,
475 + 0.20739765D-01, 0.20684374D-01, 0.50820044D+00, 0.70929073D-02,
476 + 0.31230602D+00, 0.31230602D+00, 0.20454312D-01, 0.20516855D-01,
477 + 0.20516855D-01, 0.20454312D-01, 0.48738973D+00, 0.80001390D-02,
478 + 0.30674818D+00, 0.30674818D+00, 0.20158532D-01, 0.20227071D-01,
479 + 0.20227071D-01, 0.20158532D-01, 0.46442709D+00, 0.87557943D-02,
480 + 0.30083547D+00, 0.30083547D+00, 0.19800945D-01, 0.19874398D-01,
481 + 0.19874398D-01, 0.19800945D-01, 0.43996386D+00, 0.93695568D-02/
482 DATA ((CALCPI(I,K,11),I=1,8), K=10,18) /
483 + 0.29459605D+00, 0.29459605D+00, 0.19386288D-01, 0.19463650D-01,
484 + 0.19463650D-01, 0.19386288D-01, 0.41457506D+00, 0.98515605D-02,
485 + 0.28805902D+00, 0.28805902D+00, 0.18919937D-01, 0.19000281D-01,
486 + 0.19000281D-01, 0.18919937D-01, 0.38876043D+00, 0.10212238D-01,
487 + 0.28125414D+00, 0.28125414D+00, 0.18407717D-01, 0.18490192D-01,
488 + 0.18490192D-01, 0.18407717D-01, 0.36294641D+00, 0.10462158D-01,
489 + 0.27421161D+00, 0.27421161D+00, 0.17855725D-01, 0.17939559D-01,
490 + 0.17939559D-01, 0.17855725D-01, 0.33748903D+00, 0.10611876D-01,
491 + 0.26696179D+00, 0.26696179D+00, 0.17270174D-01, 0.17354669D-01,
492 + 0.17354669D-01, 0.17270174D-01, 0.31267814D+00, 0.10671794D-01,
493 + 0.25953497D+00, 0.25953497D+00, 0.16657260D-01, 0.16741792D-01,
494 + 0.16741792D-01, 0.16657260D-01, 0.28874279D+00, 0.10652050D-01,
495 + 0.25196117D+00, 0.25196117D+00, 0.16023043D-01, 0.16107059D-01,
496 + 0.16107059D-01, 0.16023043D-01, 0.26585708D+00, 0.10562412D-01,
497 + 0.24426989D+00, 0.24426989D+00, 0.15373345D-01, 0.15456358D-01,
498 + 0.15456358D-01, 0.15373345D-01, 0.24414596D+00, 0.10412191D-01,
499 + 0.23648993D+00, 0.23648993D+00, 0.14713672D-01, 0.14795259D-01,
500 + 0.14795259D-01, 0.14713672D-01, 0.22369156D+00, 0.10210171D-01/
501 DATA ((CALCPI(I,K,11),I=1,8), K=19,20) /
502 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
503 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
504 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
505 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
506 DATA ((CALCPI(I,K,12),I=1,8), K=1,9) /
507 + 0.34023086D+00, 0.34023086D+00, 0.10652145D-01, 0.10652145D-01,
508 + 0.10652145D-01, 0.10652145D-01, 0.50042938D+00, 0.11926224D-17,
509 + 0.33524297D+00, 0.33524297D+00, 0.11089214D-01, 0.11098860D-01,
510 + 0.11098860D-01, 0.11089214D-01, 0.49097427D+00, 0.13887080D-02,
511 + 0.32977284D+00, 0.32977284D+00, 0.11429239D-01, 0.11447161D-01,
512 + 0.11447161D-01, 0.11429239D-01, 0.47678612D+00, 0.26008151D-02,
513 + 0.32385638D+00, 0.32385638D+00, 0.11675725D-01, 0.11700680D-01,
514 + 0.11700680D-01, 0.11675725D-01, 0.45892565D+00, 0.36469582D-02,
515 + 0.31752530D+00, 0.31752530D+00, 0.11832429D-01, 0.11863277D-01,
516 + 0.11863277D-01, 0.11832429D-01, 0.43823545D+00, 0.45364925D-02,
517 + 0.31081405D+00, 0.31081405D+00, 0.11905057D-01, 0.11940763D-01,
518 + 0.11940763D-01, 0.11905057D-01, 0.41549055D+00, 0.52803396D-02,
519 + 0.30375715D+00, 0.30375715D+00, 0.11899601D-01, 0.11939226D-01,
520 + 0.11939226D-01, 0.11899601D-01, 0.39137888D+00, 0.58892910D-02,
521 + 0.29638969D+00, 0.29638969D+00, 0.11822615D-01, 0.11865315D-01,
522 + 0.11865315D-01, 0.11822615D-01, 0.36649770D+00, 0.63743199D-02,
523 + 0.28874705D+00, 0.28874705D+00, 0.11681008D-01, 0.11726024D-01,
524 + 0.11726024D-01, 0.11681008D-01, 0.34135979D+00, 0.67464134D-02/
525 DATA ((CALCPI(I,K,12),I=1,8), K=10,18) /
526 + 0.28086465D+00, 0.28086465D+00, 0.11481862D-01, 0.11528519D-01,
527 + 0.11528519D-01, 0.11481862D-01, 0.31639487D+00, 0.70164284D-02,
528 + 0.27277776D+00, 0.27277776D+00, 0.11232308D-01, 0.11280009D-01,
529 + 0.11280009D-01, 0.11232308D-01, 0.29195484D+00, 0.71949907D-02,
530 + 0.26452126D+00, 0.26452126D+00, 0.10939389D-01, 0.10987608D-01,
531 + 0.10987608D-01, 0.10939389D-01, 0.26831943D+00, 0.72923849D-02,
532 + 0.25612944D+00, 0.25612944D+00, 0.10609956D-01, 0.10658235D-01,
533 + 0.10658235D-01, 0.10609956D-01, 0.24570283D+00, 0.73184702D-02,
534 + 0.24763579D+00, 0.24763579D+00, 0.10250575D-01, 0.10298520D-01,
535 + 0.10298520D-01, 0.10250575D-01, 0.22426099D+00, 0.72826065D-02,
536 + 0.23907277D+00, 0.23907277D+00, 0.98674524D-02, 0.99147235D-02,
537 + 0.99147235D-02, 0.98674524D-02, 0.20409865D+00, 0.71935922D-02,
538 + 0.23047167D+00, 0.23047167D+00, 0.94663715D-02, 0.95126842D-02,
539 + 0.95126842D-02, 0.94663715D-02, 0.18527663D+00, 0.70596202D-02,
540 + 0.22186243D+00, 0.22186243D+00, 0.90526570D-02, 0.90977743D-02,
541 + 0.90977743D-02, 0.90526570D-02, 0.16781899D+00, 0.68882501D-02,
542 + 0.21327349D+00, 0.21327349D+00, 0.86311451D-02, 0.86748738D-02,
543 + 0.86748738D-02, 0.86311451D-02, 0.15171962D+00, 0.66863917D-02/
544 DATA ((CALCPI(I,K,12),I=1,8), K=19,20) /
545 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
546 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
547 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
548 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
549 DATA ((CALCPI(I,K,13),I=1,8), K=1,9) /
550 + 0.34136380D+00, 0.34136380D+00, 0.51909346D-02, 0.51909346D-02,
551 + 0.51909346D-02, 0.51909346D-02, 0.45017673D+00, 0.22768246D-17,
552 + 0.33416865D+00, 0.33416865D+00, 0.57600044D-02, 0.57660629D-02,
553 + 0.57660629D-02, 0.57600044D-02, 0.43258651D+00, 0.10774729D-02,
554 + 0.32655596D+00, 0.32655596D+00, 0.62162263D-02, 0.62273209D-02,
555 + 0.62273209D-02, 0.62162263D-02, 0.41187469D+00, 0.19948565D-02,
556 + 0.31857365D+00, 0.31857365D+00, 0.65688772D-02, 0.65841124D-02,
557 + 0.65841124D-02, 0.65688772D-02, 0.38902390D+00, 0.27661617D-02,
558 + 0.31026278D+00, 0.31026278D+00, 0.68259178D-02, 0.68444969D-02,
559 + 0.68444969D-02, 0.68259178D-02, 0.36478314D+00, 0.34030430D-02,
560 + 0.30166582D+00, 0.30166582D+00, 0.69963705D-02, 0.70175898D-02,
561 + 0.70175898D-02, 0.69963705D-02, 0.33980354D+00, 0.39179169D-02,
562 + 0.29282433D+00, 0.29282433D+00, 0.70889405D-02, 0.71121810D-02,
563 + 0.71121810D-02, 0.70889405D-02, 0.31464086D+00, 0.43225932D-02,
564 + 0.28377922D+00, 0.28377922D+00, 0.71122362D-02, 0.71369563D-02,
565 + 0.71369563D-02, 0.71122362D-02, 0.28975167D+00, 0.46285081D-02,
566 + 0.27457061D+00, 0.27457061D+00, 0.70746759D-02, 0.71004044D-02,
567 + 0.71004044D-02, 0.70746759D-02, 0.26550532D+00, 0.48466660D-02/
568 DATA ((CALCPI(I,K,13),I=1,8), K=10,18) /
569 + 0.26523761D+00, 0.26523761D+00, 0.69843579D-02, 0.70106879D-02,
570 + 0.70106879D-02, 0.69843579D-02, 0.24218958D+00, 0.49875327D-02,
571 + 0.25581813D+00, 0.25581813D+00, 0.68489998D-02, 0.68755825D-02,
572 + 0.68755825D-02, 0.68489998D-02, 0.22001916D+00, 0.50609957D-02,
573 + 0.24634871D+00, 0.24634871D+00, 0.66758724D-02, 0.67024114D-02,
574 + 0.67024114D-02, 0.66758724D-02, 0.19914449D+00, 0.50763155D-02,
575 + 0.23686435D+00, 0.23686435D+00, 0.64717478D-02, 0.64979942D-02,
576 + 0.64979942D-02, 0.64717478D-02, 0.17966044D+00, 0.50420907D-02,
577 + 0.22739834D+00, 0.22739834D+00, 0.62428647D-02, 0.62686119D-02,
578 + 0.62686119D-02, 0.62428647D-02, 0.16161485D+00, 0.49662355D-02,
579 + 0.21798215D+00, 0.21798215D+00, 0.59949063D-02, 0.60199857D-02,
580 + 0.60199857D-02, 0.59949063D-02, 0.14501675D+00, 0.48559714D-02,
581 + 0.20864529D+00, 0.20864529D+00, 0.57329913D-02, 0.57572681D-02,
582 + 0.57572681D-02, 0.57329913D-02, 0.12984402D+00, 0.47178280D-02,
583 + 0.19941523D+00, 0.19941523D+00, 0.54616754D-02, 0.54850443D-02,
584 + 0.54850443D-02, 0.54616754D-02, 0.11605032D+00, 0.45576550D-02,
585 + 0.19031730D+00, 0.19031730D+00, 0.51849640D-02, 0.52073461D-02,
586 + 0.52073461D-02, 0.51849640D-02, 0.10357124D+00, 0.43806436D-02/
587 DATA ((CALCPI(I,K,13),I=1,8), K=19,20) /
588 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
589 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
590 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
591 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
592 DATA ((CALCPI(I,K,14),I=1,8), K=1,9) /
593 + 0.33900266D+00, 0.33900266D+00, 0.23881052D-02, 0.23881052D-02,
594 + 0.23881052D-02, 0.23881052D-02, 0.40155292D+00,-0.32526065D-18,
595 + 0.32959639D+00, 0.32959639D+00, 0.29617777D-02, 0.29653596D-02,
596 + 0.29653596D-02, 0.29617777D-02, 0.37840030D+00, 0.82753446D-03,
597 + 0.31988369D+00, 0.31988369D+00, 0.34210360D-02, 0.34275121D-02,
598 + 0.34275121D-02, 0.34210360D-02, 0.35358028D+00, 0.15147825D-02,
599 + 0.30992012D+00, 0.30992012D+00, 0.37786527D-02, 0.37874368D-02,
600 + 0.37874368D-02, 0.37786527D-02, 0.32795027D+00, 0.20773718D-02,
601 + 0.29975239D+00, 0.29975239D+00, 0.40450097D-02, 0.40555927D-02,
602 + 0.40555927D-02, 0.40450097D-02, 0.30213522D+00, 0.25278964D-02,
603 + 0.28942713D+00, 0.28942713D+00, 0.42307479D-02, 0.42426905D-02,
604 + 0.42426905D-02, 0.42307479D-02, 0.27664134D+00, 0.28790028D-02,
605 + 0.27898913D+00, 0.27898913D+00, 0.43457710D-02, 0.43586966D-02,
606 + 0.43586966D-02, 0.43457710D-02, 0.25187961D+00, 0.31424171D-02,
607 + 0.26848139D+00, 0.26848139D+00, 0.43993990D-02, 0.44129861D-02,
608 + 0.44129861D-02, 0.43993990D-02, 0.22816231D+00, 0.33291023D-02,
609 + 0.25794510D+00, 0.25794510D+00, 0.44003751D-02, 0.44143518D-02,
610 + 0.44143518D-02, 0.44003751D-02, 0.20571901D+00, 0.34492901D-02/
611 DATA ((CALCPI(I,K,14),I=1,8), K=10,18) /
612 + 0.24741943D+00, 0.24741943D+00, 0.43567934D-02, 0.43709312D-02,
613 + 0.43709312D-02, 0.43567934D-02, 0.18470551D+00, 0.35124233D-02,
614 + 0.23694143D+00, 0.23694143D+00, 0.42760971D-02, 0.42902064D-02,
615 + 0.42902064D-02, 0.42760971D-02, 0.16521479D+00, 0.35271687D-02,
616 + 0.22654590D+00, 0.22654590D+00, 0.41650665D-02, 0.41789914D-02,
617 + 0.41789914D-02, 0.41650665D-02, 0.14728752D+00, 0.35014166D-02,
618 + 0.21626529D+00, 0.21626529D+00, 0.40298181D-02, 0.40434327D-02,
619 + 0.40434327D-02, 0.40298181D-02, 0.13092160D+00, 0.34422916D-02,
620 + 0.20612961D+00, 0.20612961D+00, 0.38758163D-02, 0.38890206D-02,
621 + 0.38890206D-02, 0.38758163D-02, 0.11608121D+00, 0.33561728D-02,
622 + 0.19616637D+00, 0.19616637D+00, 0.37078918D-02, 0.37206085D-02,
623 + 0.37206085D-02, 0.37078918D-02, 0.10270492D+00, 0.32487212D-02,
624 + 0.18640053D+00, 0.18640053D+00, 0.35302685D-02, 0.35424400D-02,
625 + 0.35424400D-02, 0.35302685D-02, 0.90712788D-01, 0.31249140D-02,
626 + 0.17685449D+00, 0.17685449D+00, 0.33465970D-02, 0.33581822D-02,
627 + 0.33581822D-02, 0.33465970D-02, 0.80012495D-01, 0.29890844D-02,
628 + 0.16754808D+00, 0.16754808D+00, 0.31599922D-02, 0.31709646D-02,
629 + 0.31709646D-02, 0.31599922D-02, 0.70504500D-01, 0.28449661D-02/
630 DATA ((CALCPI(I,K,14),I=1,8), K=19,20) /
631 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
632 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
633 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
634 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
635 DATA ((CALCPI(I,K,15),I=1,8), K=1,9) /
636 + 0.33254679D+00, 0.33254679D+00, 0.10268437D-02, 0.10268437D-02,
637 + 0.10268437D-02, 0.10268437D-02, 0.35463426D+00, 0.19515639D-17,
638 + 0.32103718D+00, 0.32103718D+00, 0.15350511D-02, 0.15370048D-02,
639 + 0.15370048D-02, 0.15350511D-02, 0.32791421D+00, 0.62585318D-03,
640 + 0.30937097D+00, 0.30937097D+00, 0.19371143D-02, 0.19406061D-02,
641 + 0.19406061D-02, 0.19371143D-02, 0.30082059D+00, 0.11325320D-02,
642 + 0.29760669D+00, 0.29760669D+00, 0.22473642D-02, 0.22520477D-02,
643 + 0.22520477D-02, 0.22473642D-02, 0.27405113D+00, 0.15359343D-02,
644 + 0.28579277D+00, 0.28579277D+00, 0.24771113D-02, 0.24826921D-02,
645 + 0.24826921D-02, 0.24771113D-02, 0.24808592D+00, 0.18485763D-02,
646 + 0.27397596D+00, 0.27397596D+00, 0.26372738D-02, 0.26435029D-02,
647 + 0.26435029D-02, 0.26372738D-02, 0.22327531D+00, 0.20824543D-02,
648 + 0.26220040D+00, 0.26220040D+00, 0.27377875D-02, 0.27444559D-02,
649 + 0.27444559D-02, 0.27377875D-02, 0.19988208D+00, 0.22484764D-02,
650 + 0.25050737D+00, 0.25050737D+00, 0.27876780D-02, 0.27946119D-02,
651 + 0.27946119D-02, 0.27876780D-02, 0.17807768D+00, 0.23565372D-02,
652 + 0.23893543D+00, 0.23893543D+00, 0.27951652D-02, 0.28022210D-02,
653 + 0.28022210D-02, 0.27951652D-02, 0.15796116D+00, 0.24156353D-02/
654 DATA ((CALCPI(I,K,15),I=1,8), K=10,18) /
655 + 0.22752029D+00, 0.22752029D+00, 0.27676405D-02, 0.27747007D-02,
656 + 0.27747007D-02, 0.27676405D-02, 0.13957089D+00, 0.24338611D-02,
657 + 0.21629474D+00, 0.21629474D+00, 0.27117073D-02, 0.27186776D-02,
658 + 0.27186776D-02, 0.27117073D-02, 0.12289610D+00, 0.24184458D-02,
659 + 0.20528864D+00, 0.20528864D+00, 0.26332162D-02, 0.26400217D-02,
660 + 0.26400217D-02, 0.26332162D-02, 0.10788808D+00, 0.23758034D-02,
661 + 0.19452886D+00, 0.19452886D+00, 0.25373016D-02, 0.25438844D-02,
662 + 0.25438844D-02, 0.25373016D-02, 0.94469725D-01, 0.23115748D-02,
663 + 0.18403930D+00, 0.18403930D+00, 0.24284260D-02, 0.24347423D-02,
664 + 0.24347423D-02, 0.24284260D-02, 0.82543964D-01, 0.22306770D-02,
665 + 0.17384087D+00, 0.17384087D+00, 0.23104266D-02, 0.23164449D-02,
666 + 0.23164449D-02, 0.23104266D-02, 0.72000989D-01, 0.21373554D-02,
667 + 0.16395161D+00, 0.16395161D+00, 0.21865656D-02, 0.21922646D-02,
668 + 0.21922646D-02, 0.21865656D-02, 0.62724181D-01, 0.20352383D-02,
669 + 0.15438665D+00, 0.15438665D+00, 0.20595809D-02, 0.20649478D-02,
670 + 0.20649478D-02, 0.20595809D-02, 0.54594861D-01, 0.19273921D-02,
671 + 0.14515837D+00, 0.14515837D+00, 0.19317383D-02, 0.19367675D-02,
672 + 0.19367675D-02, 0.19317383D-02, 0.47496076D-01, 0.18163768D-02/
673 DATA ((CALCPI(I,K,15),I=1,8), K=19,20) /
674 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
675 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
676 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
677 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
678 DATA ((CALCPI(I,K,16),I=1,8), K=1,9) /
679 + 0.32176997D+00, 0.32176997D+00, 0.40737579D-03, 0.40737579D-03,
680 + 0.40737579D-03, 0.40737579D-03, 0.30950802D+00,-0.30493186D-18,
681 + 0.30835580D+00, 0.30835580D+00, 0.82070349D-03, 0.82168504D-03,
682 + 0.82168504D-03, 0.82070349D-03, 0.28083323D+00, 0.46338427D-03,
683 + 0.29496614D+00, 0.29496614D+00, 0.11419643D-02, 0.11436993D-02,
684 + 0.11436993D-02, 0.11419643D-02, 0.25291374D+00, 0.82871451D-03,
685 + 0.28165761D+00, 0.28165761D+00, 0.13855121D-02, 0.13878144D-02,
686 + 0.13878144D-02, 0.13855121D-02, 0.22626249D+00, 0.11111238D-02,
687 + 0.26847618D+00, 0.26847618D+00, 0.15624583D-02, 0.15651727D-02,
688 + 0.15651727D-02, 0.15624583D-02, 0.20120250D+00, 0.13223185D-02,
689 + 0.25546461D+00, 0.25546461D+00, 0.16830218D-02, 0.16860197D-02,
690 + 0.16860197D-02, 0.16830218D-02, 0.17792501D+00, 0.14730419D-02,
691 + 0.24266237D+00, 0.24266237D+00, 0.17563436D-02, 0.17595194D-02,
692 + 0.17595194D-02, 0.17563436D-02, 0.15654724D+00, 0.15729050D-02,
693 + 0.23010523D+00, 0.23010523D+00, 0.17904961D-02, 0.17937637D-02,
694 + 0.17937637D-02, 0.17904961D-02, 0.13710931D+00, 0.16303955D-02,
695 + 0.21782553D+00, 0.21782553D+00, 0.17926444D-02, 0.17959346D-02,
696 + 0.17959346D-02, 0.17926444D-02, 0.11959379D+00, 0.16530474D-02/
697 DATA ((CALCPI(I,K,16),I=1,8), K=10,18) /
698 + 0.20585208D+00, 0.20585208D+00, 0.17690681D-02, 0.17723261D-02,
699 + 0.17723261D-02, 0.17690681D-02, 0.10393872D+00, 0.16474683D-02,
700 + 0.19421024D+00, 0.19421024D+00, 0.17252349D-02, 0.17284179D-02,
701 + 0.17284179D-02, 0.17252349D-02, 0.90049126D-01, 0.16194169D-02,
702 + 0.18292194D+00, 0.18292194D+00, 0.16658653D-02, 0.16689407D-02,
703 + 0.16689407D-02, 0.16658653D-02, 0.77807700D-01, 0.15738726D-02,
704 + 0.17200575D+00, 0.17200575D+00, 0.15949942D-02, 0.15979379D-02,
705 + 0.15979379D-02, 0.15949942D-02, 0.67083270D-01, 0.15151004D-02,
706 + 0.16147698D+00, 0.16147698D+00, 0.15160351D-02, 0.15188303D-02,
707 + 0.15188303D-02, 0.15160351D-02, 0.57737942D-01, 0.14467182D-02,
708 + 0.15134775D+00, 0.15134775D+00, 0.14318425D-02, 0.14344781D-02,
709 + 0.14344781D-02, 0.14318425D-02, 0.49632713D-01, 0.13717621D-02,
710 + 0.14162717D+00, 0.14162717D+00, 0.13447724D-02, 0.13472423D-02,
711 + 0.13472423D-02, 0.13447724D-02, 0.42631818D-01, 0.12927493D-02,
712 + 0.13232145D+00, 0.13232145D+00, 0.12567414D-02, 0.12590432D-02,
713 + 0.12590432D-02, 0.12567414D-02, 0.36605948D-01, 0.12117391D-02,
714 + 0.12343408D+00, 0.12343408D+00, 0.11692815D-02, 0.11714161D-02,
715 + 0.11714161D-02, 0.11692815D-02, 0.31434284D-01, 0.11303904D-02/
716 DATA ((CALCPI(I,K,16),I=1,8), K=19,20) /
717 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
718 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
719 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
720 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
721 DATA ((CALCPI(I,K,17),I=1,8), K=1,9) /
722 + 0.30666522D+00, 0.30666522D+00, 0.14660455D-03, 0.14660455D-03,
723 + 0.14660455D-03, 0.14660455D-03, 0.26627487D+00, 0.00000000D+00,
724 + 0.29162048D+00, 0.29162048D+00, 0.46084555D-03, 0.46129484D-03,
725 + 0.46129484D-03, 0.46084555D-03, 0.23698931D+00, 0.33349774D-03,
726 + 0.27680572D+00, 0.27680572D+00, 0.69973234D-03, 0.70051788D-03,
727 + 0.70051788D-03, 0.69973234D-03, 0.20941648D+00, 0.58916459D-03,
728 + 0.26227037D+00, 0.26227037D+00, 0.87652031D-03, 0.87755172D-03,
729 + 0.87755172D-03, 0.87652031D-03, 0.18386799D+00, 0.78060308D-03,
730 + 0.24805358D+00, 0.24805358D+00, 0.10013164D-02, 0.10025198D-02,
731 + 0.10025198D-02, 0.10013164D-02, 0.16050388D+00, 0.91818310D-03,
732 + 0.23418982D+00, 0.23418982D+00, 0.10830110D-02, 0.10843263D-02,
733 + 0.10843263D-02, 0.10830110D-02, 0.13935998D+00, 0.10110233D-02,
734 + 0.22070989D+00, 0.22070989D+00, 0.11294484D-02, 0.11308274D-02,
735 + 0.11308274D-02, 0.11294484D-02, 0.12041707D+00, 0.10671712D-02,
736 + 0.20764028D+00, 0.20764028D+00, 0.11473733D-02, 0.11487776D-02,
737 + 0.11487776D-02, 0.11473733D-02, 0.10359781D+00, 0.10935493D-02,
738 + 0.19500357D+00, 0.19500357D+00, 0.11426309D-02, 0.11440303D-02,
739 + 0.11440303D-02, 0.11426309D-02, 0.88785455D-01, 0.10961590D-02/
740 DATA ((CALCPI(I,K,17),I=1,8), K=10,18) /
741 + 0.18281847D+00, 0.18281847D+00, 0.11202227D-02, 0.11215942D-02,
742 + 0.11215942D-02, 0.11202227D-02, 0.75837338D-01, 0.10801391D-02,
743 + 0.17109992D+00, 0.17109992D+00, 0.10843932D-02, 0.10857193D-02,
744 + 0.10857193D-02, 0.10843932D-02, 0.64594818D-01, 0.10498550D-02,
745 + 0.15985929D+00, 0.15985929D+00, 0.10387136D-02, 0.10399817D-02,
746 + 0.10399817D-02, 0.10387136D-02, 0.54892435D-01, 0.10089842D-02,
747 + 0.14910448D+00, 0.14910448D+00, 0.98615555D-03, 0.98735689D-03,
748 + 0.98735689D-03, 0.98615555D-03, 0.46564629D-01, 0.96059177D-03,
749 + 0.13884016D+00, 0.13884016D+00, 0.92916273D-03, 0.93029172D-03,
750 + 0.93029172D-03, 0.92916273D-03, 0.39450894D-01, 0.90720346D-03,
751 + 0.12906792D+00, 0.12906792D+00, 0.86971779D-03, 0.87077138D-03,
752 + 0.87077138D-03, 0.86971779D-03, 0.33399614D-01, 0.85087411D-03,
753 + 0.11978651D+00, 0.11978651D+00, 0.80940421D-03, 0.81038137D-03,
754 + 0.81038137D-03, 0.80940421D-03, 0.28270305D-01, 0.79325052D-03,
755 + 0.11099203D+00, 0.11099203D+00, 0.74946350D-03, 0.75036480D-03,
756 + 0.75036480D-03, 0.74946350D-03, 0.23934948D-01, 0.73562974D-03,
757 + 0.10267820D+00, 0.10267820D+00, 0.69084749D-03, 0.69167470D-03,
758 + 0.69167470D-03, 0.69084749D-03, 0.20278901D-01, 0.67901228D-03/
759 DATA ((CALCPI(I,K,17),I=1,8), K=19,20) /
760 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
761 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
762 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
763 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
764 DATA ((CALCPI(I,K,18),I=1,8), K=1,9) /
765 + 0.28734139D+00, 0.28734139D+00, 0.46770052D-04, 0.46770052D-04,
766 + 0.46770052D-04, 0.46770052D-04, 0.22505275D+00,-0.30704944D-18,
767 + 0.27100656D+00, 0.27100656D+00, 0.27133627D-03, 0.27152066D-03,
768 + 0.27152066D-03, 0.27133627D-03, 0.19630108D+00, 0.23112086D-03,
769 + 0.25512465D+00, 0.25512465D+00, 0.43763189D-03, 0.43795093D-03,
770 + 0.43795093D-03, 0.43763189D-03, 0.17003963D+00, 0.40306695D-03,
771 + 0.23973226D+00, 0.23973226D+00, 0.55705926D-03, 0.55747395D-03,
772 + 0.55747395D-03, 0.55705926D-03, 0.14636982D+00, 0.52737673D-03,
773 + 0.22485705D+00, 0.22485705D+00, 0.63822452D-03, 0.63870357D-03,
774 + 0.63870357D-03, 0.63822452D-03, 0.12529003D+00, 0.61275852D-03,
775 + 0.21052083D+00, 0.21052083D+00, 0.68834581D-03, 0.68886420D-03,
776 + 0.68886420D-03, 0.68834581D-03, 0.10669065D+00, 0.66651780D-03,
777 + 0.19674164D+00, 0.19674164D+00, 0.71372079D-03, 0.71425884D-03,
778 + 0.71425884D-03, 0.71372079D-03, 0.90429831D-01, 0.69502936D-03,
779 + 0.18353298D+00, 0.18353298D+00, 0.71963153D-03, 0.72017396D-03,
780 + 0.72017396D-03, 0.71963153D-03, 0.76331198D-01, 0.70364208D-03,
781 + 0.17090433D+00, 0.17090433D+00, 0.71054453D-03, 0.71107968D-03,
782 + 0.71107968D-03, 0.71054453D-03, 0.64199063D-01, 0.69688051D-03/
783 DATA ((CALCPI(I,K,18),I=1,8), K=10,18) /
784 + 0.15886127D+00, 0.15886127D+00, 0.69018951D-03, 0.69070869D-03,
785 + 0.69070869D-03, 0.69018951D-03, 0.53831081D-01, 0.67852490D-03,
786 + 0.14740574D+00, 0.14740574D+00, 0.66165007D-03, 0.66214700D-03,
787 + 0.66214700D-03, 0.66165007D-03, 0.45025818D-01, 0.65170278D-03,
788 + 0.13653626D+00, 0.13653626D+00, 0.62745152D-03, 0.62792190D-03,
789 + 0.62792190D-03, 0.62745152D-03, 0.37589423D-01, 0.61897771D-03,
790 + 0.12624821D+00, 0.12624821D+00, 0.58963464D-03, 0.59007572D-03,
791 + 0.59007572D-03, 0.58963464D-03, 0.31340230D-01, 0.58242370D-03,
792 + 0.11653405D+00, 0.11653405D+00, 0.54982595D-03, 0.55023624D-03,
793 + 0.55023624D-03, 0.54982595D-03, 0.26111312D-01, 0.54369618D-03,
794 + 0.10738366D+00, 0.10738366D+00, 0.50930001D-03, 0.50967898D-03,
795 + 0.50967898D-03, 0.50930001D-03, 0.21751783D-01, 0.50409480D-03,
796 + 0.98784558D-01, 0.98784558D-01, 0.46903555D-03, 0.46938343D-03,
797 + 0.46938343D-03, 0.46903555D-03, 0.18127816D-01, 0.46462009D-03,
798 + 0.90722201D-01, 0.90722201D-01, 0.42976572D-03, 0.43008330D-03,
799 + 0.43008330D-03, 0.42976572D-03, 0.15122287D-01, 0.42602406D-03,
800 + 0.83180274D-01, 0.83180274D-01, 0.39202129D-03, 0.39230975D-03,
801 + 0.39230975D-03, 0.39202129D-03, 0.12633363D-01, 0.38885385D-03/
802 DATA ((CALCPI(I,K,18),I=1,8), K=19,20) /
803 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
804 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
805 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
806 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
807 DATA ((CALCPI(I,K,19),I=1,8), K=1,9) /
808 + 0.26395866D+00, 0.26395866D+00, 0.12806974D-04, 0.12806974D-04,
809 + 0.12806974D-04, 0.12806974D-04, 0.18598241D+00, 0.11536589D-17,
810 + 0.24673919D+00, 0.24673919D+00, 0.16309166D-03, 0.16316063D-03,
811 + 0.16316063D-03, 0.16309166D-03, 0.15875377D+00, 0.15219856D-03,
812 + 0.23020567D+00, 0.23020567D+00, 0.27106979D-03, 0.27118771D-03,
813 + 0.27118771D-03, 0.27106979D-03, 0.13460231D+00, 0.26180772D-03,
814 + 0.21437565D+00, 0.21437565D+00, 0.34585977D-03, 0.34601127D-03,
815 + 0.34601127D-03, 0.34585977D-03, 0.11341881D+00, 0.33799152D-03,
816 + 0.19926040D+00, 0.19926040D+00, 0.39430293D-03, 0.39447595D-03,
817 + 0.39447595D-03, 0.39430293D-03, 0.95047207D-01, 0.38762528D-03,
818 + 0.18486454D+00, 0.18486454D+00, 0.42184129D-03, 0.42202640D-03,
819 + 0.42202640D-03, 0.42184129D-03, 0.79249293D-01, 0.41617954D-03,
820 + 0.17118932D+00, 0.17118932D+00, 0.43318564D-03, 0.43337558D-03,
821 + 0.43337558D-03, 0.43318564D-03, 0.65780106D-01, 0.42839007D-03,
822 + 0.15823168D+00, 0.15823168D+00, 0.43218780D-03, 0.43237712D-03,
823 + 0.43237712D-03, 0.43218780D-03, 0.54386231D-01, 0.42813014D-03,
824 + 0.14598490D+00, 0.14598490D+00, 0.42202134D-03, 0.42220599D-03,
825 + 0.42220599D-03, 0.42202134D-03, 0.44815762D-01, 0.41859167D-03/
826 DATA ((CALCPI(I,K,19),I=1,8), K=10,18) /
827 + 0.13443884D+00, 0.13443884D+00, 0.40527228D-03, 0.40544938D-03,
828 + 0.40544938D-03, 0.40527228D-03, 0.36828764D-01, 0.40237655D-03,
829 + 0.12358024D+00, 0.12358024D+00, 0.38401842D-03, 0.38418600D-03,
830 + 0.38418600D-03, 0.38401842D-03, 0.30202058D-01, 0.38157616D-03,
831 + 0.11339309D+00, 0.11339309D+00, 0.35991167D-03, 0.36006849D-03,
832 + 0.36006849D-03, 0.35991167D-03, 0.24732435D-01, 0.35785413D-03,
833 + 0.10385892D+00, 0.10385892D+00, 0.33424401D-03, 0.33438939D-03,
834 + 0.33438939D-03, 0.33424401D-03, 0.20237777D-01, 0.33251249D-03,
835 + 0.94957156D-01, 0.94957156D-01, 0.30800812D-03, 0.30814180D-03,
836 + 0.30814180D-03, 0.30800812D-03, 0.16558265D-01, 0.30655254D-03,
837 + 0.86665412D-01, 0.86665412D-01, 0.28195020D-03, 0.28207228D-03,
838 + 0.28207228D-03, 0.28195020D-03, 0.13555658D-01, 0.28072794D-03,
839 + 0.78959834D-01, 0.78959834D-01, 0.25661472D-03, 0.25672549D-03,
840 + 0.25672549D-03, 0.25661472D-03, 0.11110861D-01, 0.25558947D-03,
841 + 0.71815395D-01, 0.71815395D-01, 0.23238415D-03, 0.23248412D-03,
842 + 0.23248412D-03, 0.23238415D-03, 0.91230074D-02, 0.23152509D-03,
843 + 0.65206193D-01, 0.65206193D-01, 0.20951288D-03, 0.20960264D-03,
844 + 0.20960264D-03, 0.20951288D-03, 0.75082135D-02, 0.20879382D-03/
845 DATA ((CALCPI(I,K,19),I=1,8), K=19,20) /
846 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
847 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
848 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
849 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
850 DATA ((CALCPI(I,K,20),I=1,8), K=1,9) /
851 + 0.23669127D+00, 0.23669127D+00, 0.28711941D-05, 0.28711941D-05,
852 + 0.28711941D-05, 0.28711941D-05, 0.14923539D+00, 0.11132925D-17,
853 + 0.21906335D+00, 0.21906335D+00, 0.95742296D-04, 0.95763866D-04,
854 + 0.95763866D-04, 0.95742296D-04, 0.12439138D+00, 0.93328989D-04,
855 + 0.20235545D+00, 0.20235545D+00, 0.16022374D-03, 0.16026014D-03,
856 + 0.16026014D-03, 0.16022374D-03, 0.10300729D+00, 0.15819578D-03,
857 + 0.18655938D+00, 0.18655938D+00, 0.20299406D-03, 0.20304025D-03,
858 + 0.20304025D-03, 0.20299406D-03, 0.84772200D-01, 0.20129147D-03,
859 + 0.17166484D+00, 0.17166484D+00, 0.22907545D-03, 0.22912757D-03,
860 + 0.22912757D-03, 0.22907545D-03, 0.69390365D-01, 0.22764750D-03,
861 + 0.15765477D+00, 0.15765477D+00, 0.24221238D-03, 0.24226746D-03,
862 + 0.24226746D-03, 0.24221238D-03, 0.56518865D-01, 0.24101595D-03,
863 + 0.14450970D+00, 0.14450970D+00, 0.24563633D-03, 0.24569217D-03,
864 + 0.24569217D-03, 0.24563633D-03, 0.45834373D-01, 0.24463494D-03,
865 + 0.13220682D+00, 0.13220682D+00, 0.24192611D-03, 0.24198109D-03,
866 + 0.24198109D-03, 0.24192611D-03, 0.37031801D-01, 0.24108886D-03,
867 + 0.12072057D+00, 0.12072057D+00, 0.23314480D-03, 0.23319779D-03,
868 + 0.23319779D-03, 0.23314480D-03, 0.29828119D-01, 0.23244556D-03/
869 DATA ((CALCPI(I,K,20),I=1,8), K=10,18) /
870 + 0.11002305D+00, 0.11002305D+00, 0.22093141D-03, 0.22098162D-03,
871 + 0.22098162D-03, 0.22093141D-03, 0.23969255D-01, 0.22034808D-03,
872 + 0.10008435D+00, 0.10008435D+00, 0.20656140D-03, 0.20660833D-03,
873 + 0.20660833D-03, 0.20656140D-03, 0.19229493D-01, 0.20607532D-03,
874 + 0.90872956D-01, 0.90872956D-01, 0.19101519D-03, 0.19105859D-03,
875 + 0.19105859D-03, 0.19101519D-03, 0.15413119D-01, 0.19061061D-03,
876 + 0.82356092D-01, 0.82356092D-01, 0.17503155D-03, 0.17507129D-03,
877 + 0.17507129D-03, 0.17503155D-03, 0.12353238D-01, 0.17469519D-03,
878 + 0.74500077D-01, 0.74500077D-01, 0.15915224D-03, 0.15918834D-03,
879 + 0.15918834D-03, 0.15915224D-03, 0.99074948D-02, 0.15887290D-03,
880 + 0.67270643D-01, 0.67270643D-01, 0.14376248D-03, 0.14379504D-03,
881 + 0.14379504D-03, 0.14376248D-03, 0.79568139D-02, 0.14353077D-03,
882 + 0.60633259D-01, 0.60633259D-01, 0.12912489D-03, 0.12915408D-03,
883 + 0.12915408D-03, 0.12912489D-03, 0.64039421D-02, 0.12893290D-03,
884 + 0.54553407D-01, 0.54553407D-01, 0.11540590D-03, 0.11543192D-03,
885 + 0.11543192D-03, 0.11540590D-03, 0.51687923D-02, 0.11524699D-03,
886 + 0.48996851D-01, 0.48996851D-01, 0.10269847D-03, 0.10272155D-03,
887 + 0.10272155D-03, 0.10269847D-03, 0.41856346D-02, 0.10256709D-03/
888 DATA ((CALCPI(I,K,20),I=1,8), K=19,20) /
889 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
890 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
891 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
892 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
893 DATA ((CALCPI(I,K,21),I=1,8), K=1,9) /
894 + 0.20570850D+00, 0.20570850D+00, 0.48978368D-06, 0.48978368D-06,
895 + 0.48978368D-06, 0.48978368D-06, 0.11502730D+00,-0.42248415D-18,
896 + 0.18823309D+00, 0.18823309D+00, 0.51928064D-04, 0.51933292D-04,
897 + 0.51933292D-04, 0.51928064D-04, 0.93318280D-01, 0.51521842D-04,
898 + 0.17190066D+00, 0.17190066D+00, 0.86309571D-04, 0.86318273D-04,
899 + 0.86318273D-04, 0.86309571D-04, 0.75230196D-01, 0.85972687D-04,
900 + 0.15667015D+00, 0.15667015D+00, 0.10796120D-03, 0.10797209D-03,
901 + 0.10797209D-03, 0.10796120D-03, 0.60271024D-01, 0.10768207D-03,
902 + 0.14250440D+00, 0.14250440D+00, 0.12019981D-03, 0.12021194D-03,
903 + 0.12021194D-03, 0.12019981D-03, 0.48029010D-01, 0.11996879D-03,
904 + 0.12936043D+00, 0.12936043D+00, 0.12529323D-03, 0.12530589D-03,
905 + 0.12530589D-03, 0.12529323D-03, 0.38087955D-01, 0.12510223D-03,
906 + 0.11719466D+00, 0.11719466D+00, 0.12521077D-03, 0.12522344D-03,
907 + 0.12522344D-03, 0.12521077D-03, 0.30076498D-01, 0.12505303D-03,
908 + 0.10596194D+00, 0.10596194D+00, 0.12149382D-03, 0.12150615D-03,
909 + 0.12150615D-03, 0.12149382D-03, 0.23668193D-01, 0.12136369D-03,
910 + 0.95616090D-01, 0.95616090D-01, 0.11533162D-03, 0.11534335D-03,
911 + 0.11534335D-03, 0.11533162D-03, 0.18573882D-01, 0.11522439D-03/
912 DATA ((CALCPI(I,K,21),I=1,8), K=10,18) /
913 + 0.86110368D-01, 0.86110368D-01, 0.10764370D-03, 0.10765467D-03,
914 + 0.10765467D-03, 0.10764370D-03, 0.14547257D-01, 0.10755544D-03,
915 + 0.77397783D-01, 0.77397783D-01, 0.99122127D-04, 0.99132248D-04,
916 + 0.99132248D-04, 0.99122127D-04, 0.11381933D-01, 0.99049573D-04,
917 + 0.69431467D-01, 0.69431467D-01, 0.90277093D-04, 0.90286330D-04,
918 + 0.90286330D-04, 0.90277093D-04, 0.89038321D-02, 0.90217518D-04,
919 + 0.62164993D-01, 0.62164993D-01, 0.81474568D-04, 0.81482918D-04,
920 + 0.81482918D-04, 0.81474568D-04, 0.69694791D-02, 0.81425707D-04,
921 + 0.55552672D-01, 0.55552672D-01, 0.72968726D-04, 0.72976215D-04,
922 + 0.72976215D-04, 0.72968726D-04, 0.54642984D-02, 0.72928701D-04,
923 + 0.49549820D-01, 0.49549820D-01, 0.64926352D-04, 0.64933021D-04,
924 + 0.64933021D-04, 0.64926352D-04, 0.42953903D-02, 0.64893603D-04,
925 + 0.44112993D-01, 0.44112993D-01, 0.57447545D-04, 0.57453446D-04,
926 + 0.57453446D-04, 0.57447545D-04, 0.33874173D-02, 0.57420781D-04,
927 + 0.39200192D-01, 0.39200192D-01, 0.50584221D-04, 0.50589413D-04,
928 + 0.50589413D-04, 0.50584221D-04, 0.26818876D-02, 0.50562372D-04,
929 + 0.34771031D-01, 0.34771031D-01, 0.44352940D-04, 0.44357486D-04,
930 + 0.44357486D-04, 0.44352940D-04, 0.21333837D-02, 0.44335125D-04/
931 DATA ((CALCPI(I,K,21),I=1,8), K=19,20) /
932 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
933 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
934 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
935 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
936 DATA ((CALCPI(I,K,22),I=1,8), K=1,9) /
937 + 0.17116626D+00, 0.17116626D+00, 0.56231149D-07, 0.56231149D-07,
938 + 0.56231149D-07, 0.56231149D-07, 0.83639872D-01, 0.76929286D-18,
939 + 0.15451487D+00, 0.15451487D+00, 0.24052761D-04, 0.24053650D-04,
940 + 0.24053650D-04, 0.24052761D-04, 0.65710784D-01, 0.24006843D-04,
941 + 0.13920030D+00, 0.13920030D+00, 0.39474321D-04, 0.39475779D-04,
942 + 0.39475779D-04, 0.39474321D-04, 0.51319235D-01, 0.39436821D-04,
943 + 0.12514136D+00, 0.12514136D+00, 0.48604732D-04, 0.48606532D-04,
944 + 0.48606532D-04, 0.48604732D-04, 0.39828486D-01, 0.48574135D-04,
945 + 0.11226909D+00, 0.11226909D+00, 0.53305228D-04, 0.53307205D-04,
946 + 0.53307205D-04, 0.53305228D-04, 0.30744208D-01, 0.53280294D-04,
947 + 0.10051127D+00, 0.10051127D+00, 0.54693175D-04, 0.54695209D-04,
948 + 0.54695209D-04, 0.54693175D-04, 0.23620324D-01, 0.54672877D-04,
949 + 0.89797814D-01, 0.89797814D-01, 0.53768172D-04, 0.53770179D-04,
950 + 0.53770179D-04, 0.53768172D-04, 0.18072063D-01, 0.53751666D-04,
951 + 0.80059998D-01, 0.80059998D-01, 0.51307172D-04, 0.51309096D-04,
952 + 0.51309096D-04, 0.51307172D-04, 0.13781706D-01, 0.51293766D-04,
953 + 0.71230720D-01, 0.71230720D-01, 0.47883546D-04, 0.47885351D-04,
954 + 0.47885351D-04, 0.47883546D-04, 0.10486468D-01, 0.47872670D-04/
955 DATA ((CALCPI(I,K,22),I=1,8), K=10,18) /
956 + 0.63244892D-01, 0.63244892D-01, 0.43928693D-04, 0.43930358D-04,
957 + 0.43930358D-04, 0.43928693D-04, 0.79685257D-02, 0.43919881D-04,
958 + 0.56039652D-01, 0.56039652D-01, 0.39753609D-04, 0.39755124D-04,
959 + 0.39755124D-04, 0.39753609D-04, 0.60516870D-02, 0.39746478D-04,
960 + 0.49554577D-01, 0.49554577D-01, 0.35578482D-04, 0.35579846D-04,
961 + 0.35579846D-04, 0.35578482D-04, 0.45992228D-02, 0.35572719D-04,
962 + 0.43731867D-01, 0.43731867D-01, 0.31551726D-04, 0.31552941D-04,
963 + 0.31552941D-04, 0.31551726D-04, 0.35024657D-02, 0.31547074D-04,
964 + 0.38516480D-01, 0.38516480D-01, 0.27766366D-04, 0.27767441D-04,
965 + 0.27767441D-04, 0.27766366D-04, 0.26745429D-02, 0.27762615D-04,
966 + 0.33856245D-01, 0.33856245D-01, 0.24276734D-04, 0.24277678D-04,
967 + 0.24277678D-04, 0.24276734D-04, 0.20499269D-02, 0.24273714D-04,
968 + 0.29701933D-01, 0.29701933D-01, 0.21108366D-04, 0.21109189D-04,
969 + 0.21109189D-04, 0.21108366D-04, 0.15791350D-02, 0.21105937D-04,
970 + 0.26007296D-01, 0.26007296D-01, 0.18265891D-04, 0.18266605D-04,
971 + 0.18266605D-04, 0.18265891D-04, 0.12234050D-02, 0.18263940D-04,
972 + 0.22729072D-01, 0.22729072D-01, 0.15740611D-04, 0.15741227D-04,
973 + 0.15741227D-04, 0.15740611D-04, 0.95351049D-03, 0.15739045D-04/
974 DATA ((CALCPI(I,K,22),I=1,8), K=19,20) /
975 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
976 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
977 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
978 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
979 DATA ((CALCPI(I,K,23),I=1,8), K=1,9) /
980 + 0.13320499D+00, 0.13320499D+00, 0.34520002D-08, 0.34520002D-08,
981 + 0.34520002D-08, 0.34520002D-08, 0.55462664D-01, 0.20849354D-17,
982 + 0.11820636D+00, 0.11820636D+00, 0.80765048D-05, 0.80765896D-05,
983 + 0.80765896D-05, 0.80765048D-05, 0.41845937D-01, 0.80737395D-05,
984 + 0.10467977D+00, 0.10467977D+00, 0.13146773D-04, 0.13146909D-04,
985 + 0.13146909D-04, 0.13146773D-04, 0.31407747D-01, 0.13144557D-04,
986 + 0.92497697D-01, 0.92497697D-01, 0.15958984D-04, 0.15959147D-04,
987 + 0.15959147D-04, 0.15958984D-04, 0.23428565D-01, 0.15957210D-04,
988 + 0.81555755D-01, 0.81555755D-01, 0.17284448D-04, 0.17284623D-04,
989 + 0.17284623D-04, 0.17284448D-04, 0.17377487D-01, 0.17283030D-04,
990 + 0.71751325D-01, 0.71751325D-01, 0.17478843D-04, 0.17479019D-04,
991 + 0.17479019D-04, 0.17478843D-04, 0.12828131D-01, 0.17477710D-04,
992 + 0.62987915D-01, 0.62987915D-01, 0.16901157D-04, 0.16901326D-04,
993 + 0.16901326D-04, 0.16901157D-04, 0.94346127D-02, 0.16900253D-04,
994 + 0.55174530D-01, 0.55174530D-01, 0.15846998D-04, 0.15847157D-04,
995 + 0.15847157D-04, 0.15846998D-04, 0.69202211D-02, 0.15846278D-04,
996 + 0.48225545D-01, 0.48225545D-01, 0.14516072D-04, 0.14516218D-04,
997 + 0.14516218D-04, 0.14516072D-04, 0.50647158D-02, 0.14515499D-04/
998 DATA ((CALCPI(I,K,23),I=1,8), K=10,18) /
999 + 0.42060787D-01, 0.42060787D-01, 0.13059574D-04, 0.13059706D-04,
1000 + 0.13059706D-04, 0.13059574D-04, 0.37042547D-02, 0.13059118D-04,
1001 + 0.36605408D-01, 0.36605408D-01, 0.11583180D-04, 0.11583299D-04,
1002 + 0.11583299D-04, 0.11583180D-04, 0.27122284D-02, 0.11582819D-04,
1003 + 0.31789796D-01, 0.31789796D-01, 0.10154824D-04, 0.10154929D-04,
1004 + 0.10154929D-04, 0.10154824D-04, 0.19892978D-02, 0.10154537D-04,
1005 + 0.27549461D-01, 0.27549461D-01, 0.88180035D-05, 0.88180962D-05,
1006 + 0.88180962D-05, 0.88180035D-05, 0.14633796D-02, 0.88177766D-05,
1007 + 0.23824875D-01, 0.23824875D-01, 0.75968481D-05, 0.75969293D-05,
1008 + 0.75969293D-05, 0.75968481D-05, 0.10818880D-02, 0.75966687D-05,
1009 + 0.20561299D-01, 0.20561299D-01, 0.65012925D-05, 0.65013633D-05,
1010 + 0.65013633D-05, 0.65012925D-05, 0.80466039D-03, 0.65011509D-05,
1011 + 0.17708592D-01, 0.17708592D-01, 0.55324326D-05, 0.55324941D-05,
1012 + 0.55324941D-05, 0.55324326D-05, 0.60243422D-03, 0.55323209D-05,
1013 + 0.15220996D-01, 0.15220996D-01, 0.46854441D-05, 0.46854973D-05,
1014 + 0.46854973D-05, 0.46854441D-05, 0.45457813D-03, 0.46853561D-05,
1015 + 0.13056907D-01, 0.13056907D-01, 0.39517963D-05, 0.39518422D-05,
1016 + 0.39518422D-05, 0.39517963D-05, 0.34596439D-03, 0.39517270D-05/
1017 DATA ((CALCPI(I,K,23),I=1,8), K=19,20) /
1018 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
1019 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
1020 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
1021 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
1022 DATA ((CALCPI(I,K,24),I=1,8), K=1,9) /
1023 + 0.91950370D-01, 0.91950370D-01, 0.67606324D-10, 0.67606316D-10,
1024 + 0.67606316D-10, 0.67606324D-10, 0.31084406D-01, 0.20406171D-17,
1025 + 0.79684852D-01, 0.79684852D-01, 0.91576626D-06, 0.91575163D-06,
1026 + 0.91575163D-06, 0.91576626D-06, 0.22168331D-01, 0.91571349D-06,
1027 + 0.68912819D-01, 0.68912819D-01, 0.16482411D-05, 0.16482201D-05,
1028 + 0.16482201D-05, 0.16482411D-05, 0.15754789D-01, 0.16481998D-05,
1029 + 0.59458278D-01, 0.59458278D-01, 0.20800602D-05, 0.20800381D-05,
1030 + 0.20800381D-05, 0.20800602D-05, 0.11132200D-01, 0.20800281D-05,
1031 + 0.51182817D-01, 0.51182817D-01, 0.23331048D-05, 0.23330847D-05,
1032 + 0.23330847D-05, 0.23331048D-05, 0.78169614D-02, 0.23330797D-05,
1033 + 0.43957520D-01, 0.43957520D-01, 0.24064376D-05, 0.24064213D-05,
1034 + 0.24064213D-05, 0.24064376D-05, 0.54678262D-02, 0.24064181D-05,
1035 + 0.37664877D-01, 0.37664877D-01, 0.23357250D-05, 0.23357134D-05,
1036 + 0.23357134D-05, 0.23357250D-05, 0.38086643D-02, 0.23357098D-05,
1037 + 0.32198611D-01, 0.32198611D-01, 0.21824796D-05, 0.21824729D-05,
1038 + 0.21824729D-05, 0.21824796D-05, 0.26461108D-02, 0.21824678D-05,
1039 + 0.27462349D-01, 0.27462349D-01, 0.19814881D-05, 0.19814861D-05,
1040 + 0.19814861D-05, 0.19814881D-05, 0.18381056D-02, 0.19814790D-05/
1041 DATA ((CALCPI(I,K,24),I=1,8), K=10,18) /
1042 + 0.23369143D-01, 0.23369143D-01, 0.17576456D-05, 0.17576477D-05,
1043 + 0.17576477D-05, 0.17576456D-05, 0.12771408D-02, 0.17576385D-05,
1044 + 0.19840784D-01, 0.19840784D-01, 0.15318565D-05, 0.15318622D-05,
1045 + 0.15318622D-05, 0.15318565D-05, 0.88883501D-03, 0.15318510D-05,
1046 + 0.16807128D-01, 0.16807128D-01, 0.13168008D-05, 0.13168094D-05,
1047 + 0.13168094D-05, 0.13168008D-05, 0.62129100D-03, 0.13167965D-05,
1048 + 0.14205480D-01, 0.14205480D-01, 0.11194302D-05, 0.11194410D-05,
1049 + 0.11194410D-05, 0.11194302D-05, 0.43682179D-03, 0.11194269D-05,
1050 + 0.11979987D-01, 0.11979987D-01, 0.94337958D-06, 0.94339195D-06,
1051 + 0.94339195D-06, 0.94337958D-06, 0.30928510D-03, 0.94337707D-06,
1052 + 0.10081062D-01, 0.10081062D-01, 0.78959599D-06, 0.78960936D-06,
1053 + 0.78960936D-06, 0.78959599D-06, 0.22096292D-03, 0.78959406D-06,
1054 + 0.84648297D-02, 0.84648297D-02, 0.65732294D-06, 0.65733682D-06,
1055 + 0.65733682D-06, 0.65732294D-06, 0.15950703D-03, 0.65732145D-06,
1056 + 0.70926015D-02, 0.70926015D-02, 0.54494166D-06, 0.54495564D-06,
1057 + 0.54495564D-06, 0.54494166D-06, 0.11639689D-03, 0.54494052D-06,
1058 + 0.59303800D-02, 0.59303800D-02, 0.45036043D-06, 0.45037420D-06,
1059 + 0.45037420D-06, 0.45036043D-06, 0.85900136D-04, 0.45035956D-06/
1060 DATA ((CALCPI(I,K,24),I=1,8), K=19,20) /
1061 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
1062 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
1063 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
1064 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
1065 DATA ((CALCPI(I,K,25),I=1,8), K=1,9) /
1066 + 0.47514881D-01, 0.47514881D-01, 0.81282336D-13, 0.81281770D-13,
1067 + 0.81281770D-13, 0.81282336D-13, 0.11552384D-01,-0.44113689D-19,
1068 + 0.39558775D-01, 0.39558775D-01,-0.51902289D-06,-0.51897352D-06,
1069 + -0.51897352D-06,-0.51902289D-06, 0.74756386D-02,-0.51902295D-06,
1070 + 0.32873176D-01, 0.32873176D-01,-0.60102710D-06,-0.60094582D-06,
1071 + -0.60094582D-06,-0.60102710D-06, 0.48497831D-02,-0.60102715D-06,
1072 + 0.27249678D-01, 0.27249678D-01,-0.53440316D-06,-0.53430266D-06,
1073 + -0.53430266D-06,-0.53440316D-06, 0.31483696D-02,-0.53440320D-06,
1074 + 0.22530902D-01, 0.22530902D-01,-0.41313606D-06,-0.41302571D-06,
1075 + -0.41302571D-06,-0.41313606D-06, 0.20170920D-02,-0.41313609D-06,
1076 + 0.18583085D-01, 0.18583085D-01,-0.28292681D-06,-0.28281350D-06,
1077 + -0.28281350D-06,-0.28292681D-06, 0.12854537D-02,-0.28292683D-06,
1078 + 0.15288645D-01, 0.15288645D-01,-0.17073729D-06,-0.17062586D-06,
1079 + -0.17062586D-06,-0.17073729D-06, 0.82038244D-03,-0.17073731D-06,
1080 + 0.12546762D-01, 0.12546762D-01,-0.80730183D-07,-0.80623900D-07,
1081 + -0.80623900D-07,-0.80730183D-07, 0.52220283D-03,-0.80730194D-07,
1082 + 0.10270913D-01, 0.10270913D-01,-0.12645504D-07,-0.12546449D-07,
1083 + -0.12546449D-07,-0.12645504D-07, 0.33191852D-03,-0.12645512D-07/
1084 DATA ((CALCPI(I,K,25),I=1,8), K=10,18) /
1085 + 0.83869722D-02, 0.83869722D-02, 0.35232591D-07, 0.35323239D-07,
1086 + 0.35323239D-07, 0.35232591D-07, 0.21186226D-03, 0.35232585D-07,
1087 + 0.68316662D-02, 0.68316662D-02, 0.66871278D-07, 0.66953006D-07,
1088 + 0.66953006D-07, 0.66871278D-07, 0.13598500D-03, 0.66871273D-07,
1089 + 0.55511338D-02, 0.55511338D-02, 0.86095598D-07, 0.86168367D-07,
1090 + 0.86168367D-07, 0.86095598D-07, 0.87862421D-04, 0.86095595D-07,
1091 + 0.44996671D-02, 0.44996671D-02, 0.95920834D-07, 0.95984932D-07,
1092 + 0.95984932D-07, 0.95920834D-07, 0.57373579D-04, 0.95920831D-07,
1093 + 0.36385999D-02, 0.36385999D-02, 0.98939067D-07, 0.98994999D-07,
1094 + 0.98994999D-07, 0.98939067D-07, 0.37963821D-04, 0.98939065D-07,
1095 + 0.29353270D-02, 0.29353270D-02, 0.97257041D-07, 0.97305442D-07,
1096 + 0.97305442D-07, 0.97257041D-07, 0.25480951D-04, 0.97257039D-07,
1097 + 0.23624426D-02, 0.23624426D-02, 0.92456560D-07, 0.92498132D-07,
1098 + 0.92498132D-07, 0.92456560D-07, 0.17376552D-04, 0.92456559D-07,
1099 + 0.18969857D-02, 0.18969857D-02, 0.85730273D-07, 0.85765738D-07,
1100 + 0.85765738D-07, 0.85730273D-07, 0.12050384D-04, 0.85730272D-07,
1101 + 0.15197821D-02, 0.15197821D-02, 0.77968995D-07, 0.77999062D-07,
1102 + 0.77999062D-07, 0.77968995D-07, 0.84914737D-05, 0.77968995D-07/
1103 DATA ((CALCPI(I,K,25),I=1,8), K=19,20) /
1104 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
1105 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
1106 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
1107 + 0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
1108C----------------------------------------------------------------------
1109 DATA ISTART/0/
1110 SAVE ISTART,OWLAM2,Q02PI
1111C
1112 IF (ISTART.EQ.0) THEN
1113 ISTART=1
1114 DO 11 K=1,25
1115 DO 11 I=1,20
1116 DO 11 M=1,8
1117 11 CALCPIO(M,I,K) = CALCPI(M,I,K)
1118 OWLAM=PARPI(1)
1119 OWLAM2=OWLAM**2
1120 Q02PI=PARPI(39)
1121 Q2MAX=PARPI(40)
1122 ENDIF
1123C
1124C the conventions are : q(1)=x*u, q(2)=x*d, q(3)=x*str, q(4)=x*usea,
1125C q(5)=x*dsea, q(6)=x*charm, q(7)=x*gluon
1126C
1127 X = DX
1128 Q = DQ
1129 Q2 = Q*Q
1130 IDQ2=2
1131 SB=0.
1132 IF(Q2-Q02PI) 1,1,2
1133 2 IF(IDQ2-1) 1,1,3
1134 3 SB= LOG( LOG( MAX(Q02PI,Q2)/OWLAM2)/ LOG(Q02PI/OWLAM2))
1135 1 CALL AURPI2(1,0,X,SB,XPDF(1))
1136 CALL AURPI2(2,0,X,SB,XPDF(2))
1137 CALL AURPI2(3,0,X,SB,XPDF(3))
1138 CALL AURPI2(4,0,X,SB,XPDF(4))
1139 CALL AURPI2(5,0,X,SB,XPDF(5))
1140 CALL AURPI2(8,0,X,SB,XPDF(6))
1141 CALL AURPI2(7,0,X,SB,XPDF(7))
1142C
1143 DUPV=XPDF(1) - XPDF(4)
1144 DDNV=XPDF(2) - XPDF(5)
1145 DUSEA=XPDF(4)
1146 DDSEA=XPDF(5)
1147 DSTR=XPDF(3)
1148 DCHM=XPDF(6)
1149 DGL =XPDF(7)
1150C
1151 RETURN
1152 END