]> git.uio.no Git - u/mrichter/AliRoot.git/blob - PDF/spdf/abfkw2.F
negative indexes allowed
[u/mrichter/AliRoot.git] / PDF / spdf / abfkw2.F
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/
21 C----------------------------------------------------------------------
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/
1108 C----------------------------------------------------------------------
1109        DATA ISTART/0/
1110        SAVE ISTART,OWLAM2,Q02PI
1111 C
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
1123 C
1124 C the conventions are : q(1)=x*u, q(2)=x*d, q(3)=x*str, q(4)=x*usea,
1125 C                       q(5)=x*dsea, q(6)=x*charm, q(7)=x*gluon
1126 C
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))
1142 C
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)
1150 C
1151       RETURN
1152       END