Coding conventions

[u/mrichter/AliRoot.git] / PDF / spdf / sfctq23.F

*
* $Id$
*
* $Log$
* Revision 1.4  2000/09/18 10:02:36  hristov
* Makefile added to PDF8
*
* Revision 1.1.1.2  1996/10/30 08:29:33  cernlib
* Version 7.04
*
* Revision 1.1.1.1  1996/04/12 15:29:33  plothow
* Version 7.01
*
*
#include "pdf/pilot.h"
C-----------------------------------------------------------------------
      SUBROUTINE SFCTQ23(DX,DSCAL,
     +                     DUPV,DDNV,DSEA,DSEAD,DSTR,DCHM,DBOT,DTOP,DGL)
C
C   ********************************************************************
C   *                                                                  *
C   *        Interface to CTEQset of structure functions               *
C   *                                                                  *
C   *        Author:    H. Plothow-Besch (CERN-PPE)                    *
C   *                                                                  *
C   ********************************************************************
C
C :::::::::::: Structure functions from Morfin and Wu-Ki Tung
C :::::::::::: Lambda = 0.208 GeV (MSb - singular gluon) (CTEQ2MS)
C
#if defined(CERNLIB_DOUBLE)
      DOUBLE PRECISION
#endif
#if defined(CERNLIB_SINGLE)
      REAL
#endif
     +       DX,DSCAL,DUPV,DDNV,DSEA,DSEAD,DSTR,DCHM,DBOT,DTOP,DGL
      REAL   X, Q
      PARAMETER (ISET=2)
#include "pdf/w50511.inc"
C
      X = DX
      Q = DSCAL
C
C     generate the individual structure fcn calls
C
      UPV = Ctq2OPd(ISET, 1,X,Q,IFL)
      DUPV = UPV
      DNV = Ctq2OPd(ISET, 2,X,Q,IFL)
      DDNV = DNV
      SEAU = Ctq2OPd(ISET,-1,X,Q,IFL)
      DSEA = SEAU
      SEAD = Ctq2OPd(ISET,-2,X,Q,IFL)
      DSEAD = SEAD
      STR = Ctq2OPd(ISET,-3,X,Q,IFL)
      DSTR = STR
      CHM = Ctq2OPd(ISET,-4,X,Q,IFL)
      DCHM = CHM
      BOT = Ctq2OPd(ISET,-5,X,Q,IFL)
      DBOT = BOT
      TOP = 0.
      IF (DSCAL.GT.TMAS) TOP = Ctq2OPd(ISET,-6,X,Q,IFL)
      DTOP = TOP
      GL = Ctq2OPd(ISET, 0,X,Q,IFL)
      DGL = GL
C
      RETURN
      END