PDF/spdf/ctq2ppd.F

   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