5 * Revision 1.1.1.3 2000/04/19 13:30:02 mclareni
8 * Revision 1.2 2000/04/09 10:06:56 plothow
9 * new version 8.00 - MRS 98 and 99, CTEQ5 and GRV98 added
11 * Revision 1.1.1.1 1996/04/12 15:30:17 plothow
15 #include "pdf/pilot.h"
18 PARAMETER (LFLPRT = 2)
19 #include "pdf/w5051p2.inc"
20 #include "pdf/w5051p7.inc"
21 #include "pdf/w50510.inc"
22 #include "pdf/w50511.inc"
23 #include "pdf/w50512.inc"
24 #include "pdf/w50515.inc"
25 #include "pdf/w50516.inc"
26 #include "pdf/w50517.inc"
27 #include "pdf/w505110.inc"
28 #include "pdf/w505120.inc"
29 #include "pdf/w505122.inc"
31 #include "pdf/expdp.inc"
34 REAL X, Q, UPV, DNV, USEA, DSEA, SSEA, CSEA, BOT, TOP, GLUE
36 CHARACTER*20 PARM(NCHDIM)
37 #include "pdf/expdp.inc"
40 DATA NXX,XX/5, .002,.140,.700,0.01,0.05/
41 DATA NQ2,Q2/5, 1.30,24.0,33.5,200.0,200.0/
51 C First call to PDFSET to initialize COMMON/W505120/
56 C Test new PDFLIB versions (4.0 or bigger)
59 IF(NPGSMX(KPTY,KGRO) .EQ. 0) GOTO 400
60 DO ISET=-9,NPGSMX(KPTY,KGRO)
61 IF((KPTY*KGRO).NE.1 .AND. ISET.LE.0) GOTO 300
62 WRITE (N6,*) ' ==========================================
63 +============================================'
65 IF(ISET.LE.-8) IFLSET = 0
73 C... test different NSET writings
74 IF(ISET.GT. 5) PARM(3) = 'NSET'
75 IF(ISET.GT.10) PARM(3) = 'nset'
76 IF(ISET.GT.15) PARM(3) = 'NSet'
77 IF(ISET.GT.20) PARM(3) = 'nsET'
79 C... test different NSET writings and INVALID parameter settings
83 ELSEIF (ISET.EQ.-9) THEN
85 ELSEIF (ISET.EQ.-8) THEN
89 ELSEIF (ISET.EQ.-7) THEN
93 ELSEIF (ISET.EQ.-6) THEN
97 ELSEIF (ISET.EQ.-5) THEN
100 CALL PDFSET(PARM,VAL)
101 ELSEIF (ISET.EQ.-4) THEN
108 CALL PDFSET(PARM,VAL)
109 ELSEIF (ISET.EQ.-3) THEN
116 CALL PDFSET(PARM,VAL)
117 ELSEIF (ISET.EQ.-2) THEN
124 CALL PDFSET(PARM,VAL)
125 ELSEIF (ISET.EQ.-1) THEN
132 CALL PDFSET(PARM,VAL)
137 WRITE(N6,4000) NPTYPE,NGROUP,NSET,
138 + SFNAME(NPTYPE,NGROUP,NSET),MODECR(NPTYPE,NGROUP,NSET)
139 WRITE(N6,4001) QCDL4,QCDL5,DQ,ALF
140 4000 FORMAT(/,' Nptype = ',I1,' Ngroup = ',I2,
141 + ' Nset = ',I3,' Name = "',A8,
143 4001 FORMAT(1H ,' QCDL4, QCDL5 = ',2F7.4,
144 + ' Scale = ',F8.4,' alpha(s) = ',F7.5)
152 CALL PFTOPDG(DX,DQ,DXPDF)
178 SFFP=CH1*(FFPU+FFPUB+FFPC+FFPCB)+CH2*(FFPD+FFPDB+FFPS+FFPSB)
179 SFFN=CH1*(FFNU+FFNUB+FFNC+FFNCB)+CH2*(FFND+FFNDB+FFNS+FFNSB)
181 IF (SFFP .NE. 0.) FACNP=SFFN/SFFP
191 WRITE(N6,3000)X,Q2(IX),FACNP,UP,UV,US,DP,DV,GL,F2P
192 3000 FORMAT(1X,' X:',F6.3,2X,'Q2:',F5.1,
193 1 3X,'N/P:',F7.4,2X,'UP:',F6.3,
194 1 1X,'UV:',F6.3,1X,'US:',F6.3,1X,'DP:',F6.3,1X,'DV:',F6.3,
195 2 1X,'GL:',F6.3,2X,'F2P:',F7.4)
202 WRITE(N6,4000) NPTYPE,NGROUP,NSET,
203 + SFNAME(NPTYPE,NGROUP,NSET),MODECR(NPTYPE,NGROUP,NSET)
204 WRITE(N6,4001) QCDL4,QCDL5,DQ,ALF