| 21886bb6 |
1 | * |
| 2 | * $Id$ |
| 3 | * |
| 4 | * $Log$ |
| 5 | * Revision 1.1.1.2 1996/10/30 08:27:43 cernlib |
| 6 | * Version 7.04 |
| 7 | * |
| 8 | * Revision 1.1.1.1 1996/04/12 15:29:09 plothow |
| 9 | * Version 7.01 |
| 10 | * |
| 11 | * |
| 12 | #include "pdf/pilot.h" |
| 13 | FUNCTION Ctq2OPd (Iset, Iprtn, XX, QQ, Irt) |
| 14 | |
| 15 | C Returns xf(x,Q) -- the momentum fraction distribution !! |
| 16 | C Returns valence and sea rather than combined flavor distr. |
| 17 | |
| 18 | C Iset : PDF set label |
| 19 | |
| 20 | C Iprtn : Parton label: 2, 1 = d_ and u_ valence |
| 21 | C 0 = gluon |
| 22 | C -1, ... -6 = u, d, s, c, b, t sea quarks |
| 23 | |
| 24 | C XX : Bjorken-x |
| 25 | C QQ : scale parameter "Q" |
| 26 | C Irt : Return code |
| 27 | C 0 : no error |
| 28 | C 1 : parametrization is slightly negative; reset to 0.0. |
| 29 | C (This condition happens rarely -- only for large x where the |
| 30 | C absolute value of the parton distribution is extremely small.) |
| 31 | |
| 32 | |
| 33 | C IMPLICIT DOUBLE PRECISION (A-H, O-Z) |
| 34 | C+SEQ, IMPDP. |
| 35 | |
| 36 | PARAMETER (D0=0D0, D1=1D0, D2=2D0, D3=3D0, D4=4D0, D10=1D1) |
| 37 | PARAMETER (Nex = 5, MxFl = 6, Npn = 3, Nst = 30, Nexpt=20) |
| 38 | Parameter (Nst4 = Nst*4) |
| 39 | |
| 40 | #include "pdf/expdp.inc" |
| 41 | + P012, P34, P5, PDFD |
| 42 | REAL PDF |
| 43 | |
| 44 | DIMENSION |
| 45 | > Iord(Nst), Isch(Nst), Nqrk(Nst),Alm(Nst) |
| 46 | > , Vlm(4:6,Nst), Qms(4:6, Nst) |
| 47 | > , Xmn(Nst), Qmn(Nst), Qmx(Nst), Nexp(Nexpt) |
| 48 | > , Mex(Nst), Mpn(Nst), ExpN(Nexpt, Nst), ExpNor(Nexpt) |
| 49 | |
| 50 | C Run le26 - CTEQ2M |
| 51 | C |
| 52 | DATA |
| 53 | > Isch(1), Iord(1), Nqrk(1), Alm(1) / 1, 2, 6, .213 / |
| 54 | > ,(Vlm(I,1), I=4,6) / .213, .139, .053 / |
| 55 | > ,(Qms(I,1), I=4,6) / 1.60, 5.00, 180.0 / |
| 56 | > ,Xmn(1), Qmn(1), Qmx(1) / 1.E-5, 1.60, 1.E3 / |
| 57 | > ,Mex(1), Mpn(1), Nexp(1) / 5, 3, 10 / |
| 58 | > ,(ExpN(I,1), I=1,10 ) |
| 59 | > / 0.990,1.012,1.022,0.980,1.062,0.870,0.843 |
| 60 | > ,0.815,0.974,1.029 / |
| 61 | |
| 62 | C Run sa17 - CTEQ2MS |
| 63 | C |
| 64 | DATA |
| 65 | > Isch(2), Iord(2), Nqrk(2), Alm(2) / 1, 2, 6, .208 / |
| 66 | > ,(Vlm(I,2), I=4,6) / .208, .135, .051 / |
| 67 | > ,(Qms(I,2), I=4,6) / 1.60, 5.00, 180.0 / |
| 68 | > ,Xmn(2), Qmn(2), Qmx(2) / 1.E-5, 1.60, 1.E3 / |
| 69 | > ,Mex(2), Mpn(2), Nexp(2) / 5, 3, 10 / |
| 70 | > ,(ExpN(I,2), I=1,10 ) |
| 71 | > / 0.992,1.017,1.023,0.982,1.079,0.879,0.845 |
| 72 | > ,0.814,0.984,1.036 / |
| 73 | |
| 74 | C Run fa06 - CTEQ2MF |
| 75 | C |
| 76 | DATA |
| 77 | > Isch(3), Iord(3), Nqrk(3), Alm(3) / 1, 2, 6, .208 / |
| 78 | > ,(Vlm(I,3), I=4,6) / .208, .135, .051 / |
| 79 | > ,(Qms(I,3), I=4,6) / 1.60, 5.00, 180.0 / |
| 80 | > ,Xmn(3), Qmn(3), Qmx(3) / 1.E-5, 1.60, 1.E3 / |
| 81 | > ,Mex(3), Mpn(3), Nexp(3) / 5, 3, 10 / |
| 82 | > ,(ExpN(I,3), I=1,10 ) |
| 83 | > / 0.989,1.014,1.021,0.977,1.099,0.812,0.789 |
| 84 | > ,0.822,0.909,0.950 / |
| 85 | |
| 86 | C Run ll25 - CTEQ2ML |
| 87 | C |
| 88 | DATA |
| 89 | > Isch(4), Iord(4), Nqrk(4), Alm(4) / 1, 2, 6, .322 / |
| 90 | > ,(Vlm(I,4), I=4,6) / .322, .220, .088 / |
| 91 | > ,(Qms(I,4), I=4,6) / 1.60, 5.00, 180.0 / |
| 92 | > ,Xmn(4), Qmn(4), Qmx(4) / 1.E-5, 1.60, 1.E3 / |
| 93 | > ,Mex(4), Mpn(4), Nexp(4) / 5, 3, 10 / |
| 94 | > ,(ExpN(I,4), I=1,10 ) |
| 95 | > / 0.990,1.009,1.022,0.978,1.050,1.027,0.955 |
| 96 | > ,0.848,0.985,1.053 / |
| 97 | |
| 98 | C Run lo24 - CTEQ2L |
| 99 | C |
| 100 | DATA |
| 101 | > Isch(5), Iord(5), Nqrk(5), Alm(5) / 0, 1, 6, .190 / |
| 102 | > ,(Vlm(I,5), I=4,6) / .190, .143, .072 / |
| 103 | > ,(Qms(I,5), I=4,6) / 1.60, 5.00, 180.0 / |
| 104 | > ,Xmn(5), Qmn(5), Qmx(5) / 1.E-5, 1.60, 1.E3 / |
| 105 | > ,Mex(5), Mpn(5), Nexp(5) / 5, 3, 10 / |
| 106 | > ,(ExpN(I,5), I=1,10 ) |
| 107 | > / 0.990,1.018,1.022,0.982,0.725,0.931,0.855 |
| 108 | > ,0.795,0.965,1.007 / |
| 109 | |
| 110 | C Run da06 - CTEQ2D |
| 111 | C |
| 112 | DATA |
| 113 | > Isch(6), Iord(6), Nqrk(6), Alm(6) / 2, 2, 6, .235 / |
| 114 | > ,(Vlm(I,6), I=4,6) / .235, .155, .060 / |
| 115 | > ,(Qms(I,6), I=4,6) / 1.60, 5.00, 180.0 / |
| 116 | > ,Xmn(6), Qmn(6), Qmx(6) / 1.E-5, 1.60, 1.E3 / |
| 117 | > ,Mex(6), Mpn(6), Nexp(6) / 5, 3, 10 / |
| 118 | > ,(ExpN(I,6), I=1,10 ) |
| 119 | > / 0.989,1.019,1.021,0.978,0.961,0.958,0.904 |
| 120 | > ,0.857,0.965,1.022 / |
| 121 | |
| 122 | Data Ist, Lp, Qsto, Aln2 / 0, -10, 1.2345, 0.6931 / |
| 123 | |
| 124 | SAVE Ist, Lp, Qsto, Alam, SBL, SB, SB2, SB3 |
| 125 | |
| 126 | X = XX |
| 127 | Irt = 0 |
| 128 | if(Iset.eq.Ist .and. Qsto.eq.QQ) then |
| 129 | C if only change is in x: |
| 130 | if (Iprtn.eq.Lp) goto 100 |
| 131 | C if change in flv is within "light" partons: |
| 132 | if (Iprtn.ge.-3 .and. Lp.ge.-3) goto 501 |
| 133 | endif |
| 134 | |
| 135 | Ip = abs(Iprtn) |
| 136 | C Set up Qi for SB |
| 137 | If (Ip .GE. 4) then |
| 138 | If (QQ .LE. Qms(Ip, Iset)) Then |
| 139 | Ctq2OPd = 0.0 |
| 140 | Return |
| 141 | Endif |
| 142 | Qi = Qms(ip, Iset) |
| 143 | Else |
| 144 | Qi = Qmn(Iset) |
| 145 | Endif |
| 146 | C Use "standard lambda" of parametrization program |
| 147 | Alam = Alm (Iset) |
| 148 | |
| 149 | SBL = LOG(QQ/Alam) / LOG(Qi/Alam) |
| 150 | SB = LOG (SBL) |
| 151 | SB2 = SB*SB |
| 152 | SB3 = SB2*SB |
| 153 | |
| 154 | 501 Iflv = 3 - Iprtn |
| 155 | |
| 156 | Goto (1,2,3,4,5,6, 311) Iset |
| 157 | |
| 158 | 1 Goto(11,12,13,14,15,16,17,18,19)Iflv |
| 159 | C Ifl = 2 |
| 160 | 11 A0=Exp( 0.2140E+00+0.1896E+01*SB -0.4966E+01*SB2+0.2354E+01*SB3) |
| 161 | A1= 0.5210E+00-0.4886E+00*SB +0.1091E+01*SB2-0.5217E+00*SB3 |
| 162 | A2= 0.3178E+01+0.3404E+01*SB -0.6563E+01*SB2+0.3287E+01*SB3 |
| 163 | A3=-0.8550E+00+0.1898E+01*SB -0.4134E+01*SB2+0.2052E+01*SB3 |
| 164 | A4= 0.1821E+01-0.5544E+01*SB +0.9966E+01*SB2-0.4592E+01*SB3 |
| 165 | A5= 0.0000E+00-0.2554E+01*SB +0.6574E+01*SB2-0.3145E+01*SB3 |
| 166 | goto 100 |
| 167 | C Ifl = 1 |
| 168 | 12 A0=Exp(-0.1314E+01-0.3117E-01*SB +0.1706E+00*SB2-0.1956E+00*SB3) |
| 169 | A1= 0.2780E+00+0.2593E-01*SB +0.1416E-01*SB2-0.1153E-01*SB3 |
| 170 | A2= 0.3672E+01+0.5254E+00*SB +0.4441E-01*SB2-0.1651E-02*SB3 |
| 171 | A3= 0.2963E+02-0.1997E+02*SB +0.5682E+01*SB2+0.6001E+00*SB3 |
| 172 | A4= 0.8070E+00-0.2918E+00*SB +0.5498E-01*SB2-0.7655E-02*SB3 |
| 173 | A5= 0.0000E+00+0.2347E+00*SB -0.1332E-01*SB2+0.2772E-01*SB3 |
| 174 | goto 100 |
| 175 | C Ifl = 0 |
| 176 | 13 A0=Exp(-0.1060E+00-0.1488E+01*SB -0.2040E+00*SB2+0.4899E-01*SB3) |
| 177 | A1=-0.2580E+00+0.1448E+00*SB -0.1548E+00*SB2+0.4470E-01*SB3 |
| 178 | A2= 0.5195E+01+0.9257E+00*SB -0.1613E+00*SB2-0.5502E-01*SB3 |
| 179 | A3= 0.5126E+01+0.2061E+01*SB -0.2862E+01*SB2+0.6837E+00*SB3 |
| 180 | A4= 0.1118E+01+0.2549E+00*SB -0.5194E+00*SB2+0.1831E+00*SB3 |
| 181 | A5= 0.0000E+00+0.2507E+01*SB -0.9946E+00*SB2+0.2547E+00*SB3 |
| 182 | goto 100 |
| 183 | C Ifl = -1 |
| 184 | 14 A0=Exp(-0.1141E+00-0.3683E+01*SB +0.4654E+01*SB2-0.2191E+01*SB3) |
| 185 | A1=-0.4269E+00-0.4266E-01*SB -0.1957E+00*SB2+0.1341E+00*SB3 |
| 186 | A2= 0.9960E+01-0.3425E+01*SB +0.4851E+01*SB2-0.2131E+01*SB3 |
| 187 | A3=-0.6381E+00+0.9855E+00*SB -0.1629E+01*SB2+0.7407E+00*SB3 |
| 188 | A4= 0.3091E-03-0.1846E+00*SB +0.2512E+00*SB2-0.1031E+00*SB3 |
| 189 | A5=-0.1310E+01+0.1654E+01*SB -0.1672E+01*SB2+0.8606E+00*SB3 |
| 190 | goto 100 |
| 191 | C Ifl = -2 |
| 192 | 15 A0=Exp(-0.2025E+01-0.4256E+01*SB +0.4815E+01*SB2-0.1907E+01*SB3) |
| 193 | A1=-0.2622E+00+0.2118E+00*SB -0.2479E+00*SB2+0.6284E-01*SB3 |
| 194 | A2= 0.9113E+01-0.2913E+01*SB +0.5533E+01*SB2-0.2438E+01*SB3 |
| 195 | A3= 0.2999E+02+0.2693E+02*SB -0.2000E+02*SB2-0.4426E+00*SB3 |
| 196 | A4= 0.1524E+01-0.9878E+00*SB +0.1582E+01*SB2-0.7228E+00*SB3 |
| 197 | A5=-0.2241E+00+0.3366E+01*SB -0.3198E+01*SB2+0.1113E+01*SB3 |
| 198 | goto 100 |
| 199 | C Ifl = -3 |
| 200 | 16 A0=Exp(-0.1855E+01-0.5294E+01*SB +0.7825E+00*SB2-0.5312E-01*SB3) |
| 201 | A1=-0.4000E-02-0.1379E+01*SB +0.1215E+01*SB2-0.6000E-02*SB3 |
| 202 | A2= 0.6870E+01+0.1088E+01*SB -0.5497E+00*SB2-0.1920E+00*SB3 |
| 203 | A3= 0.0000E+00+0.4829E+01*SB +0.8971E+01*SB2+0.2999E+02*SB3 |
| 204 | A4= 0.1000E-02+0.3239E+00*SB -0.2326E+00*SB2-0.2922E+00*SB3 |
| 205 | A5= 0.0000E+00-0.2450E+01*SB +0.3630E+01*SB2-0.1137E+01*SB3 |
| 206 | goto 100 |
| 207 | C Ifl = -4 |
| 208 | 17 A0=SB** 0.4376E+00*Exp(-0.6933E+01+0.2995E+01*SB -0.2304E+01*SB2) |
| 209 | A1=-0.4228E+00+0.1222E+01*SB -0.1971E+01*SB2+0.8951E+00*SB3 |
| 210 | A2= 0.6926E+01-0.1031E+02*SB +0.1468E+02*SB2-0.6065E+01*SB3 |
| 211 | A3= 0.8573E+01-0.1998E+02*SB +0.1965E+02*SB2-0.4859E+01*SB3 |
| 212 | A4= 0.9299E-01-0.4951E+00*SB +0.6898E+00*SB2-0.1069E+00*SB3 |
| 213 | A5=-0.5687E-01+0.5526E+01*SB -0.7217E+01*SB2+0.3569E+01*SB3 |
| 214 | goto 100 |
| 215 | C Ifl = -5 |
| 216 | 18 A0=SB** 0.2274E+00*Exp(-0.2564E+01-0.7256E+01*SB +0.4891E+01*SB2) |
| 217 | A1=-0.1423E+00-0.5363E+00*SB +0.9826E+00*SB2-0.6096E+00*SB3 |
| 218 | A2= 0.6238E+01-0.2785E+01*SB +0.5089E+00*SB2+0.2333E+01*SB3 |
| 219 | A3=-0.1029E+01+0.4684E+00*SB +0.3126E+01*SB2-0.2893E+01*SB3 |
| 220 | A4= 0.1729E-01+0.7293E-01*SB -0.5278E+00*SB2+0.3967E+00*SB3 |
| 221 | A5= 0.6208E+00+0.3942E+01*SB -0.3736E+01*SB2+0.1039E+01*SB3 |
| 222 | goto 100 |
| 223 | C Ifl = -6 |
| 224 | 19 A0=SB** 0.6663E+00*Exp(-0.7572E+01+0.4202E+01*SB -0.7464E+01*SB2) |
| 225 | A1=-0.1951E+00-0.1041E+01*SB +0.6688E+01*SB2-0.9976E+01*SB3 |
| 226 | A2= 0.5920E+01-0.7929E+01*SB +0.2997E+02*SB2+0.2991E+02*SB3 |
| 227 | A3= 0.1616E+01-0.9668E+01*SB +0.2995E+02*SB2+0.3000E+02*SB3 |
| 228 | A4= 0.4595E+00-0.3761E+01*SB +0.5663E+01*SB2+0.2066E+01*SB3 |
| 229 | A5= 0.2664E+01+0.1067E+00*SB -0.9717E+01*SB2+0.9980E+01*SB3 |
| 230 | goto 100 |
| 231 | |
| 232 | 2 Goto(21,22,23,24,25,26,27,28,29)Iflv |
| 233 | C Ifl = 2 |
| 234 | 21 A0=Exp( 0.2790E+00+0.1526E+01*SB -0.4306E+01*SB2+0.2060E+01*SB3) |
| 235 | A1= 0.5380E+00-0.4199E+00*SB +0.9582E+00*SB2-0.4624E+00*SB3 |
| 236 | A2= 0.3259E+01+0.3236E+01*SB -0.6130E+01*SB2+0.3073E+01*SB3 |
| 237 | A3=-0.8390E+00+0.2092E+01*SB -0.4528E+01*SB2+0.2251E+01*SB3 |
| 238 | A4= 0.1847E+01-0.5377E+01*SB +0.9994E+01*SB2-0.4699E+01*SB3 |
| 239 | A5= 0.0000E+00-0.2097E+01*SB +0.5742E+01*SB2-0.2774E+01*SB3 |
| 240 | goto 100 |
| 241 | C Ifl = 1 |
| 242 | 22 A0=Exp(-0.1318E+01+0.2387E-01*SB +0.4416E-01*SB2-0.1337E+00*SB3) |
| 243 | A1= 0.2760E+00+0.4302E-01*SB -0.1992E-01*SB2+0.2945E-02*SB3 |
| 244 | A2= 0.3660E+01+0.5319E+00*SB +0.2772E-01*SB2+0.1030E-01*SB3 |
| 245 | A3= 0.2910E+02-0.2000E+02*SB +0.6756E+01*SB2-0.1017E-01*SB3 |
| 246 | A4= 0.8010E+00-0.2702E+00*SB +0.1663E-01*SB2+0.8604E-02*SB3 |
| 247 | A5= 0.0000E+00+0.2822E+00*SB -0.1064E+00*SB2+0.6349E-01*SB3 |
| 248 | goto 100 |
| 249 | C Ifl = 0 |
| 250 | 23 A0=Exp(-0.1623E+01-0.1156E+01*SB +0.7802E+00*SB2-0.3593E-01*SB3) |
| 251 | A1=-0.5000E+00-0.2854E-01*SB +0.2212E+00*SB2-0.9403E-01*SB3 |
| 252 | A2= 0.3821E+01+0.8254E+00*SB +0.5031E+00*SB2-0.3367E+00*SB3 |
| 253 | A3= 0.5809E+01-0.3042E+01*SB -0.2726E+01*SB2+0.1371E+01*SB3 |
| 254 | A4= 0.4500E+00-0.6445E+00*SB +0.3019E+00*SB2-0.5642E-01*SB3 |
| 255 | A5= 0.0000E+00+0.8833E+00*SB -0.8495E-01*SB2+0.1565E+00*SB3 |
| 256 | goto 100 |
| 257 | C Ifl = -1 |
| 258 | 24 A0=Exp(-0.2693E-01-0.5523E+00*SB -0.1741E+01*SB2+0.6360E+00*SB3) |
| 259 | A1=-0.8634E+00+0.2368E+00*SB +0.8493E-01*SB2-0.6031E-01*SB3 |
| 260 | A2= 0.1200E+02+0.2398E+01*SB -0.6554E+01*SB2+0.2441E+01*SB3 |
| 261 | A3= 0.2528E+02+0.3000E+02*SB +0.2826E+02*SB2-0.2000E+02*SB3 |
| 262 | A4= 0.4000E+01-0.1239E+01*SB +0.2344E+01*SB2-0.5239E+00*SB3 |
| 263 | A5=-0.3466E+01+0.1731E+01*SB +0.1155E+01*SB2-0.5539E+00*SB3 |
| 264 | goto 100 |
| 265 | C Ifl = -2 |
| 266 | 25 A0=Exp(-0.1881E+01-0.2404E+01*SB +0.2946E+01*SB2-0.1553E+01*SB3) |
| 267 | A1=-0.6715E+00+0.4590E+00*SB -0.4233E+00*SB2+0.1678E+00*SB3 |
| 268 | A2= 0.9530E+01-0.2395E+01*SB +0.4102E+01*SB2-0.1921E+01*SB3 |
| 269 | A3= 0.1828E+02+0.3982E+00*SB -0.7874E+01*SB2+0.5451E+01*SB3 |
| 270 | A4= 0.1019E+01-0.1596E+00*SB +0.5041E+00*SB2-0.2586E+00*SB3 |
| 271 | A5=-0.1804E+01+0.3477E+01*SB -0.3082E+01*SB2+0.1389E+01*SB3 |
| 272 | goto 100 |
| 273 | C Ifl = -3 |
| 274 | 26 A0=Exp(-0.3180E+01-0.4460E+01*SB +0.2338E+01*SB2-0.6414E+00*SB3) |
| 275 | A1=-0.3873E+00-0.4457E+00*SB +0.1845E+00*SB2+0.2391E-01*SB3 |
| 276 | A2= 0.9287E+01-0.3698E+01*SB +0.2525E+01*SB2-0.8127E+00*SB3 |
| 277 | A3= 0.8968E+01+0.2853E+02*SB -0.9043E+01*SB2-0.1270E+01*SB3 |
| 278 | A4= 0.4911E-02-0.1111E+00*SB +0.6279E+00*SB2-0.3286E+00*SB3 |
| 279 | A5=-0.1892E+01+0.1020E+01*SB +0.1200E+01*SB2-0.5246E+00*SB3 |
| 280 | goto 100 |
| 281 | C Ifl = -4 |
| 282 | 27 A0=SB** 0.5061E+00*Exp(-0.6545E+01+0.2128E+01*SB -0.2682E+01*SB2) |
| 283 | A1=-0.6422E+00+0.1348E+01*SB -0.2112E+01*SB2+0.1023E+01*SB3 |
| 284 | A2= 0.8633E+01-0.1098E+02*SB +0.1297E+02*SB2-0.5481E+01*SB3 |
| 285 | A3= 0.1094E+02-0.1999E+02*SB +0.1569E+02*SB2+0.1640E+01*SB3 |
| 286 | A4= 0.6101E-01-0.4030E+00*SB +0.6944E+00*SB2-0.1071E+00*SB3 |
| 287 | A5=-0.1031E+01+0.6066E+01*SB -0.7177E+01*SB2+0.3941E+01*SB3 |
| 288 | goto 100 |
| 289 | C Ifl = -5 |
| 290 | 28 A0=SB** 0.9548E+00*Exp(-0.3509E+01-0.4622E+01*SB +0.2158E+01*SB2) |
| 291 | A1=-0.3843E+00+0.3255E-01*SB +0.1015E+01*SB2-0.8073E+00*SB3 |
| 292 | A2= 0.7119E+01-0.6982E+01*SB +0.1057E+02*SB2-0.5042E+01*SB3 |
| 293 | A3=-0.2830E+00+0.6828E+00*SB +0.3744E+01*SB2-0.3577E+01*SB3 |
| 294 | A4=-0.7753E-02-0.6807E-01*SB -0.5997E+00*SB2+0.5024E+00*SB3 |
| 295 | A5= 0.5160E+00+0.2897E+01*SB -0.1816E+01*SB2+0.5683E+00*SB3 |
| 296 | goto 100 |
| 297 | C Ifl = -6 |
| 298 | 29 A0=SB** 0.6174E+00*Exp(-0.7663E+01+0.7381E+01*SB +0.1136E+01*SB2) |
| 299 | A1=-0.8717E-01-0.2683E+01*SB +0.9103E+01*SB2-0.7527E+01*SB3 |
| 300 | A2= 0.6726E+01-0.9855E+01*SB +0.2906E+02*SB2+0.2793E+02*SB3 |
| 301 | A3= 0.1893E+01-0.1341E+02*SB +0.1923E+02*SB2-0.4600E+01*SB3 |
| 302 | A4=-0.2237E+00+0.1519E+01*SB +0.3735E+00*SB2-0.7050E+01*SB3 |
| 303 | A5= 0.1873E+01+0.2195E+01*SB -0.8159E+01*SB2+0.6019E+01*SB3 |
| 304 | goto 100 |
| 305 | |
| 306 | 3 Goto(31,32,33,34,35,36,37,38,39)Iflv |
| 307 | C Ifl = 2 |
| 308 | 31 A0=Exp( 0.1670E+00+0.1075E+01*SB -0.3383E+01*SB2+0.1614E+01*SB3) |
| 309 | A1= 0.5080E+00-0.3170E+00*SB +0.7692E+00*SB2-0.3742E+00*SB3 |
| 310 | A2= 0.3243E+01+0.4730E+01*SB -0.8685E+01*SB2+0.4224E+01*SB3 |
| 311 | A3=-0.8320E+00+0.5206E+01*SB -0.1064E+02*SB2+0.5154E+01*SB3 |
| 312 | A4= 0.2189E+01-0.1541E+01*SB +0.3218E+01*SB2-0.1573E+01*SB3 |
| 313 | A5= 0.0000E+00-0.1513E+01*SB +0.4587E+01*SB2-0.2226E+01*SB3 |
| 314 | goto 100 |
| 315 | C Ifl = 1 |
| 316 | 32 A0=Exp(-0.1344E+01-0.1748E-01*SB +0.1351E+00*SB2-0.1933E+00*SB3) |
| 317 | A1= 0.2760E+00+0.3213E-01*SB +0.4936E-02*SB2-0.9453E-02*SB3 |
| 318 | A2= 0.3660E+01+0.5253E+00*SB +0.4304E-01*SB2+0.5047E-03*SB3 |
| 319 | A3= 0.2982E+02-0.2000E+02*SB +0.5957E+01*SB2+0.6646E+00*SB3 |
| 320 | A4= 0.7950E+00-0.2861E+00*SB +0.5038E-01*SB2-0.5960E-02*SB3 |
| 321 | A5= 0.0000E+00+0.2441E+00*SB -0.2216E-01*SB2+0.2776E-01*SB3 |
| 322 | goto 100 |
| 323 | C Ifl = 0 |
| 324 | 33 A0=Exp( 0.1116E+01-0.3529E+01*SB +0.1456E+01*SB2-0.1419E+00*SB3) |
| 325 | A1= 0.0000E+00-0.5181E+00*SB +0.5405E+00*SB2-0.2205E+00*SB3 |
| 326 | A2= 0.6529E+01-0.1407E+01*SB +0.5685E-01*SB2+0.3105E+00*SB3 |
| 327 | A3= 0.2643E+01-0.5811E+01*SB +0.4343E+01*SB2-0.1257E+01*SB3 |
| 328 | A4= 0.2221E+01-0.3929E+01*SB +0.2208E+01*SB2-0.3916E+00*SB3 |
| 329 | A5= 0.0000E+00+0.2176E+01*SB -0.8498E+00*SB2+0.2041E+00*SB3 |
| 330 | goto 100 |
| 331 | C Ifl = -1 |
| 332 | 34 A0=Exp(-0.6345E+00-0.3473E+01*SB +0.3619E+01*SB2-0.1521E+01*SB3) |
| 333 | A1=-0.2746E-01-0.2354E+00*SB -0.9563E-01*SB2+0.9076E-01*SB3 |
| 334 | A2= 0.9349E+01-0.2734E+01*SB +0.4205E+01*SB2-0.1868E+01*SB3 |
| 335 | A3=-0.4148E+00+0.1275E+01*SB -0.1762E+01*SB2+0.6678E+00*SB3 |
| 336 | A4=-0.6941E-02-0.3389E-01*SB +0.2910E-01*SB2-0.7428E-03*SB3 |
| 337 | A5=-0.2328E+00+0.1099E+01*SB -0.1257E+01*SB2+0.6759E+00*SB3 |
| 338 | goto 100 |
| 339 | C Ifl = -2 |
| 340 | 35 A0=Exp(-0.1043E+01-0.4115E+01*SB +0.4084E+01*SB2-0.1564E+01*SB3) |
| 341 | A1=-0.4823E-01+0.5433E-02*SB -0.3178E+00*SB2+0.1413E+00*SB3 |
| 342 | A2= 0.1038E+02-0.3667E+01*SB +0.4441E+01*SB2-0.1471E+01*SB3 |
| 343 | A3= 0.3729E+02-0.1991E+02*SB +0.2565E+01*SB2+0.3049E+01*SB3 |
| 344 | A4= 0.1831E+01-0.1349E+01*SB +0.1416E+01*SB2-0.4872E+00*SB3 |
| 345 | A5=-0.3255E+00+0.3095E+01*SB -0.3465E+01*SB2+0.1337E+01*SB3 |
| 346 | goto 100 |
| 347 | C Ifl = -3 |
| 348 | 36 A0=Exp(-0.2887E+01-0.4406E+01*SB +0.1392E+01*SB2-0.2654E+00*SB3) |
| 349 | A1=-0.1081E+00-0.1614E+00*SB -0.1005E+00*SB2+0.5118E-01*SB3 |
| 350 | A2= 0.7105E+01-0.1613E+01*SB +0.1555E+01*SB2-0.4143E+00*SB3 |
| 351 | A3= 0.9845E+00+0.7054E+01*SB +0.9067E+01*SB2+0.4381E+00*SB3 |
| 352 | A4= 0.4347E-01-0.4131E-01*SB +0.4727E+00*SB2-0.2627E+00*SB3 |
| 353 | A5= 0.1497E+00+0.7286E+00*SB +0.1079E+01*SB2-0.6934E+00*SB3 |
| 354 | goto 100 |
| 355 | C Ifl = -4 |
| 356 | 37 A0=SB** 0.5386E+00*Exp(-0.6447E+01+0.1532E+01*SB -0.2501E+01*SB2) |
| 357 | A1=-0.2340E+00+0.1158E+01*SB -0.2199E+01*SB2+0.1110E+01*SB3 |
| 358 | A2= 0.9582E+01-0.1232E+02*SB +0.1417E+02*SB2-0.5780E+01*SB3 |
| 359 | A3= 0.1136E+02-0.1977E+02*SB +0.1412E+02*SB2+0.4260E+01*SB3 |
| 360 | A4= 0.4897E-01-0.3630E+00*SB +0.6748E+00*SB2-0.8826E-01*SB3 |
| 361 | A5= 0.1638E+00+0.5938E+01*SB -0.7697E+01*SB2+0.4287E+01*SB3 |
| 362 | goto 100 |
| 363 | C Ifl = -5 |
| 364 | 38 A0=SB** 0.6821E+00*Exp(-0.3091E+01-0.6402E+01*SB +0.3642E+01*SB2) |
| 365 | A1=-0.6025E-01-0.3351E+00*SB +0.6665E+00*SB2-0.4328E+00*SB3 |
| 366 | A2= 0.7107E+01-0.2327E+01*SB +0.1482E+01*SB2+0.6384E+00*SB3 |
| 367 | A3=-0.9249E+00+0.1607E+01*SB +0.3734E+01*SB2-0.3985E+01*SB3 |
| 368 | A4= 0.5565E-01+0.5489E-01*SB -0.5636E+00*SB2+0.4372E+00*SB3 |
| 369 | A5= 0.1252E+01+0.2915E+01*SB -0.2690E+01*SB2+0.8944E+00*SB3 |
| 370 | goto 100 |
| 371 | C Ifl = -6 |
| 372 | 39 A0=SB** 0.5589E+00*Exp(-0.8150E+01+0.5395E+01*SB -0.6026E+00*SB2) |
| 373 | A1=-0.1062E-01-0.2427E+01*SB +0.8547E+01*SB2-0.6318E+01*SB3 |
| 374 | A2= 0.6721E+01-0.8373E+01*SB +0.2805E+02*SB2+0.2977E+02*SB3 |
| 375 | A3= 0.1926E+01-0.5905E+01*SB +0.2880E+01*SB2-0.1004E+02*SB3 |
| 376 | A4=-0.1701E+00+0.1258E+01*SB +0.4202E+01*SB2-0.9911E+01*SB3 |
| 377 | A5= 0.2429E+01+0.2037E+01*SB -0.2878E+01*SB2+0.3858E+01*SB3 |
| 378 | goto 100 |
| 379 | |
| 380 | 4 Goto(41,42,43,44,45,46,47,48,49)Iflv |
| 381 | C Ifl = 2 |
| 382 | 41 A0=Exp( 0.3760E+00+0.8053E+00*SB -0.2513E+01*SB2+0.1042E+01*SB3) |
| 383 | A1= 0.5650E+00-0.2477E+00*SB +0.5156E+00*SB2-0.2174E+00*SB3 |
| 384 | A2= 0.3464E+01+0.4958E+01*SB -0.8038E+01*SB2+0.3416E+01*SB3 |
| 385 | A3=-0.5850E+00+0.7780E+01*SB -0.1388E+02*SB2+0.5883E+01*SB3 |
| 386 | A4= 0.2322E+01-0.1325E+01*SB +0.2338E+01*SB2-0.9912E+00*SB3 |
| 387 | A5= 0.0000E+00-0.1051E+01*SB +0.3148E+01*SB2-0.1331E+01*SB3 |
| 388 | goto 100 |
| 389 | C Ifl = 1 |
| 390 | 42 A0=Exp(-0.1324E+01+0.4762E-01*SB -0.9509E-01*SB2-0.1792E-01*SB3) |
| 391 | A1= 0.2890E+00+0.6535E-01*SB -0.3082E-01*SB2+0.1996E-03*SB3 |
| 392 | A2= 0.3580E+01+0.5494E+00*SB -0.6856E-01*SB2+0.6323E-01*SB3 |
| 393 | A3= 0.3021E+02-0.2000E+02*SB +0.7011E+01*SB2-0.7697E+00*SB3 |
| 394 | A4= 0.7990E+00-0.2426E+00*SB +0.8071E-03*SB2+0.4460E-02*SB3 |
| 395 | A5= 0.0000E+00+0.3810E+00*SB -0.8815E-01*SB2-0.3724E-02*SB3 |
| 396 | goto 100 |
| 397 | C Ifl = 0 |
| 398 | 43 A0=Exp(-0.1920E+00-0.1130E+01*SB -0.2375E+00*SB2+0.4022E-01*SB3) |
| 399 | A1=-0.2120E+00+0.6095E-01*SB -0.9243E-01*SB2+0.2470E-01*SB3 |
| 400 | A2= 0.4549E+01+0.9281E+00*SB +0.4614E-01*SB2-0.4467E-01*SB3 |
| 401 | A3= 0.1197E+02-0.3295E+01*SB -0.1085E+01*SB2+0.4984E+00*SB3 |
| 402 | A4= 0.1616E+01-0.1152E+00*SB -0.3057E+00*SB2+0.1216E+00*SB3 |
| 403 | A5= 0.0000E+00+0.2254E+01*SB -0.8923E+00*SB2+0.2191E+00*SB3 |
| 404 | goto 100 |
| 405 | C Ifl = -1 |
| 406 | 44 A0=Exp(-0.1409E+01-0.8736E+00*SB -0.3396E+00*SB2+0.1329E+00*SB3) |
| 407 | A1=-0.2732E+00-0.7161E-01*SB +0.8535E-01*SB2-0.4319E-01*SB3 |
| 408 | A2= 0.9948E+01-0.2571E+01*SB +0.5435E+00*SB2+0.1048E+00*SB3 |
| 409 | A3= 0.3191E+02-0.6723E+01*SB -0.2000E+02*SB2+0.9788E+01*SB3 |
| 410 | A4= 0.3999E+01+0.2953E+00*SB -0.9248E+00*SB2+0.5909E+00*SB3 |
| 411 | A5=-0.6366E+00+0.6347E+00*SB +0.5611E+00*SB2-0.2659E+00*SB3 |
| 412 | goto 100 |
| 413 | C Ifl = -2 |
| 414 | 45 A0=Exp(-0.1935E+01-0.3640E+01*SB +0.3418E+01*SB2-0.1162E+01*SB3) |
| 415 | A1=-0.2333E+00+0.3532E+00*SB -0.4772E+00*SB2+0.1515E+00*SB3 |
| 416 | A2= 0.1009E+02-0.4327E+01*SB +0.4977E+01*SB2-0.1588E+01*SB3 |
| 417 | A3= 0.4998E+02-0.1658E+02*SB -0.1299E+01*SB2+0.1680E+01*SB3 |
| 418 | A4= 0.1717E+01-0.8842E+00*SB +0.8121E+00*SB2-0.2623E+00*SB3 |
| 419 | A5=-0.1852E+00+0.3591E+01*SB -0.3439E+01*SB2+0.1119E+01*SB3 |
| 420 | goto 100 |
| 421 | C Ifl = -3 |
| 422 | 46 A0=Exp(-0.1864E+01-0.3708E+01*SB -0.5756E+00*SB2+0.4375E+00*SB3) |
| 423 | A1= 0.1000E-02-0.1191E+01*SB +0.1129E+01*SB2-0.1823E+00*SB3 |
| 424 | A2= 0.6898E+01+0.7131E+00*SB -0.1816E+01*SB2+0.7422E+00*SB3 |
| 425 | A3= 0.0000E+00+0.3461E+01*SB +0.5716E+01*SB2+0.1099E+02*SB3 |
| 426 | A4= 0.1000E-02+0.2346E+00*SB -0.2442E+00*SB2-0.8945E-01*SB3 |
| 427 | A5= 0.0000E+00-0.2307E+01*SB +0.3654E+01*SB2-0.1229E+01*SB3 |
| 428 | goto 100 |
| 429 | C Ifl = -4 |
| 430 | 47 A0=SB** 0.5944E+00*Exp(-0.6580E+01+0.1186E+01*SB -0.8149E+00*SB2) |
| 431 | A1=-0.4012E+00+0.1013E+01*SB -0.1345E+01*SB2+0.4858E+00*SB3 |
| 432 | A2= 0.5745E+01-0.1001E+02*SB +0.1332E+02*SB2-0.4777E+01*SB3 |
| 433 | A3= 0.9501E+01-0.2000E+02*SB +0.2158E+02*SB2-0.7159E+01*SB3 |
| 434 | A4= 0.1330E+00-0.5430E+00*SB +0.6810E+00*SB2-0.2441E+00*SB3 |
| 435 | A5= 0.7665E-01+0.4725E+01*SB -0.5315E+01*SB2+0.1948E+01*SB3 |
| 436 | goto 100 |
| 437 | C Ifl = -5 |
| 438 | 48 A0=SB** 0.1171E+00*Exp(-0.4100E+01-0.4162E+01*SB +0.2139E+01*SB2) |
| 439 | A1=-0.1266E+00-0.3860E+00*SB +0.5611E+00*SB2-0.2887E+00*SB3 |
| 440 | A2= 0.4416E+01-0.2970E+00*SB +0.8367E-01*SB2+0.5844E+00*SB3 |
| 441 | A3=-0.1007E+01+0.7322E+00*SB +0.4878E+01*SB2-0.3640E+01*SB3 |
| 442 | A4= 0.1202E-01+0.2368E+00*SB -0.6270E+00*SB2+0.3810E+00*SB3 |
| 443 | A5= 0.1117E+01+0.2026E+01*SB -0.1649E+01*SB2+0.5582E+00*SB3 |
| 444 | goto 100 |
| 445 | C Ifl = -6 |
| 446 | 49 A0=SB** 0.4778E+00*Exp(-0.8672E+01+0.7464E+01*SB -0.2952E+01*SB2) |
| 447 | A1=-0.2503E-01-0.2713E+01*SB +0.9002E+01*SB2-0.6336E+01*SB3 |
| 448 | A2= 0.6266E+01-0.1191E+02*SB +0.2982E+02*SB2-0.7350E+01*SB3 |
| 449 | A3= 0.2402E+01-0.7116E+01*SB -0.9642E+01*SB2+0.2415E+02*SB3 |
| 450 | A4=-0.2056E+00+0.1537E+01*SB +0.2850E+01*SB2-0.9850E+01*SB3 |
| 451 | A5= 0.2237E+01+0.2495E+01*SB -0.6424E+01*SB2+0.9998E+01*SB3 |
| 452 | goto 100 |
| 453 | |
| 454 | 5 Goto(51,52,53,54,55,56,57,58,59)Iflv |
| 455 | C Ifl = 2 |
| 456 | 51 A0=Exp( 0.8000E-01+0.7364E+00*SB -0.2714E+01*SB2+0.1311E+01*SB3) |
| 457 | A1= 0.4930E+00-0.2001E+00*SB +0.5784E+00*SB2-0.2915E+00*SB3 |
| 458 | A2= 0.3001E+01+0.3538E+01*SB -0.6155E+01*SB2+0.3083E+01*SB3 |
| 459 | A3=-0.1000E+01+0.3871E+01*SB -0.8334E+01*SB2+0.4219E+01*SB3 |
| 460 | A4= 0.2986E+01+0.1597E+01*SB -0.3368E+01*SB2+0.1644E+01*SB3 |
| 461 | A5= 0.0000E+00-0.9256E+00*SB +0.3570E+01*SB2-0.1777E+01*SB3 |
| 462 | goto 100 |
| 463 | C Ifl = 1 |
| 464 | 52 A0=Exp(-0.1806E+01-0.6672E-01*SB -0.2605E+00*SB2+0.2341E-01*SB3) |
| 465 | A1= 0.1750E+00+0.3872E-01*SB -0.2189E-01*SB2+0.1415E-01*SB3 |
| 466 | A2= 0.3322E+01+0.7786E+00*SB -0.2902E+00*SB2+0.1517E+00*SB3 |
| 467 | A3= 0.4414E+02-0.1987E+02*SB +0.2597E+01*SB2+0.2670E+01*SB3 |
| 468 | A4= 0.9610E+00-0.2864E+00*SB -0.5524E-01*SB2+0.6229E-01*SB3 |
| 469 | A5= 0.0000E+00+0.2658E+00*SB -0.4728E-02*SB2+0.6048E-01*SB3 |
| 470 | goto 100 |
| 471 | C Ifl = 0 |
| 472 | 53 A0=Exp(-0.6510E+00-0.1128E+01*SB -0.6239E-01*SB2-0.8838E-01*SB3) |
| 473 | A1=-0.2590E+00+0.1822E+00*SB -0.2682E+00*SB2+0.9422E-01*SB3 |
| 474 | A2= 0.4607E+01+0.7792E+00*SB +0.8937E+00*SB2-0.5553E+00*SB3 |
| 475 | A3= 0.1627E+02-0.1114E+02*SB +0.4928E+01*SB2-0.1715E+01*SB3 |
| 476 | A4= 0.1236E+01+0.1945E+00*SB -0.3297E+00*SB2+0.6489E-01*SB3 |
| 477 | A5= 0.0000E+00+0.3346E+01*SB -0.2337E+01*SB2+0.7850E+00*SB3 |
| 478 | goto 100 |
| 479 | C Ifl = -1 |
| 480 | 54 A0=Exp(-0.1508E+01-0.5560E+00*SB -0.3523E+00*SB2+0.6562E-01*SB3) |
| 481 | A1=-0.3223E+00+0.2095E-01*SB -0.2049E-02*SB2-0.3475E-01*SB3 |
| 482 | A2= 0.9469E+01-0.3923E+01*SB +0.4333E+01*SB2-0.1654E+01*SB3 |
| 483 | A3= 0.1646E+02-0.1082E+02*SB +0.8941E+01*SB2-0.5494E+01*SB3 |
| 484 | A4= 0.2908E+01+0.2162E+01*SB -0.3233E+01*SB2+0.1267E+01*SB3 |
| 485 | A5=-0.5819E+00+0.3914E+00*SB +0.6460E+00*SB2-0.3239E+00*SB3 |
| 486 | goto 100 |
| 487 | C Ifl = -2 |
| 488 | 55 A0=Exp(-0.1951E+01-0.3435E+01*SB +0.3424E+01*SB2-0.1249E+01*SB3) |
| 489 | A1=-0.2942E+00+0.4408E+00*SB -0.5453E+00*SB2+0.1552E+00*SB3 |
| 490 | A2= 0.9782E+01-0.3454E+01*SB +0.4510E+01*SB2-0.1649E+01*SB3 |
| 491 | A3= 0.4999E+02-0.1993E+02*SB -0.2039E+01*SB2+0.5694E+00*SB3 |
| 492 | A4= 0.1938E+01-0.1351E+01*SB +0.1386E+01*SB2-0.5324E+00*SB3 |
| 493 | A5=-0.2410E+00+0.3434E+01*SB -0.3334E+01*SB2+0.1067E+01*SB3 |
| 494 | goto 100 |
| 495 | C Ifl = -3 |
| 496 | 56 A0=Exp(-0.1804E+01-0.4381E+01*SB -0.3699E+00*SB2+0.3878E+00*SB3) |
| 497 | A1=-0.1000E-02-0.9334E+00*SB +0.7156E+00*SB2-0.2029E+00*SB3 |
| 498 | A2= 0.6896E+01+0.2462E+01*SB -0.2885E+01*SB2+0.8701E+00*SB3 |
| 499 | A3= 0.0000E+00+0.5589E+01*SB +0.1047E+02*SB2+0.3000E+02*SB3 |
| 500 | A4= 0.1000E-02-0.5600E-02*SB +0.5618E-02*SB2+0.6598E-02*SB3 |
| 501 | A5= 0.0000E+00-0.3151E+01*SB +0.4025E+01*SB2-0.1232E+01*SB3 |
| 502 | goto 100 |
| 503 | C Ifl = -4 |
| 504 | 57 A0=SB** 0.7860E+00*Exp(-0.5041E+01-0.3357E+00*SB -0.4718E+00*SB2) |
| 505 | A1=-0.4989E+00+0.9571E+00*SB -0.1359E+01*SB2+0.5384E+00*SB3 |
| 506 | A2= 0.5986E+01-0.8541E+01*SB +0.1274E+02*SB2-0.5275E+01*SB3 |
| 507 | A3= 0.8121E+01-0.1753E+02*SB +0.2194E+02*SB2-0.8538E+01*SB3 |
| 508 | A4= 0.9290E-01-0.4390E+00*SB +0.6162E+00*SB2-0.2231E+00*SB3 |
| 509 | A5=-0.1257E+01+0.5677E+01*SB -0.5977E+01*SB2+0.2387E+01*SB3 |
| 510 | goto 100 |
| 511 | C Ifl = -5 |
| 512 | 58 A0=SB** 0.4537E+00*Exp(-0.3269E+01-0.5398E+01*SB +0.2893E+01*SB2) |
| 513 | A1=-0.1977E+00-0.4126E+00*SB +0.7058E+00*SB2-0.4038E+00*SB3 |
| 514 | A2= 0.4522E+01+0.6167E-01*SB -0.1849E+00*SB2+0.7345E+00*SB3 |
| 515 | A3=-0.1003E+01+0.1531E+01*SB +0.4515E+01*SB2-0.4368E+01*SB3 |
| 516 | A4= 0.3579E-01+0.1919E+00*SB -0.7268E+00*SB2+0.5192E+00*SB3 |
| 517 | A5= 0.5129E+00+0.2447E+01*SB -0.1989E+01*SB2+0.7529E+00*SB3 |
| 518 | goto 100 |
| 519 | C Ifl = -6 |
| 520 | 59 A0=SB** 0.7178E+00*Exp(-0.7327E+01+0.2277E+01*SB +0.3913E+01*SB2) |
| 521 | A1=-0.9842E-01-0.2362E+01*SB +0.8851E+01*SB2-0.7208E+01*SB3 |
| 522 | A2= 0.5552E+01-0.8935E+01*SB +0.2676E+02*SB2-0.1344E+02*SB3 |
| 523 | A3= 0.1593E+01-0.3505E+01*SB -0.1234E+01*SB2-0.1867E+02*SB3 |
| 524 | A4=-0.1723E+00+0.1530E+01*SB +0.2323E+01*SB2-0.9344E+01*SB3 |
| 525 | A5= 0.2081E+01+0.1939E+01*SB -0.3273E+01*SB2+0.9935E+01*SB3 |
| 526 | goto 100 |
| 527 | |
| 528 | 6 Goto(61,62,63,64,65,66,67,68,69)Iflv |
| 529 | C Ifl = 2 |
| 530 | 61 A0=Exp( 0.1590E+00+0.8361E+00*SB -0.2713E+01*SB2+0.1231E+01*SB3) |
| 531 | A1= 0.5110E+00-0.2081E+00*SB +0.5090E+00*SB2-0.2365E+00*SB3 |
| 532 | A2= 0.3158E+01+0.4946E+01*SB -0.8738E+01*SB2+0.4063E+01*SB3 |
| 533 | A3=-0.6000E+00+0.7591E+01*SB -0.1479E+02*SB2+0.6851E+01*SB3 |
| 534 | A4= 0.2306E+01-0.1913E+01*SB +0.3615E+01*SB2-0.1663E+01*SB3 |
| 535 | A5= 0.0000E+00-0.9682E+00*SB +0.3247E+01*SB2-0.1503E+01*SB3 |
| 536 | goto 100 |
| 537 | C Ifl = 1 |
| 538 | 62 A0=Exp(-0.1182E+01+0.1223E+00*SB +0.6512E-01*SB2-0.1208E+00*SB3) |
| 539 | A1= 0.2540E+00+0.2297E-01*SB -0.4249E-02*SB2+0.1143E-03*SB3 |
| 540 | A2= 0.3442E+01+0.5458E+00*SB +0.3803E-01*SB2-0.1598E-02*SB3 |
| 541 | A3= 0.2545E+02-0.2000E+02*SB +0.7516E+01*SB2-0.7726E+00*SB3 |
| 542 | A4= 0.9170E+00-0.3214E+00*SB +0.4666E-01*SB2-0.5126E-02*SB3 |
| 543 | A5= 0.0000E+00+0.1379E+00*SB -0.4579E-01*SB2+0.4988E-01*SB3 |
| 544 | goto 100 |
| 545 | C Ifl = 0 |
| 546 | 63 A0=Exp(-0.3410E+00-0.9586E+00*SB -0.4561E+00*SB2+0.1041E+00*SB3) |
| 547 | A1=-0.2400E+00+0.1498E+00*SB -0.1717E+00*SB2+0.4922E-01*SB3 |
| 548 | A2= 0.4841E+01+0.8606E+00*SB +0.1611E+00*SB2-0.1785E+00*SB3 |
| 549 | A3= 0.7427E+01-0.2102E+01*SB +0.1371E+01*SB2-0.8831E+00*SB3 |
| 550 | A4= 0.9600E+00+0.3385E+00*SB -0.3549E+00*SB2+0.8291E-01*SB3 |
| 551 | A5= 0.0000E+00+0.2516E+01*SB -0.1115E+01*SB2+0.2961E+00*SB3 |
| 552 | goto 100 |
| 553 | C Ifl = -1 |
| 554 | 64 A0=Exp( 0.3709E-01-0.3427E+01*SB +0.4188E+01*SB2-0.1910E+01*SB3) |
| 555 | A1=-0.3664E+00+0.1406E+00*SB -0.1340E+00*SB2+0.2966E-01*SB3 |
| 556 | A2= 0.1093E+02-0.3215E+01*SB +0.1803E+01*SB2-0.3219E+00*SB3 |
| 557 | A3=-0.7324E+00+0.9441E+00*SB -0.1698E+01*SB2+0.7797E+00*SB3 |
| 558 | A4= 0.7923E-01-0.5296E-01*SB +0.9541E-01*SB2-0.4648E-01*SB3 |
| 559 | A5=-0.1426E+01+0.2324E+01*SB -0.1860E+01*SB2+0.7350E+00*SB3 |
| 560 | goto 100 |
| 561 | C Ifl = -2 |
| 562 | 65 A0=Exp(-0.2396E+01-0.2506E+01*SB +0.3192E+01*SB2-0.1398E+01*SB3) |
| 563 | A1=-0.2366E+00+0.3132E+00*SB -0.5087E+00*SB2+0.1892E+00*SB3 |
| 564 | A2= 0.9013E+01-0.2295E+01*SB +0.4380E+01*SB2-0.1909E+01*SB3 |
| 565 | A3= 0.4999E+02-0.1713E+02*SB -0.1184E+01*SB2+0.5404E+00*SB3 |
| 566 | A4= 0.1739E+01-0.7946E+00*SB +0.9877E+00*SB2-0.4258E+00*SB3 |
| 567 | A5= 0.3889E-01+0.2840E+01*SB -0.3389E+01*SB2+0.1358E+01*SB3 |
| 568 | goto 100 |
| 569 | C Ifl = -3 |
| 570 | 66 A0=Exp(-0.1966E+01-0.4423E+01*SB -0.2825E+00*SB2+0.4253E+00*SB3) |
| 571 | A1=-0.4000E-02-0.1327E+01*SB +0.1367E+01*SB2-0.4343E+00*SB3 |
| 572 | A2= 0.6902E+01+0.2127E+01*SB -0.3095E+01*SB2+0.1140E+01*SB3 |
| 573 | A3= 0.0000E+00+0.4494E+01*SB +0.1085E+02*SB2+0.2149E+02*SB3 |
| 574 | A4= 0.1000E-02+0.3120E+00*SB -0.3684E+00*SB2+0.1059E+00*SB3 |
| 575 | A5= 0.0000E+00-0.2796E+01*SB +0.4415E+01*SB2-0.1559E+01*SB3 |
| 576 | goto 100 |
| 577 | C Ifl = -4 |
| 578 | 67 A0=SB** 0.5299E+00*Exp(-0.6308E+01+0.2148E+01*SB -0.1356E+01*SB2) |
| 579 | A1=-0.4240E+00+0.1047E+01*SB -0.1558E+01*SB2+0.6376E+00*SB3 |
| 580 | A2= 0.6562E+01-0.9971E+01*SB +0.1510E+02*SB2-0.6222E+01*SB3 |
| 581 | A3= 0.8340E+01-0.1952E+02*SB +0.2047E+02*SB2-0.6962E+01*SB3 |
| 582 | A4= 0.9903E-01-0.4558E+00*SB +0.5885E+00*SB2-0.1425E+00*SB3 |
| 583 | A5=-0.3086E+00+0.5061E+01*SB -0.6202E+01*SB2+0.2690E+01*SB3 |
| 584 | goto 100 |
| 585 | C Ifl = -5 |
| 586 | 68 A0=SB** 0.2033E+00*Exp(-0.2686E+01-0.5308E+01*SB +0.1509E+01*SB2) |
| 587 | A1=-0.1704E+00-0.2444E+00*SB +0.3557E+00*SB2-0.2174E+00*SB3 |
| 588 | A2= 0.5074E+01+0.1516E+01*SB -0.3815E+01*SB2+0.2956E+01*SB3 |
| 589 | A3=-0.1024E+01+0.6055E+00*SB +0.1338E+01*SB2+0.1855E+01*SB3 |
| 590 | A4= 0.9655E-02+0.1397E+00*SB -0.3254E+00*SB2+0.2057E+00*SB3 |
| 591 | A5= 0.5681E+00+0.2772E+01*SB -0.1963E+01*SB2+0.4607E+00*SB3 |
| 592 | goto 100 |
| 593 | C Ifl = -6 |
| 594 | 69 A0=SB** 0.8044E+00*Exp(-0.6946E+01+0.3617E+01*SB -0.9973E+01*SB2) |
| 595 | A1=-0.1901E+00-0.8342E+00*SB +0.4663E+01*SB2-0.9999E+01*SB3 |
| 596 | A2= 0.6194E+01-0.9668E+01*SB +0.2997E+02*SB2+0.2987E+02*SB3 |
| 597 | A3= 0.1173E+01-0.8138E+01*SB +0.2924E+02*SB2-0.2860E+00*SB3 |
| 598 | A4= 0.4761E+00-0.3683E+01*SB +0.4966E+01*SB2+0.9593E+01*SB3 |
| 599 | A5= 0.2691E+01-0.9906E+00*SB -0.4504E+01*SB2+0.9999E+01*SB3 |
| 600 | goto 100 |
| 601 | |
| 602 | 311 stop 'This option is not currently supported.' |
| 603 | |
| 604 | 100 CONTINUE |
| 605 | P012 = A0 *(X**A1) *((1.-X)**A2) |
| 606 | P34 = (1.+A3*(X**A4)) |
| 607 | P5 = (log(1.+1./X))**A5 |
| 608 | C |
| 609 | PDFD = P012 * P34 * P5 |
| 610 | PDF = PDFD |
| 611 | C |
| 612 | Ctq2OPd = PDF |
| 613 | if(Ctq2OPd.lt.0.0) then |
| 614 | Ctq2OPd = 0.0 |
| 615 | Irt=1 |
| 616 | endif |
| 617 | |
| 618 | Ist = Iset |
| 619 | |
| 620 | Lp = Iprtn |
| 621 | Qsto = QQ |
| 622 | |
| 623 | Return |
| 624 | C ----------------------- |
| 625 | ENTRY Wlamd2O (Iset, Iorder, Neff) |
| 626 | |
| 627 | Iorder = IORD (Iset) |
| 628 | Wlamd2O = VLM (Neff, Iset) |
| 629 | |
| 630 | Ctq2OPd = 0.0 |
| 631 | RETURN |
| 632 | C ----------------------- |
| 633 | Entry PrCtq2O |
| 634 | > (Iset, Iordr, Ischeme, MxFlv, |
| 635 | > Alam4, Alam5, Alam6, Amas4, Amas5, Amas6, |
| 636 | > Xmin, Qini, Qmax, ExpNor) |
| 637 | |
| 638 | C Return QCD parameters and Fitting parameters |
| 639 | C associated with parton distribution set Iset. |
| 640 | C Iord : Order Of Fit |
| 641 | C Ischeme : (0, 1, 2) for (LO, MS-bar-NLO, DIS-NLO) resp. |
| 642 | C MxFlv : Maximum number of flavors included |
| 643 | C Alam_i : i = 4,5,6 Effective lambda for i-flavors |
| 644 | |
| 645 | C Amas_i : i = 4,5,6 Mass parameter for flavor i |
| 646 | C Xmin, Qini, Qmax : self explanary |
| 647 | C ExpNor(I) : Normalization factor for the experimental data set used in |
| 648 | C obtaining the best global fit for parton distributions Iset: |
| 649 | C I = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 |
| 650 | C BCDMS NMC90 NMC280 CCFR E605 WA70 E706 UA6 H1 ZEUS |
| 651 | |
| 652 | Iordr = Iord (Iset) |
| 653 | Ischeme= Isch (Iset) |
| 654 | MxFlv = Nqrk (Iset) |
| 655 | |
| 656 | Alam4 = Vlm(4,Iset) |
| 657 | Alam5 = Vlm(5,Iset) |
| 658 | Alam6 = Vlm(6,Iset) |
| 659 | |
| 660 | Amas4 = Qms(4,Iset) |
| 661 | Amas5 = Qms(5,Iset) |
| 662 | Amas6 = Qms(6,Iset) |
| 663 | |
| 664 | Xmin = Xmn (Iset) |
| 665 | Qini = Qmn (Iset) |
| 666 | Qmax = Qmx (Iset) |
| 667 | |
| 668 | Do 201 Iexp = 1, Nexp(Iset) |
| 669 | 201 ExpNor(Iexp) = ExpN(Iexp, Iset) |
| 670 | |
| 671 | Ctq2OPd = 0.0 |
| 672 | Return |
| 673 | C ************************* |
| 674 | END |
git://git.uio.no / u / mrichter / AliRoot.git / blame