* * $Id$ * * $Log$ * Revision 1.1.1.3 1996/11/07 11:52:38 cernlib * Import version 7.05 * * Revision 1.3 1996/11/07 09:35:02 pdflib * New version 7.05: cut routine Ctq4RTb in 9 pieces * * Revision 1.1.1.1 1996/04/12 15:28:56 plothow * Version 7.01 * * #include "pdf/pilot.h" SUBROUTINE ACFGP1(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DGL) C C INTERPOLATION PROGRAM WHICH INTERPOLATES THE GRID "DATAGA" AND GIVES THE C QUARK AND GLUON DISTRIBUTIONS IN THE REAL PHOTON, AS FUNCTIONS OF X AND Q2 C C THE Q2-EVOLUTION IS PERFORMED WITH BLL AP-EQUATIONS AND NF=4. A MASSIVE C CHARM DISTRIBUTION (BORROWED FROM GLUCK AND REYA) IS ALSO AVAILABLE. C C THE BOUNDARY CONDITIONS ARE SUCH THAT THE DISTRIBUTION FUNCTIONS ARE GIVEN C BY A VDM "ANSATZ" AT Q2=.25 GEV**2. C C THE PROGRAM WORKS FOR 2. GEV**2 < Q2 <5.5E+5 AND .00137 < X < .9986 C C THE DISTRIBUTIONS ARE CALCULATED IN THE MSBAR FACTORIZATION SCHEME. C C THE VALUE OF LAMBDA-MSB IS 200 MEV C C THE OUTPUT IS WRITTEN IN THE FILE 'FILEOUT': C X*U=X*U(X,Q2) C X*D= ... C X*S= ... C X*C= ... (MASSLESS CHARM WITH C(X,2.)=0) C X*CM= ... (MASSIVE CHARM WITH MC=1.5 GEV ) C X*GLU=GLUON(X,Q2)*X C C C F2 = PHOTON STRUCTURE FUNCTION WITHOUT CHARM C F2C= " " " WITH MASSIVE CHARM C #include "pdf/expdp.inc" + DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL REAL X, Q, UV, DV, UB, DB, SB, CB, BB, GL REAL Q2 DIMENSION XPDF(7),PAR(30),CALC(8,20,32),CALCO(8,20,32) COMMON/W5051I7/CALCO EXTERNAL AFCPLU DATA ZERO/0.0/ C---------------------------------------------------------------------- DATA (PAR(I),I=1,30) / + 0.2000E+00, 0.4800E+00, 0.8500E+00, 0.1000E+01, + 0.1000E+01, 0.1000E+01, 0.1200E+01, 0.7500E+01, + 0.1987E+01, 0.1900E+01, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.1000E+01, + 0.0000E+00, 0.1000E+01, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.4000E+01, 0.0000E+00, 0.1000E+01, 0.1000E+01, + 0.8000E-01, 0.2000E+01/ DATA ((CALC(I,K,1 ),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1671E-02, -0.5026E-03, + 0.6657E-03, -0.5026E-03, 0.6683E-02, 0.9320E-02, + 0.0000E+00, 0.0000E+00, 0.1669E-02, -0.4992E-03, + 0.6709E-03, -0.4992E-03, 0.7711E-02, 0.1350E-01, + 0.0000E+00, 0.0000E+00, 0.1667E-02, -0.4955E-03, + 0.6764E-03, -0.4955E-03, 0.8915E-02, 0.1832E-01, + 0.0000E+00, 0.0000E+00, 0.1665E-02, -0.4917E-03, + 0.6822E-03, -0.4917E-03, 0.1031E-01, 0.2383E-01, + 0.0000E+00, 0.0000E+00, 0.1663E-02, -0.4875E-03, + 0.6883E-03, -0.4875E-03, 0.1192E-01, 0.3008E-01, + 0.0000E+00, 0.0000E+00, 0.1661E-02, -0.4831E-03, + 0.6947E-03, -0.4831E-03, 0.1375E-01, 0.3713E-01, + 0.0000E+00, 0.0000E+00, 0.1658E-02, -0.4784E-03, + 0.7014E-03, -0.4784E-03, 0.1583E-01, 0.4503E-01, + 0.0000E+00, 0.0000E+00, 0.1655E-02, -0.4733E-03, + 0.7085E-03, -0.4733E-03, 0.1817E-01, 0.5383E-01, + 0.0000E+00, 0.0000E+00, 0.1651E-02, -0.4679E-03, + 0.7159E-03, -0.4679E-03, 0.2079E-01, 0.6360E-01/ DATA ((CALC(I,K,1 ),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, 0.1647E-02, -0.4621E-03, + 0.7238E-03, -0.4621E-03, 0.2373E-01, 0.7441E-01, + 0.0000E+00, 0.0000E+00, 0.1643E-02, -0.4558E-03, + 0.7321E-03, -0.4558E-03, 0.2700E-01, 0.8633E-01, + 0.0000E+00, 0.0000E+00, 0.1639E-02, -0.4492E-03, + 0.7409E-03, -0.4492E-03, 0.3062E-01, 0.9943E-01, + 0.0000E+00, 0.0000E+00, 0.1633E-02, -0.4420E-03, + 0.7502E-03, -0.4420E-03, 0.3464E-01, 0.1138E+00, + 0.0000E+00, 0.0000E+00, 0.1628E-02, -0.4343E-03, + 0.7600E-03, -0.4343E-03, 0.3907E-01, 0.1295E+00, + 0.0000E+00, 0.0000E+00, 0.1622E-02, -0.4260E-03, + 0.7704E-03, -0.4260E-03, 0.4394E-01, 0.1467E+00, + 0.0000E+00, 0.0000E+00, 0.1615E-02, -0.4172E-03, + 0.7815E-03, -0.4172E-03, 0.4931E-01, 0.1655E+00, + 0.0000E+00, 0.0000E+00, 0.1608E-02, -0.4076E-03, + 0.7932E-03, -0.4076E-03, 0.5520E-01, 0.1859E+00, + 0.0000E+00, 0.0000E+00, 0.1599E-02, -0.3973E-03, + 0.8057E-03, -0.3973E-03, 0.6165E-01, 0.2082E+00/ DATA ((CALC(I,K,1 ),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, 0.1591E-02, -0.3863E-03, + 0.8189E-03, -0.3863E-03, 0.6871E-01, 0.2324E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,2 ),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1675E-02, -0.4282E-03, + 0.8187E-03, -0.4282E-03, 0.6700E-02, 0.8422E-02, + 0.0000E+00, 0.0000E+00, 0.1668E-02, -0.4181E-03, + 0.8319E-03, -0.4181E-03, 0.7379E-02, 0.1103E-01, + 0.0000E+00, 0.0000E+00, 0.1661E-02, -0.4075E-03, + 0.8457E-03, -0.4075E-03, 0.8140E-02, 0.1387E-01, + 0.0000E+00, 0.0000E+00, 0.1653E-02, -0.3961E-03, + 0.8603E-03, -0.3961E-03, 0.8986E-02, 0.1696E-01, + 0.0000E+00, 0.0000E+00, 0.1644E-02, -0.3838E-03, + 0.8757E-03, -0.3838E-03, 0.9923E-02, 0.2029E-01, + 0.0000E+00, 0.0000E+00, 0.1634E-02, -0.3707E-03, + 0.8919E-03, -0.3707E-03, 0.1096E-01, 0.2388E-01, + 0.0000E+00, 0.0000E+00, 0.1623E-02, -0.3567E-03, + 0.9091E-03, -0.3567E-03, 0.1209E-01, 0.2775E-01, + 0.0000E+00, 0.0000E+00, 0.1611E-02, -0.3416E-03, + 0.9272E-03, -0.3416E-03, 0.1332E-01, 0.3192E-01, + 0.0000E+00, 0.0000E+00, 0.1598E-02, -0.3254E-03, + 0.9465E-03, -0.3254E-03, 0.1467E-01, 0.3639E-01/ DATA ((CALC(I,K,2 ),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, 0.1584E-02, -0.3080E-03, + 0.9669E-03, -0.3080E-03, 0.1614E-01, 0.4119E-01, + 0.0000E+00, 0.0000E+00, 0.1568E-02, -0.2893E-03, + 0.9885E-03, -0.2893E-03, 0.1773E-01, 0.4634E-01, + 0.0000E+00, 0.0000E+00, 0.1551E-02, -0.2691E-03, + 0.1012E-02, -0.2691E-03, 0.1946E-01, 0.5188E-01, + 0.0000E+00, 0.0000E+00, 0.1532E-02, -0.2475E-03, + 0.1036E-02, -0.2475E-03, 0.2133E-01, 0.5782E-01, + 0.0000E+00, 0.0000E+00, 0.1512E-02, -0.2242E-03, + 0.1062E-02, -0.2242E-03, 0.2335E-01, 0.6421E-01, + 0.0000E+00, 0.0000E+00, 0.1490E-02, -0.1991E-03, + 0.1090E-02, -0.1991E-03, 0.2554E-01, 0.7109E-01, + 0.0000E+00, 0.0000E+00, 0.1465E-02, -0.1721E-03, + 0.1120E-02, -0.1721E-03, 0.2790E-01, 0.7848E-01, + 0.0000E+00, 0.0000E+00, 0.1439E-02, -0.1429E-03, + 0.1151E-02, -0.1429E-03, 0.3046E-01, 0.8645E-01, + 0.0000E+00, 0.0000E+00, 0.1410E-02, -0.1115E-03, + 0.1185E-02, -0.1115E-03, 0.3321E-01, 0.9504E-01/ DATA ((CALC(I,K,2 ),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, 0.1379E-02, -0.7764E-04, + 0.1222E-02, -0.7764E-04, 0.3619E-01, 0.1043E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,3 ),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1659E-02, -0.3433E-03, + 0.9723E-03, -0.3433E-03, 0.6635E-02, 0.7843E-02, + 0.0000E+00, 0.0000E+00, 0.1645E-02, -0.3259E-03, + 0.9931E-03, -0.3259E-03, 0.7163E-02, 0.9673E-02, + 0.0000E+00, 0.0000E+00, 0.1630E-02, -0.3073E-03, + 0.1015E-02, -0.3073E-03, 0.7742E-02, 0.1161E-01, + 0.0000E+00, 0.0000E+00, 0.1613E-02, -0.2873E-03, + 0.1038E-02, -0.2873E-03, 0.8372E-02, 0.1364E-01, + 0.0000E+00, 0.0000E+00, 0.1594E-02, -0.2658E-03, + 0.1063E-02, -0.2658E-03, 0.9056E-02, 0.1579E-01, + 0.0000E+00, 0.0000E+00, 0.1574E-02, -0.2428E-03, + 0.1089E-02, -0.2428E-03, 0.9796E-02, 0.1806E-01, + 0.0000E+00, 0.0000E+00, 0.1553E-02, -0.2181E-03, + 0.1116E-02, -0.2181E-03, 0.1060E-01, 0.2046E-01, + 0.0000E+00, 0.0000E+00, 0.1529E-02, -0.1915E-03, + 0.1146E-02, -0.1915E-03, 0.1146E-01, 0.2299E-01, + 0.0000E+00, 0.0000E+00, 0.1503E-02, -0.1629E-03, + 0.1177E-02, -0.1629E-03, 0.1239E-01, 0.2567E-01/ DATA ((CALC(I,K,3 ),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, 0.1475E-02, -0.1321E-03, + 0.1210E-02, -0.1321E-03, 0.1339E-01, 0.2852E-01, + 0.0000E+00, 0.0000E+00, 0.1444E-02, -0.9894E-04, + 0.1246E-02, -0.9894E-04, 0.1446E-01, 0.3155E-01, + 0.0000E+00, 0.0000E+00, 0.1410E-02, -0.6324E-04, + 0.1284E-02, -0.6324E-04, 0.1561E-01, 0.3478E-01, + 0.0000E+00, 0.0000E+00, 0.1374E-02, -0.2477E-04, + 0.1324E-02, -0.2477E-04, 0.1686E-01, 0.3822E-01, + 0.0000E+00, 0.0000E+00, 0.1335E-02, 0.1669E-04, + 0.1368E-02, 0.1669E-04, 0.1819E-01, 0.4191E-01, + 0.0000E+00, 0.0000E+00, 0.1292E-02, 0.6140E-04, + 0.1414E-02, 0.6140E-04, 0.1962E-01, 0.4586E-01, + 0.0000E+00, 0.0000E+00, 0.1245E-02, 0.1096E-03, + 0.1464E-02, 0.1096E-03, 0.2117E-01, 0.5011E-01, + 0.0000E+00, 0.0000E+00, 0.1194E-02, 0.1617E-03, + 0.1517E-02, 0.1617E-03, 0.2283E-01, 0.5469E-01, + 0.0000E+00, 0.0000E+00, 0.1140E-02, 0.2178E-03, + 0.1575E-02, 0.2178E-03, 0.2462E-01, 0.5962E-01/ DATA ((CALC(I,K,3 ),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, 0.1080E-02, 0.2785E-03, + 0.1637E-02, 0.2785E-03, 0.2655E-01, 0.6494E-01, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,4 ),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1625E-02, -0.2514E-03, + 0.1122E-02, -0.2514E-03, 0.6500E-02, 0.7358E-02, + 0.0000E+00, 0.0000E+00, 0.1602E-02, -0.2260E-03, + 0.1150E-02, -0.2260E-03, 0.6948E-02, 0.8685E-02, + 0.0000E+00, 0.0000E+00, 0.1578E-02, -0.1988E-03, + 0.1180E-02, -0.1988E-03, 0.7432E-02, 0.1006E-01, + 0.0000E+00, 0.0000E+00, 0.1551E-02, -0.1696E-03, + 0.1212E-02, -0.1696E-03, 0.7953E-02, 0.1148E-01, + 0.0000E+00, 0.0000E+00, 0.1522E-02, -0.1382E-03, + 0.1245E-02, -0.1382E-03, 0.8512E-02, 0.1296E-01, + 0.0000E+00, 0.0000E+00, 0.1490E-02, -0.1045E-03, + 0.1281E-02, -0.1045E-03, 0.9113E-02, 0.1450E-01, + 0.0000E+00, 0.0000E+00, 0.1456E-02, -0.6824E-04, + 0.1319E-02, -0.6824E-04, 0.9757E-02, 0.1610E-01, + 0.0000E+00, 0.0000E+00, 0.1419E-02, -0.2924E-04, + 0.1360E-02, -0.2924E-04, 0.1045E-01, 0.1779E-01, + 0.0000E+00, 0.0000E+00, 0.1378E-02, 0.1275E-04, + 0.1403E-02, 0.1275E-04, 0.1119E-01, 0.1957E-01/ DATA ((CALC(I,K,4 ),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, 0.1334E-02, 0.5797E-04, + 0.1450E-02, 0.5797E-04, 0.1198E-01, 0.2145E-01, + 0.0000E+00, 0.0000E+00, 0.1287E-02, 0.1067E-03, + 0.1499E-02, 0.1067E-03, 0.1284E-01, 0.2345E-01, + 0.0000E+00, 0.0000E+00, 0.1235E-02, 0.1592E-03, + 0.1553E-02, 0.1592E-03, 0.1375E-01, 0.2557E-01, + 0.0000E+00, 0.0000E+00, 0.1179E-02, 0.2158E-03, + 0.1610E-02, 0.2158E-03, 0.1474E-01, 0.2784E-01, + 0.0000E+00, 0.0000E+00, 0.1118E-02, 0.2769E-03, + 0.1671E-02, 0.2769E-03, 0.1579E-01, 0.3028E-01, + 0.0000E+00, 0.0000E+00, 0.1052E-02, 0.3428E-03, + 0.1737E-02, 0.3428E-03, 0.1693E-01, 0.3289E-01, + 0.0000E+00, 0.0000E+00, 0.9810E-03, 0.4140E-03, + 0.1808E-02, 0.4140E-03, 0.1815E-01, 0.3571E-01, + 0.0000E+00, 0.0000E+00, 0.9038E-03, 0.4909E-03, + 0.1885E-02, 0.4909E-03, 0.1947E-01, 0.3875E-01, + 0.0000E+00, 0.0000E+00, 0.8201E-03, 0.5739E-03, + 0.1967E-02, 0.5739E-03, 0.2090E-01, 0.4204E-01/ DATA ((CALC(I,K,4 ),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, 0.7295E-03, 0.6636E-03, + 0.2056E-02, 0.6636E-03, 0.2243E-01, 0.4560E-01, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,5 ),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1572E-02, -0.1542E-03, + 0.1264E-02, -0.1542E-03, 0.6288E-02, 0.6899E-02, + 0.0000E+00, 0.0000E+00, 0.1541E-02, -0.1209E-03, + 0.1299E-02, -0.1209E-03, 0.6695E-02, 0.7867E-02, + 0.0000E+00, 0.0000E+00, 0.1507E-02, -0.8520E-04, + 0.1336E-02, -0.8520E-04, 0.7132E-02, 0.8857E-02, + 0.0000E+00, 0.0000E+00, 0.1470E-02, -0.4684E-04, + 0.1375E-02, -0.4684E-04, 0.7599E-02, 0.9865E-02, + 0.0000E+00, 0.0000E+00, 0.1432E-02, -0.5605E-05, + 0.1418E-02, -0.5605E-05, 0.8099E-02, 0.1091E-01, + 0.0000E+00, 0.0000E+00, 0.1388E-02, 0.3874E-04, + 0.1463E-02, 0.3874E-04, 0.8635E-02, 0.1199E-01, + 0.0000E+00, 0.0000E+00, 0.1341E-02, 0.8647E-04, + 0.1511E-02, 0.8647E-04, 0.9210E-02, 0.1311E-01, + 0.0000E+00, 0.0000E+00, 0.1289E-02, 0.1379E-03, + 0.1563E-02, 0.1379E-03, 0.9825E-02, 0.1429E-01, + 0.0000E+00, 0.0000E+00, 0.1234E-02, 0.1932E-03, + 0.1618E-02, 0.1932E-03, 0.1049E-01, 0.1553E-01/ DATA ((CALC(I,K,5 ),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, 0.1174E-02, 0.2529E-03, + 0.1677E-02, 0.2529E-03, 0.1119E-01, 0.1684E-01, + 0.0000E+00, 0.0000E+00, 0.1109E-02, 0.3173E-03, + 0.1741E-02, 0.3173E-03, 0.1196E-01, 0.1824E-01, + 0.0000E+00, 0.0000E+00, 0.1038E-02, 0.3867E-03, + 0.1810E-02, 0.3867E-03, 0.1277E-01, 0.1973E-01, + 0.0000E+00, 0.0000E+00, 0.9620E-03, 0.4616E-03, + 0.1883E-02, 0.4616E-03, 0.1366E-01, 0.2133E-01, + 0.0000E+00, 0.0000E+00, 0.8794E-03, 0.5425E-03, + 0.1963E-02, 0.5425E-03, 0.1461E-01, 0.2305E-01, + 0.0000E+00, 0.0000E+00, 0.7902E-03, 0.6299E-03, + 0.2048E-02, 0.6299E-03, 0.1563E-01, 0.2492E-01, + 0.0000E+00, 0.0000E+00, 0.6935E-03, 0.7243E-03, + 0.2141E-02, 0.7243E-03, 0.1674E-01, 0.2694E-01, + 0.0000E+00, 0.0000E+00, 0.5890E-03, 0.8263E-03, + 0.2240E-02, 0.8263E-03, 0.1794E-01, 0.2912E-01, + 0.0000E+00, 0.0000E+00, 0.4759E-03, 0.9366E-03, + 0.2348E-02, 0.9366E-03, 0.1924E-01, 0.3150E-01/ DATA ((CALC(I,K,5 ),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, 0.3534E-03, 0.1056E-02, + 0.2464E-02, 0.1056E-02, 0.2064E-01, 0.3408E-01, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,6 ),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1509E-02, -0.5692E-04, + 0.1395E-02, -0.5692E-04, 0.6035E-02, 0.6441E-02, + 0.0000E+00, 0.0000E+00, 0.1469E-02, -0.1625E-04, + 0.1436E-02, -0.1625E-04, 0.6426E-02, 0.7136E-02, + 0.0000E+00, 0.0000E+00, 0.1425E-02, 0.2742E-04, + 0.1480E-02, 0.2742E-04, 0.6844E-02, 0.7855E-02, + 0.0000E+00, 0.0000E+00, 0.1378E-02, 0.7439E-04, + 0.1527E-02, 0.7439E-04, 0.7291E-02, 0.8574E-02, + 0.0000E+00, 0.0000E+00, 0.1328E-02, 0.1249E-03, + 0.1577E-02, 0.1249E-03, 0.7770E-02, 0.9307E-02, + 0.0000E+00, 0.0000E+00, 0.1272E-02, 0.1793E-03, + 0.1631E-02, 0.1793E-03, 0.8285E-02, 0.1009E-01, + 0.0000E+00, 0.0000E+00, 0.1213E-02, 0.2380E-03, + 0.1689E-02, 0.2380E-03, 0.8836E-02, 0.1087E-01, + 0.0000E+00, 0.0000E+00, 0.1148E-02, 0.3012E-03, + 0.1751E-02, 0.3012E-03, 0.9427E-02, 0.1171E-01, + 0.0000E+00, 0.0000E+00, 0.1079E-02, 0.3694E-03, + 0.1817E-02, 0.3694E-03, 0.1006E-01, 0.1259E-01/ DATA ((CALC(I,K,6 ),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, 0.1003E-02, 0.4429E-03, + 0.1889E-02, 0.4429E-03, 0.1075E-01, 0.1353E-01, + 0.0000E+00, 0.0000E+00, 0.9214E-03, 0.5223E-03, + 0.1966E-02, 0.5223E-03, 0.1149E-01, 0.1454E-01, + 0.0000E+00, 0.0000E+00, 0.8329E-03, 0.6081E-03, + 0.2049E-02, 0.6081E-03, 0.1229E-01, 0.1563E-01, + 0.0000E+00, 0.0000E+00, 0.7372E-03, 0.7007E-03, + 0.2139E-02, 0.7007E-03, 0.1315E-01, 0.1681E-01, + 0.0000E+00, 0.0000E+00, 0.6337E-03, 0.8008E-03, + 0.2235E-02, 0.8008E-03, 0.1408E-01, 0.1809E-01, + 0.0000E+00, 0.0000E+00, 0.5217E-03, 0.9090E-03, + 0.2340E-02, 0.9090E-03, 0.1509E-01, 0.1948E-01, + 0.0000E+00, 0.0000E+00, 0.4005E-03, 0.1026E-02, + 0.2453E-02, 0.1026E-02, 0.1618E-01, 0.2100E-01, + 0.0000E+00, 0.0000E+00, 0.2694E-03, 0.1152E-02, + 0.2574E-02, 0.1152E-02, 0.1737E-01, 0.2265E-01, + 0.0000E+00, 0.0000E+00, 0.1276E-03, 0.1289E-02, + 0.2706E-02, 0.1289E-02, 0.1865E-01, 0.2446E-01/ DATA ((CALC(I,K,6 ),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.2581E-04, 0.1437E-02, + 0.2849E-02, 0.1437E-02, 0.2005E-01, 0.2643E-01, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,7 ),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1444E-02, 0.3563E-04, + 0.1515E-02, 0.3563E-04, 0.5775E-02, 0.5978E-02, + 0.0000E+00, 0.0000E+00, 0.1396E-02, 0.8289E-04, + 0.1562E-02, 0.8289E-04, 0.6166E-02, 0.6468E-02, + 0.0000E+00, 0.0000E+00, 0.1344E-02, 0.1337E-03, + 0.1612E-02, 0.1337E-03, 0.6585E-02, 0.6960E-02, + 0.0000E+00, 0.0000E+00, 0.1288E-02, 0.1885E-03, + 0.1665E-02, 0.1885E-03, 0.7032E-02, 0.7462E-02, + 0.0000E+00, 0.0000E+00, 0.1227E-02, 0.2475E-03, + 0.1722E-02, 0.2475E-03, 0.7513E-02, 0.7978E-02, + 0.0000E+00, 0.0000E+00, 0.1162E-02, 0.3111E-03, + 0.1784E-02, 0.3111E-03, 0.8029E-02, 0.8515E-02, + 0.0000E+00, 0.0000E+00, 0.1090E-02, 0.3797E-03, + 0.1850E-02, 0.3797E-03, 0.8584E-02, 0.9080E-02, + 0.0000E+00, 0.0000E+00, 0.1014E-02, 0.4538E-03, + 0.1921E-02, 0.4538E-03, 0.9182E-02, 0.9681E-02, + 0.0000E+00, 0.0000E+00, 0.9302E-03, 0.5338E-03, + 0.1998E-02, 0.5338E-03, 0.9828E-02, 0.1032E-01/ DATA ((CALC(I,K,7 ),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, 0.8400E-03, 0.6202E-03, + 0.2080E-02, 0.6202E-03, 0.1053E-01, 0.1102E-01, + 0.0000E+00, 0.0000E+00, 0.7424E-03, 0.7136E-03, + 0.2170E-02, 0.7136E-03, 0.1128E-01, 0.1177E-01, + 0.0000E+00, 0.0000E+00, 0.6368E-03, 0.8145E-03, + 0.2266E-02, 0.8145E-03, 0.1210E-01, 0.1260E-01, + 0.0000E+00, 0.0000E+00, 0.5226E-03, 0.9237E-03, + 0.2370E-02, 0.9237E-03, 0.1298E-01, 0.1349E-01, + 0.0000E+00, 0.0000E+00, 0.3990E-03, 0.1042E-02, + 0.2483E-02, 0.1042E-02, 0.1394E-01, 0.1448E-01, + 0.0000E+00, 0.0000E+00, 0.2653E-03, 0.1170E-02, + 0.2605E-02, 0.1170E-02, 0.1498E-01, 0.1555E-01, + 0.0000E+00, 0.0000E+00, 0.1206E-03, 0.1308E-02, + 0.2736E-02, 0.1308E-02, 0.1611E-01, 0.1673E-01, + 0.0000E+00, 0.0000E+00, -0.3584E-04, 0.1457E-02, + 0.2879E-02, 0.1457E-02, 0.1734E-01, 0.1803E-01, + 0.0000E+00, 0.0000E+00, -0.2051E-03, 0.1619E-02, + 0.3034E-02, 0.1619E-02, 0.1867E-01, 0.1945E-01/ DATA ((CALC(I,K,7 ),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.3882E-03, 0.1795E-02, + 0.3201E-02, 0.1795E-02, 0.2012E-01, 0.2101E-01, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,8 ),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1383E-02, 0.1197E-03, + 0.1623E-02, 0.1197E-03, 0.5533E-02, 0.5504E-02, + 0.0000E+00, 0.0000E+00, 0.1328E-02, 0.1726E-03, + 0.1674E-02, 0.1726E-03, 0.5934E-02, 0.5833E-02, + 0.0000E+00, 0.0000E+00, 0.1269E-02, 0.2296E-03, + 0.1728E-02, 0.2296E-03, 0.6363E-02, 0.6164E-02, + 0.0000E+00, 0.0000E+00, 0.1205E-02, 0.2912E-03, + 0.1787E-02, 0.2912E-03, 0.6823E-02, 0.6502E-02, + 0.0000E+00, 0.0000E+00, 0.1135E-02, 0.3576E-03, + 0.1850E-02, 0.3576E-03, 0.7318E-02, 0.6855E-02, + 0.0000E+00, 0.0000E+00, 0.1059E-02, 0.4293E-03, + 0.1918E-02, 0.4293E-03, 0.7851E-02, 0.7229E-02, + 0.0000E+00, 0.0000E+00, 0.9777E-03, 0.5069E-03, + 0.1991E-02, 0.5069E-03, 0.8426E-02, 0.7628E-02, + 0.0000E+00, 0.0000E+00, 0.8891E-03, 0.5907E-03, + 0.2070E-02, 0.5907E-03, 0.9048E-02, 0.8061E-02, + 0.0000E+00, 0.0000E+00, 0.7933E-03, 0.6813E-03, + 0.2156E-02, 0.6813E-03, 0.9720E-02, 0.8533E-02/ DATA ((CALC(I,K,8 ),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, 0.6895E-03, 0.7794E-03, + 0.2248E-02, 0.7794E-03, 0.1045E-01, 0.9051E-02, + 0.0000E+00, 0.0000E+00, 0.5772E-03, 0.8855E-03, + 0.2348E-02, 0.8855E-03, 0.1124E-01, 0.9622E-02, + 0.0000E+00, 0.0000E+00, 0.4557E-03, 0.1000E-02, + 0.2456E-02, 0.1000E-02, 0.1209E-01, 0.1025E-01, + 0.0000E+00, 0.0000E+00, 0.3242E-03, 0.1125E-02, + 0.2573E-02, 0.1125E-02, 0.1302E-01, 0.1095E-01, + 0.0000E+00, 0.0000E+00, 0.1819E-03, 0.1259E-02, + 0.2700E-02, 0.1259E-02, 0.1403E-01, 0.1172E-01, + 0.0000E+00, 0.0000E+00, 0.2794E-04, 0.1405E-02, + 0.2838E-02, 0.1405E-02, 0.1513E-01, 0.1257E-01, + 0.0000E+00, 0.0000E+00, -0.1386E-03, 0.1563E-02, + 0.2987E-02, 0.1563E-02, 0.1632E-01, 0.1351E-01, + 0.0000E+00, 0.0000E+00, -0.3188E-03, 0.1734E-02, + 0.3149E-02, 0.1734E-02, 0.1762E-01, 0.1455E-01, + 0.0000E+00, 0.0000E+00, -0.5136E-03, 0.1919E-02, + 0.3324E-02, 0.1919E-02, 0.1903E-01, 0.1569E-01/ DATA ((CALC(I,K,8 ),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.7244E-03, 0.2119E-02, + 0.3514E-02, 0.2119E-02, 0.2056E-01, 0.1695E-01, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,9 ),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1331E-02, 0.1925E-03, + 0.1716E-02, 0.1925E-03, 0.5326E-02, 0.5026E-02, + 0.0000E+00, 0.0000E+00, 0.1270E-02, 0.2501E-03, + 0.1771E-02, 0.2501E-03, 0.5741E-02, 0.5228E-02, + 0.0000E+00, 0.0000E+00, 0.1204E-02, 0.3123E-03, + 0.1829E-02, 0.3123E-03, 0.6186E-02, 0.5433E-02, + 0.0000E+00, 0.0000E+00, 0.1132E-02, 0.3796E-03, + 0.1891E-02, 0.3796E-03, 0.6664E-02, 0.5649E-02, + 0.0000E+00, 0.0000E+00, 0.1055E-02, 0.4523E-03, + 0.1959E-02, 0.4523E-03, 0.7180E-02, 0.5880E-02, + 0.0000E+00, 0.0000E+00, 0.9702E-03, 0.5311E-03, + 0.2032E-02, 0.5311E-03, 0.7737E-02, 0.6131E-02, + 0.0000E+00, 0.0000E+00, 0.8788E-03, 0.6164E-03, + 0.2112E-02, 0.6164E-03, 0.8339E-02, 0.6410E-02, + 0.0000E+00, 0.0000E+00, 0.7797E-03, 0.7087E-03, + 0.2197E-02, 0.7087E-03, 0.8992E-02, 0.6720E-02, + 0.0000E+00, 0.0000E+00, 0.6724E-03, 0.8088E-03, + 0.2290E-02, 0.8088E-03, 0.9699E-02, 0.7067E-02/ DATA ((CALC(I,K,9 ),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, 0.5562E-03, 0.9172E-03, + 0.2391E-02, 0.9172E-03, 0.1047E-01, 0.7458E-02, + 0.0000E+00, 0.0000E+00, 0.4303E-03, 0.1035E-02, + 0.2500E-02, 0.1035E-02, 0.1130E-01, 0.7896E-02, + 0.0000E+00, 0.0000E+00, 0.2941E-03, 0.1162E-02, + 0.2618E-02, 0.1162E-02, 0.1220E-01, 0.8388E-02, + 0.0000E+00, 0.0000E+00, 0.1467E-03, 0.1300E-02, + 0.2747E-02, 0.1300E-02, 0.1319E-01, 0.8939E-02, + 0.0000E+00, 0.0000E+00, -0.1287E-04, 0.1449E-02, + 0.2886E-02, 0.1449E-02, 0.1426E-01, 0.9555E-02, + 0.0000E+00, 0.0000E+00, -0.1855E-03, 0.1611E-02, + 0.3037E-02, 0.1611E-02, 0.1542E-01, 0.1024E-01, + 0.0000E+00, 0.0000E+00, -0.3723E-03, 0.1787E-02, + 0.3202E-02, 0.1787E-02, 0.1669E-01, 0.1100E-01, + 0.0000E+00, 0.0000E+00, -0.5743E-03, 0.1977E-02, + 0.3380E-02, 0.1977E-02, 0.1806E-01, 0.1185E-01, + 0.0000E+00, 0.0000E+00, -0.7929E-03, 0.2183E-02, + 0.3574E-02, 0.2183E-02, 0.1956E-01, 0.1278E-01/ DATA ((CALC(I,K,9 ),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.1029E-02, 0.2407E-02, + 0.3785E-02, 0.2407E-02, 0.2119E-01, 0.1381E-01, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,10),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1290E-02, 0.2527E-03, + 0.1796E-02, 0.2527E-03, 0.5161E-02, 0.4547E-02, + 0.0000E+00, 0.0000E+00, 0.1224E-02, 0.3140E-03, + 0.1852E-02, 0.3140E-03, 0.5591E-02, 0.4649E-02, + 0.0000E+00, 0.0000E+00, 0.1151E-02, 0.3804E-03, + 0.1912E-02, 0.3804E-03, 0.6053E-02, 0.4759E-02, + 0.0000E+00, 0.0000E+00, 0.1073E-02, 0.4525E-03, + 0.1978E-02, 0.4525E-03, 0.6551E-02, 0.4881E-02, + 0.0000E+00, 0.0000E+00, 0.9879E-03, 0.5306E-03, + 0.2049E-02, 0.5306E-03, 0.7090E-02, 0.5021E-02, + 0.0000E+00, 0.0000E+00, 0.8954E-03, 0.6154E-03, + 0.2126E-02, 0.6154E-03, 0.7673E-02, 0.5183E-02, + 0.0000E+00, 0.0000E+00, 0.7952E-03, 0.7074E-03, + 0.2210E-02, 0.7074E-03, 0.8304E-02, 0.5372E-02, + 0.0000E+00, 0.0000E+00, 0.6865E-03, 0.8072E-03, + 0.2301E-02, 0.8072E-03, 0.8990E-02, 0.5593E-02, + 0.0000E+00, 0.0000E+00, 0.5686E-03, 0.9156E-03, + 0.2400E-02, 0.9156E-03, 0.9734E-02, 0.5850E-02/ DATA ((CALC(I,K,10),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, 0.4410E-03, 0.1033E-02, + 0.2507E-02, 0.1033E-02, 0.1054E-01, 0.6146E-02, + 0.0000E+00, 0.0000E+00, 0.3027E-03, 0.1161E-02, + 0.2624E-02, 0.1161E-02, 0.1142E-01, 0.6487E-02, + 0.0000E+00, 0.0000E+00, 0.1529E-03, 0.1299E-02, + 0.2751E-02, 0.1299E-02, 0.1238E-01, 0.6877E-02, + 0.0000E+00, 0.0000E+00, -0.9250E-05, 0.1449E-02, + 0.2890E-02, 0.1449E-02, 0.1342E-01, 0.7320E-02, + 0.0000E+00, 0.0000E+00, -0.1847E-03, 0.1612E-02, + 0.3040E-02, 0.1612E-02, 0.1456E-01, 0.7820E-02, + 0.0000E+00, 0.0000E+00, -0.3747E-03, 0.1789E-02, + 0.3204E-02, 0.1789E-02, 0.1579E-01, 0.8381E-02, + 0.0000E+00, 0.0000E+00, -0.5801E-03, 0.1981E-02, + 0.3382E-02, 0.1981E-02, 0.1713E-01, 0.9008E-02, + 0.0000E+00, 0.0000E+00, -0.8024E-03, 0.2189E-02, + 0.3575E-02, 0.2189E-02, 0.1859E-01, 0.9705E-02, + 0.0000E+00, 0.0000E+00, -0.1043E-02, 0.2414E-02, + 0.3786E-02, 0.2414E-02, 0.2018E-01, 0.1048E-01/ DATA ((CALC(I,K,10),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.1303E-02, 0.2659E-02, + 0.4015E-02, 0.2659E-02, 0.2191E-01, 0.1133E-01, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,11),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1259E-02, 0.2998E-03, + 0.1859E-02, 0.2998E-03, 0.5037E-02, 0.4072E-02, + 0.0000E+00, 0.0000E+00, 0.1188E-02, 0.3641E-03, + 0.1916E-02, 0.3641E-03, 0.5481E-02, 0.4098E-02, + 0.0000E+00, 0.0000E+00, 0.1110E-02, 0.4339E-03, + 0.1978E-02, 0.4339E-03, 0.5959E-02, 0.4135E-02, + 0.0000E+00, 0.0000E+00, 0.1026E-02, 0.5100E-03, + 0.2046E-02, 0.5100E-03, 0.6476E-02, 0.4188E-02, + 0.0000E+00, 0.0000E+00, 0.9340E-03, 0.5927E-03, + 0.2119E-02, 0.5927E-03, 0.7036E-02, 0.4261E-02, + 0.0000E+00, 0.0000E+00, 0.8342E-03, 0.6827E-03, + 0.2200E-02, 0.6827E-03, 0.7645E-02, 0.4358E-02, + 0.0000E+00, 0.0000E+00, 0.7257E-03, 0.7807E-03, + 0.2287E-02, 0.7807E-03, 0.8305E-02, 0.4483E-02, + 0.0000E+00, 0.0000E+00, 0.6081E-03, 0.8871E-03, + 0.2382E-02, 0.8871E-03, 0.9024E-02, 0.4638E-02, + 0.0000E+00, 0.0000E+00, 0.4806E-03, 0.1003E-02, + 0.2486E-02, 0.1003E-02, 0.9807E-02, 0.4829E-02/ DATA ((CALC(I,K,11),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, 0.3422E-03, 0.1129E-02, + 0.2600E-02, 0.1129E-02, 0.1066E-01, 0.5057E-02, + 0.0000E+00, 0.0000E+00, 0.1923E-03, 0.1266E-02, + 0.2724E-02, 0.1266E-02, 0.1159E-01, 0.5327E-02, + 0.0000E+00, 0.0000E+00, 0.2992E-04, 0.1414E-02, + 0.2858E-02, 0.1414E-02, 0.1260E-01, 0.5641E-02, + 0.0000E+00, 0.0000E+00, -0.1460E-03, 0.1575E-02, + 0.3005E-02, 0.1575E-02, 0.1370E-01, 0.6002E-02, + 0.0000E+00, 0.0000E+00, -0.3363E-03, 0.1751E-02, + 0.3166E-02, 0.1751E-02, 0.1489E-01, 0.6414E-02, + 0.0000E+00, 0.0000E+00, -0.5424E-03, 0.1941E-02, + 0.3340E-02, 0.1941E-02, 0.1620E-01, 0.6879E-02, + 0.0000E+00, 0.0000E+00, -0.7653E-03, 0.2148E-02, + 0.3530E-02, 0.2148E-02, 0.1761E-01, 0.7397E-02, + 0.0000E+00, 0.0000E+00, -0.1006E-02, 0.2372E-02, + 0.3738E-02, 0.2372E-02, 0.1916E-01, 0.7983E-02, + 0.0000E+00, 0.0000E+00, -0.1267E-02, 0.2615E-02, + 0.3963E-02, 0.2615E-02, 0.2084E-01, 0.8632E-02/ DATA ((CALC(I,K,11),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.1549E-02, 0.2879E-02, + 0.4209E-02, 0.2879E-02, 0.2266E-01, 0.9347E-02, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,12),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1238E-02, 0.3350E-03, + 0.1908E-02, 0.3350E-03, 0.4950E-02, 0.3607E-02, + 0.0000E+00, 0.0000E+00, 0.1162E-02, 0.4016E-03, + 0.1965E-02, 0.4016E-03, 0.5405E-02, 0.3576E-02, + 0.0000E+00, 0.0000E+00, 0.1079E-02, 0.4744E-03, + 0.2028E-02, 0.4744E-03, 0.5897E-02, 0.3560E-02, + 0.0000E+00, 0.0000E+00, 0.9892E-03, 0.5539E-03, + 0.2097E-02, 0.5539E-03, 0.6431E-02, 0.3563E-02, + 0.0000E+00, 0.0000E+00, 0.8911E-03, 0.6407E-03, + 0.2173E-02, 0.6407E-03, 0.7012E-02, 0.3588E-02, + 0.0000E+00, 0.0000E+00, 0.7844E-03, 0.7355E-03, + 0.2255E-02, 0.7355E-03, 0.7645E-02, 0.3638E-02, + 0.0000E+00, 0.0000E+00, 0.6683E-03, 0.8388E-03, + 0.2346E-02, 0.8388E-03, 0.8335E-02, 0.3717E-02, + 0.0000E+00, 0.0000E+00, 0.5422E-03, 0.9514E-03, + 0.2445E-02, 0.9514E-03, 0.9087E-02, 0.3826E-02, + 0.0000E+00, 0.0000E+00, 0.4053E-03, 0.1074E-02, + 0.2554E-02, 0.1074E-02, 0.9907E-02, 0.3969E-02/ DATA ((CALC(I,K,12),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, 0.2568E-03, 0.1208E-02, + 0.2672E-02, 0.1208E-02, 0.1080E-01, 0.4147E-02, + 0.0000E+00, 0.0000E+00, 0.9573E-04, 0.1353E-02, + 0.2802E-02, 0.1353E-02, 0.1178E-01, 0.4364E-02, + 0.0000E+00, 0.0000E+00, -0.7882E-04, 0.1511E-02, + 0.2944E-02, 0.1511E-02, 0.1284E-01, 0.4621E-02, + 0.0000E+00, 0.0000E+00, -0.2679E-03, 0.1683E-02, + 0.3099E-02, 0.1683E-02, 0.1400E-01, 0.4920E-02, + 0.0000E+00, 0.0000E+00, -0.4726E-03, 0.1871E-02, + 0.3269E-02, 0.1871E-02, 0.1526E-01, 0.5264E-02, + 0.0000E+00, 0.0000E+00, -0.6941E-03, 0.2074E-02, + 0.3454E-02, 0.2074E-02, 0.1663E-01, 0.5654E-02, + 0.0000E+00, 0.0000E+00, -0.9338E-03, 0.2294E-02, + 0.3656E-02, 0.2294E-02, 0.1813E-01, 0.6094E-02, + 0.0000E+00, 0.0000E+00, -0.1193E-02, 0.2534E-02, + 0.3876E-02, 0.2534E-02, 0.1976E-01, 0.6586E-02, + 0.0000E+00, 0.0000E+00, -0.1473E-02, 0.2794E-02, + 0.4116E-02, 0.2794E-02, 0.2153E-01, 0.7132E-02/ DATA ((CALC(I,K,12),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.1776E-02, 0.3076E-02, + 0.4377E-02, 0.3076E-02, 0.2345E-01, 0.7735E-02, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,13),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1222E-02, 0.3601E-03, + 0.1942E-02, 0.3601E-03, 0.4890E-02, 0.3159E-02, + 0.0000E+00, 0.0000E+00, 0.1143E-02, 0.4288E-03, + 0.2000E-02, 0.4288E-03, 0.5354E-02, 0.3087E-02, + 0.0000E+00, 0.0000E+00, 0.1056E-02, 0.5042E-03, + 0.2064E-02, 0.5042E-03, 0.5859E-02, 0.3033E-02, + 0.0000E+00, 0.0000E+00, 0.9603E-03, 0.5869E-03, + 0.2134E-02, 0.5869E-03, 0.6411E-02, 0.3001E-02, + 0.0000E+00, 0.0000E+00, 0.8562E-03, 0.6777E-03, + 0.2211E-02, 0.6777E-03, 0.7013E-02, 0.2993E-02, + 0.0000E+00, 0.0000E+00, 0.7427E-03, 0.7770E-03, + 0.2297E-02, 0.7770E-03, 0.7672E-02, 0.3011E-02, + 0.0000E+00, 0.0000E+00, 0.6191E-03, 0.8857E-03, + 0.2391E-02, 0.8857E-03, 0.8391E-02, 0.3058E-02, + 0.0000E+00, 0.0000E+00, 0.4847E-03, 0.1004E-02, + 0.2494E-02, 0.1004E-02, 0.9178E-02, 0.3135E-02, + 0.0000E+00, 0.0000E+00, 0.3385E-03, 0.1134E-02, + 0.2607E-02, 0.1134E-02, 0.1004E-01, 0.3244E-02/ DATA ((CALC(I,K,13),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, 0.1798E-03, 0.1275E-02, + 0.2731E-02, 0.1275E-02, 0.1098E-01, 0.3387E-02, + 0.0000E+00, 0.0000E+00, 0.7585E-05, 0.1430E-02, + 0.2867E-02, 0.1430E-02, 0.1200E-01, 0.3564E-02, + 0.0000E+00, 0.0000E+00, -0.1791E-03, 0.1598E-02, + 0.3016E-02, 0.1598E-02, 0.1312E-01, 0.3777E-02, + 0.0000E+00, 0.0000E+00, -0.3814E-03, 0.1780E-02, + 0.3180E-02, 0.1780E-02, 0.1434E-01, 0.4029E-02, + 0.0000E+00, 0.0000E+00, -0.6005E-03, 0.1979E-02, + 0.3358E-02, 0.1979E-02, 0.1567E-01, 0.4319E-02, + 0.0000E+00, 0.0000E+00, -0.8375E-03, 0.2195E-02, + 0.3554E-02, 0.2195E-02, 0.1712E-01, 0.4650E-02, + 0.0000E+00, 0.0000E+00, -0.1094E-02, 0.2430E-02, + 0.3767E-02, 0.2430E-02, 0.1869E-01, 0.5023E-02, + 0.0000E+00, 0.0000E+00, -0.1371E-02, 0.2685E-02, + 0.4000E-02, 0.2685E-02, 0.2041E-01, 0.5440E-02, + 0.0000E+00, 0.0000E+00, -0.1671E-02, 0.2962E-02, + 0.4254E-02, 0.2962E-02, 0.2227E-01, 0.5903E-02/ DATA ((CALC(I,K,13),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.1994E-02, 0.3262E-02, + 0.4531E-02, 0.3262E-02, 0.2430E-01, 0.6414E-02, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,14),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1213E-02, 0.3767E-03, + 0.1967E-02, 0.3767E-03, 0.4853E-02, 0.2733E-02, + 0.0000E+00, 0.0000E+00, 0.1130E-02, 0.4476E-03, + 0.2025E-02, 0.4476E-03, 0.5327E-02, 0.2633E-02, + 0.0000E+00, 0.0000E+00, 0.1038E-02, 0.5258E-03, + 0.2090E-02, 0.5258E-03, 0.5847E-02, 0.2554E-02, + 0.0000E+00, 0.0000E+00, 0.9374E-03, 0.6121E-03, + 0.2162E-02, 0.6121E-03, 0.6418E-02, 0.2500E-02, + 0.0000E+00, 0.0000E+00, 0.8271E-03, 0.7071E-03, + 0.2241E-02, 0.7071E-03, 0.7044E-02, 0.2471E-02, + 0.0000E+00, 0.0000E+00, 0.7066E-03, 0.8115E-03, + 0.2330E-02, 0.8115E-03, 0.7731E-02, 0.2468E-02, + 0.0000E+00, 0.0000E+00, 0.5751E-03, 0.9260E-03, + 0.2427E-02, 0.9260E-03, 0.8484E-02, 0.2494E-02, + 0.0000E+00, 0.0000E+00, 0.4317E-03, 0.1051E-02, + 0.2535E-02, 0.1051E-02, 0.9309E-02, 0.2550E-02, + 0.0000E+00, 0.0000E+00, 0.2757E-03, 0.1189E-02, + 0.2653E-02, 0.1189E-02, 0.1021E-01, 0.2635E-02/ DATA ((CALC(I,K,14),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, 0.1062E-03, 0.1340E-02, + 0.2783E-02, 0.1340E-02, 0.1120E-01, 0.2752E-02, + 0.0000E+00, 0.0000E+00, -0.7778E-04, 0.1503E-02, + 0.2927E-02, 0.1503E-02, 0.1228E-01, 0.2900E-02, + 0.0000E+00, 0.0000E+00, -0.2773E-03, 0.1681E-02, + 0.3084E-02, 0.1681E-02, 0.1346E-01, 0.3080E-02, + 0.0000E+00, 0.0000E+00, -0.4936E-03, 0.1875E-02, + 0.3256E-02, 0.1875E-02, 0.1475E-01, 0.3294E-02, + 0.0000E+00, 0.0000E+00, -0.7277E-03, 0.2086E-02, + 0.3444E-02, 0.2086E-02, 0.1615E-01, 0.3541E-02, + 0.0000E+00, 0.0000E+00, -0.9809E-03, 0.2316E-02, + 0.3650E-02, 0.2316E-02, 0.1768E-01, 0.3824E-02, + 0.0000E+00, 0.0000E+00, -0.1255E-02, 0.2565E-02, + 0.3876E-02, 0.2565E-02, 0.1934E-01, 0.4142E-02, + 0.0000E+00, 0.0000E+00, -0.1551E-02, 0.2836E-02, + 0.4122E-02, 0.2836E-02, 0.2115E-01, 0.4497E-02, + 0.0000E+00, 0.0000E+00, -0.1870E-02, 0.3130E-02, + 0.4391E-02, 0.3130E-02, 0.2312E-01, 0.4891E-02/ DATA ((CALC(I,K,14),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.2215E-02, 0.3449E-02, + 0.4684E-02, 0.3449E-02, 0.2525E-01, 0.5325E-02, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,15),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1209E-02, 0.3877E-03, + 0.1984E-02, 0.3877E-03, 0.4835E-02, 0.2333E-02, + 0.0000E+00, 0.0000E+00, 0.1121E-02, 0.4613E-03, + 0.2044E-02, 0.4613E-03, 0.5323E-02, 0.2216E-02, + 0.0000E+00, 0.0000E+00, 0.1025E-02, 0.5430E-03, + 0.2111E-02, 0.5430E-03, 0.5862E-02, 0.2123E-02, + 0.0000E+00, 0.0000E+00, 0.9181E-03, 0.6336E-03, + 0.2185E-02, 0.6336E-03, 0.6456E-02, 0.2056E-02, + 0.0000E+00, 0.0000E+00, 0.8010E-03, 0.7337E-03, + 0.2268E-02, 0.7337E-03, 0.7112E-02, 0.2015E-02, + 0.0000E+00, 0.0000E+00, 0.6726E-03, 0.8440E-03, + 0.2361E-02, 0.8440E-03, 0.7834E-02, 0.2001E-02, + 0.0000E+00, 0.0000E+00, 0.5322E-03, 0.9654E-03, + 0.2463E-02, 0.9654E-03, 0.8627E-02, 0.2015E-02, + 0.0000E+00, 0.0000E+00, 0.3789E-03, 0.1098E-02, + 0.2576E-02, 0.1098E-02, 0.9499E-02, 0.2057E-02, + 0.0000E+00, 0.0000E+00, 0.2120E-03, 0.1244E-02, + 0.2701E-02, 0.1244E-02, 0.1045E-01, 0.2126E-02/ DATA ((CALC(I,K,15),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, 0.3050E-04, 0.1404E-02, + 0.2838E-02, 0.1404E-02, 0.1150E-01, 0.2224E-02, + 0.0000E+00, 0.0000E+00, -0.1666E-03, 0.1578E-02, + 0.2990E-02, 0.1578E-02, 0.1264E-01, 0.2350E-02, + 0.0000E+00, 0.0000E+00, -0.3804E-03, 0.1768E-02, + 0.3156E-02, 0.1768E-02, 0.1389E-01, 0.2505E-02, + 0.0000E+00, 0.0000E+00, -0.6120E-03, 0.1974E-02, + 0.3338E-02, 0.1974E-02, 0.1526E-01, 0.2688E-02, + 0.0000E+00, 0.0000E+00, -0.8626E-03, 0.2199E-02, + 0.3537E-02, 0.2199E-02, 0.1674E-01, 0.2901E-02, + 0.0000E+00, 0.0000E+00, -0.1134E-02, 0.2444E-02, + 0.3756E-02, 0.2444E-02, 0.1836E-01, 0.3143E-02, + 0.0000E+00, 0.0000E+00, -0.1426E-02, 0.2710E-02, + 0.3995E-02, 0.2710E-02, 0.2012E-01, 0.3415E-02, + 0.0000E+00, 0.0000E+00, -0.1743E-02, 0.2999E-02, + 0.4255E-02, 0.2999E-02, 0.2203E-01, 0.3719E-02, + 0.0000E+00, 0.0000E+00, -0.2084E-02, 0.3311E-02, + 0.4540E-02, 0.3311E-02, 0.2411E-01, 0.4054E-02/ DATA ((CALC(I,K,15),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.2452E-02, 0.3651E-02, + 0.4851E-02, 0.3651E-02, 0.2637E-01, 0.4422E-02, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,16),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1210E-02, 0.3954E-03, + 0.2001E-02, 0.3954E-03, 0.4839E-02, 0.1962E-02, + 0.0000E+00, 0.0000E+00, 0.1118E-02, 0.4727E-03, + 0.2063E-02, 0.4727E-03, 0.5348E-02, 0.1838E-02, + 0.0000E+00, 0.0000E+00, 0.1015E-02, 0.5591E-03, + 0.2133E-02, 0.5591E-03, 0.5913E-02, 0.1740E-02, + 0.0000E+00, 0.0000E+00, 0.9014E-03, 0.6552E-03, + 0.2212E-02, 0.6552E-03, 0.6540E-02, 0.1668E-02, + 0.0000E+00, 0.0000E+00, 0.7761E-03, 0.7618E-03, + 0.2300E-02, 0.7618E-03, 0.7234E-02, 0.1623E-02, + 0.0000E+00, 0.0000E+00, 0.6384E-03, 0.8795E-03, + 0.2397E-02, 0.8795E-03, 0.8001E-02, 0.1604E-02, + 0.0000E+00, 0.0000E+00, 0.4875E-03, 0.1009E-02, + 0.2506E-02, 0.1009E-02, 0.8845E-02, 0.1611E-02, + 0.0000E+00, 0.0000E+00, 0.3226E-03, 0.1152E-02, + 0.2626E-02, 0.1152E-02, 0.9773E-02, 0.1644E-02, + 0.0000E+00, 0.0000E+00, 0.1428E-03, 0.1308E-02, + 0.2759E-02, 0.1308E-02, 0.1079E-01, 0.1703E-02/ DATA ((CALC(I,K,16),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, -0.5268E-04, 0.1478E-02, + 0.2905E-02, 0.1478E-02, 0.1191E-01, 0.1788E-02, + 0.0000E+00, 0.0000E+00, -0.2649E-03, 0.1665E-02, + 0.3066E-02, 0.1665E-02, 0.1313E-01, 0.1897E-02, + 0.0000E+00, 0.0000E+00, -0.4950E-03, 0.1868E-02, + 0.3242E-02, 0.1868E-02, 0.1446E-01, 0.2031E-02, + 0.0000E+00, 0.0000E+00, -0.7441E-03, 0.2090E-02, + 0.3436E-02, 0.2090E-02, 0.1591E-01, 0.2190E-02, + 0.0000E+00, 0.0000E+00, -0.1013E-02, 0.2330E-02, + 0.3648E-02, 0.2330E-02, 0.1749E-01, 0.2373E-02, + 0.0000E+00, 0.0000E+00, -0.1304E-02, 0.2592E-02, + 0.3881E-02, 0.2592E-02, 0.1922E-01, 0.2582E-02, + 0.0000E+00, 0.0000E+00, -0.1618E-02, 0.2876E-02, + 0.4134E-02, 0.2876E-02, 0.2109E-01, 0.2815E-02, + 0.0000E+00, 0.0000E+00, -0.1957E-02, 0.3183E-02, + 0.4412E-02, 0.3183E-02, 0.2312E-01, 0.3073E-02, + 0.0000E+00, 0.0000E+00, -0.2322E-02, 0.3517E-02, + 0.4714E-02, 0.3517E-02, 0.2533E-01, 0.3358E-02/ DATA ((CALC(I,K,16),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.2715E-02, 0.3878E-02, + 0.5044E-02, 0.3878E-02, 0.2772E-01, 0.3670E-02, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,17),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1218E-02, 0.4021E-03, + 0.2022E-02, 0.4021E-03, 0.4871E-02, 0.1625E-02, + 0.0000E+00, 0.0000E+00, 0.1119E-02, 0.4848E-03, + 0.2089E-02, 0.4848E-03, 0.5411E-02, 0.1501E-02, + 0.0000E+00, 0.0000E+00, 0.1009E-02, 0.5774E-03, + 0.2164E-02, 0.5774E-03, 0.6014E-02, 0.1404E-02, + 0.0000E+00, 0.0000E+00, 0.8866E-03, 0.6809E-03, + 0.2248E-02, 0.6809E-03, 0.6686E-02, 0.1333E-02, + 0.0000E+00, 0.0000E+00, 0.7510E-03, 0.7958E-03, + 0.2343E-02, 0.7958E-03, 0.7431E-02, 0.1288E-02, + 0.0000E+00, 0.0000E+00, 0.6017E-03, 0.9229E-03, + 0.2448E-02, 0.9229E-03, 0.8255E-02, 0.1269E-02, + 0.0000E+00, 0.0000E+00, 0.4378E-03, 0.1063E-02, + 0.2564E-02, 0.1063E-02, 0.9164E-02, 0.1274E-02, + 0.0000E+00, 0.0000E+00, 0.2588E-03, 0.1217E-02, + 0.2693E-02, 0.1217E-02, 0.1016E-01, 0.1304E-02, + 0.0000E+00, 0.0000E+00, 0.6363E-04, 0.1385E-02, + 0.2835E-02, 0.1385E-02, 0.1126E-01, 0.1356E-02/ DATA ((CALC(I,K,17),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, -0.1485E-03, 0.1569E-02, + 0.2992E-02, 0.1569E-02, 0.1246E-01, 0.1430E-02, + 0.0000E+00, 0.0000E+00, -0.3786E-03, 0.1771E-02, + 0.3164E-02, 0.1771E-02, 0.1377E-01, 0.1526E-02, + 0.0000E+00, 0.0000E+00, -0.6278E-03, 0.1989E-02, + 0.3353E-02, 0.1989E-02, 0.1520E-01, 0.1643E-02, + 0.0000E+00, 0.0000E+00, -0.8971E-03, 0.2227E-02, + 0.3560E-02, 0.2227E-02, 0.1675E-01, 0.1781E-02, + 0.0000E+00, 0.0000E+00, -0.1188E-02, 0.2486E-02, + 0.3787E-02, 0.2486E-02, 0.1845E-01, 0.1939E-02, + 0.0000E+00, 0.0000E+00, -0.1502E-02, 0.2767E-02, + 0.4034E-02, 0.2767E-02, 0.2029E-01, 0.2118E-02, + 0.0000E+00, 0.0000E+00, -0.1839E-02, 0.3071E-02, + 0.4305E-02, 0.3071E-02, 0.2229E-01, 0.2317E-02, + 0.0000E+00, 0.0000E+00, -0.2203E-02, 0.3400E-02, + 0.4600E-02, 0.3400E-02, 0.2445E-01, 0.2536E-02, + 0.0000E+00, 0.0000E+00, -0.2595E-02, 0.3757E-02, + 0.4922E-02, 0.3757E-02, 0.2680E-01, 0.2777E-02/ DATA ((CALC(I,K,17),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.3016E-02, 0.4143E-02, + 0.5273E-02, 0.4143E-02, 0.2935E-01, 0.3040E-02, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,18),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1235E-02, 0.4100E-03, + 0.2055E-02, 0.4100E-03, 0.4941E-02, 0.1323E-02, + 0.0000E+00, 0.0000E+00, 0.1128E-02, 0.5000E-03, + 0.2128E-02, 0.5000E-03, 0.5526E-02, 0.1205E-02, + 0.0000E+00, 0.0000E+00, 0.1008E-02, 0.6013E-03, + 0.2210E-02, 0.6013E-03, 0.6182E-02, 0.1114E-02, + 0.0000E+00, 0.0000E+00, 0.8735E-03, 0.7144E-03, + 0.2302E-02, 0.7144E-03, 0.6913E-02, 0.1048E-02, + 0.0000E+00, 0.0000E+00, 0.7246E-03, 0.8401E-03, + 0.2405E-02, 0.8401E-03, 0.7724E-02, 0.1008E-02, + 0.0000E+00, 0.0000E+00, 0.5606E-03, 0.9790E-03, + 0.2519E-02, 0.9790E-03, 0.8621E-02, 0.9914E-03, + 0.0000E+00, 0.0000E+00, 0.3806E-03, 0.1132E-02, + 0.2645E-02, 0.1132E-02, 0.9609E-02, 0.9973E-03, + 0.0000E+00, 0.0000E+00, 0.1842E-03, 0.1300E-02, + 0.2784E-02, 0.1300E-02, 0.1069E-01, 0.1024E-02, + 0.0000E+00, 0.0000E+00, -0.2965E-04, 0.1483E-02, + 0.2937E-02, 0.1483E-02, 0.1188E-01, 0.1072E-02/ DATA ((CALC(I,K,18),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, -0.2617E-03, 0.1683E-02, + 0.3106E-02, 0.1683E-02, 0.1318E-01, 0.1138E-02, + 0.0000E+00, 0.0000E+00, -0.5129E-03, 0.1901E-02, + 0.3291E-02, 0.1901E-02, 0.1460E-01, 0.1223E-02, + 0.0000E+00, 0.0000E+00, -0.7843E-03, 0.2138E-02, + 0.3494E-02, 0.2138E-02, 0.1614E-01, 0.1326E-02, + 0.0000E+00, 0.0000E+00, -0.1077E-02, 0.2396E-02, + 0.3716E-02, 0.2396E-02, 0.1782E-01, 0.1445E-02, + 0.0000E+00, 0.0000E+00, -0.1392E-02, 0.2674E-02, + 0.3959E-02, 0.2674E-02, 0.1964E-01, 0.1581E-02, + 0.0000E+00, 0.0000E+00, -0.1732E-02, 0.2977E-02, + 0.4224E-02, 0.2977E-02, 0.2161E-01, 0.1734E-02, + 0.0000E+00, 0.0000E+00, -0.2096E-02, 0.3303E-02, + 0.4513E-02, 0.3303E-02, 0.2375E-01, 0.1902E-02, + 0.0000E+00, 0.0000E+00, -0.2488E-02, 0.3657E-02, + 0.4828E-02, 0.3657E-02, 0.2607E-01, 0.2087E-02, + 0.0000E+00, 0.0000E+00, -0.2909E-02, 0.4039E-02, + 0.5171E-02, 0.4039E-02, 0.2858E-01, 0.2290E-02/ DATA ((CALC(I,K,18),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.3361E-02, 0.4451E-02, + 0.5545E-02, 0.4451E-02, 0.3129E-01, 0.2510E-02, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,19),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1265E-02, 0.4210E-03, + 0.2108E-02, 0.4210E-03, 0.5061E-02, 0.1058E-02, + 0.0000E+00, 0.0000E+00, 0.1146E-02, 0.5211E-03, + 0.2189E-02, 0.5211E-03, 0.5708E-02, 0.9493E-03, + 0.0000E+00, 0.0000E+00, 0.1012E-02, 0.6336E-03, + 0.2279E-02, 0.6336E-03, 0.6433E-02, 0.8673E-03, + 0.0000E+00, 0.0000E+00, 0.8619E-03, 0.7591E-03, + 0.2380E-02, 0.7591E-03, 0.7241E-02, 0.8103E-03, + 0.0000E+00, 0.0000E+00, 0.6959E-03, 0.8982E-03, + 0.2492E-02, 0.8982E-03, 0.8136E-02, 0.7765E-03, + 0.0000E+00, 0.0000E+00, 0.5132E-03, 0.1052E-02, + 0.2617E-02, 0.1052E-02, 0.9123E-02, 0.7644E-03, + 0.0000E+00, 0.0000E+00, 0.3133E-03, 0.1220E-02, + 0.2754E-02, 0.1220E-02, 0.1021E-01, 0.7723E-03, + 0.0000E+00, 0.0000E+00, 0.9558E-04, 0.1404E-02, + 0.2905E-02, 0.1404E-02, 0.1139E-01, 0.7988E-03, + 0.0000E+00, 0.0000E+00, -0.1408E-03, 0.1605E-02, + 0.3071E-02, 0.1605E-02, 0.1269E-01, 0.8425E-03/ DATA ((CALC(I,K,19),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, -0.3964E-03, 0.1823E-02, + 0.3253E-02, 0.1823E-02, 0.1411E-01, 0.9023E-03, + 0.0000E+00, 0.0000E+00, -0.6723E-03, 0.2062E-02, + 0.3452E-02, 0.2062E-02, 0.1564E-01, 0.9773E-03, + 0.0000E+00, 0.0000E+00, -0.9694E-03, 0.2319E-02, + 0.3670E-02, 0.2319E-02, 0.1731E-01, 0.1067E-02, + 0.0000E+00, 0.0000E+00, -0.1289E-02, 0.2597E-02, + 0.3908E-02, 0.2597E-02, 0.1912E-01, 0.1169E-02, + 0.0000E+00, 0.0000E+00, -0.1632E-02, 0.2898E-02, + 0.4168E-02, 0.2898E-02, 0.2108E-01, 0.1286E-02, + 0.0000E+00, 0.0000E+00, -0.2000E-02, 0.3224E-02, + 0.4452E-02, 0.3224E-02, 0.2321E-01, 0.1415E-02, + 0.0000E+00, 0.0000E+00, -0.2394E-02, 0.3576E-02, + 0.4761E-02, 0.3576E-02, 0.2550E-01, 0.1556E-02, + 0.0000E+00, 0.0000E+00, -0.2817E-02, 0.3955E-02, + 0.5097E-02, 0.3955E-02, 0.2798E-01, 0.1712E-02, + 0.0000E+00, 0.0000E+00, -0.3269E-02, 0.4364E-02, + 0.5464E-02, 0.4364E-02, 0.3066E-01, 0.1879E-02/ DATA ((CALC(I,K,19),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.3753E-02, 0.4805E-02, + 0.5862E-02, 0.4805E-02, 0.3356E-01, 0.2062E-02, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,20),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1311E-02, 0.4367E-03, + 0.2185E-02, 0.4367E-03, 0.5245E-02, 0.8288E-03, + 0.0000E+00, 0.0000E+00, 0.1175E-02, 0.5500E-03, + 0.2275E-02, 0.5500E-03, 0.5974E-02, 0.7327E-03, + 0.0000E+00, 0.0000E+00, 0.1022E-02, 0.6769E-03, + 0.2376E-02, 0.6769E-03, 0.6787E-02, 0.6623E-03, + 0.0000E+00, 0.0000E+00, 0.8518E-03, 0.8178E-03, + 0.2487E-02, 0.8178E-03, 0.7689E-02, 0.6149E-03, + 0.0000E+00, 0.0000E+00, 0.6638E-03, 0.9733E-03, + 0.2610E-02, 0.9733E-03, 0.8685E-02, 0.5888E-03, + 0.0000E+00, 0.0000E+00, 0.4578E-03, 0.1144E-02, + 0.2746E-02, 0.1144E-02, 0.9779E-02, 0.5818E-03, + 0.0000E+00, 0.0000E+00, 0.2333E-03, 0.1331E-02, + 0.2895E-02, 0.1331E-02, 0.1098E-01, 0.5922E-03, + 0.0000E+00, 0.0000E+00, -0.1005E-04, 0.1534E-02, + 0.3059E-02, 0.1534E-02, 0.1228E-01, 0.6184E-03, + 0.0000E+00, 0.0000E+00, -0.2729E-03, 0.1755E-02, + 0.3238E-02, 0.1755E-02, 0.1370E-01, 0.6589E-03/ DATA ((CALC(I,K,20),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, -0.5560E-03, 0.1989E-02, + 0.3434E-02, 0.1989E-02, 0.1524E-01, 0.7126E-03, + 0.0000E+00, 0.0000E+00, -0.8602E-03, 0.2249E-02, + 0.3648E-02, 0.2249E-02, 0.1691E-01, 0.7783E-03, + 0.0000E+00, 0.0000E+00, -0.1186E-02, 0.2529E-02, + 0.3882E-02, 0.2529E-02, 0.1872E-01, 0.8554E-03, + 0.0000E+00, 0.0000E+00, -0.1535E-02, 0.2832E-02, + 0.4138E-02, 0.2832E-02, 0.2068E-01, 0.9431E-03, + 0.0000E+00, 0.0000E+00, -0.1909E-02, 0.3158E-02, + 0.4416E-02, 0.3158E-02, 0.2279E-01, 0.1041E-02, + 0.0000E+00, 0.0000E+00, -0.2308E-02, 0.3509E-02, + 0.4718E-02, 0.3509E-02, 0.2507E-01, 0.1149E-02, + 0.0000E+00, 0.0000E+00, -0.2734E-02, 0.3887E-02, + 0.5048E-02, 0.3887E-02, 0.2753E-01, 0.1267E-02, + 0.0000E+00, 0.0000E+00, -0.3189E-02, 0.4294E-02, + 0.5406E-02, 0.4294E-02, 0.3018E-01, 0.1394E-02, + 0.0000E+00, 0.0000E+00, -0.3675E-02, 0.4732E-02, + 0.5796E-02, 0.4732E-02, 0.3305E-01, 0.1533E-02/ DATA ((CALC(I,K,20),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.4194E-02, 0.5202E-02, + 0.6219E-02, 0.5202E-02, 0.3613E-01, 0.1682E-02, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,21),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1377E-02, 0.4583E-03, + 0.2293E-02, 0.4583E-03, 0.5506E-02, 0.6351E-03, + 0.0000E+00, 0.0000E+00, 0.1217E-02, 0.5884E-03, + 0.2394E-02, 0.5884E-03, 0.6336E-02, 0.5532E-03, + 0.0000E+00, 0.0000E+00, 0.1039E-02, 0.7330E-03, + 0.2505E-02, 0.7330E-03, 0.7257E-02, 0.4951E-03, + 0.0000E+00, 0.0000E+00, 0.8426E-03, 0.8925E-03, + 0.2628E-02, 0.8925E-03, 0.8272E-02, 0.4579E-03, + 0.0000E+00, 0.0000E+00, 0.6271E-03, 0.1068E-02, + 0.2762E-02, 0.1068E-02, 0.9385E-02, 0.4395E-03, + 0.0000E+00, 0.0000E+00, 0.3917E-03, 0.1258E-02, + 0.2909E-02, 0.1258E-02, 0.1060E-01, 0.4375E-03, + 0.0000E+00, 0.0000E+00, 0.1382E-03, 0.1466E-02, + 0.3071E-02, 0.1466E-02, 0.1192E-01, 0.4502E-03, + 0.0000E+00, 0.0000E+00, -0.1353E-03, 0.1690E-02, + 0.3247E-02, 0.1690E-02, 0.1336E-01, 0.4759E-03, + 0.0000E+00, 0.0000E+00, -0.4289E-03, 0.1933E-02, + 0.3439E-02, 0.1933E-02, 0.1492E-01, 0.5131E-03/ DATA ((CALC(I,K,21),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, -0.7432E-03, 0.2195E-02, + 0.3650E-02, 0.2195E-02, 0.1660E-01, 0.5606E-03, + 0.0000E+00, 0.0000E+00, -0.1079E-02, 0.2479E-02, + 0.3879E-02, 0.2479E-02, 0.1841E-01, 0.6176E-03, + 0.0000E+00, 0.0000E+00, -0.1437E-02, 0.2782E-02, + 0.4129E-02, 0.2782E-02, 0.2036E-01, 0.6831E-03, + 0.0000E+00, 0.0000E+00, -0.1818E-02, 0.3108E-02, + 0.4400E-02, 0.3108E-02, 0.2247E-01, 0.7567E-03, + 0.0000E+00, 0.0000E+00, -0.2224E-02, 0.3459E-02, + 0.4696E-02, 0.3459E-02, 0.2474E-01, 0.8379E-03, + 0.0000E+00, 0.0000E+00, -0.2655E-02, 0.3835E-02, + 0.5018E-02, 0.3835E-02, 0.2718E-01, 0.9261E-03, + 0.0000E+00, 0.0000E+00, -0.3114E-02, 0.4239E-02, + 0.5367E-02, 0.4239E-02, 0.2980E-01, 0.1021E-02, + 0.0000E+00, 0.0000E+00, -0.3603E-02, 0.4673E-02, + 0.5747E-02, 0.4673E-02, 0.3265E-01, 0.1132E-02, + 0.0000E+00, 0.0000E+00, -0.4122E-02, 0.5139E-02, + 0.6160E-02, 0.5139E-02, 0.3569E-01, 0.1240E-02/ DATA ((CALC(I,K,21),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.4675E-02, 0.5640E-02, + 0.6608E-02, 0.5640E-02, 0.3896E-01, 0.1357E-02, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,22),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1462E-02, 0.4873E-03, + 0.2437E-02, 0.4873E-03, 0.5848E-02, 0.4746E-03, + 0.0000E+00, 0.0000E+00, 0.1272E-02, 0.6379E-03, + 0.2548E-02, 0.6379E-03, 0.6801E-02, 0.4074E-03, + 0.0000E+00, 0.0000E+00, 0.1062E-02, 0.8038E-03, + 0.2670E-02, 0.8038E-03, 0.7849E-02, 0.3617E-03, + 0.0000E+00, 0.0000E+00, 0.8324E-03, 0.9852E-03, + 0.2803E-02, 0.9852E-03, 0.8994E-02, 0.3345E-03, + 0.0000E+00, 0.0000E+00, 0.5830E-03, 0.1182E-02, + 0.2948E-02, 0.1182E-02, 0.1024E-01, 0.3231E-03, + 0.0000E+00, 0.0000E+00, 0.3145E-03, 0.1396E-02, + 0.3106E-02, 0.1396E-02, 0.1159E-01, 0.3256E-03, + 0.0000E+00, 0.0000E+00, 0.2579E-04, 0.1627E-02, + 0.3279E-02, 0.1627E-02, 0.1305E-01, 0.3399E-03, + 0.0000E+00, 0.0000E+00, -0.2823E-03, 0.1870E-02, + 0.3467E-02, 0.1870E-02, 0.1463E-01, 0.3645E-03, + 0.0000E+00, 0.0000E+00, -0.6105E-03, 0.2145E-02, + 0.3672E-02, 0.2145E-02, 0.1632E-01, 0.3980E-03/ DATA ((CALC(I,K,22),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, -0.9591E-03, 0.2431E-02, + 0.3895E-02, 0.2431E-02, 0.1815E-01, 0.4394E-03, + 0.0000E+00, 0.0000E+00, -0.1329E-02, 0.2736E-02, + 0.4138E-02, 0.2736E-02, 0.2011E-01, 0.4877E-03, + 0.0000E+00, 0.0000E+00, -0.1721E-02, 0.3061E-02, + 0.4403E-02, 0.3061E-02, 0.2221E-01, 0.5423E-03, + 0.0000E+00, 0.0000E+00, -0.2135E-02, 0.3411E-02, + 0.4690E-02, 0.3411E-02, 0.2447E-01, 0.6027E-03, + 0.0000E+00, 0.0000E+00, -0.2574E-02, 0.3786E-02, + 0.5002E-02, 0.3786E-02, 0.2689E-01, 0.6685E-03, + 0.0000E+00, 0.0000E+00, -0.3038E-02, 0.4187E-02, + 0.5341E-02, 0.4187E-02, 0.2949E-01, 0.7397E-03, + 0.0000E+00, 0.0000E+00, -0.3530E-02, 0.4617E-02, + 0.5710E-02, 0.4617E-02, 0.3227E-01, 0.8160E-03, + 0.0000E+00, 0.0000E+00, -0.4051E-02, 0.5077E-02, + 0.6110E-02, 0.5077E-02, 0.3526E-01, 0.8979E-03, + 0.0000E+00, 0.0000E+00, -0.4603E-02, 0.5570E-02, + 0.6544E-02, 0.5570E-02, 0.3848E-01, 0.9852E-03/ DATA ((CALC(I,K,22),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.5189E-02, 0.6098E-02, + 0.7015E-02, 0.6098E-02, 0.4193E-01, 0.1079E-02, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,23),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1571E-02, 0.5237E-03, + 0.2618E-02, 0.5237E-03, 0.6284E-02, 0.3454E-03, + 0.0000E+00, 0.0000E+00, 0.1341E-02, 0.6990E-03, + 0.2739E-02, 0.6990E-03, 0.7380E-02, 0.2925E-03, + 0.0000E+00, 0.0000E+00, 0.1091E-02, 0.8897E-03, + 0.2870E-02, 0.8897E-03, 0.8573E-02, 0.2584E-03, + 0.0000E+00, 0.0000E+00, 0.8202E-03, 0.1096E-02, + 0.3012E-02, 0.1096E-02, 0.9863E-02, 0.2400E-03, + 0.0000E+00, 0.0000E+00, 0.5302E-03, 0.1318E-02, + 0.3166E-02, 0.1318E-02, 0.1125E-01, 0.2347E-03, + 0.0000E+00, 0.0000E+00, 0.2214E-03, 0.1557E-02, + 0.3334E-02, 0.1557E-02, 0.1275E-01, 0.2404E-03, + 0.0000E+00, 0.0000E+00, -0.1069E-03, 0.1728E-02, + 0.3516E-02, 0.1728E-02, 0.1435E-01, 0.2554E-03, + 0.0000E+00, 0.0000E+00, -0.4535E-03, 0.2033E-02, + 0.3715E-02, 0.2033E-02, 0.1607E-01, 0.2782E-03, + 0.0000E+00, 0.0000E+00, -0.8196E-03, 0.2349E-02, + 0.3930E-02, 0.2349E-02, 0.1791E-01, 0.3075E-03/ DATA ((CALC(I,K,23),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, -0.1205E-02, 0.2657E-02, + 0.4164E-02, 0.2657E-02, 0.1987E-01, 0.3425E-03, + 0.0000E+00, 0.0000E+00, -0.1610E-02, 0.2992E-02, + 0.4419E-02, 0.2992E-02, 0.2197E-01, 0.3825E-03, + 0.0000E+00, 0.0000E+00, -0.2036E-02, 0.3345E-02, + 0.4695E-02, 0.3345E-02, 0.2422E-01, 0.4267E-03, + 0.0000E+00, 0.0000E+00, -0.2484E-02, 0.3721E-02, + 0.4996E-02, 0.3721E-02, 0.2661E-01, 0.4750E-03, + 0.0000E+00, 0.0000E+00, -0.2955E-02, 0.4122E-02, + 0.5323E-02, 0.4122E-02, 0.2918E-01, 0.5270E-03, + 0.0000E+00, 0.0000E+00, -0.3451E-02, 0.4549E-02, + 0.5677E-02, 0.4549E-02, 0.3192E-01, 0.5827E-03, + 0.0000E+00, 0.0000E+00, -0.3973E-02, 0.5005E-02, + 0.6061E-02, 0.5005E-02, 0.3485E-01, 0.6419E-03, + 0.0000E+00, 0.0000E+00, -0.4524E-02, 0.5490E-02, + 0.6479E-02, 0.5490E-02, 0.3799E-01, 0.7056E-03, + 0.0000E+00, 0.0000E+00, -0.5105E-02, 0.6009E-02, + 0.6932E-02, 0.6009E-02, 0.4136E-01, 0.7729E-03/ DATA ((CALC(I,K,23),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.5720E-02, 0.6563E-02, + 0.7423E-02, 0.6563E-02, 0.4496E-01, 0.8437E-03, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,24),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1704E-02, 0.5681E-03, + 0.2840E-02, 0.5681E-03, 0.6817E-02, 0.2441E-03, + 0.0000E+00, 0.0000E+00, 0.1424E-02, 0.7722E-03, + 0.2968E-02, 0.7722E-03, 0.8076E-02, 0.2044E-03, + 0.0000E+00, 0.0000E+00, 0.1123E-02, 0.9913E-03, + 0.3105E-02, 0.9913E-03, 0.9428E-02, 0.1806E-03, + 0.0000E+00, 0.0000E+00, 0.8029E-03, 0.1225E-02, + 0.3254E-02, 0.1225E-02, 0.1087E-01, 0.1695E-03, + 0.0000E+00, 0.0000E+00, 0.4652E-03, 0.1475E-02, + 0.3414E-02, 0.1475E-02, 0.1242E-01, 0.1688E-03, + 0.0000E+00, 0.0000E+00, 0.1104E-03, 0.1739E-02, + 0.3588E-02, 0.1739E-02, 0.1406E-01, 0.1766E-03, + 0.0000E+00, 0.0000E+00, -0.2620E-03, 0.2020E-02, + 0.3777E-02, 0.2020E-02, 0.1581E-01, 0.1913E-03, + 0.0000E+00, 0.0000E+00, -0.6507E-03, 0.2317E-02, + 0.3982E-02, 0.2317E-02, 0.1767E-01, 0.2115E-03, + 0.0000E+00, 0.0000E+00, -0.1057E-02, 0.2631E-02, + 0.4205E-02, 0.2631E-02, 0.1964E-01, 0.2364E-03/ DATA ((CALC(I,K,24),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, -0.1480E-02, 0.2963E-02, + 0.4447E-02, 0.2963E-02, 0.2173E-01, 0.2650E-03, + 0.0000E+00, 0.0000E+00, -0.1920E-02, 0.3315E-02, + 0.4710E-02, 0.3315E-02, 0.2396E-01, 0.2970E-03, + 0.0000E+00, 0.0000E+00, -0.2380E-02, 0.3687E-02, + 0.4995E-02, 0.3687E-02, 0.2633E-01, 0.3317E-03, + 0.0000E+00, 0.0000E+00, -0.2859E-02, 0.4082E-02, + 0.5305E-02, 0.4082E-02, 0.2885E-01, 0.3689E-03, + 0.0000E+00, 0.0000E+00, -0.3359E-02, 0.4501E-02, + 0.5642E-02, 0.4501E-02, 0.3153E-01, 0.4087E-03, + 0.0000E+00, 0.0000E+00, -0.3883E-02, 0.4945E-02, + 0.6008E-02, 0.4945E-02, 0.3439E-01, 0.4508E-03, + 0.0000E+00, 0.0000E+00, -0.4431E-02, 0.5418E-02, + 0.6405E-02, 0.5418E-02, 0.3744E-01, 0.4954E-03, + 0.0000E+00, 0.0000E+00, -0.5007E-02, 0.5921E-02, + 0.6836E-02, 0.5921E-02, 0.4070E-01, 0.5424E-03, + 0.0000E+00, 0.0000E+00, -0.5612E-02, 0.6458E-02, + 0.7304E-02, 0.6458E-02, 0.4419E-01, 0.5924E-03/ DATA ((CALC(I,K,24),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.6249E-02, 0.7031E-02, + 0.7811E-02, 0.7031E-02, 0.4791E-01, 0.6447E-03, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,25),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.1862E-02, 0.6211E-03, + 0.3104E-02, 0.6211E-03, 0.7449E-02, 0.1673E-03, + 0.0000E+00, 0.0000E+00, 0.1518E-02, 0.8583E-03, + 0.3234E-02, 0.8583E-03, 0.8888E-02, 0.1390E-03, + 0.0000E+00, 0.0000E+00, 0.1156E-02, 0.1109E-02, + 0.3373E-02, 0.1109E-02, 0.1041E-01, 0.1237E-03, + 0.0000E+00, 0.0000E+00, 0.7779E-03, 0.1374E-02, + 0.3523E-02, 0.1374E-02, 0.1202E-01, 0.1183E-03, + 0.0000E+00, 0.0000E+00, 0.3848E-03, 0.1651E-02, + 0.3686E-02, 0.1651E-02, 0.1372E-01, 0.1208E-03, + 0.0000E+00, 0.0000E+00, -0.2175E-04, 0.1942E-02, + 0.3862E-02, 0.1942E-02, 0.1551E-01, 0.1295E-03, + 0.0000E+00, 0.0000E+00, -0.4424E-03, 0.2248E-02, + 0.4053E-02, 0.2248E-02, 0.1739E-01, 0.1429E-03, + 0.0000E+00, 0.0000E+00, -0.8755E-03, 0.2568E-02, + 0.4259E-02, 0.2568E-02, 0.1937E-01, 0.1600E-03, + 0.0000E+00, 0.0000E+00, -0.1322E-02, 0.2904E-02, + 0.4485E-02, 0.2904E-02, 0.2147E-01, 0.1801E-03/ DATA ((CALC(I,K,25),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, -0.1782E-02, 0.3256E-02, + 0.4730E-02, 0.3256E-02, 0.2367E-01, 0.2026E-03, + 0.0000E+00, 0.0000E+00, -0.2256E-02, 0.3627E-02, + 0.4997E-02, 0.3627E-02, 0.2601E-01, 0.2271E-03, + 0.0000E+00, 0.0000E+00, -0.2746E-02, 0.4017E-02, + 0.5288E-02, 0.4017E-02, 0.2847E-01, 0.2532E-03, + 0.0000E+00, 0.0000E+00, -0.3252E-02, 0.4427E-02, + 0.5603E-02, 0.4427E-02, 0.3109E-01, 0.2808E-03, + 0.0000E+00, 0.0000E+00, -0.3777E-02, 0.4861E-02, + 0.5946E-02, 0.4861E-02, 0.3386E-01, 0.3099E-03, + 0.0000E+00, 0.0000E+00, -0.4323E-02, 0.5320E-02, + 0.6318E-02, 0.5320E-02, 0.3680E-01, 0.3404E-03, + 0.0000E+00, 0.0000E+00, -0.4891E-02, 0.5807E-02, + 0.6723E-02, 0.5807E-02, 0.3994E-01, 0.3725E-03, + 0.0000E+00, 0.0000E+00, -0.5484E-02, 0.6323E-02, + 0.7162E-02, 0.6323E-02, 0.4327E-01, 0.4061E-03, + 0.0000E+00, 0.0000E+00, -0.6105E-02, 0.6873E-02, + 0.7639E-02, 0.6873E-02, 0.4684E-01, 0.4420E-03/ DATA ((CALC(I,K,25),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.6757E-02, 0.7458E-02, + 0.8157E-02, 0.7458E-02, 0.5065E-01, 0.4794E-03, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,26),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.2044E-02, 0.6816E-03, + 0.3407E-02, 0.6816E-03, 0.8178E-02, 0.1109E-03, + 0.0000E+00, 0.0000E+00, 0.1620E-02, 0.9567E-03, + 0.3533E-02, 0.9567E-03, 0.9813E-02, 0.9210E-04, + 0.0000E+00, 0.0000E+00, 0.1184E-02, 0.1243E-02, + 0.3668E-02, 0.1243E-02, 0.1152E-01, 0.8342E-04, + 0.0000E+00, 0.0000E+00, 0.7369E-03, 0.1539E-02, + 0.3814E-02, 0.1539E-02, 0.1329E-01, 0.8212E-04, + 0.0000E+00, 0.0000E+00, 0.2815E-03, 0.1846E-02, + 0.3973E-02, 0.1846E-02, 0.1514E-01, 0.8635E-04, + 0.0000E+00, 0.0000E+00, -0.1817E-03, 0.2164E-02, + 0.4145E-02, 0.2164E-02, 0.1707E-01, 0.9473E-04, + 0.0000E+00, 0.0000E+00, -0.6536E-03, 0.2493E-02, + 0.4333E-02, 0.2493E-02, 0.1907E-01, 0.1062E-03, + 0.0000E+00, 0.0000E+00, -0.1132E-02, 0.2834E-02, + 0.4537E-02, 0.2834E-02, 0.2117E-01, 0.1199E-03, + 0.0000E+00, 0.0000E+00, -0.1618E-02, 0.3188E-02, + 0.4760E-02, 0.3188E-02, 0.2335E-01, 0.1354E-03/ DATA ((CALC(I,K,26),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, -0.2112E-02, 0.3557E-02, + 0.5004E-02, 0.3557E-02, 0.2564E-01, 0.1521E-03, + 0.0000E+00, 0.0000E+00, -0.2615E-02, 0.3941E-02, + 0.5270E-02, 0.3941E-02, 0.2805E-01, 0.1698E-03, + 0.0000E+00, 0.0000E+00, -0.3130E-02, 0.4343E-02, + 0.5561E-02, 0.4343E-02, 0.3058E-01, 0.1885E-03, + 0.0000E+00, 0.0000E+00, -0.3657E-02, 0.4765E-02, + 0.5876E-02, 0.4765E-02, 0.3325E-01, 0.2079E-03, + 0.0000E+00, 0.0000E+00, -0.4199E-02, 0.5208E-02, + 0.6219E-02, 0.5208E-02, 0.3607E-01, 0.2281E-03, + 0.0000E+00, 0.0000E+00, -0.4758E-02, 0.5674E-02, + 0.6592E-02, 0.5674E-02, 0.3905E-01, 0.2492E-03, + 0.0000E+00, 0.0000E+00, -0.5337E-02, 0.6167E-02, + 0.6998E-02, 0.6167E-02, 0.4222E-01, 0.2712E-03, + 0.0000E+00, 0.0000E+00, -0.5938E-02, 0.6688E-02, + 0.7440E-02, 0.6688E-02, 0.4559E-01, 0.2939E-03, + 0.0000E+00, 0.0000E+00, -0.6565E-02, 0.7241E-02, + 0.7919E-02, 0.7241E-02, 0.4919E-01, 0.3184E-03/ DATA ((CALC(I,K,26),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.7220E-02, 0.7829E-02, + 0.8440E-02, 0.7829E-02, 0.5300E-01, 0.3430E-03, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,27),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.2253E-02, 0.7503E-03, + 0.3754E-02, 0.7503E-03, 0.9011E-02, 0.7114E-04, + 0.0000E+00, 0.0000E+00, 0.1729E-02, 0.1069E-02, + 0.3867E-02, 0.1069E-02, 0.1086E-01, 0.5971E-04, + 0.0000E+00, 0.0000E+00, 0.1203E-02, 0.1393E-02, + 0.3989E-02, 0.1393E-02, 0.1275E-01, 0.5573E-04, + 0.0000E+00, 0.0000E+00, 0.6762E-03, 0.1723E-02, + 0.4122E-02, 0.1723E-02, 0.1469E-01, 0.5689E-04, + 0.0000E+00, 0.0000E+00, 0.1505E-03, 0.2059E-02, + 0.4269E-02, 0.2059E-02, 0.1668E-01, 0.6163E-04, + 0.0000E+00, 0.0000E+00, -0.3748E-03, 0.2402E-02, + 0.4430E-02, 0.2402E-02, 0.1872E-01, 0.6894E-04, + 0.0000E+00, 0.0000E+00, -0.8993E-03, 0.2752E-02, + 0.4608E-02, 0.2752E-02, 0.2082E-01, 0.7800E-04, + 0.0000E+00, 0.0000E+00, -0.1421E-02, 0.3111E-02, + 0.4803E-02, 0.3111E-02, 0.2300E-01, 0.8825E-04, + 0.0000E+00, 0.0000E+00, -0.1944E-02, 0.3480E-02, + 0.5018E-02, 0.3480E-02, 0.2525E-01, 0.9931E-04/ DATA ((CALC(I,K,27),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, -0.2467E-02, 0.3859E-02, + 0.5253E-02, 0.3859E-02, 0.2759E-01, 0.1109E-03, + 0.0000E+00, 0.0000E+00, -0.2993E-02, 0.4252E-02, + 0.5513E-02, 0.4252E-02, 0.3003E-01, 0.1230E-03, + 0.0000E+00, 0.0000E+00, -0.3525E-02, 0.4659E-02, + 0.5794E-02, 0.4659E-02, 0.3258E-01, 0.1354E-03, + 0.0000E+00, 0.0000E+00, -0.4064E-02, 0.5083E-02, + 0.6104E-02, 0.5083E-02, 0.3525E-01, 0.1481E-03, + 0.0000E+00, 0.0000E+00, -0.4613E-02, 0.5527E-02, + 0.6442E-02, 0.5527E-02, 0.3807E-01, 0.1612E-03, + 0.0000E+00, 0.0000E+00, -0.5175E-02, 0.5992E-02, + 0.6810E-02, 0.5992E-02, 0.4103E-01, 0.1743E-03, + 0.0000E+00, 0.0000E+00, -0.5753E-02, 0.6482E-02, + 0.7212E-02, 0.6482E-02, 0.4419E-01, 0.1887E-03, + 0.0000E+00, 0.0000E+00, -0.6350E-02, 0.6999E-02, + 0.7648E-02, 0.6999E-02, 0.4753E-01, 0.2038E-03, + 0.0000E+00, 0.0000E+00, -0.6971E-02, 0.7547E-02, + 0.8123E-02, 0.7547E-02, 0.5107E-01, 0.2184E-03/ DATA ((CALC(I,K,27),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.7618E-02, 0.8128E-02, + 0.8639E-02, 0.8128E-02, 0.5487E-01, 0.2347E-03, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,28),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.2488E-02, 0.8293E-03, + 0.4148E-02, 0.8293E-03, 0.9954E-02, 0.4417E-04, + 0.0000E+00, 0.0000E+00, 0.1839E-02, 0.1197E-02, + 0.4234E-02, 0.1197E-02, 0.1203E-01, 0.3802E-04, + 0.0000E+00, 0.0000E+00, 0.1204E-02, 0.1564E-02, + 0.4331E-02, 0.1564E-02, 0.1411E-01, 0.3700E-04, + 0.0000E+00, 0.0000E+00, 0.5842E-03, 0.1929E-02, + 0.4442E-02, 0.1929E-02, 0.1621E-01, 0.3926E-04, + 0.0000E+00, 0.0000E+00, -0.2016E-04, 0.2294E-02, + 0.4567E-02, 0.2294E-02, 0.1832E-01, 0.4360E-04, + 0.0000E+00, 0.0000E+00, -0.6087E-03, 0.2660E-02, + 0.4709E-02, 0.2660E-02, 0.2046E-01, 0.4930E-04, + 0.0000E+00, 0.0000E+00, -0.1188E-02, 0.3027E-02, + 0.4869E-02, 0.3027E-02, 0.2264E-01, 0.5583E-04, + 0.0000E+00, 0.0000E+00, -0.1752E-02, 0.3398E-02, + 0.5048E-02, 0.3398E-02, 0.2485E-01, 0.6283E-04, + 0.0000E+00, 0.0000E+00, -0.2306E-02, 0.3774E-02, + 0.5247E-02, 0.3774E-02, 0.2713E-01, 0.7008E-04/ DATA ((CALC(I,K,28),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, -0.2851E-02, 0.4158E-02, + 0.5469E-02, 0.4158E-02, 0.2947E-01, 0.7746E-04, + 0.0000E+00, 0.0000E+00, -0.3391E-02, 0.4551E-02, + 0.5712E-02, 0.4551E-02, 0.3189E-01, 0.8492E-04, + 0.0000E+00, 0.0000E+00, -0.3929E-02, 0.4955E-02, + 0.5983E-02, 0.4955E-02, 0.3441E-01, 0.9244E-04, + 0.0000E+00, 0.0000E+00, -0.4468E-02, 0.5373E-02, + 0.6280E-02, 0.5373E-02, 0.3703E-01, 0.1000E-03, + 0.0000E+00, 0.0000E+00, -0.5012E-02, 0.5809E-02, + 0.6603E-02, 0.5809E-02, 0.3978E-01, 0.1078E-03, + 0.0000E+00, 0.0000E+00, -0.5563E-02, 0.6263E-02, + 0.6961E-02, 0.6263E-02, 0.4268E-01, 0.1157E-03, + 0.0000E+00, 0.0000E+00, -0.6127E-02, 0.6740E-02, + 0.7351E-02, 0.6740E-02, 0.4578E-01, 0.1250E-03, + 0.0000E+00, 0.0000E+00, -0.6706E-02, 0.7240E-02, + 0.7776E-02, 0.7240E-02, 0.4901E-01, 0.1327E-03, + 0.0000E+00, 0.0000E+00, -0.7305E-02, 0.7772E-02, + 0.8238E-02, 0.7772E-02, 0.5246E-01, 0.1414E-03/ DATA ((CALC(I,K,28),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.7929E-02, 0.8336E-02, + 0.8741E-02, 0.8336E-02, 0.5613E-01, 0.1509E-03, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,29),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.2758E-02, 0.9195E-03, + 0.4597E-02, 0.9195E-03, 0.1103E-01, 0.2644E-04, + 0.0000E+00, 0.0000E+00, 0.1946E-02, 0.1347E-02, + 0.4639E-02, 0.1347E-02, 0.1337E-01, 0.2369E-04, + 0.0000E+00, 0.0000E+00, 0.1176E-02, 0.1761E-02, + 0.4696E-02, 0.1761E-02, 0.1565E-01, 0.2419E-04, + 0.0000E+00, 0.0000E+00, 0.4450E-03, 0.2163E-02, + 0.4770E-02, 0.2163E-02, 0.1790E-01, 0.2654E-04, + 0.0000E+00, 0.0000E+00, -0.2480E-03, 0.2556E-02, + 0.4863E-02, 0.2556E-02, 0.2012E-01, 0.2988E-04, + 0.0000E+00, 0.0000E+00, -0.9065E-03, 0.2942E-02, + 0.4976E-02, 0.2942E-02, 0.2232E-01, 0.3379E-04, + 0.0000E+00, 0.0000E+00, -0.1552E-02, 0.3321E-02, + 0.5109E-02, 0.3321E-02, 0.2453E-01, 0.3796E-04, + 0.0000E+00, 0.0000E+00, -0.2148E-02, 0.3698E-02, + 0.5264E-02, 0.3698E-02, 0.2674E-01, 0.4218E-04, + 0.0000E+00, 0.0000E+00, -0.2720E-02, 0.4074E-02, + 0.5442E-02, 0.4074E-02, 0.2899E-01, 0.4636E-04/ DATA ((CALC(I,K,29),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, -0.3274E-02, 0.4453E-02, + 0.5644E-02, 0.4453E-02, 0.3127E-01, 0.5047E-04, + 0.0000E+00, 0.0000E+00, -0.3813E-02, 0.4836E-02, + 0.5871E-02, 0.4836E-02, 0.3362E-01, 0.5450E-04, + 0.0000E+00, 0.0000E+00, -0.4342E-02, 0.5227E-02, + 0.6123E-02, 0.5227E-02, 0.3604E-01, 0.5847E-04, + 0.0000E+00, 0.0000E+00, -0.4865E-02, 0.5629E-02, + 0.6401E-02, 0.5629E-02, 0.3856E-01, 0.6242E-04, + 0.0000E+00, 0.0000E+00, -0.5387E-02, 0.6045E-02, + 0.6710E-02, 0.6045E-02, 0.4119E-01, 0.6640E-04, + 0.0000E+00, 0.0000E+00, -0.5912E-02, 0.6478E-02, + 0.7047E-02, 0.6478E-02, 0.4391E-01, 0.6982E-04, + 0.0000E+00, 0.0000E+00, -0.6446E-02, 0.6931E-02, + 0.7417E-02, 0.6931E-02, 0.4686E-01, 0.7454E-04, + 0.0000E+00, 0.0000E+00, -0.6991E-02, 0.7403E-02, + 0.7821E-02, 0.7403E-02, 0.4993E-01, 0.7895E-04, + 0.0000E+00, 0.0000E+00, -0.7554E-02, 0.7907E-02, + 0.8260E-02, 0.7907E-02, 0.5320E-01, 0.8349E-04/ DATA ((CALC(I,K,29),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.8139E-02, 0.8437E-02, + 0.8739E-02, 0.8437E-02, 0.5667E-01, 0.8814E-04, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,30),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.3081E-02, 0.1028E-02, + 0.5135E-02, 0.1028E-02, 0.1233E-01, 0.1487E-04, + 0.0000E+00, 0.0000E+00, 0.2042E-02, 0.1531E-02, + 0.5103E-02, 0.1531E-02, 0.1497E-01, 0.1400E-04, + 0.0000E+00, 0.0000E+00, 0.1093E-02, 0.2002E-02, + 0.5097E-02, 0.2002E-02, 0.1747E-01, 0.1495E-04, + 0.0000E+00, 0.0000E+00, 0.2250E-03, 0.2446E-02, + 0.5116E-02, 0.2446E-02, 0.1987E-01, 0.1677E-04, + 0.0000E+00, 0.0000E+00, -0.5694E-03, 0.2867E-02, + 0.5162E-02, 0.2867E-02, 0.2217E-01, 0.1889E-04, + 0.0000E+00, 0.0000E+00, -0.1300E-02, 0.3267E-02, + 0.5234E-02, 0.3267E-02, 0.2440E-01, 0.2112E-04, + 0.0000E+00, 0.0000E+00, -0.1974E-02, 0.3653E-02, + 0.5332E-02, 0.3653E-02, 0.2659E-01, 0.2333E-04, + 0.0000E+00, 0.0000E+00, -0.2597E-02, 0.4027E-02, + 0.5456E-02, 0.4027E-02, 0.2875E-01, 0.2542E-04, + 0.0000E+00, 0.0000E+00, -0.3183E-02, 0.4393E-02, + 0.5606E-02, 0.4393E-02, 0.3090E-01, 0.2739E-04/ DATA ((CALC(I,K,30),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, -0.3733E-02, 0.4756E-02, + 0.5782E-02, 0.4756E-02, 0.3307E-01, 0.2923E-04, + 0.0000E+00, 0.0000E+00, -0.4257E-02, 0.5118E-02, + 0.5984E-02, 0.5118E-02, 0.3527E-01, 0.3097E-04, + 0.0000E+00, 0.0000E+00, -0.4761E-02, 0.5485E-02, + 0.6207E-02, 0.5485E-02, 0.3753E-01, 0.3263E-04, + 0.0000E+00, 0.0000E+00, -0.5252E-02, 0.5858E-02, + 0.6461E-02, 0.5858E-02, 0.3986E-01, 0.3426E-04, + 0.0000E+00, 0.0000E+00, -0.5735E-02, 0.6241E-02, + 0.6747E-02, 0.6241E-02, 0.4229E-01, 0.3587E-04, + 0.0000E+00, 0.0000E+00, -0.6217E-02, 0.6640E-02, + 0.7059E-02, 0.6640E-02, 0.4481E-01, 0.3728E-04, + 0.0000E+00, 0.0000E+00, -0.6703E-02, 0.7055E-02, + 0.7402E-02, 0.7055E-02, 0.4751E-01, 0.3921E-04, + 0.0000E+00, 0.0000E+00, -0.7198E-02, 0.7475E-02, + 0.7777E-02, 0.7475E-02, 0.5034E-01, 0.4109E-04, + 0.0000E+00, 0.0000E+00, -0.7706E-02, 0.7940E-02, + 0.8185E-02, 0.7940E-02, 0.5336E-01, 0.4273E-04/ DATA ((CALC(I,K,30),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.8242E-02, 0.8430E-02, + 0.8630E-02, 0.8430E-02, 0.5653E-01, 0.4452E-04, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,31),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.3507E-02, 0.1169E-02, + 0.5847E-02, 0.1169E-02, 0.1403E-01, 0.7133E-05, + 0.0000E+00, 0.0000E+00, 0.2113E-02, 0.1787E-02, + 0.5686E-02, 0.1787E-02, 0.1711E-01, 0.7043E-05, + 0.0000E+00, 0.0000E+00, 0.9042E-03, 0.2337E-02, + 0.5577E-02, 0.2337E-02, 0.1990E-01, 0.7782E-05, + 0.0000E+00, 0.0000E+00, -0.1506E-03, 0.2833E-02, + 0.5513E-02, 0.2833E-02, 0.2246E-01, 0.8803E-05, + 0.0000E+00, 0.0000E+00, -0.1071E-02, 0.3282E-02, + 0.5492E-02, 0.3282E-02, 0.2482E-01, 0.9785E-05, + 0.0000E+00, 0.0000E+00, -0.1889E-02, 0.3693E-02, + 0.5509E-02, 0.3693E-02, 0.2704E-01, 0.1070E-04, + 0.0000E+00, 0.0000E+00, -0.2617E-02, 0.4075E-02, + 0.5562E-02, 0.4075E-02, 0.2915E-01, 0.1153E-04, + 0.0000E+00, 0.0000E+00, -0.3237E-02, 0.4435E-02, + 0.5648E-02, 0.4435E-02, 0.3118E-01, 0.1225E-04, + 0.0000E+00, 0.0000E+00, -0.3810E-02, 0.4777E-02, + 0.5764E-02, 0.4777E-02, 0.3317E-01, 0.1286E-04/ DATA ((CALC(I,K,31),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, -0.4327E-02, 0.5110E-02, + 0.5910E-02, 0.5110E-02, 0.3513E-01, 0.1339E-04, + 0.0000E+00, 0.0000E+00, -0.4804E-02, 0.5436E-02, + 0.6084E-02, 0.5436E-02, 0.3711E-01, 0.1387E-04, + 0.0000E+00, 0.0000E+00, -0.5254E-02, 0.5762E-02, + 0.6277E-02, 0.5762E-02, 0.3911E-01, 0.1430E-04, + 0.0000E+00, 0.0000E+00, -0.5684E-02, 0.6092E-02, + 0.6501E-02, 0.6092E-02, 0.4117E-01, 0.1470E-04, + 0.0000E+00, 0.0000E+00, -0.6102E-02, 0.6428E-02, + 0.6764E-02, 0.6428E-02, 0.4330E-01, 0.1511E-04, + 0.0000E+00, 0.0000E+00, -0.6516E-02, 0.6777E-02, + 0.7036E-02, 0.6777E-02, 0.4552E-01, 0.1546E-04, + 0.0000E+00, 0.0000E+00, -0.6931E-02, 0.7140E-02, + 0.7343E-02, 0.7140E-02, 0.4759E-01, 0.1532E-04, + 0.0000E+00, 0.0000E+00, -0.7349E-02, 0.7518E-02, + 0.7679E-02, 0.7518E-02, 0.5011E-01, 0.1594E-04, + 0.0000E+00, 0.0000E+00, -0.7783E-02, 0.7926E-02, + 0.8045E-02, 0.7926E-02, 0.5287E-01, 0.1677E-04/ DATA ((CALC(I,K,31),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.8237E-02, 0.8375E-02, + 0.8443E-02, 0.8375E-02, 0.5570E-01, 0.1729E-04, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ DATA ((CALC(I,K,32),I=1,8), K=1,9) / + 0.0000E+00, 0.0000E+00, 0.4201E-02, 0.1401E-02, + 0.7002E-02, 0.1401E-02, 0.1681E-01, 0.1987E-05, + 0.0000E+00, 0.0000E+00, 0.2060E-02, 0.2256E-02, + 0.6566E-02, 0.2256E-02, 0.2079E-01, 0.2022E-05, + 0.0000E+00, 0.0000E+00, 0.3736E-03, 0.2950E-02, + 0.6263E-02, 0.2950E-02, 0.2408E-01, 0.2281E-05, + 0.0000E+00, 0.0000E+00, -0.9719E-03, 0.3522E-02, + 0.6064E-02, 0.3522E-02, 0.2687E-01, 0.2574E-05, + 0.0000E+00, 0.0000E+00, -0.2051E-02, 0.4003E-02, + 0.5947E-02, 0.4003E-02, 0.2928E-01, 0.2776E-05, + 0.0000E+00, 0.0000E+00, -0.2860E-02, 0.4414E-02, + 0.5898E-02, 0.4414E-02, 0.3140E-01, 0.2928E-05, + 0.0000E+00, 0.0000E+00, -0.3754E-02, 0.4773E-02, + 0.5902E-02, 0.4773E-02, 0.3332E-01, 0.3058E-05, + 0.0000E+00, 0.0000E+00, -0.4339E-02, 0.5094E-02, + 0.5952E-02, 0.5094E-02, 0.3510E-01, 0.3146E-05, + 0.0000E+00, 0.0000E+00, -0.4796E-02, 0.5389E-02, + 0.6041E-02, 0.5389E-02, 0.3679E-01, 0.3195E-05/ DATA ((CALC(I,K,32),I=1,8), K=10,18) / + 0.0000E+00, 0.0000E+00, -0.5223E-02, 0.5666E-02, + 0.6163E-02, 0.5666E-02, 0.3843E-01, 0.3228E-05, + 0.0000E+00, 0.0000E+00, -0.5601E-02, 0.5933E-02, + 0.6312E-02, 0.5933E-02, 0.4005E-01, 0.3250E-05, + 0.0000E+00, 0.0000E+00, -0.5946E-02, 0.6197E-02, + 0.6300E-02, 0.6197E-02, 0.4168E-01, 0.3262E-05, + 0.0000E+00, 0.0000E+00, -0.6272E-02, 0.6461E-02, + 0.6473E-02, 0.6461E-02, 0.4336E-01, 0.3272E-05, + 0.0000E+00, 0.0000E+00, -0.6588E-02, 0.6732E-02, + 0.6751E-02, 0.6732E-02, 0.4510E-01, 0.3289E-05, + 0.0000E+00, 0.0000E+00, -0.6902E-02, 0.7021E-02, + 0.7020E-02, 0.7021E-02, 0.4681E-01, 0.3209E-05, + 0.0000E+00, 0.0000E+00, -0.7219E-02, 0.7310E-02, + 0.7297E-02, 0.7310E-02, 0.4882E-01, 0.3295E-05, + 0.0000E+00, 0.0000E+00, -0.7523E-02, 0.7606E-02, + 0.7596E-02, 0.7606E-02, 0.5079E-01, 0.3324E-05, + 0.0000E+00, 0.0000E+00, -0.7851E-02, 0.7922E-02, + 0.7916E-02, 0.7922E-02, 0.5287E-01, 0.3328E-05/ DATA ((CALC(I,K,32),I=1,8), K=19,20) / + 0.0000E+00, 0.0000E+00, -0.8208E-02, 0.8263E-02, + 0.8256E-02, 0.8263E-02, 0.5513E-01, 0.3346E-05, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, + 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00/ C---------------------------------------------------------------------- DATA ISTART/0/ SAVE ISTART,OWLAM2,Q02,FLAV, /W5051I7/ C IF (ISTART.EQ.0) THEN ISTART=1 DO 10 K=1,32 DO 10 I=1,20 DO 10 M=1,8 10 CALCO(M,I,K) = CALC(M,I,K) OWLAM=PAR(1) OWLAM2=OWLAM**2 Q02=PAR(30) FLAV=PAR(25) DELTA=PAR(29) CALL WATE32 ENDIF C X = DX Q = DQ Q2 = Q*Q IDQ2=2 SB=0. IF(Q2-Q02) 1,1,2 2 IF(IDQ2-1) 1,1,3 3 SB= LOG( LOG( MAX(Q02,Q2)/OWLAM2)/ LOG(Q02/OWLAM2)) 1 CONTINUE CALL AURGAM(8,0,X,SB,XPDF(7)) CALL AURGAM(7,0,X,SB,SING) CALL AURGAM(4,0,X,SB,DPLUSNS) CALL AURGAM(3,0,X,SB,CPLUSNS) CALL AURGAM(5,0,X,SB,UPLUSNS) CALL AURGAM(6,0,X,SB,SPLUSNS) XPDF(3) = CPLUSNS XPDF(4) = DPLUSNS XPDF(5) = UPLUSNS XPDF(6) = SPLUSNS XPDF(1) = SING C ADD = XPDF(1)/FLAV UPLUS=XPDF(5)+ADD DPLUS=-XPDF(4)+ADD SPLUS=-XPDF(6)+ADD CPLUS=-XPDF(3)+ADD UB=UPLUS*0.5 UV=UB DB=DPLUS*0.5 DV=DB SB=SPLUS*0.5 CB=CPLUS*0.5 SING=XPDF(1) GLU=XPDF(7) GL=GLU * DUV=MAX(ZERO,UV) DDV=MAX(ZERO,DV) DUB=MAX(ZERO,UB) DDB=MAX(ZERO,DB) DSB=MAX(ZERO,SB) DCB=MAX(ZERO,CB) DGL=MAX(ZERO,GL) C RETURN END