]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PDF/tpdf/tpdfupv.F
Removing obsolete module PDF
[u/mrichter/AliRoot.git] / PDF / tpdf / tpdfupv.F
diff --git a/PDF/tpdf/tpdfupv.F b/PDF/tpdf/tpdfupv.F
deleted file mode 100644 (file)
index 6827a2b..0000000
+++ /dev/null
@@ -1,101 +0,0 @@
-*
-* $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