PDF/spdf/sfsas13.F

   1 C-----------------------------------------------------------------------
   2 #include "pdf/pilot.h"
   3       SUBROUTINE SFSAS13(DX,DQ2,DP2,
   4      +                     DUPV,DDNV,DSEA,DSEAD,DSTR,DCHM,DBOT,DTOP,DGL)
   5 C
   6 C   ********************************************************************
   7 C   *                                                                  *
   8 C   *        Interface to SASset of structure functions                *
   9 C   *                                                                  *
  10 C   *        Author:    H. Plothow-Besch (CERN-PPE)                    *
  11 C   *                                                                  *
  12 C   ********************************************************************
  13 C
  14 C :::::::::::: Structure functions from the SAS group version 2
  15 C :::::::::::: Lambda = 0.200 GeV, Q**2 = 4 GeV**2 (DIS)
  16 C
  17 #if defined(CERNLIB_DOUBLE)
  18       DOUBLE PRECISION
  19 #endif
  20 #if defined(CERNLIB_SINGLE)
  21       REAL
  22 #endif
  23      +          DX,DQ2,DP2,
  24      +          DUPV,DDNV,DSEA,DSEAD,DSTR,DCHM,DBOT,DTOP,DGL
  25       DIMENSION XPDFGM(-6:6)
  26       REAL X, Q, Q2, P2, F2GAM, XPDFGM
  27       PARAMETER (ISET=3)
  28 #include "pdf/w50511.inc"
  29 C
  30       X  = DX
  31       Q  = SQRT(DQ2)
  32       Q2 = DQ2
  33       P2 = DP2
  34 C
  35 C     generate the individual structure fcn calls
  36 C
  37       CALL SASGAM1(ISET,X,Q2,P2,F2GAM,XPDFGM)
  38       UPV = XPDFGM(2)
  39       DUPV = UPV
  40       DNV = XPDFGM(1)
  41       DDNV = DNV
  42       SEAU = XPDFGM(-2)
  43       DSEA = SEAU
  44       SEAD = XPDFGM(-1)
  45       DSEAD = SEAD
  46       STR = XPDFGM(-3)
  47       DSTR = STR
  48       CHM = XPDFGM(-4)
  49       DCHM = CHM
  50       BOT = XPDFGM(-5)
  51       DBOT = BOT
  52       TOP = 0.
  53 C      IF (DSCAL.GT.TMAS) TOP = XPDFGM(6)
  54       DTOP = TOP
  55       GL = XPDFGM(0)
  56       DGL = GL
  57 C
  58       RETURN
  59       END