]> git.uio.no Git - u/mrichter/AliRoot.git/blob - PDF/spdf/abfkw3.F
added the delete of EMCAL object posted in the folder when new file is opened
[u/mrichter/AliRoot.git] / PDF / spdf / abfkw3.F
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/
21 C----------------------------------------------------------------------
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/
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 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))
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