PDF/spdf/strdo2.F

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