*
* $Id$
*
* $Log$
* Revision 1.1.1.2  1996/10/30 08:32:57  cernlib
* Version 7.04
*
* Revision 1.1.1.1  1996/04/12 15:30:18  plothow
* Version 7.01
*
*
#include "pdf/pilot.h"
      PROGRAM PDFUPV
C
      PARAMETER (LFLPRT = 2)
#include "pdf/w5051p2.inc"
#include "pdf/w5051p7.inc"
#include "pdf/w50510.inc"
#include "pdf/w50511.inc"
#include "pdf/w505110.inc"
#include "pdf/w505120.inc"
#include "pdf/w505122.inc"
C   define HBOOK settings
      PARAMETER (NHBMEM = 500000)
      COMMON /PAWC/ HMEM(NHBMEM)
      PARAMETER (NB=900, ID=  0)
C   define DOUBLE PRECISION variables for calling sequences
#include "pdf/expdp.inc"
     +       DX,DSCALE,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL,
     +       ALF,ALPHAS2
      REAL    X, SCALE, UPV, DNV, USEA, DSEA, STR, CHM, BOT, TOP, GL
#include "pdf/w50516.inc"
      CHARACTER*20 PARM(NCHDIM)
#include "pdf/expdp.inc"
     +       VAL(NCHDIM)
C
      DATA SCALE/80.140/
      DATA X00/0.002/, DX0/0.001/
C
      N6 = L6
      IFLPRT = LFLPRT
C
      CALL HLIMIT(NHBMEM)
C
      DSCALE=SCALE
C
C First call to PDFSET to initialize COMMON/W505120/
      PARM(1) = 'Init0'
      VAL(1)  = 0.D0
      CALL PDFSET(PARM,VAL)
C   loop over all existing sets of Nucleon structure functions (SF)
      NPTYPE = 1
      NSF = 0
      DO IGR = 1,NGRMAX
         IF(NPGSMX(NPTYPE,IGR).EQ.0) GOTO 30
         DO ISET = 1,NPGSMX(NPTYPE,IGR)
            WRITE (N6,*) ' =============================================
     +================================================================='
            NSF = NSF + 1
C   book histograms for each set of SF separately
            CALL HBOOK1(ID+NSF,'U Valence$',NB,0.0015,0.9015,0.)
C   force label printing for each set of SF (not only 1st)
            FIRST = .TRUE.
C   select and set parameters
            PARM(1) = 'Nptype'
            VAL(1)  = NPTYPE
            PARM(2) = 'Ngroup'
            VAL(2)  = IGR
            PARM(3) = 'Nset'
            VAL(3)  = ISET
            CALL PDFSET(PARM,VAL)
C   loop over all x bins
            DO 10 I=1,NB
               X= X00 +  (I-1)*DX0
               DX=X
               IF(X.LT.0.0015 .OR. X.GT.0.9015) GOTO 10
               CALL STRUCTM(DX,DSCALE,
     +                    DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
               UPV=DUPV
              IF(X.GT.0.499 .AND. X.LE.0.500) WRITE(N6,1000) X,SCALE,UPV
 1000          FORMAT(/,'  X= ',F6.4,'  Q= ',F6.3,'  UPV= ',F8.4)
               CALL HF1(ID+NSF,X,UPV)
 10         CONTINUE
C   get alpha(s) for selected set of SF at Q = SCALE
            ALF = ALPHAS2(DSCALE)
           WRITE(N6,4000) NPTYPE,NGROUP,NSET,SFNAME(NPTYPE,NGROUP,NSET),
     +                    MODECR(NPTYPE,NGROUP,NSET),NSF
            WRITE(N6,4001) DSCALE,ALF
 4000       FORMAT(/,'  Nptype = ',I1,'  Ngroup = ',I1,
     +               '  Nset = ',I2,'  Name  = "',A8,
     +               '"  CrMode = ',I3,'  HBId = ',I3)
 4001       FORMAT(1H ,' Scale  = ',F8.4,'  alpha(s) = ',F6.4)
C   get error summary for each set of SF
            CALL PDFSTA
   20       CONTINUE
         ENDDO
   30    CONTINUE
      ENDDO
C
      STOP
      END