]>
Commit | Line | Data |
---|---|---|
21886bb6 | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.2 1996/10/30 08:31:01 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 STRDO1(X,SCALE,UPV,DNV,SEA,STR,CHM,GL) | |
15 | C :::::::::::: DUKE OWENS SET 1 :::::::::::::::::::::::::::: | |
16 | #include "pdf/impdp.inc" | |
17 | #include "pdf/expdp.inc" | |
18 | + F(5),A(6,5),B1(3,6,5) | |
19 | DATA Q0,QL1/2.D0,.2D0/ | |
20 | DATA B1/3.D0,0.D0,0.D0,.419D0,.004383D0,-.007412D0, | |
21 | &3.46D0,.72432D0,-.065998D0,4.4D0,-4.8644D0,1.3274D0, | |
22 | &6*0.D0,1.D0, | |
23 | &0.D0,0.D0,.763D0,-.23696D0,.025836D0,4.D0,.62664D0,-.019163D0, | |
24 | &0.D0,-.42068D0,.032809D0,6*0.D0,1.265D0,-1.1323D0,.29268D0, | |
25 | &0.D0,-.37162D0,-.028977D0,8.05D0,1.5877D0,-.15291D0, | |
26 | &0.D0,6.3059D0,-.27342D0,0.D0,-10.543D0,-3.1674D0, | |
27 | &0.D0,14.698D0,9.798D0,0.D0,.13479D0,-.074693D0, | |
28 | &-.0355D0,-.22237D0,-.057685D0,6.3494D0,3.2649D0,-.90945D0, | |
29 | &0.D0,-3.0331D0,1.5042D0,0.D0,17.431D0,-11.255D0, | |
30 | &0.D0,-17.861D0,15.571D0,1.564D0,-1.7112D0,.63751D0, | |
31 | &0.D0,-.94892D0,.32505D0,6.D0,1.4345D0,-1.0485D0, | |
32 | &9.D0,-7.1858D0,.25494D0,0.D0,-16.457D0,10.947D0, | |
33 | &0.D0,15.261D0,-10.085D0/ | |
34 | WN=1.D0 | |
35 | S= LOG( LOG( MAX(Q0,SCALE)/QL1)/ LOG(Q0/QL1)) | |
36 | DO 10 I=1,5 | |
37 | DO 10 J=1,6 | |
38 | 10 A(J,I)=B1(1,J,I)+S*(B1(2,J,I)+S*B1(3,J,I)) | |
39 | DO 40 I=1,5 | |
40 | 40 F(I)=A(1,I)*X**A(2,I)*(WN-X)**A(3,I)*(WN+X* | |
41 | & (A(4,I)+X*(A(5,I)+X*A(6,I)))) | |
42 | DO 50 I=1,2 | |
43 | AA=WN+A(2,I)+A(3,I) | |
44 | 50 F(I)=F(I)*GAMFUN(AA)/((WN+A(2,I)*A(4,I)/AA)*GAMFUN(A(2,I)) | |
45 | &*GAMFUN(WN+A(3,I))) | |
46 | UPV=F(1)-F(2) | |
47 | DNV=F(2) | |
48 | SEA=F(3)/6.D0 | |
49 | STR=SEA | |
50 | CHM=F(4) | |
51 | GL =F(5) | |
52 | RETURN | |
53 | END |