PDF/spdf/sfsas21.F

   1 C-----------------------------------------------------------------------
   2 #include "pdf/pilot.h"
   3       SUBROUTINE SFSAS21(DX,DQ2,DP2,IP2,
   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 = 0.36 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=1)
  28 #include "pdf/w50511.inc"
  29 C
  30       X  = DX
  31       Q  = SQRT(DQ2)
  32       Q2 = DQ2
  33       P2 = DP2
  34       IP = IP2
  35       IF(IP.LT.0 .OR. IP.GT.7) IP=0
  36 C
  37 C     generate the individual structure fcn calls
  38 C
  39       CALL SASGAM2(ISET,X,Q2,P2,IP,F2GAM,XPDFGM)
  40       UPV = XPDFGM(2)
  41       DUPV = UPV
  42       DNV = XPDFGM(1)
  43       DDNV = DNV
  44       SEAU = XPDFGM(-2)
  45       DSEA = SEAU
  46       SEAD = XPDFGM(-1)
  47       DSEAD = SEAD
  48       STR = XPDFGM(-3)
  49       DSTR = STR
  50       CHM = XPDFGM(-4)
  51       DCHM = CHM
  52       BOT = XPDFGM(-5)
  53       DBOT = BOT
  54       TOP = 0.
  55 C      IF (DSCAL.GT.TMAS) TOP = XPDFGM(6)
  56       DTOP = TOP
  57       GL = XPDFGM(0)
  58       DGL = GL
  59 C
  60       RETURN
  61       END