]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PDF/spdf/sfmsqa.F
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
This commit was generated by cvs2svn to compensate for changes in r1907,
[u/mrichter/AliRoot.git] / PDF / spdf / sfmsqa.F
CommitLineData
21886bb6 1*
2* $Id$
3*
4* $Log$
5* Revision 1.1.1.1 1996/10/30 08:32:13 cernlib
6* Version 7.04
7*
8* Revision 1.2 1996/10/28 15:09:53 pdflib
9* PDFLIB new version 7.04
10*
11* Revision 1.1.1.1 1996/04/12 15:30:16 plothow
12* Version 7.01
13*
14*
15#include "pdf/pilot.h"
16 SUBROUTINE SFMSQA(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,GLU)
17C
18C THIS IS THE NEW "A" FIT -- May 1994 -- low Q^2 range
19C
20#include "pdf/impdp.inc"
21C
22 PARAMETER(NX=47)
23 PARAMETER(NTENTH=21)
24 DIMENSION F(8,NX,8),G(8),XX(NX),N0(8)
25 DATA ((F(I,1,K),I=1,5),F(7,1,K),F(6,1,K),F(8,1,K),K=1,7)/
26 + 0.00232, 0.00543, 4.23812, 1.92427, 0.00000,
27 + 0.00000, 0.53025, 1.91803,
28 + 0.00262, 0.00566, 6.42705, 1.98441, 0.00000,
29 + 0.00000, 0.62266, 1.97884,
30 + 0.00289, 0.00586, 9.43898, 2.07363, 0.00000,
31 + 0.00000, 0.73934, 2.06869,
32 + 0.00315, 0.00605, 13.16286, 2.19098, 0.00000,
33 + 0.00000, 0.88050, 2.18663,
34 + 0.00339, 0.00622, 17.54780, 2.32948, 0.00000,
35 + 0.00000, 1.04003, 2.32568,
36 + 0.00362, 0.00638, 22.55753, 2.48780, 0.05952,
37 + 0.00000, 1.21723, 2.48453,
38 + 0.00383, 0.00653, 28.15842, 2.66466, 0.26782,
39 + 0.00000, 1.41125, 2.66188 /
40 DATA ((F(I,2,K),I=1,5),F(7,2,K),F(6,2,K),F(8,2,K),K=1,7)/
41 + 0.00364, 0.00691, 3.22632, 1.55076, 0.00000,
42 + 0.00000, 0.42116, 1.54806,
43 + 0.00407, 0.00722, 5.15995, 1.60401, 0.00000,
44 + 0.00000, 0.50030, 1.60175,
45 + 0.00446, 0.00749, 7.71045, 1.67947, 0.00000,
46 + 0.00000, 0.59775, 1.67760,
47 + 0.00483, 0.00774, 10.78281, 1.77603, 0.00000,
48 + 0.00000, 0.71335, 1.77454,
49 + 0.00516, 0.00797, 14.34431, 1.88827, 0.00000,
50 + 0.00000, 0.84238, 1.88712,
51 + 0.00548, 0.00819, 18.35565, 2.01500, 0.04882,
52 + 0.00000, 0.98416, 2.01418,
53 + 0.00579, 0.00839, 22.78315, 2.15515, 0.21426,
54 + 0.00000, 1.13797, 2.15464 /
55 DATA ((F(I,3,K),I=1,5),F(7,3,K),F(6,3,K),F(8,3,K),K=1,7)/
56 + 0.00558, 0.00879, 2.45661, 1.24947, 0.00000,
57 + 0.00000, 0.33447, 1.24914,
58 + 0.00619, 0.00921, 4.14420, 1.29621, 0.00000,
59 + 0.00000, 0.40201, 1.29616,
60 + 0.00675, 0.00958, 6.30016, 1.35987, 0.00000,
61 + 0.00000, 0.48332, 1.36007,
62 + 0.00727, 0.00991, 8.83462, 1.43927, 0.00000,
63 + 0.00000, 0.57796, 1.43970,
64 + 0.00774, 0.01022, 11.72666, 1.53017, 0.00000,
65 + 0.00000, 0.68229, 1.53081,
66 + 0.00820, 0.01052, 14.93678, 1.63155, 0.04004,
67 + 0.00000, 0.79569, 1.63239,
68 + 0.00862, 0.01079, 18.43334, 1.74251, 0.17137,
69 + 0.00000, 0.91755, 1.74354 /
70 DATA ((F(I,4,K),I=1,5),F(7,4,K),F(6,4,K),F(8,4,K),K=1,7)/
71 + 0.00710, 0.01013, 2.09495, 1.10097, 0.00000,
72 + 0.00000, 0.29229, 1.10164,
73 + 0.00785, 0.01062, 3.64657, 1.14412, 0.00000,
74 + 0.00000, 0.35376, 1.14500,
75 + 0.00854, 0.01106, 5.59929, 1.20168, 0.00000,
76 + 0.00000, 0.42685, 1.20276,
77 + 0.00917, 0.01146, 7.86366, 1.27246, 0.00000,
78 + 0.00000, 0.51103, 1.27372,
79 + 0.00976, 0.01183, 10.42364, 1.35279, 0.00000,
80 + 0.00000, 0.60316, 1.35421,
81 + 0.01031, 0.01218, 13.24045, 1.44172, 0.03565,
82 + 0.00000, 0.70264, 1.44329,
83 + 0.01083, 0.01250, 16.28428, 1.53846, 0.15036,
84 + 0.00000, 0.80894, 1.54017 /
85 DATA ((F(I,5,K),I=1,5),F(7,5,K),F(6,5,K),F(8,5,K),K=1,7)/
86 + 0.00840, 0.01120, 1.87136, 1.00635, 0.00000,
87 + 0.00000, 0.26564, 1.00760,
88 + 0.00927, 0.01175, 3.33089, 1.04705, 0.00000,
89 + 0.00000, 0.32310, 1.04848,
90 + 0.01006, 0.01225, 5.15064, 1.10062, 0.00000,
91 + 0.00000, 0.39086, 1.10222,
92 + 0.01079, 0.01271, 7.24090, 1.16583, 0.00000,
93 + 0.00000, 0.46831, 1.16759,
94 + 0.01147, 0.01312, 9.58839, 1.23939, 0.00000,
95 + 0.00000, 0.55265, 1.24128,
96 + 0.01210, 0.01351, 12.15519, 1.32041, 0.03284,
97 + 0.00000, 0.64329, 1.32243,
98 + 0.01271, 0.01388, 14.91292, 1.40814, 0.13702,
99 + 0.00000, 0.73975, 1.41028 /
100 DATA ((F(I,6,K),I=1,5),F(7,6,K),F(6,6,K),F(8,6,K),K=1,7)/
101 + 0.00955, 0.01211, 1.71468, 0.93853, 0.00000,
102 + 0.00000, 0.24666, 0.94017,
103 + 0.01052, 0.01272, 3.10549, 0.97740, 0.00000,
104 + 0.00000, 0.30117, 0.97921,
105 + 0.01141, 0.01327, 4.82815, 1.02804, 0.00000,
106 + 0.00000, 0.36505, 1.03000,
107 + 0.01223, 0.01377, 6.79254, 1.08923, 0.00000,
108 + 0.00000, 0.43767, 1.09133,
109 + 0.01298, 0.01423, 8.98726, 1.15792, 0.00000,
110 + 0.00000, 0.51641, 1.16015,
111 + 0.01369, 0.01466, 11.37514, 1.23328, 0.03081,
112 + 0.00000, 0.60073, 1.23562,
113 + 0.01437, 0.01506, 13.92903, 1.31459, 0.12748,
114 + 0.00000, 0.69017, 1.31704 /
115 DATA ((F(I,7,K),I=1,5),F(7,7,K),F(6,7,K),F(8,7,K),K=1,7)/
116 + 0.01411, 0.01545, 1.30792, 0.75531, 0.00000,
117 + 0.00000, 0.19600, 0.75794,
118 + 0.01550, 0.01627, 2.50161, 0.78885, 0.00000,
119 + 0.00000, 0.24221, 0.79165,
120 + 0.01677, 0.01701, 3.95367, 0.83129, 0.00000,
121 + 0.00000, 0.29537, 0.83423,
122 + 0.01792, 0.01768, 5.57292, 0.88143, 0.00000,
123 + 0.00000, 0.35475, 0.88450,
124 + 0.01899, 0.01830, 7.35237, 0.93689, 0.00000,
125 + 0.00000, 0.41833, 0.94007,
126 + 0.01999, 0.01887, 9.25794, 0.99695, 0.02529,
127 + 0.00000, 0.48564, 1.00024,
128 + 0.02093, 0.01941, 11.26643, 1.06105, 0.10183,
129 + 0.00000, 0.55629, 1.06444 /
130 DATA ((F(I,8,K),I=1,5),F(7,8,K),F(6,8,K),F(8,8,K),K=1,7)/
131 + 0.02070, 0.01979, 0.99986, 0.60732, 0.00000,
132 + 0.00000, 0.15592, 0.61082,
133 + 0.02268, 0.02090, 2.02268, 0.63609, 0.00000,
134 + 0.00000, 0.19503, 0.63980,
135 + 0.02447, 0.02190, 3.24642, 0.67153, 0.00000,
136 + 0.00000, 0.23921, 0.67542,
137 + 0.02610, 0.02280, 4.58011, 0.71248, 0.00000,
138 + 0.00000, 0.28769, 0.71654,
139 + 0.02760, 0.02362, 6.02017, 0.75711, 0.00000,
140 + 0.00000, 0.33893, 0.76132,
141 + 0.02900, 0.02439, 7.53631, 0.80481, 0.02079,
142 + 0.00000, 0.39253, 0.80915,
143 + 0.03032, 0.02511, 9.10952, 0.85512, 0.08122,
144 + 0.00000, 0.44819, 0.85959 /
145 DATA ((F(I,9,K),I=1,5),F(7,9,K),F(6,9,K),F(8,9,K),K=1,7)/
146 + 0.02584, 0.02295, 0.85607, 0.53431, 0.00000,
147 + 0.00000, 0.13654, 0.53837,
148 + 0.02827, 0.02427, 1.79177, 0.56056, 0.00000,
149 + 0.00000, 0.17198, 0.56487,
150 + 0.03046, 0.02546, 2.89947, 0.59238, 0.00000,
151 + 0.00000, 0.21160, 0.59693,
152 + 0.03245, 0.02653, 4.08932, 0.62869, 0.00000,
153 + 0.00000, 0.25460, 0.63343,
154 + 0.03428, 0.02751, 5.35969, 0.66790, 0.00000,
155 + 0.00000, 0.29971, 0.67282,
156 + 0.03599, 0.02842, 6.68285, 0.70946, 0.01856,
157 + 0.00000, 0.34653, 0.71456,
158 + 0.03759, 0.02928, 8.04229, 0.75299, 0.07107,
159 + 0.00000, 0.39483, 0.75825 /
160 DATA ((F(I,10,K),I=1,5),F(7,10,K),F(6,10,K),F(8,10,K),K=1,7)/
161 + 0.03025, 0.02555, 0.76776, 0.48778, 0.00000,
162 + 0.00000, 0.12435, 0.49229,
163 + 0.03306, 0.02705, 1.64767, 0.51235, 0.00000,
164 + 0.00000, 0.15741, 0.51717,
165 + 0.03558, 0.02839, 2.68028, 0.54182, 0.00000,
166 + 0.00000, 0.19406, 0.54690,
167 + 0.03787, 0.02960, 3.77709, 0.57510, 0.00000,
168 + 0.00000, 0.23353, 0.58042,
169 + 0.03997, 0.03071, 4.93803, 0.61081, 0.00000,
170 + 0.00000, 0.27468, 0.61634,
171 + 0.04193, 0.03174, 6.13732, 0.64844, 0.01714,
172 + 0.00000, 0.31716, 0.65417,
173 + 0.04376, 0.03271, 7.36026, 0.68763, 0.06460,
174 + 0.00000, 0.36076, 0.69355 /
175 DATA ((F(I,11,K),I=1,5),F(7,11,K),F(6,11,K),F(8,11,K),K=1,7)/
176 + 0.03425, 0.02784, 0.70475, 0.45443, 0.00000,
177 + 0.00000, 0.11569, 0.45936,
178 + 0.03738, 0.02950, 1.54754, 0.47783, 0.00000,
179 + 0.00000, 0.14703, 0.48308,
180 + 0.04019, 0.03099, 2.52696, 0.50558, 0.00000,
181 + 0.00000, 0.18150, 0.51113,
182 + 0.04274, 0.03233, 3.55628, 0.53668, 0.00000,
183 + 0.00000, 0.21838, 0.54247,
184 + 0.04507, 0.03355, 4.63730, 0.56984, 0.00000,
185 + 0.00000, 0.25663, 0.57587,
186 + 0.04724, 0.03469, 5.74604, 0.60460, 0.01613,
187 + 0.00000, 0.29594, 0.61084,
188 + 0.04928, 0.03574, 6.86937, 0.64065, 0.05990,
189 + 0.00000, 0.33610, 0.64709 /
190 DATA ((F(I,12,K),I=1,5),F(7,12,K),F(6,12,K),F(8,12,K),K=1,7)/
191 + 0.05034, 0.03672, 0.54955, 0.36447, 0.00000,
192 + 0.00000, 0.09298, 0.37107,
193 + 0.05473, 0.03900, 1.28832, 0.38445, 0.00000,
194 + 0.00000, 0.11948, 0.39149,
195 + 0.05865, 0.04103, 2.11904, 0.40739, 0.00000,
196 + 0.00000, 0.14795, 0.41483,
197 + 0.06217, 0.04285, 2.96136, 0.43239, 0.00000,
198 + 0.00000, 0.17770, 0.44019,
199 + 0.06539, 0.04450, 3.82264, 0.45855, 0.00000,
200 + 0.00000, 0.20805, 0.46667,
201 + 0.06837, 0.04602, 4.68465, 0.48551, 0.01341,
202 + 0.00000, 0.23875, 0.49393,
203 + 0.07115, 0.04744, 5.53899, 0.51305, 0.04724,
204 + 0.00000, 0.26968, 0.52174 /
205 DATA ((F(I,13,K),I=1,5),F(7,13,K),F(6,13,K),F(8,13,K),K=1,7)/
206 + 0.07424, 0.04960, 0.45006, 0.29268, 0.00000,
207 + 0.00000, 0.07588, 0.30160,
208 + 0.08032, 0.05272, 1.10456, 0.30967, 0.00000,
209 + 0.00000, 0.09822, 0.31921,
210 + 0.08571, 0.05548, 1.80790, 0.32847, 0.00000,
211 + 0.00000, 0.12157, 0.33855,
212 + 0.09051, 0.05792, 2.49068, 0.34832, 0.00000,
213 + 0.00000, 0.14534, 0.35889,
214 + 0.09488, 0.06014, 3.16635, 0.36864, 0.00000,
215 + 0.00000, 0.16911, 0.37964,
216 + 0.09889, 0.06216, 3.82305, 0.38917, 0.01126,
217 + 0.00000, 0.19272, 0.40057,
218 + 0.10262, 0.06404, 4.45704, 0.40977, 0.03700,
219 + 0.00000, 0.21612, 0.42153 /
220 DATA ((F(I,14,K),I=1,5),F(7,14,K),F(6,14,K),F(8,14,K),K=1,7)/
221 + 0.09337, 0.05991, 0.42083, 0.25795, 0.00000,
222 + 0.00000, 0.06825, 0.26877,
223 + 0.10067, 0.06365, 1.03073, 0.27332, 0.00000,
224 + 0.00000, 0.08842, 0.28489,
225 + 0.10710, 0.06692, 1.66542, 0.28992, 0.00000,
226 + 0.00000, 0.10912, 0.30215,
227 + 0.11281, 0.06980, 2.26414, 0.30711, 0.00000,
228 + 0.00000, 0.12983, 0.31991,
229 + 0.11797, 0.07239, 2.84358, 0.32444, 0.00000,
230 + 0.00000, 0.15026, 0.33776,
231 + 0.12269, 0.07475, 3.39582, 0.34172, 0.01025,
232 + 0.00000, 0.17033, 0.35551,
233 + 0.12706, 0.07692, 3.91963, 0.35886, 0.03197,
234 + -0.00001, 0.18999, 0.37308 /
235 DATA ((F(I,15,K),I=1,5),F(7,15,K),F(6,15,K),F(8,15,K),K=1,7)/
236 + 0.11001, 0.06899, 0.41558, 0.23634, 0.00000,
237 + 0.00000, 0.06388, 0.24884,
238 + 0.11827, 0.07320, 0.99416, 0.25060, 0.00000,
239 + 0.00000, 0.08260, 0.26394,
240 + 0.12553, 0.07686, 1.58125, 0.26572, 0.00000,
241 + 0.00000, 0.10153, 0.27981,
242 + 0.13194, 0.08007, 2.12318, 0.28114, 0.00000,
243 + 0.00000, 0.12023, 0.29587,
244 + 0.13771, 0.08294, 2.63867, 0.29650, 0.00000,
245 + 0.00000, 0.13849, 0.31181,
246 + 0.14297, 0.08555, 3.12247, 0.31167, 0.00963,
247 + 0.00000, 0.15625, 0.32750,
248 + 0.14782, 0.08793, 3.57509, 0.32656, 0.02878,
249 + -0.00001, 0.17350, 0.34287 /
250 DATA ((F(I,16,K),I=1,5),F(7,16,K),F(6,16,K),F(8,16,K),K=1,7)/
251 + 0.12504, 0.07729, 0.42079, 0.22117, 0.00000,
252 + 0.00000, 0.06107, 0.23523,
253 + 0.13412, 0.08187, 0.97466, 0.23457, 0.00000,
254 + 0.00000, 0.07870, 0.24957,
255 + 0.14205, 0.08584, 1.52498, 0.24857, 0.00000,
256 + 0.00000, 0.09633, 0.26439,
257 + 0.14903, 0.08930, 2.02410, 0.26267, 0.00000,
258 + 0.00000, 0.11355, 0.27919,
259 + 0.15529, 0.09239, 2.49209, 0.27658, 0.00000,
260 + 0.00000, 0.13022, 0.29372,
261 + 0.16099, 0.09517, 2.92573, 0.29019, 0.00921,
262 + 0.00000, 0.14631, 0.30789,
263 + 0.16621, 0.09771, 3.32673, 0.30346, 0.02650,
264 + 0.00000, 0.16182, 0.32167 /
265 DATA ((F(I,17,K),I=1,5),F(7,17,K),F(6,17,K),F(8,17,K),K=1,7)/
266 + 0.18747, 0.11240, 0.49179, 0.18182, 0.00000,
267 + 0.00000, 0.05541, 0.20302,
268 + 0.19924, 0.11804, 0.96223, 0.19246, 0.00000,
269 + 0.00000, 0.06983, 0.21485,
270 + 0.20936, 0.12282, 1.39560, 0.20301, 0.00000,
271 + 0.00000, 0.08362, 0.22640,
272 + 0.21810, 0.12689, 1.76565, 0.21315, 0.00000,
273 + 0.00000, 0.09656, 0.23738,
274 + 0.22585, 0.13046, 2.09508, 0.22280, 0.00000,
275 + 0.00000, 0.10870, 0.24776,
276 + 0.23280, 0.13361, 2.38648, 0.23194, 0.00815,
277 + 0.00000, 0.12008, 0.25755,
278 + 0.23909, 0.13644, 2.64474, 0.24060, 0.02040,
279 + 0.00000, 0.13077, 0.26678 /
280 DATA ((F(I,18,K),I=1,5),F(7,18,K),F(6,18,K),F(8,18,K),K=1,7)/
281 + 0.28533, 0.16728, 0.65201, 0.15062, 0.00000,
282 + 0.00000, 0.05345, 0.18428,
283 + 0.29917, 0.17291, 1.01515, 0.15795, 0.00000,
284 + 0.00000, 0.06450, 0.19272,
285 + 0.31073, 0.17747, 1.31959, 0.16473, 0.00000,
286 + 0.00000, 0.07446, 0.20037,
287 + 0.32045, 0.18119, 1.56052, 0.17084, 0.00000,
288 + 0.00000, 0.08333, 0.20718,
289 + 0.32885, 0.18432, 1.76050, 0.17637, 0.00000,
290 + 0.00000, 0.09129, 0.21328,
291 + 0.33619, 0.18698, 1.92619, 0.18137, 0.00727,
292 + 0.00000, 0.09847, 0.21874,
293 + 0.34271, 0.18927, 2.06404, 0.18589, 0.01537,
294 + -0.00001, 0.10497, 0.22364 /
295 DATA ((F(I,19,K),I=1,5),F(7,19,K),F(6,19,K),F(8,19,K),K=1,7)/
296 + 0.36729, 0.21019, 0.77868, 0.13188, 0.00000,
297 + 0.00000, 0.05237, 0.17581,
298 + 0.38087, 0.21442, 1.06122, 0.13689, 0.00000,
299 + 0.00000, 0.06115, 0.18142,
300 + 0.39186, 0.21764, 1.28184, 0.14126, 0.00000,
301 + 0.00000, 0.06872, 0.18618,
302 + 0.40081, 0.22007, 1.44554, 0.14496, 0.00000,
303 + 0.00000, 0.07516, 0.19012,
304 + 0.40832, 0.22197, 1.57317, 0.14815, 0.00000,
305 + 0.00000, 0.08074, 0.19344,
306 + 0.41470, 0.22344, 1.67235, 0.15087, 0.00664,
307 + 0.00000, 0.08561, 0.19624,
308 + 0.42019, 0.22459, 1.74940, 0.15320, 0.01264,
309 + 0.00000, 0.08988, 0.19859 /
310 DATA ((F(I,20,K),I=1,5),F(7,20,K),F(6,20,K),F(8,20,K),K=1,7)/
311 + 0.43943, 0.24421, 0.87132, 0.11596, 0.00000,
312 + 0.00000, 0.05027, 0.16762,
313 + 0.45119, 0.24638, 1.08735, 0.11925, 0.00000,
314 + 0.00000, 0.05727, 0.17074,
315 + 0.46038, 0.24769, 1.24495, 0.12190, 0.00000,
316 + 0.00000, 0.06305, 0.17318,
317 + 0.46749, 0.24839, 1.35369, 0.12398, 0.00000,
318 + 0.00000, 0.06779, 0.17492,
319 + 0.47318, 0.24870, 1.43230, 0.12564, 0.00068,
320 + 0.00000, 0.07176, 0.17621,
321 + 0.47780, 0.24870, 1.48820, 0.12692, 0.00602,
322 + 0.00000, 0.07509, 0.17712,
323 + 0.48154, 0.24849, 1.52710, 0.12792, 0.01066,
324 + -0.00001, 0.07792, 0.17772 /
325 DATA ((F(I,21,K),I=1,5),F(7,21,K),F(6,21,K),F(8,21,K),K=1,7)/
326 + 0.50349, 0.27060, 0.93597, 0.10110, 0.00000,
327 + 0.00000, 0.04723, 0.15793,
328 + 0.51234, 0.27032, 1.09539, 0.10314, 0.00000,
329 + 0.00000, 0.05277, 0.15902,
330 + 0.51876, 0.26952, 1.20287, 0.10461, 0.00000,
331 + 0.00000, 0.05718, 0.15957,
332 + 0.52327, 0.26834, 1.26983, 0.10560, 0.00000,
333 + 0.00000, 0.06063, 0.15964,
334 + 0.52654, 0.26694, 1.31262, 0.10625, 0.00107,
335 + 0.00000, 0.06341, 0.15943,
336 + 0.52884, 0.26542, 1.33802, 0.10663, 0.00537,
337 + 0.00000, 0.06566, 0.15898,
338 + 0.53041, 0.26380, 1.35090, 0.10680, 0.00904,
339 + 0.00000, 0.06749, 0.15835 /
340 DATA ((F(I,22,K),I=1,5),F(7,22,K),F(6,22,K),F(8,22,K),K=1,7)/
341 + 0.57288, 0.29429, 0.98361, 0.08377, 0.00000,
342 + 0.00000, 0.04242, 0.14370,
343 + 0.57693, 0.29082, 1.08372, 0.08476, 0.00000,
344 + 0.00000, 0.04654, 0.14280,
345 + 0.57901, 0.28728, 1.14164, 0.08529, 0.00000,
346 + 0.00000, 0.04963, 0.14161,
347 + 0.57959, 0.28371, 1.16898, 0.08543, 0.00000,
348 + 0.00000, 0.05191, 0.14021,
349 + 0.57926, 0.28021, 1.17918, 0.08531, 0.00123,
350 + 0.00000, 0.05364, 0.13871,
351 + 0.57826, 0.27680, 1.17779, 0.08501, 0.00457,
352 + 0.00000, 0.05495, 0.13713,
353 + 0.57673, 0.27349, 1.16857, 0.08458, 0.00737,
354 + -0.00001, 0.05594, 0.13551 /
355 DATA ((F(I,23,K),I=1,5),F(7,23,K),F(6,23,K),F(8,23,K),K=1,7)/
356 + 0.63048, 0.30915, 1.00111, 0.06818, 0.00000,
357 + 0.00000, 0.03692, 0.12784,
358 + 0.62891, 0.30248, 1.05264, 0.06852, 0.00000,
359 + 0.00000, 0.03999, 0.12557,
360 + 0.62606, 0.29626, 1.07218, 0.06847, 0.00000,
361 + 0.00000, 0.04216, 0.12327,
362 + 0.62227, 0.29042, 1.07015, 0.06812, 0.00000,
363 + 0.00000, 0.04363, 0.12094,
364 + 0.61801, 0.28497, 1.05720, 0.06761, 0.00119,
365 + 0.00000, 0.04466, 0.11866,
366 + 0.61343, 0.27986, 1.03756, 0.06698, 0.00380,
367 + 0.00000, 0.04536, 0.11642,
368 + 0.60864, 0.27503, 1.01392, 0.06628, 0.00596,
369 + -0.00003, 0.04583, 0.11422 /
370 DATA ((F(I,24,K),I=1,5),F(7,24,K),F(6,24,K),F(8,24,K),K=1,7)/
371 + 0.67676, 0.31587, 0.99514, 0.05418, 0.00000,
372 + 0.00000, 0.03156, 0.11143,
373 + 0.66899, 0.30635, 1.00701, 0.05425, 0.00000,
374 + 0.00000, 0.03376, 0.10833,
375 + 0.66082, 0.29778, 0.99689, 0.05396, 0.00000,
376 + 0.00000, 0.03522, 0.10534,
377 + 0.65240, 0.28999, 0.97320, 0.05343, 0.00000,
378 + 0.00000, 0.03611, 0.10251,
379 + 0.64407, 0.28287, 0.94399, 0.05278, 0.00109,
380 + 0.00000, 0.03666, 0.09982,
381 + 0.63589, 0.27630, 0.91213, 0.05205, 0.00312,
382 + 0.00000, 0.03697, 0.09728,
383 + 0.62786, 0.27019, 0.87934, 0.05126, 0.00479,
384 + -0.00001, 0.03709, 0.09488 /
385 DATA ((F(I,25,K),I=1,5),F(7,25,K),F(6,25,K),F(8,25,K),K=1,7)/
386 + 0.71132, 0.31625, 0.97091, 0.04237, 0.00000,
387 + 0.00000, 0.02619, 0.09523,
388 + 0.69723, 0.30412, 0.95082, 0.04228, 0.00000,
389 + 0.00000, 0.02783, 0.09170,
390 + 0.68374, 0.29348, 0.91809, 0.04190, 0.00000,
391 + 0.00000, 0.02882, 0.08844,
392 + 0.67077, 0.28400, 0.87878, 0.04130, 0.00000,
393 + 0.00000, 0.02936, 0.08542,
394 + 0.65848, 0.27548, 0.83847, 0.04063, 0.00090,
395 + 0.00000, 0.02962, 0.08261,
396 + 0.64682, 0.26772, 0.79881, 0.03991, 0.00251,
397 + 0.00000, 0.02968, 0.08002,
398 + 0.63571, 0.26057, 0.76061, 0.03916, 0.00380,
399 + -0.00002, 0.02961, 0.07759 /
400 DATA ((F(I,26,K),I=1,5),F(7,26,K),F(6,26,K),F(8,26,K),K=1,7)/
401 + 0.73462, 0.31108, 0.93107, 0.03255, 0.00000,
402 + 0.00000, 0.02129, 0.07983,
403 + 0.71430, 0.29676, 0.88680, 0.03241, 0.00000,
404 + 0.00000, 0.02250, 0.07623,
405 + 0.69565, 0.28441, 0.83766, 0.03201, 0.00000,
406 + 0.00000, 0.02319, 0.07296,
407 + 0.67837, 0.27354, 0.78764, 0.03144, 0.00000,
408 + 0.00000, 0.02348, 0.07003,
409 + 0.66238, 0.26391, 0.74029, 0.03083, 0.00072,
410 + 0.00000, 0.02358, 0.06730,
411 + 0.64749, 0.25524, 0.69616, 0.03017, 0.00199,
412 + 0.00000, 0.02352, 0.06481,
413 + 0.63353, 0.24733, 0.65531, 0.02949, 0.00299,
414 + -0.00002, 0.02335, 0.06253 /
415 DATA ((F(I,27,K),I=1,5),F(7,27,K),F(6,27,K),F(8,27,K),K=1,7)/
416 + 0.74766, 0.30066, 0.87914, 0.02401, 0.00000,
417 + 0.00000, 0.01723, 0.06581,
418 + 0.72123, 0.28478, 0.81691, 0.02404, 0.00000,
419 + 0.00000, 0.01807, 0.06235,
420 + 0.69764, 0.27126, 0.75668, 0.02381, 0.00000,
421 + 0.00000, 0.01849, 0.05925,
422 + 0.67627, 0.25949, 0.70017, 0.02341, 0.00000,
423 + 0.00000, 0.01860, 0.05648,
424 + 0.65688, 0.24911, 0.64913, 0.02292, 0.00062,
425 + 0.00000, 0.01855, 0.05399,
426 + 0.63909, 0.23985, 0.60318, 0.02241, 0.00160,
427 + 0.00000, 0.01840, 0.05173,
428 + 0.62259, 0.23148, 0.56179, 0.02188, 0.00236,
429 + 0.00002, 0.01819, 0.04966 /
430 DATA ((F(I,28,K),I=1,5),F(7,28,K),F(6,28,K),F(8,28,K),K=1,7)/
431 + 0.74936, 0.28783, 0.82378, 0.01828, 0.00000,
432 + 0.00000, 0.01309, 0.05326,
433 + 0.71761, 0.27055, 0.74720, 0.01821, 0.00000,
434 + 0.00000, 0.01381, 0.05010,
435 + 0.68980, 0.25599, 0.67896, 0.01795, 0.00000,
436 + 0.00000, 0.01415, 0.04731,
437 + 0.66498, 0.24349, 0.61850, 0.01757, 0.00000,
438 + 0.00000, 0.01424, 0.04482,
439 + 0.64273, 0.23257, 0.56585, 0.01714, 0.00028,
440 + 0.00000, 0.01418, 0.04263,
441 + 0.62253, 0.22290, 0.51975, 0.01669, 0.00107,
442 + 0.00000, 0.01404, 0.04065,
443 + 0.60398, 0.21422, 0.47915, 0.01623, 0.00169,
444 + -0.00014, 0.01382, 0.03885 /
445 DATA ((F(I,29,K),I=1,5),F(7,29,K),F(6,29,K),F(8,29,K),K=1,7)/
446 + 0.74266, 0.27081, 0.75762, 0.01278, 0.00000,
447 + 0.00000, 0.01026, 0.04237,
448 + 0.70594, 0.25285, 0.67378, 0.01289, 0.00000,
449 + 0.00000, 0.01077, 0.03958,
450 + 0.67427, 0.23786, 0.60220, 0.01279, 0.00000,
451 + 0.00000, 0.01097, 0.03715,
452 + 0.64640, 0.22509, 0.54096, 0.01257, 0.00000,
453 + 0.00000, 0.01097, 0.03501,
454 + 0.62168, 0.21401, 0.48891, 0.01227, 0.00029,
455 + 0.00000, 0.01087, 0.03314,
456 + 0.59943, 0.20426, 0.44429, 0.01198, 0.00089,
457 + 0.00000, 0.01072, 0.03143,
458 + 0.57919, 0.19555, 0.40568, 0.01166, 0.00134,
459 + -0.00002, 0.01052, 0.02990 /
460 DATA ((F(I,30,K),I=1,5),F(7,30,K),F(6,30,K),F(8,30,K),K=1,7)/
461 + 0.72685, 0.25224, 0.69187, 0.00899, 0.00000,
462 + 0.00000, 0.00768, 0.03313,
463 + 0.68604, 0.23389, 0.60286, 0.00915, 0.00000,
464 + 0.00000, 0.00808, 0.03076,
465 + 0.65127, 0.21871, 0.52998, 0.00910, 0.00000,
466 + 0.00000, 0.00824, 0.02870,
467 + 0.62098, 0.20591, 0.46960, 0.00897, 0.00000,
468 + 0.00000, 0.00824, 0.02690,
469 + 0.59438, 0.19486, 0.41941, 0.00877, 0.00019,
470 + 0.00000, 0.00814, 0.02534,
471 + 0.57066, 0.18518, 0.37721, 0.00852, 0.00065,
472 + 0.00000, 0.00800, 0.02398,
473 + 0.54919, 0.17662, 0.34122, 0.00829, 0.00100,
474 + -0.00004, 0.00784, 0.02270 /
475 DATA ((F(I,31,K),I=1,5),F(7,31,K),F(6,31,K),F(8,31,K),K=1,7)/
476 + 0.70337, 0.23219, 0.62540, 0.00613, 0.00000,
477 + 0.00000, 0.00563, 0.02548,
478 + 0.65924, 0.21385, 0.53420, 0.00632, 0.00000,
479 + 0.00000, 0.00597, 0.02348,
480 + 0.62207, 0.19880, 0.46215, 0.00631, 0.00000,
481 + 0.00000, 0.00609, 0.02181,
482 + 0.59005, 0.18618, 0.40404, 0.00623, 0.00000,
483 + 0.00000, 0.00610, 0.02033,
484 + 0.56211, 0.17538, 0.35676, 0.00609, 0.00015,
485 + 0.00000, 0.00602, 0.01909,
486 + 0.53737, 0.16598, 0.31763, 0.00593, 0.00050,
487 + 0.00000, 0.00592, 0.01796,
488 + 0.51514, 0.15768, 0.28475, 0.00577, 0.00076,
489 + 0.00000, 0.00578, 0.01695 /
490 DATA ((F(I,32,K),I=1,5),F(7,32,K),F(6,32,K),F(8,32,K),K=1,7)/
491 + 0.67325, 0.21095, 0.55633, 0.00375, 0.00000,
492 + 0.00000, 0.00420, 0.01922,
493 + 0.62664, 0.19307, 0.46726, 0.00403, 0.00000,
494 + 0.00000, 0.00444, 0.01763,
495 + 0.58778, 0.17849, 0.39852, 0.00414, 0.00000,
496 + 0.00000, 0.00451, 0.01629,
497 + 0.55458, 0.16637, 0.34417, 0.00414, 0.00000,
498 + 0.00000, 0.00449, 0.01513,
499 + 0.52586, 0.15603, 0.30062, 0.00408, 0.00014,
500 + 0.00000, 0.00443, 0.01412,
501 + 0.50060, 0.14709, 0.26506, 0.00399, 0.00039,
502 + 0.00000, 0.00433, 0.01324,
503 + 0.47803, 0.13922, 0.23556, 0.00390, 0.00060,
504 + 0.00004, 0.00422, 0.01243 /
505 DATA ((F(I,33,K),I=1,5),F(7,33,K),F(6,33,K),F(8,33,K),K=1,7)/
506 + 0.63725, 0.19005, 0.49256, 0.00254, 0.00000,
507 + 0.00000, 0.00286, 0.01421,
508 + 0.58912, 0.17279, 0.40616, 0.00275, 0.00000,
509 + 0.00000, 0.00308, 0.01299,
510 + 0.54935, 0.15884, 0.34118, 0.00283, 0.00000,
511 + 0.00000, 0.00317, 0.01190,
512 + 0.51567, 0.14726, 0.29090, 0.00283, 0.00000,
513 + 0.00000, 0.00317, 0.01104,
514 + 0.48672, 0.13748, 0.25128, 0.00280, 0.00007,
515 + 0.00000, 0.00313, 0.01025,
516 + 0.46139, 0.12906, 0.21942, 0.00275, 0.00026,
517 + 0.00000, 0.00307, 0.00957,
518 + 0.43891, 0.12168, 0.19331, 0.00267, 0.00040,
519 + 0.00001, 0.00298, 0.00902 /
520 DATA ((F(I,34,K),I=1,5),F(7,34,K),F(6,34,K),F(8,34,K),K=1,7)/
521 + 0.59662, 0.16877, 0.42790, 0.00132, 0.00000,
522 + 0.00000, 0.00206, 0.01033,
523 + 0.54786, 0.15249, 0.34752, 0.00159, 0.00000,
524 + 0.00000, 0.00222, 0.00936,
525 + 0.50791, 0.13943, 0.28809, 0.00172, 0.00000,
526 + 0.00000, 0.00228, 0.00856,
527 + 0.47435, 0.12867, 0.24287, 0.00177, 0.00000,
528 + 0.00000, 0.00226, 0.00790,
529 + 0.44564, 0.11963, 0.20769, 0.00179, 0.00006,
530 + 0.00000, 0.00221, 0.00731,
531 + 0.42073, 0.11187, 0.17973, 0.00177, 0.00019,
532 + 0.00000, 0.00216, 0.00682,
533 + 0.39872, 0.10510, 0.15704, 0.00175, 0.00030,
534 + 0.00003, 0.00210, 0.00636 /
535 DATA ((F(I,35,K),I=1,5),F(7,35,K),F(6,35,K),F(8,35,K),K=1,7)/
536 + 0.55241, 0.14825, 0.36867, 0.00060, 0.00000,
537 + 0.00000, 0.00141, 0.00731,
538 + 0.50389, 0.13309, 0.29445, 0.00085, 0.00000,
539 + 0.00000, 0.00154, 0.00662,
540 + 0.46440, 0.12102, 0.24076, 0.00100, 0.00000,
541 + 0.00000, 0.00158, 0.00601,
542 + 0.43149, 0.11114, 0.20061, 0.00108, 0.00000,
543 + 0.00000, 0.00156, 0.00553,
544 + 0.40355, 0.10285, 0.16982, 0.00109, 0.00005,
545 + 0.00000, 0.00153, 0.00513,
546 + 0.37940, 0.09579, 0.14560, 0.00109, 0.00015,
547 + 0.00000, 0.00150, 0.00476,
548 + 0.35819, 0.08968, 0.12614, 0.00109, 0.00022,
549 + 0.00004, 0.00145, 0.00442 /
550 DATA ((F(I,36,K),I=1,5),F(7,36,K),F(6,36,K),F(8,36,K),K=1,7)/
551 + 0.50593, 0.12897, 0.31457, 0.00039, 0.00000,
552 + 0.00000, 0.00086, 0.00508,
553 + 0.45836, 0.11502, 0.24747, 0.00057, 0.00000,
554 + 0.00000, 0.00097, 0.00457,
555 + 0.41999, 0.10394, 0.19955, 0.00063, 0.00000,
556 + 0.00000, 0.00101, 0.00418,
557 + 0.38820, 0.09495, 0.16422, 0.00069, 0.00000,
558 + 0.00000, 0.00103, 0.00380,
559 + 0.36140, 0.08747, 0.13748, 0.00069, 0.00003,
560 + 0.00000, 0.00103, 0.00351,
561 + 0.33834, 0.08112, 0.11669, 0.00067, 0.00011,
562 + 0.00000, 0.00101, 0.00325,
563 + 0.31817, 0.07563, 0.10021, 0.00068, 0.00016,
564 + 0.00003, 0.00099, 0.00300 /
565 DATA ((F(I,37,K),I=1,5),F(7,37,K),F(6,37,K),F(8,37,K),K=1,7)/
566 + 0.45796, 0.11074, 0.26561, 0.00026, 0.00000,
567 + 0.00000, 0.00049, 0.00348,
568 + 0.41212, 0.09810, 0.20521, 0.00035, 0.00000,
569 + 0.00000, 0.00057, 0.00314,
570 + 0.37540, 0.08816, 0.16317, 0.00045, 0.00000,
571 + 0.00000, 0.00063, 0.00283,
572 + 0.34528, 0.08014, 0.13264, 0.00046, 0.00000,
573 + 0.00000, 0.00065, 0.00255,
574 + 0.31994, 0.07349, 0.10991, 0.00047, 0.00000,
575 + 0.00000, 0.00065, 0.00236,
576 + 0.29823, 0.06789, 0.09246, 0.00049, 0.00004,
577 + 0.00000, 0.00063, 0.00217,
578 + 0.27940, 0.06308, 0.07877, 0.00049, 0.00008,
579 + -0.00002, 0.00062, 0.00199 /
580 DATA ((F(I,38,K),I=1,5),F(7,38,K),F(6,38,K),F(8,38,K),K=1,7)/
581 + 0.40933, 0.09346, 0.22084, -0.00006, 0.00000,
582 + 0.00000, 0.00036, 0.00231,
583 + 0.36580, 0.08224, 0.16749, 0.00007, 0.00000,
584 + 0.00000, 0.00043, 0.00205,
585 + 0.33123, 0.07348, 0.13118, 0.00012, 0.00000,
586 + 0.00000, 0.00046, 0.00186,
587 + 0.30297, 0.06647, 0.10540, 0.00020, 0.00000,
588 + 0.00000, 0.00046, 0.00166,
589 + 0.27945, 0.06068, 0.08641, 0.00020, 0.00002,
590 + 0.00000, 0.00045, 0.00154,
591 + 0.25938, 0.05583, 0.07204, 0.00024, 0.00006,
592 + 0.00000, 0.00044, 0.00140,
593 + 0.24204, 0.05168, 0.06086, 0.00025, 0.00009,
594 + 0.00002, 0.00043, 0.00128 /
595 DATA ((F(I,39,K),I=1,5),F(7,39,K),F(6,39,K),F(8,39,K),K=1,7)/
596 + 0.36094, 0.07776, 0.17814, -0.00018, 0.00000,
597 + 0.00000, 0.00027, 0.00147,
598 + 0.32041, 0.06800, 0.13388, -0.00008, 0.00000,
599 + 0.00000, 0.00030, 0.00129,
600 + 0.28840, 0.06043, 0.10374, 0.00001, 0.00000,
601 + 0.00000, 0.00031, 0.00114,
602 + 0.26245, 0.05438, 0.08246, 0.00004, 0.00001,
603 + 0.00000, 0.00031, 0.00106,
604 + 0.24089, 0.04944, 0.06697, 0.00008, 0.00004,
605 + 0.00000, 0.00030, 0.00095,
606 + 0.22264, 0.04529, 0.05530, 0.00010, 0.00005,
607 + 0.00000, 0.00029, 0.00088,
608 + 0.20691, 0.04176, 0.04632, 0.00013, 0.00006,
609 + 0.00004, 0.00028, 0.00080 /
610 DATA ((F(I,40,K),I=1,5),F(7,40,K),F(6,40,K),F(8,40,K),K=1,7)/
611 + 0.31447, 0.06406, 0.14592, -0.00003, 0.00000,
612 + 0.00000, 0.00008, 0.00093,
613 + 0.27711, 0.05560, 0.10695, 0.00001, 0.00000,
614 + 0.00000, 0.00012, 0.00081,
615 + 0.24784, 0.04911, 0.08140, 0.00005, 0.00000,
616 + 0.00000, 0.00014, 0.00072,
617 + 0.22428, 0.04394, 0.06381, 0.00007, 0.00000,
618 + 0.00000, 0.00016, 0.00065,
619 + 0.20488, 0.03975, 0.05122, 0.00008, 0.00000,
620 + 0.00000, 0.00016, 0.00059,
621 + 0.18855, 0.03625, 0.04190, 0.00009, 0.00001,
622 + 0.00000, 0.00016, 0.00054,
623 + 0.17455, 0.03329, 0.03478, 0.00009, 0.00002,
624 + 0.00000, 0.00015, 0.00050 /
625 DATA ((F(I,41,K),I=1,5),F(7,41,K),F(6,41,K),F(8,41,K),K=1,7)/
626 + 0.26913, 0.05145, 0.11355, -0.00016, 0.00000,
627 + 0.00000, 0.00013, 0.00052,
628 + 0.23552, 0.04437, 0.08237, -0.00011, 0.00000,
629 + 0.00000, 0.00013, 0.00046,
630 + 0.20936, 0.03897, 0.06203, -0.00006, 0.00000,
631 + 0.00000, 0.00012, 0.00042,
632 + 0.18845, 0.03470, 0.04812, -0.00003, 0.00001,
633 + 0.00000, 0.00013, 0.00038,
634 + 0.17129, 0.03123, 0.03825, 0.00000, 0.00002,
635 + 0.00000, 0.00012, 0.00034,
636 + 0.15693, 0.02837, 0.03100, 0.00000, 0.00003,
637 + 0.00000, 0.00012, 0.00032,
638 + 0.14466, 0.02595, 0.02551, 0.00002, 0.00004,
639 + 0.00003, 0.00011, 0.00029 /
640 DATA ((F(I,42,K),I=1,5),F(7,42,K),F(6,42,K),F(8,42,K),K=1,7)/
641 + 0.18815, 0.03174, 0.06747, 0.00003, 0.00000,
642 + 0.00000, 0.00001, 0.00017,
643 + 0.16218, 0.02694, 0.04718, 0.00000, 0.00000,
644 + 0.00000, 0.00002, 0.00017,
645 + 0.14226, 0.02335, 0.03442, 0.00000, 0.00000,
646 + 0.00000, 0.00002, 0.00014,
647 + 0.12656, 0.02054, 0.02599, 0.00000, 0.00000,
648 + 0.00000, 0.00003, 0.00013,
649 + 0.11382, 0.01830, 0.02017, 0.00002, 0.00000,
650 + 0.00000, 0.00003, 0.00011,
651 + 0.10327, 0.01646, 0.01600, 0.00002, 0.00000,
652 + 0.00000, 0.00003, 0.00010,
653 + 0.09437, 0.01492, 0.01292, 0.00001, 0.00001,
654 + 0.00000, 0.00003, 0.00009 /
655 DATA ((F(I,43,K),I=1,5),F(7,43,K),F(6,43,K),F(8,43,K),K=1,7)/
656 + 0.12097, 0.01770, 0.03606, 0.00003, 0.00000,
657 + 0.00000, -0.00001, 0.00005,
658 + 0.10255, 0.01479, 0.02425, 0.00001, 0.00000,
659 + 0.00000, 0.00000, 0.00005,
660 + 0.08869, 0.01263, 0.01712, 0.00001, 0.00000,
661 + 0.00000, 0.00000, 0.00004,
662 + 0.07790, 0.01098, 0.01257, 0.00001, 0.00000,
663 + 0.00000, 0.00000, 0.00003,
664 + 0.06927, 0.00967, 0.00952, 0.00001, 0.00000,
665 + 0.00000, 0.00000, 0.00003,
666 + 0.06220, 0.00861, 0.00739, 0.00002, 0.00000,
667 + 0.00000, 0.00000, 0.00002,
668 + 0.05631, 0.00773, 0.00585, 0.00001, 0.00000,
669 + 0.00000, 0.00000, 0.00002 /
670 DATA ((F(I,44,K),I=1,5),F(7,44,K),F(6,44,K),F(8,44,K),K=1,7)/
671 + 0.06972, 0.00867, 0.01668, 0.00001, 0.00000,
672 + 0.00000, -0.00001, 0.00002,
673 + 0.05802, 0.00712, 0.01080, 0.00001, 0.00000,
674 + 0.00000, 0.00000, 0.00002,
675 + 0.04939, 0.00599, 0.00737, 0.00001, 0.00000,
676 + 0.00000, 0.00000, 0.00001,
677 + 0.04278, 0.00513, 0.00525, 0.00001, 0.00000,
678 + 0.00000, -0.00001, 0.00001,
679 + 0.03757, 0.00446, 0.00387, 0.00001, 0.00000,
680 + 0.00000, 0.00000, 0.00001,
681 + 0.03336, 0.00393, 0.00294, 0.00000, 0.00000,
682 + 0.00000, 0.00000, 0.00001,
683 + 0.02987, 0.00349, 0.00227, 0.00000, 0.00000,
684 + 0.00000, 0.00000, 0.00001 /
685 DATA ((F(I,45,K),I=1,5),F(7,45,K),F(6,45,K),F(8,45,K),K=1,7)/
686 + 0.03421, 0.00350, 0.00626, -0.00001, 0.00001,
687 + 0.00000, 0.00001, 0.00000,
688 + 0.02789, 0.00281, 0.00388, 0.00000, 0.00000,
689 + 0.00000, 0.00000, 0.00000,
690 + 0.02332, 0.00233, 0.00255, 0.00000, 0.00000,
691 + 0.00000, 0.00000, 0.00000,
692 + 0.01987, 0.00196, 0.00176, 0.00000, 0.00000,
693 + 0.00000, 0.00000, 0.00000,
694 + 0.01721, 0.00169, 0.00126, 0.00000, 0.00000,
695 + 0.00000, 0.00000, 0.00000,
696 + 0.01508, 0.00146, 0.00093, 0.00000, 0.00000,
697 + 0.00000, 0.00000, 0.00000,
698 + 0.01334, 0.00129, 0.00070, 0.00000, 0.00000,
699 + 0.00000, 0.00000, 0.00000 /
700 DATA ((F(I,46,K),I=1,5),F(7,46,K),F(6,46,K),F(8,46,K),K=1,7)/
701 + 0.00335, 0.00020, 0.00037, 0.00000, 0.00000,
702 + 0.00000, 0.00000, 0.00000,
703 + 0.00257, 0.00015, 0.00017, 0.00000, 0.00000,
704 + 0.00000, 0.00000, 0.00000,
705 + 0.00203, 0.00012, 0.00010, 0.00000, 0.00000,
706 + 0.00000, 0.00000, 0.00000,
707 + 0.00165, 0.00009, 0.00006, 0.00000, 0.00000,
708 + 0.00000, 0.00000, 0.00000,
709 + 0.00137, 0.00008, 0.00004, 0.00000, 0.00000,
710 + 0.00000, 0.00000, 0.00000,
711 + 0.00116, 0.00006, 0.00003, 0.00000, 0.00000,
712 + 0.00000, 0.00000, 0.00000,
713 + 0.00099, 0.00005, 0.00002, 0.00000, 0.00000,
714 + 0.00000, 0.00000, 0.00000 /
715 DATA XX/1.d-5,2.d-5,4.d-5,6.d-5,8.d-5,
716 . 1.D-4,2.D-4,4.D-4,6.D-4,8.D-4,
717 . 1.D-3,2.D-3,4.D-3,6.D-3,8.D-3,
718 . 1.D-2,2.D-2,4.D-2,6.D-2,8.D-2,
719 . .1D0,.125D0,.15D0,.175D0,.2D0,.225D0,.25D0,.275D0,
720 . .3D0,.325D0,.35D0,.375D0,.4D0,.425D0,.45D0,.475D0,
721 . .5D0,.525D0,.55D0,.575D0,.6D0,.65D0,.7D0,.75D0,
722 . .8D0,.9D0,1.D0/
723 DATA XMIN,XMAX,QSQMIN,QSQMAX/1.D-5,1.D0,0.625D0,5.D0/
724 DATA N0/2,5,5,9,0,0,9,9/
725 DATA ZEROD/0.D0/,PONED/0.1D0/,ONED/1.D0/,ONEDO/1.1D0/,TWOD/2.D0/
726 DATA INIT/0/
727C
728 xsave=x
729C
730 IF(INIT.NE.0) GOTO 10
731 INIT=1
732 DO 20 N=1,nx-1
733 DO 20 M=1,7
734C READ(34,50)F(1,N,M),F(2,N,M),F(3,N,M),F(4,N,M),F(5,N,M),F(7,N,M),
735C . F(6,N,M),F(8,N,M)
736C 1=UV 2=DV 3=GLUE 4=UBAR 5=CBAR 7=BBAR 6=SBAR 8=DBAR
737 DO 25 I=1,8
738 25 F(I,N,M)=F(I,N,M)/(ONED-XX(N))**N0(I)
739 20 CONTINUE
740 DO 31 J=1,NTENTH-1
741 XX(J)= LOG10(XX(J))+ONEDO
742 DO 31 I=1,8
743 IF(I.EQ.7 .or. I.EQ.5) GO TO 31
744 DO 30 K=1,7
745 30 F(I,J,K)= LOG(F(I,J,K))*F(I,NTENTH,K)/ LOG(F(I,NTENTH,K))
746 31 CONTINUE
747 50 FORMAT(8F10.5)
748 DO 40 I=1,8
749 DO 40 M=1,7
750 40 F(I,NX,M)=ZEROD
751 10 CONTINUE
752 IF(X.LT.XMIN) X=XMIN
753 IF(X.GT.XMAX) X=XMAX
754 QSQ=SCALE**2
755 IF(QSQ.LT.QSQMIN) QSQ=QSQMIN
756 IF(QSQ.GT.QSQMAX) QSQ=QSQMAX
757 XXX=X
758 IF(X.LT.PONED) XXX= LOG10(X)+ONEDO
759 N=0
760 70 N=N+1
761 IF(XXX.GT.XX(N+1)) GOTO 70
762 A=(XXX-XX(N))/(XX(N+1)-XX(N))
763 RM= LOG(QSQ/QSQMIN)/ LOG(TWOD)*TWOD
764 B=RM-AINT(RM)
765 M=1+ INT(RM)
766 DO 60 I=1,8
767 G(I)= (ONED-A)*(ONED-B)*F(I,N,M)+(ONED-A)*B*F(I,N,M+1)
768 . + A*(ONED-B)*F(I,N+1,M) + A*B*F(I,N+1,M+1)
769 IF(N.GE.NTENTH) GOTO 65
770 IF(I.EQ.7 .OR. I.EQ.5) GOTO 65
771 FAC=(1.D0-B)*F(I,NTENTH,M)+B*F(I,NTENTH,M+1)
772 G(I)=FAC**(G(I)/FAC)
773 65 CONTINUE
774 G(I)=G(I)*(ONED-X)**N0(I)
775 60 CONTINUE
776 UPV=G(1)
777 DNV=G(2)
778 USEA=G(4)
779 DSEA=G(8)
780 STR=G(6)
781 CHM=G(5)
782 GLU=G(3)
783 BOT=G(7)
784C
785 x=xsave
786C
787 RETURN
788 END