5 * Revision 1.1.1.2 1996/10/30 08:32:57 cernlib
8 * Revision 1.1.1.1 1996/04/12 15:30:18 plothow
12 #include "pdf/pilot.h"
15 PARAMETER (LFLPRT = 2)
16 #include "pdf/w5051p2.inc"
17 #include "pdf/w5051p7.inc"
18 #include "pdf/w50510.inc"
19 #include "pdf/w50511.inc"
20 #include "pdf/w505110.inc"
21 #include "pdf/w505120.inc"
22 #include "pdf/w505122.inc"
23 C define HBOOK settings
24 PARAMETER (NHBMEM = 500000)
25 COMMON /PAWC/ HMEM(NHBMEM)
26 PARAMETER (NB=900, ID= 0)
27 C define DOUBLE PRECISION variables for calling sequences
28 #include "pdf/expdp.inc"
29 + DX,DSCALE,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL,
31 REAL X, SCALE, UPV, DNV, USEA, DSEA, STR, CHM, BOT, TOP, GL
32 #include "pdf/w50516.inc"
33 CHARACTER*20 PARM(NCHDIM)
34 #include "pdf/expdp.inc"
38 DATA X00/0.002/, DX0/0.001/
47 C First call to PDFSET to initialize COMMON/W505120/
51 C loop over all existing sets of Nucleon structure functions (SF)
55 IF(NPGSMX(NPTYPE,IGR).EQ.0) GOTO 30
56 DO ISET = 1,NPGSMX(NPTYPE,IGR)
57 WRITE (N6,*) ' =============================================
58 +================================================================='
60 C book histograms for each set of SF separately
61 CALL HBOOK1(ID+NSF,'U Valence$',NB,0.0015,0.9015,0.)
62 C force label printing for each set of SF (not only 1st)
64 C select and set parameters
72 C loop over all x bins
76 IF(X.LT.0.0015 .OR. X.GT.0.9015) GOTO 10
77 CALL STRUCTM(DX,DSCALE,
78 + DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
80 IF(X.GT.0.499 .AND. X.LE.0.500) WRITE(N6,1000) X,SCALE,UPV
81 1000 FORMAT(/,' X= ',F6.4,' Q= ',F6.3,' UPV= ',F8.4)
82 CALL HF1(ID+NSF,X,UPV)
84 C get alpha(s) for selected set of SF at Q = SCALE
86 WRITE(N6,4000) NPTYPE,NGROUP,NSET,SFNAME(NPTYPE,NGROUP,NSET),
87 + MODECR(NPTYPE,NGROUP,NSET),NSF
88 WRITE(N6,4001) DSCALE,ALF
89 4000 FORMAT(/,' Nptype = ',I1,' Ngroup = ',I1,
90 + ' Nset = ',I2,' Name = "',A8,
91 + '" CrMode = ',I3,' HBId = ',I3)
92 4001 FORMAT(1H ,' Scale = ',F8.4,' alpha(s) = ',F6.4)
93 C get error summary for each set of SF