]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PDF/spdf/acfgp2.F
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Added SetUserDecay routine. When a particle decays the standard MC decay
[u/mrichter/AliRoot.git] / PDF / spdf / acfgp2.F
CommitLineData
fe4da5cc 1*
2* $Id$
3*
4* $Log$
5* Revision 1.1.1.3 1996/11/07 11:52:35 cernlib
6* Import version 7.05
7*
8* Revision 1.3 1996/11/07 09:35:09 pdflib
9* New version 7.05: cut routine Ctq4RTb in 9 pieces
10*
11* Revision 1.1.1.1 1996/04/12 15:28:57 plothow
12* Version 7.01
13*
14*
15#include "pdf/pilot.h"
16 SUBROUTINE ACFGP2(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DGL)
17C
18C INTERPOLATION PROGRAM WHICH INTERPOLATES THE GRID "DATAGA" AND GIVES THE
19C QUARK AND GLUON DISTRIBUTIONS IN THE REAL PHOTON, AS FUNCTIONS OF X AND Q2
20C
21C THE Q2-EVOLUTION IS PERFORMED WITH BLL AP-EQUATIONS AND NF=4. A MASSIVE
22C CHARM DISTRIBUTION (BORROWED FROM GLUCK AND REYA) IS ALSO AVAILABLE.
23C
24C THE BOUNDARY CONDITIONS ARE SUCH THAT THE DISTRIBUTION FUNCTIONS ARE GIVEN
25C BY A VDM "ANSATZ" AT Q2=.25 GEV**2.
26C
27C THE PROGRAM WORKS FOR 2. GEV**2 < Q2 <5.5E+5 AND .00137 < X < .9986
28C
29C THE DISTRIBUTIONS ARE CALCULATED IN THE MSBAR FACTORIZATION SCHEME.
30C
31C THE VALUE OF LAMBDA-MSB IS 200 MEV
32C
33C THE OUTPUT IS WRITTEN IN THE FILE 'FILEOUT':
34C X*U=X*U(X,Q2)
35C X*D= ...
36C X*S= ...
37C X*C= ... (MASSLESS CHARM WITH C(X,2.)=0)
38C X*CM= ... (MASSIVE CHARM WITH MC=1.5 GEV )
39C X*GLU=GLUON(X,Q2)*X
40C
41C
42C F2 = PHOTON STRUCTURE FUNCTION WITHOUT CHARM
43C F2C= " " " WITH MASSIVE CHARM
44C
45#include "pdf/expdp.inc"
46 + DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL
47 REAL X, Q, UV, DV, UB, DB, SB, CB, BB, GL
48 REAL Q2
49 DIMENSION XPDF(7),PAR(30),CALC(8,20,32),CALCO(8,20,32)
50 COMMON/W5051I7/CALCO
51 EXTERNAL AFCPLU
52 DATA ZERO/0.0/
53C----------------------------------------------------------------------
54 DATA (PAR(I),I=1,30) /
55 + 0.2000E+00, 0.4800E+00, 0.8500E+00, 0.1000E+01,
56 + 0.1000E+01, 0.1000E+01, 0.1200E+01, 0.7500E+01,
57 + 0.1987E+01, 0.1900E+01, 0.0000E+00, 0.0000E+00,
58 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.1000E+01,
59 + 0.0000E+00, 0.1000E+01, 0.0000E+00, 0.0000E+00,
60 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
61 + 0.4000E+01, 0.0000E+00, 0.1000E+01, 0.1000E+01,
62 + 0.8000E-01, 0.2000E+01/
63 DATA ((CALC(I,K,1 ),I=1,8), K=1,9) /
64 + 0.0000E+00, 0.0000E+00, 0.1671E-02, -0.5026E-03,
65 + 0.6657E-03, -0.5026E-03, 0.6683E-02, 0.9320E-02,
66 + 0.0000E+00, 0.0000E+00, 0.1669E-02, -0.4992E-03,
67 + 0.6709E-03, -0.4992E-03, 0.7711E-02, 0.1350E-01,
68 + 0.0000E+00, 0.0000E+00, 0.1667E-02, -0.4955E-03,
69 + 0.6764E-03, -0.4955E-03, 0.8915E-02, 0.1832E-01,
70 + 0.0000E+00, 0.0000E+00, 0.1665E-02, -0.4917E-03,
71 + 0.6822E-03, -0.4917E-03, 0.1031E-01, 0.2383E-01,
72 + 0.0000E+00, 0.0000E+00, 0.1663E-02, -0.4875E-03,
73 + 0.6883E-03, -0.4875E-03, 0.1192E-01, 0.3008E-01,
74 + 0.0000E+00, 0.0000E+00, 0.1661E-02, -0.4831E-03,
75 + 0.6947E-03, -0.4831E-03, 0.1375E-01, 0.3713E-01,
76 + 0.0000E+00, 0.0000E+00, 0.1658E-02, -0.4784E-03,
77 + 0.7014E-03, -0.4784E-03, 0.1583E-01, 0.4503E-01,
78 + 0.0000E+00, 0.0000E+00, 0.1655E-02, -0.4733E-03,
79 + 0.7085E-03, -0.4733E-03, 0.1817E-01, 0.5383E-01,
80 + 0.0000E+00, 0.0000E+00, 0.1651E-02, -0.4679E-03,
81 + 0.7159E-03, -0.4679E-03, 0.2079E-01, 0.6360E-01/
82 DATA ((CALC(I,K,1 ),I=1,8), K=10,18) /
83 + 0.0000E+00, 0.0000E+00, 0.1647E-02, -0.4621E-03,
84 + 0.7238E-03, -0.4621E-03, 0.2373E-01, 0.7441E-01,
85 + 0.0000E+00, 0.0000E+00, 0.1643E-02, -0.4558E-03,
86 + 0.7321E-03, -0.4558E-03, 0.2700E-01, 0.8633E-01,
87 + 0.0000E+00, 0.0000E+00, 0.1639E-02, -0.4492E-03,
88 + 0.7409E-03, -0.4492E-03, 0.3062E-01, 0.9943E-01,
89 + 0.0000E+00, 0.0000E+00, 0.1633E-02, -0.4420E-03,
90 + 0.7502E-03, -0.4420E-03, 0.3464E-01, 0.1138E+00,
91 + 0.0000E+00, 0.0000E+00, 0.1628E-02, -0.4343E-03,
92 + 0.7600E-03, -0.4343E-03, 0.3907E-01, 0.1295E+00,
93 + 0.0000E+00, 0.0000E+00, 0.1622E-02, -0.4260E-03,
94 + 0.7704E-03, -0.4260E-03, 0.4394E-01, 0.1467E+00,
95 + 0.0000E+00, 0.0000E+00, 0.1615E-02, -0.4172E-03,
96 + 0.7815E-03, -0.4172E-03, 0.4931E-01, 0.1655E+00,
97 + 0.0000E+00, 0.0000E+00, 0.1608E-02, -0.4076E-03,
98 + 0.7932E-03, -0.4076E-03, 0.5520E-01, 0.1859E+00,
99 + 0.0000E+00, 0.0000E+00, 0.1599E-02, -0.3973E-03,
100 + 0.8057E-03, -0.3973E-03, 0.6165E-01, 0.2082E+00/
101 DATA ((CALC(I,K,1 ),I=1,8), K=19,20) /
102 + 0.0000E+00, 0.0000E+00, 0.1591E-02, -0.3863E-03,
103 + 0.8189E-03, -0.3863E-03, 0.6871E-01, 0.2324E+00,
104 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
105 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
106 DATA ((CALC(I,K,2 ),I=1,8), K=1,9) /
107 + 0.0000E+00, 0.0000E+00, 0.1675E-02, -0.4282E-03,
108 + 0.8187E-03, -0.4282E-03, 0.6700E-02, 0.8422E-02,
109 + 0.0000E+00, 0.0000E+00, 0.1668E-02, -0.4181E-03,
110 + 0.8319E-03, -0.4181E-03, 0.7379E-02, 0.1103E-01,
111 + 0.0000E+00, 0.0000E+00, 0.1661E-02, -0.4075E-03,
112 + 0.8457E-03, -0.4075E-03, 0.8140E-02, 0.1387E-01,
113 + 0.0000E+00, 0.0000E+00, 0.1653E-02, -0.3961E-03,
114 + 0.8603E-03, -0.3961E-03, 0.8986E-02, 0.1696E-01,
115 + 0.0000E+00, 0.0000E+00, 0.1644E-02, -0.3838E-03,
116 + 0.8757E-03, -0.3838E-03, 0.9923E-02, 0.2029E-01,
117 + 0.0000E+00, 0.0000E+00, 0.1634E-02, -0.3707E-03,
118 + 0.8919E-03, -0.3707E-03, 0.1096E-01, 0.2388E-01,
119 + 0.0000E+00, 0.0000E+00, 0.1623E-02, -0.3567E-03,
120 + 0.9091E-03, -0.3567E-03, 0.1209E-01, 0.2775E-01,
121 + 0.0000E+00, 0.0000E+00, 0.1611E-02, -0.3416E-03,
122 + 0.9272E-03, -0.3416E-03, 0.1332E-01, 0.3192E-01,
123 + 0.0000E+00, 0.0000E+00, 0.1598E-02, -0.3254E-03,
124 + 0.9465E-03, -0.3254E-03, 0.1467E-01, 0.3639E-01/
125 DATA ((CALC(I,K,2 ),I=1,8), K=10,18) /
126 + 0.0000E+00, 0.0000E+00, 0.1584E-02, -0.3080E-03,
127 + 0.9669E-03, -0.3080E-03, 0.1614E-01, 0.4119E-01,
128 + 0.0000E+00, 0.0000E+00, 0.1568E-02, -0.2893E-03,
129 + 0.9885E-03, -0.2893E-03, 0.1773E-01, 0.4634E-01,
130 + 0.0000E+00, 0.0000E+00, 0.1551E-02, -0.2691E-03,
131 + 0.1012E-02, -0.2691E-03, 0.1946E-01, 0.5188E-01,
132 + 0.0000E+00, 0.0000E+00, 0.1532E-02, -0.2475E-03,
133 + 0.1036E-02, -0.2475E-03, 0.2133E-01, 0.5782E-01,
134 + 0.0000E+00, 0.0000E+00, 0.1512E-02, -0.2242E-03,
135 + 0.1062E-02, -0.2242E-03, 0.2335E-01, 0.6421E-01,
136 + 0.0000E+00, 0.0000E+00, 0.1490E-02, -0.1991E-03,
137 + 0.1090E-02, -0.1991E-03, 0.2554E-01, 0.7109E-01,
138 + 0.0000E+00, 0.0000E+00, 0.1465E-02, -0.1721E-03,
139 + 0.1120E-02, -0.1721E-03, 0.2790E-01, 0.7848E-01,
140 + 0.0000E+00, 0.0000E+00, 0.1439E-02, -0.1429E-03,
141 + 0.1151E-02, -0.1429E-03, 0.3046E-01, 0.8645E-01,
142 + 0.0000E+00, 0.0000E+00, 0.1410E-02, -0.1115E-03,
143 + 0.1185E-02, -0.1115E-03, 0.3321E-01, 0.9504E-01/
144 DATA ((CALC(I,K,2 ),I=1,8), K=19,20) /
145 + 0.0000E+00, 0.0000E+00, 0.1379E-02, -0.7764E-04,
146 + 0.1222E-02, -0.7764E-04, 0.3619E-01, 0.1043E+00,
147 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
148 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
149 DATA ((CALC(I,K,3 ),I=1,8), K=1,9) /
150 + 0.0000E+00, 0.0000E+00, 0.1659E-02, -0.3433E-03,
151 + 0.9723E-03, -0.3433E-03, 0.6635E-02, 0.7843E-02,
152 + 0.0000E+00, 0.0000E+00, 0.1645E-02, -0.3259E-03,
153 + 0.9931E-03, -0.3259E-03, 0.7163E-02, 0.9673E-02,
154 + 0.0000E+00, 0.0000E+00, 0.1630E-02, -0.3073E-03,
155 + 0.1015E-02, -0.3073E-03, 0.7742E-02, 0.1161E-01,
156 + 0.0000E+00, 0.0000E+00, 0.1613E-02, -0.2873E-03,
157 + 0.1038E-02, -0.2873E-03, 0.8372E-02, 0.1364E-01,
158 + 0.0000E+00, 0.0000E+00, 0.1594E-02, -0.2658E-03,
159 + 0.1063E-02, -0.2658E-03, 0.9056E-02, 0.1579E-01,
160 + 0.0000E+00, 0.0000E+00, 0.1574E-02, -0.2428E-03,
161 + 0.1089E-02, -0.2428E-03, 0.9796E-02, 0.1806E-01,
162 + 0.0000E+00, 0.0000E+00, 0.1553E-02, -0.2181E-03,
163 + 0.1116E-02, -0.2181E-03, 0.1060E-01, 0.2046E-01,
164 + 0.0000E+00, 0.0000E+00, 0.1529E-02, -0.1915E-03,
165 + 0.1146E-02, -0.1915E-03, 0.1146E-01, 0.2299E-01,
166 + 0.0000E+00, 0.0000E+00, 0.1503E-02, -0.1629E-03,
167 + 0.1177E-02, -0.1629E-03, 0.1239E-01, 0.2567E-01/
168 DATA ((CALC(I,K,3 ),I=1,8), K=10,18) /
169 + 0.0000E+00, 0.0000E+00, 0.1475E-02, -0.1321E-03,
170 + 0.1210E-02, -0.1321E-03, 0.1339E-01, 0.2852E-01,
171 + 0.0000E+00, 0.0000E+00, 0.1444E-02, -0.9894E-04,
172 + 0.1246E-02, -0.9894E-04, 0.1446E-01, 0.3155E-01,
173 + 0.0000E+00, 0.0000E+00, 0.1410E-02, -0.6324E-04,
174 + 0.1284E-02, -0.6324E-04, 0.1561E-01, 0.3478E-01,
175 + 0.0000E+00, 0.0000E+00, 0.1374E-02, -0.2477E-04,
176 + 0.1324E-02, -0.2477E-04, 0.1686E-01, 0.3822E-01,
177 + 0.0000E+00, 0.0000E+00, 0.1335E-02, 0.1669E-04,
178 + 0.1368E-02, 0.1669E-04, 0.1819E-01, 0.4191E-01,
179 + 0.0000E+00, 0.0000E+00, 0.1292E-02, 0.6140E-04,
180 + 0.1414E-02, 0.6140E-04, 0.1962E-01, 0.4586E-01,
181 + 0.0000E+00, 0.0000E+00, 0.1245E-02, 0.1096E-03,
182 + 0.1464E-02, 0.1096E-03, 0.2117E-01, 0.5011E-01,
183 + 0.0000E+00, 0.0000E+00, 0.1194E-02, 0.1617E-03,
184 + 0.1517E-02, 0.1617E-03, 0.2283E-01, 0.5469E-01,
185 + 0.0000E+00, 0.0000E+00, 0.1140E-02, 0.2178E-03,
186 + 0.1575E-02, 0.2178E-03, 0.2462E-01, 0.5962E-01/
187 DATA ((CALC(I,K,3 ),I=1,8), K=19,20) /
188 + 0.0000E+00, 0.0000E+00, 0.1080E-02, 0.2785E-03,
189 + 0.1637E-02, 0.2785E-03, 0.2655E-01, 0.6494E-01,
190 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
191 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
192 DATA ((CALC(I,K,4 ),I=1,8), K=1,9) /
193 + 0.0000E+00, 0.0000E+00, 0.1625E-02, -0.2514E-03,
194 + 0.1122E-02, -0.2514E-03, 0.6500E-02, 0.7358E-02,
195 + 0.0000E+00, 0.0000E+00, 0.1602E-02, -0.2260E-03,
196 + 0.1150E-02, -0.2260E-03, 0.6948E-02, 0.8685E-02,
197 + 0.0000E+00, 0.0000E+00, 0.1578E-02, -0.1988E-03,
198 + 0.1180E-02, -0.1988E-03, 0.7432E-02, 0.1006E-01,
199 + 0.0000E+00, 0.0000E+00, 0.1551E-02, -0.1696E-03,
200 + 0.1212E-02, -0.1696E-03, 0.7953E-02, 0.1148E-01,
201 + 0.0000E+00, 0.0000E+00, 0.1522E-02, -0.1382E-03,
202 + 0.1245E-02, -0.1382E-03, 0.8512E-02, 0.1296E-01,
203 + 0.0000E+00, 0.0000E+00, 0.1490E-02, -0.1045E-03,
204 + 0.1281E-02, -0.1045E-03, 0.9113E-02, 0.1450E-01,
205 + 0.0000E+00, 0.0000E+00, 0.1456E-02, -0.6824E-04,
206 + 0.1319E-02, -0.6824E-04, 0.9757E-02, 0.1610E-01,
207 + 0.0000E+00, 0.0000E+00, 0.1419E-02, -0.2924E-04,
208 + 0.1360E-02, -0.2924E-04, 0.1045E-01, 0.1779E-01,
209 + 0.0000E+00, 0.0000E+00, 0.1378E-02, 0.1275E-04,
210 + 0.1403E-02, 0.1275E-04, 0.1119E-01, 0.1957E-01/
211 DATA ((CALC(I,K,4 ),I=1,8), K=10,18) /
212 + 0.0000E+00, 0.0000E+00, 0.1334E-02, 0.5797E-04,
213 + 0.1450E-02, 0.5797E-04, 0.1198E-01, 0.2145E-01,
214 + 0.0000E+00, 0.0000E+00, 0.1287E-02, 0.1067E-03,
215 + 0.1499E-02, 0.1067E-03, 0.1284E-01, 0.2345E-01,
216 + 0.0000E+00, 0.0000E+00, 0.1235E-02, 0.1592E-03,
217 + 0.1553E-02, 0.1592E-03, 0.1375E-01, 0.2557E-01,
218 + 0.0000E+00, 0.0000E+00, 0.1179E-02, 0.2158E-03,
219 + 0.1610E-02, 0.2158E-03, 0.1474E-01, 0.2784E-01,
220 + 0.0000E+00, 0.0000E+00, 0.1118E-02, 0.2769E-03,
221 + 0.1671E-02, 0.2769E-03, 0.1579E-01, 0.3028E-01,
222 + 0.0000E+00, 0.0000E+00, 0.1052E-02, 0.3428E-03,
223 + 0.1737E-02, 0.3428E-03, 0.1693E-01, 0.3289E-01,
224 + 0.0000E+00, 0.0000E+00, 0.9810E-03, 0.4140E-03,
225 + 0.1808E-02, 0.4140E-03, 0.1815E-01, 0.3571E-01,
226 + 0.0000E+00, 0.0000E+00, 0.9038E-03, 0.4909E-03,
227 + 0.1885E-02, 0.4909E-03, 0.1947E-01, 0.3875E-01,
228 + 0.0000E+00, 0.0000E+00, 0.8201E-03, 0.5739E-03,
229 + 0.1967E-02, 0.5739E-03, 0.2090E-01, 0.4204E-01/
230 DATA ((CALC(I,K,4 ),I=1,8), K=19,20) /
231 + 0.0000E+00, 0.0000E+00, 0.7295E-03, 0.6636E-03,
232 + 0.2056E-02, 0.6636E-03, 0.2243E-01, 0.4560E-01,
233 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
234 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
235 DATA ((CALC(I,K,5 ),I=1,8), K=1,9) /
236 + 0.0000E+00, 0.0000E+00, 0.1572E-02, -0.1542E-03,
237 + 0.1264E-02, -0.1542E-03, 0.6288E-02, 0.6899E-02,
238 + 0.0000E+00, 0.0000E+00, 0.1541E-02, -0.1209E-03,
239 + 0.1299E-02, -0.1209E-03, 0.6695E-02, 0.7867E-02,
240 + 0.0000E+00, 0.0000E+00, 0.1507E-02, -0.8520E-04,
241 + 0.1336E-02, -0.8520E-04, 0.7132E-02, 0.8857E-02,
242 + 0.0000E+00, 0.0000E+00, 0.1470E-02, -0.4684E-04,
243 + 0.1375E-02, -0.4684E-04, 0.7599E-02, 0.9865E-02,
244 + 0.0000E+00, 0.0000E+00, 0.1432E-02, -0.5605E-05,
245 + 0.1418E-02, -0.5605E-05, 0.8099E-02, 0.1091E-01,
246 + 0.0000E+00, 0.0000E+00, 0.1388E-02, 0.3874E-04,
247 + 0.1463E-02, 0.3874E-04, 0.8635E-02, 0.1199E-01,
248 + 0.0000E+00, 0.0000E+00, 0.1341E-02, 0.8647E-04,
249 + 0.1511E-02, 0.8647E-04, 0.9210E-02, 0.1311E-01,
250 + 0.0000E+00, 0.0000E+00, 0.1289E-02, 0.1379E-03,
251 + 0.1563E-02, 0.1379E-03, 0.9825E-02, 0.1429E-01,
252 + 0.0000E+00, 0.0000E+00, 0.1234E-02, 0.1932E-03,
253 + 0.1618E-02, 0.1932E-03, 0.1049E-01, 0.1553E-01/
254 DATA ((CALC(I,K,5 ),I=1,8), K=10,18) /
255 + 0.0000E+00, 0.0000E+00, 0.1174E-02, 0.2529E-03,
256 + 0.1677E-02, 0.2529E-03, 0.1119E-01, 0.1684E-01,
257 + 0.0000E+00, 0.0000E+00, 0.1109E-02, 0.3173E-03,
258 + 0.1741E-02, 0.3173E-03, 0.1196E-01, 0.1824E-01,
259 + 0.0000E+00, 0.0000E+00, 0.1038E-02, 0.3867E-03,
260 + 0.1810E-02, 0.3867E-03, 0.1277E-01, 0.1973E-01,
261 + 0.0000E+00, 0.0000E+00, 0.9620E-03, 0.4616E-03,
262 + 0.1883E-02, 0.4616E-03, 0.1366E-01, 0.2133E-01,
263 + 0.0000E+00, 0.0000E+00, 0.8794E-03, 0.5425E-03,
264 + 0.1963E-02, 0.5425E-03, 0.1461E-01, 0.2305E-01,
265 + 0.0000E+00, 0.0000E+00, 0.7902E-03, 0.6299E-03,
266 + 0.2048E-02, 0.6299E-03, 0.1563E-01, 0.2492E-01,
267 + 0.0000E+00, 0.0000E+00, 0.6935E-03, 0.7243E-03,
268 + 0.2141E-02, 0.7243E-03, 0.1674E-01, 0.2694E-01,
269 + 0.0000E+00, 0.0000E+00, 0.5890E-03, 0.8263E-03,
270 + 0.2240E-02, 0.8263E-03, 0.1794E-01, 0.2912E-01,
271 + 0.0000E+00, 0.0000E+00, 0.4759E-03, 0.9366E-03,
272 + 0.2348E-02, 0.9366E-03, 0.1924E-01, 0.3150E-01/
273 DATA ((CALC(I,K,5 ),I=1,8), K=19,20) /
274 + 0.0000E+00, 0.0000E+00, 0.3534E-03, 0.1056E-02,
275 + 0.2464E-02, 0.1056E-02, 0.2064E-01, 0.3408E-01,
276 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
277 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
278 DATA ((CALC(I,K,6 ),I=1,8), K=1,9) /
279 + 0.0000E+00, 0.0000E+00, 0.1509E-02, -0.5692E-04,
280 + 0.1395E-02, -0.5692E-04, 0.6035E-02, 0.6441E-02,
281 + 0.0000E+00, 0.0000E+00, 0.1469E-02, -0.1625E-04,
282 + 0.1436E-02, -0.1625E-04, 0.6426E-02, 0.7136E-02,
283 + 0.0000E+00, 0.0000E+00, 0.1425E-02, 0.2742E-04,
284 + 0.1480E-02, 0.2742E-04, 0.6844E-02, 0.7855E-02,
285 + 0.0000E+00, 0.0000E+00, 0.1378E-02, 0.7439E-04,
286 + 0.1527E-02, 0.7439E-04, 0.7291E-02, 0.8574E-02,
287 + 0.0000E+00, 0.0000E+00, 0.1328E-02, 0.1249E-03,
288 + 0.1577E-02, 0.1249E-03, 0.7770E-02, 0.9307E-02,
289 + 0.0000E+00, 0.0000E+00, 0.1272E-02, 0.1793E-03,
290 + 0.1631E-02, 0.1793E-03, 0.8285E-02, 0.1009E-01,
291 + 0.0000E+00, 0.0000E+00, 0.1213E-02, 0.2380E-03,
292 + 0.1689E-02, 0.2380E-03, 0.8836E-02, 0.1087E-01,
293 + 0.0000E+00, 0.0000E+00, 0.1148E-02, 0.3012E-03,
294 + 0.1751E-02, 0.3012E-03, 0.9427E-02, 0.1171E-01,
295 + 0.0000E+00, 0.0000E+00, 0.1079E-02, 0.3694E-03,
296 + 0.1817E-02, 0.3694E-03, 0.1006E-01, 0.1259E-01/
297 DATA ((CALC(I,K,6 ),I=1,8), K=10,18) /
298 + 0.0000E+00, 0.0000E+00, 0.1003E-02, 0.4429E-03,
299 + 0.1889E-02, 0.4429E-03, 0.1075E-01, 0.1353E-01,
300 + 0.0000E+00, 0.0000E+00, 0.9214E-03, 0.5223E-03,
301 + 0.1966E-02, 0.5223E-03, 0.1149E-01, 0.1454E-01,
302 + 0.0000E+00, 0.0000E+00, 0.8329E-03, 0.6081E-03,
303 + 0.2049E-02, 0.6081E-03, 0.1229E-01, 0.1563E-01,
304 + 0.0000E+00, 0.0000E+00, 0.7372E-03, 0.7007E-03,
305 + 0.2139E-02, 0.7007E-03, 0.1315E-01, 0.1681E-01,
306 + 0.0000E+00, 0.0000E+00, 0.6337E-03, 0.8008E-03,
307 + 0.2235E-02, 0.8008E-03, 0.1408E-01, 0.1809E-01,
308 + 0.0000E+00, 0.0000E+00, 0.5217E-03, 0.9090E-03,
309 + 0.2340E-02, 0.9090E-03, 0.1509E-01, 0.1948E-01,
310 + 0.0000E+00, 0.0000E+00, 0.4005E-03, 0.1026E-02,
311 + 0.2453E-02, 0.1026E-02, 0.1618E-01, 0.2100E-01,
312 + 0.0000E+00, 0.0000E+00, 0.2694E-03, 0.1152E-02,
313 + 0.2574E-02, 0.1152E-02, 0.1737E-01, 0.2265E-01,
314 + 0.0000E+00, 0.0000E+00, 0.1276E-03, 0.1289E-02,
315 + 0.2706E-02, 0.1289E-02, 0.1865E-01, 0.2446E-01/
316 DATA ((CALC(I,K,6 ),I=1,8), K=19,20) /
317 + 0.0000E+00, 0.0000E+00, -0.2581E-04, 0.1437E-02,
318 + 0.2849E-02, 0.1437E-02, 0.2005E-01, 0.2643E-01,
319 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
320 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
321 DATA ((CALC(I,K,7 ),I=1,8), K=1,9) /
322 + 0.0000E+00, 0.0000E+00, 0.1444E-02, 0.3563E-04,
323 + 0.1515E-02, 0.3563E-04, 0.5775E-02, 0.5978E-02,
324 + 0.0000E+00, 0.0000E+00, 0.1396E-02, 0.8289E-04,
325 + 0.1562E-02, 0.8289E-04, 0.6166E-02, 0.6468E-02,
326 + 0.0000E+00, 0.0000E+00, 0.1344E-02, 0.1337E-03,
327 + 0.1612E-02, 0.1337E-03, 0.6585E-02, 0.6960E-02,
328 + 0.0000E+00, 0.0000E+00, 0.1288E-02, 0.1885E-03,
329 + 0.1665E-02, 0.1885E-03, 0.7032E-02, 0.7462E-02,
330 + 0.0000E+00, 0.0000E+00, 0.1227E-02, 0.2475E-03,
331 + 0.1722E-02, 0.2475E-03, 0.7513E-02, 0.7978E-02,
332 + 0.0000E+00, 0.0000E+00, 0.1162E-02, 0.3111E-03,
333 + 0.1784E-02, 0.3111E-03, 0.8029E-02, 0.8515E-02,
334 + 0.0000E+00, 0.0000E+00, 0.1090E-02, 0.3797E-03,
335 + 0.1850E-02, 0.3797E-03, 0.8584E-02, 0.9080E-02,
336 + 0.0000E+00, 0.0000E+00, 0.1014E-02, 0.4538E-03,
337 + 0.1921E-02, 0.4538E-03, 0.9182E-02, 0.9681E-02,
338 + 0.0000E+00, 0.0000E+00, 0.9302E-03, 0.5338E-03,
339 + 0.1998E-02, 0.5338E-03, 0.9828E-02, 0.1032E-01/
340 DATA ((CALC(I,K,7 ),I=1,8), K=10,18) /
341 + 0.0000E+00, 0.0000E+00, 0.8400E-03, 0.6202E-03,
342 + 0.2080E-02, 0.6202E-03, 0.1053E-01, 0.1102E-01,
343 + 0.0000E+00, 0.0000E+00, 0.7424E-03, 0.7136E-03,
344 + 0.2170E-02, 0.7136E-03, 0.1128E-01, 0.1177E-01,
345 + 0.0000E+00, 0.0000E+00, 0.6368E-03, 0.8145E-03,
346 + 0.2266E-02, 0.8145E-03, 0.1210E-01, 0.1260E-01,
347 + 0.0000E+00, 0.0000E+00, 0.5226E-03, 0.9237E-03,
348 + 0.2370E-02, 0.9237E-03, 0.1298E-01, 0.1349E-01,
349 + 0.0000E+00, 0.0000E+00, 0.3990E-03, 0.1042E-02,
350 + 0.2483E-02, 0.1042E-02, 0.1394E-01, 0.1448E-01,
351 + 0.0000E+00, 0.0000E+00, 0.2653E-03, 0.1170E-02,
352 + 0.2605E-02, 0.1170E-02, 0.1498E-01, 0.1555E-01,
353 + 0.0000E+00, 0.0000E+00, 0.1206E-03, 0.1308E-02,
354 + 0.2736E-02, 0.1308E-02, 0.1611E-01, 0.1673E-01,
355 + 0.0000E+00, 0.0000E+00, -0.3584E-04, 0.1457E-02,
356 + 0.2879E-02, 0.1457E-02, 0.1734E-01, 0.1803E-01,
357 + 0.0000E+00, 0.0000E+00, -0.2051E-03, 0.1619E-02,
358 + 0.3034E-02, 0.1619E-02, 0.1867E-01, 0.1945E-01/
359 DATA ((CALC(I,K,7 ),I=1,8), K=19,20) /
360 + 0.0000E+00, 0.0000E+00, -0.3882E-03, 0.1795E-02,
361 + 0.3201E-02, 0.1795E-02, 0.2012E-01, 0.2101E-01,
362 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
363 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
364 DATA ((CALC(I,K,8 ),I=1,8), K=1,9) /
365 + 0.0000E+00, 0.0000E+00, 0.1383E-02, 0.1197E-03,
366 + 0.1623E-02, 0.1197E-03, 0.5533E-02, 0.5504E-02,
367 + 0.0000E+00, 0.0000E+00, 0.1328E-02, 0.1726E-03,
368 + 0.1674E-02, 0.1726E-03, 0.5934E-02, 0.5833E-02,
369 + 0.0000E+00, 0.0000E+00, 0.1269E-02, 0.2296E-03,
370 + 0.1728E-02, 0.2296E-03, 0.6363E-02, 0.6164E-02,
371 + 0.0000E+00, 0.0000E+00, 0.1205E-02, 0.2912E-03,
372 + 0.1787E-02, 0.2912E-03, 0.6823E-02, 0.6502E-02,
373 + 0.0000E+00, 0.0000E+00, 0.1135E-02, 0.3576E-03,
374 + 0.1850E-02, 0.3576E-03, 0.7318E-02, 0.6855E-02,
375 + 0.0000E+00, 0.0000E+00, 0.1059E-02, 0.4293E-03,
376 + 0.1918E-02, 0.4293E-03, 0.7851E-02, 0.7229E-02,
377 + 0.0000E+00, 0.0000E+00, 0.9777E-03, 0.5069E-03,
378 + 0.1991E-02, 0.5069E-03, 0.8426E-02, 0.7628E-02,
379 + 0.0000E+00, 0.0000E+00, 0.8891E-03, 0.5907E-03,
380 + 0.2070E-02, 0.5907E-03, 0.9048E-02, 0.8061E-02,
381 + 0.0000E+00, 0.0000E+00, 0.7933E-03, 0.6813E-03,
382 + 0.2156E-02, 0.6813E-03, 0.9720E-02, 0.8533E-02/
383 DATA ((CALC(I,K,8 ),I=1,8), K=10,18) /
384 + 0.0000E+00, 0.0000E+00, 0.6895E-03, 0.7794E-03,
385 + 0.2248E-02, 0.7794E-03, 0.1045E-01, 0.9051E-02,
386 + 0.0000E+00, 0.0000E+00, 0.5772E-03, 0.8855E-03,
387 + 0.2348E-02, 0.8855E-03, 0.1124E-01, 0.9622E-02,
388 + 0.0000E+00, 0.0000E+00, 0.4557E-03, 0.1000E-02,
389 + 0.2456E-02, 0.1000E-02, 0.1209E-01, 0.1025E-01,
390 + 0.0000E+00, 0.0000E+00, 0.3242E-03, 0.1125E-02,
391 + 0.2573E-02, 0.1125E-02, 0.1302E-01, 0.1095E-01,
392 + 0.0000E+00, 0.0000E+00, 0.1819E-03, 0.1259E-02,
393 + 0.2700E-02, 0.1259E-02, 0.1403E-01, 0.1172E-01,
394 + 0.0000E+00, 0.0000E+00, 0.2794E-04, 0.1405E-02,
395 + 0.2838E-02, 0.1405E-02, 0.1513E-01, 0.1257E-01,
396 + 0.0000E+00, 0.0000E+00, -0.1386E-03, 0.1563E-02,
397 + 0.2987E-02, 0.1563E-02, 0.1632E-01, 0.1351E-01,
398 + 0.0000E+00, 0.0000E+00, -0.3188E-03, 0.1734E-02,
399 + 0.3149E-02, 0.1734E-02, 0.1762E-01, 0.1455E-01,
400 + 0.0000E+00, 0.0000E+00, -0.5136E-03, 0.1919E-02,
401 + 0.3324E-02, 0.1919E-02, 0.1903E-01, 0.1569E-01/
402 DATA ((CALC(I,K,8 ),I=1,8), K=19,20) /
403 + 0.0000E+00, 0.0000E+00, -0.7244E-03, 0.2119E-02,
404 + 0.3514E-02, 0.2119E-02, 0.2056E-01, 0.1695E-01,
405 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
406 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
407 DATA ((CALC(I,K,9 ),I=1,8), K=1,9) /
408 + 0.0000E+00, 0.0000E+00, 0.1331E-02, 0.1925E-03,
409 + 0.1716E-02, 0.1925E-03, 0.5326E-02, 0.5026E-02,
410 + 0.0000E+00, 0.0000E+00, 0.1270E-02, 0.2501E-03,
411 + 0.1771E-02, 0.2501E-03, 0.5741E-02, 0.5228E-02,
412 + 0.0000E+00, 0.0000E+00, 0.1204E-02, 0.3123E-03,
413 + 0.1829E-02, 0.3123E-03, 0.6186E-02, 0.5433E-02,
414 + 0.0000E+00, 0.0000E+00, 0.1132E-02, 0.3796E-03,
415 + 0.1891E-02, 0.3796E-03, 0.6664E-02, 0.5649E-02,
416 + 0.0000E+00, 0.0000E+00, 0.1055E-02, 0.4523E-03,
417 + 0.1959E-02, 0.4523E-03, 0.7180E-02, 0.5880E-02,
418 + 0.0000E+00, 0.0000E+00, 0.9702E-03, 0.5311E-03,
419 + 0.2032E-02, 0.5311E-03, 0.7737E-02, 0.6131E-02,
420 + 0.0000E+00, 0.0000E+00, 0.8788E-03, 0.6164E-03,
421 + 0.2112E-02, 0.6164E-03, 0.8339E-02, 0.6410E-02,
422 + 0.0000E+00, 0.0000E+00, 0.7797E-03, 0.7087E-03,
423 + 0.2197E-02, 0.7087E-03, 0.8992E-02, 0.6720E-02,
424 + 0.0000E+00, 0.0000E+00, 0.6724E-03, 0.8088E-03,
425 + 0.2290E-02, 0.8088E-03, 0.9699E-02, 0.7067E-02/
426 DATA ((CALC(I,K,9 ),I=1,8), K=10,18) /
427 + 0.0000E+00, 0.0000E+00, 0.5562E-03, 0.9172E-03,
428 + 0.2391E-02, 0.9172E-03, 0.1047E-01, 0.7458E-02,
429 + 0.0000E+00, 0.0000E+00, 0.4303E-03, 0.1035E-02,
430 + 0.2500E-02, 0.1035E-02, 0.1130E-01, 0.7896E-02,
431 + 0.0000E+00, 0.0000E+00, 0.2941E-03, 0.1162E-02,
432 + 0.2618E-02, 0.1162E-02, 0.1220E-01, 0.8388E-02,
433 + 0.0000E+00, 0.0000E+00, 0.1467E-03, 0.1300E-02,
434 + 0.2747E-02, 0.1300E-02, 0.1319E-01, 0.8939E-02,
435 + 0.0000E+00, 0.0000E+00, -0.1287E-04, 0.1449E-02,
436 + 0.2886E-02, 0.1449E-02, 0.1426E-01, 0.9555E-02,
437 + 0.0000E+00, 0.0000E+00, -0.1855E-03, 0.1611E-02,
438 + 0.3037E-02, 0.1611E-02, 0.1542E-01, 0.1024E-01,
439 + 0.0000E+00, 0.0000E+00, -0.3723E-03, 0.1787E-02,
440 + 0.3202E-02, 0.1787E-02, 0.1669E-01, 0.1100E-01,
441 + 0.0000E+00, 0.0000E+00, -0.5743E-03, 0.1977E-02,
442 + 0.3380E-02, 0.1977E-02, 0.1806E-01, 0.1185E-01,
443 + 0.0000E+00, 0.0000E+00, -0.7929E-03, 0.2183E-02,
444 + 0.3574E-02, 0.2183E-02, 0.1956E-01, 0.1278E-01/
445 DATA ((CALC(I,K,9 ),I=1,8), K=19,20) /
446 + 0.0000E+00, 0.0000E+00, -0.1029E-02, 0.2407E-02,
447 + 0.3785E-02, 0.2407E-02, 0.2119E-01, 0.1381E-01,
448 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
449 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
450 DATA ((CALC(I,K,10),I=1,8), K=1,9) /
451 + 0.0000E+00, 0.0000E+00, 0.1290E-02, 0.2527E-03,
452 + 0.1796E-02, 0.2527E-03, 0.5161E-02, 0.4547E-02,
453 + 0.0000E+00, 0.0000E+00, 0.1224E-02, 0.3140E-03,
454 + 0.1852E-02, 0.3140E-03, 0.5591E-02, 0.4649E-02,
455 + 0.0000E+00, 0.0000E+00, 0.1151E-02, 0.3804E-03,
456 + 0.1912E-02, 0.3804E-03, 0.6053E-02, 0.4759E-02,
457 + 0.0000E+00, 0.0000E+00, 0.1073E-02, 0.4525E-03,
458 + 0.1978E-02, 0.4525E-03, 0.6551E-02, 0.4881E-02,
459 + 0.0000E+00, 0.0000E+00, 0.9879E-03, 0.5306E-03,
460 + 0.2049E-02, 0.5306E-03, 0.7090E-02, 0.5021E-02,
461 + 0.0000E+00, 0.0000E+00, 0.8954E-03, 0.6154E-03,
462 + 0.2126E-02, 0.6154E-03, 0.7673E-02, 0.5183E-02,
463 + 0.0000E+00, 0.0000E+00, 0.7952E-03, 0.7074E-03,
464 + 0.2210E-02, 0.7074E-03, 0.8304E-02, 0.5372E-02,
465 + 0.0000E+00, 0.0000E+00, 0.6865E-03, 0.8072E-03,
466 + 0.2301E-02, 0.8072E-03, 0.8990E-02, 0.5593E-02,
467 + 0.0000E+00, 0.0000E+00, 0.5686E-03, 0.9156E-03,
468 + 0.2400E-02, 0.9156E-03, 0.9734E-02, 0.5850E-02/
469 DATA ((CALC(I,K,10),I=1,8), K=10,18) /
470 + 0.0000E+00, 0.0000E+00, 0.4410E-03, 0.1033E-02,
471 + 0.2507E-02, 0.1033E-02, 0.1054E-01, 0.6146E-02,
472 + 0.0000E+00, 0.0000E+00, 0.3027E-03, 0.1161E-02,
473 + 0.2624E-02, 0.1161E-02, 0.1142E-01, 0.6487E-02,
474 + 0.0000E+00, 0.0000E+00, 0.1529E-03, 0.1299E-02,
475 + 0.2751E-02, 0.1299E-02, 0.1238E-01, 0.6877E-02,
476 + 0.0000E+00, 0.0000E+00, -0.9250E-05, 0.1449E-02,
477 + 0.2890E-02, 0.1449E-02, 0.1342E-01, 0.7320E-02,
478 + 0.0000E+00, 0.0000E+00, -0.1847E-03, 0.1612E-02,
479 + 0.3040E-02, 0.1612E-02, 0.1456E-01, 0.7820E-02,
480 + 0.0000E+00, 0.0000E+00, -0.3747E-03, 0.1789E-02,
481 + 0.3204E-02, 0.1789E-02, 0.1579E-01, 0.8381E-02,
482 + 0.0000E+00, 0.0000E+00, -0.5801E-03, 0.1981E-02,
483 + 0.3382E-02, 0.1981E-02, 0.1713E-01, 0.9008E-02,
484 + 0.0000E+00, 0.0000E+00, -0.8024E-03, 0.2189E-02,
485 + 0.3575E-02, 0.2189E-02, 0.1859E-01, 0.9705E-02,
486 + 0.0000E+00, 0.0000E+00, -0.1043E-02, 0.2414E-02,
487 + 0.3786E-02, 0.2414E-02, 0.2018E-01, 0.1048E-01/
488 DATA ((CALC(I,K,10),I=1,8), K=19,20) /
489 + 0.0000E+00, 0.0000E+00, -0.1303E-02, 0.2659E-02,
490 + 0.4015E-02, 0.2659E-02, 0.2191E-01, 0.1133E-01,
491 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
492 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
493 DATA ((CALC(I,K,11),I=1,8), K=1,9) /
494 + 0.0000E+00, 0.0000E+00, 0.1259E-02, 0.2998E-03,
495 + 0.1859E-02, 0.2998E-03, 0.5037E-02, 0.4072E-02,
496 + 0.0000E+00, 0.0000E+00, 0.1188E-02, 0.3641E-03,
497 + 0.1916E-02, 0.3641E-03, 0.5481E-02, 0.4098E-02,
498 + 0.0000E+00, 0.0000E+00, 0.1110E-02, 0.4339E-03,
499 + 0.1978E-02, 0.4339E-03, 0.5959E-02, 0.4135E-02,
500 + 0.0000E+00, 0.0000E+00, 0.1026E-02, 0.5100E-03,
501 + 0.2046E-02, 0.5100E-03, 0.6476E-02, 0.4188E-02,
502 + 0.0000E+00, 0.0000E+00, 0.9340E-03, 0.5927E-03,
503 + 0.2119E-02, 0.5927E-03, 0.7036E-02, 0.4261E-02,
504 + 0.0000E+00, 0.0000E+00, 0.8342E-03, 0.6827E-03,
505 + 0.2200E-02, 0.6827E-03, 0.7645E-02, 0.4358E-02,
506 + 0.0000E+00, 0.0000E+00, 0.7257E-03, 0.7807E-03,
507 + 0.2287E-02, 0.7807E-03, 0.8305E-02, 0.4483E-02,
508 + 0.0000E+00, 0.0000E+00, 0.6081E-03, 0.8871E-03,
509 + 0.2382E-02, 0.8871E-03, 0.9024E-02, 0.4638E-02,
510 + 0.0000E+00, 0.0000E+00, 0.4806E-03, 0.1003E-02,
511 + 0.2486E-02, 0.1003E-02, 0.9807E-02, 0.4829E-02/
512 DATA ((CALC(I,K,11),I=1,8), K=10,18) /
513 + 0.0000E+00, 0.0000E+00, 0.3422E-03, 0.1129E-02,
514 + 0.2600E-02, 0.1129E-02, 0.1066E-01, 0.5057E-02,
515 + 0.0000E+00, 0.0000E+00, 0.1923E-03, 0.1266E-02,
516 + 0.2724E-02, 0.1266E-02, 0.1159E-01, 0.5327E-02,
517 + 0.0000E+00, 0.0000E+00, 0.2992E-04, 0.1414E-02,
518 + 0.2858E-02, 0.1414E-02, 0.1260E-01, 0.5641E-02,
519 + 0.0000E+00, 0.0000E+00, -0.1460E-03, 0.1575E-02,
520 + 0.3005E-02, 0.1575E-02, 0.1370E-01, 0.6002E-02,
521 + 0.0000E+00, 0.0000E+00, -0.3363E-03, 0.1751E-02,
522 + 0.3166E-02, 0.1751E-02, 0.1489E-01, 0.6414E-02,
523 + 0.0000E+00, 0.0000E+00, -0.5424E-03, 0.1941E-02,
524 + 0.3340E-02, 0.1941E-02, 0.1620E-01, 0.6879E-02,
525 + 0.0000E+00, 0.0000E+00, -0.7653E-03, 0.2148E-02,
526 + 0.3530E-02, 0.2148E-02, 0.1761E-01, 0.7397E-02,
527 + 0.0000E+00, 0.0000E+00, -0.1006E-02, 0.2372E-02,
528 + 0.3738E-02, 0.2372E-02, 0.1916E-01, 0.7983E-02,
529 + 0.0000E+00, 0.0000E+00, -0.1267E-02, 0.2615E-02,
530 + 0.3963E-02, 0.2615E-02, 0.2084E-01, 0.8632E-02/
531 DATA ((CALC(I,K,11),I=1,8), K=19,20) /
532 + 0.0000E+00, 0.0000E+00, -0.1549E-02, 0.2879E-02,
533 + 0.4209E-02, 0.2879E-02, 0.2266E-01, 0.9347E-02,
534 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
535 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
536 DATA ((CALC(I,K,12),I=1,8), K=1,9) /
537 + 0.0000E+00, 0.0000E+00, 0.1238E-02, 0.3350E-03,
538 + 0.1908E-02, 0.3350E-03, 0.4950E-02, 0.3607E-02,
539 + 0.0000E+00, 0.0000E+00, 0.1162E-02, 0.4016E-03,
540 + 0.1965E-02, 0.4016E-03, 0.5405E-02, 0.3576E-02,
541 + 0.0000E+00, 0.0000E+00, 0.1079E-02, 0.4744E-03,
542 + 0.2028E-02, 0.4744E-03, 0.5897E-02, 0.3560E-02,
543 + 0.0000E+00, 0.0000E+00, 0.9892E-03, 0.5539E-03,
544 + 0.2097E-02, 0.5539E-03, 0.6431E-02, 0.3563E-02,
545 + 0.0000E+00, 0.0000E+00, 0.8911E-03, 0.6407E-03,
546 + 0.2173E-02, 0.6407E-03, 0.7012E-02, 0.3588E-02,
547 + 0.0000E+00, 0.0000E+00, 0.7844E-03, 0.7355E-03,
548 + 0.2255E-02, 0.7355E-03, 0.7645E-02, 0.3638E-02,
549 + 0.0000E+00, 0.0000E+00, 0.6683E-03, 0.8388E-03,
550 + 0.2346E-02, 0.8388E-03, 0.8335E-02, 0.3717E-02,
551 + 0.0000E+00, 0.0000E+00, 0.5422E-03, 0.9514E-03,
552 + 0.2445E-02, 0.9514E-03, 0.9087E-02, 0.3826E-02,
553 + 0.0000E+00, 0.0000E+00, 0.4053E-03, 0.1074E-02,
554 + 0.2554E-02, 0.1074E-02, 0.9907E-02, 0.3969E-02/
555 DATA ((CALC(I,K,12),I=1,8), K=10,18) /
556 + 0.0000E+00, 0.0000E+00, 0.2568E-03, 0.1208E-02,
557 + 0.2672E-02, 0.1208E-02, 0.1080E-01, 0.4147E-02,
558 + 0.0000E+00, 0.0000E+00, 0.9573E-04, 0.1353E-02,
559 + 0.2802E-02, 0.1353E-02, 0.1178E-01, 0.4364E-02,
560 + 0.0000E+00, 0.0000E+00, -0.7882E-04, 0.1511E-02,
561 + 0.2944E-02, 0.1511E-02, 0.1284E-01, 0.4621E-02,
562 + 0.0000E+00, 0.0000E+00, -0.2679E-03, 0.1683E-02,
563 + 0.3099E-02, 0.1683E-02, 0.1400E-01, 0.4920E-02,
564 + 0.0000E+00, 0.0000E+00, -0.4726E-03, 0.1871E-02,
565 + 0.3269E-02, 0.1871E-02, 0.1526E-01, 0.5264E-02,
566 + 0.0000E+00, 0.0000E+00, -0.6941E-03, 0.2074E-02,
567 + 0.3454E-02, 0.2074E-02, 0.1663E-01, 0.5654E-02,
568 + 0.0000E+00, 0.0000E+00, -0.9338E-03, 0.2294E-02,
569 + 0.3656E-02, 0.2294E-02, 0.1813E-01, 0.6094E-02,
570 + 0.0000E+00, 0.0000E+00, -0.1193E-02, 0.2534E-02,
571 + 0.3876E-02, 0.2534E-02, 0.1976E-01, 0.6586E-02,
572 + 0.0000E+00, 0.0000E+00, -0.1473E-02, 0.2794E-02,
573 + 0.4116E-02, 0.2794E-02, 0.2153E-01, 0.7132E-02/
574 DATA ((CALC(I,K,12),I=1,8), K=19,20) /
575 + 0.0000E+00, 0.0000E+00, -0.1776E-02, 0.3076E-02,
576 + 0.4377E-02, 0.3076E-02, 0.2345E-01, 0.7735E-02,
577 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
578 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
579 DATA ((CALC(I,K,13),I=1,8), K=1,9) /
580 + 0.0000E+00, 0.0000E+00, 0.1222E-02, 0.3601E-03,
581 + 0.1942E-02, 0.3601E-03, 0.4890E-02, 0.3159E-02,
582 + 0.0000E+00, 0.0000E+00, 0.1143E-02, 0.4288E-03,
583 + 0.2000E-02, 0.4288E-03, 0.5354E-02, 0.3087E-02,
584 + 0.0000E+00, 0.0000E+00, 0.1056E-02, 0.5042E-03,
585 + 0.2064E-02, 0.5042E-03, 0.5859E-02, 0.3033E-02,
586 + 0.0000E+00, 0.0000E+00, 0.9603E-03, 0.5869E-03,
587 + 0.2134E-02, 0.5869E-03, 0.6411E-02, 0.3001E-02,
588 + 0.0000E+00, 0.0000E+00, 0.8562E-03, 0.6777E-03,
589 + 0.2211E-02, 0.6777E-03, 0.7013E-02, 0.2993E-02,
590 + 0.0000E+00, 0.0000E+00, 0.7427E-03, 0.7770E-03,
591 + 0.2297E-02, 0.7770E-03, 0.7672E-02, 0.3011E-02,
592 + 0.0000E+00, 0.0000E+00, 0.6191E-03, 0.8857E-03,
593 + 0.2391E-02, 0.8857E-03, 0.8391E-02, 0.3058E-02,
594 + 0.0000E+00, 0.0000E+00, 0.4847E-03, 0.1004E-02,
595 + 0.2494E-02, 0.1004E-02, 0.9178E-02, 0.3135E-02,
596 + 0.0000E+00, 0.0000E+00, 0.3385E-03, 0.1134E-02,
597 + 0.2607E-02, 0.1134E-02, 0.1004E-01, 0.3244E-02/
598 DATA ((CALC(I,K,13),I=1,8), K=10,18) /
599 + 0.0000E+00, 0.0000E+00, 0.1798E-03, 0.1275E-02,
600 + 0.2731E-02, 0.1275E-02, 0.1098E-01, 0.3387E-02,
601 + 0.0000E+00, 0.0000E+00, 0.7585E-05, 0.1430E-02,
602 + 0.2867E-02, 0.1430E-02, 0.1200E-01, 0.3564E-02,
603 + 0.0000E+00, 0.0000E+00, -0.1791E-03, 0.1598E-02,
604 + 0.3016E-02, 0.1598E-02, 0.1312E-01, 0.3777E-02,
605 + 0.0000E+00, 0.0000E+00, -0.3814E-03, 0.1780E-02,
606 + 0.3180E-02, 0.1780E-02, 0.1434E-01, 0.4029E-02,
607 + 0.0000E+00, 0.0000E+00, -0.6005E-03, 0.1979E-02,
608 + 0.3358E-02, 0.1979E-02, 0.1567E-01, 0.4319E-02,
609 + 0.0000E+00, 0.0000E+00, -0.8375E-03, 0.2195E-02,
610 + 0.3554E-02, 0.2195E-02, 0.1712E-01, 0.4650E-02,
611 + 0.0000E+00, 0.0000E+00, -0.1094E-02, 0.2430E-02,
612 + 0.3767E-02, 0.2430E-02, 0.1869E-01, 0.5023E-02,
613 + 0.0000E+00, 0.0000E+00, -0.1371E-02, 0.2685E-02,
614 + 0.4000E-02, 0.2685E-02, 0.2041E-01, 0.5440E-02,
615 + 0.0000E+00, 0.0000E+00, -0.1671E-02, 0.2962E-02,
616 + 0.4254E-02, 0.2962E-02, 0.2227E-01, 0.5903E-02/
617 DATA ((CALC(I,K,13),I=1,8), K=19,20) /
618 + 0.0000E+00, 0.0000E+00, -0.1994E-02, 0.3262E-02,
619 + 0.4531E-02, 0.3262E-02, 0.2430E-01, 0.6414E-02,
620 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
621 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
622 DATA ((CALC(I,K,14),I=1,8), K=1,9) /
623 + 0.0000E+00, 0.0000E+00, 0.1213E-02, 0.3767E-03,
624 + 0.1967E-02, 0.3767E-03, 0.4853E-02, 0.2733E-02,
625 + 0.0000E+00, 0.0000E+00, 0.1130E-02, 0.4476E-03,
626 + 0.2025E-02, 0.4476E-03, 0.5327E-02, 0.2633E-02,
627 + 0.0000E+00, 0.0000E+00, 0.1038E-02, 0.5258E-03,
628 + 0.2090E-02, 0.5258E-03, 0.5847E-02, 0.2554E-02,
629 + 0.0000E+00, 0.0000E+00, 0.9374E-03, 0.6121E-03,
630 + 0.2162E-02, 0.6121E-03, 0.6418E-02, 0.2500E-02,
631 + 0.0000E+00, 0.0000E+00, 0.8271E-03, 0.7071E-03,
632 + 0.2241E-02, 0.7071E-03, 0.7044E-02, 0.2471E-02,
633 + 0.0000E+00, 0.0000E+00, 0.7066E-03, 0.8115E-03,
634 + 0.2330E-02, 0.8115E-03, 0.7731E-02, 0.2468E-02,
635 + 0.0000E+00, 0.0000E+00, 0.5751E-03, 0.9260E-03,
636 + 0.2427E-02, 0.9260E-03, 0.8484E-02, 0.2494E-02,
637 + 0.0000E+00, 0.0000E+00, 0.4317E-03, 0.1051E-02,
638 + 0.2535E-02, 0.1051E-02, 0.9309E-02, 0.2550E-02,
639 + 0.0000E+00, 0.0000E+00, 0.2757E-03, 0.1189E-02,
640 + 0.2653E-02, 0.1189E-02, 0.1021E-01, 0.2635E-02/
641 DATA ((CALC(I,K,14),I=1,8), K=10,18) /
642 + 0.0000E+00, 0.0000E+00, 0.1062E-03, 0.1340E-02,
643 + 0.2783E-02, 0.1340E-02, 0.1120E-01, 0.2752E-02,
644 + 0.0000E+00, 0.0000E+00, -0.7778E-04, 0.1503E-02,
645 + 0.2927E-02, 0.1503E-02, 0.1228E-01, 0.2900E-02,
646 + 0.0000E+00, 0.0000E+00, -0.2773E-03, 0.1681E-02,
647 + 0.3084E-02, 0.1681E-02, 0.1346E-01, 0.3080E-02,
648 + 0.0000E+00, 0.0000E+00, -0.4936E-03, 0.1875E-02,
649 + 0.3256E-02, 0.1875E-02, 0.1475E-01, 0.3294E-02,
650 + 0.0000E+00, 0.0000E+00, -0.7277E-03, 0.2086E-02,
651 + 0.3444E-02, 0.2086E-02, 0.1615E-01, 0.3541E-02,
652 + 0.0000E+00, 0.0000E+00, -0.9809E-03, 0.2316E-02,
653 + 0.3650E-02, 0.2316E-02, 0.1768E-01, 0.3824E-02,
654 + 0.0000E+00, 0.0000E+00, -0.1255E-02, 0.2565E-02,
655 + 0.3876E-02, 0.2565E-02, 0.1934E-01, 0.4142E-02,
656 + 0.0000E+00, 0.0000E+00, -0.1551E-02, 0.2836E-02,
657 + 0.4122E-02, 0.2836E-02, 0.2115E-01, 0.4497E-02,
658 + 0.0000E+00, 0.0000E+00, -0.1870E-02, 0.3130E-02,
659 + 0.4391E-02, 0.3130E-02, 0.2312E-01, 0.4891E-02/
660 DATA ((CALC(I,K,14),I=1,8), K=19,20) /
661 + 0.0000E+00, 0.0000E+00, -0.2215E-02, 0.3449E-02,
662 + 0.4684E-02, 0.3449E-02, 0.2525E-01, 0.5325E-02,
663 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
664 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
665 DATA ((CALC(I,K,15),I=1,8), K=1,9) /
666 + 0.0000E+00, 0.0000E+00, 0.1209E-02, 0.3877E-03,
667 + 0.1984E-02, 0.3877E-03, 0.4835E-02, 0.2333E-02,
668 + 0.0000E+00, 0.0000E+00, 0.1121E-02, 0.4613E-03,
669 + 0.2044E-02, 0.4613E-03, 0.5323E-02, 0.2216E-02,
670 + 0.0000E+00, 0.0000E+00, 0.1025E-02, 0.5430E-03,
671 + 0.2111E-02, 0.5430E-03, 0.5862E-02, 0.2123E-02,
672 + 0.0000E+00, 0.0000E+00, 0.9181E-03, 0.6336E-03,
673 + 0.2185E-02, 0.6336E-03, 0.6456E-02, 0.2056E-02,
674 + 0.0000E+00, 0.0000E+00, 0.8010E-03, 0.7337E-03,
675 + 0.2268E-02, 0.7337E-03, 0.7112E-02, 0.2015E-02,
676 + 0.0000E+00, 0.0000E+00, 0.6726E-03, 0.8440E-03,
677 + 0.2361E-02, 0.8440E-03, 0.7834E-02, 0.2001E-02,
678 + 0.0000E+00, 0.0000E+00, 0.5322E-03, 0.9654E-03,
679 + 0.2463E-02, 0.9654E-03, 0.8627E-02, 0.2015E-02,
680 + 0.0000E+00, 0.0000E+00, 0.3789E-03, 0.1098E-02,
681 + 0.2576E-02, 0.1098E-02, 0.9499E-02, 0.2057E-02,
682 + 0.0000E+00, 0.0000E+00, 0.2120E-03, 0.1244E-02,
683 + 0.2701E-02, 0.1244E-02, 0.1045E-01, 0.2126E-02/
684 DATA ((CALC(I,K,15),I=1,8), K=10,18) /
685 + 0.0000E+00, 0.0000E+00, 0.3050E-04, 0.1404E-02,
686 + 0.2838E-02, 0.1404E-02, 0.1150E-01, 0.2224E-02,
687 + 0.0000E+00, 0.0000E+00, -0.1666E-03, 0.1578E-02,
688 + 0.2990E-02, 0.1578E-02, 0.1264E-01, 0.2350E-02,
689 + 0.0000E+00, 0.0000E+00, -0.3804E-03, 0.1768E-02,
690 + 0.3156E-02, 0.1768E-02, 0.1389E-01, 0.2505E-02,
691 + 0.0000E+00, 0.0000E+00, -0.6120E-03, 0.1974E-02,
692 + 0.3338E-02, 0.1974E-02, 0.1526E-01, 0.2688E-02,
693 + 0.0000E+00, 0.0000E+00, -0.8626E-03, 0.2199E-02,
694 + 0.3537E-02, 0.2199E-02, 0.1674E-01, 0.2901E-02,
695 + 0.0000E+00, 0.0000E+00, -0.1134E-02, 0.2444E-02,
696 + 0.3756E-02, 0.2444E-02, 0.1836E-01, 0.3143E-02,
697 + 0.0000E+00, 0.0000E+00, -0.1426E-02, 0.2710E-02,
698 + 0.3995E-02, 0.2710E-02, 0.2012E-01, 0.3415E-02,
699 + 0.0000E+00, 0.0000E+00, -0.1743E-02, 0.2999E-02,
700 + 0.4255E-02, 0.2999E-02, 0.2203E-01, 0.3719E-02,
701 + 0.0000E+00, 0.0000E+00, -0.2084E-02, 0.3311E-02,
702 + 0.4540E-02, 0.3311E-02, 0.2411E-01, 0.4054E-02/
703 DATA ((CALC(I,K,15),I=1,8), K=19,20) /
704 + 0.0000E+00, 0.0000E+00, -0.2452E-02, 0.3651E-02,
705 + 0.4851E-02, 0.3651E-02, 0.2637E-01, 0.4422E-02,
706 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
707 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
708 DATA ((CALC(I,K,16),I=1,8), K=1,9) /
709 + 0.0000E+00, 0.0000E+00, 0.1210E-02, 0.3954E-03,
710 + 0.2001E-02, 0.3954E-03, 0.4839E-02, 0.1962E-02,
711 + 0.0000E+00, 0.0000E+00, 0.1118E-02, 0.4727E-03,
712 + 0.2063E-02, 0.4727E-03, 0.5348E-02, 0.1838E-02,
713 + 0.0000E+00, 0.0000E+00, 0.1015E-02, 0.5591E-03,
714 + 0.2133E-02, 0.5591E-03, 0.5913E-02, 0.1740E-02,
715 + 0.0000E+00, 0.0000E+00, 0.9014E-03, 0.6552E-03,
716 + 0.2212E-02, 0.6552E-03, 0.6540E-02, 0.1668E-02,
717 + 0.0000E+00, 0.0000E+00, 0.7761E-03, 0.7618E-03,
718 + 0.2300E-02, 0.7618E-03, 0.7234E-02, 0.1623E-02,
719 + 0.0000E+00, 0.0000E+00, 0.6384E-03, 0.8795E-03,
720 + 0.2397E-02, 0.8795E-03, 0.8001E-02, 0.1604E-02,
721 + 0.0000E+00, 0.0000E+00, 0.4875E-03, 0.1009E-02,
722 + 0.2506E-02, 0.1009E-02, 0.8845E-02, 0.1611E-02,
723 + 0.0000E+00, 0.0000E+00, 0.3226E-03, 0.1152E-02,
724 + 0.2626E-02, 0.1152E-02, 0.9773E-02, 0.1644E-02,
725 + 0.0000E+00, 0.0000E+00, 0.1428E-03, 0.1308E-02,
726 + 0.2759E-02, 0.1308E-02, 0.1079E-01, 0.1703E-02/
727 DATA ((CALC(I,K,16),I=1,8), K=10,18) /
728 + 0.0000E+00, 0.0000E+00, -0.5268E-04, 0.1478E-02,
729 + 0.2905E-02, 0.1478E-02, 0.1191E-01, 0.1788E-02,
730 + 0.0000E+00, 0.0000E+00, -0.2649E-03, 0.1665E-02,
731 + 0.3066E-02, 0.1665E-02, 0.1313E-01, 0.1897E-02,
732 + 0.0000E+00, 0.0000E+00, -0.4950E-03, 0.1868E-02,
733 + 0.3242E-02, 0.1868E-02, 0.1446E-01, 0.2031E-02,
734 + 0.0000E+00, 0.0000E+00, -0.7441E-03, 0.2090E-02,
735 + 0.3436E-02, 0.2090E-02, 0.1591E-01, 0.2190E-02,
736 + 0.0000E+00, 0.0000E+00, -0.1013E-02, 0.2330E-02,
737 + 0.3648E-02, 0.2330E-02, 0.1749E-01, 0.2373E-02,
738 + 0.0000E+00, 0.0000E+00, -0.1304E-02, 0.2592E-02,
739 + 0.3881E-02, 0.2592E-02, 0.1922E-01, 0.2582E-02,
740 + 0.0000E+00, 0.0000E+00, -0.1618E-02, 0.2876E-02,
741 + 0.4134E-02, 0.2876E-02, 0.2109E-01, 0.2815E-02,
742 + 0.0000E+00, 0.0000E+00, -0.1957E-02, 0.3183E-02,
743 + 0.4412E-02, 0.3183E-02, 0.2312E-01, 0.3073E-02,
744 + 0.0000E+00, 0.0000E+00, -0.2322E-02, 0.3517E-02,
745 + 0.4714E-02, 0.3517E-02, 0.2533E-01, 0.3358E-02/
746 DATA ((CALC(I,K,16),I=1,8), K=19,20) /
747 + 0.0000E+00, 0.0000E+00, -0.2715E-02, 0.3878E-02,
748 + 0.5044E-02, 0.3878E-02, 0.2772E-01, 0.3670E-02,
749 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
750 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
751 DATA ((CALC(I,K,17),I=1,8), K=1,9) /
752 + 0.0000E+00, 0.0000E+00, 0.1218E-02, 0.4021E-03,
753 + 0.2022E-02, 0.4021E-03, 0.4871E-02, 0.1625E-02,
754 + 0.0000E+00, 0.0000E+00, 0.1119E-02, 0.4848E-03,
755 + 0.2089E-02, 0.4848E-03, 0.5411E-02, 0.1501E-02,
756 + 0.0000E+00, 0.0000E+00, 0.1009E-02, 0.5774E-03,
757 + 0.2164E-02, 0.5774E-03, 0.6014E-02, 0.1404E-02,
758 + 0.0000E+00, 0.0000E+00, 0.8866E-03, 0.6809E-03,
759 + 0.2248E-02, 0.6809E-03, 0.6686E-02, 0.1333E-02,
760 + 0.0000E+00, 0.0000E+00, 0.7510E-03, 0.7958E-03,
761 + 0.2343E-02, 0.7958E-03, 0.7431E-02, 0.1288E-02,
762 + 0.0000E+00, 0.0000E+00, 0.6017E-03, 0.9229E-03,
763 + 0.2448E-02, 0.9229E-03, 0.8255E-02, 0.1269E-02,
764 + 0.0000E+00, 0.0000E+00, 0.4378E-03, 0.1063E-02,
765 + 0.2564E-02, 0.1063E-02, 0.9164E-02, 0.1274E-02,
766 + 0.0000E+00, 0.0000E+00, 0.2588E-03, 0.1217E-02,
767 + 0.2693E-02, 0.1217E-02, 0.1016E-01, 0.1304E-02,
768 + 0.0000E+00, 0.0000E+00, 0.6363E-04, 0.1385E-02,
769 + 0.2835E-02, 0.1385E-02, 0.1126E-01, 0.1356E-02/
770 DATA ((CALC(I,K,17),I=1,8), K=10,18) /
771 + 0.0000E+00, 0.0000E+00, -0.1485E-03, 0.1569E-02,
772 + 0.2992E-02, 0.1569E-02, 0.1246E-01, 0.1430E-02,
773 + 0.0000E+00, 0.0000E+00, -0.3786E-03, 0.1771E-02,
774 + 0.3164E-02, 0.1771E-02, 0.1377E-01, 0.1526E-02,
775 + 0.0000E+00, 0.0000E+00, -0.6278E-03, 0.1989E-02,
776 + 0.3353E-02, 0.1989E-02, 0.1520E-01, 0.1643E-02,
777 + 0.0000E+00, 0.0000E+00, -0.8971E-03, 0.2227E-02,
778 + 0.3560E-02, 0.2227E-02, 0.1675E-01, 0.1781E-02,
779 + 0.0000E+00, 0.0000E+00, -0.1188E-02, 0.2486E-02,
780 + 0.3787E-02, 0.2486E-02, 0.1845E-01, 0.1939E-02,
781 + 0.0000E+00, 0.0000E+00, -0.1502E-02, 0.2767E-02,
782 + 0.4034E-02, 0.2767E-02, 0.2029E-01, 0.2118E-02,
783 + 0.0000E+00, 0.0000E+00, -0.1839E-02, 0.3071E-02,
784 + 0.4305E-02, 0.3071E-02, 0.2229E-01, 0.2317E-02,
785 + 0.0000E+00, 0.0000E+00, -0.2203E-02, 0.3400E-02,
786 + 0.4600E-02, 0.3400E-02, 0.2445E-01, 0.2536E-02,
787 + 0.0000E+00, 0.0000E+00, -0.2595E-02, 0.3757E-02,
788 + 0.4922E-02, 0.3757E-02, 0.2680E-01, 0.2777E-02/
789 DATA ((CALC(I,K,17),I=1,8), K=19,20) /
790 + 0.0000E+00, 0.0000E+00, -0.3016E-02, 0.4143E-02,
791 + 0.5273E-02, 0.4143E-02, 0.2935E-01, 0.3040E-02,
792 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
793 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
794 DATA ((CALC(I,K,18),I=1,8), K=1,9) /
795 + 0.0000E+00, 0.0000E+00, 0.1235E-02, 0.4100E-03,
796 + 0.2055E-02, 0.4100E-03, 0.4941E-02, 0.1323E-02,
797 + 0.0000E+00, 0.0000E+00, 0.1128E-02, 0.5000E-03,
798 + 0.2128E-02, 0.5000E-03, 0.5526E-02, 0.1205E-02,
799 + 0.0000E+00, 0.0000E+00, 0.1008E-02, 0.6013E-03,
800 + 0.2210E-02, 0.6013E-03, 0.6182E-02, 0.1114E-02,
801 + 0.0000E+00, 0.0000E+00, 0.8735E-03, 0.7144E-03,
802 + 0.2302E-02, 0.7144E-03, 0.6913E-02, 0.1048E-02,
803 + 0.0000E+00, 0.0000E+00, 0.7246E-03, 0.8401E-03,
804 + 0.2405E-02, 0.8401E-03, 0.7724E-02, 0.1008E-02,
805 + 0.0000E+00, 0.0000E+00, 0.5606E-03, 0.9790E-03,
806 + 0.2519E-02, 0.9790E-03, 0.8621E-02, 0.9914E-03,
807 + 0.0000E+00, 0.0000E+00, 0.3806E-03, 0.1132E-02,
808 + 0.2645E-02, 0.1132E-02, 0.9609E-02, 0.9973E-03,
809 + 0.0000E+00, 0.0000E+00, 0.1842E-03, 0.1300E-02,
810 + 0.2784E-02, 0.1300E-02, 0.1069E-01, 0.1024E-02,
811 + 0.0000E+00, 0.0000E+00, -0.2965E-04, 0.1483E-02,
812 + 0.2937E-02, 0.1483E-02, 0.1188E-01, 0.1072E-02/
813 DATA ((CALC(I,K,18),I=1,8), K=10,18) /
814 + 0.0000E+00, 0.0000E+00, -0.2617E-03, 0.1683E-02,
815 + 0.3106E-02, 0.1683E-02, 0.1318E-01, 0.1138E-02,
816 + 0.0000E+00, 0.0000E+00, -0.5129E-03, 0.1901E-02,
817 + 0.3291E-02, 0.1901E-02, 0.1460E-01, 0.1223E-02,
818 + 0.0000E+00, 0.0000E+00, -0.7843E-03, 0.2138E-02,
819 + 0.3494E-02, 0.2138E-02, 0.1614E-01, 0.1326E-02,
820 + 0.0000E+00, 0.0000E+00, -0.1077E-02, 0.2396E-02,
821 + 0.3716E-02, 0.2396E-02, 0.1782E-01, 0.1445E-02,
822 + 0.0000E+00, 0.0000E+00, -0.1392E-02, 0.2674E-02,
823 + 0.3959E-02, 0.2674E-02, 0.1964E-01, 0.1581E-02,
824 + 0.0000E+00, 0.0000E+00, -0.1732E-02, 0.2977E-02,
825 + 0.4224E-02, 0.2977E-02, 0.2161E-01, 0.1734E-02,
826 + 0.0000E+00, 0.0000E+00, -0.2096E-02, 0.3303E-02,
827 + 0.4513E-02, 0.3303E-02, 0.2375E-01, 0.1902E-02,
828 + 0.0000E+00, 0.0000E+00, -0.2488E-02, 0.3657E-02,
829 + 0.4828E-02, 0.3657E-02, 0.2607E-01, 0.2087E-02,
830 + 0.0000E+00, 0.0000E+00, -0.2909E-02, 0.4039E-02,
831 + 0.5171E-02, 0.4039E-02, 0.2858E-01, 0.2290E-02/
832 DATA ((CALC(I,K,18),I=1,8), K=19,20) /
833 + 0.0000E+00, 0.0000E+00, -0.3361E-02, 0.4451E-02,
834 + 0.5545E-02, 0.4451E-02, 0.3129E-01, 0.2510E-02,
835 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
836 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
837 DATA ((CALC(I,K,19),I=1,8), K=1,9) /
838 + 0.0000E+00, 0.0000E+00, 0.1265E-02, 0.4210E-03,
839 + 0.2108E-02, 0.4210E-03, 0.5061E-02, 0.1058E-02,
840 + 0.0000E+00, 0.0000E+00, 0.1146E-02, 0.5211E-03,
841 + 0.2189E-02, 0.5211E-03, 0.5708E-02, 0.9493E-03,
842 + 0.0000E+00, 0.0000E+00, 0.1012E-02, 0.6336E-03,
843 + 0.2279E-02, 0.6336E-03, 0.6433E-02, 0.8673E-03,
844 + 0.0000E+00, 0.0000E+00, 0.8619E-03, 0.7591E-03,
845 + 0.2380E-02, 0.7591E-03, 0.7241E-02, 0.8103E-03,
846 + 0.0000E+00, 0.0000E+00, 0.6959E-03, 0.8982E-03,
847 + 0.2492E-02, 0.8982E-03, 0.8136E-02, 0.7765E-03,
848 + 0.0000E+00, 0.0000E+00, 0.5132E-03, 0.1052E-02,
849 + 0.2617E-02, 0.1052E-02, 0.9123E-02, 0.7644E-03,
850 + 0.0000E+00, 0.0000E+00, 0.3133E-03, 0.1220E-02,
851 + 0.2754E-02, 0.1220E-02, 0.1021E-01, 0.7723E-03,
852 + 0.0000E+00, 0.0000E+00, 0.9558E-04, 0.1404E-02,
853 + 0.2905E-02, 0.1404E-02, 0.1139E-01, 0.7988E-03,
854 + 0.0000E+00, 0.0000E+00, -0.1408E-03, 0.1605E-02,
855 + 0.3071E-02, 0.1605E-02, 0.1269E-01, 0.8425E-03/
856 DATA ((CALC(I,K,19),I=1,8), K=10,18) /
857 + 0.0000E+00, 0.0000E+00, -0.3964E-03, 0.1823E-02,
858 + 0.3253E-02, 0.1823E-02, 0.1411E-01, 0.9023E-03,
859 + 0.0000E+00, 0.0000E+00, -0.6723E-03, 0.2062E-02,
860 + 0.3452E-02, 0.2062E-02, 0.1564E-01, 0.9773E-03,
861 + 0.0000E+00, 0.0000E+00, -0.9694E-03, 0.2319E-02,
862 + 0.3670E-02, 0.2319E-02, 0.1731E-01, 0.1067E-02,
863 + 0.0000E+00, 0.0000E+00, -0.1289E-02, 0.2597E-02,
864 + 0.3908E-02, 0.2597E-02, 0.1912E-01, 0.1169E-02,
865 + 0.0000E+00, 0.0000E+00, -0.1632E-02, 0.2898E-02,
866 + 0.4168E-02, 0.2898E-02, 0.2108E-01, 0.1286E-02,
867 + 0.0000E+00, 0.0000E+00, -0.2000E-02, 0.3224E-02,
868 + 0.4452E-02, 0.3224E-02, 0.2321E-01, 0.1415E-02,
869 + 0.0000E+00, 0.0000E+00, -0.2394E-02, 0.3576E-02,
870 + 0.4761E-02, 0.3576E-02, 0.2550E-01, 0.1556E-02,
871 + 0.0000E+00, 0.0000E+00, -0.2817E-02, 0.3955E-02,
872 + 0.5097E-02, 0.3955E-02, 0.2798E-01, 0.1712E-02,
873 + 0.0000E+00, 0.0000E+00, -0.3269E-02, 0.4364E-02,
874 + 0.5464E-02, 0.4364E-02, 0.3066E-01, 0.1879E-02/
875 DATA ((CALC(I,K,19),I=1,8), K=19,20) /
876 + 0.0000E+00, 0.0000E+00, -0.3753E-02, 0.4805E-02,
877 + 0.5862E-02, 0.4805E-02, 0.3356E-01, 0.2062E-02,
878 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
879 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
880 DATA ((CALC(I,K,20),I=1,8), K=1,9) /
881 + 0.0000E+00, 0.0000E+00, 0.1311E-02, 0.4367E-03,
882 + 0.2185E-02, 0.4367E-03, 0.5245E-02, 0.8288E-03,
883 + 0.0000E+00, 0.0000E+00, 0.1175E-02, 0.5500E-03,
884 + 0.2275E-02, 0.5500E-03, 0.5974E-02, 0.7327E-03,
885 + 0.0000E+00, 0.0000E+00, 0.1022E-02, 0.6769E-03,
886 + 0.2376E-02, 0.6769E-03, 0.6787E-02, 0.6623E-03,
887 + 0.0000E+00, 0.0000E+00, 0.8518E-03, 0.8178E-03,
888 + 0.2487E-02, 0.8178E-03, 0.7689E-02, 0.6149E-03,
889 + 0.0000E+00, 0.0000E+00, 0.6638E-03, 0.9733E-03,
890 + 0.2610E-02, 0.9733E-03, 0.8685E-02, 0.5888E-03,
891 + 0.0000E+00, 0.0000E+00, 0.4578E-03, 0.1144E-02,
892 + 0.2746E-02, 0.1144E-02, 0.9779E-02, 0.5818E-03,
893 + 0.0000E+00, 0.0000E+00, 0.2333E-03, 0.1331E-02,
894 + 0.2895E-02, 0.1331E-02, 0.1098E-01, 0.5922E-03,
895 + 0.0000E+00, 0.0000E+00, -0.1005E-04, 0.1534E-02,
896 + 0.3059E-02, 0.1534E-02, 0.1228E-01, 0.6184E-03,
897 + 0.0000E+00, 0.0000E+00, -0.2729E-03, 0.1755E-02,
898 + 0.3238E-02, 0.1755E-02, 0.1370E-01, 0.6589E-03/
899 DATA ((CALC(I,K,20),I=1,8), K=10,18) /
900 + 0.0000E+00, 0.0000E+00, -0.5560E-03, 0.1989E-02,
901 + 0.3434E-02, 0.1989E-02, 0.1524E-01, 0.7126E-03,
902 + 0.0000E+00, 0.0000E+00, -0.8602E-03, 0.2249E-02,
903 + 0.3648E-02, 0.2249E-02, 0.1691E-01, 0.7783E-03,
904 + 0.0000E+00, 0.0000E+00, -0.1186E-02, 0.2529E-02,
905 + 0.3882E-02, 0.2529E-02, 0.1872E-01, 0.8554E-03,
906 + 0.0000E+00, 0.0000E+00, -0.1535E-02, 0.2832E-02,
907 + 0.4138E-02, 0.2832E-02, 0.2068E-01, 0.9431E-03,
908 + 0.0000E+00, 0.0000E+00, -0.1909E-02, 0.3158E-02,
909 + 0.4416E-02, 0.3158E-02, 0.2279E-01, 0.1041E-02,
910 + 0.0000E+00, 0.0000E+00, -0.2308E-02, 0.3509E-02,
911 + 0.4718E-02, 0.3509E-02, 0.2507E-01, 0.1149E-02,
912 + 0.0000E+00, 0.0000E+00, -0.2734E-02, 0.3887E-02,
913 + 0.5048E-02, 0.3887E-02, 0.2753E-01, 0.1267E-02,
914 + 0.0000E+00, 0.0000E+00, -0.3189E-02, 0.4294E-02,
915 + 0.5406E-02, 0.4294E-02, 0.3018E-01, 0.1394E-02,
916 + 0.0000E+00, 0.0000E+00, -0.3675E-02, 0.4732E-02,
917 + 0.5796E-02, 0.4732E-02, 0.3305E-01, 0.1533E-02/
918 DATA ((CALC(I,K,20),I=1,8), K=19,20) /
919 + 0.0000E+00, 0.0000E+00, -0.4194E-02, 0.5202E-02,
920 + 0.6219E-02, 0.5202E-02, 0.3613E-01, 0.1682E-02,
921 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
922 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
923 DATA ((CALC(I,K,21),I=1,8), K=1,9) /
924 + 0.0000E+00, 0.0000E+00, 0.1377E-02, 0.4583E-03,
925 + 0.2293E-02, 0.4583E-03, 0.5506E-02, 0.6351E-03,
926 + 0.0000E+00, 0.0000E+00, 0.1217E-02, 0.5884E-03,
927 + 0.2394E-02, 0.5884E-03, 0.6336E-02, 0.5532E-03,
928 + 0.0000E+00, 0.0000E+00, 0.1039E-02, 0.7330E-03,
929 + 0.2505E-02, 0.7330E-03, 0.7257E-02, 0.4951E-03,
930 + 0.0000E+00, 0.0000E+00, 0.8426E-03, 0.8925E-03,
931 + 0.2628E-02, 0.8925E-03, 0.8272E-02, 0.4579E-03,
932 + 0.0000E+00, 0.0000E+00, 0.6271E-03, 0.1068E-02,
933 + 0.2762E-02, 0.1068E-02, 0.9385E-02, 0.4395E-03,
934 + 0.0000E+00, 0.0000E+00, 0.3917E-03, 0.1258E-02,
935 + 0.2909E-02, 0.1258E-02, 0.1060E-01, 0.4375E-03,
936 + 0.0000E+00, 0.0000E+00, 0.1382E-03, 0.1466E-02,
937 + 0.3071E-02, 0.1466E-02, 0.1192E-01, 0.4502E-03,
938 + 0.0000E+00, 0.0000E+00, -0.1353E-03, 0.1690E-02,
939 + 0.3247E-02, 0.1690E-02, 0.1336E-01, 0.4759E-03,
940 + 0.0000E+00, 0.0000E+00, -0.4289E-03, 0.1933E-02,
941 + 0.3439E-02, 0.1933E-02, 0.1492E-01, 0.5131E-03/
942 DATA ((CALC(I,K,21),I=1,8), K=10,18) /
943 + 0.0000E+00, 0.0000E+00, -0.7432E-03, 0.2195E-02,
944 + 0.3650E-02, 0.2195E-02, 0.1660E-01, 0.5606E-03,
945 + 0.0000E+00, 0.0000E+00, -0.1079E-02, 0.2479E-02,
946 + 0.3879E-02, 0.2479E-02, 0.1841E-01, 0.6176E-03,
947 + 0.0000E+00, 0.0000E+00, -0.1437E-02, 0.2782E-02,
948 + 0.4129E-02, 0.2782E-02, 0.2036E-01, 0.6831E-03,
949 + 0.0000E+00, 0.0000E+00, -0.1818E-02, 0.3108E-02,
950 + 0.4400E-02, 0.3108E-02, 0.2247E-01, 0.7567E-03,
951 + 0.0000E+00, 0.0000E+00, -0.2224E-02, 0.3459E-02,
952 + 0.4696E-02, 0.3459E-02, 0.2474E-01, 0.8379E-03,
953 + 0.0000E+00, 0.0000E+00, -0.2655E-02, 0.3835E-02,
954 + 0.5018E-02, 0.3835E-02, 0.2718E-01, 0.9261E-03,
955 + 0.0000E+00, 0.0000E+00, -0.3114E-02, 0.4239E-02,
956 + 0.5367E-02, 0.4239E-02, 0.2980E-01, 0.1021E-02,
957 + 0.0000E+00, 0.0000E+00, -0.3603E-02, 0.4673E-02,
958 + 0.5747E-02, 0.4673E-02, 0.3265E-01, 0.1132E-02,
959 + 0.0000E+00, 0.0000E+00, -0.4122E-02, 0.5139E-02,
960 + 0.6160E-02, 0.5139E-02, 0.3569E-01, 0.1240E-02/
961 DATA ((CALC(I,K,21),I=1,8), K=19,20) /
962 + 0.0000E+00, 0.0000E+00, -0.4675E-02, 0.5640E-02,
963 + 0.6608E-02, 0.5640E-02, 0.3896E-01, 0.1357E-02,
964 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
965 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
966 DATA ((CALC(I,K,22),I=1,8), K=1,9) /
967 + 0.0000E+00, 0.0000E+00, 0.1462E-02, 0.4873E-03,
968 + 0.2437E-02, 0.4873E-03, 0.5848E-02, 0.4746E-03,
969 + 0.0000E+00, 0.0000E+00, 0.1272E-02, 0.6379E-03,
970 + 0.2548E-02, 0.6379E-03, 0.6801E-02, 0.4074E-03,
971 + 0.0000E+00, 0.0000E+00, 0.1062E-02, 0.8038E-03,
972 + 0.2670E-02, 0.8038E-03, 0.7849E-02, 0.3617E-03,
973 + 0.0000E+00, 0.0000E+00, 0.8324E-03, 0.9852E-03,
974 + 0.2803E-02, 0.9852E-03, 0.8994E-02, 0.3345E-03,
975 + 0.0000E+00, 0.0000E+00, 0.5830E-03, 0.1182E-02,
976 + 0.2948E-02, 0.1182E-02, 0.1024E-01, 0.3231E-03,
977 + 0.0000E+00, 0.0000E+00, 0.3145E-03, 0.1396E-02,
978 + 0.3106E-02, 0.1396E-02, 0.1159E-01, 0.3256E-03,
979 + 0.0000E+00, 0.0000E+00, 0.2579E-04, 0.1627E-02,
980 + 0.3279E-02, 0.1627E-02, 0.1305E-01, 0.3399E-03,
981 + 0.0000E+00, 0.0000E+00, -0.2823E-03, 0.1870E-02,
982 + 0.3467E-02, 0.1870E-02, 0.1463E-01, 0.3645E-03,
983 + 0.0000E+00, 0.0000E+00, -0.6105E-03, 0.2145E-02,
984 + 0.3672E-02, 0.2145E-02, 0.1632E-01, 0.3980E-03/
985 DATA ((CALC(I,K,22),I=1,8), K=10,18) /
986 + 0.0000E+00, 0.0000E+00, -0.9591E-03, 0.2431E-02,
987 + 0.3895E-02, 0.2431E-02, 0.1815E-01, 0.4394E-03,
988 + 0.0000E+00, 0.0000E+00, -0.1329E-02, 0.2736E-02,
989 + 0.4138E-02, 0.2736E-02, 0.2011E-01, 0.4877E-03,
990 + 0.0000E+00, 0.0000E+00, -0.1721E-02, 0.3061E-02,
991 + 0.4403E-02, 0.3061E-02, 0.2221E-01, 0.5423E-03,
992 + 0.0000E+00, 0.0000E+00, -0.2135E-02, 0.3411E-02,
993 + 0.4690E-02, 0.3411E-02, 0.2447E-01, 0.6027E-03,
994 + 0.0000E+00, 0.0000E+00, -0.2574E-02, 0.3786E-02,
995 + 0.5002E-02, 0.3786E-02, 0.2689E-01, 0.6685E-03,
996 + 0.0000E+00, 0.0000E+00, -0.3038E-02, 0.4187E-02,
997 + 0.5341E-02, 0.4187E-02, 0.2949E-01, 0.7397E-03,
998 + 0.0000E+00, 0.0000E+00, -0.3530E-02, 0.4617E-02,
999 + 0.5710E-02, 0.4617E-02, 0.3227E-01, 0.8160E-03,
1000 + 0.0000E+00, 0.0000E+00, -0.4051E-02, 0.5077E-02,
1001 + 0.6110E-02, 0.5077E-02, 0.3526E-01, 0.8979E-03,
1002 + 0.0000E+00, 0.0000E+00, -0.4603E-02, 0.5570E-02,
1003 + 0.6544E-02, 0.5570E-02, 0.3848E-01, 0.9852E-03/
1004 DATA ((CALC(I,K,22),I=1,8), K=19,20) /
1005 + 0.0000E+00, 0.0000E+00, -0.5189E-02, 0.6098E-02,
1006 + 0.7015E-02, 0.6098E-02, 0.4193E-01, 0.1079E-02,
1007 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
1008 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
1009 DATA ((CALC(I,K,23),I=1,8), K=1,9) /
1010 + 0.0000E+00, 0.0000E+00, 0.1571E-02, 0.5237E-03,
1011 + 0.2618E-02, 0.5237E-03, 0.6284E-02, 0.3454E-03,
1012 + 0.0000E+00, 0.0000E+00, 0.1341E-02, 0.6990E-03,
1013 + 0.2739E-02, 0.6990E-03, 0.7380E-02, 0.2925E-03,
1014 + 0.0000E+00, 0.0000E+00, 0.1091E-02, 0.8897E-03,
1015 + 0.2870E-02, 0.8897E-03, 0.8573E-02, 0.2584E-03,
1016 + 0.0000E+00, 0.0000E+00, 0.8202E-03, 0.1096E-02,
1017 + 0.3012E-02, 0.1096E-02, 0.9863E-02, 0.2400E-03,
1018 + 0.0000E+00, 0.0000E+00, 0.5302E-03, 0.1318E-02,
1019 + 0.3166E-02, 0.1318E-02, 0.1125E-01, 0.2347E-03,
1020 + 0.0000E+00, 0.0000E+00, 0.2214E-03, 0.1557E-02,
1021 + 0.3334E-02, 0.1557E-02, 0.1275E-01, 0.2404E-03,
1022 + 0.0000E+00, 0.0000E+00, -0.1069E-03, 0.1728E-02,
1023 + 0.3516E-02, 0.1728E-02, 0.1435E-01, 0.2554E-03,
1024 + 0.0000E+00, 0.0000E+00, -0.4535E-03, 0.2033E-02,
1025 + 0.3715E-02, 0.2033E-02, 0.1607E-01, 0.2782E-03,
1026 + 0.0000E+00, 0.0000E+00, -0.8196E-03, 0.2349E-02,
1027 + 0.3930E-02, 0.2349E-02, 0.1791E-01, 0.3075E-03/
1028 DATA ((CALC(I,K,23),I=1,8), K=10,18) /
1029 + 0.0000E+00, 0.0000E+00, -0.1205E-02, 0.2657E-02,
1030 + 0.4164E-02, 0.2657E-02, 0.1987E-01, 0.3425E-03,
1031 + 0.0000E+00, 0.0000E+00, -0.1610E-02, 0.2992E-02,
1032 + 0.4419E-02, 0.2992E-02, 0.2197E-01, 0.3825E-03,
1033 + 0.0000E+00, 0.0000E+00, -0.2036E-02, 0.3345E-02,
1034 + 0.4695E-02, 0.3345E-02, 0.2422E-01, 0.4267E-03,
1035 + 0.0000E+00, 0.0000E+00, -0.2484E-02, 0.3721E-02,
1036 + 0.4996E-02, 0.3721E-02, 0.2661E-01, 0.4750E-03,
1037 + 0.0000E+00, 0.0000E+00, -0.2955E-02, 0.4122E-02,
1038 + 0.5323E-02, 0.4122E-02, 0.2918E-01, 0.5270E-03,
1039 + 0.0000E+00, 0.0000E+00, -0.3451E-02, 0.4549E-02,
1040 + 0.5677E-02, 0.4549E-02, 0.3192E-01, 0.5827E-03,
1041 + 0.0000E+00, 0.0000E+00, -0.3973E-02, 0.5005E-02,
1042 + 0.6061E-02, 0.5005E-02, 0.3485E-01, 0.6419E-03,
1043 + 0.0000E+00, 0.0000E+00, -0.4524E-02, 0.5490E-02,
1044 + 0.6479E-02, 0.5490E-02, 0.3799E-01, 0.7056E-03,
1045 + 0.0000E+00, 0.0000E+00, -0.5105E-02, 0.6009E-02,
1046 + 0.6932E-02, 0.6009E-02, 0.4136E-01, 0.7729E-03/
1047 DATA ((CALC(I,K,23),I=1,8), K=19,20) /
1048 + 0.0000E+00, 0.0000E+00, -0.5720E-02, 0.6563E-02,
1049 + 0.7423E-02, 0.6563E-02, 0.4496E-01, 0.8437E-03,
1050 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
1051 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
1052 DATA ((CALC(I,K,24),I=1,8), K=1,9) /
1053 + 0.0000E+00, 0.0000E+00, 0.1704E-02, 0.5681E-03,
1054 + 0.2840E-02, 0.5681E-03, 0.6817E-02, 0.2441E-03,
1055 + 0.0000E+00, 0.0000E+00, 0.1424E-02, 0.7722E-03,
1056 + 0.2968E-02, 0.7722E-03, 0.8076E-02, 0.2044E-03,
1057 + 0.0000E+00, 0.0000E+00, 0.1123E-02, 0.9913E-03,
1058 + 0.3105E-02, 0.9913E-03, 0.9428E-02, 0.1806E-03,
1059 + 0.0000E+00, 0.0000E+00, 0.8029E-03, 0.1225E-02,
1060 + 0.3254E-02, 0.1225E-02, 0.1087E-01, 0.1695E-03,
1061 + 0.0000E+00, 0.0000E+00, 0.4652E-03, 0.1475E-02,
1062 + 0.3414E-02, 0.1475E-02, 0.1242E-01, 0.1688E-03,
1063 + 0.0000E+00, 0.0000E+00, 0.1104E-03, 0.1739E-02,
1064 + 0.3588E-02, 0.1739E-02, 0.1406E-01, 0.1766E-03,
1065 + 0.0000E+00, 0.0000E+00, -0.2620E-03, 0.2020E-02,
1066 + 0.3777E-02, 0.2020E-02, 0.1581E-01, 0.1913E-03,
1067 + 0.0000E+00, 0.0000E+00, -0.6507E-03, 0.2317E-02,
1068 + 0.3982E-02, 0.2317E-02, 0.1767E-01, 0.2115E-03,
1069 + 0.0000E+00, 0.0000E+00, -0.1057E-02, 0.2631E-02,
1070 + 0.4205E-02, 0.2631E-02, 0.1964E-01, 0.2364E-03/
1071 DATA ((CALC(I,K,24),I=1,8), K=10,18) /
1072 + 0.0000E+00, 0.0000E+00, -0.1480E-02, 0.2963E-02,
1073 + 0.4447E-02, 0.2963E-02, 0.2173E-01, 0.2650E-03,
1074 + 0.0000E+00, 0.0000E+00, -0.1920E-02, 0.3315E-02,
1075 + 0.4710E-02, 0.3315E-02, 0.2396E-01, 0.2970E-03,
1076 + 0.0000E+00, 0.0000E+00, -0.2380E-02, 0.3687E-02,
1077 + 0.4995E-02, 0.3687E-02, 0.2633E-01, 0.3317E-03,
1078 + 0.0000E+00, 0.0000E+00, -0.2859E-02, 0.4082E-02,
1079 + 0.5305E-02, 0.4082E-02, 0.2885E-01, 0.3689E-03,
1080 + 0.0000E+00, 0.0000E+00, -0.3359E-02, 0.4501E-02,
1081 + 0.5642E-02, 0.4501E-02, 0.3153E-01, 0.4087E-03,
1082 + 0.0000E+00, 0.0000E+00, -0.3883E-02, 0.4945E-02,
1083 + 0.6008E-02, 0.4945E-02, 0.3439E-01, 0.4508E-03,
1084 + 0.0000E+00, 0.0000E+00, -0.4431E-02, 0.5418E-02,
1085 + 0.6405E-02, 0.5418E-02, 0.3744E-01, 0.4954E-03,
1086 + 0.0000E+00, 0.0000E+00, -0.5007E-02, 0.5921E-02,
1087 + 0.6836E-02, 0.5921E-02, 0.4070E-01, 0.5424E-03,
1088 + 0.0000E+00, 0.0000E+00, -0.5612E-02, 0.6458E-02,
1089 + 0.7304E-02, 0.6458E-02, 0.4419E-01, 0.5924E-03/
1090 DATA ((CALC(I,K,24),I=1,8), K=19,20) /
1091 + 0.0000E+00, 0.0000E+00, -0.6249E-02, 0.7031E-02,
1092 + 0.7811E-02, 0.7031E-02, 0.4791E-01, 0.6447E-03,
1093 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
1094 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
1095 DATA ((CALC(I,K,25),I=1,8), K=1,9) /
1096 + 0.0000E+00, 0.0000E+00, 0.1862E-02, 0.6211E-03,
1097 + 0.3104E-02, 0.6211E-03, 0.7449E-02, 0.1673E-03,
1098 + 0.0000E+00, 0.0000E+00, 0.1518E-02, 0.8583E-03,
1099 + 0.3234E-02, 0.8583E-03, 0.8888E-02, 0.1390E-03,
1100 + 0.0000E+00, 0.0000E+00, 0.1156E-02, 0.1109E-02,
1101 + 0.3373E-02, 0.1109E-02, 0.1041E-01, 0.1237E-03,
1102 + 0.0000E+00, 0.0000E+00, 0.7779E-03, 0.1374E-02,
1103 + 0.3523E-02, 0.1374E-02, 0.1202E-01, 0.1183E-03,
1104 + 0.0000E+00, 0.0000E+00, 0.3848E-03, 0.1651E-02,
1105 + 0.3686E-02, 0.1651E-02, 0.1372E-01, 0.1208E-03,
1106 + 0.0000E+00, 0.0000E+00, -0.2175E-04, 0.1942E-02,
1107 + 0.3862E-02, 0.1942E-02, 0.1551E-01, 0.1295E-03,
1108 + 0.0000E+00, 0.0000E+00, -0.4424E-03, 0.2248E-02,
1109 + 0.4053E-02, 0.2248E-02, 0.1739E-01, 0.1429E-03,
1110 + 0.0000E+00, 0.0000E+00, -0.8755E-03, 0.2568E-02,
1111 + 0.4259E-02, 0.2568E-02, 0.1937E-01, 0.1600E-03,
1112 + 0.0000E+00, 0.0000E+00, -0.1322E-02, 0.2904E-02,
1113 + 0.4485E-02, 0.2904E-02, 0.2147E-01, 0.1801E-03/
1114 DATA ((CALC(I,K,25),I=1,8), K=10,18) /
1115 + 0.0000E+00, 0.0000E+00, -0.1782E-02, 0.3256E-02,
1116 + 0.4730E-02, 0.3256E-02, 0.2367E-01, 0.2026E-03,
1117 + 0.0000E+00, 0.0000E+00, -0.2256E-02, 0.3627E-02,
1118 + 0.4997E-02, 0.3627E-02, 0.2601E-01, 0.2271E-03,
1119 + 0.0000E+00, 0.0000E+00, -0.2746E-02, 0.4017E-02,
1120 + 0.5288E-02, 0.4017E-02, 0.2847E-01, 0.2532E-03,
1121 + 0.0000E+00, 0.0000E+00, -0.3252E-02, 0.4427E-02,
1122 + 0.5603E-02, 0.4427E-02, 0.3109E-01, 0.2808E-03,
1123 + 0.0000E+00, 0.0000E+00, -0.3777E-02, 0.4861E-02,
1124 + 0.5946E-02, 0.4861E-02, 0.3386E-01, 0.3099E-03,
1125 + 0.0000E+00, 0.0000E+00, -0.4323E-02, 0.5320E-02,
1126 + 0.6318E-02, 0.5320E-02, 0.3680E-01, 0.3404E-03,
1127 + 0.0000E+00, 0.0000E+00, -0.4891E-02, 0.5807E-02,
1128 + 0.6723E-02, 0.5807E-02, 0.3994E-01, 0.3725E-03,
1129 + 0.0000E+00, 0.0000E+00, -0.5484E-02, 0.6323E-02,
1130 + 0.7162E-02, 0.6323E-02, 0.4327E-01, 0.4061E-03,
1131 + 0.0000E+00, 0.0000E+00, -0.6105E-02, 0.6873E-02,
1132 + 0.7639E-02, 0.6873E-02, 0.4684E-01, 0.4420E-03/
1133 DATA ((CALC(I,K,25),I=1,8), K=19,20) /
1134 + 0.0000E+00, 0.0000E+00, -0.6757E-02, 0.7458E-02,
1135 + 0.8157E-02, 0.7458E-02, 0.5065E-01, 0.4794E-03,
1136 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
1137 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
1138 DATA ((CALC(I,K,26),I=1,8), K=1,9) /
1139 + 0.0000E+00, 0.0000E+00, 0.2044E-02, 0.6816E-03,
1140 + 0.3407E-02, 0.6816E-03, 0.8178E-02, 0.1109E-03,
1141 + 0.0000E+00, 0.0000E+00, 0.1620E-02, 0.9567E-03,
1142 + 0.3533E-02, 0.9567E-03, 0.9813E-02, 0.9210E-04,
1143 + 0.0000E+00, 0.0000E+00, 0.1184E-02, 0.1243E-02,
1144 + 0.3668E-02, 0.1243E-02, 0.1152E-01, 0.8342E-04,
1145 + 0.0000E+00, 0.0000E+00, 0.7369E-03, 0.1539E-02,
1146 + 0.3814E-02, 0.1539E-02, 0.1329E-01, 0.8212E-04,
1147 + 0.0000E+00, 0.0000E+00, 0.2815E-03, 0.1846E-02,
1148 + 0.3973E-02, 0.1846E-02, 0.1514E-01, 0.8635E-04,
1149 + 0.0000E+00, 0.0000E+00, -0.1817E-03, 0.2164E-02,
1150 + 0.4145E-02, 0.2164E-02, 0.1707E-01, 0.9473E-04,
1151 + 0.0000E+00, 0.0000E+00, -0.6536E-03, 0.2493E-02,
1152 + 0.4333E-02, 0.2493E-02, 0.1907E-01, 0.1062E-03,
1153 + 0.0000E+00, 0.0000E+00, -0.1132E-02, 0.2834E-02,
1154 + 0.4537E-02, 0.2834E-02, 0.2117E-01, 0.1199E-03,
1155 + 0.0000E+00, 0.0000E+00, -0.1618E-02, 0.3188E-02,
1156 + 0.4760E-02, 0.3188E-02, 0.2335E-01, 0.1354E-03/
1157 DATA ((CALC(I,K,26),I=1,8), K=10,18) /
1158 + 0.0000E+00, 0.0000E+00, -0.2112E-02, 0.3557E-02,
1159 + 0.5004E-02, 0.3557E-02, 0.2564E-01, 0.1521E-03,
1160 + 0.0000E+00, 0.0000E+00, -0.2615E-02, 0.3941E-02,
1161 + 0.5270E-02, 0.3941E-02, 0.2805E-01, 0.1698E-03,
1162 + 0.0000E+00, 0.0000E+00, -0.3130E-02, 0.4343E-02,
1163 + 0.5561E-02, 0.4343E-02, 0.3058E-01, 0.1885E-03,
1164 + 0.0000E+00, 0.0000E+00, -0.3657E-02, 0.4765E-02,
1165 + 0.5876E-02, 0.4765E-02, 0.3325E-01, 0.2079E-03,
1166 + 0.0000E+00, 0.0000E+00, -0.4199E-02, 0.5208E-02,
1167 + 0.6219E-02, 0.5208E-02, 0.3607E-01, 0.2281E-03,
1168 + 0.0000E+00, 0.0000E+00, -0.4758E-02, 0.5674E-02,
1169 + 0.6592E-02, 0.5674E-02, 0.3905E-01, 0.2492E-03,
1170 + 0.0000E+00, 0.0000E+00, -0.5337E-02, 0.6167E-02,
1171 + 0.6998E-02, 0.6167E-02, 0.4222E-01, 0.2712E-03,
1172 + 0.0000E+00, 0.0000E+00, -0.5938E-02, 0.6688E-02,
1173 + 0.7440E-02, 0.6688E-02, 0.4559E-01, 0.2939E-03,
1174 + 0.0000E+00, 0.0000E+00, -0.6565E-02, 0.7241E-02,
1175 + 0.7919E-02, 0.7241E-02, 0.4919E-01, 0.3184E-03/
1176 DATA ((CALC(I,K,26),I=1,8), K=19,20) /
1177 + 0.0000E+00, 0.0000E+00, -0.7220E-02, 0.7829E-02,
1178 + 0.8440E-02, 0.7829E-02, 0.5300E-01, 0.3430E-03,
1179 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
1180 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
1181 DATA ((CALC(I,K,27),I=1,8), K=1,9) /
1182 + 0.0000E+00, 0.0000E+00, 0.2253E-02, 0.7503E-03,
1183 + 0.3754E-02, 0.7503E-03, 0.9011E-02, 0.7114E-04,
1184 + 0.0000E+00, 0.0000E+00, 0.1729E-02, 0.1069E-02,
1185 + 0.3867E-02, 0.1069E-02, 0.1086E-01, 0.5971E-04,
1186 + 0.0000E+00, 0.0000E+00, 0.1203E-02, 0.1393E-02,
1187 + 0.3989E-02, 0.1393E-02, 0.1275E-01, 0.5573E-04,
1188 + 0.0000E+00, 0.0000E+00, 0.6762E-03, 0.1723E-02,
1189 + 0.4122E-02, 0.1723E-02, 0.1469E-01, 0.5689E-04,
1190 + 0.0000E+00, 0.0000E+00, 0.1505E-03, 0.2059E-02,
1191 + 0.4269E-02, 0.2059E-02, 0.1668E-01, 0.6163E-04,
1192 + 0.0000E+00, 0.0000E+00, -0.3748E-03, 0.2402E-02,
1193 + 0.4430E-02, 0.2402E-02, 0.1872E-01, 0.6894E-04,
1194 + 0.0000E+00, 0.0000E+00, -0.8993E-03, 0.2752E-02,
1195 + 0.4608E-02, 0.2752E-02, 0.2082E-01, 0.7800E-04,
1196 + 0.0000E+00, 0.0000E+00, -0.1421E-02, 0.3111E-02,
1197 + 0.4803E-02, 0.3111E-02, 0.2300E-01, 0.8825E-04,
1198 + 0.0000E+00, 0.0000E+00, -0.1944E-02, 0.3480E-02,
1199 + 0.5018E-02, 0.3480E-02, 0.2525E-01, 0.9931E-04/
1200 DATA ((CALC(I,K,27),I=1,8), K=10,18) /
1201 + 0.0000E+00, 0.0000E+00, -0.2467E-02, 0.3859E-02,
1202 + 0.5253E-02, 0.3859E-02, 0.2759E-01, 0.1109E-03,
1203 + 0.0000E+00, 0.0000E+00, -0.2993E-02, 0.4252E-02,
1204 + 0.5513E-02, 0.4252E-02, 0.3003E-01, 0.1230E-03,
1205 + 0.0000E+00, 0.0000E+00, -0.3525E-02, 0.4659E-02,
1206 + 0.5794E-02, 0.4659E-02, 0.3258E-01, 0.1354E-03,
1207 + 0.0000E+00, 0.0000E+00, -0.4064E-02, 0.5083E-02,
1208 + 0.6104E-02, 0.5083E-02, 0.3525E-01, 0.1481E-03,
1209 + 0.0000E+00, 0.0000E+00, -0.4613E-02, 0.5527E-02,
1210 + 0.6442E-02, 0.5527E-02, 0.3807E-01, 0.1612E-03,
1211 + 0.0000E+00, 0.0000E+00, -0.5175E-02, 0.5992E-02,
1212 + 0.6810E-02, 0.5992E-02, 0.4103E-01, 0.1743E-03,
1213 + 0.0000E+00, 0.0000E+00, -0.5753E-02, 0.6482E-02,
1214 + 0.7212E-02, 0.6482E-02, 0.4419E-01, 0.1887E-03,
1215 + 0.0000E+00, 0.0000E+00, -0.6350E-02, 0.6999E-02,
1216 + 0.7648E-02, 0.6999E-02, 0.4753E-01, 0.2038E-03,
1217 + 0.0000E+00, 0.0000E+00, -0.6971E-02, 0.7547E-02,
1218 + 0.8123E-02, 0.7547E-02, 0.5107E-01, 0.2184E-03/
1219 DATA ((CALC(I,K,27),I=1,8), K=19,20) /
1220 + 0.0000E+00, 0.0000E+00, -0.7618E-02, 0.8128E-02,
1221 + 0.8639E-02, 0.8128E-02, 0.5487E-01, 0.2347E-03,
1222 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
1223 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
1224 DATA ((CALC(I,K,28),I=1,8), K=1,9) /
1225 + 0.0000E+00, 0.0000E+00, 0.2488E-02, 0.8293E-03,
1226 + 0.4148E-02, 0.8293E-03, 0.9954E-02, 0.4417E-04,
1227 + 0.0000E+00, 0.0000E+00, 0.1839E-02, 0.1197E-02,
1228 + 0.4234E-02, 0.1197E-02, 0.1203E-01, 0.3802E-04,
1229 + 0.0000E+00, 0.0000E+00, 0.1204E-02, 0.1564E-02,
1230 + 0.4331E-02, 0.1564E-02, 0.1411E-01, 0.3700E-04,
1231 + 0.0000E+00, 0.0000E+00, 0.5842E-03, 0.1929E-02,
1232 + 0.4442E-02, 0.1929E-02, 0.1621E-01, 0.3926E-04,
1233 + 0.0000E+00, 0.0000E+00, -0.2016E-04, 0.2294E-02,
1234 + 0.4567E-02, 0.2294E-02, 0.1832E-01, 0.4360E-04,
1235 + 0.0000E+00, 0.0000E+00, -0.6087E-03, 0.2660E-02,
1236 + 0.4709E-02, 0.2660E-02, 0.2046E-01, 0.4930E-04,
1237 + 0.0000E+00, 0.0000E+00, -0.1188E-02, 0.3027E-02,
1238 + 0.4869E-02, 0.3027E-02, 0.2264E-01, 0.5583E-04,
1239 + 0.0000E+00, 0.0000E+00, -0.1752E-02, 0.3398E-02,
1240 + 0.5048E-02, 0.3398E-02, 0.2485E-01, 0.6283E-04,
1241 + 0.0000E+00, 0.0000E+00, -0.2306E-02, 0.3774E-02,
1242 + 0.5247E-02, 0.3774E-02, 0.2713E-01, 0.7008E-04/
1243 DATA ((CALC(I,K,28),I=1,8), K=10,18) /
1244 + 0.0000E+00, 0.0000E+00, -0.2851E-02, 0.4158E-02,
1245 + 0.5469E-02, 0.4158E-02, 0.2947E-01, 0.7746E-04,
1246 + 0.0000E+00, 0.0000E+00, -0.3391E-02, 0.4551E-02,
1247 + 0.5712E-02, 0.4551E-02, 0.3189E-01, 0.8492E-04,
1248 + 0.0000E+00, 0.0000E+00, -0.3929E-02, 0.4955E-02,
1249 + 0.5983E-02, 0.4955E-02, 0.3441E-01, 0.9244E-04,
1250 + 0.0000E+00, 0.0000E+00, -0.4468E-02, 0.5373E-02,
1251 + 0.6280E-02, 0.5373E-02, 0.3703E-01, 0.1000E-03,
1252 + 0.0000E+00, 0.0000E+00, -0.5012E-02, 0.5809E-02,
1253 + 0.6603E-02, 0.5809E-02, 0.3978E-01, 0.1078E-03,
1254 + 0.0000E+00, 0.0000E+00, -0.5563E-02, 0.6263E-02,
1255 + 0.6961E-02, 0.6263E-02, 0.4268E-01, 0.1157E-03,
1256 + 0.0000E+00, 0.0000E+00, -0.6127E-02, 0.6740E-02,
1257 + 0.7351E-02, 0.6740E-02, 0.4578E-01, 0.1250E-03,
1258 + 0.0000E+00, 0.0000E+00, -0.6706E-02, 0.7240E-02,
1259 + 0.7776E-02, 0.7240E-02, 0.4901E-01, 0.1327E-03,
1260 + 0.0000E+00, 0.0000E+00, -0.7305E-02, 0.7772E-02,
1261 + 0.8238E-02, 0.7772E-02, 0.5246E-01, 0.1414E-03/
1262 DATA ((CALC(I,K,28),I=1,8), K=19,20) /
1263 + 0.0000E+00, 0.0000E+00, -0.7929E-02, 0.8336E-02,
1264 + 0.8741E-02, 0.8336E-02, 0.5613E-01, 0.1509E-03,
1265 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
1266 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
1267 DATA ((CALC(I,K,29),I=1,8), K=1,9) /
1268 + 0.0000E+00, 0.0000E+00, 0.2758E-02, 0.9195E-03,
1269 + 0.4597E-02, 0.9195E-03, 0.1103E-01, 0.2644E-04,
1270 + 0.0000E+00, 0.0000E+00, 0.1946E-02, 0.1347E-02,
1271 + 0.4639E-02, 0.1347E-02, 0.1337E-01, 0.2369E-04,
1272 + 0.0000E+00, 0.0000E+00, 0.1176E-02, 0.1761E-02,
1273 + 0.4696E-02, 0.1761E-02, 0.1565E-01, 0.2419E-04,
1274 + 0.0000E+00, 0.0000E+00, 0.4450E-03, 0.2163E-02,
1275 + 0.4770E-02, 0.2163E-02, 0.1790E-01, 0.2654E-04,
1276 + 0.0000E+00, 0.0000E+00, -0.2480E-03, 0.2556E-02,
1277 + 0.4863E-02, 0.2556E-02, 0.2012E-01, 0.2988E-04,
1278 + 0.0000E+00, 0.0000E+00, -0.9065E-03, 0.2942E-02,
1279 + 0.4976E-02, 0.2942E-02, 0.2232E-01, 0.3379E-04,
1280 + 0.0000E+00, 0.0000E+00, -0.1552E-02, 0.3321E-02,
1281 + 0.5109E-02, 0.3321E-02, 0.2453E-01, 0.3796E-04,
1282 + 0.0000E+00, 0.0000E+00, -0.2148E-02, 0.3698E-02,
1283 + 0.5264E-02, 0.3698E-02, 0.2674E-01, 0.4218E-04,
1284 + 0.0000E+00, 0.0000E+00, -0.2720E-02, 0.4074E-02,
1285 + 0.5442E-02, 0.4074E-02, 0.2899E-01, 0.4636E-04/
1286 DATA ((CALC(I,K,29),I=1,8), K=10,18) /
1287 + 0.0000E+00, 0.0000E+00, -0.3274E-02, 0.4453E-02,
1288 + 0.5644E-02, 0.4453E-02, 0.3127E-01, 0.5047E-04,
1289 + 0.0000E+00, 0.0000E+00, -0.3813E-02, 0.4836E-02,
1290 + 0.5871E-02, 0.4836E-02, 0.3362E-01, 0.5450E-04,
1291 + 0.0000E+00, 0.0000E+00, -0.4342E-02, 0.5227E-02,
1292 + 0.6123E-02, 0.5227E-02, 0.3604E-01, 0.5847E-04,
1293 + 0.0000E+00, 0.0000E+00, -0.4865E-02, 0.5629E-02,
1294 + 0.6401E-02, 0.5629E-02, 0.3856E-01, 0.6242E-04,
1295 + 0.0000E+00, 0.0000E+00, -0.5387E-02, 0.6045E-02,
1296 + 0.6710E-02, 0.6045E-02, 0.4119E-01, 0.6640E-04,
1297 + 0.0000E+00, 0.0000E+00, -0.5912E-02, 0.6478E-02,
1298 + 0.7047E-02, 0.6478E-02, 0.4391E-01, 0.6982E-04,
1299 + 0.0000E+00, 0.0000E+00, -0.6446E-02, 0.6931E-02,
1300 + 0.7417E-02, 0.6931E-02, 0.4686E-01, 0.7454E-04,
1301 + 0.0000E+00, 0.0000E+00, -0.6991E-02, 0.7403E-02,
1302 + 0.7821E-02, 0.7403E-02, 0.4993E-01, 0.7895E-04,
1303 + 0.0000E+00, 0.0000E+00, -0.7554E-02, 0.7907E-02,
1304 + 0.8260E-02, 0.7907E-02, 0.5320E-01, 0.8349E-04/
1305 DATA ((CALC(I,K,29),I=1,8), K=19,20) /
1306 + 0.0000E+00, 0.0000E+00, -0.8139E-02, 0.8437E-02,
1307 + 0.8739E-02, 0.8437E-02, 0.5667E-01, 0.8814E-04,
1308 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
1309 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
1310 DATA ((CALC(I,K,30),I=1,8), K=1,9) /
1311 + 0.0000E+00, 0.0000E+00, 0.3081E-02, 0.1028E-02,
1312 + 0.5135E-02, 0.1028E-02, 0.1233E-01, 0.1487E-04,
1313 + 0.0000E+00, 0.0000E+00, 0.2042E-02, 0.1531E-02,
1314 + 0.5103E-02, 0.1531E-02, 0.1497E-01, 0.1400E-04,
1315 + 0.0000E+00, 0.0000E+00, 0.1093E-02, 0.2002E-02,
1316 + 0.5097E-02, 0.2002E-02, 0.1747E-01, 0.1495E-04,
1317 + 0.0000E+00, 0.0000E+00, 0.2250E-03, 0.2446E-02,
1318 + 0.5116E-02, 0.2446E-02, 0.1987E-01, 0.1677E-04,
1319 + 0.0000E+00, 0.0000E+00, -0.5694E-03, 0.2867E-02,
1320 + 0.5162E-02, 0.2867E-02, 0.2217E-01, 0.1889E-04,
1321 + 0.0000E+00, 0.0000E+00, -0.1300E-02, 0.3267E-02,
1322 + 0.5234E-02, 0.3267E-02, 0.2440E-01, 0.2112E-04,
1323 + 0.0000E+00, 0.0000E+00, -0.1974E-02, 0.3653E-02,
1324 + 0.5332E-02, 0.3653E-02, 0.2659E-01, 0.2333E-04,
1325 + 0.0000E+00, 0.0000E+00, -0.2597E-02, 0.4027E-02,
1326 + 0.5456E-02, 0.4027E-02, 0.2875E-01, 0.2542E-04,
1327 + 0.0000E+00, 0.0000E+00, -0.3183E-02, 0.4393E-02,
1328 + 0.5606E-02, 0.4393E-02, 0.3090E-01, 0.2739E-04/
1329 DATA ((CALC(I,K,30),I=1,8), K=10,18) /
1330 + 0.0000E+00, 0.0000E+00, -0.3733E-02, 0.4756E-02,
1331 + 0.5782E-02, 0.4756E-02, 0.3307E-01, 0.2923E-04,
1332 + 0.0000E+00, 0.0000E+00, -0.4257E-02, 0.5118E-02,
1333 + 0.5984E-02, 0.5118E-02, 0.3527E-01, 0.3097E-04,
1334 + 0.0000E+00, 0.0000E+00, -0.4761E-02, 0.5485E-02,
1335 + 0.6207E-02, 0.5485E-02, 0.3753E-01, 0.3263E-04,
1336 + 0.0000E+00, 0.0000E+00, -0.5252E-02, 0.5858E-02,
1337 + 0.6461E-02, 0.5858E-02, 0.3986E-01, 0.3426E-04,
1338 + 0.0000E+00, 0.0000E+00, -0.5735E-02, 0.6241E-02,
1339 + 0.6747E-02, 0.6241E-02, 0.4229E-01, 0.3587E-04,
1340 + 0.0000E+00, 0.0000E+00, -0.6217E-02, 0.6640E-02,
1341 + 0.7059E-02, 0.6640E-02, 0.4481E-01, 0.3728E-04,
1342 + 0.0000E+00, 0.0000E+00, -0.6703E-02, 0.7055E-02,
1343 + 0.7402E-02, 0.7055E-02, 0.4751E-01, 0.3921E-04,
1344 + 0.0000E+00, 0.0000E+00, -0.7198E-02, 0.7475E-02,
1345 + 0.7777E-02, 0.7475E-02, 0.5034E-01, 0.4109E-04,
1346 + 0.0000E+00, 0.0000E+00, -0.7706E-02, 0.7940E-02,
1347 + 0.8185E-02, 0.7940E-02, 0.5336E-01, 0.4273E-04/
1348 DATA ((CALC(I,K,30),I=1,8), K=19,20) /
1349 + 0.0000E+00, 0.0000E+00, -0.8242E-02, 0.8430E-02,
1350 + 0.8630E-02, 0.8430E-02, 0.5653E-01, 0.4452E-04,
1351 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
1352 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
1353 DATA ((CALC(I,K,31),I=1,8), K=1,9) /
1354 + 0.0000E+00, 0.0000E+00, 0.3507E-02, 0.1169E-02,
1355 + 0.5847E-02, 0.1169E-02, 0.1403E-01, 0.7133E-05,
1356 + 0.0000E+00, 0.0000E+00, 0.2113E-02, 0.1787E-02,
1357 + 0.5686E-02, 0.1787E-02, 0.1711E-01, 0.7043E-05,
1358 + 0.0000E+00, 0.0000E+00, 0.9042E-03, 0.2337E-02,
1359 + 0.5577E-02, 0.2337E-02, 0.1990E-01, 0.7782E-05,
1360 + 0.0000E+00, 0.0000E+00, -0.1506E-03, 0.2833E-02,
1361 + 0.5513E-02, 0.2833E-02, 0.2246E-01, 0.8803E-05,
1362 + 0.0000E+00, 0.0000E+00, -0.1071E-02, 0.3282E-02,
1363 + 0.5492E-02, 0.3282E-02, 0.2482E-01, 0.9785E-05,
1364 + 0.0000E+00, 0.0000E+00, -0.1889E-02, 0.3693E-02,
1365 + 0.5509E-02, 0.3693E-02, 0.2704E-01, 0.1070E-04,
1366 + 0.0000E+00, 0.0000E+00, -0.2617E-02, 0.4075E-02,
1367 + 0.5562E-02, 0.4075E-02, 0.2915E-01, 0.1153E-04,
1368 + 0.0000E+00, 0.0000E+00, -0.3237E-02, 0.4435E-02,
1369 + 0.5648E-02, 0.4435E-02, 0.3118E-01, 0.1225E-04,
1370 + 0.0000E+00, 0.0000E+00, -0.3810E-02, 0.4777E-02,
1371 + 0.5764E-02, 0.4777E-02, 0.3317E-01, 0.1286E-04/
1372 DATA ((CALC(I,K,31),I=1,8), K=10,18) /
1373 + 0.0000E+00, 0.0000E+00, -0.4327E-02, 0.5110E-02,
1374 + 0.5910E-02, 0.5110E-02, 0.3513E-01, 0.1339E-04,
1375 + 0.0000E+00, 0.0000E+00, -0.4804E-02, 0.5436E-02,
1376 + 0.6084E-02, 0.5436E-02, 0.3711E-01, 0.1387E-04,
1377 + 0.0000E+00, 0.0000E+00, -0.5254E-02, 0.5762E-02,
1378 + 0.6277E-02, 0.5762E-02, 0.3911E-01, 0.1430E-04,
1379 + 0.0000E+00, 0.0000E+00, -0.5684E-02, 0.6092E-02,
1380 + 0.6501E-02, 0.6092E-02, 0.4117E-01, 0.1470E-04,
1381 + 0.0000E+00, 0.0000E+00, -0.6102E-02, 0.6428E-02,
1382 + 0.6764E-02, 0.6428E-02, 0.4330E-01, 0.1511E-04,
1383 + 0.0000E+00, 0.0000E+00, -0.6516E-02, 0.6777E-02,
1384 + 0.7036E-02, 0.6777E-02, 0.4552E-01, 0.1546E-04,
1385 + 0.0000E+00, 0.0000E+00, -0.6931E-02, 0.7140E-02,
1386 + 0.7343E-02, 0.7140E-02, 0.4759E-01, 0.1532E-04,
1387 + 0.0000E+00, 0.0000E+00, -0.7349E-02, 0.7518E-02,
1388 + 0.7679E-02, 0.7518E-02, 0.5011E-01, 0.1594E-04,
1389 + 0.0000E+00, 0.0000E+00, -0.7783E-02, 0.7926E-02,
1390 + 0.8045E-02, 0.7926E-02, 0.5287E-01, 0.1677E-04/
1391 DATA ((CALC(I,K,31),I=1,8), K=19,20) /
1392 + 0.0000E+00, 0.0000E+00, -0.8237E-02, 0.8375E-02,
1393 + 0.8443E-02, 0.8375E-02, 0.5570E-01, 0.1729E-04,
1394 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
1395 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
1396 DATA ((CALC(I,K,32),I=1,8), K=1,9) /
1397 + 0.0000E+00, 0.0000E+00, 0.4201E-02, 0.1401E-02,
1398 + 0.7002E-02, 0.1401E-02, 0.1681E-01, 0.1987E-05,
1399 + 0.0000E+00, 0.0000E+00, 0.2060E-02, 0.2256E-02,
1400 + 0.6566E-02, 0.2256E-02, 0.2079E-01, 0.2022E-05,
1401 + 0.0000E+00, 0.0000E+00, 0.3736E-03, 0.2950E-02,
1402 + 0.6263E-02, 0.2950E-02, 0.2408E-01, 0.2281E-05,
1403 + 0.0000E+00, 0.0000E+00, -0.9719E-03, 0.3522E-02,
1404 + 0.6064E-02, 0.3522E-02, 0.2687E-01, 0.2574E-05,
1405 + 0.0000E+00, 0.0000E+00, -0.2051E-02, 0.4003E-02,
1406 + 0.5947E-02, 0.4003E-02, 0.2928E-01, 0.2776E-05,
1407 + 0.0000E+00, 0.0000E+00, -0.2860E-02, 0.4414E-02,
1408 + 0.5898E-02, 0.4414E-02, 0.3140E-01, 0.2928E-05,
1409 + 0.0000E+00, 0.0000E+00, -0.3754E-02, 0.4773E-02,
1410 + 0.5902E-02, 0.4773E-02, 0.3332E-01, 0.3058E-05,
1411 + 0.0000E+00, 0.0000E+00, -0.4339E-02, 0.5094E-02,
1412 + 0.5952E-02, 0.5094E-02, 0.3510E-01, 0.3146E-05,
1413 + 0.0000E+00, 0.0000E+00, -0.4796E-02, 0.5389E-02,
1414 + 0.6041E-02, 0.5389E-02, 0.3679E-01, 0.3195E-05/
1415 DATA ((CALC(I,K,32),I=1,8), K=10,18) /
1416 + 0.0000E+00, 0.0000E+00, -0.5223E-02, 0.5666E-02,
1417 + 0.6163E-02, 0.5666E-02, 0.3843E-01, 0.3228E-05,
1418 + 0.0000E+00, 0.0000E+00, -0.5601E-02, 0.5933E-02,
1419 + 0.6312E-02, 0.5933E-02, 0.4005E-01, 0.3250E-05,
1420 + 0.0000E+00, 0.0000E+00, -0.5946E-02, 0.6197E-02,
1421 + 0.6300E-02, 0.6197E-02, 0.4168E-01, 0.3262E-05,
1422 + 0.0000E+00, 0.0000E+00, -0.6272E-02, 0.6461E-02,
1423 + 0.6473E-02, 0.6461E-02, 0.4336E-01, 0.3272E-05,
1424 + 0.0000E+00, 0.0000E+00, -0.6588E-02, 0.6732E-02,
1425 + 0.6751E-02, 0.6732E-02, 0.4510E-01, 0.3289E-05,
1426 + 0.0000E+00, 0.0000E+00, -0.6902E-02, 0.7021E-02,
1427 + 0.7020E-02, 0.7021E-02, 0.4681E-01, 0.3209E-05,
1428 + 0.0000E+00, 0.0000E+00, -0.7219E-02, 0.7310E-02,
1429 + 0.7297E-02, 0.7310E-02, 0.4882E-01, 0.3295E-05,
1430 + 0.0000E+00, 0.0000E+00, -0.7523E-02, 0.7606E-02,
1431 + 0.7596E-02, 0.7606E-02, 0.5079E-01, 0.3324E-05,
1432 + 0.0000E+00, 0.0000E+00, -0.7851E-02, 0.7922E-02,
1433 + 0.7916E-02, 0.7922E-02, 0.5287E-01, 0.3328E-05/
1434 DATA ((CALC(I,K,32),I=1,8), K=19,20) /
1435 + 0.0000E+00, 0.0000E+00, -0.8208E-02, 0.8263E-02,
1436 + 0.8256E-02, 0.8263E-02, 0.5513E-01, 0.3346E-05,
1437 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,
1438 + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/
1439C----------------------------------------------------------------------
1440 DATA ISTART/0/
1441 SAVE ISTART,OWLAM2,Q02,FLAV, /W5051I7/
1442C
1443 IF (ISTART.EQ.0) THEN
1444 ISTART=1
1445 DO 10 K=1,32
1446 DO 10 I=1,20
1447 DO 10 M=1,8
1448 10 CALCO(M,I,K) = CALC(M,I,K)
1449 OWLAM=PAR(1)
1450 OWLAM2=OWLAM**2
1451 Q02=PAR(30)
1452 FLAV=PAR(25)
1453 DELTA=PAR(29)
1454 CALL WATE32
1455 ENDIF
1456C
1457 X = DX
1458 Q = DQ
1459 Q2 = Q*Q
1460 IDQ2=2
1461 SB=0.
1462 IF(Q2-Q02) 1,1,2
1463 2 IF(IDQ2-1) 1,1,3
1464 3 SB= LOG( LOG( MAX(Q02,Q2)/OWLAM2)/ LOG(Q02/OWLAM2))
1465 1 CONTINUE
1466 CALL AURGAM(8,0,X,SB,XPDF(7))
1467 CALL AURGAM(7,0,X,SB,SING)
1468 CALL AURGAM(4,0,X,SB,DPLUSNS)
1469 CALL AURGAM(3,0,X,SB,CPLUSNS)
1470 CALL AURGAM(5,0,X,SB,UPLUSNS)
1471 CALL AURGAM(6,0,X,SB,SPLUSNS)
1472 XPDF(3) = CPLUSNS
1473 XPDF(4) = DPLUSNS
1474 XPDF(5) = UPLUSNS
1475 XPDF(6) = SPLUSNS
1476 XPDF(1) = SING
1477C
1478 ADD = XPDF(1)/FLAV
1479 UPLUS=XPDF(5)+ADD
1480 DPLUS=-XPDF(4)+ADD
1481 SPLUS=-XPDF(6)+ADD
1482 CPLUS=-XPDF(3)+ADD
1483 UB=UPLUS*0.5
1484 UV=UB
1485 DB=DPLUS*0.5
1486 DV=DB
1487 SB=SPLUS*0.5
1488 CB=CPLUS*0.5
1489 SING=XPDF(1)
1490 GLU=XPDF(7)
1491 GL=GLU
1492C... get parton density with massive charm
1493 CPLUM=AFCPLU(X,Q2)
1494 CB=CPLUM*0.5
1495*
1496 DUV=MAX(ZERO,UV)
1497 DDV=MAX(ZERO,DV)
1498 DUB=MAX(ZERO,UB)
1499 DDB=MAX(ZERO,DB)
1500 DSB=MAX(ZERO,SB)
1501 DCB=MAX(ZERO,CB)
1502 DGL=MAX(ZERO,GL)
1503C
1504 RETURN
1505 END