]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PDF/spdf/abfkw3.F
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Removing warnings (icc), adding more detailed description
[u/mrichter/AliRoot.git] / PDF / spdf / abfkw3.F
CommitLineData
21886bb6 1*
2* $Id$
3*
4* $Log$
5* Revision 1.1.1.2 1996/10/30 08:27:28 cernlib
6* Version 7.04
7*
8* Revision 1.1.1.1 1996/04/12 15:28:54 plothow
9* Version 7.01
10*
11*
12#include "pdf/pilot.h"
13 SUBROUTINE ABFKW3(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/W5051I6/CALCPIO
19 REAL XPDF(7)
20 DATA ZEROD/0.D0/
21C----------------------------------------------------------------------
22 DATA (PARPI(I),I=1,40) /
23 + 0.28100002D+00, 0.38845206D+00, 0.69999999D+00, 0.00000000D+00,
24 + 0.00000000D+00, 0.38845206D+00, 0.69999999D+00, 0.00000000D+00,
25 + 0.56699985D+00, 0.53000002D+01, 0.19346480D+01, 0.00000000D+00,
26 + 0.25959997D+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.12821481D+00, 0.12821481D+00, 0.94000201D-01, 0.94000201D-01,
35 + 0.94000201D-01, 0.94000201D-01, 0.19296297D+01, 0.34694470D-17,
36 + 0.16706269D+00, 0.16706269D+00, 0.13168872D+00, 0.13167343D+00,
37 + 0.13167343D+00, 0.13168872D+00, 0.29529459D+01, 0.37834649D-01,
38 + 0.21189126D+00, 0.21189126D+00, 0.17534196D+00, 0.17531252D+00,
39 + 0.17531252D+00, 0.17534196D+00, 0.41492404D+01, 0.81641090D-01,
40 + 0.26300280D+00, 0.26300280D+00, 0.22526172D+00, 0.22521926D+00,
41 + 0.22521926D+00, 0.22526172D+00, 0.55198687D+01, 0.13172154D+00,
42 + 0.32065997D+00, 0.32065997D+00, 0.28171031D+00, 0.28165595D+00,
43 + 0.28165595D+00, 0.28171031D+00, 0.70633861D+01, 0.18833868D+00,
44 + 0.38510573D+00, 0.38510573D+00, 0.34493036D+00, 0.34486519D+00,
45 + 0.34486519D+00, 0.34493036D+00, 0.87758667D+01, 0.25173554D+00,
46 + 0.45655511D+00, 0.45655511D+00, 0.41513665D+00, 0.41506174D+00,
47 + 0.41506174D+00, 0.41513665D+00, 0.10651241D+02, 0.32212729D+00,
48 + 0.53518537D+00, 0.53518537D+00, 0.49250627D+00, 0.49242266D+00,
49 + 0.49242266D+00, 0.49250627D+00, 0.12681165D+02, 0.39969144D+00,
50 + 0.62112973D+00, 0.62112973D+00, 0.57717230D+00, 0.57708100D+00,
51 + 0.57708100D+00, 0.57717230D+00, 0.14855182D+02, 0.48456148D+00/
52 DATA ((CALCPI(I,K,1 ),I=1,8), K=10,18) /
53 + 0.71446944D+00, 0.71446944D+00, 0.66921594D+00, 0.66911794D+00,
54 + 0.66911794D+00, 0.66921594D+00, 0.17160841D+02, 0.57681905D+00,
55 + 0.81522779D+00, 0.81522779D+00, 0.76866048D+00, 0.76855674D+00,
56 + 0.76855674D+00, 0.76866048D+00, 0.19583877D+02, 0.67648786D+00,
57 + 0.92336459D+00, 0.92336459D+00, 0.87546584D+00, 0.87535726D+00,
58 + 0.87535726D+00, 0.87546584D+00, 0.22108423D+02, 0.78352823D+00,
59 + 0.10387719D+01, 0.10387719D+01, 0.98952429D+00, 0.98941176D+00,
60 + 0.98941176D+00, 0.98952429D+00, 0.24717248D+02, 0.89783290D+00,
61 + 0.11612713D+01, 0.11612713D+01, 0.11106575D+01, 0.11105419D+01,
62 + 0.11105419D+01, 0.11106575D+01, 0.27392024D+02, 0.10192239D+01,
63 + 0.12906118D+01, 0.12906118D+01, 0.12386150D+01, 0.12384970D+01,
64 + 0.12384970D+01, 0.12386150D+01, 0.30113602D+02, 0.11474512D+01,
65 + 0.14264699D+01, 0.14264699D+01, 0.13730736D+01, 0.13729542D+01,
66 + 0.13729542D+01, 0.13730736D+01, 0.32862311D+02, 0.12821921D+01,
67 + 0.15684509D+01, 0.15684509D+01, 0.15136392D+01, 0.15135189D+01,
68 + 0.15135189D+01, 0.15136392D+01, 0.35618247D+02, 0.14230528D+01,
69 + 0.17160913D+01, 0.17160913D+01, 0.16598486D+01, 0.16597283D+01,
70 + 0.16597283D+01, 0.16598486D+01, 0.38361571D+02, 0.15695705D+01/
71 DATA ((CALCPI(I,K,1 ),I=1,8), K=19,20) /
72 + 0.18688627D+01, 0.18688627D+01, 0.18111744D+01, 0.18110546D+01,
73 + 0.18110546D+01, 0.18111744D+01, 0.41072781D+02, 0.17212181D+01,
74 + 0.20261771D+01, 0.20261771D+01, 0.19670292D+01, 0.19669105D+01,
75 + 0.19669105D+01, 0.19670292D+01, 0.43732985D+02, 0.18774087D+01/
76 DATA ((CALCPI(I,K,2 ),I=1,8), K=1,9) /
77 + 0.13825770D+00, 0.13825770D+00, 0.93502573D-01, 0.93502573D-01,
78 + 0.93502573D-01, 0.93502573D-01, 0.19246194D+01, 0.17347235D-17,
79 + 0.17173978D+00, 0.17173978D+00, 0.12552391D+00, 0.12550690D+00,
80 + 0.12550690D+00, 0.12552391D+00, 0.27841085D+01, 0.32253347D-01,
81 + 0.20927357D+00, 0.20927357D+00, 0.16158395D+00, 0.16155147D+00,
82 + 0.16155147D+00, 0.16158395D+00, 0.37516700D+01, 0.68555971D-01,
83 + 0.25097815D+00, 0.25097815D+00, 0.20180138D+00, 0.20175493D+00,
84 + 0.20175493D+00, 0.20180138D+00, 0.48230626D+01, 0.10902704D+00,
85 + 0.29693944D+00, 0.29693944D+00, 0.24626185D+00, 0.24620286D+00,
86 + 0.24620286D+00, 0.24626185D+00, 0.59928410D+01, 0.15375265D+00,
87 + 0.34721071D+00, 0.34721071D+00, 0.29501848D+00, 0.29494834D+00,
88 + 0.29494834D+00, 0.29501848D+00, 0.72538108D+01, 0.20278636D+00,
89 + 0.40181628D+00, 0.40181628D+00, 0.34809551D+00, 0.34801554D+00,
90 + 0.34801554D+00, 0.34809551D+00, 0.85975660D+01, 0.25615288D+00,
91 + 0.46074351D+00, 0.46074351D+00, 0.40548033D+00, 0.40539180D+00,
92 + 0.40539180D+00, 0.40548033D+00, 0.10014442D+02, 0.31384006D+00,
93 + 0.52394162D+00, 0.52394162D+00, 0.46712230D+00, 0.46702642D+00,
94 + 0.46702642D+00, 0.46712230D+00, 0.11493727D+02, 0.37579770D+00/
95 DATA ((CALCPI(I,K,2 ),I=1,8), K=10,18) /
96 + 0.59131859D+00, 0.59131859D+00, 0.53292963D+00, 0.53282756D+00,
97 + 0.53282756D+00, 0.53292963D+00, 0.13023825D+02, 0.44193449D+00,
98 + 0.66273983D+00, 0.66273983D+00, 0.60276805D+00, 0.60266088D+00,
99 + 0.60266088D+00, 0.60276805D+00, 0.14592450D+02, 0.51211659D+00,
100 + 0.73802729D+00, 0.73802729D+00, 0.67645999D+00, 0.67634876D+00,
101 + 0.67634876D+00, 0.67645999D+00, 0.16186834D+02, 0.58616688D+00,
102 + 0.81695970D+00, 0.81695970D+00, 0.75378471D+00, 0.75367041D+00,
103 + 0.75367041D+00, 0.75378471D+00, 0.17793936D+02, 0.66386502D+00,
104 + 0.89927349D+00, 0.89927349D+00, 0.83447932D+00, 0.83436285D+00,
105 + 0.83436285D+00, 0.83447932D+00, 0.19400649D+02, 0.74494851D+00,
106 + 0.98466465D+00, 0.98466465D+00, 0.91824059D+00, 0.91812282D+00,
107 + 0.91812282D+00, 0.91824059D+00, 0.20994011D+02, 0.82911454D+00,
108 + 0.10727915D+01, 0.10727915D+01, 0.10047277D+01, 0.10046095D+01,
109 + 0.10046095D+01, 0.10047277D+01, 0.22561398D+02, 0.91602264D+00,
110 + 0.11632781D+01, 0.11632781D+01, 0.10935658D+01, 0.10934478D+01,
111 + 0.10934478D+01, 0.10935658D+01, 0.24090706D+02, 0.10052982D+01,
112 + 0.12557185D+01, 0.12557185D+01, 0.11843498D+01, 0.11842328D+01,
113 + 0.11842328D+01, 0.11843498D+01, 0.25570518D+02, 0.10965366D+01/
114 DATA ((CALCPI(I,K,2 ),I=1,8), K=19,20) /
115 + 0.13496811D+01, 0.13496811D+01, 0.12766495D+01, 0.12765341D+01,
116 + 0.12765341D+01, 0.12766495D+01, 0.26990245D+02, 0.11893079D+01,
117 + 0.14447144D+01, 0.14447144D+01, 0.13700146D+01, 0.13699015D+01,
118 + 0.13699015D+01, 0.13700146D+01, 0.28340247D+02, 0.12831619D+01/
119 DATA ((CALCPI(I,K,3 ),I=1,8), K=1,9) /
120 + 0.15101546D+00, 0.15101546D+00, 0.92513729D-01, 0.92513729D-01,
121 + 0.92513729D-01, 0.92513729D-01, 0.19146227D+01, 0.52041704D-17,
122 + 0.17916929D+00, 0.17916929D+00, 0.11886813D+00, 0.11885069D+00,
123 + 0.11885069D+00, 0.11886813D+00, 0.26079692D+01, 0.26723688D-01,
124 + 0.20995742D+00, 0.20995742D+00, 0.14785133D+00, 0.14781832D+00,
125 + 0.14781832D+00, 0.14785133D+00, 0.33618996D+01, 0.56091119D-01,
126 + 0.24336338D+00, 0.24336338D+00, 0.17944661D+00, 0.17939979D+00,
127 + 0.17939979D+00, 0.17944661D+00, 0.41695510D+01, 0.88086045D-01,
128 + 0.27935265D+00, 0.27935265D+00, 0.21361940D+00, 0.21356042D+00,
129 + 0.21356042D+00, 0.21361940D+00, 0.50240311D+01, 0.12267432D+00,
130 + 0.31785866D+00, 0.31785866D+00, 0.25030321D+00, 0.25023363D+00,
131 + 0.25023363D+00, 0.25030321D+00, 0.59174672D+01, 0.15978989D+00,
132 + 0.35879420D+00, 0.35879420D+00, 0.28941106D+00, 0.28933234D+00,
133 + 0.28933234D+00, 0.28941106D+00, 0.68415748D+01, 0.19934620D+00,
134 + 0.40204568D+00, 0.40204568D+00, 0.33082974D+00, 0.33074327D+00,
135 + 0.33074327D+00, 0.33082974D+00, 0.77876678D+01, 0.24123046D+00,
136 + 0.44747474D+00, 0.44747474D+00, 0.37442143D+00, 0.37432850D+00,
137 + 0.37432850D+00, 0.37442143D+00, 0.87468404D+01, 0.28530523D+00/
138 DATA ((CALCPI(I,K,3 ),I=1,8), K=10,18) /
139 + 0.49491821D+00, 0.49491821D+00, 0.42002365D+00, 0.41992544D+00,
140 + 0.41992544D+00, 0.42002365D+00, 0.97101188D+01, 0.33140839D+00,
141 + 0.54418909D+00, 0.54418909D+00, 0.46745020D+00, 0.46734786D+00,
142 + 0.46734786D+00, 0.46745020D+00, 0.10668609D+02, 0.37935413D+00,
143 + 0.59507796D+00, 0.59507796D+00, 0.51649267D+00, 0.51638722D+00,
144 + 0.51638722D+00, 0.51649267D+00, 0.11613643D+02, 0.42893431D+00,
145 + 0.64735488D+00, 0.64735488D+00, 0.56692221D+00, 0.56681463D+00,
146 + 0.56681463D+00, 0.56692221D+00, 0.12536911D+02, 0.47992040D+00,
147 + 0.70077172D+00, 0.70077172D+00, 0.61849196D+00, 0.61838313D+00,
148 + 0.61838313D+00, 0.61849196D+00, 0.13430581D+02, 0.53206577D+00,
149 + 0.75506487D+00, 0.75506487D+00, 0.67093970D+00, 0.67083044D+00,
150 + 0.67083044D+00, 0.67093970D+00, 0.14287405D+02, 0.58510841D+00,
151 + 0.80995831D+00, 0.80995831D+00, 0.72399093D+00, 0.72388200D+00,
152 + 0.72388200D+00, 0.72399093D+00, 0.15100808D+02, 0.63877400D+00,
153 + 0.86516694D+00, 0.86516694D+00, 0.77736219D+00, 0.77725427D+00,
154 + 0.77725427D+00, 0.77736219D+00, 0.15864957D+02, 0.69277918D+00,
155 + 0.92040002D+00, 0.92040002D+00, 0.83076449D+00, 0.83065822D+00,
156 + 0.83065822D+00, 0.83076449D+00, 0.16574813D+02, 0.74683507D+00/
157 DATA ((CALCPI(I,K,3 ),I=1,8), K=19,20) /
158 + 0.97536478D+00, 0.97536478D+00, 0.88390694D+00, 0.88380289D+00,
159 + 0.88380289D+00, 0.88390694D+00, 0.17226163D+02, 0.80065078D+00,
160 + 0.10297700D+01, 0.10297700D+01, 0.93650026D+00, 0.93639894D+00,
161 + 0.93639894D+00, 0.93650026D+00, 0.17815637D+02, 0.85393701D+00/
162 DATA ((CALCPI(I,K,4 ),I=1,8), K=1,9) /
163 + 0.16692415D+00, 0.16692415D+00, 0.90561489D-01, 0.90561489D-01,
164 + 0.90561489D-01, 0.90561489D-01, 0.18947253D+01, 0.17347235D-17,
165 + 0.19026586D+00, 0.19026586D+00, 0.11175541D+00, 0.11173886D+00,
166 + 0.11173886D+00, 0.11175541D+00, 0.24332533D+01, 0.21778920D-01,
167 + 0.21511552D+00, 0.21511552D+00, 0.13446678D+00, 0.13443561D+00,
168 + 0.13443561D+00, 0.13446678D+00, 0.29966805D+01, 0.45094164D-01,
169 + 0.24137560D+00, 0.24137560D+00, 0.15859806D+00, 0.15855405D+00,
170 + 0.15855405D+00, 0.15859806D+00, 0.35779189D+01, 0.69848543D-01,
171 + 0.26894901D+00, 0.26894901D+00, 0.18405246D+00, 0.18399726D+00,
172 + 0.18399726D+00, 0.18405246D+00, 0.41706395D+01, 0.95945524D-01,
173 + 0.29771728D+00, 0.29771728D+00, 0.21071199D+00, 0.21064713D+00,
174 + 0.21064713D+00, 0.21071199D+00, 0.47682891D+01, 0.12326738D+00,
175 + 0.32755295D+00, 0.32755295D+00, 0.23844989D+00, 0.23837677D+00,
176 + 0.23837677D+00, 0.23844989D+00, 0.53644818D+01, 0.15168758D+00,
177 + 0.35831595D+00, 0.35831595D+00, 0.26712702D+00, 0.26704695D+00,
178 + 0.26704695D+00, 0.26712702D+00, 0.59530343D+01, 0.18106719D+00,
179 + 0.38985562D+00, 0.38985562D+00, 0.29659383D+00, 0.29650800D+00,
180 + 0.29650800D+00, 0.29659383D+00, 0.65280608D+01, 0.21125682D+00/
181 DATA ((CALCPI(I,K,4 ),I=1,8), K=10,18) /
182 + 0.42201175D+00, 0.42201175D+00, 0.32669141D+00, 0.32660093D+00,
183 + 0.32660093D+00, 0.32669141D+00, 0.70840600D+01, 0.24209768D+00,
184 + 0.45461604D+00, 0.45461604D+00, 0.35725295D+00, 0.35715883D+00,
185 + 0.35715883D+00, 0.35725295D+00, 0.76159841D+01, 0.27342304D+00,
186 + 0.48749384D+00, 0.48749384D+00, 0.38810547D+00, 0.38800864D+00,
187 + 0.38800864D+00, 0.38810547D+00, 0.81192978D+01, 0.30505992D+00,
188 + 0.52046602D+00, 0.52046602D+00, 0.41907170D+00, 0.41897300D+00,
189 + 0.41897300D+00, 0.41907170D+00, 0.85900243D+01, 0.33683103D+00,
190 + 0.55335107D+00, 0.55335107D+00, 0.44997207D+00, 0.44987228D+00,
191 + 0.44987228D+00, 0.44997207D+00, 0.90247770D+01, 0.36855671D+00,
192 + 0.58596713D+00, 0.58596713D+00, 0.48062687D+00, 0.48052671D+00,
193 + 0.48052671D+00, 0.48062687D+00, 0.94207781D+01, 0.40005710D+00,
194 + 0.61813422D+00, 0.61813422D+00, 0.51085839D+00, 0.51075848D+00,
195 + 0.51075848D+00, 0.51085839D+00, 0.97758647D+01, 0.43115428D+00,
196 + 0.64967637D+00, 0.64967637D+00, 0.54049301D+00, 0.54039394D+00,
197 + 0.54039394D+00, 0.54049301D+00, 0.10088482D+02, 0.46167434D+00,
198 + 0.68042364D+00, 0.68042364D+00, 0.56936325D+00, 0.56926555D+00,
199 + 0.56926555D+00, 0.56936325D+00, 0.10357667D+02, 0.49144950D+00/
200 DATA ((CALCPI(I,K,4 ),I=1,8), K=19,20) /
201 + 0.71021408D+00, 0.71021408D+00, 0.59730972D+00, 0.59721388D+00,
202 + 0.59721388D+00, 0.59730972D+00, 0.10583019D+02, 0.52031996D+00,
203 + 0.73889551D+00, 0.73889551D+00, 0.62418286D+00, 0.62408930D+00,
204 + 0.62408930D+00, 0.62418286D+00, 0.10764666D+02, 0.54813569D+00/
205 DATA ((CALCPI(I,K,5 ),I=1,8), K=1,9) /
206 + 0.18615009D+00, 0.18615009D+00, 0.86757238D-01, 0.86757238D-01,
207 + 0.86757238D-01, 0.86757238D-01, 0.18553132D+01, 0.34694470D-17,
208 + 0.20495309D+00, 0.20495309D+00, 0.10312192D+00, 0.10310608D+00,
209 + 0.10310608D+00, 0.10312192D+00, 0.22468147D+01, 0.17268707D-01,
210 + 0.22435737D+00, 0.22435737D+00, 0.12012473D+00, 0.12009485D+00,
211 + 0.12009485D+00, 0.12012473D+00, 0.26379097D+01, 0.35194729D-01,
212 + 0.24424792D+00, 0.24424792D+00, 0.13765097D+00, 0.13760865D+00,
213 + 0.13760865D+00, 0.13765097D+00, 0.30233411D+01, 0.53663379D-01,
214 + 0.26452068D+00, 0.26452068D+00, 0.15559727D+00, 0.15554401D+00,
215 + 0.15554401D+00, 0.15559727D+00, 0.33988506D+01, 0.72571273D-01,
216 + 0.28506247D+00, 0.28506247D+00, 0.17385140D+00, 0.17378858D+00,
217 + 0.17378858D+00, 0.17385140D+00, 0.37604487D+01, 0.91806124D-01,
218 + 0.30575952D+00, 0.30575952D+00, 0.19230084D+00, 0.19222972D+00,
219 + 0.19222972D+00, 0.19230084D+00, 0.41045496D+01, 0.11125528D+00,
220 + 0.32649595D+00, 0.32649595D+00, 0.21083122D+00, 0.21075297D+00,
221 + 0.21075297D+00, 0.21083122D+00, 0.44280073D+01, 0.13080421D+00,
222 + 0.34715510D+00, 0.34715510D+00, 0.22932760D+00, 0.22924328D+00,
223 + 0.22924328D+00, 0.22932760D+00, 0.47281239D+01, 0.15033770D+00/
224 DATA ((CALCPI(I,K,5 ),I=1,8), K=10,18) /
225 + 0.36762059D+00, 0.36762059D+00, 0.24767556D+00, 0.24758618D+00,
226 + 0.24758618D+00, 0.24767556D+00, 0.50026696D+01, 0.16974106D+00,
227 + 0.38777749D+00, 0.38777749D+00, 0.26576234D+00, 0.26566879D+00,
228 + 0.26566879D+00, 0.26576234D+00, 0.52498855D+01, 0.18890110D+00,
229 + 0.40751352D+00, 0.40751352D+00, 0.28347797D+00, 0.28338110D+00,
230 + 0.28338110D+00, 0.28347797D+00, 0.54684785D+01, 0.20770743D+00,
231 + 0.42672029D+00, 0.42672029D+00, 0.30071651D+00, 0.30061709D+00,
232 + 0.30061709D+00, 0.30071651D+00, 0.56576093D+01, 0.22605359D+00,
233 + 0.44529442D+00, 0.44529442D+00, 0.31737719D+00, 0.31727594D+00,
234 + 0.31727594D+00, 0.31737719D+00, 0.58168719D+01, 0.24383821D+00,
235 + 0.46313871D+00, 0.46313871D+00, 0.33336551D+00, 0.33326309D+00,
236 + 0.33326309D+00, 0.33336551D+00, 0.59462671D+01, 0.26096612D+00,
237 + 0.48016314D+00, 0.48016314D+00, 0.34859422D+00, 0.34849124D+00,
238 + 0.34849124D+00, 0.34859422D+00, 0.60461703D+01, 0.27734934D+00,
239 + 0.49628577D+00, 0.49628577D+00, 0.36298420D+00, 0.36288122D+00,
240 + 0.36288122D+00, 0.36298420D+00, 0.61172950D+01, 0.29290795D+00,
241 + 0.51143347D+00, 0.51143347D+00, 0.37646518D+00, 0.37636272D+00,
242 + 0.37636272D+00, 0.37646518D+00, 0.61606536D+01, 0.30757082D+00/
243 DATA ((CALCPI(I,K,5 ),I=1,8), K=19,20) /
244 + 0.52554254D+00, 0.52554254D+00, 0.38897630D+00, 0.38887487D+00,
245 + 0.38887487D+00, 0.38897630D+00, 0.61775152D+01, 0.32127618D+00,
246 + 0.53855912D+00, 0.53855912D+00, 0.40046654D+00, 0.40036658D+00,
247 + 0.40036658D+00, 0.40046654D+00, 0.61693635D+01, 0.33397202D+00/
248 DATA ((CALCPI(I,K,6 ),I=1,8), K=1,9) /
249 + 0.20815692D+00, 0.20815692D+00, 0.79537344D-01, 0.79537344D-01,
250 + 0.79537344D-01, 0.79537344D-01, 0.17780108D+01, 0.00000000D+00,
251 + 0.22237127D+00, 0.22237127D+00, 0.91233284D-01, 0.91217009D-01,
252 + 0.91217009D-01, 0.91233284D-01, 0.20285821D+01, 0.13016954D-01,
253 + 0.23652036D+00, 0.23652036D+00, 0.10294598D+00, 0.10291520D+00,
254 + 0.10291520D+00, 0.10294598D+00, 0.22645799D+01, 0.26060969D-01,
255 + 0.25051540D+00, 0.25051540D+00, 0.11458697D+00, 0.11454329D+00,
256 + 0.11454329D+00, 0.11458697D+00, 0.24835026D+01, 0.39043296D-01,
257 + 0.26428040D+00, 0.26428040D+00, 0.12608111D+00, 0.12602602D+00,
258 + 0.12602602D+00, 0.12608111D+00, 0.26836091D+01, 0.51888460D-01,
259 + 0.27773968D+00, 0.27773968D+00, 0.13735389D+00, 0.13728875D+00,
260 + 0.13728875D+00, 0.13735389D+00, 0.28635744D+01, 0.64521505D-01,
261 + 0.29082109D+00, 0.29082109D+00, 0.14833460D+00, 0.14826068D+00,
262 + 0.14826068D+00, 0.14833460D+00, 0.30224311D+01, 0.76871221D-01,
263 + 0.30345607D+00, 0.30345607D+00, 0.15895645D+00, 0.15887492D+00,
264 + 0.15887492D+00, 0.15895645D+00, 0.31595930D+01, 0.88870204D-01,
265 + 0.31558008D+00, 0.31558008D+00, 0.16915683D+00, 0.16906879D+00,
266 + 0.16906879D+00, 0.16915683D+00, 0.32748177D+01, 0.10045518D+00/
267 DATA ((CALCPI(I,K,6 ),I=1,8), K=10,18) /
268 + 0.32713323D+00, 0.32713323D+00, 0.17887803D+00, 0.17878449D+00,
269 + 0.17878449D+00, 0.17887803D+00, 0.33681856D+01, 0.11156764D+00,
270 + 0.33806080D+00, 0.33806080D+00, 0.18806760D+00, 0.18796950D+00,
271 + 0.18796950D+00, 0.18806760D+00, 0.34400699D+01, 0.12215431D+00,
272 + 0.34831370D+00, 0.34831370D+00, 0.19667883D+00, 0.19657707D+00,
273 + 0.19657707D+00, 0.19667883D+00, 0.34911049D+01, 0.13216757D+00,
274 + 0.35784889D+00, 0.35784889D+00, 0.20467114D+00, 0.20456653D+00,
275 + 0.20456653D+00, 0.20467114D+00, 0.35221527D+01, 0.14156582D+00,
276 + 0.36662969D+00, 0.36662969D+00, 0.21201031D+00, 0.21190364D+00,
277 + 0.21190364D+00, 0.21201031D+00, 0.35342677D+01, 0.15031382D+00,
278 + 0.37462601D+00, 0.37462601D+00, 0.21866872D+00, 0.21856070D+00,
279 + 0.21856070D+00, 0.21866872D+00, 0.35286620D+01, 0.15838283D+00,
280 + 0.38181449D+00, 0.38181449D+00, 0.22462542D+00, 0.22451674D+00,
281 + 0.22451674D+00, 0.22462542D+00, 0.35066703D+01, 0.16575074D+00,
282 + 0.38817854D+00, 0.38817854D+00, 0.22986612D+00, 0.22975741D+00,
283 + 0.22975741D+00, 0.22986612D+00, 0.34697162D+01, 0.17240205D+00,
284 + 0.39370821D+00, 0.39370821D+00, 0.23438309D+00, 0.23427495D+00,
285 + 0.23427495D+00, 0.23438309D+00, 0.34192800D+01, 0.17832780D+00/
286 DATA ((CALCPI(I,K,6 ),I=1,8), K=19,20) /
287 + 0.39840004D+00, 0.39840004D+00, 0.23817494D+00, 0.23806792D+00,
288 + 0.23806792D+00, 0.23817494D+00, 0.33568697D+01, 0.18352533D+00,
289 + 0.40225680D+00, 0.40225680D+00, 0.24124632D+00, 0.24114093D+00,
290 + 0.24114093D+00, 0.24124632D+00, 0.32839933D+01, 0.18799799D+00/
291 DATA ((CALCPI(I,K,7 ),I=1,8), K=1,9) /
292 + 0.23100293D+00, 0.23100293D+00, 0.66556988D-01, 0.66556988D-01,
293 + 0.66556988D-01, 0.66556988D-01, 0.16294238D+01, 0.17347235D-17,
294 + 0.24034847D+00, 0.24034847D+00, 0.73819110D-01, 0.73806365D-01,
295 + 0.73806365D-01, 0.73819110D-01, 0.17504465D+01, 0.89976789D-02,
296 + 0.24921003D+00, 0.24921003D+00, 0.80745149D-01, 0.80721143D-01,
297 + 0.80721143D-01, 0.80745149D-01, 0.18533173D+01, 0.17645027D-01,
298 + 0.25754928D+00, 0.25754928D+00, 0.87295945D-01, 0.87262040D-01,
299 + 0.87262040D-01, 0.87295945D-01, 0.19380990D+01, 0.25902814D-01,
300 + 0.26533543D+00, 0.26533543D+00, 0.93440694D-01, 0.93398162D-01,
301 + 0.93398162D-01, 0.93440694D-01, 0.20051235D+01, 0.33739992D-01,
302 + 0.27254263D+00, 0.27254263D+00, 0.99153735D-01, 0.99103759D-01,
303 + 0.99103759D-01, 0.99153735D-01, 0.20550180D+01, 0.41130553D-01,
304 + 0.27914856D+00, 0.27914856D+00, 0.10441325D+00, 0.10435693D+00,
305 + 0.10435693D+00, 0.10441325D+00, 0.20885678D+01, 0.48052202D-01,
306 + 0.28513538D+00, 0.28513538D+00, 0.10920207D+00, 0.10914044D+00,
307 + 0.10914044D+00, 0.10920207D+00, 0.21067258D+01, 0.54487217D-01,
308 + 0.29048947D+00, 0.29048947D+00, 0.11350739D+00, 0.11344140D+00,
309 + 0.11344140D+00, 0.11350739D+00, 0.21105715D+01, 0.60422105D-01/
310 DATA ((CALCPI(I,K,7 ),I=1,8), K=10,18) /
311 + 0.29520153D+00, 0.29520153D+00, 0.11732076D+00, 0.11725131D+00,
312 + 0.11725131D+00, 0.11732076D+00, 0.21012807D+01, 0.65847705D-01,
313 + 0.29926662D+00, 0.29926662D+00, 0.12063808D+00, 0.12056601D+00,
314 + 0.12056601D+00, 0.12063808D+00, 0.20800963D+01, 0.70759107D-01,
315 + 0.30268404D+00, 0.30268404D+00, 0.12345944D+00, 0.12338552D+00,
316 + 0.12338552D+00, 0.12345944D+00, 0.20483002D+01, 0.75155547D-01,
317 + 0.30545725D+00, 0.30545725D+00, 0.12578895D+00, 0.12571390D+00,
318 + 0.12571390D+00, 0.12578895D+00, 0.20071875D+01, 0.79040253D-01,
319 + 0.30759371D+00, 0.30759371D+00, 0.12763456D+00, 0.12755904D+00,
320 + 0.12755904D+00, 0.12763456D+00, 0.19580434D+01, 0.82420247D-01,
321 + 0.30910463D+00, 0.30910463D+00, 0.12900776D+00, 0.12893238D+00,
322 + 0.12893238D+00, 0.12900776D+00, 0.19021229D+01, 0.85306089D-01,
323 + 0.31000475D+00, 0.31000475D+00, 0.12992330D+00, 0.12984862D+00,
324 + 0.12984862D+00, 0.12992330D+00, 0.18406327D+01, 0.87711587D-01,
325 + 0.31031198D+00, 0.31031198D+00, 0.13039884D+00, 0.13032538D+00,
326 + 0.13032538D+00, 0.13039884D+00, 0.17747172D+01, 0.89653466D-01,
327 + 0.31004713D+00, 0.31004713D+00, 0.13045458D+00, 0.13038281D+00,
328 + 0.13038281D+00, 0.13045458D+00, 0.17054465D+01, 0.91151017D-01/
329 DATA ((CALCPI(I,K,7 ),I=1,8), K=19,20) /
330 + 0.30923347D+00, 0.30923347D+00, 0.13011289D+00, 0.13004323D+00,
331 + 0.13004323D+00, 0.13011289D+00, 0.16338078D+01, 0.92225716D-01,
332 + 0.30789643D+00, 0.30789643D+00, 0.12939792D+00, 0.12933075D+00,
333 + 0.12933075D+00, 0.12939792D+00, 0.15606997D+01, 0.92900840D-01/
334 DATA ((CALCPI(I,K,8 ),I=1,8), K=1,9) /
335 + 0.24434354D+00, 0.24434354D+00, 0.54065089D-01, 0.54065089D-01,
336 + 0.54065089D-01, 0.54065089D-01, 0.14716872D+01, 0.34694470D-17,
337 + 0.25029537D+00, 0.25029537D+00, 0.58741459D-01, 0.58731575D-01,
338 + 0.58731575D-01, 0.58741459D-01, 0.15201810D+01, 0.65477534D-02,
339 + 0.25563519D+00, 0.25563519D+00, 0.63002654D-01, 0.62984361D-01,
340 + 0.62984361D-01, 0.63002654D-01, 0.15531386D+01, 0.12642397D-01,
341 + 0.26035998D+00, 0.26035998D+00, 0.66842795D-01, 0.66817425D-01,
342 + 0.66817425D-01, 0.66842795D-01, 0.15717256D+01, 0.18278793D-01,
343 + 0.26446844D+00, 0.26446844D+00, 0.70258812D-01, 0.70227586D-01,
344 + 0.70227586D-01, 0.70258812D-01, 0.15770653D+01, 0.23454263D-01,
345 + 0.26796426D+00, 0.26796426D+00, 0.73252646D-01, 0.73216680D-01,
346 + 0.73216680D-01, 0.73252646D-01, 0.15704171D+01, 0.28171043D-01,
347 + 0.27085297D+00, 0.27085297D+00, 0.75828350D-01, 0.75788660D-01,
348 + 0.75788660D-01, 0.75828350D-01, 0.15530551D+01, 0.32433334D-01,
349 + 0.27314313D+00, 0.27314313D+00, 0.77993030D-01, 0.77950544D-01,
350 + 0.77950544D-01, 0.77993030D-01, 0.15262694D+01, 0.36248289D-01,
351 + 0.27484590D+00, 0.27484590D+00, 0.79756386D-01, 0.79711946D-01,
352 + 0.79711946D-01, 0.79756386D-01, 0.14913414D+01, 0.39625566D-01/
353 DATA ((CALCPI(I,K,8 ),I=1,8), K=10,18) /
354 + 0.27597500D+00, 0.27597500D+00, 0.81130588D-01, 0.81084956D-01,
355 + 0.81084956D-01, 0.81130588D-01, 0.14495215D+01, 0.42577226D-01,
356 + 0.27654656D+00, 0.27654656D+00, 0.82130068D-01, 0.82083931D-01,
357 + 0.82083931D-01, 0.82130068D-01, 0.14020127D+01, 0.45117543D-01,
358 + 0.27657897D+00, 0.27657897D+00, 0.82771279D-01, 0.82725254D-01,
359 + 0.82725254D-01, 0.82771279D-01, 0.13499555D+01, 0.47262777D-01,
360 + 0.27609267D+00, 0.27609267D+00, 0.83072447D-01, 0.83027081D-01,
361 + 0.83027081D-01, 0.83072447D-01, 0.12944158D+01, 0.49030944D-01,
362 + 0.27510991D+00, 0.27510991D+00, 0.83053293D-01, 0.83009072D-01,
363 + 0.83009072D-01, 0.83053293D-01, 0.12363767D+01, 0.50441551D-01,
364 + 0.27365451D+00, 0.27365451D+00, 0.82734758D-01, 0.82692107D-01,
365 + 0.82692107D-01, 0.82734758D-01, 0.11767315D+01, 0.51515331D-01,
366 + 0.27175161D+00, 0.27175161D+00, 0.82138704D-01, 0.82097992D-01,
367 + 0.82097992D-01, 0.82138704D-01, 0.11162806D+01, 0.52273956D-01,
368 + 0.26942735D+00, 0.26942735D+00, 0.81287631D-01, 0.81249172D-01,
369 + 0.81249172D-01, 0.81287631D-01, 0.10557294D+01, 0.52739765D-01,
370 + 0.26670869D+00, 0.26670869D+00, 0.80204389D-01, 0.80168447D-01,
371 + 0.80168447D-01, 0.80204389D-01, 0.99569003D+00, 0.52935485D-01/
372 DATA ((CALCPI(I,K,8 ),I=1,8), K=19,20) /
373 + 0.26362306D+00, 0.26362306D+00, 0.78911912D-01, 0.78878704D-01,
374 + 0.78878704D-01, 0.78911912D-01, 0.93668316D+00, 0.52883968D-01,
375 + 0.26019814D+00, 0.26019814D+00, 0.77432959D-01, 0.77402657D-01,
376 + 0.77402657D-01, 0.77432959D-01, 0.87914238D+00, 0.52607944D-01/
377 DATA ((CALCPI(I,K,9 ),I=1,8), K=1,9) /
378 + 0.25393598D+00, 0.25393598D+00, 0.39934850D-01, 0.39934850D-01,
379 + 0.39934850D-01, 0.39934850D-01, 0.12687427D+01, 0.00000000D+00,
380 + 0.25658251D+00, 0.25658251D+00, 0.42620330D-01, 0.42615827D-01,
381 + 0.42615827D-01, 0.42620330D-01, 0.12634749D+01, 0.44683915D-02,
382 + 0.25860356D+00, 0.25860356D+00, 0.44916190D-01, 0.44908124D-01,
383 + 0.44908124D-01, 0.44916190D-01, 0.12476484D+01, 0.84881719D-02,
384 + 0.26002316D+00, 0.26002316D+00, 0.46840163D-01, 0.46829358D-01,
385 + 0.46829358D-01, 0.46840163D-01, 0.12228611D+01, 0.12079175D-01,
386 + 0.26086160D+00, 0.26086160D+00, 0.48407661D-01, 0.48394850D-01,
387 + 0.48394850D-01, 0.48407661D-01, 0.11904608D+01, 0.15258316D-01,
388 + 0.26114245D+00, 0.26114245D+00, 0.49637207D-01, 0.49623039D-01,
389 + 0.49623039D-01, 0.49637207D-01, 0.11517779D+01, 0.18045530D-01,
390 + 0.26088927D+00, 0.26088927D+00, 0.50547429D-01, 0.50532475D-01,
391 + 0.50532475D-01, 0.50547429D-01, 0.11080570D+01, 0.20460680D-01,
392 + 0.26012660D+00, 0.26012660D+00, 0.51157841D-01, 0.51142601D-01,
393 + 0.51142601D-01, 0.51157841D-01, 0.10604534D+01, 0.22524392D-01,
394 + 0.25887973D+00, 0.25887973D+00, 0.51488547D-01, 0.51473455D-01,
395 + 0.51473455D-01, 0.51488547D-01, 0.10100286D+01, 0.24257785D-01/
396 DATA ((CALCPI(I,K,9 ),I=1,8), K=10,18) /
397 + 0.25717463D+00, 0.25717463D+00, 0.51560086D-01, 0.51545514D-01,
398 + 0.51545514D-01, 0.51560086D-01, 0.95774231D+00, 0.25682336D-01,
399 + 0.25503783D+00, 0.25503783D+00, 0.51393280D-01, 0.51379545D-01,
400 + 0.51379545D-01, 0.51393280D-01, 0.90445031D+00, 0.26819749D-01,
401 + 0.25249633D+00, 0.25249633D+00, 0.51009049D-01, 0.50996415D-01,
402 + 0.50996415D-01, 0.51009049D-01, 0.85090359D+00, 0.27691789D-01,
403 + 0.24957740D+00, 0.24957740D+00, 0.50428241D-01, 0.50416926D-01,
404 + 0.50416926D-01, 0.50428241D-01, 0.79774993D+00, 0.28320122D-01,
405 + 0.24630848D+00, 0.24630848D+00, 0.49671453D-01, 0.49661632D-01,
406 + 0.49661632D-01, 0.49671453D-01, 0.74553722D+00, 0.28726156D-01,
407 + 0.24271699D+00, 0.24271699D+00, 0.48758866D-01, 0.48750672D-01,
408 + 0.48750672D-01, 0.48758866D-01, 0.69471827D+00, 0.28930879D-01,
409 + 0.23883021D+00, 0.23883021D+00, 0.47710084D-01, 0.47703616D-01,
410 + 0.47703616D-01, 0.47710084D-01, 0.64565687D+00, 0.28954713D-01,
411 + 0.23467506D+00, 0.23467506D+00, 0.46543992D-01, 0.46539315D-01,
412 + 0.46539315D-01, 0.46543992D-01, 0.59863493D+00, 0.28817373D-01,
413 + 0.23027803D+00, 0.23027803D+00, 0.45278628D-01, 0.45275779D-01,
414 + 0.45275779D-01, 0.45278628D-01, 0.55386011D+00, 0.28537751D-01/
415 DATA ((CALCPI(I,K,9 ),I=1,8), K=19,20) /
416 + 0.22566501D+00, 0.22566501D+00, 0.43931080D-01, 0.43930068D-01,
417 + 0.43930068D-01, 0.43931080D-01, 0.51147390D+00, 0.28133812D-01,
418 + 0.22086116D+00, 0.22086116D+00, 0.42517397D-01, 0.42518208D-01,
419 + 0.42518208D-01, 0.42517397D-01, 0.47155983D+00, 0.27622508D-01/
420 DATA ((CALCPI(I,K,10),I=1,8), K=1,9) /
421 + 0.25832112D+00, 0.25832112D+00, 0.28960672D-01, 0.28960672D-01,
422 + 0.28960672D-01, 0.28960672D-01, 0.10839859D+01, 0.23852448D-17,
423 + 0.25846394D+00, 0.25846394D+00, 0.30539011D-01, 0.30537508D-01,
424 + 0.30537508D-01, 0.30539011D-01, 0.10505497D+01, 0.31277892D-02,
425 + 0.25803268D+00, 0.25803268D+00, 0.31792680D-01, 0.31790208D-01,
426 + 0.31790208D-01, 0.31792680D-01, 0.10110277D+01, 0.58632227D-02,
427 + 0.25706457D+00, 0.25706457D+00, 0.32748875D-01, 0.32745880D-01,
428 + 0.32745880D-01, 0.32748875D-01, 0.96693666D+00, 0.82367171D-02,
429 + 0.25559005D+00, 0.25559005D+00, 0.33429793D-01, 0.33426657D-01,
430 + 0.33426657D-01, 0.33429793D-01, 0.91947643D+00, 0.10272905D-01,
431 + 0.25364116D+00, 0.25364116D+00, 0.33858861D-01, 0.33855906D-01,
432 + 0.33855906D-01, 0.33858861D-01, 0.86974126D+00, 0.11997545D-01,
433 + 0.25124881D+00, 0.25124881D+00, 0.34058406D-01, 0.34055898D-01,
434 + 0.34055898D-01, 0.34058406D-01, 0.81870325D+00, 0.13435082D-01,
435 + 0.24844356D+00, 0.24844356D+00, 0.34050212D-01, 0.34048367D-01,
436 + 0.34048367D-01, 0.34050212D-01, 0.76720685D+00, 0.14609274D-01,
437 + 0.24525556D+00, 0.24525556D+00, 0.33855439D-01, 0.33854431D-01,
438 + 0.33854431D-01, 0.33855439D-01, 0.71597830D+00, 0.15543132D-01/
439 DATA ((CALCPI(I,K,10),I=1,8), K=10,18) /
440 + 0.24171455D+00, 0.24171455D+00, 0.33494520D-01, 0.33494484D-01,
441 + 0.33494484D-01, 0.33494520D-01, 0.66562790D+00, 0.16258838D-01,
442 + 0.23784979D+00, 0.23784979D+00, 0.32987116D-01, 0.32988149D-01,
443 + 0.32988149D-01, 0.32987116D-01, 0.61665720D+00, 0.16777720D-01,
444 + 0.23369006D+00, 0.23369006D+00, 0.32352031D-01, 0.32354201D-01,
445 + 0.32354201D-01, 0.32352031D-01, 0.56946665D+00, 0.17120189D-01,
446 + 0.22926354D+00, 0.22926354D+00, 0.31607147D-01, 0.31610494D-01,
447 + 0.31610494D-01, 0.31607147D-01, 0.52436404D+00, 0.17305679D-01,
448 + 0.22459780D+00, 0.22459780D+00, 0.30769357D-01, 0.30773896D-01,
449 + 0.30773896D-01, 0.30769357D-01, 0.48157360D+00, 0.17352600D-01,
450 + 0.21971966D+00, 0.21971966D+00, 0.29854510D-01, 0.29860234D-01,
451 + 0.29860234D-01, 0.29854510D-01, 0.44124553D+00, 0.17278283D-01,
452 + 0.21465516D+00, 0.21465516D+00, 0.28877363D-01, 0.28884249D-01,
453 + 0.28884249D-01, 0.28877363D-01, 0.40346553D+00, 0.17098949D-01,
454 + 0.20942947D+00, 0.20942947D+00, 0.27851556D-01, 0.27859564D-01,
455 + 0.27859564D-01, 0.27851556D-01, 0.36826407D+00, 0.16829679D-01,
456 + 0.20406685D+00, 0.20406685D+00, 0.26789589D-01, 0.26798665D-01,
457 + 0.26798665D-01, 0.26789589D-01, 0.33562547D+00, 0.16484402D-01/
458 DATA ((CALCPI(I,K,10),I=1,8), K=19,20) /
459 + 0.19859055D+00, 0.19859055D+00, 0.25702823D-01, 0.25712905D-01,
460 + 0.25712905D-01, 0.25702823D-01, 0.30549636D+00, 0.16075896D-01,
461 + 0.19302282D+00, 0.19302282D+00, 0.24601491D-01, 0.24612506D-01,
462 + 0.24612506D-01, 0.24601491D-01, 0.27779349D+00, 0.15615799D-01/
463 DATA ((CALCPI(I,K,11),I=1,8), K=1,9) /
464 + 0.25965109D+00, 0.25965109D+00, 0.20571046D-01, 0.20571046D-01,
465 + 0.20571046D-01, 0.20571046D-01, 0.91677136D+00, 0.23852448D-17,
466 + 0.25773527D+00, 0.25773527D+00, 0.21504001D-01, 0.21504079D-01,
467 + 0.21504079D-01, 0.21504001D-01, 0.86884765D+00, 0.22017070D-02,
468 + 0.25532059D+00, 0.25532059D+00, 0.22180457D-01, 0.22180859D-01,
469 + 0.22180859D-01, 0.22180457D-01, 0.81844202D+00, 0.40795342D-02,
470 + 0.25245013D+00, 0.25245013D+00, 0.22629676D-01, 0.22630592D-01,
471 + 0.22630592D-01, 0.22629676D-01, 0.76678664D+00, 0.56666584D-02,
472 + 0.24915877D+00, 0.24915877D+00, 0.22874877D-01, 0.22876456D-01,
473 + 0.22876456D-01, 0.22874877D-01, 0.71480295D+00, 0.69893088D-02,
474 + 0.24548165D+00, 0.24548165D+00, 0.22939078D-01, 0.22941430D-01,
475 + 0.22941430D-01, 0.22939078D-01, 0.66326481D+00, 0.80733850D-02,
476 + 0.24145225D+00, 0.24145225D+00, 0.22843625D-01, 0.22846828D-01,
477 + 0.22846828D-01, 0.22843625D-01, 0.61282237D+00, 0.89428821D-02,
478 + 0.23710290D+00, 0.23710290D+00, 0.22608528D-01, 0.22612631D-01,
479 + 0.22612631D-01, 0.22608528D-01, 0.56399578D+00, 0.96202862D-02,
480 + 0.23246494D+00, 0.23246494D+00, 0.22252557D-01, 0.22257586D-01,
481 + 0.22257586D-01, 0.22252557D-01, 0.51719258D+00, 0.10126694D-01/
482 DATA ((CALCPI(I,K,11),I=1,8), K=10,18) /
483 + 0.22756863D+00, 0.22756863D+00, 0.21793214D-01, 0.21799172D-01,
484 + 0.21799172D-01, 0.21793214D-01, 0.47271709D+00, 0.10481803D-01,
485 + 0.22244329D+00, 0.22244329D+00, 0.21246762D-01, 0.21253633D-01,
486 + 0.21253633D-01, 0.21246762D-01, 0.43078194D+00, 0.10703969D-01,
487 + 0.21711722D+00, 0.21711722D+00, 0.20628233D-01, 0.20635988D-01,
488 + 0.20635988D-01, 0.20628233D-01, 0.39151968D+00, 0.10810222D-01,
489 + 0.21161773D+00, 0.21161773D+00, 0.19951446D-01, 0.19960041D-01,
490 + 0.19960041D-01, 0.19951446D-01, 0.35499405D+00, 0.10816298D-01,
491 + 0.20597110D+00, 0.20597110D+00, 0.19229015D-01, 0.19238398D-01,
492 + 0.19238398D-01, 0.19229015D-01, 0.32121090D+00, 0.10736660D-01,
493 + 0.20020258D+00, 0.20020258D+00, 0.18472376D-01, 0.18482486D-01,
494 + 0.18482486D-01, 0.18472376D-01, 0.29012844D+00, 0.10584531D-01,
495 + 0.19433630D+00, 0.19433630D+00, 0.17691814D-01, 0.17702584D-01,
496 + 0.17702584D-01, 0.17691814D-01, 0.26166677D+00, 0.10371926D-01,
497 + 0.18839533D+00, 0.18839533D+00, 0.16896497D-01, 0.16907855D-01,
498 + 0.16907855D-01, 0.16896497D-01, 0.23571664D+00, 0.10109695D-01,
499 + 0.18240157D+00, 0.18240157D+00, 0.16094523D-01, 0.16106395D-01,
500 + 0.16106395D-01, 0.16094523D-01, 0.21214719D+00, 0.98075669D-02/
501 DATA ((CALCPI(I,K,11),I=1,8), K=19,20) /
502 + 0.17637578D+00, 0.17637578D+00, 0.15292969D-01, 0.15305279D-01,
503 + 0.15305279D-01, 0.15292969D-01, 0.19081278D+00, 0.94742091D-02,
504 + 0.17033758D+00, 0.17033758D+00, 0.14497950D-01, 0.14510623D-01,
505 + 0.14510623D-01, 0.14497950D-01, 0.17155890D+00, 0.91172811D-02/
506 DATA ((CALCPI(I,K,12),I=1,8), K=1,9) /
507 + 0.25879811D+00, 0.25879811D+00, 0.14270896D-01, 0.14270896D-01,
508 + 0.14270896D-01, 0.14270896D-01, 0.76643707D+00, 0.11926224D-17,
509 + 0.25510266D+00, 0.25510266D+00, 0.14821666D-01, 0.14822448D-01,
510 + 0.14822448D-01, 0.14821666D-01, 0.71227877D+00, 0.15418105D-02,
511 + 0.25099427D+00, 0.25099427D+00, 0.15176761D-01, 0.15178394D-01,
512 + 0.15178394D-01, 0.15176761D-01, 0.65839243D+00, 0.28266503D-02,
513 + 0.24651783D+00, 0.24651783D+00, 0.15363638D-01, 0.15366162D-01,
514 + 0.15366162D-01, 0.15363638D-01, 0.60565244D+00, 0.38861289D-02,
515 + 0.24170962D+00, 0.24170962D+00, 0.15403683D-01, 0.15407113D-01,
516 + 0.15407113D-01, 0.15403683D-01, 0.55466879D+00, 0.47448491D-02,
517 + 0.23660526D+00, 0.23660526D+00, 0.15317134D-01, 0.15321464D-01,
518 + 0.15321464D-01, 0.15317134D-01, 0.50589349D+00, 0.54261089D-02,
519 + 0.23123850D+00, 0.23123850D+00, 0.15122407D-01, 0.15127614D-01,
520 + 0.15127614D-01, 0.15122407D-01, 0.45967243D+00, 0.59511441D-02,
521 + 0.22564155D+00, 0.22564155D+00, 0.14836254D-01, 0.14842303D-01,
522 + 0.14842303D-01, 0.14836254D-01, 0.41623954D+00, 0.63393360D-02,
523 + 0.21984527D+00, 0.21984527D+00, 0.14473973D-01, 0.14480816D-01,
524 + 0.14480816D-01, 0.14473973D-01, 0.37573706D+00, 0.66084403D-02/
525 DATA ((CALCPI(I,K,12),I=1,8), K=10,18) /
526 + 0.21387917D+00, 0.21387917D+00, 0.14049427D-01, 0.14057009D-01,
527 + 0.14057009D-01, 0.14049427D-01, 0.33822838D+00, 0.67746355D-02,
528 + 0.20777152D+00, 0.20777152D+00, 0.13575140D-01, 0.13583397D-01,
529 + 0.13583397D-01, 0.13575140D-01, 0.30371095D+00, 0.68526308D-02,
530 + 0.20154932D+00, 0.20154932D+00, 0.13062362D-01, 0.13071226D-01,
531 + 0.13071226D-01, 0.13062362D-01, 0.27212883D+00, 0.68557464D-02,
532 + 0.19523836D+00, 0.19523836D+00, 0.12521126D-01, 0.12530527D-01,
533 + 0.12530527D-01, 0.12521126D-01, 0.24338378D+00, 0.67959882D-02,
534 + 0.18886319D+00, 0.18886319D+00, 0.11960319D-01, 0.11970185D-01,
535 + 0.11970185D-01, 0.11960319D-01, 0.21734539D+00, 0.66841218D-02,
536 + 0.18244713D+00, 0.18244713D+00, 0.11387747D-01, 0.11398003D-01,
537 + 0.11398003D-01, 0.11387747D-01, 0.19386016D+00, 0.65297451D-02,
538 + 0.17601222D+00, 0.17601222D+00, 0.10810194D-01, 0.10820769D-01,
539 + 0.10820769D-01, 0.10810194D-01, 0.17275924D+00, 0.63413614D-02,
540 + 0.16957930D+00, 0.16957930D+00, 0.10233498D-01, 0.10244321D-01,
541 + 0.10244321D-01, 0.10233498D-01, 0.15386517D+00, 0.61264529D-02,
542 + 0.16316789D+00, 0.16316789D+00, 0.96626210D-02, 0.96736224D-02,
543 + 0.96736224D-02, 0.96626210D-02, 0.13699740D+00, 0.58915553D-02/
544 DATA ((CALCPI(I,K,12),I=1,8), K=19,20) /
545 + 0.15679627D+00, 0.15679627D+00, 0.91017156D-02, 0.91128305D-02,
546 + 0.91128305D-02, 0.91017156D-02, 0.12197683D+00, 0.56423328D-02,
547 + 0.15048145D+00, 0.15048145D+00, 0.85542025D-02, 0.85653690D-02,
548 + 0.85653690D-02, 0.85542025D-02, 0.10862933D+00, 0.53836526D-02/
549 DATA ((CALCPI(I,K,13),I=1,8), K=1,9) /
550 + 0.25612786D+00, 0.25612786D+00, 0.96354581D-02, 0.96354581D-02,
551 + 0.96354581D-02, 0.96354581D-02, 0.63230272D+00,-0.10842022D-18,
552 + 0.25084444D+00, 0.25084444D+00, 0.99594084D-02, 0.99604949D-02,
553 + 0.99604949D-02, 0.99594084D-02, 0.57717980D+00, 0.10657702D-02,
554 + 0.24524055D+00, 0.24524055D+00, 0.10138463D-01, 0.10140599D-01,
555 + 0.10140599D-01, 0.10138463D-01, 0.52432560D+00, 0.19345699D-02,
556 + 0.23936064D+00, 0.23936064D+00, 0.10196373D-01, 0.10199504D-01,
557 + 0.10199504D-01, 0.10196373D-01, 0.47424860D+00, 0.26341541D-02,
558 + 0.23324074D+00, 0.23324074D+00, 0.10151392D-01, 0.10155456D-01,
559 + 0.10155456D-01, 0.10151392D-01, 0.42725979D+00, 0.31858901D-02,
560 + 0.22691557D+00, 0.22691557D+00, 0.10020105D-01, 0.10025031D-01,
561 + 0.10025031D-01, 0.10020105D-01, 0.38352407D+00, 0.36092957D-02,
562 + 0.22041799D+00, 0.22041799D+00, 0.98173862D-02, 0.98230987D-02,
563 + 0.98230987D-02, 0.98173862D-02, 0.34312781D+00, 0.39219477D-02,
564 + 0.21377912D+00, 0.21377912D+00, 0.95564254D-02, 0.95628451D-02,
565 + 0.95628451D-02, 0.95564254D-02, 0.30607374D+00, 0.41395380D-02,
566 + 0.20702852D+00, 0.20702852D+00, 0.92489818D-02, 0.92560275D-02,
567 + 0.92560275D-02, 0.92489818D-02, 0.27230038D+00, 0.42761561D-02/
568 DATA ((CALCPI(I,K,13),I=1,8), K=10,18) /
569 + 0.20019427D+00, 0.20019427D+00, 0.89054557D-02, 0.89130459D-02,
570 + 0.89130459D-02, 0.89054557D-02, 0.24169604D+00, 0.43443792D-02,
571 + 0.19330303D+00, 0.19330303D+00, 0.85350013D-02, 0.85430549D-02,
572 + 0.85430549D-02, 0.85350013D-02, 0.21411055D+00, 0.43554021D-02,
573 + 0.18638002D+00, 0.18638002D+00, 0.81456233D-02, 0.81540613D-02,
574 + 0.81540613D-02, 0.81456233D-02, 0.18936651D+00, 0.43191475D-02,
575 + 0.17944908D+00, 0.17944908D+00, 0.77442653D-02, 0.77530111D-02,
576 + 0.77530111D-02, 0.77442653D-02, 0.16726861D+00, 0.42443660D-02,
577 + 0.17253267D+00, 0.17253267D+00, 0.73368952D-02, 0.73458755D-02,
578 + 0.73458755D-02, 0.73368952D-02, 0.14761164D+00, 0.41387306D-02,
579 + 0.16565185D+00, 0.16565185D+00, 0.69285874D-02, 0.69377329D-02,
580 + 0.69377329D-02, 0.69285874D-02, 0.13018709D+00, 0.40089273D-02,
581 + 0.15882630D+00, 0.15882630D+00, 0.65236030D-02, 0.65328486D-02,
582 + 0.65328486D-02, 0.65236030D-02, 0.11478856D+00, 0.38607412D-02,
583 + 0.15207430D+00, 0.15207430D+00, 0.61254672D-02, 0.61347526D-02,
584 + 0.61347526D-02, 0.61254672D-02, 0.10121588D+00, 0.36991402D-02,
585 + 0.14541274D+00, 0.14541274D+00, 0.57370450D-02, 0.57463147D-02,
586 + 0.57463147D-02, 0.57370450D-02, 0.89278292D-01, 0.35283545D-02/
587 DATA ((CALCPI(I,K,13),I=1,8), K=19,20) /
588 + 0.13885713D+00, 0.13885713D+00, 0.53606133D-02, 0.53698169D-02,
589 + 0.53698169D-02, 0.53606133D-02, 0.78796647D-01, 0.33519532D-02,
590 + 0.13242156D+00, 0.13242156D+00, 0.49979302D-02, 0.50070221D-02,
591 + 0.50070221D-02, 0.49979302D-02, 0.69604844D-01, 0.31729167D-02/
592 DATA ((CALCPI(I,K,14),I=1,8), K=1,9) /
593 + 0.25177353D+00, 0.25177353D+00, 0.63042492D-02, 0.63042492D-02,
594 + 0.63042492D-02, 0.63042492D-02, 0.51366824D+00,-0.32526065D-18,
595 + 0.24504651D+00, 0.24504651D+00, 0.64940417D-02, 0.64952158D-02,
596 + 0.64952158D-02, 0.64940417D-02, 0.46100530D+00, 0.72215456D-03,
597 + 0.23809737D+00, 0.23809737D+00, 0.65797224D-02, 0.65819651D-02,
598 + 0.65819651D-02, 0.65797224D-02, 0.41195408D+00, 0.12984830D-02,
599 + 0.23096882D+00, 0.23096882D+00, 0.65805932D-02, 0.65837976D-02,
600 + 0.65837976D-02, 0.65805932D-02, 0.36668482D+00, 0.17518714D-02,
601 + 0.22369566D+00, 0.22369566D+00, 0.65113670D-02, 0.65154268D-02,
602 + 0.65154268D-02, 0.65113670D-02, 0.32524568D+00, 0.20998136D-02,
603 + 0.21631089D+00, 0.21631089D+00, 0.63848855D-02, 0.63896963D-02,
604 + 0.63896963D-02, 0.63848855D-02, 0.28756533D+00, 0.23577495D-02,
605 + 0.20884573D+00, 0.20884573D+00, 0.62125425D-02, 0.62180031D-02,
606 + 0.62180031D-02, 0.62125425D-02, 0.25352523D+00, 0.25394573D-02,
607 + 0.20132955D+00, 0.20132955D+00, 0.60042098D-02, 0.60102228D-02,
608 + 0.60102228D-02, 0.60042098D-02, 0.22295563D+00, 0.26570028D-02,
609 + 0.19379003D+00, 0.19379003D+00, 0.57684967D-02, 0.57749692D-02,
610 + 0.57749692D-02, 0.57684967D-02, 0.19565154D+00, 0.27210248D-02/
611 DATA ((CALCPI(I,K,14),I=1,8), K=10,18) /
612 + 0.18625325D+00, 0.18625325D+00, 0.55128456D-02, 0.55196899D-02,
613 + 0.55196899D-02, 0.55128456D-02, 0.17138562D+00, 0.27408469D-02,
614 + 0.17874369D+00, 0.17874369D+00, 0.52436469D-02, 0.52507806D-02,
615 + 0.52507806D-02, 0.52436469D-02, 0.14991744D+00, 0.27246091D-02,
616 + 0.17128429D+00, 0.17128429D+00, 0.49663453D-02, 0.49736921D-02,
617 + 0.49736921D-02, 0.49663453D-02, 0.13100207D+00, 0.26793859D-02,
618 + 0.16389643D+00, 0.16389643D+00, 0.46855324D-02, 0.46930216D-02,
619 + 0.46930216D-02, 0.46855324D-02, 0.11439664D+00, 0.26112898D-02,
620 + 0.15659995D+00, 0.15659995D+00, 0.44050347D-02, 0.44126020D-02,
621 + 0.44126020D-02, 0.44050347D-02, 0.99865574D-01, 0.25255689D-02,
622 + 0.14941313D+00, 0.14941313D+00, 0.41279969D-02, 0.41355839D-02,
623 + 0.41355839D-02, 0.41279969D-02, 0.87184474D-01, 0.24266972D-02,
624 + 0.14235275D+00, 0.14235275D+00, 0.38569581D-02, 0.38645125D-02,
625 + 0.38645125D-02, 0.38569581D-02, 0.76142868D-01, 0.23184578D-02,
626 + 0.13543401D+00, 0.13543401D+00, 0.35939252D-02, 0.36014005D-02,
627 + 0.36014005D-02, 0.35939252D-02, 0.66545972D-01, 0.22040221D-02,
628 + 0.12867062D+00, 0.12867062D+00, 0.33404406D-02, 0.33477961D-02,
629 + 0.33477961D-02, 0.33404406D-02, 0.58215662D-01, 0.20860219D-02/
630 DATA ((CALCPI(I,K,14),I=1,8), K=19,20) /
631 + 0.12207475D+00, 0.12207475D+00, 0.30976447D-02, 0.31048451D-02,
632 + 0.31048451D-02, 0.30976447D-02, 0.50990753D-01, 0.19666172D-02,
633 + 0.11565711D+00, 0.11565711D+00, 0.28663337D-02, 0.28733488D-02,
634 + 0.28733488D-02, 0.28663337D-02, 0.44726750D-01, 0.18475573D-02/
635 DATA ((CALCPI(I,K,15),I=1,8), K=1,9) /
636 + 0.24575242D+00, 0.24575242D+00, 0.39751539D-02, 0.39751539D-02,
637 + 0.39751539D-02, 0.39751539D-02, 0.40981099D+00, 0.65052130D-18,
638 + 0.23770237D+00, 0.23770237D+00, 0.40858460D-02, 0.40869562D-02,
639 + 0.40869562D-02, 0.40858460D-02, 0.36179506D+00, 0.47613270D-03,
640 + 0.22953528D+00, 0.22953528D+00, 0.41242411D-02, 0.41263228D-02,
641 + 0.41263228D-02, 0.41242411D-02, 0.31817903D+00, 0.84836616D-03,
642 + 0.22129122D+00, 0.22129122D+00, 0.41051569D-02, 0.41080800D-02,
643 + 0.41080800D-02, 0.41051569D-02, 0.27884363D+00, 0.11345171D-02,
644 + 0.21300294D+00, 0.21300294D+00, 0.40398494D-02, 0.40434915D-02,
645 + 0.40434915D-02, 0.40398494D-02, 0.24362265D+00, 0.13481608D-02,
646 + 0.20470108D+00, 0.20470108D+00, 0.39377405D-02, 0.39419872D-02,
647 + 0.39419872D-02, 0.39377405D-02, 0.21226565D+00, 0.15008645D-02,
648 + 0.19641453D+00, 0.19641453D+00, 0.38071279D-02, 0.38118734D-02,
649 + 0.38118734D-02, 0.38071279D-02, 0.18450588D+00, 0.16028860D-02,
650 + 0.18817021D+00, 0.18817021D+00, 0.36550500D-02, 0.36601966D-02,
651 + 0.36601966D-02, 0.36550500D-02, 0.16005751D+00, 0.16630537D-02,
652 + 0.17999326D+00, 0.17999326D+00, 0.34875156D-02, 0.34929741D-02,
653 + 0.34929741D-02, 0.34875156D-02, 0.13862640D+00, 0.16890199D-02/
654 DATA ((CALCPI(I,K,15),I=1,8), K=10,18) /
655 + 0.17190702D+00, 0.17190702D+00, 0.33096082D-02, 0.33152973D-02,
656 + 0.33152973D-02, 0.33096082D-02, 0.11992046D+00, 0.16873797D-02,
657 + 0.16393313D+00, 0.16393313D+00, 0.31255879D-02, 0.31314343D-02,
658 + 0.31314343D-02, 0.31255879D-02, 0.10365558D+00, 0.16637865D-02,
659 + 0.15609148D+00, 0.15609148D+00, 0.29389907D-02, 0.29449289D-02,
660 + 0.29449289D-02, 0.29389907D-02, 0.89560899D-01, 0.16230630D-02,
661 + 0.14840025D+00, 0.14840025D+00, 0.27527137D-02, 0.27586853D-02,
662 + 0.27586853D-02, 0.27527137D-02, 0.77382365D-01, 0.15692957D-02,
663 + 0.14087590D+00, 0.14087590D+00, 0.25690932D-02, 0.25750472D-02,
664 + 0.25750472D-02, 0.25690932D-02, 0.66885018D-01, 0.15059217D-02,
665 + 0.13353322D+00, 0.13353322D+00, 0.23899777D-02, 0.23958695D-02,
666 + 0.23958695D-02, 0.23899777D-02, 0.57854248D-01, 0.14358084D-02,
667 + 0.12638528D+00, 0.12638528D+00, 0.22167932D-02, 0.22225844D-02,
668 + 0.22225844D-02, 0.22167932D-02, 0.50096238D-01, 0.13613257D-02,
669 + 0.11944349D+00, 0.11944349D+00, 0.20506038D-02, 0.20562621D-02,
670 + 0.20562621D-02, 0.20506038D-02, 0.43437702D-01, 0.12844116D-02,
671 + 0.11271763D+00, 0.11271763D+00, 0.18921664D-02, 0.18976647D-02,
672 + 0.18976647D-02, 0.18921664D-02, 0.37725089D-01, 0.12066320D-02/
673 DATA ((CALCPI(I,K,15),I=1,8), K=19,20) /
674 + 0.10621587D+00, 0.10621587D+00, 0.17419795D-02, 0.17472958D-02,
675 + 0.17472958D-02, 0.17419795D-02, 0.32823472D-01, 0.11292341D-02,
676 + 0.99944809D-01, 0.99944809D-01, 0.16003276D-02, 0.16054443D-02,
677 + 0.16054443D-02, 0.16003276D-02, 0.28615198D-01, 0.10531942D-02/
678 DATA ((CALCPI(I,K,16),I=1,8), K=1,9) /
679 + 0.23802032D+00, 0.23802032D+00, 0.23986819D-02, 0.23986819D-02,
680 + 0.23986819D-02, 0.23986819D-02, 0.31998394D+00,-0.10842022D-18,
681 + 0.22875953D+00, 0.22875953D+00, 0.24626918D-02, 0.24636621D-02,
682 + 0.24636621D-02, 0.24626918D-02, 0.27792199D+00, 0.30282026D-03,
683 + 0.21949492D+00, 0.21949492D+00, 0.24785578D-02, 0.24803491D-02,
684 + 0.24803491D-02, 0.24785578D-02, 0.24058610D+00, 0.53488192D-03,
685 + 0.21026301D+00, 0.21026301D+00, 0.24569821D-02, 0.24594604D-02,
686 + 0.24594604D-02, 0.24569821D-02, 0.20762793D+00, 0.70923900D-03,
687 + 0.20109372D+00, 0.20109372D+00, 0.24061201D-02, 0.24091636D-02,
688 + 0.24091636D-02, 0.24061201D-02, 0.17872072D+00, 0.83585721D-03,
689 + 0.19201449D+00, 0.19201449D+00, 0.23324957D-02, 0.23359947D-02,
690 + 0.23359947D-02, 0.23324957D-02, 0.15349297D+00, 0.92291695D-03,
691 + 0.18305101D+00, 0.18305101D+00, 0.22418345D-02, 0.22456909D-02,
692 + 0.22456909D-02, 0.22418345D-02, 0.13158424D+00, 0.97764917D-03,
693 + 0.17422686D+00, 0.17422686D+00, 0.21389009D-02, 0.21430272D-02,
694 + 0.21430272D-02, 0.21389009D-02, 0.11264366D+00, 0.10061823D-02,
695 + 0.16556367D+00, 0.16556367D+00, 0.20276779D-02, 0.20319965D-02,
696 + 0.20319965D-02, 0.20276779D-02, 0.96334645D-01, 0.10137412D-02/
697 DATA ((CALCPI(I,K,16),I=1,8), K=10,18) /
698 + 0.15708114D+00, 0.15708114D+00, 0.19114664D-02, 0.19159092D-02,
699 + 0.19159092D-02, 0.19114664D-02, 0.82341921D-01, 0.10047596D-02,
700 + 0.14879705D+00, 0.14879705D+00, 0.17929667D-02, 0.17974743D-02,
701 + 0.17974743D-02, 0.17929667D-02, 0.70373881D-01, 0.98297281D-03,
702 + 0.14072726D+00, 0.14072726D+00, 0.16743616D-02, 0.16788826D-02,
703 + 0.16788826D-02, 0.16743616D-02, 0.60164413D-01, 0.95151031D-03,
704 + 0.13288576D+00, 0.13288576D+00, 0.15573858D-02, 0.15618762D-02,
705 + 0.15618762D-02, 0.15573858D-02, 0.51473534D-01, 0.91297372D-03,
706 + 0.12528468D+00, 0.12528468D+00, 0.14433891D-02, 0.14478117D-02,
707 + 0.14478117D-02, 0.14433891D-02, 0.44087108D-01, 0.86950738D-03,
708 + 0.11793429D+00, 0.11793429D+00, 0.13333942D-02, 0.13377180D-02,
709 + 0.13377180D-02, 0.13333942D-02, 0.37815928D-01, 0.82286184D-03,
710 + 0.11084309D+00, 0.11084309D+00, 0.12281469D-02, 0.12323465D-02,
711 + 0.12323465D-02, 0.12281469D-02, 0.32494234D-01, 0.77445028D-03,
712 + 0.10401779D+00, 0.10401779D+00, 0.11281622D-02, 0.11322172D-02,
713 + 0.11322172D-02, 0.11281622D-02, 0.27977991D-01, 0.72539909D-03,
714 + 0.97463439D-01, 0.97463439D-01, 0.10337644D-02, 0.10376590D-02,
715 + 0.10376590D-02, 0.10337644D-02, 0.24143013D-01, 0.67659252D-03/
716 DATA ((CALCPI(I,K,16),I=1,8), K=19,20) /
717 + 0.91183419D-01, 0.91183419D-01, 0.94512214D-03, 0.94884464D-03,
718 + 0.94884464D-03, 0.94512214D-03, 0.20882905D-01, 0.62871192D-03,
719 + 0.85179574D-01, 0.85179574D-01, 0.86227955D-03, 0.86582162D-03,
720 + 0.86582162D-03, 0.86227955D-03, 0.18107040D-01, 0.58226995D-03/
721 DATA ((CALCPI(I,K,17),I=1,8), K=1,9) /
722 + 0.22849684D+00, 0.22849684D+00, 0.13723279D-02, 0.13723279D-02,
723 + 0.13723279D-02, 0.13723279D-02, 0.24341185D+00, 0.20599841D-17,
724 + 0.21814046D+00, 0.21814046D+00, 0.14085736D-02, 0.14093568D-02,
725 + 0.14093568D-02, 0.14085736D-02, 0.20795831D+00, 0.18369815D-03,
726 + 0.20790326D+00, 0.20790326D+00, 0.14144576D-02, 0.14158835D-02,
727 + 0.14158835D-02, 0.14144576D-02, 0.17718140D+00, 0.32182400D-03,
728 + 0.19781707D+00, 0.19781707D+00, 0.13971715D-02, 0.13991185D-02,
729 + 0.13991185D-02, 0.13971715D-02, 0.15057030D+00, 0.42329321D-03,
730 + 0.18790787D+00, 0.18790787D+00, 0.13622673D-02, 0.13646280D-02,
731 + 0.13646280D-02, 0.13622673D-02, 0.12769314D+00, 0.49497766D-03,
732 + 0.17819872D+00, 0.17819872D+00, 0.13139597D-02, 0.13166402D-02,
733 + 0.13166402D-02, 0.13139597D-02, 0.10811239D+00, 0.54228532D-03,
734 + 0.16871092D+00, 0.16871092D+00, 0.12559563D-02, 0.12588747D-02,
735 + 0.12588747D-02, 0.12559563D-02, 0.91423399D-01, 0.57000173D-03,
736 + 0.15946347D+00, 0.15946347D+00, 0.11912927D-02, 0.11943780D-02,
737 + 0.11943780D-02, 0.11912927D-02, 0.77254110D-01, 0.58213097D-03,
738 + 0.15047321D+00, 0.15047321D+00, 0.11224534D-02, 0.11256447D-02,
739 + 0.11256447D-02, 0.11224534D-02, 0.65263768D-01, 0.58203167D-03/
740 DATA ((CALCPI(I,K,17),I=1,8), K=10,18) /
741 + 0.14175490D+00, 0.14175490D+00, 0.10514589D-02, 0.10547040D-02,
742 + 0.10547040D-02, 0.10514589D-02, 0.55146147D-01, 0.57251426D-03,
743 + 0.13332119D+00, 0.13332119D+00, 0.97992289D-03, 0.98317770D-03,
744 + 0.98317770D-03, 0.97992289D-03, 0.46628587D-01, 0.55590776D-03,
745 + 0.12518267D+00, 0.12518267D+00, 0.90911521D-03, 0.91234284D-03,
746 + 0.91234284D-03, 0.90911521D-03, 0.39470762D-01, 0.53412985D-03,
747 + 0.11734793D+00, 0.11734793D+00, 0.84001331D-03, 0.84318330D-03,
748 + 0.84318330D-03, 0.84001331D-03, 0.33462929D-01, 0.50874553D-03,
749 + 0.10982357D+00, 0.10982357D+00, 0.77334844D-03, 0.77643606D-03,
750 + 0.77643606D-03, 0.77334844D-03, 0.28423477D-01, 0.48101891D-03,
751 + 0.10261430D+00, 0.10261430D+00, 0.70964688D-03, 0.71263245D-03,
752 + 0.71263245D-03, 0.70964688D-03, 0.24196408D-01, 0.45195960D-03,
753 + 0.95722986D-01, 0.95722986D-01, 0.64926600D-03, 0.65213430D-03,
754 + 0.65213430D-03, 0.64926600D-03, 0.20648794D-01, 0.42236335D-03,
755 + 0.89150754D-01, 0.89150754D-01, 0.59242565D-03, 0.59516536D-03,
756 + 0.59516536D-03, 0.59242565D-03, 0.17668059D-01, 0.39284740D-03,
757 + 0.82897064D-01, 0.82897064D-01, 0.53923530D-03, 0.54183849D-03,
758 + 0.54183849D-03, 0.53923530D-03, 0.15159374D-01, 0.36388128D-03/
759 DATA ((CALCPI(I,K,17),I=1,8), K=19,20) /
760 + 0.76959826D-01, 0.76959826D-01, 0.48971741D-03, 0.49217904D-03,
761 + 0.49217904D-03, 0.48971741D-03, 0.13043364D-01, 0.33581336D-03,
762 + 0.71335508D-01, 0.71335508D-01, 0.44382708D-03, 0.44614461D-03,
763 + 0.44614461D-03, 0.44382708D-03, 0.11253923D-01, 0.30889356D-03/
764 DATA ((CALCPI(I,K,18),I=1,8), K=1,9) /
765 + 0.21707483D+00, 0.21707483D+00, 0.73510510D-03, 0.73510510D-03,
766 + 0.73510510D-03, 0.73510510D-03, 0.17928739D+00,-0.35914197D-18,
767 + 0.20575033D+00, 0.20575033D+00, 0.75465779D-03, 0.75523717D-03,
768 + 0.75523717D-03, 0.75465779D-03, 0.15059863D+00, 0.10462869D-03,
769 + 0.19467944D+00, 0.19467944D+00, 0.75645398D-03, 0.75749891D-03,
770 + 0.75749891D-03, 0.75645398D-03, 0.12624355D+00, 0.18197158D-03,
771 + 0.18388760D+00, 0.18388760D+00, 0.74487900D-03, 0.74629248D-03,
772 + 0.74629248D-03, 0.74487900D-03, 0.10561678D+00, 0.23758744D-03,
773 + 0.17339526D+00, 0.17339526D+00, 0.72342403D-03, 0.72512180D-03,
774 + 0.72512180D-03, 0.72342403D-03, 0.88234668D-01, 0.27585990D-03,
775 + 0.16321947D+00, 0.16321947D+00, 0.69457725D-03, 0.69648679D-03,
776 + 0.69648679D-03, 0.69457725D-03, 0.73643508D-01, 0.30007333D-03,
777 + 0.15337545D+00, 0.15337545D+00, 0.66055268D-03, 0.66261193D-03,
778 + 0.66261193D-03, 0.66055268D-03, 0.61437212D-01, 0.31315650D-03,
779 + 0.14387589D+00, 0.14387589D+00, 0.62314382D-03, 0.62530004D-03,
780 + 0.62530004D-03, 0.62314382D-03, 0.51258263D-01, 0.31753900D-03,
781 + 0.13473123D+00, 0.13473123D+00, 0.58378747D-03, 0.58599623D-03,
782 + 0.58599623D-03, 0.58378747D-03, 0.42791125D-01, 0.31522635D-03/
783 DATA ((CALCPI(I,K,18),I=1,8), K=10,18) /
784 + 0.12594964D+00, 0.12594964D+00, 0.54363293D-03, 0.54585716D-03,
785 + 0.54585716D-03, 0.54363293D-03, 0.35761787D-01, 0.30787701D-03,
786 + 0.11753706D+00, 0.11753706D+00, 0.50357781D-03, 0.50578694D-03,
787 + 0.50578694D-03, 0.50357781D-03, 0.29934302D-01, 0.29684490D-03,
788 + 0.10949734D+00, 0.10949734D+00, 0.46430950D-03, 0.46647869D-03,
789 + 0.46647869D-03, 0.46430950D-03, 0.25106893D-01, 0.28322667D-03,
790 + 0.10183227D+00, 0.10183227D+00, 0.42633993D-03, 0.42844940D-03,
791 + 0.42844940D-03, 0.42633993D-03, 0.21108414D-01, 0.26790152D-03,
792 + 0.94541657D-01, 0.94541657D-01, 0.39003592D-03, 0.39207023D-03,
793 + 0.39207023D-03, 0.39003592D-03, 0.17794878D-01, 0.25156582D-03,
794 + 0.87623499D-01, 0.87623499D-01, 0.35564601D-03, 0.35759352D-03,
795 + 0.35759352D-03, 0.35564601D-03, 0.15045868D-01, 0.23476386D-03,
796 + 0.81074038D-01, 0.81074038D-01, 0.32332339D-03, 0.32517571D-03,
797 + 0.32517571D-03, 0.32332339D-03, 0.12761147D-01, 0.21791416D-03,
798 + 0.74887922D-01, 0.74887922D-01, 0.29314573D-03, 0.29489728D-03,
799 + 0.29489728D-03, 0.29314573D-03, 0.10857858D-01, 0.20133223D-03,
800 + 0.69058330D-01, 0.69058330D-01, 0.26513197D-03, 0.26677950D-03,
801 + 0.26677950D-03, 0.26513197D-03, 0.92678860D-02, 0.18525014D-03/
802 DATA ((CALCPI(I,K,18),I=1,8), K=19,20) /
803 + 0.63577119D-01, 0.63577119D-01, 0.23925617D-03, 0.24079840D-03,
804 + 0.24079840D-03, 0.23925617D-03, 0.79353281D-02, 0.16983269D-03,
805 + 0.58434971D-01, 0.58434971D-01, 0.21545905D-03, 0.21689632D-03,
806 + 0.21689632D-03, 0.21545905D-03, 0.68144746D-02, 0.15519119D-03/
807 DATA ((CALCPI(I,K,19),I=1,8), K=1,9) /
808 + 0.20361916D+00, 0.20361916D+00, 0.36223666D-03, 0.36223666D-03,
809 + 0.36223666D-03, 0.36223666D-03, 0.12676601D+00, 0.19718927D-17,
810 + 0.19147582D+00, 0.19147582D+00, 0.37174717D-03, 0.37217687D-03,
811 + 0.37217687D-03, 0.37174717D-03, 0.10461111D+00, 0.54642443D-04,
812 + 0.17973344D+00, 0.17973344D+00, 0.37197592D-03, 0.37274080D-03,
813 + 0.37274080D-03, 0.37197592D-03, 0.86231484D-01, 0.94521885D-04,
814 + 0.16840854D+00, 0.16840854D+00, 0.36523130D-03, 0.36625293D-03,
815 + 0.36625293D-03, 0.36523130D-03, 0.70993300D-01, 0.12267888D-03,
816 + 0.15751378D+00, 0.15751378D+00, 0.35348077D-03, 0.35469262D-03,
817 + 0.35469262D-03, 0.35348077D-03, 0.58410377D-01, 0.14163351D-03,
818 + 0.14705787D+00, 0.14705787D+00, 0.33801272D-03, 0.33935887D-03,
819 + 0.33935887D-03, 0.33801272D-03, 0.48055083D-01, 0.15315995D-03,
820 + 0.13704767D+00, 0.13704767D+00, 0.32000812D-03, 0.32144197D-03,
821 + 0.32144197D-03, 0.32000812D-03, 0.39554499D-01, 0.15887226D-03,
822 + 0.12748741D+00, 0.12748741D+00, 0.30042465D-03, 0.30190764D-03,
823 + 0.30190764D-03, 0.30042465D-03, 0.32592378D-01, 0.16010956D-03,
824 + 0.11837899D+00, 0.11837899D+00, 0.28001302D-03, 0.28151364D-03,
825 + 0.28151364D-03, 0.28001302D-03, 0.26899082D-01, 0.15796012D-03/
826 DATA ((CALCPI(I,K,19),I=1,8), K=10,18) /
827 + 0.10972203D+00, 0.10972203D+00, 0.25936718D-03, 0.26085996D-03,
828 + 0.26085996D-03, 0.25936718D-03, 0.22247529D-01, 0.15331696D-03,
829 + 0.10151397D+00, 0.10151397D+00, 0.23894237D-03, 0.24040709D-03,
830 + 0.24040709D-03, 0.23894237D-03, 0.18447916D-01, 0.14690052D-03,
831 + 0.93750249D-01, 0.93750249D-01, 0.21907852D-03, 0.22049943D-03,
832 + 0.22049943D-03, 0.21907852D-03, 0.15342800D-01, 0.13928598D-03,
833 + 0.86424367D-01, 0.86424367D-01, 0.20002104D-03, 0.20138622D-03,
834 + 0.20138622D-03, 0.20002104D-03, 0.12802401D-01, 0.13092753D-03,
835 + 0.79528091D-01, 0.79528091D-01, 0.18193789D-03, 0.18323866D-03,
836 + 0.18323866D-03, 0.18193789D-03, 0.10720084D-01, 0.12217838D-03,
837 + 0.73051590D-01, 0.73051590D-01, 0.16493526D-03, 0.16616567D-03,
838 + 0.16616567D-03, 0.16493526D-03, 0.90089315D-02, 0.11330909D-03,
839 + 0.66983604D-01, 0.66983604D-01, 0.14907078D-03, 0.15022714D-03,
840 + 0.15022714D-03, 0.14907078D-03, 0.75985466D-02, 0.10452301D-03,
841 + 0.61311617D-01, 0.61311617D-01, 0.13436434D-03, 0.13544483D-03,
842 + 0.13544483D-03, 0.13436434D-03, 0.64319118D-02, 0.95969183D-04,
843 + 0.56022033D-01, 0.56022033D-01, 0.12080736D-03, 0.12181167D-03,
844 + 0.12181167D-03, 0.12080736D-03, 0.54630107D-02, 0.87753354D-04/
845 DATA ((CALCPI(I,K,19),I=1,8), K=19,20) /
846 + 0.51100347D-01, 0.51100347D-01, 0.10837037D-03, 0.10929942D-03,
847 + 0.10929942D-03, 0.10837037D-03, 0.46549386D-02, 0.79947206D-04,
848 + 0.46531326D-01, 0.46531326D-01, 0.97009053D-04, 0.97864705D-04,
849 + 0.97864705D-04, 0.97009053D-04, 0.39780773D-02, 0.72595761D-04/
850 DATA ((CALCPI(I,K,20),I=1,8), K=1,9) /
851 + 0.18795535D+00, 0.18795535D+00, 0.16002003D-03, 0.16002003D-03,
852 + 0.16002003D-03, 0.16002003D-03, 0.84958910D-01, 0.20803129D-17,
853 + 0.17517565D+00, 0.17517565D+00, 0.16366860D-03, 0.16396317D-03,
854 + 0.16396317D-03, 0.16366860D-03, 0.68804471D-01, 0.25133578D-04,
855 + 0.16295836D+00, 0.16295836D+00, 0.16331040D-03, 0.16382720D-03,
856 + 0.16382720D-03, 0.16331040D-03, 0.55727437D-01, 0.43430035D-04,
857 + 0.15130733D+00, 0.15130733D+00, 0.15986156D-03, 0.16054231D-03,
858 + 0.16054231D-03, 0.15986156D-03, 0.45126075D-01, 0.56215712D-04,
859 + 0.14022419D+00, 0.14022419D+00, 0.15425394D-03, 0.15505054D-03,
860 + 0.15505054D-03, 0.15425394D-03, 0.36554058D-01, 0.64729405D-04,
861 + 0.12970621D+00, 0.12970621D+00, 0.14702170D-03, 0.14789477D-03,
862 + 0.14789477D-03, 0.14702170D-03, 0.29642711D-01, 0.69770569D-04,
863 + 0.11974891D+00, 0.11974891D+00, 0.13869255D-03, 0.13961021D-03,
864 + 0.13961021D-03, 0.13869255D-03, 0.24077970D-01, 0.72101631D-04,
865 + 0.11034530D+00, 0.11034530D+00, 0.12970743D-03, 0.13064411D-03,
866 + 0.13064411D-03, 0.12970743D-03, 0.19602592D-01, 0.72367902D-04,
867 + 0.10148624D+00, 0.10148624D+00, 0.12040478D-03, 0.12134025D-03,
868 + 0.12134025D-03, 0.12040478D-03, 0.16004097D-01, 0.71086585D-04/
869 DATA ((CALCPI(I,K,20),I=1,8), K=10,18) /
870 + 0.93160516D-01, 0.93160516D-01, 0.11105260D-03, 0.11197118D-03,
871 + 0.11197118D-03, 0.11105260D-03, 0.13109239D-01, 0.68682310D-04,
872 + 0.85355089D-01, 0.85355089D-01, 0.10185503D-03, 0.10274477D-03,
873 + 0.10274477D-03, 0.10185503D-03, 0.10777823D-01, 0.65496416D-04,
874 + 0.78055239D-01, 0.78055239D-01, 0.92961348D-04, 0.93813482D-04,
875 + 0.93813482D-04, 0.92961348D-04, 0.88962548D-02, 0.61798506D-04,
876 + 0.71244778D-01, 0.71244778D-01, 0.84476855D-04, 0.85285199D-04,
877 + 0.85285199D-04, 0.84476855D-04, 0.73734434D-02, 0.57799327D-04,
878 + 0.64906248D-01, 0.64906248D-01, 0.76471165D-04, 0.77231676D-04,
879 + 0.77231676D-04, 0.76471165D-04, 0.61369375D-02, 0.53660941D-04,
880 + 0.59021124D-01, 0.59021124D-01, 0.68985444D-04, 0.69695809D-04,
881 + 0.69695809D-04, 0.68985444D-04, 0.51289226D-02, 0.49505510D-04,
882 + 0.53570023D-01, 0.53570023D-01, 0.62038867D-04, 0.62698166D-04,
883 + 0.62698166D-04, 0.62038867D-04, 0.43034132D-02, 0.45423131D-04,
884 + 0.48532899D-01, 0.48532899D-01, 0.55634069D-04, 0.56242481D-04,
885 + 0.56242481D-04, 0.55634069D-04, 0.36241649D-02, 0.41478495D-04,
886 + 0.43889247D-01, 0.43889247D-01, 0.49761407D-04, 0.50319962D-04,
887 + 0.50319962D-04, 0.49761407D-04, 0.30625534D-02, 0.37716203D-04/
888 DATA ((CALCPI(I,K,20),I=1,8), K=19,20) /
889 + 0.39618284D-01, 0.39618284D-01, 0.44402426D-04, 0.44912797D-04,
890 + 0.44912797D-04, 0.44402426D-04, 0.25958712D-02, 0.34165165D-04,
891 + 0.35699128D-01, 0.35699128D-01, 0.39532719D-04, 0.39997051D-04,
892 + 0.39997051D-04, 0.39532719D-04, 0.22061789D-02, 0.30842283D-04/
893 DATA ((CALCPI(I,K,21),I=1,8), K=1,9) /
894 + 0.16984505D+00, 0.16984505D+00, 0.60885465D-04, 0.60885465D-04,
895 + 0.60885465D-04, 0.60885465D-04, 0.52924211D-01, 0.88345536D-18,
896 + 0.15666146D+00, 0.15666146D+00, 0.61446509D-04, 0.61623934D-04,
897 + 0.61623934D-04, 0.61446509D-04, 0.42003073D-01, 0.93958376D-05,
898 + 0.14421618D+00, 0.14421618D+00, 0.60899758D-04, 0.61206837D-04,
899 + 0.61206837D-04, 0.60899758D-04, 0.33396444D-01, 0.16418275D-04,
900 + 0.13249492D+00, 0.13249492D+00, 0.59362264D-04, 0.59761509D-04,
901 + 0.59761509D-04, 0.59362264D-04, 0.26585821D-01, 0.21383842D-04,
902 + 0.12148386D+00, 0.12148386D+00, 0.57143167D-04, 0.57604407D-04,
903 + 0.57604407D-04, 0.57143167D-04, 0.21198459D-01, 0.24749454D-04,
904 + 0.11116464D+00, 0.11116464D+00, 0.54369512D-04, 0.54868619D-04,
905 + 0.54868619D-04, 0.54369512D-04, 0.16947054D-01, 0.26767364D-04,
906 + 0.10151771D+00, 0.10151771D+00, 0.51206680D-04, 0.51724642D-04,
907 + 0.51724642D-04, 0.51206680D-04, 0.13590929D-01, 0.27712781D-04,
908 + 0.92521541D-01, 0.92521541D-01, 0.47812039D-04, 0.48334076D-04,
909 + 0.48334076D-04, 0.47812039D-04, 0.10940195D-01, 0.27837482D-04,
910 + 0.84153018D-01, 0.84153018D-01, 0.44305192D-04, 0.44820015D-04,
911 + 0.44820015D-04, 0.44305192D-04, 0.88441370D-02, 0.27342560D-04/
912 DATA ((CALCPI(I,K,21),I=1,8), K=10,18) /
913 + 0.76387680D-01, 0.76387680D-01, 0.40783839D-04, 0.41283030D-04,
914 + 0.41283030D-04, 0.40783839D-04, 0.71824561D-02, 0.26396185D-04,
915 + 0.69199949D-01, 0.69199949D-01, 0.37324313D-04, 0.37801800D-04,
916 + 0.37801800D-04, 0.37324313D-04, 0.58607578D-02, 0.25135626D-04,
917 + 0.62563370D-01, 0.62563370D-01, 0.33982934D-04, 0.34434547D-04,
918 + 0.34434547D-04, 0.33982934D-04, 0.48056069D-02, 0.23669904D-04,
919 + 0.56450853D-01, 0.56450853D-01, 0.30799491D-04, 0.31222586D-04,
920 + 0.31222586D-04, 0.30799491D-04, 0.39594263D-02, 0.22084388D-04,
921 + 0.50834903D-01, 0.50834903D-01, 0.27800110D-04, 0.28193257D-04,
922 + 0.28193257D-04, 0.27800110D-04, 0.32771831D-02, 0.20444607D-04,
923 + 0.45687844D-01, 0.45687844D-01, 0.25000224D-04, 0.25362943D-04,
924 + 0.25362943D-04, 0.25000224D-04, 0.27241266D-02, 0.18800043D-04,
925 + 0.40982029D-01, 0.40982029D-01, 0.22406888D-04, 0.22739427D-04,
926 + 0.22739427D-04, 0.22406888D-04, 0.22733375D-02, 0.17187147D-04,
927 + 0.36690042D-01, 0.36690042D-01, 0.20020736D-04, 0.20323893D-04,
928 + 0.20323893D-04, 0.20020736D-04, 0.19037882D-02, 0.15631911D-04,
929 + 0.32784875D-01, 0.32784875D-01, 0.17837766D-04, 0.18112738D-04,
930 + 0.18112738D-04, 0.17837766D-04, 0.15991541D-02, 0.14152190D-04/
931 DATA ((CALCPI(I,K,21),I=1,8), K=19,20) /
932 + 0.29240095D-01, 0.29240095D-01, 0.15850674D-04, 0.16098936D-04,
933 + 0.16098936D-04, 0.15850674D-04, 0.13467598D-02, 0.12759488D-04,
934 + 0.26029990D-01, 0.26029990D-01, 0.14049880D-04, 0.14273088D-04,
935 + 0.14273088D-04, 0.14049880D-04, 0.11366537D-02, 0.11460395D-04/
936 DATA ((CALCPI(I,K,22),I=1,8), K=1,9) /
937 + 0.14893693D+00, 0.14893693D+00, 0.18659088D-04, 0.18659088D-04,
938 + 0.18659088D-04, 0.18659088D-04, 0.29653500D-01, 0.10452387D-17,
939 + 0.13565988D+00, 0.13565988D+00, 0.17993271D-04, 0.18088547D-04,
940 + 0.18088547D-04, 0.17993271D-04, 0.23024856D-01, 0.22765974D-05,
941 + 0.12331002D+00, 0.12331002D+00, 0.17476260D-04, 0.17638657D-04,
942 + 0.17638657D-04, 0.17476260D-04, 0.17956049D-01, 0.42411576D-05,
943 + 0.11184701D+00, 0.11184701D+00, 0.16890258D-04, 0.17098264D-04,
944 + 0.17098264D-04, 0.16890258D-04, 0.14050693D-01, 0.57552835D-05,
945 + 0.10123550D+00, 0.10123550D+00, 0.16241964D-04, 0.16478715D-04,
946 + 0.16478715D-04, 0.16241964D-04, 0.11030654D-01, 0.68838036D-05,
947 + 0.91436148D-01, 0.91436148D-01, 0.15493096D-04, 0.15745494D-04,
948 + 0.15745494D-04, 0.15493096D-04, 0.86995268D-02, 0.76364989D-05,
949 + 0.82409700D-01, 0.82409700D-01, 0.14649555D-04, 0.14907639D-04,
950 + 0.14907639D-04, 0.14649555D-04, 0.68965978D-02, 0.80610204D-05,
951 + 0.74116264D-01, 0.74116264D-01, 0.13740791D-04, 0.13997113D-04,
952 + 0.13997113D-04, 0.13740791D-04, 0.54973300D-02, 0.82221692D-05,
953 + 0.66515671D-01, 0.66515671D-01, 0.12790392D-04, 0.13039529D-04,
954 + 0.13039529D-04, 0.12790392D-04, 0.44067959D-02, 0.81735671D-05/
955 DATA ((CALCPI(I,K,22),I=1,8), K=10,18) /
956 + 0.59567782D-01, 0.59567782D-01, 0.11822215D-04, 0.12060361D-04,
957 + 0.12060361D-04, 0.11822215D-04, 0.35530492D-02, 0.79646535D-05,
958 + 0.53232735D-01, 0.53232735D-01, 0.10858081D-04, 0.11082705D-04,
959 + 0.11082705D-04, 0.10858081D-04, 0.28810796D-02, 0.76390325D-05,
960 + 0.47471183D-01, 0.47471183D-01, 0.99150661D-05, 0.10124631D-04,
961 + 0.10124631D-04, 0.99150661D-05, 0.23486757D-02, 0.72323212D-05,
962 + 0.42244523D-01, 0.42244523D-01, 0.90064951D-05, 0.92002321D-05,
963 + 0.92002321D-05, 0.90064951D-05, 0.19241141D-02, 0.67736298D-05,
964 + 0.37515095D-01, 0.37515095D-01, 0.81421094D-05, 0.83198279D-05,
965 + 0.83198279D-05, 0.81421094D-05, 0.15834187D-02, 0.62861281D-05,
966 + 0.33246370D-01, 0.33246370D-01, 0.73283836D-05, 0.74903217D-05,
967 + 0.74903217D-05, 0.73283836D-05, 0.13081901D-02, 0.57877078D-05,
968 + 0.29403104D-01, 0.29403104D-01, 0.65692553D-05, 0.67159583D-05,
969 + 0.67159583D-05, 0.65692553D-05, 0.10844362D-02, 0.52920045D-05,
970 + 0.25951472D-01, 0.25951472D-01, 0.58665925D-05, 0.59988169D-05,
971 + 0.59988169D-05, 0.58665925D-05, 0.90150889D-03, 0.48091163D-05,
972 + 0.22859180D-01, 0.22859180D-01, 0.52206265D-05, 0.53392645D-05,
973 + 0.53392645D-05, 0.52206265D-05, 0.75117828D-03, 0.43462499D-05/
974 DATA ((CALCPI(I,K,22),I=1,8), K=19,20) /
975 + 0.20095550D-01, 0.20095550D-01, 0.46303832D-05, 0.47364044D-05,
976 + 0.47364044D-05, 0.46303832D-05, 0.62705922D-03, 0.39083370D-05,
977 + 0.17631577D-01, 0.17631577D-01, 0.40939948D-05, 0.41884015D-05,
978 + 0.41884015D-05, 0.40939948D-05, 0.52418344D-03, 0.34985020D-05/
979 DATA ((CALCPI(I,K,23),I=1,8), K=1,9) /
980 + 0.12466155D+00, 0.12466155D+00, 0.40617655D-05, 0.40617655D-05,
981 + 0.40617655D-05, 0.40617655D-05, 0.14051896D-01, 0.73495990D-18,
982 + 0.11173314D+00, 0.11173314D+00, 0.32938963D-05, 0.33358710D-05,
983 + 0.33358710D-05, 0.32938963D-05, 0.10660771D-01,-0.64809452D-07,
984 + 0.99927656D-01, 0.99927656D-01, 0.29399618D-05, 0.30092196D-05,
985 + 0.30092196D-05, 0.29399618D-05, 0.81551927D-02, 0.16269877D-06,
986 + 0.89167355D-01, 0.89167355D-01, 0.27551315D-05, 0.28413834D-05,
987 + 0.28413834D-05, 0.27551315D-05, 0.62804841D-02, 0.46082146D-06,
988 + 0.79387013D-01, 0.79387013D-01, 0.26597156D-05, 0.27556416D-05,
989 + 0.27556416D-05, 0.26597156D-05, 0.48642075D-02, 0.76656066D-06,
990 + 0.70519683D-01, 0.70519683D-01, 0.25930830D-05, 0.26933197D-05,
991 + 0.26933197D-05, 0.25930830D-05, 0.37947917D-02, 0.10326492D-05,
992 + 0.62501104D-01, 0.62501104D-01, 0.25218351D-05, 0.26226000D-05,
993 + 0.26226000D-05, 0.25218351D-05, 0.29816936D-02, 0.12371630D-05,
994 + 0.55268986D-01, 0.55268986D-01, 0.24372334D-05, 0.25359336D-05,
995 + 0.25359336D-05, 0.24372334D-05, 0.23591941D-02, 0.13809056D-05,
996 + 0.48763274D-01, 0.48763274D-01, 0.23354218D-05, 0.24303250D-05,
997 + 0.24303250D-05, 0.23354218D-05, 0.18796033D-02, 0.14679651D-05/
998 DATA ((CALCPI(I,K,23),I=1,8), K=10,18) /
999 + 0.42926390D-01, 0.42926390D-01, 0.22165796D-05, 0.23065871D-05,
1000 + 0.23065871D-05, 0.22165796D-05, 0.15068345D-02, 0.15051628D-05,
1001 + 0.37703340D-01, 0.37703340D-01, 0.20844329D-05, 0.21689076D-05,
1002 + 0.21689076D-05, 0.20844329D-05, 0.12147442D-02, 0.15017730D-05,
1003 + 0.33041846D-01, 0.33041846D-01, 0.19431609D-05, 0.20217950D-05,
1004 + 0.20217950D-05, 0.19431609D-05, 0.98419305D-03, 0.14666069D-05,
1005 + 0.28892432D-01, 0.28892432D-01, 0.17967575D-05, 0.18694754D-05,
1006 + 0.18694754D-05, 0.17967575D-05, 0.80077746D-03, 0.14075252D-05,
1007 + 0.25208481D-01, 0.25208481D-01, 0.16489621D-05, 0.17158482D-05,
1008 + 0.17158482D-05, 0.16489621D-05, 0.65382251D-03, 0.13314928D-05,
1009 + 0.21946264D-01, 0.21946264D-01, 0.15029021D-05, 0.15641478D-05,
1010 + 0.15641478D-05, 0.15029021D-05, 0.53537088D-03, 0.12443253D-05,
1011 + 0.19064943D-01, 0.19064943D-01, 0.13610616D-05, 0.14169271D-05,
1012 + 0.14169271D-05, 0.13610616D-05, 0.43936516D-03, 0.11507452D-05,
1013 + 0.16526541D-01, 0.16526541D-01, 0.12253560D-05, 0.12761424D-05,
1014 + 0.12761424D-05, 0.12253560D-05, 0.36118695D-03, 0.10545294D-05,
1015 + 0.14295900D-01, 0.14295900D-01, 0.10971589D-05, 0.11431890D-05,
1016 + 0.11431890D-05, 0.10971589D-05, 0.29728882D-03, 0.95859812D-06/
1017 DATA ((CALCPI(I,K,23),I=1,8), K=19,20) /
1018 + 0.12340602D-01, 0.12340602D-01, 0.97737632D-06, 0.10189805D-05,
1019 + 0.10189805D-05, 0.97737632D-06, 0.24490074D-03, 0.86513980D-06,
1020 + 0.10630890D-01, 0.10630890D-01, 0.86653311D-06, 0.90404004D-06,
1021 + 0.90404004D-06, 0.86653311D-06, 0.20184479D-03, 0.77574186D-06/
1022 DATA ((CALCPI(I,K,24),I=1,8), K=1,9) /
1023 + 0.95962050D-01, 0.95962050D-01, 0.47362220D-06, 0.47362220D-06,
1024 + 0.47362220D-06, 0.47362220D-06, 0.49045878D-02, 0.20394965D-17,
1025 + 0.84079722D-01, 0.84079722D-01, 0.70130062D-07, 0.69988528D-07,
1026 + 0.69988528D-07, 0.70130062D-07, 0.36451349D-02,-0.31177713D-06,
1027 + 0.73503277D-01, 0.73503277D-01,-0.93949559D-07,-0.90343984D-07,
1028 + -0.90343984D-07,-0.93949559D-07, 0.27512246D-02,-0.40201316D-06,
1029 + 0.64100842D-01, 0.64100842D-01,-0.14362678D-06,-0.13436729D-06,
1030 + -0.13436729D-06,-0.14362678D-06, 0.21021321D-02,-0.39189456D-06,
1031 + 0.55766799D-01, 0.55766799D-01,-0.13315620D-06,-0.11760078D-06,
1032 + -0.11760078D-06,-0.13315620D-06, 0.16167642D-02,-0.33297146D-06,
1033 + 0.48398988D-01, 0.48398988D-01,-0.92538893D-07,-0.70742324D-07,
1034 + -0.70742323D-07,-0.92538893D-07, 0.12544540D-02,-0.25317481D-06,
1035 + 0.41902807D-01, 0.41902807D-01,-0.42348392D-07,-0.14880959D-07,
1036 + -0.14880958D-07,-0.42348392D-07, 0.98199191D-03,-0.17132676D-06,
1037 + 0.36190734D-01, 0.36190734D-01, 0.74222474D-08, 0.39687164D-07,
1038 + 0.39687164D-07, 0.74222474D-08, 0.77374477D-03,-0.95999516D-07,
1039 + 0.31181740D-01, 0.31181740D-01, 0.51602539D-07, 0.87648579D-07,
1040 + 0.87648579D-07, 0.51602539D-07, 0.61314725D-03,-0.31213106D-07/
1041 DATA ((CALCPI(I,K,24),I=1,8), K=10,18) /
1042 + 0.26801252D-01, 0.26801252D-01, 0.87520725D-07, 0.12629318D-06,
1043 + 0.12629318D-06, 0.87520725D-07, 0.48832221D-03, 0.21298934D-07,
1044 + 0.22980857D-01, 0.22980857D-01, 0.11469355D-06, 0.15517350D-06,
1045 + 0.15517350D-06, 0.11469355D-06, 0.39039186D-03, 0.61816996D-07,
1046 + 0.19658021D-01, 0.19658021D-01, 0.13365167D-06, 0.17490421D-06,
1047 + 0.17490421D-06, 0.13365167D-06, 0.31302545D-03, 0.91492748D-07,
1048 + 0.16775799D-01, 0.16775799D-01, 0.14529177D-06, 0.18649424D-06,
1049 + 0.18649424D-06, 0.14529177D-06, 0.25156265D-03, 0.11172781D-06,
1050 + 0.14282533D-01, 0.14282533D-01, 0.15076037D-06, 0.19121612D-06,
1051 + 0.19121612D-06, 0.15076037D-06, 0.20248137D-03, 0.12407887D-06,
1052 + 0.12131537D-01, 0.12131537D-01, 0.15126112D-06, 0.19040264D-06,
1053 + 0.19040264D-06, 0.15126112D-06, 0.16314146D-03, 0.13008240D-06,
1054 + 0.10280782D-01, 0.10280782D-01, 0.14791039D-06, 0.18529515D-06,
1055 + 0.18529515D-06, 0.14791039D-06, 0.13152577D-03, 0.13112458D-06,
1056 + 0.86925746D-02, 0.86925746D-02, 0.14170656D-06, 0.17700791D-06,
1057 + 0.17700791D-06, 0.14170656D-06, 0.10606210D-03, 0.12842218D-06,
1058 + 0.73332400D-02, 0.73332400D-02, 0.13351233D-06, 0.16650745D-06,
1059 + 0.16650745D-06, 0.13351233D-06, 0.85523934D-04, 0.12301445D-06/
1060 DATA ((CALCPI(I,K,24),I=1,8), K=19,20) /
1061 + 0.61728066D-02, 0.61728066D-02, 0.12404355D-06, 0.15459991D-06,
1062 + 0.15459991D-06, 0.12404355D-06, 0.68945956D-04, 0.11575973D-06,
1063 + 0.51847001D-02, 0.51847001D-02, 0.11387807D-06, 0.14193947D-06,
1064 + 0.14193947D-06, 0.11387807D-06, 0.55558553D-04, 0.10735073D-06/
1065 DATA ((CALCPI(I,K,25),I=1,8), K=1,9) /
1066 + 0.60325058D-01, 0.60325058D-01, 0.12021897D-07, 0.12021897D-07,
1067 + 0.12021897D-07, 0.12021897D-07, 0.81120147D-03,-0.11270876D-17,
1068 + 0.50834153D-01, 0.50834153D-01,-0.36274860D-07, 0.32872280D-07,
1069 + 0.32872280D-07,-0.36274860D-07, 0.61140969D-03,-0.45571468D-07,
1070 + 0.42744651D-01, 0.42744651D-01,-0.34258286D-07, 0.80910456D-07,
1071 + 0.80910456D-07,-0.34258286D-07, 0.46844302D-03,-0.41456923D-07,
1072 + 0.35848234D-01, 0.35848234D-01,-0.10842241D-07, 0.13313926D-06,
1073 + 0.13313926D-06,-0.10842241D-07, 0.36321200D-03,-0.16412313D-07,
1074 + 0.29985251D-01, 0.29985251D-01, 0.18938907D-07, 0.17867746D-06,
1075 + 0.17867746D-06, 0.18938907D-07, 0.28205770D-03, 0.14636623D-07,
1076 + 0.25015626D-01, 0.25015626D-01, 0.49079739D-07, 0.21480938D-06,
1077 + 0.21480938D-06, 0.49079739D-07, 0.21978589D-03, 0.45760655D-07,
1078 + 0.20814681D-01, 0.20814681D-01, 0.75514551D-07, 0.24016718D-06,
1079 + 0.24016718D-06, 0.75514551D-07, 0.17211389D-03, 0.72957264D-07,
1080 + 0.17273565D-01, 0.17273565D-01, 0.96144340D-07, 0.25477275D-06,
1081 + 0.25477275D-06, 0.96144340D-07, 0.13490269D-03, 0.94176960D-07,
1082 + 0.14297126D-01, 0.14297126D-01, 0.11060699D-06, 0.25992244D-06,
1083 + 0.25992244D-06, 0.11060699D-06, 0.10566281D-03, 0.10909573D-06/
1084 DATA ((CALCPI(I,K,25),I=1,8), K=10,18) /
1085 + 0.11802457D-01, 0.11802457D-01, 0.11923453D-06, 0.25722643D-06,
1086 + 0.25722643D-06, 0.11923453D-06, 0.82766725D-04, 0.11807544D-06,
1087 + 0.97175665D-02, 0.97175665D-02, 0.12275219D-06, 0.24837617D-06,
1088 + 0.24837617D-06, 0.12275219D-06, 0.64789615D-04, 0.12186462D-06,
1089 + 0.79801385D-02, 0.79801385D-02, 0.12209600D-06, 0.23502257D-06,
1090 + 0.23502257D-06, 0.12209600D-06, 0.50635791D-04, 0.12141745D-06,
1091 + 0.65364096D-02, 0.65364096D-02, 0.11820595D-06, 0.21861936D-06,
1092 + 0.21861936D-06, 0.11820595D-06, 0.39514631D-04, 0.11768805D-06,
1093 + 0.53401547D-02, 0.53401547D-02, 0.11195715D-06, 0.20039654D-06,
1094 + 0.20039654D-06, 0.11195715D-06, 0.30792202D-04, 0.11156250D-06,
1095 + 0.43517661D-02, 0.43517661D-02, 0.10412636D-06, 0.18136202D-06,
1096 + 0.18136202D-06, 0.10412636D-06, 0.23955915D-04, 0.10382614D-06,
1097 + 0.35374248D-02, 0.35374248D-02, 0.95360916D-07, 0.16229946D-06,
1098 + 0.16229946D-06, 0.95360916D-07, 0.18607293D-04, 0.95132905D-07,
1099 + 0.28683557D-02, 0.28683557D-02, 0.86178300D-07, 0.14379144D-06,
1100 + 0.14379144D-06, 0.86178300D-07, 0.14431662D-04, 0.86005414D-07,
1101 + 0.23201619D-02, 0.23201619D-02, 0.76978650D-07, 0.12625102D-06,
1102 + 0.12625102D-06, 0.76978650D-07, 0.11177135D-04, 0.76847775D-07/
1103 DATA ((CALCPI(I,K,25),I=1,8), K=19,20) /
1104 + 0.18722314D-02, 0.18722314D-02, 0.68056829D-07, 0.10994897D-06,
1105 + 0.10994897D-06, 0.68056829D-07, 0.86446436D-05, 0.67957915D-07,
1106 + 0.15072110D-02, 0.15072110D-02, 0.59617793D-07, 0.95041092D-07,
1107 + 0.95041092D-07, 0.59617793D-07, 0.66774538D-05, 0.59543153D-07/
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 AURPI3(1,0,X,SB,XPDF(1))
1136 CALL AURPI3(2,0,X,SB,XPDF(2))
1137 CALL AURPI3(3,0,X,SB,XPDF(3))
1138 CALL AURPI3(4,0,X,SB,XPDF(4))
1139 CALL AURPI3(5,0,X,SB,XPDF(5))
1140 CALL AURPI3(8,0,X,SB,XPDF(6))
1141 CALL AURPI3(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