PDF/spdf/pdfsta.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.3  2000/04/19 13:28:20  mclareni
   6 * Import version 8.04
   7 *
   8 * Revision 1.2  2000/04/09 10:06:26  plothow
   9 * new version 8.00 - MRS 98 and 99, CTEQ5 and GRV98 added
  10 *
  11 * Revision 1.1.1.1  1996/04/12 15:29:25  plothow
  12 * Version 7.01
  13 *
  14 *
  15 #include "pdf/pilot.h"
  16       SUBROUTINE PDFSTA
  17 C
  18 C **********************************************************************
  19 C *                                                                    *
  20 C *   Routine to print status of parameters in the COMMON/W50511/,     *
  21 C *                                                COMMON/W50512/,     *
  22 C *                                                COMMON/W50513/      *
  23 C *                                                COMMON/W50514/      *
  24 C *                                                COMMON/W505110/     *
  25 C *                                            and COMMON/W505122/     *
  26 C *                                                                    *
  27 C *   Author:  H. Plothow-Besch                                        *
  28 C *            CERN-PPE, CH - 1211 Geneva 23, Switzerland              *
  29 C *                                                                    *
  30 C **********************************************************************
  31 C
  32 #include "pdf/impdp.inc"
  33 C
  34 #include "pdf/w5051p2.inc"
  35 #include "pdf/w5051p7.inc"
  36 #include "pdf/w50510.inc"
  37 #include "pdf/w50511.inc"
  38 #include "pdf/w50512.inc"
  39 #include "pdf/w50513.inc"
  40 #include "pdf/w50514.inc"
  41 #include "pdf/w50515.inc"
  42 #include "pdf/w50517.inc"
  43 #include "pdf/w505110.inc"
  44 #include "pdf/w505122.inc"
  45       DATA ZEROD/0.D0/
  46 C
  47       IFLSTA = 1
  48 C.
  49 C print values
  50       IF(N6.LE.0) N6 = L6
  51       WRITE(N6,2000)
  52       WRITE(N6,2001) NPTYPE,NGROUP,NSET,SFNAME(NPTYPE,NGROUP,NSET),
  53      +               MODECR(NPTYPE,NGROUP,NSET)
  54       WRITE(N6,2002) NFL,LO,TMAS
  55       WRITE(N6,2003) QCDL4,QCDL5
  56       WRITE(N6,2004) XMIN,XMAX
  57       WRITE(N6,2005) Q2MIN,Q2MAX
  58  2000 FORMAT(/,'  PDFLIB : Summary from PDFSTA')
  59  2001 FORMAT(1H ,' Nptype = ',I1,'  Ngroup = ',I2,'  Nset = ',I3,
  60      +         '  Name  = "',A8,'"  CrMode = ',I3)
  61  2002 FORMAT(1H ,' Nfl    = ',I2,',  LO = ',I1,
  62      +           ',  Tmas = ',F7.2,' GeV/c**2')
  63  2003 FORMAT(1H ,' QCDL4  = ',F7.4,' GeV,  QCDL5 = ',F7.4,' GeV')
  64  2004 FORMAT(1H ,' Xmin   = ',1X,E8.2,',    Xmax = ',E11.5)
  65  2005 FORMAT(1H ,' Q2min  =',F7.3,' (GeV/c)**2,  Q2max = ',E8.2,
  66      +          ' (GeV/c)**2')
  67 C.
  68       WRITE(N6,3000)
  69       IF(WXMIN .GT. ZEROD) THEN
  70          PERC = 100.D0*WXMIN/WTXMIN
  71          WRITE(N6,3001) WXMIN,WTXMIN,PERC,XMIN
  72       ENDIF
  73       IF(WXMAX .GT. ZEROD) THEN
  74          PERC = 100.D0*WXMAX/WTXMAX
  75          WRITE(N6,3002) WXMAX,WTXMAX,PERC,XMAX
  76       ENDIF
  77       IF(WQ2MIN .GT. ZEROD) THEN
  78          PERC = 100.D0*WQ2MIN/WTQ2MIN
  79          WRITE(N6,3003) WQ2MIN,WTQ2MIN,PERC,Q2MIN
  80       ENDIF
  81       IF(WQ2MAX .GT. ZEROD) THEN
  82          PERC = 100.D0*WQ2MAX/WTQ2MAX
  83          WRITE(N6,3004) WQ2MAX,WTQ2MAX,PERC,Q2MAX
  84       ENDIF
  85       IF((WXMIN+WXMAX+WQ2MIN+WQ2MAX) .EQ. ZEROD) WRITE(N6,3005)
  86  3000 FORMAT(' ')
  87  3001 FORMAT(1H ,' PDFSTA:    X values ',E9.3,' times out of ',E9.3,
  88      +         ' ( = ',E9.3,' %) smaller allowed minimum of ',E9.3)
  89  3002 FORMAT(1H ,' PDFSTA:    X values ',E9.3,' times out of ',E9.3,
  90      +         ' ( = ',E9.3,' %) bigger  allowed maximum of ',E11.5)
  91  3003 FORMAT(1H ,' PDFSTA: Q**2 values ',E9.3,' times out of ',E9.3,
  92      +         ' ( = ',E9.3,' %) smaller allowed minimum of ',E9.3,
  93      +         ' (GeV/c)**2')
  94  3004 FORMAT(1H ,' PDFSTA: Q**2 values ',E9.3,' times out of ',E9.3,
  95      +         ' ( = ',E9.3,' %) bigger  allowed maximum of ',E9.3,
  96      +         ' (GeV/c)**2')
  97  3005 FORMAT(1H ,' PDFSTA:    NO errors occured')
  98 C.
  99       RETURN
 100       END