PDF/spdf/sfmsqa.F

   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)
  17 C
  18 C     THIS IS THE NEW  "A" FIT -- May 1994 -- low Q^2 range
  19 C
  20 #include "pdf/impdp.inc"
  21 C
  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/
 727 C
 728       xsave=x
 729 C
 730       IF(INIT.NE.0) GOTO 10
 731       INIT=1
 732       DO 20 N=1,nx-1
 733       DO 20 M=1,7
 734 C     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),
 735 C    .  F(6,N,M),F(8,N,M)
 736 C 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)
 784 C
 785       x=xsave
 786 C
 787       RETURN
 788       END