PDF/spdf/srctq23.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.4  2000/09/18 10:02:37  hristov
   6 * Makefile added to PDF8
   7 *
   8 * Revision 1.1.1.2  1996/10/30 08:30:55  cernlib
   9 * Version 7.04
  10 *
  11 * Revision 1.1.1.1  1996/04/12 15:29:48  plothow
  12 * Version 7.01
  13 *
  14 *
  15 #include "pdf/pilot.h"
  16 C-----------------------------------------------------------------------
  17       SUBROUTINE SRCTQ23(DX,DSCAL,
  18      +                     DUPV,DDNV,DSEA,DSEAD,DSTR,DCHM,DBOT,DTOP,DGL)
  19 C
  20 C   ********************************************************************
  21 C   *                                                                  *
  22 C   *        Interface to CTEQset of structure functions               *
  23 C   *                                                                  *
  24 C   *        Author:    H. Plothow-Besch (CERN-PPE)                    *
  25 C   *                                                                  *
  26 C   ********************************************************************
  27 C
  28 C :::::::::::: Structure functions from Morfin and Wu-Ki Tung
  29 C :::::::::::: Lambda = 0.208 GeV (MSb - singular gluon) (CTEQ2MS)
  30 C
  31 #if defined(CERNLIB_DOUBLE)
  32       DOUBLE PRECISION
  33 #endif
  34 #if defined(CERNLIB_SINGLE)
  35       REAL
  36 #endif
  37      +       DX,DSCAL,DUPV,DDNV,DSEA,DSEAD,DSTR,DCHM,DBOT,DTOP,DGL
  38       REAL   X, Q
  39       PARAMETER (ISET=2)
  40 #include "pdf/w50511.inc"
  41 C
  42       X = DX
  43       Q = DSCAL
  44 C
  45 C     generate the individual structure fcn calls
  46 C
  47       UPV = Ctq2pPd(ISET, 1,X,Q,IFL)
  48       DUPV = UPV
  49       DNV = Ctq2pPd(ISET, 2,X,Q,IFL)
  50       DDNV = DNV
  51       SEAU = Ctq2pPd(ISET,-1,X,Q,IFL)
  52       DSEA = SEAU
  53       SEAD = Ctq2pPd(ISET,-2,X,Q,IFL)
  54       DSEAD = SEAD
  55       STR = Ctq2pPd(ISET,-3,X,Q,IFL)
  56       DSTR = STR
  57       CHM = Ctq2pPd(ISET,-4,X,Q,IFL)
  58       DCHM = CHM
  59       BOT = Ctq2pPd(ISET,-5,X,Q,IFL)
  60       DBOT = BOT
  61       TOP = 0.
  62       IF (DSCAL.GT.TMAS) TOP = Ctq2pPd(ISET,-6,X,Q,IFL)
  63       DTOP = TOP
  64       GL = Ctq2pPd(ISET, 0,X,Q,IFL)
  65       DGL = GL
  66 C
  67       RETURN
  68       END