* * $Id$ * * $Log$ * Revision 1.1.1.1 1996/10/30 08:32:13 cernlib * Version 7.04 * * Revision 1.2 1996/10/28 15:09:53 pdflib * PDFLIB new version 7.04 * * Revision 1.1.1.1 1996/04/12 15:30:16 plothow * Version 7.01 * * #include "pdf/pilot.h" SUBROUTINE SFMSQA(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,GLU) C C THIS IS THE NEW "A" FIT -- May 1994 -- low Q^2 range 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)/ + 0.00232, 0.00543, 4.23812, 1.92427, 0.00000, + 0.00000, 0.53025, 1.91803, + 0.00262, 0.00566, 6.42705, 1.98441, 0.00000, + 0.00000, 0.62266, 1.97884, + 0.00289, 0.00586, 9.43898, 2.07363, 0.00000, + 0.00000, 0.73934, 2.06869, + 0.00315, 0.00605, 13.16286, 2.19098, 0.00000, + 0.00000, 0.88050, 2.18663, + 0.00339, 0.00622, 17.54780, 2.32948, 0.00000, + 0.00000, 1.04003, 2.32568, + 0.00362, 0.00638, 22.55753, 2.48780, 0.05952, + 0.00000, 1.21723, 2.48453, + 0.00383, 0.00653, 28.15842, 2.66466, 0.26782, + 0.00000, 1.41125, 2.66188 / DATA ((F(I,2,K),I=1,5),F(7,2,K),F(6,2,K),F(8,2,K),K=1,7)/ + 0.00364, 0.00691, 3.22632, 1.55076, 0.00000, + 0.00000, 0.42116, 1.54806, + 0.00407, 0.00722, 5.15995, 1.60401, 0.00000, + 0.00000, 0.50030, 1.60175, + 0.00446, 0.00749, 7.71045, 1.67947, 0.00000, + 0.00000, 0.59775, 1.67760, + 0.00483, 0.00774, 10.78281, 1.77603, 0.00000, + 0.00000, 0.71335, 1.77454, + 0.00516, 0.00797, 14.34431, 1.88827, 0.00000, + 0.00000, 0.84238, 1.88712, + 0.00548, 0.00819, 18.35565, 2.01500, 0.04882, + 0.00000, 0.98416, 2.01418, + 0.00579, 0.00839, 22.78315, 2.15515, 0.21426, + 0.00000, 1.13797, 2.15464 / DATA ((F(I,3,K),I=1,5),F(7,3,K),F(6,3,K),F(8,3,K),K=1,7)/ + 0.00558, 0.00879, 2.45661, 1.24947, 0.00000, + 0.00000, 0.33447, 1.24914, + 0.00619, 0.00921, 4.14420, 1.29621, 0.00000, + 0.00000, 0.40201, 1.29616, + 0.00675, 0.00958, 6.30016, 1.35987, 0.00000, + 0.00000, 0.48332, 1.36007, + 0.00727, 0.00991, 8.83462, 1.43927, 0.00000, + 0.00000, 0.57796, 1.43970, + 0.00774, 0.01022, 11.72666, 1.53017, 0.00000, + 0.00000, 0.68229, 1.53081, + 0.00820, 0.01052, 14.93678, 1.63155, 0.04004, + 0.00000, 0.79569, 1.63239, + 0.00862, 0.01079, 18.43334, 1.74251, 0.17137, + 0.00000, 0.91755, 1.74354 / DATA ((F(I,4,K),I=1,5),F(7,4,K),F(6,4,K),F(8,4,K),K=1,7)/ + 0.00710, 0.01013, 2.09495, 1.10097, 0.00000, + 0.00000, 0.29229, 1.10164, + 0.00785, 0.01062, 3.64657, 1.14412, 0.00000, + 0.00000, 0.35376, 1.14500, + 0.00854, 0.01106, 5.59929, 1.20168, 0.00000, + 0.00000, 0.42685, 1.20276, + 0.00917, 0.01146, 7.86366, 1.27246, 0.00000, + 0.00000, 0.51103, 1.27372, + 0.00976, 0.01183, 10.42364, 1.35279, 0.00000, + 0.00000, 0.60316, 1.35421, + 0.01031, 0.01218, 13.24045, 1.44172, 0.03565, + 0.00000, 0.70264, 1.44329, + 0.01083, 0.01250, 16.28428, 1.53846, 0.15036, + 0.00000, 0.80894, 1.54017 / DATA ((F(I,5,K),I=1,5),F(7,5,K),F(6,5,K),F(8,5,K),K=1,7)/ + 0.00840, 0.01120, 1.87136, 1.00635, 0.00000, + 0.00000, 0.26564, 1.00760, + 0.00927, 0.01175, 3.33089, 1.04705, 0.00000, + 0.00000, 0.32310, 1.04848, + 0.01006, 0.01225, 5.15064, 1.10062, 0.00000, + 0.00000, 0.39086, 1.10222, + 0.01079, 0.01271, 7.24090, 1.16583, 0.00000, + 0.00000, 0.46831, 1.16759, + 0.01147, 0.01312, 9.58839, 1.23939, 0.00000, + 0.00000, 0.55265, 1.24128, + 0.01210, 0.01351, 12.15519, 1.32041, 0.03284, + 0.00000, 0.64329, 1.32243, + 0.01271, 0.01388, 14.91292, 1.40814, 0.13702, + 0.00000, 0.73975, 1.41028 / DATA ((F(I,6,K),I=1,5),F(7,6,K),F(6,6,K),F(8,6,K),K=1,7)/ + 0.00955, 0.01211, 1.71468, 0.93853, 0.00000, + 0.00000, 0.24666, 0.94017, + 0.01052, 0.01272, 3.10549, 0.97740, 0.00000, + 0.00000, 0.30117, 0.97921, + 0.01141, 0.01327, 4.82815, 1.02804, 0.00000, + 0.00000, 0.36505, 1.03000, + 0.01223, 0.01377, 6.79254, 1.08923, 0.00000, + 0.00000, 0.43767, 1.09133, + 0.01298, 0.01423, 8.98726, 1.15792, 0.00000, + 0.00000, 0.51641, 1.16015, + 0.01369, 0.01466, 11.37514, 1.23328, 0.03081, + 0.00000, 0.60073, 1.23562, + 0.01437, 0.01506, 13.92903, 1.31459, 0.12748, + 0.00000, 0.69017, 1.31704 / DATA ((F(I,7,K),I=1,5),F(7,7,K),F(6,7,K),F(8,7,K),K=1,7)/ + 0.01411, 0.01545, 1.30792, 0.75531, 0.00000, + 0.00000, 0.19600, 0.75794, + 0.01550, 0.01627, 2.50161, 0.78885, 0.00000, + 0.00000, 0.24221, 0.79165, + 0.01677, 0.01701, 3.95367, 0.83129, 0.00000, + 0.00000, 0.29537, 0.83423, + 0.01792, 0.01768, 5.57292, 0.88143, 0.00000, + 0.00000, 0.35475, 0.88450, + 0.01899, 0.01830, 7.35237, 0.93689, 0.00000, + 0.00000, 0.41833, 0.94007, + 0.01999, 0.01887, 9.25794, 0.99695, 0.02529, + 0.00000, 0.48564, 1.00024, + 0.02093, 0.01941, 11.26643, 1.06105, 0.10183, + 0.00000, 0.55629, 1.06444 / DATA ((F(I,8,K),I=1,5),F(7,8,K),F(6,8,K),F(8,8,K),K=1,7)/ + 0.02070, 0.01979, 0.99986, 0.60732, 0.00000, + 0.00000, 0.15592, 0.61082, + 0.02268, 0.02090, 2.02268, 0.63609, 0.00000, + 0.00000, 0.19503, 0.63980, + 0.02447, 0.02190, 3.24642, 0.67153, 0.00000, + 0.00000, 0.23921, 0.67542, + 0.02610, 0.02280, 4.58011, 0.71248, 0.00000, + 0.00000, 0.28769, 0.71654, + 0.02760, 0.02362, 6.02017, 0.75711, 0.00000, + 0.00000, 0.33893, 0.76132, + 0.02900, 0.02439, 7.53631, 0.80481, 0.02079, + 0.00000, 0.39253, 0.80915, + 0.03032, 0.02511, 9.10952, 0.85512, 0.08122, + 0.00000, 0.44819, 0.85959 / DATA ((F(I,9,K),I=1,5),F(7,9,K),F(6,9,K),F(8,9,K),K=1,7)/ + 0.02584, 0.02295, 0.85607, 0.53431, 0.00000, + 0.00000, 0.13654, 0.53837, + 0.02827, 0.02427, 1.79177, 0.56056, 0.00000, + 0.00000, 0.17198, 0.56487, + 0.03046, 0.02546, 2.89947, 0.59238, 0.00000, + 0.00000, 0.21160, 0.59693, + 0.03245, 0.02653, 4.08932, 0.62869, 0.00000, + 0.00000, 0.25460, 0.63343, + 0.03428, 0.02751, 5.35969, 0.66790, 0.00000, + 0.00000, 0.29971, 0.67282, + 0.03599, 0.02842, 6.68285, 0.70946, 0.01856, + 0.00000, 0.34653, 0.71456, + 0.03759, 0.02928, 8.04229, 0.75299, 0.07107, + 0.00000, 0.39483, 0.75825 / DATA ((F(I,10,K),I=1,5),F(7,10,K),F(6,10,K),F(8,10,K),K=1,7)/ + 0.03025, 0.02555, 0.76776, 0.48778, 0.00000, + 0.00000, 0.12435, 0.49229, + 0.03306, 0.02705, 1.64767, 0.51235, 0.00000, + 0.00000, 0.15741, 0.51717, + 0.03558, 0.02839, 2.68028, 0.54182, 0.00000, + 0.00000, 0.19406, 0.54690, + 0.03787, 0.02960, 3.77709, 0.57510, 0.00000, + 0.00000, 0.23353, 0.58042, + 0.03997, 0.03071, 4.93803, 0.61081, 0.00000, + 0.00000, 0.27468, 0.61634, + 0.04193, 0.03174, 6.13732, 0.64844, 0.01714, + 0.00000, 0.31716, 0.65417, + 0.04376, 0.03271, 7.36026, 0.68763, 0.06460, + 0.00000, 0.36076, 0.69355 / DATA ((F(I,11,K),I=1,5),F(7,11,K),F(6,11,K),F(8,11,K),K=1,7)/ + 0.03425, 0.02784, 0.70475, 0.45443, 0.00000, + 0.00000, 0.11569, 0.45936, + 0.03738, 0.02950, 1.54754, 0.47783, 0.00000, + 0.00000, 0.14703, 0.48308, + 0.04019, 0.03099, 2.52696, 0.50558, 0.00000, + 0.00000, 0.18150, 0.51113, + 0.04274, 0.03233, 3.55628, 0.53668, 0.00000, + 0.00000, 0.21838, 0.54247, + 0.04507, 0.03355, 4.63730, 0.56984, 0.00000, + 0.00000, 0.25663, 0.57587, + 0.04724, 0.03469, 5.74604, 0.60460, 0.01613, + 0.00000, 0.29594, 0.61084, + 0.04928, 0.03574, 6.86937, 0.64065, 0.05990, + 0.00000, 0.33610, 0.64709 / DATA ((F(I,12,K),I=1,5),F(7,12,K),F(6,12,K),F(8,12,K),K=1,7)/ + 0.05034, 0.03672, 0.54955, 0.36447, 0.00000, + 0.00000, 0.09298, 0.37107, + 0.05473, 0.03900, 1.28832, 0.38445, 0.00000, + 0.00000, 0.11948, 0.39149, + 0.05865, 0.04103, 2.11904, 0.40739, 0.00000, + 0.00000, 0.14795, 0.41483, + 0.06217, 0.04285, 2.96136, 0.43239, 0.00000, + 0.00000, 0.17770, 0.44019, + 0.06539, 0.04450, 3.82264, 0.45855, 0.00000, + 0.00000, 0.20805, 0.46667, + 0.06837, 0.04602, 4.68465, 0.48551, 0.01341, + 0.00000, 0.23875, 0.49393, + 0.07115, 0.04744, 5.53899, 0.51305, 0.04724, + 0.00000, 0.26968, 0.52174 / DATA ((F(I,13,K),I=1,5),F(7,13,K),F(6,13,K),F(8,13,K),K=1,7)/ + 0.07424, 0.04960, 0.45006, 0.29268, 0.00000, + 0.00000, 0.07588, 0.30160, + 0.08032, 0.05272, 1.10456, 0.30967, 0.00000, + 0.00000, 0.09822, 0.31921, + 0.08571, 0.05548, 1.80790, 0.32847, 0.00000, + 0.00000, 0.12157, 0.33855, + 0.09051, 0.05792, 2.49068, 0.34832, 0.00000, + 0.00000, 0.14534, 0.35889, + 0.09488, 0.06014, 3.16635, 0.36864, 0.00000, + 0.00000, 0.16911, 0.37964, + 0.09889, 0.06216, 3.82305, 0.38917, 0.01126, + 0.00000, 0.19272, 0.40057, + 0.10262, 0.06404, 4.45704, 0.40977, 0.03700, + 0.00000, 0.21612, 0.42153 / DATA ((F(I,14,K),I=1,5),F(7,14,K),F(6,14,K),F(8,14,K),K=1,7)/ + 0.09337, 0.05991, 0.42083, 0.25795, 0.00000, + 0.00000, 0.06825, 0.26877, + 0.10067, 0.06365, 1.03073, 0.27332, 0.00000, + 0.00000, 0.08842, 0.28489, + 0.10710, 0.06692, 1.66542, 0.28992, 0.00000, + 0.00000, 0.10912, 0.30215, + 0.11281, 0.06980, 2.26414, 0.30711, 0.00000, + 0.00000, 0.12983, 0.31991, + 0.11797, 0.07239, 2.84358, 0.32444, 0.00000, + 0.00000, 0.15026, 0.33776, + 0.12269, 0.07475, 3.39582, 0.34172, 0.01025, + 0.00000, 0.17033, 0.35551, + 0.12706, 0.07692, 3.91963, 0.35886, 0.03197, + -0.00001, 0.18999, 0.37308 / DATA ((F(I,15,K),I=1,5),F(7,15,K),F(6,15,K),F(8,15,K),K=1,7)/ + 0.11001, 0.06899, 0.41558, 0.23634, 0.00000, + 0.00000, 0.06388, 0.24884, + 0.11827, 0.07320, 0.99416, 0.25060, 0.00000, + 0.00000, 0.08260, 0.26394, + 0.12553, 0.07686, 1.58125, 0.26572, 0.00000, + 0.00000, 0.10153, 0.27981, + 0.13194, 0.08007, 2.12318, 0.28114, 0.00000, + 0.00000, 0.12023, 0.29587, + 0.13771, 0.08294, 2.63867, 0.29650, 0.00000, + 0.00000, 0.13849, 0.31181, + 0.14297, 0.08555, 3.12247, 0.31167, 0.00963, + 0.00000, 0.15625, 0.32750, + 0.14782, 0.08793, 3.57509, 0.32656, 0.02878, + -0.00001, 0.17350, 0.34287 / DATA ((F(I,16,K),I=1,5),F(7,16,K),F(6,16,K),F(8,16,K),K=1,7)/ + 0.12504, 0.07729, 0.42079, 0.22117, 0.00000, + 0.00000, 0.06107, 0.23523, + 0.13412, 0.08187, 0.97466, 0.23457, 0.00000, + 0.00000, 0.07870, 0.24957, + 0.14205, 0.08584, 1.52498, 0.24857, 0.00000, + 0.00000, 0.09633, 0.26439, + 0.14903, 0.08930, 2.02410, 0.26267, 0.00000, + 0.00000, 0.11355, 0.27919, + 0.15529, 0.09239, 2.49209, 0.27658, 0.00000, + 0.00000, 0.13022, 0.29372, + 0.16099, 0.09517, 2.92573, 0.29019, 0.00921, + 0.00000, 0.14631, 0.30789, + 0.16621, 0.09771, 3.32673, 0.30346, 0.02650, + 0.00000, 0.16182, 0.32167 / DATA ((F(I,17,K),I=1,5),F(7,17,K),F(6,17,K),F(8,17,K),K=1,7)/ + 0.18747, 0.11240, 0.49179, 0.18182, 0.00000, + 0.00000, 0.05541, 0.20302, + 0.19924, 0.11804, 0.96223, 0.19246, 0.00000, + 0.00000, 0.06983, 0.21485, + 0.20936, 0.12282, 1.39560, 0.20301, 0.00000, + 0.00000, 0.08362, 0.22640, + 0.21810, 0.12689, 1.76565, 0.21315, 0.00000, + 0.00000, 0.09656, 0.23738, + 0.22585, 0.13046, 2.09508, 0.22280, 0.00000, + 0.00000, 0.10870, 0.24776, + 0.23280, 0.13361, 2.38648, 0.23194, 0.00815, + 0.00000, 0.12008, 0.25755, + 0.23909, 0.13644, 2.64474, 0.24060, 0.02040, + 0.00000, 0.13077, 0.26678 / DATA ((F(I,18,K),I=1,5),F(7,18,K),F(6,18,K),F(8,18,K),K=1,7)/ + 0.28533, 0.16728, 0.65201, 0.15062, 0.00000, + 0.00000, 0.05345, 0.18428, + 0.29917, 0.17291, 1.01515, 0.15795, 0.00000, + 0.00000, 0.06450, 0.19272, + 0.31073, 0.17747, 1.31959, 0.16473, 0.00000, + 0.00000, 0.07446, 0.20037, + 0.32045, 0.18119, 1.56052, 0.17084, 0.00000, + 0.00000, 0.08333, 0.20718, + 0.32885, 0.18432, 1.76050, 0.17637, 0.00000, + 0.00000, 0.09129, 0.21328, + 0.33619, 0.18698, 1.92619, 0.18137, 0.00727, + 0.00000, 0.09847, 0.21874, + 0.34271, 0.18927, 2.06404, 0.18589, 0.01537, + -0.00001, 0.10497, 0.22364 / DATA ((F(I,19,K),I=1,5),F(7,19,K),F(6,19,K),F(8,19,K),K=1,7)/ + 0.36729, 0.21019, 0.77868, 0.13188, 0.00000, + 0.00000, 0.05237, 0.17581, + 0.38087, 0.21442, 1.06122, 0.13689, 0.00000, + 0.00000, 0.06115, 0.18142, + 0.39186, 0.21764, 1.28184, 0.14126, 0.00000, + 0.00000, 0.06872, 0.18618, + 0.40081, 0.22007, 1.44554, 0.14496, 0.00000, + 0.00000, 0.07516, 0.19012, + 0.40832, 0.22197, 1.57317, 0.14815, 0.00000, + 0.00000, 0.08074, 0.19344, + 0.41470, 0.22344, 1.67235, 0.15087, 0.00664, + 0.00000, 0.08561, 0.19624, + 0.42019, 0.22459, 1.74940, 0.15320, 0.01264, + 0.00000, 0.08988, 0.19859 / DATA ((F(I,20,K),I=1,5),F(7,20,K),F(6,20,K),F(8,20,K),K=1,7)/ + 0.43943, 0.24421, 0.87132, 0.11596, 0.00000, + 0.00000, 0.05027, 0.16762, + 0.45119, 0.24638, 1.08735, 0.11925, 0.00000, + 0.00000, 0.05727, 0.17074, + 0.46038, 0.24769, 1.24495, 0.12190, 0.00000, + 0.00000, 0.06305, 0.17318, + 0.46749, 0.24839, 1.35369, 0.12398, 0.00000, + 0.00000, 0.06779, 0.17492, + 0.47318, 0.24870, 1.43230, 0.12564, 0.00068, + 0.00000, 0.07176, 0.17621, + 0.47780, 0.24870, 1.48820, 0.12692, 0.00602, + 0.00000, 0.07509, 0.17712, + 0.48154, 0.24849, 1.52710, 0.12792, 0.01066, + -0.00001, 0.07792, 0.17772 / DATA ((F(I,21,K),I=1,5),F(7,21,K),F(6,21,K),F(8,21,K),K=1,7)/ + 0.50349, 0.27060, 0.93597, 0.10110, 0.00000, + 0.00000, 0.04723, 0.15793, + 0.51234, 0.27032, 1.09539, 0.10314, 0.00000, + 0.00000, 0.05277, 0.15902, + 0.51876, 0.26952, 1.20287, 0.10461, 0.00000, + 0.00000, 0.05718, 0.15957, + 0.52327, 0.26834, 1.26983, 0.10560, 0.00000, + 0.00000, 0.06063, 0.15964, + 0.52654, 0.26694, 1.31262, 0.10625, 0.00107, + 0.00000, 0.06341, 0.15943, + 0.52884, 0.26542, 1.33802, 0.10663, 0.00537, + 0.00000, 0.06566, 0.15898, + 0.53041, 0.26380, 1.35090, 0.10680, 0.00904, + 0.00000, 0.06749, 0.15835 / DATA ((F(I,22,K),I=1,5),F(7,22,K),F(6,22,K),F(8,22,K),K=1,7)/ + 0.57288, 0.29429, 0.98361, 0.08377, 0.00000, + 0.00000, 0.04242, 0.14370, + 0.57693, 0.29082, 1.08372, 0.08476, 0.00000, + 0.00000, 0.04654, 0.14280, + 0.57901, 0.28728, 1.14164, 0.08529, 0.00000, + 0.00000, 0.04963, 0.14161, + 0.57959, 0.28371, 1.16898, 0.08543, 0.00000, + 0.00000, 0.05191, 0.14021, + 0.57926, 0.28021, 1.17918, 0.08531, 0.00123, + 0.00000, 0.05364, 0.13871, + 0.57826, 0.27680, 1.17779, 0.08501, 0.00457, + 0.00000, 0.05495, 0.13713, + 0.57673, 0.27349, 1.16857, 0.08458, 0.00737, + -0.00001, 0.05594, 0.13551 / DATA ((F(I,23,K),I=1,5),F(7,23,K),F(6,23,K),F(8,23,K),K=1,7)/ + 0.63048, 0.30915, 1.00111, 0.06818, 0.00000, + 0.00000, 0.03692, 0.12784, + 0.62891, 0.30248, 1.05264, 0.06852, 0.00000, + 0.00000, 0.03999, 0.12557, + 0.62606, 0.29626, 1.07218, 0.06847, 0.00000, + 0.00000, 0.04216, 0.12327, + 0.62227, 0.29042, 1.07015, 0.06812, 0.00000, + 0.00000, 0.04363, 0.12094, + 0.61801, 0.28497, 1.05720, 0.06761, 0.00119, + 0.00000, 0.04466, 0.11866, + 0.61343, 0.27986, 1.03756, 0.06698, 0.00380, + 0.00000, 0.04536, 0.11642, + 0.60864, 0.27503, 1.01392, 0.06628, 0.00596, + -0.00003, 0.04583, 0.11422 / DATA ((F(I,24,K),I=1,5),F(7,24,K),F(6,24,K),F(8,24,K),K=1,7)/ + 0.67676, 0.31587, 0.99514, 0.05418, 0.00000, + 0.00000, 0.03156, 0.11143, + 0.66899, 0.30635, 1.00701, 0.05425, 0.00000, + 0.00000, 0.03376, 0.10833, + 0.66082, 0.29778, 0.99689, 0.05396, 0.00000, + 0.00000, 0.03522, 0.10534, + 0.65240, 0.28999, 0.97320, 0.05343, 0.00000, + 0.00000, 0.03611, 0.10251, + 0.64407, 0.28287, 0.94399, 0.05278, 0.00109, + 0.00000, 0.03666, 0.09982, + 0.63589, 0.27630, 0.91213, 0.05205, 0.00312, + 0.00000, 0.03697, 0.09728, + 0.62786, 0.27019, 0.87934, 0.05126, 0.00479, + -0.00001, 0.03709, 0.09488 / DATA ((F(I,25,K),I=1,5),F(7,25,K),F(6,25,K),F(8,25,K),K=1,7)/ + 0.71132, 0.31625, 0.97091, 0.04237, 0.00000, + 0.00000, 0.02619, 0.09523, + 0.69723, 0.30412, 0.95082, 0.04228, 0.00000, + 0.00000, 0.02783, 0.09170, + 0.68374, 0.29348, 0.91809, 0.04190, 0.00000, + 0.00000, 0.02882, 0.08844, + 0.67077, 0.28400, 0.87878, 0.04130, 0.00000, + 0.00000, 0.02936, 0.08542, + 0.65848, 0.27548, 0.83847, 0.04063, 0.00090, + 0.00000, 0.02962, 0.08261, + 0.64682, 0.26772, 0.79881, 0.03991, 0.00251, + 0.00000, 0.02968, 0.08002, + 0.63571, 0.26057, 0.76061, 0.03916, 0.00380, + -0.00002, 0.02961, 0.07759 / DATA ((F(I,26,K),I=1,5),F(7,26,K),F(6,26,K),F(8,26,K),K=1,7)/ + 0.73462, 0.31108, 0.93107, 0.03255, 0.00000, + 0.00000, 0.02129, 0.07983, + 0.71430, 0.29676, 0.88680, 0.03241, 0.00000, + 0.00000, 0.02250, 0.07623, + 0.69565, 0.28441, 0.83766, 0.03201, 0.00000, + 0.00000, 0.02319, 0.07296, + 0.67837, 0.27354, 0.78764, 0.03144, 0.00000, + 0.00000, 0.02348, 0.07003, + 0.66238, 0.26391, 0.74029, 0.03083, 0.00072, + 0.00000, 0.02358, 0.06730, + 0.64749, 0.25524, 0.69616, 0.03017, 0.00199, + 0.00000, 0.02352, 0.06481, + 0.63353, 0.24733, 0.65531, 0.02949, 0.00299, + -0.00002, 0.02335, 0.06253 / DATA ((F(I,27,K),I=1,5),F(7,27,K),F(6,27,K),F(8,27,K),K=1,7)/ + 0.74766, 0.30066, 0.87914, 0.02401, 0.00000, + 0.00000, 0.01723, 0.06581, + 0.72123, 0.28478, 0.81691, 0.02404, 0.00000, + 0.00000, 0.01807, 0.06235, + 0.69764, 0.27126, 0.75668, 0.02381, 0.00000, + 0.00000, 0.01849, 0.05925, + 0.67627, 0.25949, 0.70017, 0.02341, 0.00000, + 0.00000, 0.01860, 0.05648, + 0.65688, 0.24911, 0.64913, 0.02292, 0.00062, + 0.00000, 0.01855, 0.05399, + 0.63909, 0.23985, 0.60318, 0.02241, 0.00160, + 0.00000, 0.01840, 0.05173, + 0.62259, 0.23148, 0.56179, 0.02188, 0.00236, + 0.00002, 0.01819, 0.04966 / DATA ((F(I,28,K),I=1,5),F(7,28,K),F(6,28,K),F(8,28,K),K=1,7)/ + 0.74936, 0.28783, 0.82378, 0.01828, 0.00000, + 0.00000, 0.01309, 0.05326, + 0.71761, 0.27055, 0.74720, 0.01821, 0.00000, + 0.00000, 0.01381, 0.05010, + 0.68980, 0.25599, 0.67896, 0.01795, 0.00000, + 0.00000, 0.01415, 0.04731, + 0.66498, 0.24349, 0.61850, 0.01757, 0.00000, + 0.00000, 0.01424, 0.04482, + 0.64273, 0.23257, 0.56585, 0.01714, 0.00028, + 0.00000, 0.01418, 0.04263, + 0.62253, 0.22290, 0.51975, 0.01669, 0.00107, + 0.00000, 0.01404, 0.04065, + 0.60398, 0.21422, 0.47915, 0.01623, 0.00169, + -0.00014, 0.01382, 0.03885 / DATA ((F(I,29,K),I=1,5),F(7,29,K),F(6,29,K),F(8,29,K),K=1,7)/ + 0.74266, 0.27081, 0.75762, 0.01278, 0.00000, + 0.00000, 0.01026, 0.04237, + 0.70594, 0.25285, 0.67378, 0.01289, 0.00000, + 0.00000, 0.01077, 0.03958, + 0.67427, 0.23786, 0.60220, 0.01279, 0.00000, + 0.00000, 0.01097, 0.03715, + 0.64640, 0.22509, 0.54096, 0.01257, 0.00000, + 0.00000, 0.01097, 0.03501, + 0.62168, 0.21401, 0.48891, 0.01227, 0.00029, + 0.00000, 0.01087, 0.03314, + 0.59943, 0.20426, 0.44429, 0.01198, 0.00089, + 0.00000, 0.01072, 0.03143, + 0.57919, 0.19555, 0.40568, 0.01166, 0.00134, + -0.00002, 0.01052, 0.02990 / DATA ((F(I,30,K),I=1,5),F(7,30,K),F(6,30,K),F(8,30,K),K=1,7)/ + 0.72685, 0.25224, 0.69187, 0.00899, 0.00000, + 0.00000, 0.00768, 0.03313, + 0.68604, 0.23389, 0.60286, 0.00915, 0.00000, + 0.00000, 0.00808, 0.03076, + 0.65127, 0.21871, 0.52998, 0.00910, 0.00000, + 0.00000, 0.00824, 0.02870, + 0.62098, 0.20591, 0.46960, 0.00897, 0.00000, + 0.00000, 0.00824, 0.02690, + 0.59438, 0.19486, 0.41941, 0.00877, 0.00019, + 0.00000, 0.00814, 0.02534, + 0.57066, 0.18518, 0.37721, 0.00852, 0.00065, + 0.00000, 0.00800, 0.02398, + 0.54919, 0.17662, 0.34122, 0.00829, 0.00100, + -0.00004, 0.00784, 0.02270 / DATA ((F(I,31,K),I=1,5),F(7,31,K),F(6,31,K),F(8,31,K),K=1,7)/ + 0.70337, 0.23219, 0.62540, 0.00613, 0.00000, + 0.00000, 0.00563, 0.02548, + 0.65924, 0.21385, 0.53420, 0.00632, 0.00000, + 0.00000, 0.00597, 0.02348, + 0.62207, 0.19880, 0.46215, 0.00631, 0.00000, + 0.00000, 0.00609, 0.02181, + 0.59005, 0.18618, 0.40404, 0.00623, 0.00000, + 0.00000, 0.00610, 0.02033, + 0.56211, 0.17538, 0.35676, 0.00609, 0.00015, + 0.00000, 0.00602, 0.01909, + 0.53737, 0.16598, 0.31763, 0.00593, 0.00050, + 0.00000, 0.00592, 0.01796, + 0.51514, 0.15768, 0.28475, 0.00577, 0.00076, + 0.00000, 0.00578, 0.01695 / DATA ((F(I,32,K),I=1,5),F(7,32,K),F(6,32,K),F(8,32,K),K=1,7)/ + 0.67325, 0.21095, 0.55633, 0.00375, 0.00000, + 0.00000, 0.00420, 0.01922, + 0.62664, 0.19307, 0.46726, 0.00403, 0.00000, + 0.00000, 0.00444, 0.01763, + 0.58778, 0.17849, 0.39852, 0.00414, 0.00000, + 0.00000, 0.00451, 0.01629, + 0.55458, 0.16637, 0.34417, 0.00414, 0.00000, + 0.00000, 0.00449, 0.01513, + 0.52586, 0.15603, 0.30062, 0.00408, 0.00014, + 0.00000, 0.00443, 0.01412, + 0.50060, 0.14709, 0.26506, 0.00399, 0.00039, + 0.00000, 0.00433, 0.01324, + 0.47803, 0.13922, 0.23556, 0.00390, 0.00060, + 0.00004, 0.00422, 0.01243 / DATA ((F(I,33,K),I=1,5),F(7,33,K),F(6,33,K),F(8,33,K),K=1,7)/ + 0.63725, 0.19005, 0.49256, 0.00254, 0.00000, + 0.00000, 0.00286, 0.01421, + 0.58912, 0.17279, 0.40616, 0.00275, 0.00000, + 0.00000, 0.00308, 0.01299, + 0.54935, 0.15884, 0.34118, 0.00283, 0.00000, + 0.00000, 0.00317, 0.01190, + 0.51567, 0.14726, 0.29090, 0.00283, 0.00000, + 0.00000, 0.00317, 0.01104, + 0.48672, 0.13748, 0.25128, 0.00280, 0.00007, + 0.00000, 0.00313, 0.01025, + 0.46139, 0.12906, 0.21942, 0.00275, 0.00026, + 0.00000, 0.00307, 0.00957, + 0.43891, 0.12168, 0.19331, 0.00267, 0.00040, + 0.00001, 0.00298, 0.00902 / DATA ((F(I,34,K),I=1,5),F(7,34,K),F(6,34,K),F(8,34,K),K=1,7)/ + 0.59662, 0.16877, 0.42790, 0.00132, 0.00000, + 0.00000, 0.00206, 0.01033, + 0.54786, 0.15249, 0.34752, 0.00159, 0.00000, + 0.00000, 0.00222, 0.00936, + 0.50791, 0.13943, 0.28809, 0.00172, 0.00000, + 0.00000, 0.00228, 0.00856, + 0.47435, 0.12867, 0.24287, 0.00177, 0.00000, + 0.00000, 0.00226, 0.00790, + 0.44564, 0.11963, 0.20769, 0.00179, 0.00006, + 0.00000, 0.00221, 0.00731, + 0.42073, 0.11187, 0.17973, 0.00177, 0.00019, + 0.00000, 0.00216, 0.00682, + 0.39872, 0.10510, 0.15704, 0.00175, 0.00030, + 0.00003, 0.00210, 0.00636 / DATA ((F(I,35,K),I=1,5),F(7,35,K),F(6,35,K),F(8,35,K),K=1,7)/ + 0.55241, 0.14825, 0.36867, 0.00060, 0.00000, + 0.00000, 0.00141, 0.00731, + 0.50389, 0.13309, 0.29445, 0.00085, 0.00000, + 0.00000, 0.00154, 0.00662, + 0.46440, 0.12102, 0.24076, 0.00100, 0.00000, + 0.00000, 0.00158, 0.00601, + 0.43149, 0.11114, 0.20061, 0.00108, 0.00000, + 0.00000, 0.00156, 0.00553, + 0.40355, 0.10285, 0.16982, 0.00109, 0.00005, + 0.00000, 0.00153, 0.00513, + 0.37940, 0.09579, 0.14560, 0.00109, 0.00015, + 0.00000, 0.00150, 0.00476, + 0.35819, 0.08968, 0.12614, 0.00109, 0.00022, + 0.00004, 0.00145, 0.00442 / DATA ((F(I,36,K),I=1,5),F(7,36,K),F(6,36,K),F(8,36,K),K=1,7)/ + 0.50593, 0.12897, 0.31457, 0.00039, 0.00000, + 0.00000, 0.00086, 0.00508, + 0.45836, 0.11502, 0.24747, 0.00057, 0.00000, + 0.00000, 0.00097, 0.00457, + 0.41999, 0.10394, 0.19955, 0.00063, 0.00000, + 0.00000, 0.00101, 0.00418, + 0.38820, 0.09495, 0.16422, 0.00069, 0.00000, + 0.00000, 0.00103, 0.00380, + 0.36140, 0.08747, 0.13748, 0.00069, 0.00003, + 0.00000, 0.00103, 0.00351, + 0.33834, 0.08112, 0.11669, 0.00067, 0.00011, + 0.00000, 0.00101, 0.00325, + 0.31817, 0.07563, 0.10021, 0.00068, 0.00016, + 0.00003, 0.00099, 0.00300 / DATA ((F(I,37,K),I=1,5),F(7,37,K),F(6,37,K),F(8,37,K),K=1,7)/ + 0.45796, 0.11074, 0.26561, 0.00026, 0.00000, + 0.00000, 0.00049, 0.00348, + 0.41212, 0.09810, 0.20521, 0.00035, 0.00000, + 0.00000, 0.00057, 0.00314, + 0.37540, 0.08816, 0.16317, 0.00045, 0.00000, + 0.00000, 0.00063, 0.00283, + 0.34528, 0.08014, 0.13264, 0.00046, 0.00000, + 0.00000, 0.00065, 0.00255, + 0.31994, 0.07349, 0.10991, 0.00047, 0.00000, + 0.00000, 0.00065, 0.00236, + 0.29823, 0.06789, 0.09246, 0.00049, 0.00004, + 0.00000, 0.00063, 0.00217, + 0.27940, 0.06308, 0.07877, 0.00049, 0.00008, + -0.00002, 0.00062, 0.00199 / DATA ((F(I,38,K),I=1,5),F(7,38,K),F(6,38,K),F(8,38,K),K=1,7)/ + 0.40933, 0.09346, 0.22084, -0.00006, 0.00000, + 0.00000, 0.00036, 0.00231, + 0.36580, 0.08224, 0.16749, 0.00007, 0.00000, + 0.00000, 0.00043, 0.00205, + 0.33123, 0.07348, 0.13118, 0.00012, 0.00000, + 0.00000, 0.00046, 0.00186, + 0.30297, 0.06647, 0.10540, 0.00020, 0.00000, + 0.00000, 0.00046, 0.00166, + 0.27945, 0.06068, 0.08641, 0.00020, 0.00002, + 0.00000, 0.00045, 0.00154, + 0.25938, 0.05583, 0.07204, 0.00024, 0.00006, + 0.00000, 0.00044, 0.00140, + 0.24204, 0.05168, 0.06086, 0.00025, 0.00009, + 0.00002, 0.00043, 0.00128 / DATA ((F(I,39,K),I=1,5),F(7,39,K),F(6,39,K),F(8,39,K),K=1,7)/ + 0.36094, 0.07776, 0.17814, -0.00018, 0.00000, + 0.00000, 0.00027, 0.00147, + 0.32041, 0.06800, 0.13388, -0.00008, 0.00000, + 0.00000, 0.00030, 0.00129, + 0.28840, 0.06043, 0.10374, 0.00001, 0.00000, + 0.00000, 0.00031, 0.00114, + 0.26245, 0.05438, 0.08246, 0.00004, 0.00001, + 0.00000, 0.00031, 0.00106, + 0.24089, 0.04944, 0.06697, 0.00008, 0.00004, + 0.00000, 0.00030, 0.00095, + 0.22264, 0.04529, 0.05530, 0.00010, 0.00005, + 0.00000, 0.00029, 0.00088, + 0.20691, 0.04176, 0.04632, 0.00013, 0.00006, + 0.00004, 0.00028, 0.00080 / DATA ((F(I,40,K),I=1,5),F(7,40,K),F(6,40,K),F(8,40,K),K=1,7)/ + 0.31447, 0.06406, 0.14592, -0.00003, 0.00000, + 0.00000, 0.00008, 0.00093, + 0.27711, 0.05560, 0.10695, 0.00001, 0.00000, + 0.00000, 0.00012, 0.00081, + 0.24784, 0.04911, 0.08140, 0.00005, 0.00000, + 0.00000, 0.00014, 0.00072, + 0.22428, 0.04394, 0.06381, 0.00007, 0.00000, + 0.00000, 0.00016, 0.00065, + 0.20488, 0.03975, 0.05122, 0.00008, 0.00000, + 0.00000, 0.00016, 0.00059, + 0.18855, 0.03625, 0.04190, 0.00009, 0.00001, + 0.00000, 0.00016, 0.00054, + 0.17455, 0.03329, 0.03478, 0.00009, 0.00002, + 0.00000, 0.00015, 0.00050 / DATA ((F(I,41,K),I=1,5),F(7,41,K),F(6,41,K),F(8,41,K),K=1,7)/ + 0.26913, 0.05145, 0.11355, -0.00016, 0.00000, + 0.00000, 0.00013, 0.00052, + 0.23552, 0.04437, 0.08237, -0.00011, 0.00000, + 0.00000, 0.00013, 0.00046, + 0.20936, 0.03897, 0.06203, -0.00006, 0.00000, + 0.00000, 0.00012, 0.00042, + 0.18845, 0.03470, 0.04812, -0.00003, 0.00001, + 0.00000, 0.00013, 0.00038, + 0.17129, 0.03123, 0.03825, 0.00000, 0.00002, + 0.00000, 0.00012, 0.00034, + 0.15693, 0.02837, 0.03100, 0.00000, 0.00003, + 0.00000, 0.00012, 0.00032, + 0.14466, 0.02595, 0.02551, 0.00002, 0.00004, + 0.00003, 0.00011, 0.00029 / DATA ((F(I,42,K),I=1,5),F(7,42,K),F(6,42,K),F(8,42,K),K=1,7)/ + 0.18815, 0.03174, 0.06747, 0.00003, 0.00000, + 0.00000, 0.00001, 0.00017, + 0.16218, 0.02694, 0.04718, 0.00000, 0.00000, + 0.00000, 0.00002, 0.00017, + 0.14226, 0.02335, 0.03442, 0.00000, 0.00000, + 0.00000, 0.00002, 0.00014, + 0.12656, 0.02054, 0.02599, 0.00000, 0.00000, + 0.00000, 0.00003, 0.00013, + 0.11382, 0.01830, 0.02017, 0.00002, 0.00000, + 0.00000, 0.00003, 0.00011, + 0.10327, 0.01646, 0.01600, 0.00002, 0.00000, + 0.00000, 0.00003, 0.00010, + 0.09437, 0.01492, 0.01292, 0.00001, 0.00001, + 0.00000, 0.00003, 0.00009 / DATA ((F(I,43,K),I=1,5),F(7,43,K),F(6,43,K),F(8,43,K),K=1,7)/ + 0.12097, 0.01770, 0.03606, 0.00003, 0.00000, + 0.00000, -0.00001, 0.00005, + 0.10255, 0.01479, 0.02425, 0.00001, 0.00000, + 0.00000, 0.00000, 0.00005, + 0.08869, 0.01263, 0.01712, 0.00001, 0.00000, + 0.00000, 0.00000, 0.00004, + 0.07790, 0.01098, 0.01257, 0.00001, 0.00000, + 0.00000, 0.00000, 0.00003, + 0.06927, 0.00967, 0.00952, 0.00001, 0.00000, + 0.00000, 0.00000, 0.00003, + 0.06220, 0.00861, 0.00739, 0.00002, 0.00000, + 0.00000, 0.00000, 0.00002, + 0.05631, 0.00773, 0.00585, 0.00001, 0.00000, + 0.00000, 0.00000, 0.00002 / DATA ((F(I,44,K),I=1,5),F(7,44,K),F(6,44,K),F(8,44,K),K=1,7)/ + 0.06972, 0.00867, 0.01668, 0.00001, 0.00000, + 0.00000, -0.00001, 0.00002, + 0.05802, 0.00712, 0.01080, 0.00001, 0.00000, + 0.00000, 0.00000, 0.00002, + 0.04939, 0.00599, 0.00737, 0.00001, 0.00000, + 0.00000, 0.00000, 0.00001, + 0.04278, 0.00513, 0.00525, 0.00001, 0.00000, + 0.00000, -0.00001, 0.00001, + 0.03757, 0.00446, 0.00387, 0.00001, 0.00000, + 0.00000, 0.00000, 0.00001, + 0.03336, 0.00393, 0.00294, 0.00000, 0.00000, + 0.00000, 0.00000, 0.00001, + 0.02987, 0.00349, 0.00227, 0.00000, 0.00000, + 0.00000, 0.00000, 0.00001 / DATA ((F(I,45,K),I=1,5),F(7,45,K),F(6,45,K),F(8,45,K),K=1,7)/ + 0.03421, 0.00350, 0.00626, -0.00001, 0.00001, + 0.00000, 0.00001, 0.00000, + 0.02789, 0.00281, 0.00388, 0.00000, 0.00000, + 0.00000, 0.00000, 0.00000, + 0.02332, 0.00233, 0.00255, 0.00000, 0.00000, + 0.00000, 0.00000, 0.00000, + 0.01987, 0.00196, 0.00176, 0.00000, 0.00000, + 0.00000, 0.00000, 0.00000, + 0.01721, 0.00169, 0.00126, 0.00000, 0.00000, + 0.00000, 0.00000, 0.00000, + 0.01508, 0.00146, 0.00093, 0.00000, 0.00000, + 0.00000, 0.00000, 0.00000, + 0.01334, 0.00129, 0.00070, 0.00000, 0.00000, + 0.00000, 0.00000, 0.00000 / DATA ((F(I,46,K),I=1,5),F(7,46,K),F(6,46,K),F(8,46,K),K=1,7)/ + 0.00335, 0.00020, 0.00037, 0.00000, 0.00000, + 0.00000, 0.00000, 0.00000, + 0.00257, 0.00015, 0.00017, 0.00000, 0.00000, + 0.00000, 0.00000, 0.00000, + 0.00203, 0.00012, 0.00010, 0.00000, 0.00000, + 0.00000, 0.00000, 0.00000, + 0.00165, 0.00009, 0.00006, 0.00000, 0.00000, + 0.00000, 0.00000, 0.00000, + 0.00137, 0.00008, 0.00004, 0.00000, 0.00000, + 0.00000, 0.00000, 0.00000, + 0.00116, 0.00006, 0.00003, 0.00000, 0.00000, + 0.00000, 0.00000, 0.00000, + 0.00099, 0.00005, 0.00002, 0.00000, 0.00000, + 0.00000, 0.00000, 0.00000 / 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 xsave=x C IF(INIT.NE.0) GOTO 10 INIT=1 DO 20 N=1,nx-1 DO 20 M=1,7 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), 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=(1.D0-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 x=xsave C RETURN END