]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 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 |