PDF/spdf/grv94di.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:28:34  cernlib
   6 * Version 7.04
   7 *
   8 * Revision 1.1.1.1  1996/04/12 15:29:22  plothow
   9 * Version 7.01
  10 *
  11 *
  12 #include "pdf/pilot.h"
  13 *
  14 *...NLO PARAMETRIZATION (DIS) :
  15 *
  16        SUBROUTINE GRV94DI (ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZTB,ZGL)
  17        IMPLICIT REAL (A - Y)
  18 C+SEQ, IMPDP.
  19 #include "pdf/expdp.inc"
  20      +      ZX,ZQ,ZUV,ZDV,ZDEL,ZUDB,ZUB,ZDB,ZSB,ZCB,ZBB,ZTB,ZGL
  21        REAL X, Q
  22        X = ZX
  23        Q = ZQ
  24        MU2  = 0.34
  25        LAM2 = 0.248 * 0.248
  26        Q2 = Q*Q
  27        ARG1 = MAX(MU2,Q2)/LAM2
  28        ARG2 = MU2/LAM2
  29        S  =  LOG ( LOG( ARG1  ) /  LOG( ARG2   ))
  30        DS =  SQRT (S)
  31        S2 = S * S
  32        S3 = S2 * S
  33 *...UV :
  34        NU  =  2.484 + 0.116 * S + 0.093 * S2
  35        AKU =  0.563 - 0.025 * S
  36        BKU =  0.054 + 0.154 * S
  37        AU  = -0.326 - 0.058 * S - 0.135 * S2
  38        BU  = -3.322 + 8.259 * S - 3.119 * S2 + 0.291 * S3
  39        CU  =  11.52 - 12.99 * S + 3.161 * S2
  40        DU  =  2.808 + 1.400 * S - 0.557 * S2 + 0.119 * S3
  41        UV  = GR94FV (X, NU, AKU, BKU, AU, BU, CU, DU)
  42        ZUV = UV
  43 *...DV :
  44        ND  =  0.156 - 0.017 * S
  45        AKD =  0.299 - 0.022 * S
  46        BKD =  0.259 - 0.015 * S
  47        AD  =  3.445 + 1.278 * S + 0.326 * S2
  48        BD  = -6.934 + 37.45 * S - 18.95 * S2 + 1.463 * S3
  49        CD  =  55.45 - 69.92 * S + 20.78 * S2
  50        DD  =  3.577 + 1.441 * S - 0.683 * S2 + 0.179 * S3
  51        DV  = GR94FV (X, ND, AKD, BKD, AD, BD, CD, DD)
  52        ZDV = DV
  53 *...DEL :
  54        NE  =  0.099 + 0.019 * S + 0.002 * S2
  55        AKE =  0.419 - 0.013 * S
  56        BKE =  1.064 - 0.038 * S
  57        AE  = -44.00 + 98.70 * S - 14.79 * S2
  58        BE  =  28.59 - 40.94 * S - 13.66 * S2 + 2.523 * S3
  59        CE  =  84.57 - 108.8 * S + 31.52 * S2
  60        DE  =  7.469 + 2.480 * S - 0.866 * S2
  61        DEL = GR94FV (X, NE, AKE, BKE, AE, BE, CE, DE)
  62        ZDEL= DEL
  63 *...UDB :
  64        ALX =  1.215
  65        BEX =  0.466
  66        AKX =  0.326 + 0.150 * S
  67        BKX =  0.956 + 0.405 * S
  68        AGX =  0.272
  69        BGX =  3.794 - 2.359 * DS
  70        CX  =  2.014
  71        DX  =  7.941 + 0.534 * DS - 0.940 * S + 0.410 * S2
  72        EX  =  3.049 + 1.597 * S
  73        ESX =  4.396 - 4.594 * DS + 3.268 * S
  74        UDB = GR94FW(X, S, ALX, BEX, AKX, BKX, AGX, BGX, CX, DX, EX, ESX)
  75        ZUDB= UDB
  76 *...UB
  77        ZUB = (ZUDB - ZDEL)/2.D0
  78 *...DB
  79        ZDB = (ZUDB + ZDEL)/2.D0
  80 *...SB :
  81        ALS =  0.175
  82        BES =  0.344
  83        AKS =  1.415 - 0.641 * DS
  84        AS  =  0.580 - 9.763 * DS + 6.795 * S  - 0.558 * S2
  85        BS  =  5.617 + 5.709 * DS - 3.972 * S
  86        DST =  13.78 - 9.581 * S  + 5.370 * S2 - 0.996 * S3
  87        EST =  4.546 + 0.372 * S2
  88        ESS =  5.053 - 1.070 * S  + 0.805 * S2
  89        SB  = GR94FWS (X, S, ALS, BES, AKS, AS, BS, DST, EST, ESS)
  90        ZSB = SB
  91 *...GL :
  92        ALG =  1.258
  93        BEG =  1.846
  94        AKG =  2.423
  95        BKG =  2.427 + 1.311 * S  - 0.153 * S2
  96        AG  =  25.09 - 7.935 * S
  97        BG  = -14.84 - 124.3 * DS + 72.18 * S
  98        CG  =  590.3 - 173.8 * S
  99        DG  =  5.196 + 1.857 * S
 100        EG  = -1.648 + 3.988 * S  - 0.432 * S2
 101        ESG =  3.232 - 0.542 * S
 102        GL  = GR94FW (X, S, ALG, BEG, AKG, BKG, AG, BG, CG, DG, EG, ESG)
 103        ZGL = GL
 104 *...
 105        ZCB = 0.D0
 106        ZBB = 0.D0
 107        ZTB = 0.D0
 108 *...
 109        RETURN
 110        END