#include "pdf/pilot.h" SUBROUTINE SFDSQA(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,GLU) C THIS IS THE NEW "A" FIT -- May 1994 -- low Q^2 range C DIS DIS DIS DIS DIS DIS DIS DIS DIS DIS DIS DIS DIS DIS C #include "pdf/impdp.inc" C PARAMETER(NX=47) PARAMETER(NTENTH=21) DIMENSION F(8,NX,8),G(8),XX(NX),N0(8) DATA ((F(I,1,K),I=1,5),F(7,1,K),F(6,1,K),F(8,1,K),K=1,7)/ + .00353, .00627, 4.23812, 1.87623, .00000, + .00000, .51265, 1.87021, + .00377, .00645, 6.42705, 1.94786, .00000, + .00000, .61236, 1.94247, + .00399, .00660, 9.43898, 2.04451, .00000, + .00000, .73330, 2.03973, + .00420, .00675, 13.16286, 2.16589, .00000, + .00000, .87600, 2.16169, + .00440, .00688, 17.54780, 2.30621, .00000, + .00000, 1.03533, 2.30254, + .00460, .00702, 22.55753, 2.46430, .06825, + .00000, 1.21072, 2.46116, + .00477, .00714, 28.12229, 2.63824, .27055, + .00000, 1.40039, 2.63556 / DATA ((F(I,2,K),I=1,5),F(7,2,K),F(6,2,K),F(8,2,K),K=1,7)/ + .00538, .00802, 3.22632, 1.51244, .00000, + .00000, .40761, 1.50993, + .00571, .00826, 5.15995, 1.57471, .00000, + .00000, .49238, 1.57262, + .00602, .00846, 7.71045, 1.65616, .00000, + .00000, .59325, 1.65446, + .00632, .00866, 10.78281, 1.75593, .00000, + .00000, .71002, 1.75459, + .00659, .00884, 14.34431, 1.86953, .00000, + .00000, .83878, 1.86852, + .00686, .00902, 18.35565, 1.99595, .05595, + .00000, .97896, 1.99526, + .00712, .00919, 22.75467, 2.13362, .21642, + .00000, 1.12915, 2.13323 / DATA ((F(I,3,K),I=1,5),F(7,3,K),F(6,3,K),F(8,3,K),K=1,7)/ + .00806, .01027, 2.45661, 1.21901, .00000, + .00000, .32420, 1.21891, + .00853, .01059, 4.14420, 1.27292, .00000, + .00000, .39616, 1.27308, + .00896, .01087, 6.30016, 1.34140, .00000, + .00000, .48020, 1.34180, + .00937, .01113, 8.83462, 1.42328, .00000, + .00000, .57568, 1.42390, + .00975, .01137, 11.72666, 1.51512, .00000, + .00000, .67963, 1.51595, + .01013, .01162, 14.93678, 1.61607, .04591, + .00000, .79155, 1.61709, + .01048, .01184, 18.41091, 1.72486, .17321, + .00000, .91032, 1.72606 / DATA ((F(I,4,K),I=1,5),F(7,4,K),F(6,4,K),F(8,4,K),K=1,7)/ + .01014, .01188, 2.09495, 1.07443, .00000, + .00000, .28369, 1.07538, + .01071, .01225, 3.64657, 1.12390, .00000, + .00000, .34904, 1.12504, + .01125, .01259, 5.59929, 1.18568, .00000, + .00000, .42450, 1.18701, + .01174, .01289, 7.86366, 1.25856, .00000, + .00000, .50933, 1.26006, + .01220, .01319, 10.42364, 1.33960, .00000, + .00000, .60099, 1.34125, + .01266, .01348, 13.24045, 1.42800, .04089, + .00000, .69902, 1.42979, + .01308, .01374, 16.26478, 1.52271, .15196, + .00000, .80246, 1.52463 / DATA ((F(I,5,K),I=1,5),F(7,5,K),F(6,5,K),F(8,5,K),K=1,7)/ + .01193, .01319, 1.87136, .98235, .00000, + .00000, .25813, .98392, + .01258, .01359, 3.33089, 1.02884, .00000, + .00000, .31915, 1.03058, + .01318, .01397, 5.15064, 1.08626, .00000, + .00000, .38906, 1.08815, + .01374, .01433, 7.24090, 1.15331, .00000, + .00000, .46701, 1.15535, + .01428, .01465, 9.58839, 1.22740, .00000, + .00000, .55080, 1.22955, + .01479, .01496, 12.15519, 1.30782, .03767, + .00000, .63999, 1.31010, + .01529, .01527, 14.89528, 1.39358, .13842, + .00000, .73373, 1.39596 / DATA ((F(I,6,K),I=1,5),F(7,6,K),F(6,6,K),F(8,6,K),K=1,7)/ + .01351, .01430, 1.71468, .91640, .00000, + .00000, .23995, .91839, + .01423, .01475, 3.10549, .96070, .00000, + .00000, .29779, .96283, + .01490, .01516, 4.82815, 1.01490, .00000, + .00000, .36364, 1.01717, + .01552, .01555, 6.79254, 1.07773, .00000, + .00000, .43666, 1.08013, + .01611, .01591, 8.98726, 1.14683, .00000, + .00000, .51480, 1.14934, + .01668, .01626, 11.37514, 1.22153, .03529, + .00000, .59766, 1.22414, + .01724, .01658, 13.91271, 1.30091, .12850, + .00000, .68447, 1.30361 / DATA ((F(I,7,K),I=1,5),F(7,7,K),F(6,7,K),F(8,7,K),K=1,7)/ + .01973, .01843, 1.30792, .73836, .00000, + .00000, .19165, .74150, + .02071, .01901, 2.50161, .77638, .00000, + .00000, .24064, .77965, + .02164, .01956, 3.95367, .82157, .00000, + .00000, .29525, .82495, + .02250, .02006, 5.57292, .87276, .00000, + .00000, .35466, .87624, + .02331, .02054, 7.35237, .92821, .00000, + .00000, .41740, .93177, + .02411, .02099, 9.25794, .98738, .02895, + .00000, .48318, .99104, + .02486, .02143, 11.25363, 1.04960, .10236, + .00000, .55137, 1.05333 / DATA ((F(I,8,K),I=1,5),F(7,8,K),F(6,8,K),F(8,8,K),K=1,7)/ + .02860, .02387, .99986, .59496, .00000, + .00000, .15390, .59922, + .02994, .02463, 2.02268, .62745, .00000, + .00000, .19534, .63186, + .03120, .02534, 3.24642, .66488, .00000, + .00000, .24044, .66941, + .03238, .02600, 4.58011, .70627, .00000, + .00000, .28853, .71093, + .03349, .02660, 6.02017, .75042, .00000, + .00000, .33861, .75519, + .03458, .02721, 7.53631, .79695, .02377, + .00000, .39051, .80182, + .03562, .02777, 9.09952, .84533, .08133, + .00000, .44376, .85029 / DATA ((F(I,9,K),I=1,5),F(7,9,K),F(6,9,K),F(8,9,K),K=1,7)/ + .03541, .02785, .85607, .52448, .00000, + .00000, .13594, .52948, + .03702, .02872, 1.79177, .55403, .00000, + .00000, .17343, .55918, + .03853, .02954, 2.89947, .58737, .00000, + .00000, .21363, .59270, + .03994, .03030, 4.08932, .62376, .00000, + .00000, .25595, .62921, + .04128, .03102, 5.35969, .66221, .00000, + .00000, .29970, .66779, + .04259, .03172, 6.68285, .70242, .02116, + .00000, .34469, .70814, + .04384, .03239, 8.03367, .74396, .07087, + .00000, .39058, .74981 / DATA ((F(I,10,K),I=1,5),F(7,10,K),F(6,10,K),F(8,10,K),K=1,7)/ + .04114, .03109, .76776, .47972, .00000, + .00000, .12483, .48529, + .04297, .03205, 1.64767, .50719, .00000, + .00000, .15966, .51298, + .04468, .03295, 2.68028, .53784, .00000, + .00000, .19662, .54379, + .04629, .03380, 3.77709, .57092, .00000, + .00000, .23521, .57706, + .04781, .03460, 4.93803, .60569, .00000, + .00000, .27487, .61197, + .04931, .03538, 6.13732, .64187, .01945, + .00000, .31546, .64830, + .05072, .03614, 7.35252, .67905, .06418, + .00000, .35667, .68564 / DATA ((F(I,11,K),I=1,5),F(7,11,K),F(6,11,K),F(8,11,K),K=1,7)/ + .04628, .03389, .70475, .44767, .00000, + .00000, .11695, .45366, + .04828, .03494, 1.54754, .47364, .00000, + .00000, .14979, .47983, + .05016, .03594, 2.52696, .50227, .00000, + .00000, .18433, .50868, + .05194, .03687, 3.55628, .53299, .00000, + .00000, .22019, .53956, + .05361, .03774, 4.63730, .56508, .00000, + .00000, .25686, .57183, + .05525, .03861, 5.74604, .59832, .01817, + .00000, .29424, .60523, + .05683, .03941, 6.86227, .63238, .05918, + .00000, .33205, .63944 / DATA ((F(I,12,K),I=1,5),F(7,12,K),F(6,12,K),F(8,12,K),K=1,7)/ + .06634, .04460, .54955, .36164, .00000, + .00000, .09697, .36963, + .06902, .04597, 1.28832, .38302, .00000, + .00000, .12399, .39128, + .07156, .04728, 2.11904, .40586, .00000, + .00000, .15169, .41438, + .07395, .04850, 2.96136, .42982, .00000, + .00000, .17986, .43859, + .07622, .04966, 3.82264, .45449, .00000, + .00000, .20829, .46348, + .07843, .05079, 4.68465, .47971, .01465, + .00000, .23691, .48893, + .08055, .05186, 5.53362, .50526, .04585, + .00000, .26555, .51468 / DATA ((F(I,13,K),I=1,5),F(7,13,K),F(6,13,K),F(8,13,K),K=1,7)/ + .09459, .05933, .45006, .29360, .00000, + .00000, .08276, .30437, + .09814, .06112, 1.10456, .31050, .00000, + .00000, .10425, .32162, + .10155, .06285, 1.80790, .32807, .00000, + .00000, .12580, .33952, + .10473, .06446, 2.49068, .34620, .00000, + .00000, .14739, .35798, + .10777, .06600, 3.16635, .36464, .00000, + .00000, .16891, .37671, + .11072, .06748, 3.82305, .38326, .01152, + .00000, .19031, .39562, + .11353, .06890, 4.45308, .40190, .03481, + .00000, .21149, .41453 / DATA ((F(I,14,K),I=1,5),F(7,14,K),F(6,14,K),F(8,14,K),K=1,7)/ + .11598, .07030, .42083, .26052, .00000, + .00000, .07654, .27338, + .12018, .07245, 1.03073, .27483, .00000, + .00000, .09494, .28811, + .12420, .07450, 1.66542, .28956, .00000, + .00000, .11322, .30324, + .12798, .07641, 2.26414, .30469, .00000, + .00000, .13143, .31874, + .13156, .07822, 2.84358, .31997, .00000, + .00000, .14946, .33439, + .13503, .07996, 3.39582, .33531, .00979, + .00000, .16730, .35006, + .13832, .08161, 3.91637, .35057, .02918, + -.00001, .18482, .36564 / DATA ((F(I,15,K),I=1,5),F(7,15,K),F(6,15,K),F(8,15,K),K=1,7)/ + .13377, .07939, .41558, .23959, .00000, + .00000, .07283, .25410, + .13849, .08184, .99416, .25213, .00000, + .00000, .08912, .26710, + .14304, .08416, 1.58125, .26497, .00000, + .00000, .10522, .28042, + .14730, .08633, 2.12318, .27815, .00000, + .00000, .12124, .29401, + .15132, .08837, 2.63867, .29140, .00000, + .00000, .13703, .30768, + .15519, .09033, 3.12247, .30466, .00857, + .00000, .15259, .32132, + .15886, .09215, 3.57229, .31776, .02549, + -.00001, .16780, .33479 / DATA ((F(I,16,K),I=1,5),F(7,16,K),F(6,16,K),F(8,16,K),K=1,7)/ + .14903, .08698, .42079, .22390, .00000, + .00000, .06996, .23937, + .15428, .08974, .97466, .23519, .00000, + .00000, .08469, .25125, + .15929, .09235, 1.52498, .24674, .00000, + .00000, .09927, .26336, + .16398, .09477, 2.02410, .25857, .00000, + .00000, .11375, .27570, + .16840, .09703, 2.49209, .27045, .00000, + .00000, .12799, .28805, + .17264, .09917, 2.92573, .28225, .00757, + .00000, .14197, .30029, + .17663, .10116, 3.32426, .29388, .02270, + .00000, .15557, .31233 / DATA ((F(I,17,K),I=1,5),F(7,17,K),F(6,17,K),F(8,17,K),K=1,7)/ + .20755, .11696, .49179, .18309, .00000, + .00000, .06331, .20391, + .21483, .12090, .96223, .19043, .00000, + .00000, .07349, .21220, + .22167, .12450, 1.39560, .19812, .00000, + .00000, .08371, .22074, + .22795, .12773, 1.76565, .20605, .00000, + .00000, .09388, .22941, + .23378, .13067, 2.09508, .21395, .00000, + .00000, .10382, .23798, + .23925, .13335, 2.38648, .22169, .00451, + .00000, .11345, .24633, + .24430, .13580, 2.64317, .22922, .01502, + .00000, .12270, .25439 / DATA ((F(I,18,K),I=1,5),F(7,18,K),F(6,18,K),F(8,18,K),K=1,7)/ + .28677, .15844, .65201, .14832, .00000, + .00000, .05678, .17806, + .29713, .16379, 1.01515, .15204, .00000, + .00000, .06325, .18312, + .30643, .16834, 1.31959, .15631, .00000, + .00000, .06998, .18849, + .31466, .17220, 1.56052, .16082, .00000, + .00000, .07671, .19391, + .32205, .17554, 1.76050, .16531, .00000, + .00000, .08320, .19917, + .32867, .17842, 1.92619, .16961, .00128, + .00000, .08935, .20410, + .33465, .18092, 2.06321, .17369, .00833, + -.00001, .09512, .20870 / DATA ((F(I,19,K),I=1,5),F(7,19,K),F(6,19,K),F(8,19,K),K=1,7)/ + .34542, .18867, .77868, .12710, .00000, + .00000, .05139, .16412, + .35819, .19454, 1.06122, .12899, .00000, + .00000, .05624, .16744, + .36909, .19926, 1.28184, .13152, .00000, + .00000, .06139, .17101, + .37837, .20302, 1.44554, .13426, .00000, + .00000, .06644, .17453, + .38641, .20610, 1.57317, .13698, .00000, + .00000, .07123, .17784, + .39332, .20854, 1.67235, .13951, .00001, + .00000, .07566, .18081, + .39934, .21054, 1.74894, .14183, .00482, + .00000, .07971, .18346 / DATA ((F(I,20,K),I=1,5),F(7,20,K),F(6,20,K),F(8,20,K),K=1,7)/ + .39417, .21269, .87132, .11000, .00000, + .00000, .04597, .15304, + .40857, .21843, 1.08735, .11080, .00000, + .00000, .04998, .15496, + .42039, .22263, 1.24495, .11220, .00000, + .00000, .05416, .15715, + .42999, .22574, 1.35369, .11380, .00000, + .00000, .05818, .15920, + .43796, .22808, 1.43230, .11540, .00000, + .00000, .06191, .16105, + .44452, .22971, 1.48820, .11679, .00001, + .00000, .06523, .16258, + .44997, .23090, 1.52687, .11804, .00257, + -.00001, .06820, .16384 / DATA ((F(I,21,K),I=1,5),F(7,21,K),F(6,21,K),F(8,21,K),K=1,7)/ + .43740, .23195, .93597, .09467, .00000, + .00000, .04051, .14281, + .45273, .23699, 1.09539, .09482, .00000, + .00000, .04393, .14355, + .46465, .24034, 1.20287, .09550, .00000, + .00000, .04744, .14447, + .47385, .24250, 1.26983, .09636, .00000, + .00000, .05068, .14526, + .48116, .24382, 1.31262, .09719, .00000, + .00000, .05359, .14591, + .48675, .24447, 1.33802, .09787, .00001, + .00000, .05612, .14629, + .49116, .24467, 1.35083, .09841, .00096, + .00000, .05830, .14647 / DATA ((F(I,22,K),I=1,5),F(7,22,K),F(6,22,K),F(8,22,K),K=1,7)/ + .48495, .25086, .98361, .07726, .00000, + .00000, .03390, .12954, + .50004, .25430, 1.08372, .07710, .00000, + .00000, .03696, .12905, + .51094, .25600, 1.14164, .07738, .00000, + .00000, .03990, .12871, + .51874, .25651, 1.16898, .07772, .00000, + .00000, .04250, .12835, + .52440, .25628, 1.17918, .07799, .00000, + .00000, .04472, .12789, + .52821, .25543, 1.17779, .07813, .00001, + .00000, .04656, .12725, + .53076, .25425, 1.16862, .07816, .00001, + -.00001, .04809, .12650 / DATA ((F(I,23,K),I=1,5),F(7,23,K),F(6,23,K),F(8,23,K),K=1,7)/ + .52685, .26533, 1.00111, .06198, .00000, + .00000, .02725, .11577, + .54044, .26638, 1.05264, .06174, .00000, + .00000, .03024, .11435, + .54931, .26589, 1.07218, .06180, .00000, + .00000, .03289, .11316, + .55485, .26445, 1.07015, .06185, .00000, + .00000, .03506, .11194, + .55816, .26246, 1.05720, .06185, .00000, + .00000, .03682, .11070, + .55957, .26006, 1.03756, .06171, .00001, + .00000, .03820, .10935, + .55978, .25746, 1.01406, .06150, .00001, + -.00003, .03931, .10794 / DATA ((F(I,24,K),I=1,5),F(7,24,K),F(6,24,K),F(8,24,K),K=1,7)/ + .56454, .27445, .99514, .04789, .00000, + .00000, .02161, .10191, + .57507, .27304, 1.00701, .04803, .00000, + .00000, .02440, .09988, + .58075, .27034, .99689, .04819, .00000, + .00000, .02674, .09801, + .58312, .26698, .97320, .04827, .00000, + .00000, .02854, .09624, + .58342, .26328, .94399, .04821, .00000, + .00000, .02995, .09448, + .58201, .25935, .91213, .04801, .00001, + .00000, .03101, .09271, + .57958, .25541, .87953, .04770, .00001, + -.00001, .03178, .09099 / DATA ((F(I,25,K),I=1,5),F(7,25,K),F(6,25,K),F(8,25,K),K=1,7)/ + .59691, .27992, .97091, .03603, .00000, + .00000, .01614, .08820, + .60334, .27563, .95082, .03651, .00000, + .00000, .01898, .08577, + .60506, .27058, .91809, .03687, .00000, + .00000, .02114, .08357, + .60374, .26525, .87878, .03699, .00000, + .00000, .02277, .08146, + .60062, .25988, .83847, .03696, .00000, + .00000, .02394, .07943, + .59609, .25451, .79881, .03678, .00001, + .00000, .02476, .07747, + .59086, .24929, .76083, .03650, .00001, + -.00002, .02533, .07557 / DATA ((F(I,26,K),I=1,5),F(7,26,K),F(6,26,K),F(8,26,K),K=1,7)/ + .62371, .28203, .93107, .02649, .00000, + .00000, .01124, .07481, + .62519, .27467, .88680, .02720, .00000, + .00000, .01413, .07232, + .62239, .26718, .83766, .02763, .00000, + .00000, .01625, .07001, + .61707, .25983, .78764, .02781, .00000, + .00000, .01772, .06790, + .61031, .25286, .74029, .02785, .00000, + .00000, .01879, .06577, + .60253, .24619, .69616, .02770, .00001, + .00000, .01949, .06376, + .59439, .23988, .65555, .02744, .00001, + -.00002, .01993, .06188 / DATA ((F(I,27,K),I=1,5),F(7,27,K),F(6,27,K),F(8,27,K),K=1,7)/ + .64685, .27834, .87914, .01677, .00000, + .00000, .00808, .06246, + .64193, .26874, .81691, .01836, .00000, + .00000, .01072, .06006, + .63371, .25949, .75668, .01935, .00000, + .00000, .01258, .05778, + .62377, .25070, .70017, .01991, .00000, + .00000, .01381, .05566, + .61311, .24249, .64913, .02015, .00000, + .00000, .01465, .05364, + .60199, .23483, .60318, .02021, .00001, + .00000, .01518, .05174, + .59088, .22773, .56203, .02013, .00001, + .00002, .01554, .04995 / DATA ((F(I,28,K),I=1,5),F(7,28,K),F(6,28,K),F(8,28,K),K=1,7)/ + .66107, .27515, .82378, .01211, .00000, + .00000, .00374, .05101, + .65020, .26245, .74720, .01349, .00000, + .00000, .00665, .04889, + .63701, .25088, .67896, .01434, .00000, + .00000, .00861, .04683, + .62284, .24040, .61850, .01480, .00000, + .00000, .00991, .04486, + .60857, .23090, .56585, .01502, .00000, + .00000, .01076, .04303, + .59439, .22225, .51975, .01504, .00001, + .00000, .01129, .04132, + .58066, .21436, .47938, .01495, .00001, + -.00014, .01159, .03972 / DATA ((F(I,29,K),I=1,5),F(7,29,K),F(6,29,K),F(8,29,K),K=1,7)/ + .67179, .26585, .75762, .00606, .00000, + .00000, .00178, .04082, + .65387, .25121, .67378, .00798, .00000, + .00000, .00448, .03902, + .63509, .23825, .60220, .00918, .00000, + .00000, .00622, .03725, + .61643, .22673, .54096, .00988, .00000, + .00000, .00733, .03555, + .59846, .21645, .48891, .01024, .00000, + .00000, .00805, .03397, + .58120, .20725, .44429, .01048, .00001, + .00000, .00851, .03240, + .56495, .19892, .40590, .01053, .00001, + -.00002, .00877, .03101 / DATA ((F(I,30,K),I=1,5),F(7,30,K),F(6,30,K),F(8,30,K),K=1,7)/ + .67527, .25529, .69187, .00245, .00000, + .00000, .00001, .03191, + .65063, .23870, .60286, .00452, .00000, + .00000, .00236, .03056, + .62661, .22441, .52998, .00578, .00000, + .00000, .00405, .02908, + .60372, .21204, .46960, .00656, .00000, + .00000, .00510, .02766, + .58244, .20113, .41941, .00701, .00000, + .00000, .00575, .02631, + .56258, .19146, .37721, .00720, .00001, + .00000, .00615, .02511, + .54411, .18290, .34143, .00732, .00001, + -.00004, .00641, .02390 / DATA ((F(I,31,K),I=1,5),F(7,31,K),F(6,31,K),F(8,31,K),K=1,7)/ + .67251, .24277, .62540, .00010, .00000, + .00000, .00001, .02445, + .64126, .22455, .53420, .00214, .00000, + .00000, .00078, .02338, + .61229, .20924, .46215, .00332, .00000, + .00000, .00236, .02228, + .58571, .19616, .40404, .00407, .00000, + .00000, .00338, .02110, + .56145, .18486, .35676, .00452, .00000, + .00000, .00399, .02008, + .53927, .17498, .31763, .00477, .00001, + .00000, .00440, .01905, + .51898, .16628, .28494, .00492, .00001, + .00000, .00462, .01810 / DATA ((F(I,32,K),I=1,5),F(7,32,K),F(6,32,K),F(8,32,K),K=1,7)/ + .66391, .22706, .55633, .00001, .00000, + .00000, .00001, .01820, + .62628, .20811, .46726, .00001, .00000, + .00000, .00003, .01752, + .59273, .19238, .39852, .00116, .00000, + .00000, .00139, .01675, + .56270, .17922, .34417, .00202, .00000, + .00000, .00225, .01585, + .53591, .16787, .30062, .00253, .00000, + .00000, .00280, .01501, + .51187, .15809, .26506, .00284, .00001, + .00000, .00311, .01421, + .49013, .14952, .23573, .00306, .00001, + .00004, .00331, .01340 / DATA ((F(I,33,K),I=1,5),F(7,33,K),F(6,33,K),F(8,33,K),K=1,7)/ + .64879, .21198, .49256, .00001, .00000, + .00000, .00001, .01307, + .60557, .19218, .40616, .00001, .00000, + .00000, .00001, .01281, + .56813, .17607, .34118, .00036, .00000, + .00000, .00038, .01218, + .53536, .16266, .29090, .00107, .00000, + .00000, .00120, .01162, + .50661, .15134, .25128, .00154, .00000, + .00000, .00172, .01096, + .48111, .14168, .21942, .00184, .00001, + .00000, .00204, .01034, + .45837, .13328, .19346, .00198, .00001, + .00001, .00221, .00988 / DATA ((F(I,34,K),I=1,5),F(7,34,K),F(6,34,K),F(8,34,K),K=1,7)/ + .62786, .19433, .42790, .00001, .00000, + .00000, .00001, .00916, + .57984, .17452, .34752, .00001, .00000, + .00000, .00001, .00903, + .53927, .15865, .28809, .00001, .00000, + .00000, .00001, .00871, + .50443, .14561, .24287, .00016, .00000, + .00000, .00065, .00832, + .47415, .13472, .20769, .00067, .00000, + .00000, .00105, .00783, + .44781, .12546, .17973, .00096, .00001, + .00000, .00131, .00743, + .42452, .11746, .15717, .00116, .00001, + .00003, .00148, .00700 / DATA ((F(I,35,K),I=1,5),F(7,35,K),F(6,35,K),F(8,35,K),K=1,7)/ + .60128, .17617, .36867, .00001, .00000, + .00000, .00001, .00600, + .54953, .15670, .29445, .00001, .00000, + .00000, .00001, .00617, + .50651, .14136, .24076, .00001, .00000, + .00000, .00001, .00596, + .47024, .12886, .20061, .00001, .00000, + .00000, .00024, .00577, + .43920, .11844, .16982, .00005, .00000, + .00000, .00059, .00551, + .41244, .10970, .14560, .00035, .00001, + .00000, .00083, .00521, + .38903, .10226, .12625, .00054, .00001, + .00004, .00095, .00490 / DATA ((F(I,36,K),I=1,5),F(7,36,K),F(6,36,K),F(8,36,K),K=1,7)/ + .57046, .15882, .31457, .00001, .00000, + .00000, .00001, .00370, + .51571, .13983, .24747, .00001, .00000, + .00000, .00001, .00397, + .47116, .12491, .19955, .00001, .00000, + .00000, .00001, .00406, + .43400, .11296, .16422, .00001, .00000, + .00000, .00001, .00391, + .40271, .10316, .13748, .00001, .00000, + .00000, .00027, .00374, + .37598, .09498, .11669, .00001, .00001, + .00000, .00048, .00356, + .35274, .08802, .10031, .00018, .00001, + .00003, .00061, .00334 / DATA ((F(I,37,K),I=1,5),F(7,37,K),F(6,37,K),F(8,37,K),K=1,7)/ + .53556, .14173, .26561, .00001, .00000, + .00000, .00001, .00219, + .47905, .12344, .20521, .00001, .00000, + .00000, .00001, .00262, + .43373, .10934, .16317, .00001, .00000, + .00000, .00001, .00266, + .39671, .09815, .13264, .00001, .00000, + .00000, .00001, .00257, + .36558, .08902, .10991, .00001, .00000, + .00000, .00001, .00250, + .33924, .08151, .09246, .00007, .00001, + .00000, .00015, .00237, + .31668, .07519, .07885, .00019, .00001, + -.00002, .00027, .00222 / DATA ((F(I,38,K),I=1,5),F(7,38,K),F(6,38,K),F(8,38,K),K=1,7)/ + .49683, .12381, .22084, .00001, .00000, + .00000, .00001, .00113, + .43944, .10672, .16749, .00001, .00000, + .00000, .00001, .00150, + .39428, .09371, .13118, .00001, .00000, + .00000, .00001, .00168, + .35761, .08355, .10540, .00001, .00000, + .00000, .00001, .00159, + .32742, .07528, .08641, .00001, .00000, + .00000, .00001, .00161, + .30203, .06857, .07204, .00001, .00001, + .00000, .00009, .00150, + .28045, .06294, .06092, .00001, .00001, + .00002, .00018, .00141 / DATA ((F(I,39,K),I=1,5),F(7,39,K),F(6,39,K),F(8,39,K),K=1,7)/ + .45431, .10635, .17814, .00001, .00000, + .00000, .00001, .00036, + .39783, .09091, .13388, .00001, .00000, + .00000, .00001, .00072, + .35381, .07923, .10374, .00001, .00000, + .00000, .00001, .00084, + .31856, .07014, .08246, .00001, .00000, + .00000, .00001, .00096, + .28960, .06287, .06697, .00001, .00000, + .00000, .00001, .00093, + .26556, .05692, .05530, .00001, .00001, + .00000, .00004, .00092, + .24521, .05198, .04637, .00001, .00001, + .00004, .00011, .00086 / DATA ((F(I,40,K),I=1,5),F(7,40,K),F(6,40,K),F(8,40,K),K=1,7)/ + .41254, .09148, .14592, .00001, .00000, + .00000, .00001, .00003, + .35692, .07724, .10695, .00001, .00000, + .00000, .00001, .00033, + .31428, .06668, .08140, .00001, .00000, + .00000, .00001, .00047, + .28057, .05850, .06381, .00001, .00000, + .00000, .00001, .00053, + .25332, .05206, .05122, .00001, .00000, + .00000, .00001, .00056, + .23091, .04683, .04190, .00001, .00001, + .00000, .00001, .00054, + .21210, .04253, .03482, .00001, .00001, + .00000, .00001, .00053 / DATA ((F(I,41,K),I=1,5),F(7,41,K),F(6,41,K),F(8,41,K),K=1,7)/ + .36734, .07563, .11355, .00001, .00000, + .00000, .00001, .00001, + .31446, .06338, .08237, .00001, .00000, + .00000, .00001, .00001, + .27437, .05433, .06203, .00001, .00000, + .00000, .00001, .00017, + .24304, .04740, .04812, .00001, .00000, + .00000, .00001, .00024, + .21785, .04189, .03825, .00001, .00000, + .00000, .00001, .00027, + .19734, .03751, .03100, .00001, .00001, + .00000, .00001, .00030, + .18021, .03390, .02554, .00001, .00001, + .00003, .00001, .00029 / DATA ((F(I,42,K),I=1,5),F(7,42,K),F(6,42,K),F(8,42,K),K=1,7)/ + .27911, .05064, .06747, .00001, .00000, + .00000, .00001, .00001, + .23354, .04145, .04718, .00001, .00000, + .00000, .00001, .00001, + .19979, .03487, .03442, .00001, .00000, + .00000, .00001, .00001, + .17401, .02990, .02599, .00001, .00000, + .00000, .00001, .00005, + .15364, .02606, .02017, .00001, .00000, + .00000, .00001, .00006, + .13734, .02302, .01600, .00001, .00001, + .00000, .00001, .00008, + .12395, .02056, .01294, .00001, .00001, + .00000, .00001, .00008 / DATA ((F(I,43,K),I=1,5),F(7,43,K),F(6,43,K),F(8,43,K),K=1,7)/ + .19610, .03075, .03606, .00001, .00000, + .00000, .00001, .00001, + .16013, .02462, .02425, .00001, .00000, + .00000, .00001, .00001, + .13421, .02028, .01712, .00001, .00000, + .00000, .00001, .00001, + .11479, .01712, .01257, .00001, .00000, + .00000, .00001, .00001, + .09976, .01469, .00952, .00001, .00000, + .00000, .00001, .00001, + .08792, .01280, .00739, .00000, .00001, + .00000, .00001, .00001, + .07838, .01129, .00586, .00001, .00001, + .00000, .00001, .00000 / DATA ((F(I,44,K),I=1,5),F(7,44,K),F(6,44,K),F(8,44,K),K=1,7)/ + .12471, .01656, .01668, .00001, .00000, + .00000, .00001, .00001, + .09913, .01295, .01080, .00001, .00000, + .00000, .00001, .00001, + .08124, .01044, .00737, .00001, .00000, + .00000, .00001, .00001, + .06814, .00865, .00525, .00001, .00000, + .00000, .00001, .00001, + .05819, .00728, .00387, .00000, .00000, + .00000, .00001, .00000, + .05052, .00627, .00294, .00001, .00001, + .00000, .00001, .00001, + .04438, .00545, .00227, .00001, .00001, + .00000, .00001, .00001 / DATA ((F(I,45,K),I=1,5),F(7,45,K),F(6,45,K),F(8,45,K),K=1,7)/ + .06789, .00717, .00626, .00001, .00000, + .00000, .00001, .00001, + .05238, .00549, .00388, .00001, .00000, + .00000, .00001, .00001, + .04184, .00436, .00255, .00001, .00000, + .00000, .00001, .00001, + .03429, .00357, .00176, .00001, .00000, + .00000, .00001, .00001, + .02872, .00295, .00126, .00001, .00000, + .00000, .00001, .00001, + .02448, .00252, .00093, .00001, .00001, + .00000, .00001, .00001, + .02116, .00218, .00070, .00001, .00001, + .00000, .00001, .00001 / DATA ((F(I,46,K),I=1,5),F(7,46,K),F(6,46,K),F(8,46,K),K=1,7)/ + .00856, .00063, .00037, .00001, .00000, + .00000, .00001, .00001, + .00610, .00043, .00017, .00001, .00000, + .00000, .00001, .00001, + .00453, .00032, .00010, .00001, .00000, + .00000, .00001, .00001, + .00349, .00023, .00006, .00001, .00000, + .00000, .00001, .00001, + .00277, .00019, .00004, .00001, .00000, + .00000, .00001, .00001, + .00226, .00014, .00003, .00001, .00001, + .00000, .00001, .00001, + .00187, .00011, .00002, .00001, .00001, + .00000, .00001, .00001 / DATA XX/1.D-5,2.D-5,4.D-5,6.D-5,8.D-5, . 1.D-4,2.D-4,4.D-4,6.D-4,8.D-4, . 1.D-3,2.D-3,4.D-3,6.D-3,8.D-3, . 1.D-2,2.D-2,4.D-2,6.D-2,8.D-2, . .1D0,.125D0,.15D0,.175D0,.2D0,.225D0,.25D0,.275D0, . .3D0,.325D0,.35D0,.375D0,.4D0,.425D0,.45D0,.475D0, . .5D0,.525D0,.55D0,.575D0,.6D0,.65D0,.7D0,.75D0, . .8D0,.9D0,1.D0/ DATA XMIN,XMAX,QSQMIN,QSQMAX/1.D-5,1.D0,0.625D0,5.D0/ DATA N0/2,5,5,9,0,0,9,9/ DATA ZEROD/0.D0/,PONED/0.1D0/,ONED/1.D0/,ONEDO/1.1D0/,TWOD/2.D0/ DATA INIT/0/ C don't let x be altered if it's out of range!! xsave=x C IF(INIT.NE.0) GOTO 10 INIT=1 DO 20 N=1,NX-1 DO 20 M=1,7 C READ(36,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), C . F(6,N,M),F(8,N,M) C 1=UV 2=DV 3=GLUE 4=UBAR 5=CBAR 7=BBAR 6=SBAR 8=DBAR DO 25 I=1,8 25 F(I,N,M)=F(I,N,M)/(ONED-XX(N))**N0(I) 20 CONTINUE DO 31 J=1,NTENTH-1 XX(J)= LOG10(XX(J))+ONEDO DO 31 I=1,8 IF(I.EQ.7.OR.I.EQ.5) GO TO 31 DO 30 K=1,7 30 F(I,J,K)= LOG(F(I,J,K))*F(I,NTENTH,K)/ LOG(F(I,NTENTH,K)) 31 CONTINUE 50 FORMAT(8F10.5) DO 40 I=1,8 DO 40 M=1,7 40 F(I,nx,M)=ZEROD 10 CONTINUE IF(X.LT.XMIN) X=XMIN IF(X.GT.XMAX) X=XMAX QSQ=SCALE**2 IF(QSQ.LT.QSQMIN) QSQ=QSQMIN IF(QSQ.GT.QSQMAX) QSQ=QSQMAX XXX=X IF(X.LT.PONED) XXX= LOG10(X)+ONEDO N=0 70 N=N+1 IF(XXX.GT.XX(N+1)) GOTO 70 A=(XXX-XX(N))/(XX(N+1)-XX(N)) RM= LOG(QSQ/QSQMIN)/ LOG(TWOD)*TWOD B=RM-AINT(RM) M=1+ INT(RM) DO 60 I=1,8 G(I)= (ONED-A)*(ONED-B)*F(I,N,M)+(ONED-A)*B*F(I,N,M+1) . + A*(ONED-B)*F(I,N+1,M) + A*B*F(I,N+1,M+1) IF(N.GE.ntenth) GOTO 65 IF(I.EQ.7.or.i.eq.5) GOTO 65 FAC=(ONED-B)*F(I,ntenth,M)+B*F(I,ntenth,M+1) G(I)=FAC**(G(I)/FAC) 65 CONTINUE G(I)=G(I)*(ONED-X)**N0(I) 60 CONTINUE UPV=G(1) DNV=G(2) USEA=G(4) DSEA=G(8) STR=G(6) CHM=G(5) GLU=G(3) BOT=G(7) C restore x x=xsave C RETURN END