]> git.uio.no Git - u/mrichter/AliRoot.git/blob - PDF/spdf/sfabfow.F
Adding some QCD diffractive states to the PDG list
[u/mrichter/AliRoot.git] / PDF / spdf / sfabfow.F
1 *
2 * $Id$
3 *
4 * $Log$
5 * Revision 1.1.1.2  1996/10/30 08:29:21  cernlib
6 * Version 7.04
7 *
8 * Revision 1.1.1.1  1996/04/12 15:29:29  plothow
9 * Version 7.01
10 *
11 *
12 #include "pdf/pilot.h"
13       SUBROUTINE SFABFOW(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DGL)
14 #include "pdf/expdp.inc"
15      +       PAR(40),CALC(8,20,25),CALCO(8,20,25),ZEROD,
16      +       DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DGL
17       REAL    X, Q, UPV, DNV, USEA, DSEA, STR, CHM, GL
18       COMMON/W5051I3/CALCO
19       REAL XPDF(7)
20       DATA ZEROD/0.D0/
21 C ---------------------------------------------------------------------
22         DATA (PAR(I),I=1,40) /
23      +  0.23000000D+00, 0.87291000D+00, 0.30563000D+01,-0.81541000D+00,
24      +  0.00000000D+00, 0.10280000D+01, 0.53963000D+01, 0.00000000D+00,
25      +  0.11002000D+01, 0.87514000D+01, 0.21308913D+01, 0.00000000D+00,
26      +  0.40000000D+01, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
27      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
28      +  0.10000000D+01, 0.10000000D+01, 0.92300000D+00, 0.11240000D+01,
29      +  0.10000000D+01, 0.10000000D+01, 0.10000000D+01, 0.10520000D+01,
30      +  0.18200000D+01, 0.16496000D+03, 0.00000000D+00, 0.00000000D+00,
31      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
32      +  0.10000000D+01, 0.20000000D-07, 0.20000000D+01, 0.10000000D+09/
33         DATA ((CALC(I,K,1 ),I=1,8), K=1,9) /
34      +  0.20194028D+00, 0.18738953D+00, 0.18176816D+00, 0.18176816D+00,
35      +  0.18176816D+00, 0.18176816D+00, 0.21223805D+01,-0.17347235D-17,
36      +  0.24423567D+00, 0.22842866D+00, 0.22167093D+00, 0.22163287D+00,
37      +  0.22164650D+00, 0.22167093D+00, 0.33355957D+01, 0.40286953D-01,
38      +  0.29313746D+00, 0.27603823D+00, 0.26809481D+00, 0.26802280D+00,
39      +  0.26804826D+00, 0.26809481D+00, 0.47295147D+01, 0.87113536D-01,
40      +  0.34896078D+00, 0.33053340D+00, 0.32135373D+00, 0.32125188D+00,
41      +  0.32128744D+00, 0.32135373D+00, 0.63030943D+01, 0.14079460D+00,
42      +  0.41196816D+00, 0.39217673D+00, 0.38170912D+00, 0.38158152D+00,
43      +  0.38162550D+00, 0.38170912D+00, 0.80526510D+01, 0.20159242D+00,
44      +  0.48239147D+00, 0.46120006D+00, 0.44939173D+00, 0.44924244D+00,
45      +  0.44929323D+00, 0.44939173D+00, 0.99718541D+01, 0.26973861D+00,
46      +  0.56042446D+00, 0.53779719D+00, 0.52459425D+00, 0.52442730D+00,
47      +  0.52448335D+00, 0.52459425D+00, 0.12052181D+02, 0.34542675D+00,
48      +  0.64621175D+00, 0.62211276D+00, 0.60746034D+00, 0.60727971D+00,
49      +  0.60733952D+00, 0.60746034D+00, 0.14282793D+02, 0.42880139D+00,
50      +  0.73984307D+00, 0.71423660D+00, 0.69807888D+00, 0.69788851D+00,
51      +  0.69795065D+00, 0.69807888D+00, 0.16650759D+02, 0.51995229D+00/
52         DATA ((CALC(I,K,1 ),I=1,8), K=10,18) /
53      +  0.84134569D+00, 0.81419609D+00, 0.79647647D+00, 0.79628021D+00,
54      +  0.79634332D+00, 0.79647647D+00, 0.19141215D+02, 0.61890693D+00,
55      +  0.95067917D+00, 0.92195096D+00, 0.90261215D+00, 0.90241381D+00,
56      +  0.90247660D+00, 0.90261215D+00, 0.21737598D+02, 0.72562522D+00,
57      +  0.10677309D+01, 0.10373889D+01, 0.10163730D+01, 0.10161763D+01,
58      +  0.10162376D+01, 0.10163730D+01, 0.24421896D+02, 0.83999511D+00,
59      +  0.11923132D+01, 0.11603224D+01, 0.11375714D+01, 0.11373799D+01,
60      +  0.11374385D+01, 0.11375714D+01, 0.27174942D+02, 0.96182969D+00,
61      +  0.13241617D+01, 0.12904878D+01, 0.12659430D+01, 0.12657604D+01,
62      +  0.12658152D+01, 0.12659430D+01, 0.29976710D+02, 0.10908656D+01,
63      +  0.14629355D+01, 0.14275444D+01, 0.14011475D+01, 0.14009772D+01,
64      +  0.14010272D+01, 0.14011475D+01, 0.32806631D+02, 0.12267633D+01,
65      +  0.16082186D+01, 0.15710768D+01, 0.15427697D+01, 0.15426149D+01,
66      +  0.15426593D+01, 0.15427697D+01, 0.35643910D+02, 0.13691081D+01,
67      +  0.17595228D+01, 0.17205977D+01, 0.16903224D+01, 0.16901864D+01,
68      +  0.16902243D+01, 0.16903224D+01, 0.38467840D+02, 0.15174136D+01,
69      +  0.19162923D+01, 0.18755519D+01, 0.18432512D+01, 0.18431371D+01,
70      +  0.18431678D+01, 0.18432512D+01, 0.41258096D+02, 0.16711262D+01/
71         DATA ((CALC(I,K,1 ),I=1,8), K=19,20) /
72      +  0.20779090D+01, 0.20353224D+01, 0.20009399D+01, 0.20008506D+01,
73      +  0.20008735D+01, 0.20009399D+01, 0.43995019D+02, 0.18296301D+01,
74      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
75      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
76         DATA ((CALC(I,K,2 ),I=1,8), K=1,9) /
77      +  0.21583491D+00, 0.19158330D+00, 0.18018200D+00, 0.18018200D+00,
78      +  0.18018200D+00, 0.18018200D+00, 0.21138952D+01,-0.52041704D-17,
79      +  0.25255607D+00, 0.22649496D+00, 0.21327094D+00, 0.21323316D+00,
80      +  0.21324695D+00, 0.21327094D+00, 0.31115744D+01, 0.33702174D-01,
81      +  0.29378966D+00, 0.26588969D+00, 0.25077892D+00, 0.25070944D+00,
82      +  0.25073454D+00, 0.25077892D+00, 0.42176629D+01, 0.71850099D-01,
83      +  0.33964459D+00, 0.30987604D+00, 0.29281357D+00, 0.29271828D+00,
84      +  0.29275242D+00, 0.29281357D+00, 0.54262314D+01, 0.11455242D+00,
85      +  0.39018756D+00, 0.35852046D+00, 0.33944066D+00, 0.33932526D+00,
86      +  0.33936635D+00, 0.33944066D+00, 0.67303079D+01, 0.16187578D+00,
87      +  0.44544250D+00, 0.41184667D+00, 0.39068345D+00, 0.39055345D+00,
88      +  0.39059953D+00, 0.39068345D+00, 0.81210348D+01, 0.21384434D+00,
89      +  0.50539764D+00, 0.46984285D+00, 0.44652981D+00, 0.44639052D+00,
90      +  0.44643980D+00, 0.44652981D+00, 0.95883720D+01, 0.27044684D+00,
91      +  0.56999702D+00, 0.53245303D+00, 0.50692372D+00, 0.50678027D+00,
92      +  0.50683109D+00, 0.50692372D+00, 0.11121064D+02, 0.33162813D+00,
93      +  0.63914062D+00, 0.59957732D+00, 0.57176547D+00, 0.57162277D+00,
94      +  0.57167363D+00, 0.57176547D+00, 0.12706903D+02, 0.39728933D+00/
95         DATA ((CALC(I,K,2 ),I=1,8), K=10,18) /
96      +  0.71268211D+00, 0.67106968D+00, 0.64090944D+00, 0.64077222D+00,
97      +  0.64082174D+00, 0.64090944D+00, 0.14332935D+02, 0.46728562D+00,
98      +  0.79042864D+00, 0.74673766D+00, 0.71416390D+00, 0.71403665D+00,
99      +  0.71408358D+00, 0.71416390D+00, 0.15985692D+02, 0.54142599D+00,
100      +  0.87214117D+00, 0.82634277D+00, 0.79129132D+00, 0.79117834D+00,
101      +  0.79122157D+00, 0.79129132D+00, 0.17651425D+02, 0.61947365D+00,
102      +  0.95753582D+00, 0.90960182D+00, 0.87200980D+00, 0.87191517D+00,
103      +  0.87195370D+00, 0.87200980D+00, 0.19316334D+02, 0.70114732D+00,
104      +  0.10462861D+01, 0.99618915D+00, 0.95599526D+00, 0.95592283D+00,
105      +  0.95595580D+00, 0.95599526D+00, 0.20966798D+02, 0.78612350D+00,
106      +  0.11380258D+01, 0.10857397D+01, 0.10428845D+01, 0.10428379D+01,
107      +  0.10428646D+01, 0.10428845D+01, 0.22589582D+02, 0.87403953D+00,
108      +  0.12323534D+01, 0.11778530D+01, 0.11322792D+01, 0.11322618D+01,
109      +  0.11322815D+01, 0.11322792D+01, 0.24172034D+02, 0.96449744D+00,
110      +  0.13288358D+01, 0.12720973D+01, 0.12237501D+01, 0.12237652D+01,
111      +  0.12237774D+01, 0.12237501D+01, 0.25702254D+02, 0.10570685D+01,
112      +  0.14270141D+01, 0.13680153D+01, 0.13168427D+01, 0.13168932D+01,
113      +  0.13168974D+01, 0.13168427D+01, 0.27169250D+02, 0.11512983D+01/
114         DATA ((CALC(I,K,2 ),I=1,8), K=19,20) /
115      +  0.15264089D+01, 0.14651291D+01, 0.14110822D+01, 0.14111709D+01,
116      +  0.14111667D+01, 0.14110822D+01, 0.28563054D+02, 0.12467124D+01,
117      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
118      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
119         DATA ((CALC(I,K,3 ),I=1,8), K=1,9) /
120      +  0.23955354D+00, 0.20004558D+00, 0.17704643D+00, 0.17704643D+00,
121      +  0.17704643D+00, 0.17704643D+00, 0.20970011D+01, 0.52041704D-17,
122      +  0.27130160D+00, 0.22936997D+00, 0.20347723D+00, 0.20345215D+00,
123      +  0.20346114D+00, 0.20347723D+00, 0.28837840D+01, 0.27406605D-01,
124      +  0.30603895D+00, 0.26166466D+00, 0.23280614D+00, 0.23276307D+00,
125      +  0.23277886D+00, 0.23280614D+00, 0.37286621D+01, 0.57746771D-01,
126      +  0.34371025D+00, 0.29687360D+00, 0.26497877D+00, 0.26492438D+00,
127      +  0.26494503D+00, 0.26497877D+00, 0.46232124D+01, 0.90966868D-01,
128      +  0.38424325D+00, 0.33492418D+00, 0.29992417D+00, 0.29986475D+00,
129      +  0.29988858D+00, 0.29992417D+00, 0.55592695D+01, 0.12699659D+00,
130      +  0.42752842D+00, 0.37570676D+00, 0.33753472D+00, 0.33747622D+00,
131      +  0.33750172D+00, 0.33753472D+00, 0.65277378D+01, 0.16572895D+00,
132      +  0.47343421D+00, 0.41908995D+00, 0.37768129D+00, 0.37762929D+00,
133      +  0.37765516D+00, 0.37768129D+00, 0.75192617D+01, 0.20703539D+00,
134      +  0.52180091D+00, 0.46491448D+00, 0.42020721D+00, 0.42016694D+00,
135      +  0.42019206D+00, 0.42020721D+00, 0.85242535D+01, 0.25075977D+00,
136      +  0.57244319D+00, 0.51299568D+00, 0.46493073D+00, 0.46490708D+00,
137      +  0.46493048D+00, 0.46493073D+00, 0.95330873D+01, 0.29672077D+00/
138         DATA ((CALC(I,K,3 ),I=1,8), K=10,18) /
139      +  0.62515070D+00, 0.56312415D+00, 0.51164567D+00, 0.51164319D+00,
140      +  0.51166406D+00, 0.51164567D+00, 0.10536268D+02, 0.34471262D+00,
141      +  0.67968980D+00, 0.61506741D+00, 0.56012316D+00, 0.56014606D+00,
142      +  0.56016373D+00, 0.56012316D+00, 0.11524580D+02, 0.39450672D+00,
143      +  0.73580569D+00, 0.66857209D+00, 0.61011374D+00, 0.61016592D+00,
144      +  0.61017983D+00, 0.61011374D+00, 0.12489237D+02, 0.44585384D+00,
145      +  0.79322499D+00, 0.72336646D+00, 0.66134996D+00, 0.66143497D+00,
146      +  0.66144468D+00, 0.66134996D+00, 0.13421999D+02, 0.49848661D+00,
147      +  0.85165870D+00, 0.77916341D+00, 0.71354930D+00, 0.71367041D+00,
148      +  0.71367557D+00, 0.71354930D+00, 0.14315288D+02, 0.55212255D+00,
149      +  0.91080557D+00, 0.83566379D+00, 0.76641753D+00, 0.76657768D+00,
150      +  0.76657805D+00, 0.76641753D+00, 0.15162271D+02, 0.60646732D+00,
151      +  0.97035560D+00, 0.89255991D+00, 0.81965219D+00, 0.81985402D+00,
152      +  0.81984942D+00, 0.81965219D+00, 0.15956931D+02, 0.66121826D+00,
153      +  0.10299938D+01, 0.94953929D+00, 0.87294633D+00, 0.87319219D+00,
154      +  0.87318251D+00, 0.87294633D+00, 0.16694112D+02, 0.71606809D+00,
155      +  0.10894039D+01, 0.10062884D+01, 0.92599221D+00, 0.92628417D+00,
156      +  0.92626936D+00, 0.92599221D+00, 0.17369551D+02, 0.77070866D+00/
157         DATA ((CALC(I,K,3 ),I=1,8), K=19,20) /
158      +  0.11482723D+01, 0.10624964D+01, 0.97848509D+00, 0.97882493D+00,
159      +  0.97880498D+00, 0.97848509D+00, 0.17979883D+02, 0.82483467D+00,
160      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
161      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
162         DATA ((CALC(I,K,4 ),I=1,8), K=1,9) /
163      +  0.27919948D+00, 0.21681220D+00, 0.17091990D+00, 0.17091990D+00,
164      +  0.17091990D+00, 0.17091990D+00, 0.20635167D+01,-0.86736174D-17,
165      +  0.30719611D+00, 0.24175123D+00, 0.19141648D+00, 0.19141170D+00,
166      +  0.19141779D+00, 0.19141648D+00, 0.26646269D+01, 0.22020703D-01,
167      +  0.33688109D+00, 0.26837614D+00, 0.21356081D+00, 0.21355815D+00,
168      +  0.21356933D+00, 0.21356081D+00, 0.32852034D+01, 0.45728609D-01,
169      +  0.36811194D+00, 0.29654497D+00, 0.23721569D+00, 0.23722139D+00,
170      +  0.23723688D+00, 0.23721569D+00, 0.39172171D+01, 0.70986696D-01,
171      +  0.40075271D+00, 0.32612258D+00, 0.26225048D+00, 0.26227027D+00,
172      +  0.26228943D+00, 0.26225048D+00, 0.45537009D+01, 0.97664462D-01,
173      +  0.43464550D+00, 0.35695222D+00, 0.28851342D+00, 0.28855254D+00,
174      +  0.28857485D+00, 0.28851342D+00, 0.51876525D+01, 0.12561019D+00,
175      +  0.46962465D+00, 0.38886971D+00, 0.31584569D+00, 0.31590887D+00,
176      +  0.31593392D+00, 0.31584569D+00, 0.58124017D+01, 0.15466500D+00,
177      +  0.50551327D+00, 0.42169989D+00, 0.34407789D+00, 0.34416942D+00,
178      +  0.34419687D+00, 0.34407789D+00, 0.64216588D+01, 0.18465932D+00,
179      +  0.54212547D+00, 0.45525893D+00, 0.37303230D+00, 0.37315602D+00,
180      +  0.37318560D+00, 0.37303230D+00, 0.70095902D+01, 0.21541516D+00/
181         DATA ((CALC(I,K,4 ),I=1,8), K=10,18) /
182      +  0.57926783D+00, 0.48935570D+00, 0.40252430D+00, 0.40268365D+00,
183      +  0.40271512D+00, 0.40252430D+00, 0.75708958D+01, 0.24674751D+00,
184      +  0.61674122D+00, 0.52379364D+00, 0.43236416D+00, 0.43256218D+00,
185      +  0.43259533D+00, 0.43236416D+00, 0.81008614D+01, 0.27846612D+00,
186      +  0.65434284D+00, 0.55837275D+00, 0.46235905D+00, 0.46259841D+00,
187      +  0.46263305D+00, 0.46235905D+00, 0.85954001D+01, 0.31037751D+00,
188      +  0.69186835D+00, 0.59289174D+00, 0.49231516D+00, 0.49259814D+00,
189      +  0.49263410D+00, 0.49231516D+00, 0.90510786D+01, 0.34228709D+00,
190      +  0.72911419D+00, 0.62715029D+00, 0.52203989D+00, 0.52236843D+00,
191      +  0.52240554D+00, 0.52203989D+00, 0.94651297D+01, 0.37400133D+00,
192      +  0.76587978D+00, 0.66095125D+00, 0.55134404D+00, 0.55171978D+00,
193      +  0.55175786D+00, 0.55134404D+00, 0.98354508D+01, 0.40533002D+00,
194      +  0.80196976D+00, 0.69410291D+00, 0.58004405D+00, 0.58046828D+00,
195      +  0.58050714D+00, 0.58004405D+00, 0.10160591D+02, 0.43608839D+00,
196      +  0.83719619D+00, 0.72642107D+00, 0.60796400D+00, 0.60843772D+00,
197      +  0.60847716D+00, 0.60796400D+00, 0.10439724D+02, 0.46609923D+00,
198      +  0.87138042D+00, 0.75773103D+00, 0.63493754D+00, 0.63546150D+00,
199      +  0.63550128D+00, 0.63493754D+00, 0.10672617D+02, 0.49519478D+00/
200         DATA ((CALC(I,K,4 ),I=1,8), K=19,20) /
201      +  0.90435501D+00, 0.78786935D+00, 0.66080966D+00, 0.66138430D+00,
202      +  0.66142418D+00, 0.66080966D+00, 0.10859586D+02, 0.52321845D+00,
203      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
204      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
205         DATA ((CALC(I,K,5 ),I=1,8), K=1,9) /
206      +  0.34327808D+00, 0.24880963D+00, 0.15922736D+00, 0.15922736D+00,
207      +  0.15922736D+00, 0.15922736D+00, 0.19977525D+01, 0.00000000D+00,
208      +  0.36786460D+00, 0.26977928D+00, 0.17402038D+00, 0.17403996D+00,
209      +  0.17405794D+00, 0.17402038D+00, 0.24242245D+01, 0.17061951D-01,
210      +  0.39300678D+00, 0.29134137D+00, 0.18947912D+00, 0.18952311D+00,
211      +  0.18955798D+00, 0.18947912D+00, 0.28433784D+01, 0.34819499D-01,
212      +  0.41855640D+00, 0.31334873D+00, 0.20546319D+00, 0.20553578D+00,
213      +  0.20558647D+00, 0.20546319D+00, 0.32499600D+01, 0.53131698D-01,
214      +  0.44438149D+00, 0.33567084D+00, 0.22184868D+00, 0.22195352D+00,
215      +  0.22201896D+00, 0.22184868D+00, 0.36398535D+01, 0.71874004D-01,
216      +  0.47034243D+00, 0.35816994D+00, 0.23850495D+00, 0.23864521D+00,
217      +  0.23872434D+00, 0.23850495D+00, 0.40093587D+01, 0.90915002D-01,
218      +  0.49630067D+00, 0.38070976D+00, 0.25530329D+00, 0.25548169D+00,
219      +  0.25557342D+00, 0.25530329D+00, 0.43552632D+01, 0.11012505D+00,
220      +  0.52211769D+00, 0.40315441D+00, 0.27211578D+00, 0.27233459D+00,
221      +  0.27243782D+00, 0.27211578D+00, 0.46748810D+01, 0.12937513D+00,
222      +  0.54765633D+00, 0.42536971D+00, 0.28881646D+00, 0.28907758D+00,
223      +  0.28919119D+00, 0.28881646D+00, 0.49660412D+01, 0.14853807D+00/
224         DATA ((CALC(I,K,5 ),I=1,8), K=10,18) /
225      +  0.57278214D+00, 0.44722448D+00, 0.30528267D+00, 0.30558761D+00,
226      +  0.30571046D+00, 0.30528267D+00, 0.52270916D+01, 0.16748981D+00,
227      +  0.59736451D+00, 0.46859171D+00, 0.32139606D+00, 0.32174599D+00,
228      +  0.32187692D+00, 0.32139606D+00, 0.54568849D+01, 0.18611047D+00,
229      +  0.62127804D+00, 0.48934981D+00, 0.33704382D+00, 0.33743958D+00,
230      +  0.33757743D+00, 0.33704382D+00, 0.56547590D+01, 0.20428558D+00,
231      +  0.64440369D+00, 0.50938381D+00, 0.35211979D+00, 0.35256193D+00,
232      +  0.35270551D+00, 0.35211979D+00, 0.58205086D+01, 0.22190716D+00,
233      +  0.66662991D+00, 0.52858643D+00, 0.36652544D+00, 0.36701423D+00,
234      +  0.36716235D+00, 0.36652544D+00, 0.59543523D+01, 0.23887477D+00,
235      +  0.68785373D+00, 0.54685908D+00, 0.38017085D+00, 0.38070628D+00,
236      +  0.38085774D+00, 0.38017085D+00, 0.60568930D+01, 0.25509642D+00,
237      +  0.70798155D+00, 0.56411274D+00, 0.39297543D+00, 0.39355726D+00,
238      +  0.39371088D+00, 0.39297543D+00, 0.61290758D+01, 0.27048936D+00,
239      +  0.72692995D+00, 0.58026859D+00, 0.40486854D+00, 0.40549631D+00,
240      +  0.40565090D+00, 0.40486854D+00, 0.61721427D+01, 0.28498071D+00,
241      +  0.74462615D+00, 0.59525856D+00, 0.41578998D+00, 0.41646301D+00,
242      +  0.41661740D+00, 0.41578998D+00, 0.61875858D+01, 0.29850793D+00/
243         DATA ((CALC(I,K,5 ),I=1,8), K=19,20) /
244      +  0.76100847D+00, 0.60902567D+00, 0.42569018D+00, 0.42640762D+00,
245      +  0.42656065D+00, 0.42569018D+00, 0.61771012D+01, 0.31101904D+00,
246      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
247      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
248         DATA ((CALC(I,K,6 ),I=1,8), K=1,9) /
249      +  0.44106537D+00, 0.30515567D+00, 0.13794628D+00, 0.13794628D+00,
250      +  0.13794628D+00, 0.13794628D+00, 0.18709523D+01,-0.17347235D-17,
251      +  0.46119251D+00, 0.32136177D+00, 0.14716767D+00, 0.14720135D+00,
252      +  0.14723645D+00, 0.14716767D+00, 0.21299443D+01, 0.12306566D-01,
253      +  0.48095509D+00, 0.33729624D+00, 0.15641371D+00, 0.15648271D+00,
254      +  0.15654848D+00, 0.15641371D+00, 0.23687832D+01, 0.24632075D-01,
255      +  0.50024830D+00, 0.35285650D+00, 0.16558478D+00, 0.16569038D+00,
256      +  0.16578272D+00, 0.16558478D+00, 0.25855725D+01, 0.36876049D-01,
257      +  0.51898317D+00, 0.36795632D+00, 0.17459826D+00, 0.17474146D+00,
258      +  0.17485651D+00, 0.17459826D+00, 0.27791414D+01, 0.48954867D-01,
259      +  0.53707332D+00, 0.38251257D+00, 0.18337484D+00, 0.18355637D+00,
260      +  0.18369053D+00, 0.18337484D+00, 0.29487906D+01, 0.60788045D-01,
261      +  0.55443725D+00, 0.39644746D+00, 0.19184096D+00, 0.19206130D+00,
262      +  0.19221117D+00, 0.19184096D+00, 0.30941757D+01, 0.72300639D-01,
263      +  0.57099901D+00, 0.40968919D+00, 0.19992912D+00, 0.20018850D+00,
264      +  0.20035090D+00, 0.19992912D+00, 0.32153306D+01, 0.83423599D-01,
265      +  0.58668840D+00, 0.42217206D+00, 0.20757792D+00, 0.20787635D+00,
266      +  0.20804828D+00, 0.20757792D+00, 0.33126163D+01, 0.94093852D-01/
267         DATA ((CALC(I,K,6 ),I=1,8), K=10,18) /
268      +  0.60144177D+00, 0.43383715D+00, 0.21473262D+00, 0.21506991D+00,
269      +  0.21524858D+00, 0.21473262D+00, 0.33866886D+01, 0.10425484D+00,
270      +  0.61520235D+00, 0.44463269D+00, 0.22134538D+00, 0.22172112D+00,
271      +  0.22190390D+00, 0.22134538D+00, 0.34384597D+01, 0.11385675D+00,
272      +  0.62792073D+00, 0.45451441D+00, 0.22737543D+00, 0.22778906D+00,
273      +  0.22797350D+00, 0.22737543D+00, 0.34690589D+01, 0.12285682D+00,
274      +  0.63955514D+00, 0.46344576D+00, 0.23278927D+00, 0.23324005D+00,
275      +  0.23342389D+00, 0.23278927D+00, 0.34797928D+01, 0.13121941D+00,
276      +  0.65007163D+00, 0.47139806D+00, 0.23756070D+00, 0.23804772D+00,
277      +  0.23822885D+00, 0.23756070D+00, 0.34721068D+01, 0.13891615D+00,
278      +  0.65944422D+00, 0.47835053D+00, 0.24167075D+00, 0.24219297D+00,
279      +  0.24236945D+00, 0.24167075D+00, 0.34475470D+01, 0.14592585D+00,
280      +  0.66765484D+00, 0.48429022D+00, 0.24510753D+00, 0.24566378D+00,
281      +  0.24583383D+00, 0.24510753D+00, 0.34077251D+01, 0.15223439D+00,
282      +  0.67469321D+00, 0.48921180D+00, 0.24786598D+00, 0.24845499D+00,
283      +  0.24861699D+00, 0.24786598D+00, 0.33542852D+01, 0.15783446D+00,
284      +  0.68055661D+00, 0.49311730D+00, 0.24994757D+00, 0.25056793D+00,
285      +  0.25072043D+00, 0.24994757D+00, 0.32888738D+01, 0.16272529D+00/
286         DATA ((CALC(I,K,6 ),I=1,8), K=19,20) /
287      +  0.68524957D+00, 0.49601572D+00, 0.25135982D+00, 0.25201006D+00,
288      +  0.25215174D+00, 0.25135982D+00, 0.32131134D+01, 0.16691222D+00,
289      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
290      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
291         DATA ((CALC(I,K,7 ),I=1,8), K=1,9) /
292      +  0.57479616D+00, 0.38719266D+00, 0.10278833D+00, 0.10278833D+00,
293      +  0.10278833D+00, 0.10278833D+00, 0.16355532D+01,-0.34694470D-17,
294      +  0.58720174D+00, 0.39563566D+00, 0.10712829D+00, 0.10715295D+00,
295      +  0.10717037D+00, 0.10712829D+00, 0.17455238D+01, 0.78757968D-02,
296      +  0.59856678D+00, 0.40324505D+00, 0.11119185D+00, 0.11124212D+00,
297      +  0.11127289D+00, 0.11119185D+00, 0.18349430D+01, 0.15405101D-01,
298      +  0.60886221D+00, 0.40999594D+00, 0.11493931D+00, 0.11501585D+00,
299      +  0.11505634D+00, 0.11493931D+00, 0.19046215D+01, 0.22549219D-01,
300      +  0.61806541D+00, 0.41587010D+00, 0.11834085D+00, 0.11844404D+00,
301      +  0.11849098D+00, 0.11834085D+00, 0.19554890D+01, 0.29278758D-01,
302      +  0.62616072D+00, 0.42085649D+00, 0.12137335D+00, 0.12150333D+00,
303      +  0.12155379D+00, 0.12137335D+00, 0.19887275D+01, 0.35570892D-01,
304      +  0.63313652D+00, 0.42494831D+00, 0.12401823D+00, 0.12417493D+00,
305      +  0.12422629D+00, 0.12401823D+00, 0.20056144D+01, 0.41407112D-01,
306      +  0.63898674D+00, 0.42814438D+00, 0.12626237D+00, 0.12644549D+00,
307      +  0.12649543D+00, 0.12626237D+00, 0.20075283D+01, 0.46774184D-01,
308      +  0.64371052D+00, 0.43044863D+00, 0.12809747D+00, 0.12830656D+00,
309      +  0.12835302D+00, 0.12809747D+00, 0.19959097D+01, 0.51663585D-01/
310         DATA ((CALC(I,K,7 ),I=1,8), K=10,18) /
311      +  0.64731232D+00, 0.43187014D+00, 0.12951997D+00, 0.12975440D+00,
312      +  0.12979560D+00, 0.12951997D+00, 0.19722294D+01, 0.56071426D-01,
313      +  0.64980187D+00, 0.43242299D+00, 0.13053075D+00, 0.13078975D+00,
314      +  0.13082413D+00, 0.13053075D+00, 0.19379614D+01, 0.59998201D-01,
315      +  0.65119406D+00, 0.43212600D+00, 0.13113485D+00, 0.13141751D+00,
316      +  0.13144375D+00, 0.13113485D+00, 0.18945582D+01, 0.63448515D-01,
317      +  0.65150876D+00, 0.43100251D+00, 0.13134112D+00, 0.13164643D+00,
318      +  0.13166342D+00, 0.13134112D+00, 0.18434285D+01, 0.66430783D-01,
319      +  0.65077057D+00, 0.42908000D+00, 0.13116188D+00, 0.13148871D+00,
320      +  0.13149555D+00, 0.13116188D+00, 0.17859194D+01, 0.68956888D-01,
321      +  0.64900852D+00, 0.42638974D+00, 0.13061250D+00, 0.13095966D+00,
322      +  0.13095562D+00, 0.13061250D+00, 0.17233014D+01, 0.71041811D-01,
323      +  0.64625572D+00, 0.42296638D+00, 0.12971104D+00, 0.13007724D+00,
324      +  0.13006177D+00, 0.12971104D+00, 0.16567571D+01, 0.72703245D-01,
325      +  0.64254896D+00, 0.41884749D+00, 0.12847776D+00, 0.12886168D+00,
326      +  0.12883440D+00, 0.12847776D+00, 0.15873722D+01, 0.73961188D-01,
327      +  0.63792833D+00, 0.41407310D+00, 0.12693476D+00, 0.12733502D+00,
328      +  0.12729569D+00, 0.12693476D+00, 0.15161309D+01, 0.74837535D-01/
329         DATA ((CALC(I,K,7 ),I=1,8), K=19,20) /
330      +  0.63243672D+00, 0.40868521D+00, 0.12510555D+00, 0.12552074D+00,
331      +  0.12546924D+00, 0.12510555D+00, 0.14439130D+01, 0.75355673D-01,
332      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
333      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
334         DATA ((CALC(I,K,8 ),I=1,8), K=1,9) /
335      +  0.66868210D+00, 0.43706345D+00, 0.72925294D-01, 0.72925294D-01,
336      +  0.72925294D-01, 0.72925294D-01, 0.13980776D+01, 0.17347235D-17,
337      +  0.67274888D+00, 0.43768741D+00, 0.74937012D-01, 0.74969093D-01,
338      +  0.74969402D-01, 0.74937012D-01, 0.14300181D+01, 0.52887558D-02,
339      +  0.67556992D+00, 0.43742644D+00, 0.76645438D-01, 0.76708235D-01,
340      +  0.76706403D-01, 0.76645438D-01, 0.14465159D+01, 0.10164601D-01,
341      +  0.67717298D+00, 0.43630833D+00, 0.78045072D-01, 0.78137148D-01,
342      +  0.78131171D-01, 0.78045072D-01, 0.14492734D+01, 0.14625694D-01,
343      +  0.67758301D+00, 0.43435862D+00, 0.79133832D-01, 0.79253695D-01,
344      +  0.79241904D-01, 0.79133832D-01, 0.14397941D+01, 0.18672619D-01,
345      +  0.67683131D+00, 0.43160862D+00, 0.79914994D-01, 0.80061106D-01,
346      +  0.80042142D-01, 0.79914994D-01, 0.14196350D+01, 0.22311181D-01,
347      +  0.67495101D+00, 0.42809125D+00, 0.80394014D-01, 0.80564798D-01,
348      +  0.80537583D-01, 0.80394014D-01, 0.13902979D+01, 0.25549084D-01,
349      +  0.67197853D+00, 0.42384224D+00, 0.80579342D-01, 0.80773192D-01,
350      +  0.80736896D-01, 0.80579342D-01, 0.13532281D+01, 0.28396835D-01,
351      +  0.66795328D+00, 0.41889965D+00, 0.80481868D-01, 0.80697156D-01,
352      +  0.80651177D-01, 0.80481868D-01, 0.13097992D+01, 0.30867233D-01/
353         DATA ((CALC(I,K,8 ),I=1,8), K=10,18) /
354      +  0.66291757D+00, 0.41330371D+00, 0.80114657D-01, 0.80349743D-01,
355      +  0.80293679D-01, 0.80114657D-01, 0.12612964D+01, 0.32975137D-01,
356      +  0.65691643D+00, 0.40709650D+00, 0.79492644D-01, 0.79745886D-01,
357      +  0.79679516D-01, 0.79492644D-01, 0.12089080D+01, 0.34737194D-01,
358      +  0.64999738D+00, 0.40032171D+00, 0.78632316D-01, 0.78902073D-01,
359      +  0.78825336D-01, 0.78632316D-01, 0.11537173D+01, 0.36171551D-01,
360      +  0.64221020D+00, 0.39302424D+00, 0.77551397D-01, 0.77836042D-01,
361      +  0.77749018D-01, 0.77551397D-01, 0.10966991D+01, 0.37297558D-01,
362      +  0.63360659D+00, 0.38524991D+00, 0.76268533D-01, 0.76566460D-01,
363      +  0.76469352D-01, 0.76268533D-01, 0.10387177D+01, 0.38135485D-01,
364      +  0.62423984D+00, 0.37704505D+00, 0.74802994D-01, 0.75112622D-01,
365      +  0.75005743D-01, 0.74802994D-01, 0.98052794D+00, 0.38706224D-01,
366      +  0.61416450D+00, 0.36845616D+00, 0.73174378D-01, 0.73494162D-01,
367      +  0.73377914D-01, 0.73174378D-01, 0.92277789D+00, 0.39031021D-01,
368      +  0.60343597D+00, 0.35952958D+00, 0.71402341D-01, 0.71730776D-01,
369      +  0.71605639D-01, 0.71402341D-01, 0.86601365D+00, 0.39131206D-01,
370      +  0.59211017D+00, 0.35031110D+00, 0.69506346D-01, 0.69841971D-01,
371      +  0.69708492D-01, 0.69506346D-01, 0.81068532D+00, 0.39027956D-01/
372         DATA ((CALC(I,K,8 ),I=1,8), K=19,20) /
373      +  0.58024319D+00, 0.34084568D+00, 0.67505436D-01, 0.67846845D-01,
374      +  0.67705621D-01, 0.67505436D-01, 0.75715400D+00, 0.38742071D-01,
375      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
376      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
377         DATA ((CALC(I,K,9 ),I=1,8), K=1,9) /
378      +  0.73575945D+00, 0.45004095D+00, 0.44221888D-01, 0.44221888D-01,
379      +  0.44221888D-01, 0.44221888D-01, 0.11123387D+01, 0.17347235D-17,
380      +  0.72905091D+00, 0.44200376D+00, 0.44938216D-01, 0.44962187D-01,
381      +  0.44951773D-01, 0.44938216D-01, 0.10927997D+01, 0.32370744D-02,
382      +  0.72122873D+00, 0.43339858D+00, 0.45406888D-01, 0.45452988D-01,
383      +  0.45432116D-01, 0.45406888D-01, 0.10647140D+01, 0.61056719D-02,
384      +  0.71236648D+00, 0.42428172D+00, 0.45642314D-01, 0.45708755D-01,
385      +  0.45677543D-01, 0.45642314D-01, 0.10298317D+01, 0.86263445D-02,
386      +  0.70252667D+00, 0.41470109D+00, 0.45657349D-01, 0.45742385D-01,
387      +  0.45701072D-01, 0.45657349D-01, 0.98955628D+00, 0.10816695D-01,
388      +  0.69177531D+00, 0.40470676D+00, 0.45467489D-01, 0.45569428D-01,
389      +  0.45518350D-01, 0.45467489D-01, 0.94519155D+00, 0.12696783D-01,
390      +  0.68017772D+00, 0.39434784D+00, 0.45088271D-01, 0.45205471D-01,
391      +  0.45145049D-01, 0.45088271D-01, 0.89791119D+00, 0.14286361D-01,
392      +  0.66779950D+00, 0.38367297D+00, 0.44535809D-01, 0.44666684D-01,
393      +  0.44597411D-01, 0.44535809D-01, 0.84875360D+00, 0.15605498D-01,
394      +  0.65470640D+00, 0.37273021D+00, 0.43826532D-01, 0.43969552D-01,
395      +  0.43891975D-01, 0.43826532D-01, 0.79862959D+00, 0.16674360D-01/
396         DATA ((CALC(I,K,9 ),I=1,8), K=10,18) /
397      +  0.64096417D+00, 0.36156675D+00, 0.42976966D-01, 0.43130659D-01,
398      +  0.43045371D-01, 0.42976966D-01, 0.74832328D+00, 0.17513024D-01,
399      +  0.62663835D+00, 0.35022874D+00, 0.42003571D-01, 0.42166530D-01,
400      +  0.42074155D-01, 0.42003571D-01, 0.69849820D+00, 0.18141352D-01,
401      +  0.61179409D+00, 0.33876109D+00, 0.40922567D-01, 0.41093446D-01,
402      +  0.40994632D-01, 0.40922567D-01, 0.64970429D+00, 0.18578832D-01,
403      +  0.59649591D+00, 0.32720723D+00, 0.39749769D-01, 0.39927289D-01,
404      +  0.39822699D-01, 0.39749769D-01, 0.60238607D+00, 0.18844445D-01,
405      +  0.58080739D+00, 0.31560889D+00, 0.38500446D-01, 0.38683394D-01,
406      +  0.38573695D-01, 0.38500446D-01, 0.55689160D+00, 0.18956529D-01,
407      +  0.56479094D+00, 0.30400594D+00, 0.37189188D-01, 0.37376419D-01,
408      +  0.37262280D-01, 0.37189188D-01, 0.51348211D+00, 0.18932666D-01,
409      +  0.54850750D+00, 0.29243618D+00, 0.35829802D-01, 0.36020242D-01,
410      +  0.35902321D-01, 0.35829802D-01, 0.47234183D+00, 0.18789586D-01,
411      +  0.53201633D+00, 0.28093522D+00, 0.34435222D-01, 0.34627863D-01,
412      +  0.34506807D-01, 0.34435222D-01, 0.43358779D+00, 0.18543082D-01,
413      +  0.51537470D+00, 0.26953633D+00, 0.33017438D-01, 0.33211345D-01,
414      +  0.33087782D-01, 0.33017438D-01, 0.39727941D+00, 0.18207948D-01/
415         DATA ((CALC(I,K,9 ),I=1,8), K=19,20) /
416      +  0.49863771D+00, 0.25827032D+00, 0.31587457D-01, 0.31781763D-01,
417      +  0.31656298D-01, 0.31587457D-01, 0.36342766D+00, 0.17797939D-01,
418      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
419      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
420         DATA ((CALC(I,K,10),I=1,8), K=1,9) /
421      +  0.75033233D+00, 0.41928473D+00, 0.26014900D-01, 0.26014900D-01,
422      +  0.26014900D-01, 0.26014900D-01, 0.87281313D+00,-0.16696713D-16,
423      +  0.73440392D+00, 0.40568081D+00, 0.26292098D-01, 0.26309965D-01,
424      +  0.26298466D-01, 0.26292098D-01, 0.83163604D+00, 0.20306844D-02,
425      +  0.71775563D+00, 0.39196078D+00, 0.26384089D-01, 0.26417954D-01,
426      +  0.26395831D-01, 0.26384089D-01, 0.78707432D+00, 0.37712302D-02,
427      +  0.70047819D+00, 0.37817567D+00, 0.26310868D-01, 0.26358988D-01,
428      +  0.26327118D-01, 0.26310868D-01, 0.74047548D+00, 0.52484018D-02,
429      +  0.68264767D+00, 0.36436792D+00, 0.26088831D-01, 0.26149562D-01,
430      +  0.26108817D-01, 0.26088831D-01, 0.69284970D+00, 0.64838907D-02,
431      +  0.66434041D+00, 0.35057926D+00, 0.25735001D-01, 0.25806801D-01,
432      +  0.25758039D-01, 0.25735001D-01, 0.64505024D+00, 0.74997675D-02,
433      +  0.64563037D+00, 0.33684942D+00, 0.25265452D-01, 0.25346873D-01,
434      +  0.25290936D-01, 0.25265452D-01, 0.59779857D+00, 0.83167852D-02,
435      +  0.62658942D+00, 0.32321607D+00, 0.24695598D-01, 0.24785288D-01,
436      +  0.24722992D-01, 0.24695598D-01, 0.55167852D+00, 0.89547331D-02,
437      +  0.60728743D+00, 0.30971488D+00, 0.24040168D-01, 0.24136866D-01,
438      +  0.24069001D-01, 0.24040168D-01, 0.50715463D+00, 0.94324539D-02/
439         DATA ((CALC(I,K,10),I=1,8), K=10,18) /
440      +  0.58779207D+00, 0.29637937D+00, 0.23313100D-01, 0.23415631D-01,
441      +  0.23342957D-01, 0.23313100D-01, 0.46458231D+00, 0.97677659D-02,
442      +  0.56816868D+00, 0.28324079D+00, 0.22527502D-01, 0.22634777D-01,
443      +  0.22558022D-01, 0.22527502D-01, 0.42422017D+00, 0.99774566D-02,
444      +  0.54848012D+00, 0.27032807D+00, 0.21695605D-01, 0.21806618D-01,
445      +  0.21726472D-01, 0.21695605D-01, 0.38624235D+00, 0.10077255D-01,
446      +  0.52878659D+00, 0.25766771D+00, 0.20828727D-01, 0.20942551D-01,
447      +  0.20859669D-01, 0.20828727D-01, 0.35075014D+00, 0.10081812D-01,
448      +  0.50914547D+00, 0.24528374D+00, 0.19937246D-01, 0.20053032D-01,
449      +  0.19968028D-01, 0.19937246D-01, 0.31778344D+00, 0.10004694D-01,
450      +  0.48961113D+00, 0.23319766D+00, 0.19030589D-01, 0.19147561D-01,
451      +  0.19061011D-01, 0.19030589D-01, 0.28733133D+00, 0.98583776D-02,
452      +  0.47023480D+00, 0.22142843D+00, 0.18117230D-01, 0.18234685D-01,
453      +  0.18147123D-01, 0.18117230D-01, 0.25934185D+00, 0.96542623D-02,
454      +  0.45106446D+00, 0.20999249D+00, 0.17204704D-01, 0.17322007D-01,
455      +  0.17233927D-01, 0.17204704D-01, 0.23373099D+00, 0.94026882D-02,
456      +  0.43214467D+00, 0.19890378D+00, 0.16299628D-01, 0.16416209D-01,
457      +  0.16328067D-01, 0.16299628D-01, 0.21039060D+00, 0.91129667D-02/
458         DATA ((CALC(I,K,10),I=1,8), K=19,20) /
459      +  0.41351656D+00, 0.18817378D+00, 0.15407739D-01, 0.15523092D-01,
460      +  0.15435299D-01, 0.15407739D-01, 0.18919538D+00, 0.87934189D-02,
461      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
462      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
463         DATA ((CALC(I,K,11),I=1,8), K=1,9) /
464      +  0.72612470D+00, 0.36397949D+00, 0.14788752D-01, 0.14788752D-01,
465      +  0.14788752D-01, 0.14788752D-01, 0.67422733D+00,-0.34694470D-17,
466      +  0.70301699D+00, 0.34773427D+00, 0.14929544D-01, 0.14942274D-01,
467      +  0.14932700D-01, 0.14929544D-01, 0.62610341D+00, 0.12768068D-02,
468      +  0.67971451D+00, 0.33178907D+00, 0.14935458D-01, 0.14959291D-01,
469      +  0.14941229D-01, 0.14935458D-01, 0.57809221D+00, 0.23388262D-02,
470      +  0.65630177D+00, 0.31617176D+00, 0.14826093D-01, 0.14859555D-01,
471      +  0.14834012D-01, 0.14826093D-01, 0.53104444D+00, 0.32117599D-02,
472      +  0.63284895D+00, 0.30090515D+00, 0.14616981D-01, 0.14658717D-01,
473      +  0.14626640D-01, 0.14616981D-01, 0.48554748D+00, 0.39159194D-02,
474      +  0.60942378D+00, 0.28600989D+00, 0.14323061D-01, 0.14371833D-01,
475      +  0.14334102D-01, 0.14323061D-01, 0.44202991D+00, 0.44707510D-02,
476      +  0.58609074D+00, 0.27150488D+00, 0.13958017D-01, 0.14012694D-01,
477      +  0.13970128D-01, 0.13958017D-01, 0.40081736D+00, 0.48941032D-02,
478      +  0.56291077D+00, 0.25740685D+00, 0.13534380D-01, 0.13593930D-01,
479      +  0.13547291D-01, 0.13534380D-01, 0.36212733D+00, 0.52023720D-02,
480      +  0.53994148D+00, 0.24373049D+00, 0.13063647D-01, 0.13127130D-01,
481      +  0.13077122D-01, 0.13063647D-01, 0.32608954D+00, 0.54106573D-02/
482         DATA ((CALC(I,K,11),I=1,8), K=10,18) /
483      +  0.51723696D+00, 0.23048836D+00, 0.12556258D-01, 0.12622825D-01,
484      +  0.12570094D-01, 0.12556258D-01, 0.29275945D+00, 0.55327699D-02,
485      +  0.49484775D+00, 0.21769088D+00, 0.12021642D-01, 0.12090523D-01,
486      +  0.12035665D-01, 0.12021642D-01, 0.26213118D+00, 0.55812975D-02,
487      +  0.47282077D+00, 0.20534634D+00, 0.11468242D-01, 0.11538746D-01,
488      +  0.11482305D-01, 0.11468242D-01, 0.23414975D+00, 0.55676557D-02,
489      +  0.45119921D+00, 0.19346088D+00, 0.10903556D-01, 0.10975062D-01,
490      +  0.10917533D-01, 0.10903556D-01, 0.20872192D+00, 0.55021381D-02,
491      +  0.43002249D+00, 0.18203854D+00, 0.10334171D-01, 0.10406127D-01,
492      +  0.10347957D-01, 0.10334171D-01, 0.18572580D+00, 0.53939699D-02,
493      +  0.40932621D+00, 0.17108134D+00, 0.97658089D-02, 0.98377266D-02,
494      +  0.97793177D-02, 0.97658089D-02, 0.16501926D+00, 0.52513640D-02,
495      +  0.38914211D+00, 0.16058928D+00, 0.92033789D-02, 0.92748267D-02,
496      +  0.92165392D-02, 0.92033789D-02, 0.14644725D+00, 0.50815808D-02,
497      +  0.36949815D+00, 0.15056052D+00, 0.86510266D-02, 0.87216286D-02,
498      +  0.86637817D-02, 0.86510266D-02, 0.12984774D+00, 0.48909896D-02,
499      +  0.35041843D+00, 0.14099142D+00, 0.81121930D-02, 0.81816238D-02,
500      +  0.81244986D-02, 0.81121930D-02, 0.11505677D+00, 0.46851323D-02/
501         DATA ((CALC(I,K,11),I=1,8), K=19,20) /
502      +  0.33192335D+00, 0.13187666D+00, 0.75896736D-02, 0.76576547D-02,
503      +  0.76014961D-02, 0.75896736D-02, 0.10191238D+00, 0.44687892D-02,
504      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
505      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
506         DATA ((CALC(I,K,12),I=1,8), K=1,9) /
507      +  0.67334252D+00, 0.29833441D+00, 0.80855825D-02, 0.80855825D-02,
508      +  0.80855825D-02, 0.80855825D-02, 0.51162697D+00,-0.11817804D-16,
509      +  0.64536318D+00, 0.28180015D+00, 0.81879384D-02, 0.81964326D-02,
510      +  0.81895119D-02, 0.81879384D-02, 0.46430924D+00, 0.79455875D-03,
511      +  0.61774910D+00, 0.26586805D+00, 0.81936895D-02, 0.82094634D-02,
512      +  0.81965473D-02, 0.81936895D-02, 0.41926385D+00, 0.14371427D-02,
513      +  0.59056001D+00, 0.25053723D+00, 0.81196347D-02, 0.81416057D-02,
514      +  0.81235309D-02, 0.81196347D-02, 0.37688432D+00, 0.19493803D-02,
515      +  0.56384502D+00, 0.23580681D+00, 0.79788388D-02, 0.80060248D-02,
516      +  0.79835596D-02, 0.79788388D-02, 0.33739524D+00, 0.23481052D-02,
517      +  0.53764896D+00, 0.22167368D+00, 0.77831646D-02, 0.78146805D-02,
518      +  0.77885249D-02, 0.77831646D-02, 0.30089054D+00, 0.26486947D-02,
519      +  0.51201355D+00, 0.20813409D+00, 0.75432692D-02, 0.75783182D-02,
520      +  0.75491096D-02, 0.75432692D-02, 0.26739948D+00, 0.28650364D-02,
521      +  0.48697655D+00, 0.19518293D+00, 0.72685826D-02, 0.73064492D-02,
522      +  0.72747660D-02, 0.72685826D-02, 0.23688280D+00, 0.30095327D-02,
523      +  0.46257205D+00, 0.18281391D+00, 0.69674848D-02, 0.70075288D-02,
524      +  0.69738943D-02, 0.69674848D-02, 0.20924982D+00, 0.30933111D-02/
525         DATA ((CALC(I,K,12),I=1,8), K=10,18) /
526      +  0.43883036D+00, 0.17101948D+00, 0.66473444D-02, 0.66889947D-02,
527      +  0.66538804D-02, 0.66473444D-02, 0.18437152D+00, 0.31262846D-02,
528      +  0.41577800D+00, 0.15979084D+00, 0.63145934D-02, 0.63573430D-02,
529      +  0.63211718D-02, 0.63145934D-02, 0.16209093D+00, 0.31172474D-02,
530      +  0.39343772D+00, 0.14911797D+00, 0.59747963D-02, 0.60181973D-02,
531      +  0.59813468D-02, 0.59747963D-02, 0.14223264D+00, 0.30739602D-02,
532      +  0.37182850D+00, 0.13898967D+00, 0.56327140D-02, 0.56763733D-02,
533      +  0.56391785D-02, 0.56327140D-02, 0.12461070D+00, 0.30032288D-02,
534      +  0.35096559D+00, 0.12939361D+00, 0.52923700D-02, 0.53359442D-02,
535      +  0.52987006D-02, 0.52923700D-02, 0.10903509D+00, 0.29109819D-02,
536      +  0.33086057D+00, 0.12031641D+00, 0.49571134D-02, 0.50003055D-02,
537      +  0.49632718D-02, 0.49571134D-02, 0.95316920D-01, 0.28023456D-02,
538      +  0.31152147D+00, 0.11174375D+00, 0.46296827D-02, 0.46722381D-02,
539      +  0.46356387D-02, 0.46296827D-02, 0.83272498D-01, 0.26817158D-02,
540      +  0.29295284D+00, 0.10366045D+00, 0.43122686D-02, 0.43539711D-02,
541      +  0.43179989D-02, 0.43122686D-02, 0.72726380D-01, 0.25528291D-02,
542      +  0.27515593D+00, 0.96050577D-01, 0.40065739D-02, 0.40472427D-02,
543      +  0.40120614D-02, 0.40065739D-02, 0.63513533D-01, 0.24188299D-02/
544         DATA ((CALC(I,K,12),I=1,8), K=19,20) /
545      +  0.25812880D+00, 0.88897572D-01, 0.37138719D-02, 0.37533580D-02,
546      +  0.37191047D-02, 0.37138719D-02, 0.55480708D-01, 0.22823356D-02,
547      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
548      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
549         DATA ((CALC(I,K,13),I=1,8), K=1,9) /
550      +  0.60077312D+00, 0.23221598D+00, 0.42271747D-02, 0.42271747D-02,
551      +  0.42271747D-02, 0.42271747D-02, 0.38037736D+00, 0.86736174D-18,
552      +  0.57027955D+00, 0.21703986D+00, 0.43144662D-02, 0.43201956D-02,
553      +  0.43152794D-02, 0.43144662D-02, 0.33787940D+00, 0.48472802D-03,
554      +  0.54066854D+00, 0.20263136D+00, 0.43340412D-02, 0.43446000D-02,
555      +  0.43355085D-02, 0.43340412D-02, 0.29879859D+00, 0.86632144D-03,
556      +  0.51196237D+00, 0.18896360D+00, 0.42991523D-02, 0.43137484D-02,
557      +  0.43011396D-02, 0.42991523D-02, 0.26316287D+00, 0.11614835D-02,
558      +  0.48417903D+00, 0.17601483D+00, 0.42200012D-02, 0.42379246D-02,
559      +  0.42223933D-02, 0.42200012D-02, 0.23092273D+00, 0.13830959D-02,
560      +  0.45733163D+00, 0.16376189D+00, 0.41053988D-02, 0.41260179D-02,
561      +  0.41080969D-02, 0.41053988D-02, 0.20193947D+00, 0.15424386D-02,
562      +  0.43143123D+00, 0.15218235D+00, 0.39631584D-02, 0.39859122D-02,
563      +  0.39660784D-02, 0.39631584D-02, 0.17604658D+00, 0.16495813D-02,
564      +  0.40648558D+00, 0.14125372D+00, 0.38000062D-02, 0.38243982D-02,
565      +  0.38030767D-02, 0.38000062D-02, 0.15304734D+00, 0.17133078D-02,
566      +  0.38249945D+00, 0.13095353D+00, 0.36217606D-02, 0.36473530D-02,
567      +  0.36249215D-02, 0.36217606D-02, 0.13272620D+00, 0.17413214D-02/
568         DATA ((CALC(I,K,13),I=1,8), K=10,18) /
569      +  0.35947453D+00, 0.12125927D+00, 0.34333998D-02, 0.34598086D-02,
570      +  0.34366006D-02, 0.34333998D-02, 0.11485913D+00, 0.17403296D-02,
571      +  0.33740948D+00, 0.11214834D+00, 0.32391427D-02, 0.32660329D-02,
572      +  0.32423417D-02, 0.32391427D-02, 0.99220244D-01, 0.17161406D-02,
573      +  0.31629996D+00, 0.10359808D+00, 0.30425279D-02, 0.30696090D-02,
574      +  0.30456907D-02, 0.30425279D-02, 0.85587700D-01, 0.16737553D-02,
575      +  0.29613873D+00, 0.95585734D-01, 0.28464824D-02, 0.28735045D-02,
576      +  0.28495812D-02, 0.28464824D-02, 0.73748068D-01, 0.16174479D-02,
577      +  0.27691572D+00, 0.88088565D-01, 0.26533892D-02, 0.26801392D-02,
578      +  0.26564018D-02, 0.26533892D-02, 0.63499475D-01, 0.15508429D-02,
579      +  0.25861818D+00, 0.81083843D-01, 0.24651504D-02, 0.24914485D-02,
580      +  0.24680595D-02, 0.24651504D-02, 0.54653783D-01, 0.14769869D-02,
581      +  0.24123083D+00, 0.74548927D-01, 0.22832461D-02, 0.23089427D-02,
582      +  0.22860388D-02, 0.22832461D-02, 0.47037889D-01, 0.13984152D-02,
583      +  0.22473600D+00, 0.68461330D-01, 0.21087905D-02, 0.21337628D-02,
584      +  0.21114573D-02, 0.21087905D-02, 0.40494306D-01, 0.13172148D-02,
585      +  0.20911389D+00, 0.62798787D-01, 0.19425826D-02, 0.19667322D-02,
586      +  0.19451171D-02, 0.19425826D-02, 0.34881210D-01, 0.12350809D-02/
587         DATA ((CALC(I,K,13),I=1,8), K=19,20) /
588      +  0.19434270D+00, 0.57539325D-01, 0.17851541D-02, 0.18084041D-02,
589      +  0.17875526D-02, 0.17851541D-02, 0.30072058D-01, 0.11533708D-02,
590      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
591      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
592         DATA ((CALC(I,K,14),I=1,8), K=1,9) /
593      +  0.51632537D+00, 0.17190006D+00, 0.20981348D-02, 0.20981348D-02,
594      +  0.20981348D-02, 0.20981348D-02, 0.27616356D+00, 0.11926224D-17,
595      +  0.48552628D+00, 0.15904827D+00, 0.21710311D-02, 0.21749231D-02,
596      +  0.21714614D-02, 0.21710311D-02, 0.24030554D+00, 0.28736293D-03,
597      +  0.45602587D+00, 0.14700374D+00, 0.21979326D-02, 0.22050439D-02,
598      +  0.21987031D-02, 0.21979326D-02, 0.20827361D+00, 0.50772576D-03,
599      +  0.42780301D+00, 0.13572121D+00, 0.21886054D-02, 0.21983531D-02,
600      +  0.21896412D-02, 0.21886054D-02, 0.17983633D+00, 0.67310759D-03,
601      +  0.40084025D+00, 0.12516450D+00, 0.21505429D-02, 0.21624119D-02,
602      +  0.21517801D-02, 0.21505429D-02, 0.15476148D+00, 0.79274167D-03,
603      +  0.37511595D+00, 0.11529731D+00, 0.20899179D-02, 0.21034569D-02,
604      +  0.20913030D-02, 0.20899179D-02, 0.13277141D+00, 0.87438403D-03,
605      +  0.35060822D+00, 0.10608530D+00, 0.20121614D-02, 0.20269756D-02,
606      +  0.20136489D-02, 0.20121614D-02, 0.11359060D+00, 0.92490325D-03,
607      +  0.32729347D+00, 0.97495376D-01, 0.19218442D-02, 0.19375899D-02,
608      +  0.19233964D-02, 0.19218442D-02, 0.96944662D-01, 0.95016958D-03,
609      +  0.30514661D+00, 0.89495614D-01, 0.18228262D-02, 0.18392058D-02,
610      +  0.18244118D-02, 0.18228262D-02, 0.82565447D-01, 0.95522237D-03/
611         DATA ((CALC(I,K,14),I=1,8), K=10,18) /
612      +  0.28414102D+00, 0.82055168D-01, 0.17183336D-02, 0.17350908D-02,
613      +  0.17199267D-02, 0.17183336D-02, 0.70197637D-01, 0.94435964D-03,
614      +  0.26424851D+00, 0.75144175D-01, 0.16110292D-02, 0.16279447D-02,
615      +  0.16126089D-02, 0.16110292D-02, 0.59601612D-01, 0.92121971D-03,
616      +  0.24543941D+00, 0.68733711D-01, 0.15030839D-02, 0.15199719D-02,
617      +  0.15046335D-02, 0.15030839D-02, 0.50555726D-01, 0.88886214D-03,
618      +  0.22768264D+00, 0.62795758D-01, 0.13962381D-02, 0.14129425D-02,
619      +  0.13977442D-02, 0.13962381D-02, 0.42857567D-01, 0.84983733D-03,
620      +  0.21094579D+00, 0.57303184D-01, 0.12918587D-02, 0.13082502D-02,
621      +  0.12933112D-02, 0.12918587D-02, 0.36324336D-01, 0.80625106D-03,
622      +  0.19519530D+00, 0.52229746D-01, 0.11909920D-02, 0.12069649D-02,
623      +  0.11923834D-02, 0.11909920D-02, 0.30792628D-01, 0.75982316D-03,
624      +  0.18039655D+00, 0.47550088D-01, 0.10944108D-02, 0.11098805D-02,
625      +  0.10957358D-02, 0.10944108D-02, 0.26117732D-01, 0.71194080D-03,
626      +  0.16651410D+00, 0.43239757D-01, 0.10026578D-02, 0.10175579D-02,
627      +  0.10039127D-02, 0.10026578D-02, 0.22172590D-01, 0.66370657D-03,
628      +  0.15351180D+00, 0.39275213D-01, 0.91608379D-03, 0.93036438D-03,
629      +  0.91726670D-03, 0.91608379D-03, 0.18846509D-01, 0.61598147D-03/
630         DATA ((CALC(I,K,14),I=1,8), K=19,20) /
631      +  0.14135303D+00, 0.35633841D-01, 0.83488202D-03, 0.84850733D-03,
632      +  0.83599225D-03, 0.83488202D-03, 0.16043733D-01, 0.56942327D-03,
633      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
634      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
635         DATA ((CALC(I,K,15),I=1,8), K=1,9) /
636      +  0.42705580D+00, 0.12082176D+00, 0.97978438D-03, 0.97978438D-03,
637      +  0.97978438D-03, 0.97978438D-03, 0.19498989D+00,-0.25072175D-18,
638      +  0.39784903D+00, 0.11068616D+00, 0.10356646D-02, 0.10382274D-02,
639      +  0.10358873D-02, 0.10356646D-02, 0.16625863D+00, 0.16399453D-03,
640      +  0.37021828D+00, 0.10130049D+00, 0.10616397D-02, 0.10662818D-02,
641      +  0.10620353D-02, 0.10616397D-02, 0.14125776D+00, 0.28655433D-03,
642      +  0.34409999D+00, 0.92610460D-01, 0.10644378D-02, 0.10707467D-02,
643      +  0.10649653D-02, 0.10644378D-02, 0.11959872D+00, 0.37575575D-03,
644      +  0.31944260D+00, 0.84573060D-01, 0.10492410D-02, 0.10568576D-02,
645      +  0.10498663D-02, 0.10492410D-02, 0.10094846D+00, 0.43781550D-03,
646      +  0.29619219D+00, 0.77146471D-01, 0.10201214D-02, 0.10287351D-02,
647      +  0.10208158D-02, 0.10201214D-02, 0.84967494D-01, 0.47773454D-03,
648      +  0.27429679D+00, 0.70291650D-01, 0.98064902D-03, 0.98999276D-03,
649      +  0.98138896D-03, 0.98064902D-03, 0.71339424D-01, 0.49992158D-03,
650      +  0.25370497D+00, 0.63971782D-01, 0.93377413D-03, 0.94361922D-03,
651      +  0.93454012D-03, 0.93377413D-03, 0.59770974D-01, 0.50807768D-03,
652      +  0.23436584D+00, 0.58152064D-01, 0.88193337D-03, 0.89208528D-03,
653      +  0.88270955D-03, 0.88193337D-03, 0.49991630D-01, 0.50531460D-03/
654         DATA ((CALC(I,K,15),I=1,8), K=10,18) /
655      +  0.21622892D+00, 0.52799564D-01, 0.82712198D-03, 0.83741647D-03,
656      +  0.82789559D-03, 0.82712198D-03, 0.41756407D-01, 0.49423547D-03,
657      +  0.19924413D+00, 0.47883066D-01, 0.77094786D-03, 0.78124752D-03,
658      +  0.77170876D-03, 0.77094786D-03, 0.34845585D-01, 0.47699601D-03,
659      +  0.18336173D+00, 0.43372975D-01, 0.71468791D-03, 0.72487906D-03,
660      +  0.71542819D-03, 0.71468791D-03, 0.29064050D-01, 0.45536721D-03,
661      +  0.16853240D+00, 0.39241220D-01, 0.65933595D-03, 0.66932576D-03,
662      +  0.66004960D-03, 0.65933595D-03, 0.24240208D-01, 0.43078871D-03,
663      +  0.15470724D+00, 0.35461188D-01, 0.60564592D-03, 0.61535992D-03,
664      +  0.60632853D-03, 0.60564592D-03, 0.20224430D-01, 0.40441676D-03,
665      +  0.14183788D+00, 0.32007664D-01, 0.55417086D-03, 0.56355056D-03,
666      +  0.55481932D-03, 0.55417086D-03, 0.16887306D-01, 0.37716744D-03,
667      +  0.12987657D+00, 0.28856782D-01, 0.50529713D-03, 0.51429798D-03,
668      +  0.50590946D-03, 0.50529713D-03, 0.14117757D-01, 0.34975466D-03,
669      +  0.11877628D+00, 0.25985981D-01, 0.45927458D-03, 0.46786404D-03,
670      +  0.45984969D-03, 0.45927458D-03, 0.11821079D-01, 0.32272356D-03,
671      +  0.10849081D+00, 0.23373964D-01, 0.41624290D-03, 0.42439878D-03,
672      +  0.41678046D-03, 0.41624290D-03, 0.99170713D-02, 0.29647999D-03/
673         DATA ((CALC(I,K,15),I=1,8), K=19,20) /
674      +  0.98974901D-01, 0.21000664D-01, 0.37625434D-03, 0.38396327D-03,
675      +  0.37675462D-03, 0.37625434D-03, 0.83382513D-02, 0.27131583D-03,
676      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
677      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
678         DATA ((CALC(I,K,16),I=1,8), K=1,9) /
679      +  0.33903118D+00, 0.80276474D-01, 0.42548719D-03, 0.42548719D-03,
680      +  0.42548719D-03, 0.42548719D-03, 0.13318071D+00,-0.17618285D-18,
681      +  0.31288724D+00, 0.72817098D-01, 0.46413041D-03, 0.46575439D-03,
682      +  0.46424142D-03, 0.46413041D-03, 0.11125492D+00, 0.89075813D-04,
683      +  0.28844204D+00, 0.65988931D-01, 0.48442470D-03, 0.48733857D-03,
684      +  0.48462021D-03, 0.48442470D-03, 0.92647001D-01, 0.15398437D-03,
685      +  0.26559663D+00, 0.59737298D-01, 0.49067020D-03, 0.49459379D-03,
686      +  0.49092880D-03, 0.49067020D-03, 0.76899569D-01, 0.19976257D-03,
687      +  0.24427155D+00, 0.54019328D-01, 0.48622142D-03, 0.49091428D-03,
688      +  0.48652541D-03, 0.48622142D-03, 0.63645454D-01, 0.23032621D-03,
689      +  0.22438755D+00, 0.48794267D-01, 0.47358551D-03, 0.47884331D-03,
690      +  0.47392034D-03, 0.47358551D-03, 0.52541413D-01, 0.24867979D-03,
691      +  0.20586960D+00, 0.44024481D-01, 0.45496396D-03, 0.46061397D-03,
692      +  0.45531776D-03, 0.45496396D-03, 0.43278927D-01, 0.25747077D-03,
693      +  0.18864559D+00, 0.39675068D-01, 0.43215079D-03, 0.43804785D-03,
694      +  0.43251400D-03, 0.43215079D-03, 0.35585185D-01, 0.25888872D-03,
695      +  0.17264613D+00, 0.35713543D-01, 0.40659624D-03, 0.41261939D-03,
696      +  0.40696119D-03, 0.40659624D-03, 0.29218655D-01, 0.25473615D-03/
697         DATA ((CALC(I,K,16),I=1,8), K=10,18) /
698      +  0.15780432D+00, 0.32109654D-01, 0.37946663D-03, 0.38551605D-03,
699      +  0.37982732D-03, 0.37946663D-03, 0.23968575D-01, 0.24649356D-03,
700      +  0.14405560D+00, 0.28835199D-01, 0.35168151D-03, 0.35767579D-03,
701      +  0.35203327D-03, 0.35168151D-03, 0.19652525D-01, 0.23536084D-03,
702      +  0.13133761D+00, 0.25863873D-01, 0.32395319D-03, 0.32982692D-03,
703      +  0.32429252D-03, 0.32395319D-03, 0.16113788D-01, 0.22230062D-03,
704      +  0.11959016D+00, 0.23171139D-01, 0.29682056D-03, 0.30252215D-03,
705      +  0.29714489D-03, 0.29682056D-03, 0.13218820D-01, 0.20807534D-03,
706      +  0.10875513D+00, 0.20734118D-01, 0.27067861D-03, 0.27616838D-03,
707      +  0.27098617D-03, 0.27067861D-03, 0.10854644D-01, 0.19327964D-03,
708      +  0.98776501D-01, 0.18531488D-01, 0.24580455D-03, 0.25105304D-03,
709      +  0.24609422D-03, 0.24580455D-03, 0.89263130D-02, 0.17836902D-03,
710      +  0.89600324D-01, 0.16543398D-01, 0.22238049D-03, 0.22736690D-03,
711      +  0.22265165D-03, 0.22238049D-03, 0.73545566D-02, 0.16368469D-03,
712      +  0.81174729D-01, 0.14751386D-01, 0.20051276D-03, 0.20522362D-03,
713      +  0.20076523D-03, 0.20051276D-03, 0.60736319D-02, 0.14947494D-03,
714      +  0.73449942D-01, 0.13138306D-01, 0.18024833D-03, 0.18467630D-03,
715      +  0.18048226D-03, 0.18024833D-03, 0.50292731D-02, 0.13591328D-03/
716         DATA ((CALC(I,K,16),I=1,8), K=19,20) /
717      +  0.66378294D-01, 0.11688258D-01, 0.16158863D-03, 0.16573145D-03,
718      +  0.16180444D-03, 0.16158863D-03, 0.41768921D-02, 0.12311380D-03,
719      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
720      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
721         DATA ((CALC(I,K,17),I=1,8), K=1,9) /
722      +  0.25716544D+00, 0.50047417D-01, 0.16921711D-03, 0.16921711D-03,
723      +  0.16921711D-03, 0.16921711D-03, 0.87379852D-01, 0.20803129D-17,
724      +  0.23506196D+00, 0.44940973D-01, 0.19311456D-03, 0.19408476D-03,
725      +  0.19316694D-03, 0.19311456D-03, 0.71473141D-01, 0.45371071D-04,
726      +  0.21463003D+00, 0.40319879D-01, 0.20657405D-03, 0.20829675D-03,
727      +  0.20666547D-03, 0.20657405D-03, 0.58303762D-01, 0.77635734D-04,
728      +  0.19574688D+00, 0.36135696D-01, 0.21210461D-03, 0.21440070D-03,
729      +  0.21222447D-03, 0.21210461D-03, 0.47414988D-01, 0.99667201D-04,
730      +  0.17831486D+00, 0.32350818D-01, 0.21171439D-03, 0.21443283D-03,
731      +  0.21185408D-03, 0.21171439D-03, 0.38455136D-01, 0.11374899D-03,
732      +  0.16223931D+00, 0.28930202D-01, 0.20682295D-03, 0.20983764D-03,
733      +  0.20697547D-03, 0.20682295D-03, 0.31116267D-01, 0.12154118D-03,
734      +  0.14743170D+00, 0.25841808D-01, 0.19868652D-03, 0.20189297D-03,
735      +  0.19884629D-03, 0.19868652D-03, 0.25128836D-01, 0.12451427D-03,
736      +  0.13380862D+00, 0.23056348D-01, 0.18831994D-03, 0.19163212D-03,
737      +  0.18848253D-03, 0.18831994D-03, 0.20263181D-01, 0.12387123D-03,
738      +  0.12129123D+00, 0.20546938D-01, 0.17652595D-03, 0.17987386D-03,
739      +  0.17668791D-03, 0.17652595D-03, 0.16322892D-01, 0.12058072D-03/
740         DATA ((CALC(I,K,17),I=1,8), K=10,18) /
741      +  0.10980498D+00, 0.18288894D-01, 0.16394017D-03, 0.16726752D-03,
742      +  0.16409884D-03, 0.16394017D-03, 0.13141904D-01, 0.11542487D-03,
743      +  0.99279307D-01, 0.16259533D-01, 0.15105373D-03, 0.15431597D-03,
744      +  0.15120712D-03, 0.15105373D-03, 0.10580912D-01, 0.10902423D-03,
745      +  0.89647348D-01, 0.14438001D-01, 0.13823793D-03, 0.14140052D-03,
746      +  0.13838460D-03, 0.13823793D-03, 0.85237466D-02, 0.10186447D-03,
747      +  0.80845790D-01, 0.12805132D-01, 0.12576592D-03, 0.12880284D-03,
748      +  0.12590488D-03, 0.12576592D-03, 0.68742070D-02, 0.94319829D-04,
749      +  0.72814689D-01, 0.11343312D-01, 0.11383101D-03, 0.11672337D-03,
750      +  0.11396162D-03, 0.11383101D-03, 0.55531346D-02, 0.86672930D-04,
751      +  0.65497339D-01, 0.10036368D-01, 0.10256263D-03, 0.10529755D-03,
752      +  0.10268455D-03, 0.10256263D-03, 0.44958340D-02, 0.79132120D-04,
753      +  0.58840152D-01, 0.88694607D-02, 0.92039859D-04, 0.94609435D-04,
754      +  0.92152981D-04, 0.92039859D-04, 0.36496885D-02, 0.71846147D-04,
755      +  0.52792555D-01, 0.78289913D-02, 0.82302735D-04, 0.84703144D-04,
756      +  0.82407123D-04, 0.82302735D-04, 0.29720846D-02, 0.64916525D-04,
757      +  0.47306895D-01, 0.69025121D-02, 0.73361805D-04, 0.75592543D-04,
758      +  0.73457663D-04, 0.73361805D-04, 0.24287797D-02, 0.58407999D-04/
759         DATA ((CALC(I,K,17),I=1,8), K=19,20) /
760      +  0.42338349D-01, 0.60786459D-02, 0.65205807D-04, 0.67269026D-04,
761      +  0.65293443D-04, 0.65205807D-04, 0.19923875D-02, 0.52357046D-04,
762      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
763      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
764         DATA ((CALC(I,K,18),I=1,8), K=1,9) /
765      +  0.18509435D+00, 0.28948732D-01, 0.60365385D-04, 0.60365385D-04,
766      +  0.60365385D-04, 0.60365385D-04, 0.54550804D-01,-0.99102855D-19,
767      +  0.16750202D+00, 0.25723871D-01, 0.73400658D-04, 0.73935524D-04,
768      +  0.73423603D-04, 0.73400658D-04, 0.43645194D-01, 0.21221013D-04,
769      +  0.15142629D+00, 0.22839415D-01, 0.81065112D-04, 0.82004390D-04,
770      +  0.81104814D-04, 0.81065112D-04, 0.34842363D-01, 0.35978415D-04,
771      +  0.13673471D+00, 0.20257027D-01, 0.84656333D-04, 0.85894569D-04,
772      +  0.84707952D-04, 0.84656333D-04, 0.27735150D-01, 0.45730076D-04,
773      +  0.12332287D+00, 0.17947249D-01, 0.85264971D-04, 0.86714700D-04,
774      +  0.85324626D-04, 0.85264971D-04, 0.22019925D-01, 0.51686866D-04,
775      +  0.11109181D+00, 0.15883149D-01, 0.83612150D-04, 0.85201675D-04,
776      +  0.83676741D-04, 0.83612150D-04, 0.17446055D-01, 0.54673316D-04,
777      +  0.99949989D-01, 0.14040339D-01, 0.80347140D-04, 0.82018279D-04,
778      +  0.80414230D-04, 0.80347140D-04, 0.13798816D-01, 0.55430154D-04,
779      +  0.89812509D-01, 0.12396828D-01, 0.75994424D-04, 0.77700363D-04,
780      +  0.76062118D-04, 0.75994424D-04, 0.10901131D-01, 0.54561355D-04,
781      +  0.80600357D-01, 0.10932696D-01, 0.70959398D-04, 0.72663041D-04,
782      +  0.71026248D-04, 0.70959398D-04, 0.86064001D-02, 0.52541686D-04/
783         DATA ((CALC(I,K,18),I=1,8), K=10,18) /
784      +  0.72239983D-01, 0.96298964D-02, 0.65558694D-04, 0.67231142D-04,
785      +  0.65623620D-04, 0.65558694D-04, 0.67942217D-02, 0.49748407D-04,
786      +  0.64662859D-01, 0.84720691D-02, 0.60032677D-04, 0.61651919D-04,
787      +  0.60094895D-04, 0.60032677D-04, 0.53665277D-02, 0.46474873D-04,
788      +  0.57805120D-01, 0.74443801D-02, 0.54558688D-04, 0.56108474D-04,
789      +  0.54617657D-04, 0.54558688D-04, 0.42438055D-02, 0.42944754D-04,
790      +  0.51607257D-01, 0.65333812D-02, 0.49263551D-04, 0.50732414D-04,
791      +  0.49318922D-04, 0.49263551D-04, 0.33620652D-02, 0.39325376D-04,
792      +  0.46013841D-01, 0.57268858D-02, 0.44233750D-04, 0.45614158D-04,
793      +  0.44285328D-04, 0.44233750D-04, 0.26700731D-02, 0.35738622D-04,
794      +  0.40973286D-01, 0.50138581D-02, 0.39524115D-04, 0.40811745D-04,
795      +  0.39571827D-04, 0.39524115D-04, 0.21269517D-02, 0.32270251D-04,
796      +  0.36437630D-01, 0.43843135D-02, 0.35165232D-04, 0.36358345D-04,
797      +  0.35209098D-04, 0.35165232D-04, 0.17003738D-02, 0.28977854D-04,
798      +  0.32362334D-01, 0.38292285D-02, 0.31169503D-04, 0.32268408D-04,
799      +  0.31209611D-04, 0.31169503D-04, 0.13649232D-02, 0.25897389D-04,
800      +  0.28706102D-01, 0.33404582D-02, 0.27535873D-04, 0.28542471D-04,
801      +  0.27572364D-04, 0.27535873D-04, 0.11005374D-02, 0.23048338D-04/
802         DATA ((CALC(I,K,18),I=1,8), K=19,20) /
803      +  0.25430700D-01, 0.29106608D-02, 0.24253781D-04, 0.25171186D-04,
804      +  0.24286832D-04, 0.24253781D-04, 0.89149947D-03, 0.20438022D-04,
805      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
806      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
807         DATA ((CALC(I,K,19),I=1,8), K=1,9) /
808      +  0.12511596D+00, 0.15277434D-01, 0.18761940D-04, 0.18761940D-04,
809      +  0.18761940D-04, 0.18761940D-04, 0.31976696D-01, 0.11244362D-17,
810      +  0.11203183D+00, 0.13425405D-01, 0.24867409D-04, 0.25139448D-04,
811      +  0.24876929D-04, 0.24867409D-04, 0.24985115D-01, 0.88225515D-05,
812      +  0.10021714D+00, 0.11789005D-01, 0.28571124D-04, 0.29040474D-04,
813      +  0.28587419D-04, 0.28571124D-04, 0.19490648D-01, 0.14854089D-04,
814      +  0.89543890D-01, 0.10341160D-01, 0.30433528D-04, 0.31041543D-04,
815      +  0.30454492D-04, 0.30433528D-04, 0.15164455D-01, 0.18716337D-04,
816      +  0.79912217D-01, 0.90612886D-02, 0.30974160D-04, 0.31673538D-04,
817      +  0.30998138D-04, 0.30974160D-04, 0.11767860D-01, 0.20974399D-04,
818      +  0.71229380D-01, 0.79309402D-02, 0.30508076D-04, 0.31261180D-04,
819      +  0.30533769D-04, 0.30508076D-04, 0.91152641D-02, 0.21981861D-04,
820      +  0.63410450D-01, 0.69335629D-02, 0.29328990D-04, 0.30106333D-04,
821      +  0.29355400D-04, 0.29328990D-04, 0.70505733D-02, 0.22066269D-04,
822      +  0.56377807D-01, 0.60544377D-02, 0.27677772D-04, 0.28456519D-04,
823      +  0.27704142D-04, 0.27677772D-04, 0.54488746D-02, 0.21497573D-04,
824      +  0.50060242D-01, 0.52804060D-02, 0.25735893D-04, 0.26498771D-04,
825      +  0.25761661D-04, 0.25735893D-04, 0.42101156D-02, 0.20482340D-04/
826         DATA ((CALC(I,K,19),I=1,8), K=10,18) /
827      +  0.44392424D-01, 0.45997051D-02, 0.23643606D-04, 0.24377904D-04,
828      +  0.23668370D-04, 0.23643606D-04, 0.32543735D-02, 0.19182514D-04,
829      +  0.39314381D-01, 0.40018155D-02, 0.21505942D-04, 0.22202662D-04,
830      +  0.21529422D-04, 0.21505942D-04, 0.25184902D-02, 0.17721871D-04,
831      +  0.34771051D-01, 0.34773286D-02, 0.19398514D-04, 0.20051671D-04,
832      +  0.19420532D-04, 0.19398514D-04, 0.19527211D-02, 0.16192231D-04,
833      +  0.30711899D-01, 0.30178307D-02, 0.17373779D-04, 0.17979794D-04,
834      +  0.17394234D-04, 0.17373779D-04, 0.15180062D-02, 0.14660054D-04,
835      +  0.27090563D-01, 0.26158026D-02, 0.15466106D-04, 0.16023302D-04,
836      +  0.15484957D-04, 0.15466106D-04, 0.11839900D-02, 0.13171807D-04,
837      +  0.23864549D-01, 0.22645289D-02, 0.13696033D-04, 0.14204210D-04,
838      +  0.13713284D-04, 0.13696033D-04, 0.92724223D-03, 0.11758474D-04,
839      +  0.20994946D-01, 0.19580180D-02, 0.12073599D-04, 0.12533683D-04,
840      +  0.12089289D-04, 0.12073599D-04, 0.72955528D-03, 0.10439109D-04,
841      +  0.18446167D-01, 0.16909298D-02, 0.10601234D-04, 0.11014980D-04,
842      +  0.10615424D-04, 0.10601234D-04, 0.57689779D-03, 0.92239109D-05,
843      +  0.16185714D-01, 0.14585109D-02, 0.92761477D-05, 0.96459052D-05,
844      +  0.92889183D-05, 0.92761477D-05, 0.45869843D-03, 0.81167844D-05/
845         DATA ((CALC(I,K,19),I=1,8), K=19,20) /
846      +  0.14183953D-01, 0.12565363D-02, 0.80919019D-05, 0.84204158D-05,
847      +  0.81033428D-05, 0.80919019D-05, 0.36686556D-03, 0.71170484D-05,
848      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
849      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
850         DATA ((CALC(I,K,20),I=1,8), K=1,9) /
851      +  0.78207939D-01, 0.71718715D-02, 0.48685484D-05, 0.48685484D-05,
852      +  0.48685484D-05, 0.48685484D-05, 0.17260222D-01, 0.65899163D-18,
853      +  0.69231380D-01, 0.62267774D-02, 0.72000531D-05, 0.73038093D-05,
854      +  0.72035216D-05, 0.72000531D-05, 0.13139860D-01, 0.30850984D-05,
855      +  0.61227629D-01, 0.54026433D-02, 0.86663244D-05, 0.88395579D-05,
856      +  0.86721878D-05, 0.86663244D-05, 0.99947918D-02, 0.51890455D-05,
857      +  0.54085345D-01, 0.46826456D-02, 0.94402700D-05, 0.96573590D-05,
858      +  0.94477236D-05, 0.94402700D-05, 0.75850886D-02, 0.65043408D-05,
859      +  0.47718551D-01, 0.40541438D-02, 0.97259942D-05, 0.99672413D-05,
860      +  0.97344198D-05, 0.97259942D-05, 0.57407678D-02, 0.72495189D-05,
861      +  0.42049033D-01, 0.35060412D-02, 0.96322988D-05, 0.98828263D-05,
862      +  0.96412218D-05, 0.96322988D-05, 0.43380534D-02, 0.75453521D-05,
863      +  0.37005885D-01, 0.30284797D-02, 0.92685924D-05, 0.95174911D-05,
864      +  0.92776571D-05, 0.92685924D-05, 0.32745391D-02, 0.75116772D-05,
865      +  0.32525236D-01, 0.26128175D-02, 0.87291103D-05, 0.89685648D-05,
866      +  0.87380555D-05, 0.87291103D-05, 0.24704846D-02, 0.72515884D-05,
867      +  0.28549325D-01, 0.22514364D-02, 0.80827884D-05, 0.83074776D-05,
868      +  0.80914274D-05, 0.80827884D-05, 0.18643596D-02, 0.68415702D-05/
869         DATA ((CALC(I,K,20),I=1,8), K=10,18) /
870      +  0.25025915D-01, 0.19376193D-02, 0.73826785D-05, 0.75892343D-05,
871      +  0.73908832D-05, 0.73826785D-05, 0.14084418D-02, 0.63411222D-05,
872      +  0.21907749D-01, 0.16654436D-02, 0.66682955D-05, 0.68548566D-05,
873      +  0.66759835D-05, 0.66682955D-05, 0.10659816D-02, 0.57952644D-05,
874      +  0.19152084D-01, 0.14296863D-02, 0.59674032D-05, 0.61332583D-05,
875      +  0.59745275D-05, 0.59674032D-05, 0.80895634D-03, 0.52364600D-05,
876      +  0.16720268D-01, 0.12257434D-02, 0.52988039D-05, 0.54440952D-05,
877      +  0.53053441D-05, 0.52988039D-05, 0.61617187D-03, 0.46875223D-05,
878      +  0.14577375D-01, 0.10495583D-02, 0.46742455D-05, 0.47997302D-05,
879      +  0.46802014D-05, 0.46742455D-05, 0.47144499D-03, 0.41636212D-05,
880      +  0.12691871D-01, 0.89756029D-03, 0.41001724D-05, 0.42070311D-05,
881      +  0.41055580D-05, 0.41001724D-05, 0.36251748D-03, 0.36741181D-05,
882      +  0.11035317D-01, 0.76660894D-03, 0.35793348D-05, 0.36690201D-05,
883      +  0.35841745D-05, 0.35793348D-05, 0.28039012D-03, 0.32242487D-05,
884      +  0.95820983D-02, 0.65394572D-03, 0.31117463D-05, 0.31858637D-05,
885      +  0.31160713D-05, 0.31117463D-05, 0.21831921D-03, 0.28161414D-05,
886      +  0.83091737D-02, 0.55715039D-03, 0.26954812D-05, 0.27556979D-05,
887      +  0.26993268D-05, 0.26954812D-05, 0.17110569D-03, 0.24496689D-05/
888         DATA ((CALC(I,K,20),I=1,8), K=19,20) /
889      +  0.71958512D-02, 0.47410266D-03, 0.23275372D-05, 0.23755137D-05,
890      +  0.23309410D-05, 0.23275372D-05, 0.13497267D-03, 0.21233560D-05,
891      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
892      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
893         DATA ((CALC(I,K,21),I=1,8), K=1,9) /
894      +  0.44126198D-01, 0.28774181D-02, 0.98730783D-06, 0.98730783D-06,
895      +  0.98730783D-06, 0.98730783D-06, 0.83237942D-02, 0.16570082D-19,
896      +  0.38566452D-01, 0.24647091D-02, 0.16316902D-05, 0.16407939D-05,
897      +  0.16327293D-05, 0.16316902D-05, 0.61526696D-02, 0.80814262D-06,
898      +  0.33677597D-01, 0.21100556D-02, 0.20721243D-05, 0.20823724D-05,
899      +  0.20738572D-05, 0.20721243D-05, 0.45490572D-02, 0.13852264D-05,
900      +  0.29373200D-01, 0.18045320D-02, 0.23203461D-05, 0.23267903D-05,
901      +  0.23225203D-05, 0.23203461D-05, 0.33575991D-02, 0.17479170D-05,
902      +  0.25587283D-01, 0.15415184D-02, 0.24333717D-05, 0.24329421D-05,
903      +  0.24357984D-05, 0.24333717D-05, 0.24706562D-02, 0.19568236D-05,
904      +  0.22261127D-01, 0.13153392D-02, 0.24344018D-05, 0.24255218D-05,
905      +  0.24369394D-05, 0.24344018D-05, 0.18156650D-02, 0.20380617D-05,
906      +  0.19342057D-01, 0.11210040D-02, 0.23525948D-05, 0.23348337D-05,
907      +  0.23551400D-05, 0.23525948D-05, 0.13335055D-02, 0.20233070D-05,
908      +  0.16783341D-01, 0.95420544D-03, 0.22171962D-05, 0.21908372D-05,
909      +  0.22196760D-05, 0.22171962D-05, 0.97931452D-03, 0.19439162D-05,
910      +  0.14543347D-01, 0.81120410D-03, 0.20488639D-05, 0.20147058D-05,
911      +  0.20512283D-05, 0.20488639D-05, 0.71978424D-03, 0.18223200D-05/
912         DATA ((CALC(I,K,21),I=1,8), K=10,18) /
913      +  0.12584987D-01, 0.68875151D-03, 0.18635568D-05, 0.18227219D-05,
914      +  0.18657735D-05, 0.18635568D-05, 0.52992811D-03, 0.16759719D-05,
915      +  0.10875222D-01, 0.58402703D-03, 0.16735644D-05, 0.16273312D-05,
916      +  0.16756148D-05, 0.16735644D-05, 0.39128689D-03, 0.15184202D-05,
917      +  0.93846245D-02, 0.49458096D-03, 0.14873320D-05, 0.14370258D-05,
918      +  0.14892073D-05, 0.14873320D-05, 0.29005601D-03, 0.13591698D-05,
919      +  0.80869988D-02, 0.41828666D-03, 0.13103514D-05, 0.12572654D-05,
920      +  0.13120505D-05, 0.13103514D-05, 0.21597406D-03, 0.12046045D-05,
921      +  0.69590429D-02, 0.35329943D-03, 0.11460860D-05, 0.10914244D-05,
922      +  0.11476129D-05, 0.11460860D-05, 0.16173656D-03, 0.10589379D-05,
923      +  0.59800496D-02, 0.29802041D-03, 0.99635391D-06, 0.94119654D-06,
924      +  0.99771627D-06, 0.99635391D-06, 0.12200019D-03, 0.92461894D-06,
925      +  0.51316381D-02, 0.25106512D-03, 0.86169801D-06, 0.80698098D-06,
926      +  0.86290587D-06, 0.86169801D-06, 0.92683828D-04, 0.80272005D-06,
927      +  0.43975147D-02, 0.21123630D-03, 0.74193773D-06, 0.68844610D-06,
928      +  0.74300255D-06, 0.74193773D-06, 0.70911410D-04, 0.69350509D-06,
929      +  0.37632565D-02, 0.17749999D-03, 0.63643257D-06, 0.58480153D-06,
930      +  0.63736646D-06, 0.63643257D-06, 0.54725864D-04, 0.59670595D-06/
931         DATA ((CALC(I,K,21),I=1,8), K=19,20) /
932      +  0.32161163D-02, 0.14896435D-03, 0.54415857D-06, 0.49488097D-06,
933      +  0.54497380D-06, 0.54415857D-06, 0.42606524D-04, 0.51161037D-06,
934      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
935      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
936         DATA ((CALC(I,K,22),I=1,8), K=1,9) /
937      +  0.21580972D-01, 0.92209954D-03, 0.14007288D-06, 0.14007288D-06,
938      +  0.14007288D-06, 0.14007288D-06, 0.34094253D-02, 0.22811987D-19,
939      +  0.18585618D-01, 0.77752973D-03, 0.22352135D-06, 0.19058807D-06,
940      +  0.22376004D-06, 0.22352135D-06, 0.24340210D-02, 0.10846867D-06,
941      +  0.15993270D-01, 0.65539026D-03, 0.30802777D-06, 0.24852177D-06,
942      +  0.30841980D-06, 0.30802777D-06, 0.17410110D-02, 0.21351491D-06,
943      +  0.13745610D-01, 0.55186052D-03, 0.36621607D-06, 0.28633073D-06,
944      +  0.36670050D-06, 0.36621607D-06, 0.12444936D-02, 0.28864325D-06,
945      +  0.11798701D-01, 0.46414836D-03, 0.40311511D-06, 0.30870323D-06,
946      +  0.40364752D-06, 0.40311511D-06, 0.88612427D-03, 0.33951680D-06,
947      +  0.10114392D-01, 0.38992726D-03, 0.41758062D-06, 0.31368737D-06,
948      +  0.41812855D-06, 0.41758062D-06, 0.63037536D-03, 0.36549185D-06,
949      +  0.86588490D-02, 0.32717401D-03, 0.41277224D-06, 0.30369141D-06,
950      +  0.41331294D-06, 0.41277224D-06, 0.44853047D-03, 0.37015640D-06,
951      +  0.74025551D-02, 0.27417276D-03, 0.39489841D-06, 0.28417881D-06,
952      +  0.39541648D-06, 0.39489841D-06, 0.31919521D-03, 0.36007271D-06,
953      +  0.63196564D-02, 0.22946005D-03, 0.36836096D-06, 0.25883661D-06,
954      +  0.36884652D-06, 0.36836096D-06, 0.22754880D-03, 0.33993476D-06/
955         DATA ((CALC(I,K,22),I=1,8), K=10,18) /
956      +  0.53874971D-02, 0.19178564D-03, 0.33665343D-06, 0.23050806D-06,
957      +  0.33710077D-06, 0.33665343D-06, 0.16272555D-03, 0.31347867D-06,
958      +  0.45862357D-02, 0.16008248D-03, 0.30269325D-06, 0.20153905D-06,
959      +  0.30309968D-06, 0.30269325D-06, 0.11675915D-03, 0.28382289D-06,
960      +  0.38985004D-02, 0.13343992D-03, 0.26863798D-06, 0.17359831D-06,
961      +  0.26900300D-06, 0.26863798D-06, 0.84194380D-04, 0.25329154D-06,
962      +  0.33090918D-02, 0.11108101D-03, 0.23591169D-06, 0.14769895D-06,
963      +  0.23623631D-06, 0.23591169D-06, 0.61171621D-04, 0.22344662D-06,
964      +  0.28047231D-02, 0.92343339D-04, 0.20533311D-06, 0.12432414D-06,
965      +  0.20561935D-06, 0.20533311D-06, 0.44770743D-04, 0.19522106D-06,
966      +  0.23737922D-02, 0.76662798D-04, 0.17739050D-06, 0.10369546D-06,
967      +  0.17764102D-06, 0.17739050D-06, 0.32997468D-04, 0.16919749D-06,
968      +  0.20061808D-02, 0.63559531D-04, 0.15232747D-06, 0.85849969D-07,
969      +  0.15254527D-06, 0.15232747D-06, 0.24567206D-04, 0.14569754D-06,
970      +  0.16930761D-02, 0.52625754D-04, 0.13013297D-06, 0.70622292D-07,
971      +  0.13032119D-06, 0.13013297D-06, 0.18493395D-04, 0.12477449D-06,
972      +  0.14268147D-02, 0.43515450D-04, 0.11066644D-06, 0.57761959D-07,
973      +  0.11082822D-06, 0.11066644D-06, 0.14033552D-04, 0.10634084D-06/
974         DATA ((CALC(I,K,22),I=1,8), K=19,20) /
975      +  0.12007428D-02, 0.35935474D-04, 0.93752153D-07, 0.47019945D-07,
976      +  0.93890517D-07, 0.93752153D-07, 0.10740693D-04, 0.90264520D-07,
977      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
978      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
979         DATA ((CALC(I,K,23),I=1,8), K=1,9) /
980      +  0.84858156D-02, 0.20777556D-03, 0.11297044D-07, 0.11297044D-07,
981      +  0.11297044D-07, 0.11297044D-07, 0.10787630D-02, 0.81852830D-19,
982      +  0.71772078D-02, 0.17184644D-03,-0.66395807D-08,-0.41478760D-07,
983      + -0.66053772D-08,-0.66395807D-08, 0.73734465D-03,-0.15743793D-07,
984      +  0.60664390D-02, 0.14211945D-03,-0.20154037D-08,-0.57171527D-07,
985      + -0.19618372D-08,-0.20154037D-08, 0.50633255D-03,-0.93555012D-08,
986      +  0.51211524D-02, 0.11741519D-03, 0.68791051D-08,-0.58842106D-07,
987      +  0.69423207D-08, 0.68791051D-08, 0.34830947D-03, 0.96628812D-09,
988      +  0.43173411D-02, 0.96886760D-04, 0.16322381D-07,-0.53462304D-07,
989      +  0.16388920D-07, 0.16322381D-07, 0.23817074D-03, 0.11565259D-07,
990      +  0.36348173D-02, 0.79853924D-04, 0.23685523D-07,-0.45664853D-07,
991      +  0.23751105D-07, 0.23685523D-07, 0.16267318D-03, 0.19862202D-07,
992      +  0.30559156D-02, 0.65732582D-04, 0.27959697D-07,-0.38052035D-07,
993      +  0.28021657D-07, 0.27959697D-07, 0.11140415D-03, 0.24890346D-07,
994      +  0.25655227D-02, 0.54037278D-04, 0.29830311D-07,-0.31147471D-07,
995      +  0.29887152D-07, 0.29830311D-07, 0.76410313D-04, 0.27369238D-07,
996      +  0.21506722D-02, 0.44362930D-04, 0.29807241D-07,-0.25259246D-07,
997      +  0.29858244D-07, 0.29807241D-07, 0.52448899D-04, 0.27836343D-07/
998         DATA ((CALC(I,K,23),I=1,8), K=10,18) /
999      +  0.18002247D-02, 0.36370383D-04, 0.28397781D-07,-0.20440888D-07,
1000      +  0.28442756D-07, 0.28397781D-07, 0.36149121D-04, 0.26821429D-07,
1001      +  0.15046181D-02, 0.29776112D-04, 0.26176173D-07,-0.16505236D-07,
1002      +  0.26215275D-07, 0.26176173D-07, 0.25121414D-04, 0.24917013D-07,
1003      +  0.12556495D-02, 0.24343044D-04, 0.23524267D-07,-0.13316204D-07,
1004      +  0.23557864D-07, 0.23524267D-07, 0.17574399D-04, 0.22519787D-07,
1005      +  0.10462870D-02, 0.19873096D-04, 0.20716252D-07,-0.10742717D-07,
1006      +  0.20744828D-07, 0.20716252D-07, 0.12362798D-04, 0.19915994D-07,
1007      +  0.87050901D-03, 0.16200928D-04, 0.17966909D-07,-0.86435474D-08,
1008      +  0.17991001D-07, 0.17966909D-07, 0.87997065D-05, 0.17330192D-07,
1009      +  0.72316490D-03, 0.13188630D-04, 0.15391385D-07,-0.69279941D-08,
1010      +  0.15411538D-07, 0.15391385D-07, 0.63519415D-05, 0.14885453D-07,
1011      +  0.59985502D-03, 0.10721318D-04, 0.13045821D-07,-0.55309984D-08,
1012      +  0.13062561D-07, 0.13045821D-07, 0.46234869D-05, 0.12644335D-07,
1013      +  0.49682666D-03, 0.87034365D-05, 0.10964387D-07,-0.43885536D-08,
1014      +  0.10978202D-07, 0.10964387D-07, 0.33981657D-05, 0.10646197D-07,
1015      +  0.41088397D-03, 0.70556024D-05, 0.91516560D-08,-0.34539907D-08,
1016      +  0.91629915D-08, 0.91516560D-08, 0.25366900D-05, 0.88998043D-08/
1017         DATA ((CALC(I,K,23),I=1,8), K=19,20) /
1018      +  0.33930984D-03, 0.57119689D-05, 0.75919874D-08,-0.26946949D-08,
1019      +  0.76012374D-08, 0.75919874D-08, 0.19154362D-05, 0.73928945D-08,
1020      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
1021      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
1022         DATA ((CALC(I,K,24),I=1,8), K=1,9) /
1023      +  0.22702571D-02, 0.24708982D-04, 0.32502839D-09, 0.32502839D-09,
1024      +  0.32502839D-09, 0.32502839D-09, 0.21308933D-03,-0.18486130D-20,
1025      +  0.18744090D-02, 0.19903444D-04,-0.14246691D-07,-0.42350326D-08,
1026      + -0.14249211D-07,-0.14246691D-07, 0.13714357D-03,-0.14501962D-07,
1027      +  0.15469124D-02, 0.16040142D-04,-0.18652171D-07,-0.27630390D-08,
1028      + -0.18656122D-07,-0.18652171D-07, 0.89146620D-04,-0.18852845D-07,
1029      +  0.12749615D-02, 0.12914882D-04,-0.19042325D-07,-0.74608701D-10,
1030      + -0.19046993D-07,-0.19042325D-07, 0.58368090D-04,-0.19199960D-07,
1031      +  0.10492560D-02, 0.10384592D-04,-0.17684678D-07, 0.24683636D-08,
1032      + -0.17689594D-07,-0.17684678D-07, 0.37797614D-04,-0.17808318D-07,
1033      +  0.86227886D-03, 0.83402949D-05,-0.15549653D-07, 0.44773317D-08,
1034      + -0.15554502D-07,-0.15549653D-07, 0.24455882D-04,-0.15646524D-07,
1035      +  0.70756042D-03, 0.66897909D-05,-0.13333272D-07, 0.57206659D-08,
1036      + -0.13337857D-07,-0.13333272D-07, 0.15933279D-04,-0.13409081D-07,
1037      +  0.57970642D-03, 0.53585136D-05,-0.11252145D-07, 0.63275839D-08,
1038      + -0.11256354D-07,-0.11252145D-07, 0.10386156D-04,-0.11311394D-07,
1039      +  0.47420896D-03, 0.42861267D-05,-0.93772561D-08, 0.64667477D-08,
1040      + -0.93810351D-08,-0.93772561D-08, 0.67963855D-05,-0.94235001D-08/
1041         DATA ((CALC(I,K,24),I=1,8), K=10,18) /
1042      +  0.38729093D-03, 0.34233910D-05,-0.77688585D-08, 0.62450390D-08,
1043      + -0.77721932D-08,-0.77688585D-08, 0.44870343D-05,-0.78049036D-08,
1044      +  0.31579344D-03, 0.27302787D-05,-0.64116419D-08, 0.57920215D-08,
1045      + -0.64145433D-08,-0.64116419D-08, 0.29804498D-05,-0.64396984D-08,
1046      +  0.25707652D-03, 0.21742626D-05,-0.52650150D-08, 0.52214625D-08,
1047      + -0.52675098D-08,-0.52650150D-08, 0.19997594D-05,-0.52868229D-08,
1048      +  0.20893582D-03, 0.17288911D-05,-0.43056522D-08, 0.46000381D-08,
1049      + -0.43077759D-08,-0.43056522D-08, 0.13632810D-05,-0.43225792D-08,
1050      +  0.16953312D-03, 0.13726894D-05,-0.35096481D-08, 0.39740039D-08,
1051      + -0.35114402D-08,-0.35096481D-08, 0.94021476D-06,-0.35227681D-08,
1052      +  0.13733760D-03, 0.10882495D-05,-0.28493808D-08, 0.33785446D-08,
1053      + -0.28508814D-08,-0.28493808D-08, 0.65571517D-06,-0.28595357D-08,
1054      +  0.11107635D-03, 0.86147086D-06,-0.23030700D-08, 0.28329892D-08,
1055      + -0.23043179D-08,-0.23030700D-08, 0.46575548D-06,-0.23109189D-08,
1056      +  0.89692587D-04, 0.68094803D-06,-0.18535793D-08, 0.23456287D-08,
1057      + -0.18546104D-08,-0.18535793D-08, 0.33609604D-06,-0.18596374D-08,
1058      +  0.72310507D-04, 0.53747241D-06,-0.14850461D-08, 0.19196948D-08,
1059      + -0.14858934D-08,-0.14850461D-08, 0.24478891D-06,-0.14897156D-08/
1060         DATA ((CALC(I,K,24),I=1,8), K=19,20) /
1061      +  0.58205571D-04, 0.42361984D-06,-0.11839162D-08, 0.15541943D-08,
1062      + -0.11846086D-08,-0.11839162D-08, 0.18047922D-06,-0.11875105D-08,
1063      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
1064      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
1065         DATA ((CALC(I,K,25),I=1,8), K=1,9) /
1066      +  0.24429591D-03, 0.62021186D-06, 0.75420150D-12, 0.75420150D-12,
1067      +  0.75420150D-12, 0.75420150D-12, 0.13318083D-04, 0.95227098D-22,
1068      +  0.19366150D-03, 0.47939321D-06, 0.21511713D-09,-0.40057740D-08,
1069      +  0.21523678D-09, 0.21511713D-09, 0.77271079D-05, 0.21454950D-09,
1070      +  0.15355328D-03, 0.37131637D-06, 0.67681623D-09,-0.60145920D-08,
1071      +  0.67699959D-09, 0.67681623D-09, 0.45920150D-05, 0.67638808D-09,
1072      +  0.12159520D-03, 0.28752233D-06, 0.10810370D-08,-0.69046506D-08,
1073      +  0.10812489D-08, 0.10810370D-08, 0.28193096D-05, 0.10807142D-08,
1074      +  0.96111689D-04, 0.22228227D-06, 0.13333946D-08,-0.71542527D-08,
1075      +  0.13336133D-08, 0.13333946D-08, 0.16691017D-05, 0.13331517D-08,
1076      +  0.75852785D-04, 0.17167953D-06, 0.14702839D-08,-0.69732754D-08,
1077      +  0.14704954D-08, 0.14702839D-08, 0.97059807D-06, 0.14701014D-08,
1078      +  0.59767344D-04, 0.13245573D-06, 0.15134131D-08,-0.65356568D-08,
1079      +  0.15136090D-08, 0.15134131D-08, 0.59145851D-06, 0.15132761D-08,
1080      +  0.47011851D-04, 0.10205625D-06, 0.14805928D-08,-0.59668029D-08,
1081      +  0.14807691D-08, 0.14805928D-08, 0.36409669D-06, 0.14804901D-08,
1082      +  0.36914306D-04, 0.78527064D-07, 0.13973951D-08,-0.53401944D-08,
1083      +  0.13975503D-08, 0.13973951D-08, 0.21651983D-06, 0.13973183D-08/
1084         DATA ((CALC(I,K,25),I=1,8), K=10,18) /
1085      +  0.28934434D-04, 0.60339449D-07, 0.12841196D-08,-0.47044321D-08,
1086      +  0.12842538D-08, 0.12841196D-08, 0.13171744D-06, 0.12840622D-08,
1087      +  0.22638993D-04, 0.46297848D-07, 0.11551000D-08,-0.40916924D-08,
1088      +  0.11552145D-08, 0.11551000D-08, 0.84441910D-07, 0.11550572D-08,
1089      +  0.17681367D-04, 0.35472543D-07, 0.10216026D-08,-0.35203358D-08,
1090      +  0.10216991D-08, 0.10216026D-08, 0.53712308D-07, 0.10215707D-08,
1091      +  0.13784433D-04, 0.27139041D-07, 0.89105708D-09,-0.30004447D-08,
1092      +  0.89113766D-09, 0.89105708D-09, 0.33895730D-07, 0.89103334D-09,
1093      +  0.10726952D-04, 0.20732993D-07, 0.76790480D-09,-0.25365403D-08,
1094      +  0.76797148D-09, 0.76790480D-09, 0.22507340D-07, 0.76788717D-09,
1095      +  0.83326046D-05, 0.15816051D-07, 0.65502489D-09,-0.21287637D-08,
1096      +  0.65507965D-09, 0.65502489D-09, 0.15432815D-07, 0.65501182D-09,
1097      +  0.64611100D-05, 0.12047888D-07, 0.55385358D-09,-0.17747234D-08,
1098      +  0.55389824D-09, 0.55385358D-09, 0.10449599D-07, 0.55384390D-09,
1099      +  0.50010620D-05, 0.91644363D-08, 0.46463511D-09,-0.14706962D-08,
1100      +  0.46467131D-09, 0.46463511D-09, 0.71931193D-08, 0.46462796D-09,
1101      +  0.38641562D-05, 0.69612878D-08, 0.38702015D-09,-0.12120739D-08,
1102      +  0.38704934D-09, 0.38702015D-09, 0.51519281D-08, 0.38701488D-09/
1103         DATA ((CALC(I,K,25),I=1,8), K=19,20) /
1104      +  0.29805372D-05, 0.52805026D-08, 0.32033222D-09,-0.99383802D-09,
1105      +  0.32035562D-09, 0.32033222D-09, 0.37055056D-08, 0.32032833D-09,
1106      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00,
1107      +  0.00000000D+00, 0.00000000D+00, 0.00000000D+00, 0.00000000D+00/
1108 C ---------------------------------------------------------------------
1109       DATA ISTART/0/
1110       SAVE ISTART,OWLAM2,Q02PR
1111 C
1112       IF(ISTART.EQ.0) THEN
1113         ISTART=1
1114         DO 11 K=1,25
1115         DO 11 I=1,20
1116         DO 11 M=1,8
1117    11   CALCO(M,I,K) = CALC(M,I,K)
1118            OWLAM=PAR(1)
1119            OWLAM2=OWLAM**2
1120            Q02PR=PAR(39)
1121            Q2MAX=PAR(40)
1122          ENDIF
1123 C
1124 C the conventions are : q(1)=x*u, q(2)=x*d, q(3)=x*str, q(4)=x*usea,
1125 C                       q(5)=x*dsea, q(6)=x*charm, q(7)=x*gluon
1126 C
1127       X = DX
1128       Q = DQ
1129       Q2 = Q*Q
1130       IDQ2=2
1131       SB=0.
1132       IF(Q2-Q02PR) 1,1,2
1133     2 IF(IDQ2-1) 1,1,3
1134     3 SB= LOG( LOG(Q2/OWLAM2)/ LOG(Q02PR/OWLAM2))
1135     1 CONTINUE
1136       CALL AURPR(1,0,X,SB,XPDF(1))
1137       CALL AURPR(2,0,X,SB,XPDF(2))
1138       CALL AURPR(3,0,X,SB,XPDF(3))
1139       CALL AURPR(4,0,X,SB,XPDF(4))
1140       CALL AURPR(5,0,X,SB,XPDF(5))
1141       CALL AURPR(8,0,X,SB,XPDF(6))
1142       CALL AURPR(7,0,X,SB,XPDF(7))
1143 C
1144       DUPV=XPDF(1) - XPDF(4)
1145       DDNV=XPDF(2) - XPDF(5)
1146       DUSEA=XPDF(4)
1147       DDSEA=XPDF(5)
1148       DSTR=XPDF(3)
1149       DCHM=XPDF(6)
1150       DGL =XPDF(7)
1151 C
1152       RETURN
1153       END