| fe4da5cc |
1 | * |
| 2 | * $Id$ |
| 3 | * |
| 4 | * $Log$ |
| 5 | * Revision 1.1.1.2 1996/10/30 08:27:47 cernlib |
| 6 | * Version 7.04 |
| 7 | * |
| 8 | * Revision 1.2 1996/10/29 09:30:58 pdflib |
| 9 | * more LINNUX corrections in DATA statements |
| 10 | * |
| 11 | * Revision 1.1.1.1 1996/04/12 15:29:10 plothow |
| 12 | * Version 7.01 |
| 13 | * |
| 14 | * |
| 15 | #include "pdf/pilot.h" |
| 16 | FUNCTION Ctq2pPd (Iset, Iprtn, XX, QQ, Irt) |
| 17 | |
| 18 | C Returns xf(x,Q) -- the momentum fraction distribution !! |
| 19 | C Returns valence and sea rather than combined flavor distr. |
| 20 | |
| 21 | C Iset : PDF set label |
| 22 | |
| 23 | C Iprtn : Parton label: 2, 1 = d_ and u_ valence |
| 24 | C 0 = gluon |
| 25 | C -1, ... -6 = u, d, s, c, b, t sea quarks |
| 26 | |
| 27 | C XX : Bjorken-x |
| 28 | C QQ : scale parameter "Q" |
| 29 | C Irt : Return code |
| 30 | C 0 : no error |
| 31 | C 1 : parametrization is slightly negative; reset to 0.0. |
| 32 | C (This condition happens rarely -- only for large x where the |
| 33 | C absolute value of the parton distribution is extremely small.) |
| 34 | |
| 35 | |
| 36 | C IMPLICIT DOUBLE PRECISION (A-H, O-Z) |
| 37 | C+SEQ, IMPDP. |
| 38 | |
| 39 | PARAMETER (D0=0D0, D1=1D0, D2=2D0, D3=3D0, D4=4D0, D10=1D1) |
| 40 | PARAMETER (Nex = 5, MxFl = 6, Npn = 3, Nst = 30, Nexpt=20) |
| 41 | Parameter (Nst4 = Nst*4) |
| 42 | |
| 43 | #include "pdf/expdp.inc" |
| 44 | + P012, P34, P5, PDFD |
| 45 | |
| 46 | DIMENSION |
| 47 | > Iord(Nst), Isch(Nst), Nqrk(Nst),Alm(Nst) |
| 48 | > , Vlm(4:6,Nst), Qms(4:6, Nst) |
| 49 | > , Xmn(Nst), Qmn(Nst), Qmx(Nst), Nexp(Nexpt) |
| 50 | > , Mex(Nst), Mpn(Nst), ExpN(Nexpt, Nst), ExpNor(Nexpt) |
| 51 | |
| 52 | C RUN LE26 - CTEQ2M |
| 53 | C |
| 54 | DATA |
| 55 | > Isch(1), Iord(1), Nqrk(1), Alm(1) / 1, 2, 6, .213 / |
| 56 | > ,(Vlm(I,1), I=4,6) / .213, .139, .053 / |
| 57 | > ,(Qms(I,1), I=4,6) / 1.60, 5.00, 180.0 / |
| 58 | > ,Xmn(1), Qmn(1), Qmx(1) / 1.E-5, 1.60, 1.E3 / |
| 59 | > ,Mex(1), Mpn(1), Nexp(1) / 5, 3, 10 / |
| 60 | > ,(ExpN(I,1), I=1,10 ) |
| 61 | > / 0.990,1.012,1.022,0.980,1.062,0.870,0.843 |
| 62 | > ,0.815,0.974,1.029 / |
| 63 | |
| 64 | C RUN SA17 - CTEQ2MS |
| 65 | C |
| 66 | DATA |
| 67 | > Isch(2), Iord(2), Nqrk(2), Alm(2) / 1, 2, 6, .208 / |
| 68 | > ,(Vlm(I,2), I=4,6) / .208, .135, .051 / |
| 69 | > ,(Qms(I,2), I=4,6) / 1.60, 5.00, 180.0 / |
| 70 | > ,Xmn(2), Qmn(2), Qmx(2) / 1.E-5, 1.60, 1.E3 / |
| 71 | > ,Mex(2), Mpn(2), Nexp(2) / 5, 3, 10 / |
| 72 | > ,(ExpN(I,2), I=1,10 ) |
| 73 | > / 0.992,1.017,1.023,0.982,1.079,0.879,0.845 |
| 74 | > ,0.814,0.984,1.036 / |
| 75 | |
| 76 | C RUN FA06 - CTEQ2MF |
| 77 | C |
| 78 | DATA |
| 79 | > Isch(3), Iord(3), Nqrk(3), Alm(3) / 1, 2, 6, .208 / |
| 80 | > ,(Vlm(I,3), I=4,6) / .208, .135, .051 / |
| 81 | > ,(Qms(I,3), I=4,6) / 1.60, 5.00, 180.0 / |
| 82 | > ,Xmn(3), Qmn(3), Qmx(3) / 1.E-5, 1.60, 1.E3 / |
| 83 | > ,Mex(3), Mpn(3), Nexp(3) / 5, 3, 10 / |
| 84 | > ,(ExpN(I,3), I=1,10 ) |
| 85 | > / 0.989,1.014,1.021,0.977,1.099,0.812,0.789 |
| 86 | > ,0.822,0.909,0.950 / |
| 87 | |
| 88 | C RUN LL25 - CTEQ2ML |
| 89 | C |
| 90 | DATA |
| 91 | > Isch(4), Iord(4), Nqrk(4), Alm(4) / 1, 2, 6, .322 / |
| 92 | > ,(Vlm(I,4), I=4,6) / .322, .220, .088 / |
| 93 | > ,(Qms(I,4), I=4,6) / 1.60, 5.00, 180.0 / |
| 94 | > ,Xmn(4), Qmn(4), Qmx(4) / 1.E-5, 1.60, 1.E3 / |
| 95 | > ,Mex(4), Mpn(4), Nexp(4) / 5, 3, 10 / |
| 96 | > ,(ExpN(I,4), I=1,10 ) |
| 97 | > / 0.990,1.009,1.022,0.978,1.050,1.027,0.955 |
| 98 | > ,0.848,0.985,1.053 / |
| 99 | |
| 100 | C Run lo24 - CTEQ2L |
| 101 | C |
| 102 | DATA |
| 103 | > Isch(5), Iord(5), Nqrk(5), Alm(5) / 0, 1, 6, .190 / |
| 104 | > ,(Vlm(I,5), I=4,6) / .190, .143, .072 / |
| 105 | > ,(Qms(I,5), I=4,6) / 1.60, 5.00, 180.0 / |
| 106 | > ,Xmn(5), Qmn(5), Qmx(5) / 1.E-5, 1.60, 1.E3 / |
| 107 | > ,Mex(5), Mpn(5), Nexp(5) / 5, 3, 10 / |
| 108 | > ,(ExpN(I,5), I=1,10 ) |
| 109 | > / 0.990,1.018,1.022,0.982,0.725,0.931,0.855 |
| 110 | > ,0.795,0.965,1.007 / |
| 111 | |
| 112 | C Run da06 - CTEQ2D |
| 113 | C |
| 114 | DATA |
| 115 | > Isch(6), Iord(6), Nqrk(6), Alm(6) / 2, 2, 6, .235 / |
| 116 | > ,(Vlm(I,6), I=4,6) / .235, .155, .060 / |
| 117 | > ,(Qms(I,6), I=4,6) / 1.60, 5.00, 180.0 / |
| 118 | > ,Xmn(6), Qmn(6), Qmx(6) / 1.E-5, 1.60, 1.E3 / |
| 119 | > ,Mex(6), Mpn(6), Nexp(6) / 5, 3, 10 / |
| 120 | > ,(ExpN(I,6), I=1,10 ) |
| 121 | > / 0.989,1.019,1.021,0.978,0.961,0.958,0.904 |
| 122 | > ,0.857,0.965,1.022 / |
| 123 | |
| 124 | Data Ist, Lp, Qsto, Aln2 / 0, -10, 1.2345, 0.6931 / |
| 125 | |
| 126 | X = XX |
| 127 | QQ = MAX(QQ,Qmn(iset)) |
| 128 | Irt = 0 |
| 129 | if(Iset.eq.Ist .and. Qsto.eq.QQ) then |
| 130 | C if only change is in x: |
| 131 | if (Iprtn.eq.Lp) goto 100 |
| 132 | C if change in flv is within "light" partons: |
| 133 | if (Iprtn.ge.-3 .and. Lp.ge.-3) goto 501 |
| 134 | endif |
| 135 | |
| 136 | Ip = abs(Iprtn) |
| 137 | C Set up Qi for SB |
| 138 | If (Ip .GE. 4) then |
| 139 | If (QQ .LE. Qms(Ip, Iset)) Then |
| 140 | Ctq2pPd = 0.0 |
| 141 | Return |
| 142 | Endif |
| 143 | Qi = Qms(ip, Iset) |
| 144 | Else |
| 145 | Qi = Qmn(Iset) |
| 146 | Endif |
| 147 | C Use "standard lambda" of parametrization program |
| 148 | Alam = Alm (Iset) |
| 149 | |
| 150 | SBL = LOG(QQ/Alam) / LOG(Qi/Alam) |
| 151 | SB = LOG (SBL) |
| 152 | SB2 = SB*SB |
| 153 | SB3 = SB2*SB |
| 154 | |
| 155 | 501 Iflv = 3 - Iprtn |
| 156 | |
| 157 | Goto (1,2,3,4,5,6, 311) Iset |
| 158 | |
| 159 | 1 Goto(11,12,13,14,15,16,17,18,19)Iflv |
| 160 | C Ifl = 2 |
| 161 | 11 A0=Exp( 0.2143E+00+0.8417E+00*SB -0.2451E+01*SB2+0.9875E+00*SB3) |
| 162 | A1= 0.5209E+00-0.2384E+00*SB +0.5086E+00*SB2-0.2123E+00*SB3 |
| 163 | A2= 0.3178E+01+0.5258E+01*SB -0.8102E+01*SB2+0.3334E+01*SB3 |
| 164 | A3=-0.8537E+00+0.5921E+01*SB -0.1007E+02*SB2+0.4146E+01*SB3 |
| 165 | A4= 0.1821E+01+0.2822E-01*SB +0.1662E+00*SB2-0.1058E+00*SB3 |
| 166 | A5= 0.0000E+00-0.1090E+01*SB +0.3136E+01*SB2-0.1301E+01*SB3 |
| 167 | goto 100 |
| 168 | C Ifl = 1 |
| 169 | 12 A0=Exp(-0.1314E+01-0.1342E-01*SB +0.1136E+00*SB2-0.1557E+00*SB3) |
| 170 | A1= 0.2780E+00+0.2558E-01*SB +0.4467E-02*SB2-0.2472E-02*SB3 |
| 171 | A2= 0.3672E+01+0.5324E+00*SB +0.3531E-01*SB2+0.7928E-03*SB3 |
| 172 | A3= 0.2957E+02-0.2000E+02*SB +0.5929E+01*SB2+0.3390E+00*SB3 |
| 173 | A4= 0.8069E+00-0.2877E+00*SB +0.3574E-01*SB2+0.5622E-02*SB3 |
| 174 | A5= 0.0000E+00+0.2287E+00*SB -0.4052E-01*SB2+0.5589E-01*SB3 |
| 175 | goto 100 |
| 176 | C Ifl = 0 |
| 177 | 13 A0=Exp(-0.1059E+00-0.1461E+01*SB -0.2544E+00*SB2+0.4526E-01*SB3) |
| 178 | A1=-0.2578E+00+0.1385E+00*SB -0.1383E+00*SB2+0.3811E-01*SB3 |
| 179 | A2= 0.5195E+01+0.9648E+00*SB -0.2103E+00*SB2-0.6701E-01*SB3 |
| 180 | A3= 0.5131E+01+0.2151E+01*SB -0.2880E+01*SB2+0.6608E+00*SB3 |
| 181 | A4= 0.1118E+01+0.2636E+00*SB -0.5140E+00*SB2+0.1613E+00*SB3 |
| 182 | A5= 0.0000E+00+0.2456E+01*SB -0.8741E+00*SB2+0.2136E+00*SB3 |
| 183 | goto 100 |
| 184 | C Ifl = -1 |
| 185 | 14 A0=Exp(-0.2732E+00-0.3523E+01*SB +0.3657E+01*SB2-0.1415E+01*SB3) |
| 186 | A1=-0.3807E+00+0.1211E+00*SB -0.1231E+00*SB2+0.3753E-01*SB3 |
| 187 | A2= 0.9698E+01-0.2596E+01*SB +0.2412E+01*SB2-0.9257E+00*SB3 |
| 188 | A3=-0.6165E+00+0.1120E+01*SB -0.1708E+01*SB2+0.6383E+00*SB3 |
| 189 | A4= 0.7292E-01-0.1339E+00*SB +0.2104E+00*SB2-0.7987E-01*SB3 |
| 190 | A5=-0.1370E+01+0.2452E+01*SB -0.1804E+01*SB2+0.6459E+00*SB3 |
| 191 | goto 100 |
| 192 | C Ifl = -2 |
| 193 | 15 A0=Exp(-0.2319E+01-0.3182E+01*SB +0.3572E+01*SB2-0.1431E+01*SB3) |
| 194 | A1=-0.2622E+00+0.3085E+00*SB -0.4394E+00*SB2+0.1496E+00*SB3 |
| 195 | A2= 0.9481E+01-0.3627E+01*SB +0.5640E+01*SB2-0.2265E+01*SB3 |
| 196 | A3= 0.5000E+02-0.1851E+02*SB +0.2640E+01*SB2-0.6001E+00*SB3 |
| 197 | A4= 0.1566E+01-0.7375E+00*SB +0.8736E+00*SB2-0.3449E+00*SB3 |
| 198 | A5=-0.7983E-01+0.3236E+01*SB -0.3373E+01*SB2+0.1236E+01*SB3 |
| 199 | goto 100 |
| 200 | C Ifl = -3 |
| 201 | 16 A0=Exp(-0.1855E+01-0.5302E+01*SB +0.8433E+00*SB2-0.1236E+00*SB3) |
| 202 | A1=-0.4000E-02-0.1345E+01*SB +0.1192E+01*SB2-0.3039E+00*SB3 |
| 203 | A2= 0.6870E+01+0.1246E+01*SB -0.8968E+00*SB2-0.9791E-01*SB3 |
| 204 | A3= 0.0000E+00+0.4616E+01*SB +0.1026E+02*SB2+0.2844E+02*SB3 |
| 205 | A4= 0.1000E-02+0.4098E+00*SB -0.4250E+00*SB2+0.1100E+00*SB3 |
| 206 | A5= 0.0000E+00-0.2151E+01*SB +0.2991E+01*SB2-0.7717E+00*SB3 |
| 207 | goto 100 |
| 208 | C Ifl = -4 |
| 209 | 17 A0=SB** 0.7722E+00*Exp(-0.7241E+01-0.7885E-01*SB -0.1124E+01*SB2) |
| 210 | A1=-0.3971E+00+0.9132E+00*SB -0.1175E+01*SB2+0.3573E+00*SB3 |
| 211 | A2= 0.6367E+01-0.6565E+01*SB +0.8114E+01*SB2-0.2666E+01*SB3 |
| 212 | A3= 0.2878E+02-0.2000E+02*SB +0.7000E+00*SB2+0.3000E+02*SB3 |
| 213 | A4= 0.1010E+00-0.4592E+00*SB +0.5877E+00*SB2-0.1472E+00*SB3 |
| 214 | A5= 0.1749E+00+0.3875E+01*SB -0.3768E+01*SB2+0.1316E+01*SB3 |
| 215 | goto 100 |
| 216 | C Ifl = -5 |
| 217 | 18 A0=SB** 0.1299E+00*Exp(-0.4868E+01-0.4339E+01*SB +0.7080E+00*SB2) |
| 218 | A1=-0.1705E+00-0.3381E+00*SB +0.5287E+00*SB2-0.2644E+00*SB3 |
| 219 | A2= 0.5610E+01-0.1365E+01*SB +0.1835E+01*SB2-0.5655E+00*SB3 |
| 220 | A3=-0.1001E+01+0.3044E+01*SB +0.2680E+01*SB2+0.1426E+02*SB3 |
| 221 | A4= 0.3814E-02+0.3430E+00*SB -0.6926E+00*SB2+0.3486E+00*SB3 |
| 222 | A5= 0.1156E+01+0.2016E+01*SB -0.1674E+01*SB2+0.5981E+00*SB3 |
| 223 | goto 100 |
| 224 | C Ifl = -6 |
| 225 | 19 A0=SB** 0.9819E+00*Exp(-0.7859E+01+0.6819E+00*SB -0.3386E+01*SB2) |
| 226 | A1=-0.1055E+00-0.1413E+01*SB +0.3451E+01*SB2-0.2466E+01*SB3 |
| 227 | A2= 0.4055E+01+0.8107E+01*SB -0.1576E+02*SB2+0.8094E+01*SB3 |
| 228 | A3= 0.3799E+01+0.9616E+01*SB -0.1984E+02*SB2+0.2641E+02*SB3 |
| 229 | A4= 0.3619E+00-0.8627E+00*SB -0.9390E-01*SB2+0.9196E+00*SB3 |
| 230 | A5= 0.3779E+01-0.6073E+01*SB +0.9999E+01*SB2-0.4304E+01*SB3 |
| 231 | goto 100 |
| 232 | |
| 233 | 2 Goto(21,22,23,24,25,26,27,28,29)Iflv |
| 234 | C Ifl = 2 |
| 235 | 21 A0=Exp( 0.2790E+00+0.7294E+00*SB -0.2202E+01*SB2+0.8599E+00*SB3) |
| 236 | A1= 0.5380E+00-0.2261E+00*SB +0.4636E+00*SB2-0.1871E+00*SB3 |
| 237 | A2= 0.3259E+01+0.2141E+01*SB -0.2947E+01*SB2+0.1245E+01*SB3 |
| 238 | A3=-0.8390E+00+0.1448E+01*SB -0.2331E+01*SB2+0.8658E+00*SB3 |
| 239 | A4= 0.1847E+01-0.3943E+01*SB +0.5998E+01*SB2-0.2191E+01*SB3 |
| 240 | A5= 0.0000E+00-0.9719E+00*SB +0.2830E+01*SB2-0.1137E+01*SB3 |
| 241 | goto 100 |
| 242 | C Ifl = 1 |
| 243 | 22 A0=Exp(-0.1318E+01+0.2328E-01*SB +0.5179E-01*SB2-0.1305E+00*SB3) |
| 244 | A1= 0.2760E+00+0.4429E-01*SB -0.2626E-01*SB2+0.7143E-02*SB3 |
| 245 | A2= 0.3660E+01+0.5232E+00*SB +0.5491E-01*SB2-0.4115E-02*SB3 |
| 246 | A3= 0.2910E+02-0.2000E+02*SB +0.6631E+01*SB2-0.3050E-01*SB3 |
| 247 | A4= 0.8010E+00-0.2688E+00*SB +0.1051E-01*SB2+0.1195E-01*SB3 |
| 248 | A5= 0.0000E+00+0.2887E+00*SB -0.1398E+00*SB2+0.8194E-01*SB3 |
| 249 | goto 100 |
| 250 | C Ifl = 0 |
| 251 | 23 A0=Exp(-0.1623E+01-0.7232E+00*SB +0.1889E+00*SB2+0.1140E+00*SB3) |
| 252 | A1=-0.5000E+00+0.8611E-01*SB +0.2203E-01*SB2-0.1401E-01*SB3 |
| 253 | A2= 0.3821E+01+0.8976E+00*SB +0.1400E+00*SB2-0.9163E-01*SB3 |
| 254 | A3= 0.5809E+01-0.5060E+01*SB +0.3808E+00*SB2+0.2519E+00*SB3 |
| 255 | A4= 0.4500E+00-0.5121E+00*SB +0.1979E+00*SB2-0.2705E-01*SB3 |
| 256 | A5= 0.0000E+00+0.1210E+01*SB -0.2921E+00*SB2+0.1240E+00*SB3 |
| 257 | goto 100 |
| 258 | C Ifl = -1 |
| 259 | 24 A0=Exp(-0.6986E-01-0.5954E+00*SB -0.1582E+01*SB2+0.5104E+00*SB3) |
| 260 | A1=-0.8461E+00+0.2127E+00*SB +0.9425E-01*SB2-0.5264E-01*SB3 |
| 261 | A2= 0.1200E+02+0.1659E+01*SB -0.5354E+01*SB2+0.1795E+01*SB3 |
| 262 | A3= 0.2958E+02+0.3000E+02*SB +0.3000E+02*SB2-0.1965E+02*SB3 |
| 263 | A4= 0.4000E+01-0.4865E+00*SB +0.9460E+00*SB2+0.3432E+00*SB3 |
| 264 | A5=-0.3378E+01+0.1656E+01*SB +0.1123E+01*SB2-0.4667E+00*SB3 |
| 265 | goto 100 |
| 266 | C Ifl = -2 |
| 267 | 25 A0=Exp(-0.1929E+01-0.2626E+01*SB +0.2926E+01*SB2-0.1297E+01*SB3) |
| 268 | A1=-0.6627E+00+0.4561E+00*SB -0.3818E+00*SB2+0.1239E+00*SB3 |
| 269 | A2= 0.9506E+01-0.2724E+01*SB +0.4283E+01*SB2-0.1804E+01*SB3 |
| 270 | A3= 0.1897E+02+0.1642E+01*SB -0.8390E+01*SB2+0.3894E+01*SB3 |
| 271 | A4= 0.1024E+01-0.1786E+00*SB +0.4535E+00*SB2-0.2075E+00*SB3 |
| 272 | A5=-0.1746E+01+0.3572E+01*SB -0.2908E+01*SB2+0.1093E+01*SB3 |
| 273 | goto 100 |
| 274 | C Ifl = -3 |
| 275 | 26 A0=Exp(-0.4913E+00-0.6866E+01*SB +0.1432E+01*SB2-0.1749E+00*SB3) |
| 276 | A1=-0.1157E+00-0.1567E+01*SB +0.1439E+01*SB2-0.3724E+00*SB3 |
| 277 | A2= 0.7730E+01+0.9748E+00*SB -0.1157E+01*SB2-0.8358E-02*SB3 |
| 278 | A3=-0.6050E+00+0.1835E+01*SB +0.3788E+01*SB2+0.3000E+02*SB3 |
| 279 | A4= 0.1620E-08+0.4590E+00*SB -0.4070E+00*SB2+0.8900E-01*SB3 |
| 280 | A5=-0.7048E+00-0.2505E+01*SB +0.4000E+01*SB2-0.1161E+01*SB3 |
| 281 | goto 100 |
| 282 | C Ifl = -4 |
| 283 | 27 A0=SB** 0.7393E+00*Exp(-0.6518E+01-0.3998E+00*SB -0.1111E+01*SB2) |
| 284 | A1=-0.6482E+00+0.1125E+01*SB -0.1290E+01*SB2+0.3940E+00*SB3 |
| 285 | A2= 0.8487E+01-0.9235E+01*SB +0.9353E+01*SB2-0.2913E+01*SB3 |
| 286 | A3= 0.2265E+02-0.1999E+02*SB +0.4105E+01*SB2+0.2144E+02*SB3 |
| 287 | A4= 0.8990E-01-0.4372E+00*SB +0.5941E+00*SB2-0.1469E+00*SB3 |
| 288 | A5=-0.9690E+00+0.5068E+01*SB -0.4368E+01*SB2+0.1503E+01*SB3 |
| 289 | goto 100 |
| 290 | C Ifl = -5 |
| 291 | 28 A0=SB** 0.9880E+00*Exp(-0.7180E+01-0.2494E+01*SB +0.3561E-01*SB2) |
| 292 | A1=-0.4301E+00-0.2611E+00*SB +0.3914E+00*SB2-0.1638E+00*SB3 |
| 293 | A2= 0.5137E+01+0.1506E+01*SB -0.9588E+00*SB2-0.1596E+00*SB3 |
| 294 | A3= 0.1483E+02+0.2998E+02*SB +0.2357E+02*SB2-0.9353E+01*SB3 |
| 295 | A4= 0.2426E+00+0.1371E+00*SB -0.3791E+00*SB2+0.1948E+00*SB3 |
| 296 | A5= 0.1463E+01+0.1907E+00*SB +0.3557E+00*SB2+0.2097E-01*SB3 |
| 297 | goto 100 |
| 298 | C Ifl = -6 |
| 299 | 29 A0=SB** 0.1005E+01*Exp(-0.5255E+01-0.9866E-01*SB -0.2737E+01*SB2) |
| 300 | A1=-0.3140E+00-0.2055E+00*SB +0.5594E+00*SB2-0.2960E+00*SB3 |
| 301 | A2= 0.9227E+01-0.4569E+01*SB -0.9724E+01*SB2+0.1026E+02*SB3 |
| 302 | A3= 0.1131E+02-0.1972E+02*SB -0.1107E+02*SB2+0.2311E+02*SB3 |
| 303 | A4= 0.1488E+01+0.1737E+01*SB +0.4323E+01*SB2-0.9925E+01*SB3 |
| 304 | A5= 0.1895E+01-0.7350E+00*SB +0.3780E+01*SB2-0.1408E+01*SB3 |
| 305 | goto 100 |
| 306 | |
| 307 | 3 Goto(31,32,33,34,35,36,37,38,39)Iflv |
| 308 | C Ifl = 2 |
| 309 | 31 A0=Exp(-0.7913E+00-0.2789E+01*SB -0.7289E-01*SB2+0.1770E+00*SB3) |
| 310 | A1= 0.4942E+00-0.7886E-01*SB +0.9057E-01*SB2-0.5259E-01*SB3 |
| 311 | A2= 0.3727E+01+0.1089E+01*SB -0.1004E+01*SB2+0.4345E+00*SB3 |
| 312 | A3= 0.1944E+01+0.7846E+01*SB +0.7984E+01*SB2+0.5548E+01*SB3 |
| 313 | A4= 0.2940E-02+0.8428E-04*SB +0.1266E+00*SB2-0.3517E-01*SB3 |
| 314 | A5=-0.1060E+00-0.1192E-01*SB +0.1130E+01*SB2-0.4527E+00*SB3 |
| 315 | goto 100 |
| 316 | C Ifl = 1 |
| 317 | 32 A0=Exp(-0.1344E+01+0.7859E-02*SB +0.4623E-01*SB2-0.1273E+00*SB3) |
| 318 | A1= 0.2760E+00+0.4201E-01*SB -0.1795E-01*SB2+0.3212E-02*SB3 |
| 319 | A2= 0.3660E+01+0.5247E+00*SB +0.4405E-01*SB2+0.1391E-02*SB3 |
| 320 | A3= 0.2981E+02-0.2000E+02*SB +0.6566E+01*SB2+0.2479E-01*SB3 |
| 321 | A4= 0.7950E+00-0.2732E+00*SB +0.2470E-01*SB2+0.6157E-02*SB3 |
| 322 | A5= 0.0000E+00+0.2793E+00*SB -0.9197E-01*SB2+0.5953E-01*SB3 |
| 323 | goto 100 |
| 324 | C Ifl = 0 |
| 325 | 33 A0=Exp( 0.9746E+00-0.3252E+01*SB +0.1664E+01*SB2-0.6410E+00*SB3) |
| 326 | A1=-0.5271E-02-0.3198E+00*SB +0.1279E+00*SB2-0.1256E-02*SB3 |
| 327 | A2= 0.5740E+01-0.3139E+01*SB +0.3841E+01*SB2-0.1415E+01*SB3 |
| 328 | A3= 0.7161E-01-0.4363E+01*SB +0.4925E+01*SB2-0.1614E+01*SB3 |
| 329 | A4= 0.1860E+01+0.1342E+01*SB -0.2234E+01*SB2+0.1047E+01*SB3 |
| 330 | A5= 0.7409E-01+0.2390E+01*SB -0.1457E+01*SB2+0.5853E+00*SB3 |
| 331 | goto 100 |
| 332 | C Ifl = -1 |
| 333 | 34 A0=Exp(-0.8454E+00-0.3334E+01*SB +0.3591E+01*SB2-0.1485E+01*SB3) |
| 334 | A1=-0.2826E-02-0.2810E+00*SB -0.3809E-01*SB2+0.6585E-01*SB3 |
| 335 | A2= 0.9139E+01-0.2811E+01*SB +0.4730E+01*SB2-0.2157E+01*SB3 |
| 336 | A3=-0.3120E+00+0.1217E+01*SB -0.1726E+01*SB2+0.6220E+00*SB3 |
| 337 | A4= 0.1793E-01-0.4608E-01*SB +0.5294E-01*SB2-0.1709E-01*SB3 |
| 338 | A5=-0.1471E+00+0.1104E+01*SB -0.1358E+01*SB2+0.7200E+00*SB3 |
| 339 | goto 100 |
| 340 | C Ifl = -2 |
| 341 | 35 A0=Exp(-0.1398E+01-0.3536E+01*SB +0.3849E+01*SB2-0.1549E+01*SB3) |
| 342 | A1=-0.1332E-01-0.2155E-01*SB -0.3404E+00*SB2+0.1569E+00*SB3 |
| 343 | A2= 0.9981E+01-0.3499E+01*SB +0.5448E+01*SB2-0.2198E+01*SB3 |
| 344 | A3= 0.3736E+02-0.2000E+02*SB +0.6675E+01*SB2-0.7276E+00*SB3 |
| 345 | A4= 0.1705E+01-0.1013E+01*SB +0.1122E+01*SB2-0.4057E+00*SB3 |
| 346 | A5=-0.1189E-01+0.2698E+01*SB -0.3429E+01*SB2+0.1389E+01*SB3 |
| 347 | goto 100 |
| 348 | C Ifl = -3 |
| 349 | 36 A0=Exp(-0.2979E+01-0.6085E+01*SB +0.2428E+01*SB2-0.6482E+00*SB3) |
| 350 | A1=-0.1372E+00-0.1281E+00*SB +0.1587E+00*SB2-0.9637E-01*SB3 |
| 351 | A2= 0.7009E+01-0.1609E+01*SB +0.2765E+01*SB2-0.1177E+01*SB3 |
| 352 | A3= 0.1308E+01+0.9583E+01*SB +0.2360E+02*SB2+0.2999E+02*SB3 |
| 353 | A4= 0.2509E-01+0.2106E+00*SB -0.4405E+00*SB2+0.2075E+00*SB3 |
| 354 | A5=-0.2069E-01+0.1971E+01*SB -0.1615E+01*SB2+0.6039E+00*SB3 |
| 355 | goto 100 |
| 356 | C Ifl = -4 |
| 357 | 37 A0=SB** 0.8072E+00*Exp(-0.6920E+01-0.5031E+00*SB -0.9965E+00*SB2) |
| 358 | A1=-0.2118E+00+0.7930E+00*SB -0.1101E+01*SB2+0.3302E+00*SB3 |
| 359 | A2= 0.8039E+01-0.7170E+01*SB +0.8657E+01*SB2-0.2893E+01*SB3 |
| 360 | A3= 0.2926E+02-0.1993E+02*SB +0.1841E+01*SB2+0.2996E+02*SB3 |
| 361 | A4= 0.1339E+00-0.5531E+00*SB +0.6505E+00*SB2-0.1595E+00*SB3 |
| 362 | A5= 0.7439E+00+0.3307E+01*SB -0.3284E+01*SB2+0.1152E+01*SB3 |
| 363 | goto 100 |
| 364 | C Ifl = -5 |
| 365 | 38 A0=SB** 0.9925E+00*Exp(-0.2190E+01-0.3393E+01*SB -0.8631E+00*SB2) |
| 366 | A1=-0.1261E+00-0.2368E+00*SB +0.4143E+00*SB2-0.1577E+00*SB3 |
| 367 | A2= 0.4585E+01+0.5227E+01*SB -0.3248E+01*SB2-0.2599E+00*SB3 |
| 368 | A3=-0.1094E+01+0.4927E+00*SB -0.9921E+00*SB2+0.3138E+01*SB3 |
| 369 | A4= 0.1396E+00+0.2562E+00*SB +0.1844E+00*SB2-0.1599E+00*SB3 |
| 370 | A5= 0.8621E+00+0.4715E+00*SB +0.2547E+01*SB2-0.8429E+00*SB3 |
| 371 | goto 100 |
| 372 | C Ifl = -6 |
| 373 | 39 A0=SB** 0.1016E+01*Exp(-0.5397E+01-0.1979E+01*SB -0.2441E+00*SB2) |
| 374 | A1=-0.1426E+00-0.2861E+00*SB +0.7434E+00*SB2-0.5214E+00*SB3 |
| 375 | A2= 0.6363E+01+0.4028E+00*SB -0.8356E+01*SB2+0.6814E+01*SB3 |
| 376 | A3=-0.2526E+00+0.2425E+01*SB -0.1407E+02*SB2+0.3000E+02*SB3 |
| 377 | A4= 0.1125E+00-0.1089E+01*SB +0.9977E+01*SB2+0.1000E+02*SB3 |
| 378 | A5= 0.2669E+01-0.6366E+00*SB +0.4355E+01*SB2-0.2919E+01*SB3 |
| 379 | goto 100 |
| 380 | |
| 381 | 4 Goto(41,42,43,44,45,46,47,48,49)Iflv |
| 382 | C Ifl = 2 |
| 383 | 41 A0=Exp( 0.3760E+00+0.5491E+00*SB -0.1845E+01*SB2+0.6803E+00*SB3) |
| 384 | A1= 0.5650E+00-0.1953E+00*SB +0.3761E+00*SB2-0.1419E+00*SB3 |
| 385 | A2= 0.3464E+01+0.3817E+01*SB -0.5384E+01*SB2+0.2057E+01*SB3 |
| 386 | A3=-0.5850E+00+0.5566E+01*SB -0.9000E+01*SB2+0.3433E+01*SB3 |
| 387 | A4= 0.2322E+01-0.1431E+00*SB +0.3901E+00*SB2-0.1678E+00*SB3 |
| 388 | A5= 0.0000E+00-0.7370E+00*SB +0.2310E+01*SB2-0.8743E+00*SB3 |
| 389 | goto 100 |
| 390 | C Ifl = 1 |
| 391 | 42 A0=Exp(-0.1324E+01+0.1169E-01*SB +0.1969E-01*SB2-0.7583E-01*SB3) |
| 392 | A1= 0.2890E+00+0.5832E-01*SB -0.2921E-01*SB2+0.4701E-02*SB3 |
| 393 | A2= 0.3580E+01+0.5291E+00*SB -0.5662E-02*SB2+0.2746E-01*SB3 |
| 394 | A3= 0.3021E+02-0.1999E+02*SB +0.6250E+01*SB2-0.3035E+00*SB3 |
| 395 | A4= 0.7990E+00-0.2531E+00*SB +0.5556E-02*SB2+0.8272E-02*SB3 |
| 396 | A5= 0.0000E+00+0.3674E+00*SB -0.1383E+00*SB2+0.4665E-01*SB3 |
| 397 | goto 100 |
| 398 | C Ifl = 0 |
| 399 | 43 A0=Exp(-0.1920E+00-0.7015E+00*SB -0.9113E+00*SB2+0.2352E+00*SB3) |
| 400 | A1=-0.2120E+00+0.1133E-01*SB -0.1553E-01*SB2+0.2822E-02*SB3 |
| 401 | A2= 0.4549E+01+0.1250E+01*SB -0.4647E+00*SB2+0.9617E-01*SB3 |
| 402 | A3= 0.1197E+02-0.4156E+01*SB +0.1413E+00*SB2+0.1607E+00*SB3 |
| 403 | A4= 0.1616E+01+0.1082E+00*SB -0.6651E+00*SB2+0.2356E+00*SB3 |
| 404 | A5= 0.0000E+00+0.1824E+01*SB -0.2063E+00*SB2+0.1148E-01*SB3 |
| 405 | goto 100 |
| 406 | C Ifl = -1 |
| 407 | 44 A0=Exp(-0.1388E+01-0.7408E+00*SB -0.6454E+00*SB2+0.2373E+00*SB3) |
| 408 | A1=-0.2928E+00-0.1726E-01*SB +0.4033E-01*SB2-0.2514E-01*SB3 |
| 409 | A2= 0.9975E+01-0.2048E+01*SB -0.6060E+00*SB2+0.5225E+00*SB3 |
| 410 | A3= 0.2687E+02-0.4683E+01*SB -0.1999E+02*SB2+0.1188E+02*SB3 |
| 411 | A4= 0.4000E+01-0.6773E+00*SB +0.4301E+00*SB2+0.4524E+00*SB3 |
| 412 | A5=-0.7164E+00+0.7488E+00*SB +0.5766E+00*SB2-0.2609E+00*SB3 |
| 413 | goto 100 |
| 414 | C Ifl = -2 |
| 415 | 45 A0=Exp(-0.2272E+01-0.2998E+01*SB +0.3282E+01*SB2-0.1203E+01*SB3) |
| 416 | A1=-0.2062E+00+0.3320E+00*SB -0.5074E+00*SB2+0.1655E+00*SB3 |
| 417 | A2= 0.9667E+01-0.3497E+01*SB +0.5271E+01*SB2-0.1984E+01*SB3 |
| 418 | A3= 0.4996E+02-0.3241E+01*SB -0.1425E+02*SB2+0.3849E+01*SB3 |
| 419 | A4= 0.1619E+01-0.5354E+00*SB +0.5753E+00*SB2-0.2238E+00*SB3 |
| 420 | A5= 0.8755E-01+0.3195E+01*SB -0.3496E+01*SB2+0.1197E+01*SB3 |
| 421 | goto 100 |
| 422 | C Ifl = -3 |
| 423 | 46 A0=Exp(-0.1864E+01-0.5258E+01*SB +0.1034E+01*SB2-0.1550E+00*SB3) |
| 424 | A1= 0.1000E-02-0.1090E+01*SB +0.8345E+00*SB2-0.1887E+00*SB3 |
| 425 | A2= 0.6898E+01-0.4951E+00*SB +0.4279E+00*SB2-0.2727E+00*SB3 |
| 426 | A3= 0.0000E+00+0.4322E+01*SB +0.8181E+01*SB2+0.2309E+02*SB3 |
| 427 | A4= 0.1000E-02+0.3550E+00*SB -0.3220E+00*SB2+0.7294E-01*SB3 |
| 428 | A5= 0.0000E+00-0.1347E+01*SB +0.1896E+01*SB2-0.4491E+00*SB3 |
| 429 | goto 100 |
| 430 | C Ifl = -4 |
| 431 | 47 A0=SB** 0.7528E+00*Exp(-0.7684E+01+0.6791E-01*SB -0.9094E+00*SB2) |
| 432 | A1=-0.3732E+00+0.8408E+00*SB -0.1020E+01*SB2+0.3046E+00*SB3 |
| 433 | A2= 0.4984E+01-0.5534E+01*SB +0.6418E+01*SB2-0.1856E+01*SB3 |
| 434 | A3= 0.3761E+02-0.1999E+02*SB -0.3358E+01*SB2+0.2999E+02*SB3 |
| 435 | A4= 0.1161E+00-0.4680E+00*SB +0.5567E+00*SB2-0.1633E+00*SB3 |
| 436 | A5= 0.3028E+00+0.3339E+01*SB -0.3004E+01*SB2+0.9160E+00*SB3 |
| 437 | goto 100 |
| 438 | C Ifl = -5 |
| 439 | 48 A0=SB** 0.1011E+01*Exp(-0.7217E+01-0.2288E+01*SB +0.3450E+00*SB2) |
| 440 | A1=-0.1955E+00-0.3371E+00*SB +0.5111E+00*SB2-0.2210E+00*SB3 |
| 441 | A2= 0.4302E+01-0.1214E+01*SB +0.3104E+01*SB2-0.1408E+01*SB3 |
| 442 | A3= 0.1487E+02+0.1549E+02*SB +0.2875E+02*SB2-0.1922E+02*SB3 |
| 443 | A4= 0.8935E-02+0.3571E+00*SB -0.6668E+00*SB2+0.3037E+00*SB3 |
| 444 | A5= 0.1570E+01+0.7105E+00*SB -0.6070E+00*SB2+0.3796E+00*SB3 |
| 445 | goto 100 |
| 446 | C Ifl = -6 |
| 447 | 49 A0=SB** 0.9986E+00*Exp(-0.5847E+01-0.2798E+00*SB -0.9882E+00*SB2) |
| 448 | A1=-0.2154E+00-0.8282E-01*SB +0.3611E-01*SB2+0.2623E-01*SB3 |
| 449 | A2= 0.3250E+01+0.9635E+01*SB -0.1274E+02*SB2+0.4453E+01*SB3 |
| 450 | A3=-0.2594E+01+0.9097E+01*SB +0.1581E+02*SB2-0.9123E+01*SB3 |
| 451 | A4= 0.1768E+01-0.2749E+01*SB +0.9999E+01*SB2+0.9995E+01*SB3 |
| 452 | A5= 0.2521E+01-0.1802E-01*SB +0.4820E+00*SB2+0.2004E+00*SB3 |
| 453 | goto 100 |
| 454 | |
| 455 | 5 Goto(51,52,53,54,55,56,57,58,59)Iflv |
| 456 | C Ifl = 2 |
| 457 | 51 A0=Exp( 0.7248E-01+0.3941E+00*SB -0.1772E+01*SB2+0.7629E+00*SB3) |
| 458 | A1= 0.4964E+00-0.1224E+00*SB +0.3646E+00*SB2-0.1685E+00*SB3 |
| 459 | A2= 0.3000E+01+0.2780E+01*SB -0.4028E+01*SB2+0.1816E+01*SB3 |
| 460 | A3=-0.1064E+01+0.3062E+01*SB -0.5927E+01*SB2+0.2785E+01*SB3 |
| 461 | A4= 0.3193E+01+0.1499E+01*SB -0.2765E+01*SB2+0.1019E+01*SB3 |
| 462 | A5= 0.1524E-01-0.4541E+00*SB +0.2281E+01*SB2-0.1033E+01*SB3 |
| 463 | goto 100 |
| 464 | C Ifl = 1 |
| 465 | 52 A0=Exp(-0.1794E+01-0.2055E+00*SB -0.3350E-01*SB2-0.5084E-01*SB3) |
| 466 | A1= 0.1748E+00+0.4637E-01*SB -0.2048E-01*SB2+0.2596E-02*SB3 |
| 467 | A2= 0.3321E+01+0.6253E+00*SB +0.2148E-01*SB2+0.1288E-01*SB3 |
| 468 | A3= 0.4355E+02-0.2000E+02*SB +0.5486E+01*SB2+0.1536E+00*SB3 |
| 469 | A4= 0.9586E+00-0.3217E+00*SB +0.4458E-01*SB2-0.1404E-03*SB3 |
| 470 | A5=-0.6595E-02+0.3499E+00*SB -0.7048E-01*SB2+0.2619E-01*SB3 |
| 471 | goto 100 |
| 472 | C Ifl = 0 |
| 473 | 53 A0=Exp(-0.6194E+00-0.2643E+00*SB -0.1875E+01*SB2+0.6011E+00*SB3) |
| 474 | A1=-0.2600E+00+0.8704E-01*SB -0.7375E-01*SB2+0.1876E-01*SB3 |
| 475 | A2= 0.4620E+01+0.1578E+01*SB -0.8411E+00*SB2+0.1527E+00*SB3 |
| 476 | A3= 0.1604E+02-0.1230E+02*SB +0.6939E+01*SB2-0.2012E+01*SB3 |
| 477 | A4= 0.1255E+01+0.4769E+00*SB -0.9915E+00*SB2+0.3439E+00*SB3 |
| 478 | A5= 0.1116E-02+0.2409E+01*SB -0.4442E+00*SB2+0.3431E-01*SB3 |
| 479 | goto 100 |
| 480 | C Ifl = -1 |
| 481 | 54 A0=Exp(-0.1571E+01-0.1905E+00*SB -0.8672E+00*SB2+0.2070E+00*SB3) |
| 482 | A1=-0.3266E+00+0.6428E-01*SB -0.8694E-01*SB2+0.1778E-01*SB3 |
| 483 | A2= 0.8921E+01-0.5010E+00*SB -0.9658E+00*SB2+0.3893E+00*SB3 |
| 484 | A3= 0.1329E+02+0.4652E+01*SB -0.2000E+02*SB2+0.1001E+02*SB3 |
| 485 | A4= 0.3283E+01-0.3400E+00*SB -0.1957E+00*SB2+0.8063E+00*SB3 |
| 486 | A5=-0.5701E+00+0.4042E+00*SB +0.5239E+00*SB2-0.1665E+00*SB3 |
| 487 | goto 100 |
| 488 | C Ifl = -2 |
| 489 | 55 A0=Exp(-0.2281E+01-0.2768E+01*SB +0.3137E+01*SB2-0.1278E+01*SB3) |
| 490 | A1=-0.2624E+00+0.4142E+00*SB -0.5936E+00*SB2+0.1937E+00*SB3 |
| 491 | A2= 0.9438E+01-0.3179E+01*SB +0.5107E+01*SB2-0.2179E+01*SB3 |
| 492 | A3= 0.5000E+02-0.1802E+02*SB -0.7515E+01*SB2+0.2991E+01*SB3 |
| 493 | A4= 0.1809E+01-0.9121E+00*SB +0.8854E+00*SB2-0.3582E+00*SB3 |
| 494 | A5= 0.4056E-01+0.3033E+01*SB -0.3431E+01*SB2+0.1253E+01*SB3 |
| 495 | goto 100 |
| 496 | C Ifl = -3 |
| 497 | 56 A0=Exp(-0.2318E+01-0.4104E+01*SB -0.1502E+00*SB2+0.1693E+00*SB3) |
| 498 | A1=-0.2251E-01-0.1101E+01*SB +0.1037E+01*SB2-0.3290E+00*SB3 |
| 499 | A2= 0.6989E+01+0.1794E+01*SB -0.1811E+01*SB2+0.3061E+00*SB3 |
| 500 | A3= 0.7972E+00+0.7806E+01*SB +0.1869E+02*SB2+0.2999E+02*SB3 |
| 501 | A4= 0.4795E-01+0.1622E+00*SB -0.3977E+00*SB2+0.1920E+00*SB3 |
| 502 | A5=-0.5275E-01-0.2616E+01*SB +0.3076E+01*SB2-0.7425E+00*SB3 |
| 503 | goto 100 |
| 504 | C Ifl = -4 |
| 505 | 57 A0=SB** 0.8431E+00*Exp(-0.6539E+01-0.1875E+00*SB -0.1346E+01*SB2) |
| 506 | A1=-0.4970E+00+0.9062E+00*SB -0.1169E+01*SB2+0.3703E+00*SB3 |
| 507 | A2= 0.4939E+01-0.2995E+01*SB +0.4483E+01*SB2-0.1704E+01*SB3 |
| 508 | A3= 0.3113E+02-0.1997E+02*SB +0.1540E+01*SB2+0.3000E+02*SB3 |
| 509 | A4= 0.1349E+00-0.5418E+00*SB +0.6142E+00*SB2-0.1360E+00*SB3 |
| 510 | A5=-0.8590E+00+0.3956E+01*SB -0.3612E+01*SB2+0.1401E+01*SB3 |
| 511 | goto 100 |
| 512 | C Ifl = -5 |
| 513 | 58 A0=SB** 0.2639E-01*Exp(-0.2099E+01-0.2681E+01*SB +0.2925E+00*SB2) |
| 514 | A1=-0.2243E+00-0.5343E-01*SB -0.1953E-01*SB2+0.1586E-01*SB3 |
| 515 | A2= 0.4294E+01+0.1102E+01*SB -0.1822E+00*SB2-0.2481E+00*SB3 |
| 516 | A3=-0.9998E+00+0.8275E-01*SB +0.5494E+00*SB2-0.1982E+00*SB3 |
| 517 | A4= 0.5904E-04+0.9222E-01*SB -0.9293E-01*SB2+0.9159E-01*SB3 |
| 518 | A5= 0.2657E+00+0.1770E+01*SB -0.7111E+00*SB2+0.2525E+00*SB3 |
| 519 | goto 100 |
| 520 | C Ifl = -6 |
| 521 | 59 A0=SB** 0.1009E+01*Exp(-0.7032E+01+0.4562E+01*SB -0.9081E+01*SB2) |
| 522 | A1=-0.1412E+00-0.5076E+00*SB +0.9513E+00*SB2-0.4326E+00*SB3 |
| 523 | A2= 0.5385E+01+0.3023E+01*SB -0.1162E+02*SB2+0.7006E+01*SB3 |
| 524 | A3= 0.4997E+01-0.1600E+02*SB +0.1342E+02*SB2+0.1197E+02*SB3 |
| 525 | A4= 0.5825E+00+0.3994E+00*SB -0.1255E+01*SB2+0.6486E+00*SB3 |
| 526 | A5= 0.3365E+01-0.4026E+01*SB +0.8385E+01*SB2-0.2260E+01*SB3 |
| 527 | goto 100 |
| 528 | |
| 529 | 6 Goto(61,62,63,64,65,66,67,68,69)Iflv |
| 530 | C Ifl = 2 |
| 531 | 61 A0=Exp( 0.1590E+00+0.5580E+00*SB -0.1838E+01*SB2+0.7018E+00*SB3) |
| 532 | A1= 0.5110E+00-0.1625E+00*SB +0.3547E+00*SB2-0.1412E+00*SB3 |
| 533 | A2= 0.3158E+01+0.3962E+01*SB -0.5866E+01*SB2+0.2375E+01*SB3 |
| 534 | A3=-0.6000E+00+0.6144E+01*SB -0.1056E+02*SB2+0.4345E+01*SB3 |
| 535 | A4= 0.2306E+01-0.4669E-01*SB +0.2711E+00*SB2-0.1640E+00*SB3 |
| 536 | A5= 0.0000E+00-0.6638E+00*SB +0.2239E+01*SB2-0.8843E+00*SB3 |
| 537 | goto 100 |
| 538 | C Ifl = 1 |
| 539 | 62 A0=Exp(-0.1182E+01+0.1449E+00*SB +0.2753E-01*SB2-0.1009E+00*SB3) |
| 540 | A1= 0.2540E+00+0.2686E-01*SB -0.1546E-01*SB2+0.5396E-02*SB3 |
| 541 | A2= 0.3442E+01+0.5576E+00*SB +0.1937E-01*SB2+0.6696E-02*SB3 |
| 542 | A3= 0.2545E+02-0.2000E+02*SB +0.7355E+01*SB2-0.7058E+00*SB3 |
| 543 | A4= 0.9170E+00-0.3090E+00*SB +0.1705E-01*SB2+0.8534E-02*SB3 |
| 544 | A5= 0.0000E+00+0.1449E+00*SB -0.7821E-01*SB2+0.6405E-01*SB3 |
| 545 | goto 100 |
| 546 | C Ifl = 0 |
| 547 | 63 A0=Exp(-0.3410E+00-0.9613E+00*SB -0.4969E+00*SB2+0.9360E-01*SB3) |
| 548 | A1=-0.2400E+00+0.1473E+00*SB -0.1593E+00*SB2+0.4538E-01*SB3 |
| 549 | A2= 0.4841E+01+0.9311E+00*SB +0.1601E-03*SB2-0.1331E+00*SB3 |
| 550 | A3= 0.7427E+01-0.1397E+01*SB +0.1489E+00*SB2-0.2848E+00*SB3 |
| 551 | A4= 0.9600E+00+0.3697E+00*SB -0.4246E+00*SB2+0.1032E+00*SB3 |
| 552 | A5= 0.0000E+00+0.2484E+01*SB -0.9908E+00*SB2+0.2568E+00*SB3 |
| 553 | goto 100 |
| 554 | C Ifl = -1 |
| 555 | 64 A0=Exp( 0.1176E+00-0.3418E+01*SB +0.3529E+01*SB2-0.1367E+01*SB3) |
| 556 | A1=-0.3654E+00+0.1914E+00*SB -0.2192E+00*SB2+0.6933E-01*SB3 |
| 557 | A2= 0.1099E+02-0.4281E+01*SB +0.3729E+01*SB2-0.1254E+01*SB3 |
| 558 | A3=-0.7514E+00+0.7696E+00*SB -0.1134E+01*SB2+0.4245E+00*SB3 |
| 559 | A4= 0.7690E-01-0.6558E-01*SB +0.8726E-01*SB2-0.3345E-01*SB3 |
| 560 | A5=-0.1447E+01+0.2617E+01*SB -0.2094E+01*SB2+0.7536E+00*SB3 |
| 561 | goto 100 |
| 562 | C Ifl = -2 |
| 563 | 65 A0=Exp(-0.2412E+01-0.2522E+01*SB +0.3126E+01*SB2-0.1305E+01*SB3) |
| 564 | A1=-0.2353E+00+0.3118E+00*SB -0.4864E+00*SB2+0.1689E+00*SB3 |
| 565 | A2= 0.9017E+01-0.2437E+01*SB +0.4659E+01*SB2-0.2044E+01*SB3 |
| 566 | A3= 0.5000E+02-0.1158E+02*SB -0.9260E+01*SB2+0.2847E+01*SB3 |
| 567 | A4= 0.1726E+01-0.6849E+00*SB +0.7864E+00*SB2-0.3300E+00*SB3 |
| 568 | A5= 0.5080E-01+0.2858E+01*SB -0.3297E+01*SB2+0.1246E+01*SB3 |
| 569 | goto 100 |
| 570 | C Ifl = -3 |
| 571 | 66 A0=Exp(-0.1966E+01-0.4405E+01*SB +0.2436E+00*SB2+0.4576E-01*SB3) |
| 572 | A1=-0.4000E-02-0.1229E+01*SB +0.1118E+01*SB2-0.2988E+00*SB3 |
| 573 | A2= 0.6902E+01+0.1266E+01*SB -0.1068E+01*SB2+0.3062E-01*SB3 |
| 574 | A3= 0.0000E+00+0.3987E+01*SB +0.9389E+01*SB2+0.1881E+02*SB3 |
| 575 | A4= 0.1000E-02+0.3528E+00*SB -0.4201E+00*SB2+0.1248E+00*SB3 |
| 576 | A5= 0.0000E+00-0.2149E+01*SB +0.2925E+01*SB2-0.7609E+00*SB3 |
| 577 | goto 100 |
| 578 | C Ifl = -4 |
| 579 | 67 A0=SB** 0.7561E+00*Exp(-0.6960E+01+0.5634E-01*SB -0.1170E+01*SB2) |
| 580 | A1=-0.4232E+00+0.9269E+00*SB -0.1161E+01*SB2+0.3470E+00*SB3 |
| 581 | A2= 0.6057E+01-0.5790E+01*SB +0.7352E+01*SB2-0.2435E+01*SB3 |
| 582 | A3= 0.2941E+02-0.1999E+02*SB -0.8345E+00*SB2+0.3000E+02*SB3 |
| 583 | A4= 0.1069E+00-0.4620E+00*SB +0.5614E+00*SB2-0.1336E+00*SB3 |
| 584 | A5=-0.1865E+00+0.3953E+01*SB -0.3791E+01*SB2+0.1315E+01*SB3 |
| 585 | goto 100 |
| 586 | C Ifl = -5 |
| 587 | 68 A0=SB** 0.5661E-02*Exp(-0.2123E+01-0.3026E+01*SB +0.1912E+00*SB2) |
| 588 | A1=-0.2011E+00-0.1338E-01*SB -0.3974E-01*SB2+0.1948E-01*SB3 |
| 589 | A2= 0.4906E+01+0.1740E+01*SB -0.1387E+01*SB2+0.1263E+00*SB3 |
| 590 | A3=-0.1000E+01+0.5767E-01*SB +0.6377E+00*SB2+0.4736E-01*SB3 |
| 591 | A4= 0.5927E-04+0.1039E+00*SB -0.9797E-01*SB2+0.6881E-01*SB3 |
| 592 | A5= 0.4017E+00+0.1981E+01*SB -0.7758E+00*SB2+0.2916E+00*SB3 |
| 593 | goto 100 |
| 594 | C Ifl = -6 |
| 595 | 69 A0=SB** 0.1008E+01*Exp(-0.7211E+01+0.3273E+01*SB -0.6979E+01*SB2) |
| 596 | A1=-0.1026E+00-0.4948E+00*SB +0.1188E+01*SB2-0.8016E+00*SB3 |
| 597 | A2= 0.5397E+01+0.2135E+01*SB -0.9531E+01*SB2+0.6115E+01*SB3 |
| 598 | A3= 0.4966E+01-0.1111E+02*SB +0.4732E+01*SB2+0.1568E+02*SB3 |
| 599 | A4= 0.5345E+00-0.1935E+00*SB +0.5816E+00*SB2-0.6794E+00*SB3 |
| 600 | A5= 0.3569E+01-0.3477E+01*SB +0.8756E+01*SB2-0.4139E+01*SB3 |
| 601 | goto 100 |
| 602 | |
| 603 | 311 stop 'This option is not currently supported.' |
| 604 | |
| 605 | 100 CONTINUE |
| 606 | P012 = A0 *(x**A1) *((1.-x)**A2) |
| 607 | P34 = (1.+A3*(x**A4)) |
| 608 | P5 = (log(1.+1./x))**A5 |
| 609 | C |
| 610 | PDFD = P012 * P34 * P5 |
| 611 | PDF = PDFD |
| 612 | C |
| 613 | Ctq2pPd = PDF |
| 614 | if(Ctq2pPd.lt.0.0) then |
| 615 | Ctq2pPd = 0.0 |
| 616 | Irt=1 |
| 617 | endif |
| 618 | |
| 619 | Ist = Iset |
| 620 | |
| 621 | Lp = Iprtn |
| 622 | Qsto = QQ |
| 623 | |
| 624 | Return |
| 625 | C ----------------------- |
| 626 | ENTRY Wlamd2p (Iset, Iorder, Neff) |
| 627 | |
| 628 | Iorder = IORD (Iset) |
| 629 | Wlamd2p = VLM (Neff, Iset) |
| 630 | |
| 631 | Ctq2pPd = 0.0 |
| 632 | RETURN |
| 633 | C ----------------------- |
| 634 | Entry PrCtq2p |
| 635 | > (Iset, Iordr, Ischeme, MxFlv, |
| 636 | > Alam4, Alam5, Alam6, Amas4, Amas5, Amas6, |
| 637 | > Xmin, Qini, Qmax, ExpNor) |
| 638 | |
| 639 | C Return QCD parameters and Fitting parameters |
| 640 | C associated with parton distribution set Iset. |
| 641 | C Iord : Order Of Fit |
| 642 | C Ischeme : (0, 1, 2) for (LO, MS-bar-NLO, DIS-NLO) resp. |
| 643 | C MxFlv : Maximum number of flavors included |
| 644 | C Alam_i : i = 4,5,6 Effective lambda for i-flavors |
| 645 | |
| 646 | C Amas_i : i = 4,5,6 Mass parameter for flavor i |
| 647 | C Xmin, Qini, Qmax : self explanary |
| 648 | C ExpNor(I) : Normalization factor for the experimental data set used in |
| 649 | C obtaining the best global fit for parton distributions Iset: |
| 650 | C I = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 |
| 651 | C BCDMS NMC90 NMC280 CCFR E605 WA70 E706 UA6 H1 ZEUS |
| 652 | |
| 653 | Iordr = Iord (Iset) |
| 654 | Ischeme= Isch (Iset) |
| 655 | MxFlv = Nqrk (Iset) |
| 656 | |
| 657 | Alam4 = Vlm(4,Iset) |
| 658 | Alam5 = Vlm(5,Iset) |
| 659 | Alam6 = Vlm(6,Iset) |
| 660 | |
| 661 | Amas4 = Qms(4,Iset) |
| 662 | Amas5 = Qms(5,Iset) |
| 663 | Amas6 = Qms(6,Iset) |
| 664 | |
| 665 | Xmin = Xmn (Iset) |
| 666 | Qini = Qmn (Iset) |
| 667 | Qmax = Qmx (Iset) |
| 668 | |
| 669 | Do 201 Iexp = 1, Nexp(Iset) |
| 670 | 201 ExpNor(Iexp) = ExpN(Iexp, Iset) |
| 671 | |
| 672 | Ctq2pPd = 0.0 |
| 673 | Return |
| 674 | C ************************* |
| 675 | END |
git://git.uio.no / u / mrichter / AliRoot.git / blame