PDF/spdf/strdo2.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:31:02  cernlib
   9 * Version 7.04
  10 *
  11 * Revision 1.1.1.1  1996/04/12 15:29:49  plothow
  12 * Version 7.01
  13 *
  14 *
  15 #include "pdf/pilot.h"
  16 C-----------------------------------------------------------------------
  17       SUBROUTINE STRDO2(X,SCALE,UPV,DNV,SEA,STR,CHM,GL)
  18 C :::::::::::: DUKE OWENS SET 2 ::::::::::::::::::::::::::::
  19 C
  20 #include "pdf/impdp.inc"
  21 #if defined(CERNLIB_DOUBLE)
  22       DOUBLE PRECISION
  23 #endif
  24 #if defined(CERNLIB_SINGLE)
  25       REAL
  26 #endif
  27      +       F(5),A(6,5),B2(3,6,5)
  28       DATA Q0,QL2/2.D0,.4D0/
  29       DATA B2/3.D0,0.D0,0.D0,.3743D0,.013946D0,-.00031695D0,
  30      &3.329D0,.75343D0,-.076125D0,6.032D0,-6.2153D0,1.5561D0,
  31      &6*0.D0,1.D0,0.D0,
  32      &0.D0,.7608D0,-.2317D0,.023232D0,3.83D0,.62746D0,-.019155D0,
  33      &0.D0,-.41843D0,.035972D0,6*0.D0,1.6714D0,-1.9168D0,.58175D0,
  34      &0.D0,-.27307D0,-.16392D0,9.145D0,.53045D0,-.76271D0,
  35      &0.D0,15.665D0,-2.8341D0,0.D0,-100.63D0,44.658D0,
  36      &0.D0,223.24D0,-116.76D0,0.D0,.067368D0,-.030574D0,
  37      &-.11989D0,-.23293D0,-.023273D0,3.5087D0,3.6554D0,-.45313D0,
  38      &0.D0,-.47369D0,.35793D0,0.D0,9.5041D0,-5.4303D0,
  39      &0.D0,-16.563D0,15.524D0,.8789D0,-.97093D0,.43388D0,
  40      &0.D0,-1.1612D0,.4759D0,4.D0,1.2271D0,-.25369D0,
  41      &9.D0,-5.6354D0,-.81747D0,0.D0,-7.5438D0,5.5034D0,
  42      &0.D0,-.59649D0,.12611D0/
  43       WN=1.D0
  44       S= LOG( LOG( MAX(Q0,SCALE)/QL2)/ LOG(Q0/QL2))
  45       DO 10 I=1,5
  46          DO 10 J=1,6
  47    10 A(J,I)=B2(1,J,I)+S*(B2(2,J,I)+S*B2(3,J,I))
  48       DO 20 I=1,5
  49    20 F(I)=A(1,I)*X**A(2,I)*(WN-X)**A(3,I)*(WN+X* (A(4,I)+X*(A(5,I)+X*A
  50      +(6,I))))
  51       DO 30 I=1,2
  52          AA=WN+A(2,I)+A(3,I)
  53    30 F(I)=F(I)*GAMFUN(AA)/((WN+A(2,I)*A(4,I)/AA)*GAMFUN(A(2,I))
  54      +*GAMFUN(WN+A(3,I)))
  55       UPV=F(1)-F(2)
  56       DNV=F(2)
  57       SEA=F(3)/6.D0
  58       STR=SEA
  59       CHM=F(4)
  60       GL =F(5)
  61       RETURN
  62       END