PDF/spdf/grvho.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:28:42  cernlib
   6 * Version 7.04
   7 *
   8 * Revision 1.1.1.1  1996/04/12 15:29:23  plothow
   9 * Version 7.01
  10 *
  11 *
  12 #include "pdf/pilot.h"
  13 C
  14 C...HO PARAMATRIZATIONS :
  15        SUBROUTINE GRVHO (ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL)
  16        IMPLICIT REAL (A - Y)
  17 #include "pdf/expdp.inc"
  18      +        ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL
  19        REAL  X, Q
  20        X = ZX
  21        Q = ZQ
  22        MU2  = 0.3
  23        LAM2 = 0.248 * 0.248
  24        Q2 = Q*Q
  25        S  = ALOG (ALOG( MAX(MU2,Q2)/LAM2) / ALOG(MU2/LAM2))
  26        DS = SQRT (S)
  27        S2 = S * S
  28        S3 = S2 * S
  29 C...X * (UV + DV) :
  30        NUD  = 0.330 + 0.151 * S - 0.059 * S2 + 0.027 * S3
  31        AKUD = 0.285
  32        AGUD = -2.28 + 15.73 * S -  4.58 * S2
  33        BUD  =  56.7 -  53.6 * S + 11.21 * S2
  34        DUD  =  3.17 +  1.17 * S -  0.47 * S2 +  0.09 * S3
  35        UDV  = GRVFV (X, NUD, AKUD, AGUD, BUD, DUD)
  36 C...X * DV :
  37        ND  = 0.459 + 0.315 * DS + 0.515 * S
  38        AKD = 0.624              - 0.031 * S
  39        AGD =  8.13 -  6.77 * DS +  0.46 * S
  40        BD  =  6.59 - 12.83 * DS +  5.65 * S
  41        DD  =  3.98              +  1.04 * S  -  0.34 * S2
  42        DV  = GRVFV (X, ND, AKD, AGD, BD, DD)
  43        ZDV = DV
  44 C...X * UV :
  45        UV = UDV - DV
  46        ZUV = UV
  47 C...X * G :
  48        ALG =  1.128
  49        BEG =  1.575
  50        AKG =  0.323 + 1.653 * S
  51        BKG =  0.811 + 2.044 * S
  52        AGG =   0.0  + 1.963 * S - 0.519 * S2
  53        BGG =  0.078 +  6.24 * S
  54        CG  =  30.77 - 24.19 * S
  55        DG  =  3.188 + 0.720 * S
  56        EG  = -0.881 + 2.687 * S
  57        ESG =  2.466
  58        GL  = GRVFW (X, S, ALG, BEG, AKG, BKG, AGG, BGG, CG, DG, EG, ESG)
  59        ZGL = GL
  60 C...X * UBAR = X * DBAR :
  61        ALU =  0.594
  62        BEU =  0.614
  63        AKU =  0.636 - 0.084 * S
  64        BKU =   0.0
  65        AGU =  1.121 - 0.193 * S
  66        BGU =  0.751 - 0.785 * S
  67        CU  =   8.57 - 1.763 * S
  68        DU  =  10.22 + 0.668 * S
  69        EU  =  3.784 + 1.280 * S
  70        ESU =  1.808 + 0.980 * S
  71        UDB = GRVFW (X, S, ALU, BEU, AKU, BKU, AGU, BGU, CU, DU, EU, ESU)
  72        ZUDB = UDB
  73 C...X * SBAR = X * S :
  74        SS  =   0.0
  75        ALS =  0.756
  76        BES =  0.101
  77        AKS =  2.942 - 1.016 * S
  78        AGS =  -4.60 + 1.167 * S
  79        BS  =   9.31 - 1.324 * S
  80        DS  =  11.49 - 1.198 * S + 0.053 * S2
  81        EST =  2.630 + 1.729 * S
  82        ESS =   8.12
  83        SB  = GRVFWS (X, S, SS, ALS, BES, AKS, AGS, BS, DS, EST, ESS)
  84        ZSB = SB
  85 C...X * CBAR = X * C :
  86        SC  =  0.820
  87        ALC =   0.98
  88        BEC =   0.0
  89        AKC = -0.625 - 0.523 * S
  90        AGC =   0.0
  91        BC  =  1.896 + 1.616 * S
  92        DC  =   4.12 + 0.683 * S
  93        EC  =   4.36 + 1.328 * S
  94        ESC =  0.677 + 0.679 * S
  95        CB  = GRVFWS (X, S, SC, ALC, BEC, AKC, AGC, BC, DC, EC, ESC)
  96        ZCB = CB
  97 C...X * BBAR = X * B :
  98        SBO =  1.297
  99        ALB =   0.99
 100        BEB =   0.0
 101        AKB =   0.0  - 0.193 * S
 102        AGB =   0.0
 103        BBO =   0.0
 104        DB  =  3.447 + 0.927 * S
 105        EB  =   4.68 + 1.259 * S
 106        ESB =  1.892 + 2.199 * S
 107        BB  = GRVFWS (X, S, SBO, ALB, BEB, AKB, AGB, BBO, DB, EB, ESB)
 108        ZBB = BB
 109 C...X * TBAR = X * T :
 110        TB=0.
 111        ZTB = TB
 112        RETURN
 113        END