*
* $Id: pdfdoc.dat,v 1.1.1.1 1996/03/08 17:48:06 mclareni Exp $
*
* $Log: pdfdoc.dat,v $
* Revision 1.1.1.1  1996/03/08 17:48:06  mclareni
* Pdf
*
*
1***********************************************************************
 *                                                                     *
 *   CERN                                                              *
 *                                                                     *
 *   W5051 PDFLIB                                         1995.03.15   *
 *                                                                     *
 *   ---------------------------------------------------------------   *
 *                                                                     *
 *      PDFLIB: Nucleon, Pion and Photon Parton Density Functions      *
 *                      and alpha(s) Calculations                      *
 *                                                                     *
 *                           User's Manual                             *
 *                                                                     *
 *                           Version 6.06                              *
 *                                                                     *
 *                                                                     *
 *                     H. Plothow-Besch/CERN-PPE                       *
 *                                                                     *
 *                         CH-1211 Geneva 23                           *
 *            BITNET/EARN address: PLOTHOW@CERNVM.CERN.CH              *
 *                    Tel.: +41 - 22 - 767 39 80                       *
 *                                                                     *
 *   A complete manual in PostScript format exists on a separate file  *
 *     Requests can be addressed to the CERN Program Library Office    *
 *                                                                     *
 *                   Copyright by H. Plothow-Besch                     *
 *                                                                     *
 *                            Abstract                                 *
 *                                                                     *
 *       This document describes an integrated package of Parton Den-  *
 *  sity Functions called PDFLIB which has been added to the CERN      *
 *  Program Library Pool W999 and is labelled as W5051. In this pack-  *
 *  age all the different sets of parton density functions of the      *
 *  NUCLEON, the PION and the PHOTON which are available today (about  *
 *  100 different sets) have been put together. All these sets have    *
 *  been combined in a consistent way such that they all have similar  *
 *  calling sequences and no external data files have to be read in    *
 *  anymore. A default set has been prepared, although those prefering *
 *  their own set or wanting to test a new one may do so within the    *
 *  package. The selection mode, the program parameters as well as     *
 *  the possibilities to modify the defaults and to control errors     *
 *  occured during execution are described.                            *
 *                                                                     *
 *       The package offers as well a program to calculate the strong  *
 *  coupling constant alpha(s) to first or second order (default).     *
 *  The correct Lambda(QCD) associated to the selected set of struc-   *
 *  ture functions and the number of allowed flavours with respect to  *
 *  the given Q**2 is automatically used in the caculation. For those  *
 *  structure functions where the evolution has been performed to      *
 *  leading order only, an alpha(s) value fixed to first order is re-  *
 *  turned.                                                            *
 *                                                                     *
 *       The package is written in standard FORTRAN-77 and is avai-    *
 *  lable on all CERN system areas.                                    *
 *                                                                     *
 *        Submitter:   H. Plothow-Besch/CERN-PPE                       *
 *        Language:    FORTRAN 77                                      *
 *        Library:     PDFLIB                                          *
 *                                                                     *
 ***********************************************************************
1***********************************************************************
 *                                                                     *
 *                                                                     *
 *   _______________________________________________________________   *
 *   |                                                             |   *
 *   |  PDFLIB - User's Manual                                     |   *
 *   |                                                             |   *
 *   |  CERN Program Library entry W5051                           |   *
 *   |                                                             |   *
 *   |  Copyright CERN, Geneva 1990, 1991, 1992, 1993, 1994, 1995  |   *
 *   |                                                             |   *
 *   |  Copyright and any other appropriate legal protection of    |   *
 *   |  these computer programs and associated documentation re-   |   *
 *   |  served in all countries of the world.                      |   *
 *   |                                                             |   *
 *   |  These programs or documentation may not be reproduced      |   *
 *   |  and/or redistributed by any method without prior written   |   *
 *   |  consent of the Director-General of CERN or his delegate.   |   *
 *   |                                                             |   *
 *   |  Permission for the scientific usage of any programs des-   |   *
 *   |  cribed herein is granted apriori to those scientific in-   |   *
 *   |  stitutes associated with the CERN experimental program or  |   *
 *   |  with whom CERN has concluded a scientific collaboration    |   *
 *   |  agreement.                                                 |   *
 *   |                                                             |   *
 *   |  Commercial utilisation requires explicit a priori per-     |   *
 *   |  mission from CERN and will be subjected to payment of a    |   *
 *   |  licence fee.                                               |   *
 *   |                                                             |   *
 *   |  These statements are without prejudice to any rights that  |   *
 *   |  third parties may hold in these programs.                  |   *
 *   |                                                             |   *
 *   |                                                             |   *
 *   |  Requests for information should be addressed to:           |   *
 *   |                                                             |   *
 *   |        CERN Program Library Office                          |   *
 *   |        CERN-DD Division                                     |   *
 *   |        CH-1211 Geneva 23                                    |   *
 *   |        Switzerland                                          |   *
 *   |        Tel. +41 22 767 4951                                 |   *
 *   |        Fax. +41 22 767 7155                                 |   *
 *   |        Electronic Mail address:                             |   *
 *   |                                                             |   *
 *   |        EARN/Bitnet:   CERNLIB@CERNVM                        |   *
 *   |        DECnet:        VXCERN::CERNLIB ( node 22.190 )       |   *
 *   |        Internet:      CERNLIB@CERNVM.CERN.CH                |   *
 *   |_____________________________________________________________|   *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 ******************************** - 1 - ********************************
1***********************************************************************
 *                                                                     *
 *      Table of Contents                                              *
 *                                                               page  *
 *                                                                     *
 *      PROGRAM SUMMARY .......................................... 3   *
 *                                                                     *
 *  0.  Major Differences compared to earlier Versions             4   *
 *      0.1  The actual Version 6.06 wrt Version 5.00 ............ 4   *
 *      0.2  Version 5.00 wrt Version 4.00 ....................... 4   *
 *      0.3  Version 4.00 wrt Version 3.00 ....................... 4   *
 *      0.4  Version 3.00 wrt Version 2.00 ....................... 6   *
 *      0.5  Version 2.00 wrt Version 1.00 ....................... 7   *
 *                                                                     *
 *  1.  Introduction                                               8   *
 *                                                                     *
 *  2.  The Structure Functions Package PDFLIB                    10   *
 *      2.1  Parameter Setting and Selection of sets ............ 10   *
 *      2.2  The Calling Sequence: PDFLIB Format ................ 12   *
 *      2.3  The Calling Sequence: PDG Format ................... 14   *
 *      2.4  Error Handling ..................................... 14   *
 *      2.5  Features ........................................... 15   *
 *      2.6  PION Structure Functions ........................... 16   *
 *      2.7  PHOTON Structure Functions ......................... 17   *
 *      2.8  How to find PDFLIB ................................. 18   *
 *                                                                     *
 *  3.  alpha(s) Calculation                                      18   *
 *                                                                     *
 *  4.  Conclusions                                               19   *
 *                                                                     *
 *  5.  Acknowledgements                                          21   *
 *                                                                     *
 *  6.  Appendix                                                  22   *
 *      A : Example of how to use PDFLIB ........................ 22   *
 *      B : List of subprograms and COMMON blocks in PDFLIB            *
 *          which are relevant to the user ...................... 25   *
 *      C : Complete list of subprograms in PDFLIB .............. 26   *
 *                                                                     *
 *  7.  References                                                32   *
 *                                                                     *
 *  8.  List of Tables                                            34   *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 ******************************** - 2 - ********************************
1***********************************************************************
 *                                                                     *
 *  PROGRAM SUMMARY                                                    *
 *                                                                     *
 *                                                                     *
 *  Titel of Program : PDFLIB, version 6.06                            *
 *                                                                     *
 *  Catalogue number : W5051 in the CERN Computer Program Library      *
 *                                                                     *
 *  Program obtainable from : CERN Computer Program Library Office     *
 *                                                                     *
 *  Computer for which the program is designed and others on which it  *
 *  has been tested:                                                   *
 *  Computer : IBM 3090, VAX, IBM-RISC, DECS, SUN, APOLLO, HP-UX,      *
 *             Silicon Graphics, NEXT, CDC, CRAY, IBM-PC               *
 *                                                                     *
 *  Operating systems : VM/CMS, VAX/VMS, ULTRIX                        *
 *                                                                     *
 *  Programming language used : FORTRAN 77, program is available also  *
 *                              in the PATCHY and the CMZ format       *
 *                                                                     *
 *  Number of lines in distributed progam,                             *
 *  including test data, etc : 99 000                                  *
 *                                                                     *
 *  Memory required,                                                   *
 *  including testprograms and some documentation : 5 700 000 bytes    *
 *                                                                     *
 *  No. of bits in a word : 64                                         *
 *                                                                     *
 *  Peripherals used : line printer (optional)                         *
 *                                                                     *
 *  Keywords : parton density functions, structure functions, W5051,   *
 *             nucleon-, proton-, pion-, photon density parametrisa-   *
 *             tions, Altarelli-Parisi equation, strong coupling       *
 *             constant alpha(s)                                       *
 *                                                                     *
 *  Nature of physical problem :                                       *
 *  All theoretical calculations of cross sections involving the       *
 *  hadronic structure of the incoming particles, i.e. protons, pions  *
 *  or photons, use the parton density at a given longitudinal momen-  *
 *  tum X of the parton and at a given momentum transfer Q**2. There   *
 *  is a variety of parametrisations available, but each of them has   *
 *  its own structure for input and/or output.                         *
 *                                                                     *
 *  Method of solution :                                               *
 *  Collection of all available parton density parametrisations and    *
 *  application in an identical format for all of them, regardless of  *
 *  the author's convention or of the particle type.                   *
 *                                                                     *
 *  References :                                                       *
 *  [1] H. Plothow-Besch, 'PDFLIB: Structure Functions and alpha(s)    *
 *      Calculations', User's Manual - Version 1.00, W5051 PDFLIB,     *
 *      1991.03.21, CERN-PPE.                                          *
 *  [2] H. Plothow-Besch, 'PDFLIB: a library of all available parton   *
 *      density functions of the nucleon, the pion and the photon and  *
 *      the corresponding alpha(s) calculations',                      *
 *      Comp. Phys. Comm. 75 (1993) 396-416.                           *
 *                                                                     *
 *                                                                     *
 ******************************** - 3 - ********************************
1***********************************************************************
 *                                                                     *
 *  0.  Major Differences compared to earlier Versions                 *
 *                                                                     *
 *  0.1  The actual Version 6.06 wrt Version 5.00:                     *
 *                                                                     *
 *     The following major changes compared to version 5.00 of PDFLIB  *
 *  have been implemented:                                             *
 *                                                                     *
 *  Several NUCLEON structure function sets have been added:           *
 *    - the NLL sets A' and G of Martin, Roberts and Stirling (MRS-Ap) *
 *      and (MRS-G) of February 1995 with grid data as well as in its  *
 *      parametrized version (MRS-Ap-F) and (MRS-G-Fit),               *
 *    - the version 3 of the CTEQ group of October 1994 as NLL parame- *
 *      trisation in the MS_bar as well as in the DIS renormalisation  *
 *      scheme and as LO (CTEQ3M, CTEQ3D, CTEQ3L),                     *
 *    - the new sets of GRV of November 1994 as NLL in the MS_bar as   *
 *      well as in the DIS renormalisation scheme and as LO GRV        *
 *      (GRV94-MS, GRV94DI, GRV94-LO),                                 *
 *                                                                     *
 *  a NEW group of PHOTON structure function sets have been added:     *
 *    - the six NLL sets of the WHIT group (WHIT1-G to WHIT6-G).       *
 *                                                                     *
 *  Please note that as new the default set the MRS set (G) has been   *
 *  chosen (NPTYPE = 1, NGROUP = 3, NSET = 41).                        *
 *                                                                     *
 *     More details can be found on the next pages of this updated     *
 *  User's Manual.                                                     *
 *                                                                     *
 *     The response from the Physics community has been very satis-    *
 *  factory, which is encouraging us to keep up-to-date with the       *
 *  latest developments in this area. The author will be pleased to    *
 *  learn about new parton density functions, as well as to receive    *
 *  suggestions to improve both, the usage and the documentation.      *
 *                                                                     *
 *                                                                     *
 *  0.2  Version 5.00 wrt Version 4.00:                                *
 *                                                                     *
 *     The following major changes compared to version 4.00 of PDFLIB  *
 *  have been implemented:                                             *
 *                                                                     *
 *  NUCLEON structure function sets have been added:                   *
 *    - the NLL set A of Martin, Roberts and Stirling (MRS-A) and the  *
 *      parametrized version (MRS-A-Fit)                               *
 *  PHOTON structure function sets have been added:                    *
 *    - the NLL set of Aurenche, Fontannaz and Guillet (AFG-G).        *
 *                                                                     *
 *                                                                     *
 *  0.3  Version 4.00 wrt Version 3.00:                                *
 *                                                                     *
 *     The following major changes compared to version 3.00 of PDFLIB  *
 *  have been implemented:                                             *
 *                                                                     *
 *  NUCLEON structure function sets have been updated or added:        *
 *    - the updated NLL sets of Martin, Roberts and Stirling (MRS) of  *
 *      November 1992, the sets S0', D0' and D-' in the MSbar as well  *
 *      as the sets S0', D0' and D-' in the DIS renormalisation        *
 *      scheme. Please note that the previous S0, D0 and D- sets have  *
 *                                                                     *
 ******************************** - 4 - ********************************
1***********************************************************************
 *                                                                     *
 *      been retracted by the authors. The new MRS(H) set in both      *
 *      schemes have also been added,                                  *
 *    - the NLL and LO sets of the CTEQ collaboration (CTEQ) in its    *
 *      parametrised form of January 1993. These are the sets 1M, 1MS  *
 *      and 1ML in the MSbar, the set 1D in the DIS renormalisation    *
 *      scheme and the LO set 1L,                                      *
 *    - the new NLL sets A and B of Berger and Meng (BM) of            *
 *      February 1993 with special gluon distributions, and finally    *
 *    - the new NLL and LO sets of the CTEQ collaboration (CTEQ) in    *
 *      its parametrised form of August 1993. These are the sets 2M,   *
 *      2MS, 2MF and 2ML in the MSbar, the set 2D in the DIS renorma-  *
 *      lisation scheme and the LO set 2L. The revised CTEQ2 version,  *
 *      CTEQ2p, of November 1993 is recommended by the authors.        *
 *                                                                     *
 *  Please note that the default set has been put to the GRV set HO    *
 *  (NPTYPE = 1, NGROUP = 5, NSET = 3 / former MODE = 72).             *
 *  The MOST DRASTICAL change of the package has been, however, the    *
 *  NEW FORMAT of PARAMETER SETTINGS. Instead of using only ONE para-  *
 *  meter (MODE) to select a set of parton density functions (PDF),    *
 *  each PDF set is now identified by THREE parameters. These are      *
 *  NPTYPE, NGROUP and NSET, for PARTICLE TYPE (Nucleons, Pions or     *
 *  Photons) AUTHOR GROUP (f.ex. MRS, CTEQ, GRV, etc) and the PARAME-  *
 *  TRISATION SET within the group. Therefore the user should provide  *
 *  THREE parameters (NPTYPE, NGROUP, NSET) to the subroutine PDFSET,  *
 *  each time a PDF set shall be selected. But to make life nice and   *
 *  easy, it is also possible to select a PDF set by transferring only *
 *  ONE parameter to PDFSET as before: this parameter should then      *
 *  identify the NAME of the author group in an unequivocal way; the   *
 *  selection of the PDF set within the author group is made by the    *
 *  corresponding VALUE value (f.ex. PARM(1) = 'MRS' with VALUE(1) =   *
 *  29.0D0 would select the 29th set of the Nucleon PDFs of MRS, and   *
 *  f.ex. PARM(1) = 'GRVph' with VALUE(1) = 1.0D0 would select the     *
 *  first set of the Photon PDFs of GRV). As a suggestion, the NAME for*
 *  each PDF set can be found in the NEW COMMON block /W505110/SFNAME, *
 *  where SFNAME is a CHARACTER*8 array with dimension SFNAME(NPTYMX,  *
 *  NGRMAX,NSETMX), and NPTYMX, NGRMAX, NSETMX are defined by the      *
 *  parameter sequence, +SEQ, W5051P2. In case the user wants to se-   *
 *  lect the DEFAULT, it is sufficient to transfer as the ONLY parame- *
 *  ter PARM(1) = 'Nset' with VALUE(1) = 0.0D0. This change of format  *
 *  had been necessary to take into account the inflationary growing   *
 *  number of PDFs until now and to foresee already the possibility of *
 *  adding new PDFs in a more flexible way. Please note that in case   *
 *  of any INCONSISTANCY of parameter settings, it is always chosen    *
 *  the DEFAULT settings !
 *  Please also note that the internal COMMON block /W50511/ has been  *
 *  modified to /W50511/ NPTYPE,NGROUP,NSET,MODE,NFL,LO,TMAS.          *
 *                                                                     *
 *  Full backward compatibility with previous versions is ensured !!   *
 *  This means that old programs using PDF sets up to version 3.00 of  *
 *  PDFLIB should run without ANY change. To ensure this, four new     *
 *  COMMON blocks have been added:                                     *
 *    - /W50519/ NEWVER :  to flag the NEW/OLD version with the        *
 *      LOGICAL variable NEWVER,                                       *
 *    - /W505120/ NPGSMX(NPTYMX,NGRMAX),NSETFL(NPTYMX,NGRMAX,NSETMX) : *
 *      to provide the actual maximum number of parametrisation sets   *
 *      per particle type and author group, and to flag the retracted  *
 *                                                                     *
 ******************************** - 5 - ********************************
1***********************************************************************
 *                                                                     *
 *      PDF sets (NSETFL(I,J,K) = 0),                                  *
 *    - /W505121/ NPTYCR(MODEMX),NGROMX(MODEMX),NSETCR(MODEMAX) :      *
 *      to provide the cross reference to the NPTYPE, NGROUP and NSET  *
 *      value (NEW format) for a given MODE value (OLD format),        *
 *    and                                                              *
 *    - /W505122/ MODECR(NPTYMX,NGRMAX,NSETMX) :  to provide the       *
 *      cross reference to the MODE value (OLD version) for a given    *
 *      NPTYPE, NGROUP and NSET value (NEW format). In case of the     *
 *      above added PDF sets NO MODE value is available; in that case  *
 *      MODECR(I,J,K) = -1 is returned.                                *
 *                                                                     *
 *                                                                     *
 *  0.4  Version 3.00 wrt Version 2.00                                 *
 *                                                                     *
 *     The following major changes compared to version 2.00 of PDFLIB  *
 *  had been implemented:                                              *
 *                                                                     *
 *  NUCLEON structure function sets had been updated or added:         *
 *    - the NLL set of Aurenche et al. (ABFOW),                        *
 *    - the updated NLL set of Kwiecinski, Martin, Roberts and         *
 *      Stirling (KMRS) set B0 with u_bar not equal to d_bar           *
 *      of April 1992,                                                 *
 *    - the new NLL sets of Martin, Roberts and Stirling (MRS)         *
 *      set S0, D0 and D- with u_bar not equal to d_bar of April 1992. *
 *  Please note that the default set had been put to the MRS set S0    *
 *  (MODE = 45).                                                       *
 *                                                                     *
 *  PION structure function sets had been added:                       *
 *    - the new LO and NLL sets of Gl"uck, Reya and Vogt (GRV-P)       *
 *      of August 1991,                                                *
 *    - the new sets 2 and 3 of Aurenche et al. (ABFKW-P).             *
 *  Please also note that the ABFKW set had changed MODE number from   *
 *  110 -> 180, as well as the SMRS sets from 120ff -> 130ff, and      *
 *  that earlier inconsistencies in the parton density definitions     *
 *  had been corrected.                                                *
 *                                                                     *
 *  For the first time also PHOTON structure function sets had been    *
 *  added to the package. These sets were the following:               *
 *    - the old Duke and Owens (DO-G) LO and NLL sets,                 *
 *    - the old LO sets 1 to 4 of Drees and Grassie (DG-G),            *
 *    - the new LO sets 1 to 3 of Abramowicz, Charchula and Levy       *
 *      (LAC-G),                                                       *
 *    - the new LO sets 1 and 2 and the NLL set of Gordon and Storrow  *
 *      (GS-G) of July 1991,                                           *
 *    - the new LO and NLL sets, as well as the leading term of the    *
 *      NLL calculation of Gl"uck, Reya and Vogt (GRV-G)               *
 *      of December 1991,                                              *
 *    - the new NLL sets of Aurenche et al. (ACFGP-G) of April 1992.   *
 *                                                                     *
 *  Some improvements in the usage of the package had been made:       *
 *    - more protections had been implemented and inconsistencies      *
 *      had been corrected to make the package more robust and         *
 *      available on ALL system areas, and finally                     *
 *    - the subroutine STRUCTM had been added. STRUCTM is basically    *
 *      the same routine as the old steering subroutine STRUCTF apart  *
 *      from ONE additional parameter in the argument list which had   *
 *                                                                     *
 ******************************** - 6 - ********************************
1***********************************************************************
 *                                                                     *
 *      been added to transfer separately the parton densities for     *
 *      u_bar and d_bar. The argument list of the subroutine STRUCTF   *
 *      was as before; in case of STRUCTF sea = (u_bar + d_bar)/2 is   *
 *      transferred. The subroutine PFTOPDG, which is the sterring     *
 *      routine in PDG format, was unchanged; this routine always      *
 *      transferred u_bar and d_bar separately.                        *
 *                                                                     *
 *                                                                     *
 *  0.5  Version 2.00 wrt Version 1.00                                 *
 *                                                                     *
 *     The following major changes compared to version 1.00 of PDFLIB  *
 *  had been implemented:                                              *
 *                                                                     *
 *  New parton density functions of the NUCLEON had been added:        *
 *    - the old LO set of Buras and Gaemers of 1978 (BEBC) for com-    *
 *      pleteness,                                                     *
 *    - all the LO sets of Diemoz, Ferroni, Longo and Martinelli (DFLM)*
 *      including their parametrisation for the NLL average fit,       *
 *    - the new LO set 1.1 of Owens (DO 1.1) of June 1991,             *
 *    and finally                                                      *
 *    - the new LO and NLL sets of Gl"uck, Reya and Vogt (GRV) of      *
 *      May 1991.                                                      *
 *                                                                     *
 *  PION structure function sets had been added to the package:        *
 *    - the old Owens (OW-P) sets 1 and 2,                             *
 *    - the newer Aurenche et al. (ABFKW-P) set and                    *
 *    - the new Sutton, Martin, Roberts and Stirling (SMRS-P) sets 1,  *
 *      2 and 3.                                                       *
 *                                                                     *
 *  Some improvements in the usage of the package had been made:       *
 *    - COMMON block names had been changed to avoid interference      *
 *      to those defined by the user (/IMODE/ -> /W50511/,             *
 *      /ALAMB/ -> /W50512/),                                          *
 *    - three new COMMON blocks had been added:                        *
 *      a) /W50513/ XMIN,XMAX,Q2MIN,Q2MAX :                            *
 *         to provide the minimum and maximum values of X and Q**2 for *
 *         each set of structure functions,                            *
 *      b) /W50510/ IFLPRT :                                           *
 *         to provide a print flag at different levels to print        *
 *         b.1) the variables of the three COMMON blocks               *
 *              /W50511/, /W50512/ and /W50513/ at job initialisation  *
 *         b.2) error messages during job execution                    *
 *      c) /W50514W/ PDFWGT :                                          *
 *         to take into account the event weight for error counting    *
 *    - a subroutine PDFSTA had been added to allow at job termination *
 *      for a summary how often X and Q**2 limits have been exceeded,  *
 *    - the MODE number of some structure functions had been changed   *
 *      to keep consistensy,                                           *
 *    - the number of parameters which can be set through the sub-     *
 *      routine PDFSET had been extended to 20,                        *
 *    - more protections had been implemented and inconsistencies      *
 *      had been corrected to make the package more robust,            *
 *    and finally                                                      *
 *    - a subroutine PFTOPDG had been provided as an interface for     *
 *      those users who prefer the flavour code convention of the PDG  *
 *      instead of the PDFLIB format.                                  *
 *                                                                     *
 ******************************** - 7 - ********************************
1***********************************************************************
 *                                                                     *
 *  1.  Introduction                                                   *
 *                                                                     *
 *                                                                     *
 *     The number of useful sets of parton density functions of the    *
 *  NUCLEON available today is about 70 and is still growing. The ol-  *
 *  der sets are leading order evolutions only, like                   *
 *                                                                     *
 *   -  Buras and Gaemers (BEBC) [1],                                  *
 *   -  Owens and Reya (OR) [2],                                       *
 *   -  Baier, Engels and Petersson (BEP) [3],                         *
 *   -  Gl"uck, Hoffmann and Reya (GHR) [4],                           *
 *   -  Duke and Owens (DO) sets 1 and 2 [5],                          *
 *   -  Eichten, Hinchliffe, Lane and Quigg (EHLQ) sets 1 and 2 [6],   *
 *   -  Diemoz, Ferroni, Longo and Martinelli (DFLM) with              *
 *      Lambda_4(QCD) = 200 MeV for soft and hard valence quark and    *
 *      soft and hard gluon distributions [7], and finally             *
 *   -  the new LO paramerisation of Owens set 1.1 (DO 1.1) [8].       *
 *                                                                     *
 *     Most of the sets, essentially the more recent ones, are         *
 *  next-to-leading order evolutions, like                             *
 *                                                                     *
 *   -  Diemoz, Ferroni, Longo and Martinelli (DFLM) with              *
 *      Lambda_4(QCD) = 160, 260, 360 MeV [7],                         *
 *   -  Martin, Roberts and Stirling (MRS) sets 1 to 3, sets E and B   *
 *      and sets E' and B' [9],                                        *
 *   -  Harriman, Martin, Roberts and Stirling (HMRS) sets E, E+, E-   *
 *      and B [10],                                                    *
 *   -  Kwiecinski, Martin, Roberts and Stirling (KMRS) [11] sets B0   *
 *      and B-, and for the B- set with different radii for shadowing, *
 *   -  Martin, Roberts and Stirling (MRS-B Lambda) [11] for the B0    *
 *      set also with different values of Lambda_4(QCD) = 135, 160,    *
 *      200 and 235 MeV,                                               *
 *   -  Martin, Roberts and Stirling (MRS) [12] sets S0, D0 and D-     *
 *      as well as the upgraded KMRS set B0, all having u_bar not      *
 *      equal to d_bar,                                                *
 *   -  Morfin and Tung (MT) sets 1 to 6 [13],                         *
 *   -  Gl"uck, Reya and Vogt (GRV) [14] including their new para-     *
 *      metrisation from August 1991,                                  *
 *   -  Aurenche et al. (ABFOW) [15],                                  *
 *   -  the upgraded Martin, Roberts and Stirling (MRS) [16] sets      *
 *      S0', D0' and D-' in the MSbar renormalisation scheme and MRS(H)*
 *   -  the upgraded Martin, Roberts and Stirling (MRS) [17] sets      *
 *      S0', D0' and D-' in the DIS renormalisation scheme and MRS(H), *
 *    - the sets 1M, 1MS and 1ML in the MSbar and the set 1D in the    *
 *      DIS renormalisation scheme of the CTEQ collaboration (CTEQ)    *
 *      [18] in its parametrised form,                                 *
 *    - the new sets A and B of Berger and Meng (BM) [19] with special *
 *      gluon distributions, and                                       *
 *    - the new sets 2M, 2MS, 2MF and 2ML in the MSbar and the set 2D  *
 *      in the DIS renormalisation scheme of the CTEQ collaboration    *
 *      (CTEQ) [20] in its parametrised form which superseed already   *
 *      the previous CTEQ sets 1M, 1MS, 1ML and 1D. The CTEQ collabo-  *
 *      ration provides also parametrisations to leading order, the    *
 *      sets 1L [18] and 2L [20].                                      *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 ******************************** - 8 - ********************************
1***********************************************************************
 *                                                                     *
 *     There is also a limited number of parton density functions of   *
 *  the PION available. These are the old set to leading order evolu-  *
 *  tion of                                                            *
 *                                                                     *
 *   -  Owens (OW-P) [21], sets 1 and 2,                               *
 *                                                                     *
 *  and the newer and the very new next-to-leading order evolutions of *
 *                                                                     *
 *   -  Aurenche et al. (ABFKW-P) [22], sets 1, 2 and 3,               *
 *   -  Sutton, Martin, Roberts and Stirling (SMRS-P) [23], sets 1, 2  *
 *      and 3 and                                                      *
 *   -  Gl"uck, Reya and Vogt (GRV-P) [24], the LO and NLL sets.       *
 *                                                                     *
 *     In addition there are also a few parton density functions of    *
 *  the PHOTON available. These are the old sets to leading order evo- *
 *  lution of                                                          *
 *                                                                     *
 *   -  Duke and Owens (DO-G) (and the NLL) [25] using an asymptotic   *
 *      solution of the Altarelli-Parisi equation,                     *
 *   -  Drees and Grassie (DG-G) sets 1 to 4 [26] using a full solu-   *
 *      tion of the Altarelli-Parisi equation,                         *
 *                                                                     *
 *  the new leading order evolution of                                 *
 *                                                                     *
 *   -  Abramowicz, Charchula and Levy (LAC-G) [27] sets 1 to 3 using  *
 *      a full solution of the Altarelli-Parisi equation,              *
 *                                                                     *
 *  and the new next-to-leading order evolution of                     *
 *                                                                     *
 *   -  Gordon and Storrow (GS-G) [28], the LO sets 1 and 2 and the    *
 *      NLL set,                                                       *
 *   -  Gl"uck, Reya and Vogt (GRV-G) [29], the LO, NLL and the        *
 *      leading term of the NLL sets, and                              *
 *   -  Aurenche et al. (ACFGP-G) [30] sets 1 and 2, without and with  *
 *      massive charm.                                                 *
 *                                                                     *
 *     The different sets differ also in the renormalisation scheme    *
 *  they use. The DFLM series, the sets 1 to 6 of the MT series, the   *
 *  new MRS sets S0', D0' and D-' and the new CTEQ sets 1D and 2D are  *
 *  performed in the DIS renormalisation scheme [31]. If the LO evo-   *
 *  lutions of these authors are ignored, all the other sets are       *
 *  performed in the MS_bar renormalisation scheme. All sets are made  *
 *  using four flavours in the initial state (NF = 4). The value of    *
 *  the QCD scale factor, Lambda_4(QCD), ranges from 45 to 500 MeV.    *
 *                                                                     *
 *     We have put together all these different sets of parton density *
 *  functions in one single package [32]. We have modified the sets    *
 *  such that no external data files for the grids have to be read in  *
 *  anymore. All these structure function sets have been combined in a *
 *  consistent way such that they all appear in an IDENTICAL STRUCTURE *
 *  to the user. The selection is made via THREE parameters, NPTYPE,   *
 *  NGROUP and NSET, which identify a PDF set by its PARTICLE TYPE     *
 *  (Nucleon: NPTYPE = 1, Pion: NPTYPE = 2, Photon: NPTYPE = 3), its   *
 *  AUTHOR GROUP (f.ex. MRS: NGROUP = 3, CTEQ: NGROUP = 4,             *
 *  GRV: NGROUP = 5) and its PDF SET WITHIN THE GROUP (NSET = 1,2,3    *
 *  etc.). It is also possible to select a PDF set by transferring     *
 *                                                                     *
 ******************************** - 9 - ********************************
1***********************************************************************
 *                                                                     *
 *  only ONE parameter to PDFSET as before: this parameter should then *
 *  identify the NAME of the author group in an unequivocal way; the   *
 *  selection of the PDF set within the author group is made by the    *
 *  corresponding VALUE value (f.ex. PARM(1) = 'MRS' with VALUE(1) =   *
 *  29.0D0 would select the 29th set of the Nucleon PDFs of MRS, and   *
 *  f.ex. PARM(1) = 'GRVph' with VALUE(1) = 1.0D0 would select the     *
 *  first set of the Photon PDFs of GRV). As a suggestion, the NAME for*
 *  each PDF set can be found in the new COMMON block /W505110/SFNAME, *
 *  where SFNAME is a CHARACTER*8 array with dimension SFNAME(NPTYMX,  *
 *  NGRMAX,NSETMX), and NPTYMX, NGRMAX, NSETMX are defined by the      *
 *  parameter sequence, +SEQ, W5051P2. The parameters should be,       *
 *  NPTYPE, NGROUP, NSET or the 'NAME' of the PDF set with their VALUE *
 *  values should be set with a call to the subroutine PDFSET at the   *
 *  initialisation phase. A simple SUBROUTINE call, which is IDENTICAL *
 *  for all applications (nucleons, pions, photons) returns the parton *
 *  densities for all partons (u, d, s, c, b, t, gluon and their anti- *
 *  quarks) at a given X value (where X is the fraction of the longi-  *
 *  tudinal momentum carried by the parton) and the Q-scale SCALE (in  *
 *  GeV). A default set has been prepared, in which case it is suffi-  *
 *  cient to transfer as the ONLY parameter PARM(1) = 'Nset' with      *
 *  VALUE(1) = 0.0D0, although those preferring their own private set  *
 *  or wanting to test a new one may do so within the package. Error   *
 *  control can be obtained by setting a print flag to obtain output   *
 *  either during execution and/or as a summary at job termination via *
 *  a call to the subroutine PDFSTA. In addition, the different parton *
 *  density sets have been modified such that the library can be used  *
 *  on all the different computer systems known today. The library has *
 *  been tested on the different system areas to ensure identical re-  *
 *  sults within the given machine precision.                          *
 *                                                                     *
 *     Within the same package a program is also provided to calculate *
 *  the strong coupling constant alpha(s) to second order (default)    *
 *  or to first order (by user's choice). The correct Lambda(QCD) as-  *
 *  sociated to the selected set of structure functions and the number *
 *  of allowed flavours with respect to the given Q**2 is automati-    *
 *  cally used in the calculation. For those structure functions where *
 *  the evolution has been performed to leading order only, the        *
 *  alpha(s) value to first order is returned.                         *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *  2.  The Structure Functions Package PDFLIB                         *
 *                                                                     *
 *                                                                     *
 *     In the following the use of and the access to the PDFLIB pack-  *
 *  age is described.                                                  *
 *                                                                     *
 *  2.1  Parameter Setting and Selection of Sets                       *
 *                                                                     *
 *                                                                     *
 *     The user may access his/her preferred set of structure func-    *
 *  tions by setting the THREE parameters, NPTYPE, NGROUP and NSET,    *
 *  which identify each set of parton density functions via a call     *
 *  to the subroutine                                                  *
 *                                                                     *
 *       CALL PDFSET(PARM,VALUE)                                       *
 *                                                                     *
 ******************************* - 10 - ********************************
1***********************************************************************
 *                                                                     *
 *  at the initialization phase of his/her MAIN program, or use the    *
 *  package as a 'black box' with the default values (see below). The  *
 *  arguments PARM and VALUE are vectors of dimension 20 and have the  *
 *  following meaning                                                  *
 *                                                                     *
 *       PARM(I)  = character*20 variable, which defines in any order  *
 *                  the variables                                      *
 *                  'NPTYPE', 'NGROUP', 'NSET', 'MODE', 'INIT0',       *
 *                  'NFL', 'LO', 'TMAS', 'QCDL4', 'QCDL5' and          *
 *                  'XMIN', 'XMAX', 'Q2MIN', 'Q2MAX'                   *
 *       VALUE(I) = the corresponding numerical value of the variable  *
 *                  PARM(I)                                            *
 *                  (TMAS, QCDL4, QCDL5, XMIN, XMAX, Q2MIN and Q2MAX   *
 *                  are DOUBLE PRECISION variables),                   *
 *                                                                     *
 *  where                                                              *
 *                                                                     *
 *       NPTYPE = number of particle type ranging from 1 to 3          *
 *                (Nucleons: NPTYPE = 1, Pions: NPTYPE = 2 and         *
 *                 Photons: NPTYPE = 3) (Default: NPTYPE = 1)          *
 *       NGROUP = number of author group ranging from 1 to 7           *
 *                (Default: NGROUP = 5)                                *
 *       NSET   = number of a selected structure function set within   *
 *                the author group ranging from 1 to 34                *
 *                (Default: NSET = 3)                                  *
 *      (MODE   = number of a selected structure function set ranging  *
 *                from 0 to 281 - OLD format !!)                       *
 *                (Default: MODE = 45)                                 *
 *       INIT0  = in case of PARM(1) = 'INIT0' PDFSET fills as the     *
 *                only action the COMMON blocks /W505110/, /W505120/   *
 *                /W505121/ and /W505122/                              *
 *       NFL    = desired number of flavours in the alpha(s) calcula-  *
 *                tion ranging from 3 to 6                             *
 *                (Default: NFL = 5)                                   *
 *       LO     = order of alpha(s) calculation; if LO = 1, alpha(s)   *
 *                to first order only                                  *
 *                (Default: LO = 2)                                    *
 *       TMAS   = the user defined value of the top-quark mass in      *
 *                GeV/c**2 (optional)                                  *
 *                (Default: TMAS = 100.0D0)                            *
 *       QCDL4  = QCD scale, Lambda_4(QCD), in GeV for four flavours   *
 *       QCDL5  = QCD scale, Lambda_5(QCD), in GeV for five flavours   *
 *                corresponding to QCDL4                               *
 *                                                                     *
 *  and                                                                *
 *                                                                     *
 *       XMIN   = minimum allowed X value                              *
 *       XMAX   = maximum allowed X value                              *
 *       Q2MIN  = minimum allowed Q**2 value (in (GeV/c)**2)           *
 *       Q2MAX  = maximum allowed Q**2 value (in (GeV/c)**2)           *
 *                                                                     *
 *  for each set of structure functions.                               *
 *                                                                     *
 *     If the user wants to select the default settings, it is suffi-  *
 *  cient to give PARM(1) = 'Nset' with VALUE(1) = 0.0D0 as the ONLY   *
 *  parameter to PDFSET. Please note that a PDF set may also be iden-  *
 *                                                                     *
 ******************************* - 11 - ********************************
1***********************************************************************
 *                                                                     *
 *  tified by the NAME of the author group provided in an unequivocal  *
 *  way. The selection of the PDF set within the author group is, in   *
 *  that case, given by the corresponding VALUE value. The NAME for    *
 *  each PDF set can be found in the COMMON block /W505110/SFNAME,     *
 *  where SFNAME is a CHARACTER*8 array with dimension SFNAME(NPTYMX,  *
 *  NGRMAX,NSETMX), and NPTYMX, NGRMAX, NSETMX are defined by the para-*
 *  meter sequence, +SEQ, W5051P2. Please also note that the variable  *
 *  names can be transferred in lower, upper or mixed character modes. *
 *                                                                     *
 *     Please also note that PDFSET can be called as often as the user *
 *  likes. In order to redefine the parameters to select other sets of *
 *  structure functions which the user wants to investigate, it might  *
 *  be necessary to call PDFSET in an alternating way, but it is al-   *
 *  ways mandatory to transfer either the THREE parameters, NPTYPE,    *
 *  NGROUP, NSET, or the 'NAME' of the author group with their corres- *
 *  ponding VALUE values. The subroutine PDFSET fills the internal     *
 *  COMMON blocks                                                      *
 *                                                                     *
 *       COMMON/W50511/ NPTYPE,NGROUP,NSET,MODE,NFL,LO,TMAS            *
 *       COMMON/W50512/ QCDL4,QCDL5                                    *
 *       COMMON/W50513/ XMIN,XMAX,Q2MIN,Q2MAX                          *
 *                                                                     *
 *  at the time the routine is called. Please note that in case of a   *
 *  multiple call to PDFSET with fewer parameters redefined than in a  *
 *  preceeding call, always the last parameters are kept in memory.    *
 *  All variables of the three COMMON blocks can be automatically      *
 *  printed at job initialisation by setting the print flag IFLPRT in  *
 *  the COMMON block                                                   *
 *                                                                     *
 *       COMMON/W50510/ IFLPRT                                         *
 *                                                                     *
 *  to IFLPRT = 2.                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *  2.2  The Calling Sequence: PDFLIB Format                           *
 *                                                                     *
 *                                                                     *
 *     The main steering routine for a set of structure functions is   *
 *  accessed as follows:                                               *
 *                                                                     *
 *       CALL STRUCTM(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)    *
 *                                                                     *
 *     Please note that all variables are defined as DOUBLE PRECISION. *
 *  The user has to provide the following INPUTs:                      *
 *                                                                     *
 *       X      = X value of parton                                    *
 *       SCALE  = QCD scale in GeV                                     *
 *                                                                     *
 *     The subroutine STRUCTM returns the following OUTPUT:            *
 *                                                                     *
 *       UPV    = up valence quark                                     *
 *       DNV    = down valence quark                                   *
 *       USEA   = sea (up_bar)                                         *
 *       DSEA   = sea (down_bar)                                       *
 *       STR    = strange quark                                        *
 *                                                                     *
 ******************************* - 12 - ********************************
1***********************************************************************
 *                                                                     *
 *       CHM    = charm quark                                          *
 *       BOT    = bottom quark                                         *
 *       TOP    = top quark                                            *
 *       GL     = gluon                                                *
 *                                                                     *
 *     This steering routine has separate arguments for u_bar and      *
 *  d_bar-quarks to take into account that some of the new sets of     *
 *  nucleon density functions provide separately the densities for     *
 *  u_bar and d_bar. In case u_bar is not given separately from d_bar  *
 *  it is set USEA = DSEA.                                             *
 *                                                                     *
 *  Please note that STRUCTM returns  X * parton distribution          *
 *  function !                                                         *
 *                                                                     *
 *     An older version of the main steering routine
 *                                                                     *
 *       CALL STRUCTF(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,TOP,GL)          *
 *                                                                     *
 *  is still available. Also here all variables are defined as         *
 *  DOUBLE PRECISION. Please note that STRUCTF has only ONE argument   *
 *  for the sea-quark, which gives u_bar = d_bar or (u_bar + d_bar)/2. *
 *  The user has to provide as above the following INPUTs:             *
 *                                                                     *
 *       X      = X value of parton                                    *
 *       SCALE  = QCD scale in GeV                                     *
 *                                                                     *
 *     The subroutine STRUCTF returns the following OUTPUT:            *
 *                                                                     *
 *       UPV    = up valence quark                                     *
 *       DNV    = down valence quark                                   *
 *       SEA    = sea (up_bar = down_bar or (up_bar + down_bar)/2)     *
 *       STR    = strange quark                                        *
 *       CHM    = charm quark                                          *
 *       BOT    = bottom quark                                         *
 *       TOP    = top quark                                            *
 *       GL     = gluon                                                *
 *                                                                     *
 *     Please note that also STRUCTF returns  X * parton distribution  *
 *  function !                                                         *
 *                                                                     *
 *     A list of available sets of NUCLEON structure functions in the  *
 *  package is given in Table 1.                                       *
 *                                                                     *
 *     The recommended set of structure functions from the different   *
 *  authors is indicated with a flash. If NSET is set to zero          *
 *  (or if any of the three parameters, NPTYPE, NGROUP and NSET, is    *
 *  undefined) a default set is selected which is a nucleon parton     *
 *  dennsity set, the set of MRS (G) with Lambda_4(QCD) = 255 MeV      *
 *  (NPTYPE = 1, NGROUP = 3, NSET = 41).                               *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 ******************************* - 13 - ********************************
1***********************************************************************
 *                                                                     *
 *  2.3  The Calling Sequence: PDG Format                              *
 *                                                                     *
 *                                                                     *
 *     To allow the use of the package with the flavour code conven-   *
 *  tion of the Particle Data Group (PDG) an interface has been writ-  *
 *  ten which translates the PDFLIB format (see Section 2.2) into the  *
 *  PDG format. Instead of calling the subroutine STRUCTM the user     *
 *  preferring the flavour code convention of the PDG accesses the     *
 *  package via a call to the subroutine                               *
 *                                                                     *
 *       CALL PFTOPDG(X,SCALE,DXPDF).                                  *
 *                                                                     *
 *     All variables are defined as DOUBLE PRECISION as before, and    *
 *  X and SCALE are INPUTs provided by the user having the same        *
 *  meaning as above (see Section 2.2).                                *
 *                                                                     *
 *     The subroutine PFTOPDG outputs a vector DXPDF(-6:6) of which    *
 *  its variables have the following meaning:                          *
 *                                                                     *
 *       DXPDF(0)  = gluon                                             *
 *       DXPDF(1)  = sum of down valence + down_bar quark              *
 *       DXPDF(2)  = sum of up valence + up_bar quark                  *
 *       DXPDF(3)  = strange quark                                     *
 *       DXPDF(4)  = charm quark                                       *
 *       DXPDF(5)  = bottom quark                                      *
 *       DXPDF(6)  = top quark                                         *
 *  and  DXPDF(-1) to DXPDF(-6) are the corresponding antiquarks.      *
 *                                                                     *
 *     DXPDF(-1) = DXPDF(-2) stands for down_bar = up_bar or           *
 *  (down_bar + up_bar)/2 in almost all parametrisations. In case      *
 *  where u_bar is not equal to d_bar (like for the retracted MRS S0,  *
 *  D0 and D- sets, the updated MRS S0', D0' and D-' in MS_bar and DIS *
 *  renormalisation scheme as well as the upgraded KMRS set B0 and one *
 *  of the MT series) DXPDF(-1) = d_bar and DXPDF(-2) = u_bar. In all  *
 *  sets is DXPDF(3) = DXPDF(-3), DXPDF(4) = DXPDF(-4),                *
 *  DXPDF(5) = DXPDF(-5) and DXPDF(6) = DXPDF(-6) so far.              *
 *                                                                     *
 *     Please note that also PFTOPDG returns X * parton distribution   *
 *  function !                                                         *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *  2.4  Error Handling                                                *
 *                                                                     *
 *                                                                     *
 *     A note of caution should be addressed here. All structure       *
 *  functions are limited in X and Q**2 which range for most of the    *
 *  sets from                                                          *
 *                                                                     *
 *       4 - 10            < Q**2 < 10**6 - 10**8 (GeV/c)**2           *
 *       10**{-5} - 10**{-4} < X < 1.                                  *
 *                                                                     *
 *     If the user wants to extend the X or Q**2 ranges, it is pos-    *
 *  sible for some of the structure function sets, but the result      *
 *  should be looked at with great caution. If in doubt, please check  *
 *  with the authors of the structure function set. The program is     *
 *                                                                     *
 ******************************* - 14 - ********************************
1***********************************************************************
 *                                                                     *
 *  protected against calculations of Q**2 values below Q2MIN (in this *
 *  case Q**2 = Q2MIN) and of X values in unphysical regions (X < 0 or *
 *  X > 1). The execution of the program is stopped in the later case. *
 *                                                                     *
 *     A COMMON block can be accessed                                  *
 *                                                                     *
 *       COMMON/W50513/ XMIN,XMAX,Q2MIN,Q2MAX                          *
 *                                                                     *
 *  where                                                              *
 *                                                                     *
 *       XMIN   = minimum allowed X value                              *
 *       XMAX   = maximum allowed X value                              *
 *       Q2MIN  = minimum allowed Q**2 value (in (GeV/c)**2)           *
 *       Q2MAX  = maximum allowed Q**2 value (in (GeV/c)**2)           *
 *                                                                     *
 *  is given for each set of structure functions. To control how often *
 *  the X or Q**2 ranges have been exceeded during execution of the    *
 *  user's program a call to the subroutine                            *
 *                                                                     *
 *       CALL PDFSTA                                                   *
 *                                                                     *
 *  at the termination phase of his/her MAIN program allows to print   *
 *  a summary of these errors. If the COMMON block                     *
 *                                                                     *
 *       COMMON/W50514W/ PDFWGT                                        *
 *                                                                     *
 *  where                                                              *
 *                                                                     *
 *       PDFWGT = weight                                               *
 *                                                                     *
 *  is filled by the user on an event-to-event basis the statistics    *
 *  in the error summary takes into account this weight.               *
 *                                                                     *
 *     The print flag in the COMMON/W50510/IFLPRT may be set to        *
 *  IFLPRT = 3 to print an error message EACH TIME a limit is exeeded  *
 *  during job execution. Please note that setting this flag may pro-  *
 *  duce an ENORMOUS AMOUNT of output!                                 *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *  2.5  Features                                                      *
 *                                                                     *
 *                                                                     *
 *     If NPTYPE = 1, NGROUP = 1 and NSET = 1, the user has the possi- *
 *  bility to use his/her own private set of structure functions       *
 *  (i.e. CDHS effective structure functions could be inseted here,    *
 *  or the user could implement a new set of structure functions for   *
 *  testing purposes). To insert the private set the DUMMY subroutine  *
 *                                                                     *
 *    SUBROUTINE STRPRIV(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL) *
 *                                                                     *
 *  should be activated by the user. All variables have to be defined  *
 *  as DOUBLE PRECISION. The variables have the same meaning as des-   *
 *  cribed for the subroutine STRUCTM in Section 2.2 (for explanations *
 *  see there).                                                        *
 *                                                                     *
 *                                                                     *
 ******************************* - 15 - ********************************
1***********************************************************************
 *                                                                     *
 *     For the calculation of the strong coupling constant alpha(s)    *
 *  (see Section 3) with the user's private set of structure functions *
 *  the user should provide a value for the variables QCDL4 and QCDL5  *
 *  via the subroutine PDFSET at the initialization phase. In fact,    *
 *  because only the parameter QCDL5 is used in the alpha(s) calcul-   *
 *  ation, only this definition is mandatory.                          *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *  2.6  PION Structure Functions                                      *
 *                                                                     *
 *                                                                     *
 *     All the available sets of structure functions of the PION have  *
 *  been added to the package as well:                                 *
 *                                                                     *
 *                                                                     *
 *   -  the old Owens (OW-P) sets 1 and 2 [21],                        *
 *   -  the newer sets 1 to 3 of Aurenche et al. (ABFKW-P) [22],       *
 *   -  the very new Sutton, Martin, Roberts and Stirling (SMRS-P)     *
 *      sets 1 to 3 [23] and                                           *
 *   -  the very new Gl"uck, Reya and Vogt (GRV-P) [24] sets of LO     *
 *      and NLL calculations.                                          *
 *                                                                     *
 *     The calling sequence to the pion set of structure functions is  *
 *  kept identical to those described in Section 2.2 for the nucleon   *
 *  structure fuctions. It can be either                               *
 *                                                                     *
 *       CALL STRUCTM(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)    *
 *                                                                     *
 *  or                                                                 *
 *                                                                     *
 *       CALL PFTOPDG(X,SCALE,DXPDF)                                   *
 *                                                                     *
 *  or the old version of the main steering routine                    *
 *                                                                     *
 *       CALL STRUCTF(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,TOP,GL)          *
 *                                                                     *
 *  where, as before, the INPUTs are to be provided by the user in     *
 *  DOUBLE PRECISION. The OUTPUT has the same definitions as before    *
 *  (see Section 2.2). In case of the pion structure functions it is   *
 *  also returned  X * parton distribution function with               *
 *  UPV = DNV, SEA = STR (USEA = DSEA = STR) and TOP = 0               *
 *  (DXPDF(1) = DXPDF(2), DXPDF(3) = DXPDF(-1) = DXPDF(-2) and         *
 *  DXPDF(6) = 0).                                                     *
 *                                                                     *
 *     A list of available sets of PION structure functions in the     *
 *  package is given in Table 2.                                       *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 ******************************* - 16 - ********************************
1***********************************************************************
 *                                                                     *
 *  2.7  PHOTON Structure Functions                                    *
 *                                                                     *
 *                                                                     *
 *     In the previous version, all the available sets of structure    *
 *  functions of the PHOTON have been added to the package as well,    *
 *                                                                     *
 *   -  the old Duke and Owens (DO-G) LO and NLL sets [25],            *
 *   -  the old Drees and Grassie (DG-G) sets 1 to 4 [26],             *
 *   -  the new Abramowicz, Charchula and Levy (LAC-G) [27]            *
 *      sets 1 to 3,                                                   *
 *   -  the new Gordon and Storrow (GS-G) [28] sets, the LO sets 1     *
 *      and 2 and the NLL set,                                         *
 *   -  the new Gl"uck, Reya and Vogt (GRV-G) [29] sets, the LO, NLL   *
 *      and the leading term of the NLL calculations, and              *
 *   -  the new Aurenche et al. (ACFGP-G) [30] NLL sets 1 and 2,       *
 *      without and with massive charm.                                *
 *                                                                     *
 *     The calling sequence to the photon set of structure functions   *
 *  is also kept identical to those described in Section 2.2 for the   *
 *  nucleon structure fuctions. It can be either                       *
 *                                                                     *
 *       CALL STRUCTM(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)    *
 *                                                                     *
 *  or                                                                 *
 *                                                                     *
 *       CALL PFTOPDG(X,SCALE,DXPDF)                                   *
 *                                                                     *
 *  or the old version of the main steering routine                    *
 *                                                                     *
 *       CALL STRUCTF(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,TOP,GL)          *
 *                                                                     *
 *  where, as before, the INPUTs are to be provided by the user in     *
 *  DOUBLE PRECISION. The OUTPUT has the same definitions as before    *
 *  (see Section 2.2). In case of the photon structure functions it    *
 *  is also returned  X * parton distribution function (note that the  *
 *  alpha(QED) has to be taken care of by the user !) with UPV = USEA, *
 *  DNV = DSEA and TOP = 0 (DXPDF(1) = 2 * DXPDF(-1), DXPDF(2) =       *
 *  2 * DXPDF(-2), as a result of the definitions - see Section 2.3 -, *
 *  and DXPDF(6) = 0), while all the other quarks are set identical to *
 *  their antiquark densities.                                         *
 *                                                                     *
 *     A list of available sets of PHOTON structure functions in the   *
 *  package is given in Table 3.                                       *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 ******************************* - 17 - ********************************
1***********************************************************************
 *                                                                     *
 *  2.8  How to find PDFLIB                                            *
 *                                                                     *
 *                                                                     *
 *     At CERN the PDFLIB library and corresponding documentation are  *
 *  available on all systems in the standard areas.                    *
 *                                                                     *
 *    System    Library                     Documentation              *
 *                                                                     *
 *    CERNVM    PDFLIB TXTLIB Q             PDFLIB DOC Q               *
 *    VXCERN    CERN:[PRO.LIB]PDFLIB.OLB    CERN:[PRO.DOC]PDFLIB.PS    *
 *    Unix      /cern/pro/lib/libpdflib.a   /cern/pro/doc/pdf.ps       *
 *    systems                                                          *
 *                                                                     *
 *     The package may be accessed via the CERNLIB command:            *
 *                                                                     *
 *       CERNLIB PDFLIB                                                *
 *                                                                     *
 *                                                                     *
 *     Whenever you use this program, please give reference to the     *
 *  following papers:                                                  *
 *                                                                     *
 *       H. Plothow-Besch, 'PDFLIB: a library of all available parton  *
 *       density functions of the nucleon, the pion and the photon     *
 *       and the corresponding alpha(s) calculations',                 *
 *       Comp. Phys. Comm. 75 (1993) 396-416.                          *
 *                                                                     *
 *   and/or                                                            *
 *                                                                     *
 *       H. Plothow-Besch, 'PDFLIB: Structure Functions and            *
 *       alpha_s Calculations', User's Manual - Version 1.00,          *
 *       W5051 PDFLIB, 1991.03.21, CERN-PPE.                           *
 *                                                                     *
 *   and/or                                                            *
 *                                                                     *
 *       H. Plothow-Besch, 'PDFLIB: Nucleon, Pion and Photon Parton    *
 *       Density Functions and alpha_s Calculations',                  *
 *       User's Manual - Version 6.06, W5051 PDFLIB, 1995.03.15,       *
 *       CERN-PPE.                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *  3.  alpha(s) Calculation                                           *
 *                                                                     *
 *                                                                     *
 *     Within the same package a program is provided to calculate the  *
 *  strong coupling constant alpha(s) to second order as a function    *
 *  of Lambda(QCD) of five flavours and the desired number of fla-     *
 *  vours (NFL) for the selected set of structure functions, which     *
 *  fixes Lambda(QCD). The formula on which the calculations are       *
 *  based upon can be found in Ref. [33]. The same three parameters,   *
 *  NPTYPE, NGROUP and NSET, which select a structure function set,    *
 *  is used to steer the calculation of the alpha(s) value at a given  *
 *  scale from the Lambda(QCD) defined in the selected PDF set.        *
 *                                                                     *
 *      The value of alpha(s) is matched at the thresholds q = m_q.    *
 *  When invoked with NFL < 0, it chooses NFL as the number of fla-    *
 *                                                                     *
 ******************************* - 18 - ********************************
1***********************************************************************
 *                                                                     *
 *  vours for which the masses are less then q. For the quark masses   *
 *  where thresholds are changed the following values have been used:  *
 *  m_charm = 1.5 GeV/c**2, m_bottom = 4.75 GeV/c**2 and               *
 *    m_top = 100 GeV/c**2.                                            *
 *                                                                     *
 *     The alpha(s) value can be obtained via a call to:               *
 *                                                                     *
 *       FUNCTION ALPHAS2(SCALE)                                       *
 *                                                                     *
 *     The user has to provide the following INPUT:                    *
 *                                                                     *
 *       SCALE  =  QCD scale in GeV                                    *
 *                                                                     *
 *     The function ALPHAS2 returns the following OUTPUT:              *
 *                                                                     *
 *       ALPHAS2 =  alpha strong to second order,                      *
 *                                                                     *
 *  if LO not equal to one. For those structure functions, for which   *
 *  the evolution is done to leading order only, alpha(s) to first     *
 *  order is returned.                                                 *
 *                                                                     *
 *     The same subroutine                                             *
 *                                                                     *
 *       SUBROUTINE PDFSET(PARM,VALUE)                                 *
 *                                                                     *
 *  as described in Section 2.1 should be used to change the default   *
 *  values of the variables NPTYPE, NGROUP, NSET, NFL, LO, TMAS,       *
 *  QCDL4, QCDL5, XMIN, XMAX, Q2MIN and Q2MAX.                         *
 *                                                                     *
 *     The internal COMMON block                                       *
 *                                                                     *
 *       COMMON/W50512/ QCDL4,QCDL5                                    *
 *                                                                     *
 *  where                                                              *
 *                                                                     *
 *       QCDL4  = QCD scale, Lambda_4(QCD), in GeV for four flavours   *
 *       QCDL5  = QCD scale, Lambda_5(QCD), in GeV for five flavours   *
 *                corresponding to QCDL4                               *
 *                                                                     *
 *  provides the actual value of Lambda_4(QCD) and Lambda_5(QCD) used  *
 *  in the alpha(s) calculation for four and five flavours, respectiv- *
 *  ly, for each set of structure functions.                           *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *  4.  Conclusions                                                    *
 *                                                                     *
 *                                                                     *
 *     In the Long Write-up as an example, the NUCLEON structure       *
 *  function distributions at the scale Q**2 = m_W**2 for u and d va-  *
 *  lence quarks as well as for the sea and gluon contributions are    *
 *  shown in Figs. 1 to 4 respectively, as a function of the parton X. *
 *  The average X values for CERN and FNAL as well as for LHC and SSC  *
 *  energies at that scale are indicated. Please, note the logarithmic *
 *  scale for the gluon distribution. The d/u ratio of the valence     *
 *  quarks is shown in Fig. 5. It can be seen from all these figures   *
 *                                                                     *
 ******************************* - 19 - ********************************
1***********************************************************************
 *                                                                     *
 *  these figures that there is a large spread in shape for the dif-   *
 *  ferent sets in all contributions leading to different results in   *
 *  the cross section calculations. It should be pointed out that the  *
 *  spread from the different parametrisations in the SAME renormali-  *
 *  sation scheme is MUCH LARGER than the spread from parametrisations *
 *  in different renormalisation schemes. This is true for ALL parton  *
 *  densities, and in particular at small X values !                   *
 *                                                                     *
 *     A comparison to recent data is shown in Fig. 6, where the ratio *
 *  F_2(n)/F_2(p) at Q**2 = m_W**2 is displayed. This ratio mainly re- *
 *  flects the d/u ratio of valence quarks which is the dominant un-   *
 *  certainty of the W and Z boson cross section ratio in pp or pp_bar *
 *  reactions at lower energies. In fact, the older sets of structure  *
 *  functions, namely the set of Owens and Reya [2] and the set of     *
 *  Baier et al. [3] have only been implemented for completeness.      *
 *  These sets should not be used for cross section calculations of    *
 *  any hard processes at high energies anymore. From Fig. 6 we con-   *
 *  clude that also the other older sets of structure functions,       *
 *  namely the two sets of Duke and Owens [5] and the two sets of      *
 *  Eichten et al. [6] should be used with some care, because they do  *
 *  not fit the recent low energy deep-inelastic lepton-nucleon data   *
 *  from NMC and BCDMS [34]. This is not surprising because these data *
 *  were not available when the sets have been made. The new set of    *
 *  DO 1.1 superseeds the old DO sets 1 and 2 but still gives very     *
 *  limited results in the low X region, even though data are avai-    *
 *  lable there. (Please note that the last data point at the lowest   *
 *  X value shown in Fig. 6 is out of range in Q**2 for all sets of    *
 *  structure functions and that the point before last is very much at *
 *  the limit of the allowed range).                                   *
 *                                                                     *
 *     In Fig. 7 the F_2(p) distribution of the recent data from Ref.  *
 *  [35] at Q**2 = 5 GeV**2, together with the very recent data at     *
 *  Q**2 = 15 GeV**2 from the H1 and the ZEUS experiments [36] at the  *
 *  the ep-collider HERA is shown as a function of X. Overlayed on the *
 *  same Figure are a few theoretical predictions of more recent par-  *
 *  ton parametrisations [12,13,14,15,16,20] at Q**2 = 15 GeV**2.      *
 *  Please note that most of these parametrisations are NOT real pre-  *
 *  dictions because the data are fully or partially used in the fits. *
 *  Apart from Ref. [13] set B2, from Ref. [14] set HO, from Ref. [16] *
 *  set D-' and (H) and from Ref. [20] the new set 2pM, most of these  *
 *  more recent parametrisations fail to describe well enough the low  *
 *  X behaviour of the present data.                                   *
 *                                                                     *
 *     From Figs. 6 and 7 we conclude that the recent sets of parton   *
 *  densities, nameley the GRV set HO [14], but also the old MT set    *
 *  B2 [13] - both sets are real predictions -, and the updated MRS    *
 *  set (H) [16], could be preferred for all theoretical predictions   *
 *  involving structure functions, because they fit the present data   *
 *  best which, at the moment, still suffer from large uncertainties.  *
 *                                                                     *
 *     As an example, the u_bar, strange and gluon densities for the   *
 *  different sets of structure functions of the PHOTON are shown in   *
 *  Figs. 7, 8 and 9 of the Long Write-up, respectively, as a function *
 *  of X at the scale Q**2 = 100 (GeV/c)**2. In Fig. 10 a comparison   *
 *  of F_2(gamma)/alpha(QED) as a function of X at Q**2 = 100 (GeV/c)**2
 *  for the different photon structure function sets is made. The      *
 *                                                                     *
 ******************************* - 20 - ********************************
1***********************************************************************
 *                                                                     *
 *  charm quark density has been taken into account. The theoretical   *
 *  predictions are compared with experimental data from the JADE      *
 *  experiment [37] at the same Q**2 value. Because of the large       *
 *  experimental uncertainties no distinction between the              *
 *  different parton density sets can be made yet.                     *
 *                                                                     *
 *     Coming new sets of structure functions could and should be      *
 *  easily implemented in this package. Authors of new sets are kindly *
 *  requested to provide us with the relevant information.             *
 *                                                                     *
 *     Please return any problems, questions, suggestions for im-      *
 *  provements to the author of the package                            *
 *  (e-mail address: PLOTHOW@CERNVM.CERN.CH).                          *
 *                                                                     *
 *     Requests for the Long Write-up of the User's Manual can be ad-  *
 *  dressed either to the author or to the CERN Program Library Office.*
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *  5.  Acknowledgements                                               *
 *                                                                     *
 *                                                                     *
 *     The author of the package would like to thank                   *
 *  M.A. Marquina/CERN-CN for help and support making this package     *
 *  available to the Physics community as part of the CERN Program     *
 *  Library. In addition we give many thanks to G. Folger/CERN-CN and  *
 *  I. McLaren/CERN-CN for thoroughly testing out the package on the   *
 *  different system areas.                                            *
 *                                                                     *
 *     We also would like to thank J.M. Marraffino/FNAL for the help-  *
 *  ful comments and the correction set which he supplied.             *
 *                                                                     *
 *     And finally we acknowledge K. Charchula/DESY who provided the   *
 *  source code of some of the photon structure functions which had    *
 *  previously been implemented in the package.                        *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 ******************************* - 21 - ********************************
1***********************************************************************
 *                                                                     *
 *  Appendix A                                                         *
 *                                                                     *
 *                                                                     *
 *     As an example how to use the package, the FORTRAN code to ob-   *
 *  tain Fig. 1 of the Long Write-up which shows the parton density    *
 *  function distributions for the up valence quark as a function of   *
 *  the parton X for all sets of structure functions is given below:   *
 *                                                                     *
 *      PROGRAM PDFUPV                                                 *
 *  C  define the maximum number of PDF sets, the name and the cross reference
 *      PARAMETER (NPTYMX = 3, NGRMAX = 8, NSETMX = 41)                *
 *      COMMON /W505120/ NPGSMX(NPTYMX,NGRMAX),NSETFL(NPTYMX,NGRMAX,NSETMX)
 *      CHARACTER*8 SFNAME(NPTYMX,NGRMAX,NSETMX)                       *
 *      COMMON /W505110/ SFNAME                                        *
 *      COMMON /W505122/ MODECR(NPTYMX,NGRMAX,NSETMX)                  *
 *  C  define HBOOK settings                                           *
 *      PARAMETER (NHBMEM = 500000)                                    *
 *      COMMON /PAWC/ HMEM(NHBMEM)                                     *
 *      PARAMETER (NB=900, ID=100)                                     *
 *  C  define DOUBLE PRECISION variables for calling sequence to STRUCTM
 *      DOUBLE PRECISION DX,DSCALE,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL
 *      DOUBLE PRECISION ALF,ALPHAS2                                   *
 *      REAL    X, SCALE, UPV, DNV, USEA, DSEA, STR, CHM, BOT, TOP, GL *
 *      COMMON/W50516/ FIRST                                           *
 *      LOGICAL FIRST                                                  *
 *      CHARACTER*20 PARM(20)                                          *
 *      DOUBLE PRECISION VAL(20)                                       *
 *      DATA SCALE/80.140/                                             *
 *      DATA X00/0.002/,DX0/0.001/,XLOW/0.0015/,XUP/0.9015/            *
 *  C                                                                  *
 *      CALL HLIMIT(NHBMEM)                                            *
 *      DSCALE=SCALE                                                   *
 *  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                                                     *
 *      NHB = 0                                                        *
 *      DO 20 IGR = 1,NGRMAX                                           *
 *         IF(NPGSMX(NPTYPE,IGR).EQ.0) GOTO 20                         *
 *         DO ISET = 1,NPGSMX(NPTYPE,IGR)                              *
 *  C  book histograms for each set of SF separately                   *
 *         NHB = NHB + 1                                               *
 *         CALL HBOOK1(ID+NHB,'U Valence quark ',NB,XLOW,XUP,0.)       *
 *  C  force label printing for each set of SF (not only the 1st)      *
 *         FIRST = .TRUE.                                              *
 *  C  define 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)                                       *
 *                                                                     *
 *                                                                     *
 ******************************* - 22 - ********************************
1***********************************************************************
 *                                                                     *
 *  Example of Fig. 1, cont'd                                          *
 *                                                                     *
 *  C  loop over all X bins                                            *
 *        DO 10 I=1,NB                                                 *
 *           X = X00 + (I-1)*DX0                                       *
 *           DX = X                                                    *
 *           IF(X.LT.XLOW .OR. X.GT.XUP) GOTO 10                       *
 *           CALL STRUCTM(DX,DSCALE,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
 *           UPV=DUPV                                                  *
 *           IF(X.GT.0.499 .AND. X.LE.0.500) WRITE(6,1000) X,SCALE,UPV *
 *  1000     FORMAT(/,'  X= ',F6.4,'  Q= ',F6.3,'  UPV= ',8.4)         *
 *           CALL HF1(ID+NHB,X,UPV)                                    *
 *  10    CONTINUE                                                     *
 *  C  get alpha(s) for selected set of SF at Q = SCALE                *
 *        ALF = ALPHAS2(DSCALE)                                        *
 *        WRITE(6,4000) NPTYPE,IGR,ISET,SFNAME(NPTYPE,IGR,ISET),       *
 *   +                  MODECR(NPTYPE,IGR,ISET),NSF                    *
 *        WRITE(6,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             *
 *  C  get error summary for each set of SF                            *
 *        CALL PDFSTA                                                  *
 *     ENDDO                                                           *
 *  20 CONTINUE                                                        *
 *     STOP                                                            *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 ******************************* - 23 - ********************************
1***********************************************************************
 *                                                                     *
 *     This code will produce the following output:                    *
 *                                                                     *
 *                                                                     *
 *         1. HBOOK Output :                                           *
 *                                                                     *
 *     one histogram for each set of structure functions. The histo-   *
 *  gram output has then been accessed, and each histogram has been    *
 *  superimposed on one single frame using the PAW package [38] to     *
 *  obtain Fig. 1 of the Long Write-up.                                *
 *                                                                     *
 *                                                                     *
 *         2. Print Output (Example is given for NPTYPE = 1,           *
 *  NGROUP = 3, NSET = 41 only :                                       *
 *                                                                     *
 * ___________________________________________________________________ *
 *|                                                                   |*
 *| 1*****  PDFLIB Version:  6.06  Released on  950315  at  17.05  in |*
 *| the CERN Program Library     W5051  *****                         |*
 *|  *****  Library compiled on      950315  at        2337   *****   |*
 *|                                                                   |*
 *| Nucleon PDFs:  Ngroup = 3,  Nset = 41,  for MRS Set (G) (L255-MSb)|*
 *| ------------------------------------------------------------------|*
 *|                                             Structure Functions   |*
 *|                                             -------------------   |*
 *|                                                                   |*
 *| X= 0.5000  Q= 80.140  UPV=   0.1176                               |*
 *|                                                                   |*
 *| Nptype = 1  Ngroup = 3  Nset = 41 Name = "MRS-G   "  CrMode =  -1  HBId = 60
 *| Scale  =  80.1400  alpha(s) =  0.1154                             |*
 *|                                                                   |*
 *| PDFLIB : Summary from PDFSTA                                      |*
 *| Nptype = 1  Ngroup = 3  Nset = 41  Name = "MRS-G   "  CrMode =  -1|*
 *| Nfl    = -5  LO = 2  Tmas = 180.00 GeV/c**2                       |*
 *| QCDL4  = 0.2550 GeV,  QCDL5 =   0.1708 GeV                        |*
 *| Xmin   = 0.10E-04,  Xmax =  0.99999E+00,  Q2min =  5.000 (GeV/c)**2,
 *|   Q2max =  0.27E+10 (GeV/c)**2                                    |*
 *|                                                                   |*
 *| PDFSTA:    NO errors occured                                      |*
 *|___________________________________________________________________|*
 *                                                                     *
 *                                                                     *
 *     Apart from the two top lines the print output shown above is    *
 *  repeated for each set of structure functions.                      *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 ******************************* - 24 - ********************************
1***********************************************************************
 *                                                                     *
 *  Appendix B                                                         *
 *                                                                     *
 *                                                                     *
 *     A list of subroutines, functions, COMMON blocks and parameters  *
 *  of RELEVANCE to the user which are used in PDFLIB is given below:  *
 *                                                                     *
 *                                                                     *
 *     List of relevant subprograms and COMMON blocks in PDFLIB:       *
 *                                                                     *
 *  _________________________________________________________________  *
 *  |                                                               |  *
 *  |  (S = Subroutine, F = Function, C = COMMON, P = Parameter)    |  *
 *  |                                                               |  *
 *  | ------------------------------------------------------------- |  *
 *  |                                                               |  *
 *  |  S  PDFSET   : to set all parameters for PDFLIB               |  *
 *  |  S  STRUCTF  : to access PDFLIB with output in PDFLIB format  |  *
 *  |  S  STRUCTM  : to access PDFLIB with output in PDFLIB format  |  *
 *  |                (with u_bar not equal to d_bar)                |  *
 *  |  S  PFTOPDG  : to access PDFLIB with output in PDG format     |  *
 *  |  S  PDFSTA   : to print summary of parameters and error statistics
 *  |  F  ALPHAS2  : to calculate alpha(s) to second order QCD      |  *
 *  |                                                               |  *
 *  |  C  W50510   : IFLPRT - to set print flag                     |  *
 *  |  C  W50511   : NPTYPE,NGROUP,NSET,MODE,NFL,LO,TMAS            |  *
 *  |                - to select a structure function set, to       |  *
 *  |                define number of flavours, flag for leading    |  *
 *  |                order and value for the mass of the top quark  |  *
 *  |  C  W50512   : QCDL4,QCDL5 - contains Lambda_4(QCD) and       |  *
 *  |                Lambda_5(QCD) (in GeV)                         |  *
 *  |  C  W50513   : XMIN,XMAX,Q2MIN,Q2MAX - contains minimum and   |  *
 *  |                maximum values of X and Q**2                   |  *
 *  |  C  W50516   : FIRST - to force printing name of structure    |  *
 *  |                function set in case of multiple choices       |  *
 *  |  C  W50519   : NEWVER - flags NEW/OLD version of PDFLIB format|  *
 *  |  C  W505110  : SFNAME - contains NAME (CHARACTER*8 array with |  *
 *  |                dimension SFNAME(NPTYMX,NGRMAX,NSETMX)) for    |  *
 *  |                each set of PDFs                               |  *
 *  |  C  W505120  : NPGSMX(NPTYMX,NGRMAX),NSETFL(NPTYMX,NGRMAX,NSETMX)*
 *  |                - contains maximum number of stucture functions|  *
 *  |                per particle type and author group, and FLAG   |  *
 *  |                for retracted PDF set (NSETFL(I,J,K) = 0),     |  *
 *  |  C  W505121  : NPTYCR(MODEMX),NGROCR(MODEMX),NSETCR(MODEMX)   |  *
 *  |                - cross reference to the NPTYPE, NGROUP and    |  *
 *  |                NSET value (NEW format) for a given MODE       |  *
 *  |                value (OLD format)                             |  *
 *  |  C  W505122  : MODECR(NPTYMX,NGRMAX,NSETMX) - cross reference |  *
 *  |                to the MODE value (OLD format) for a given     |  *
 *  |                NPTYPE, NGROUP and NSET value (NEW format)     |  *
 *  |  P  W5051P2  : MODEMX - maximum number of structure function  |  *
 *  |                sets in the OLD PDFLIB format                  |  *
 *  |              : NPTYMX,NGRMAX,NSETMX - maximum number of       |  *
 *  |                structure function sets in the NEW PDFLIB      |  *
 *  |                format for particle type, author group number  |  *
 *  |                and number of sets within an author group      |  *
 *  |  P  W5051P7  : L6 - to define logical print unit number       |  *
 *  |_______________________________________________________________|  *
 ******************************* - 25 - ********************************
1***********************************************************************
 *                                                                     *
 *  Appendix C                                                         *
 *                                                                     *
 *                                                                     *
 *     A complete list of subroutines, functions, COMMON blocks and    *
 *  parameters which are used in PDFLIB is given below:                *
 *                                                                     *
 *                                                                     *
 *     Complete list of subprograms and COMMON blocks in PDFLIB:       *
 *                       (in alphabetical order)                       *
 *                                                                     *
 *  _________________________________________________________________  *
 *  |                                                               |  *
 *  |  (S = Subroutine, F = Function, C = COMMON, P = Parameter)    |  *
 *  |                                                               |  *
 *  | ------------------------------------------------------------- |  *
 *  |                                                               |  *
 *  |  S  ABFKW1   : ABFKW SF set 1 in the pion (NP=2,NG=6,NS=1)    |  *
 *  |  S  ABFKW2   : ABFKW SF set 2 in the pion (NP=2,NG=6,NS=2)    |  *
 *  |  S  ABFKW3   : ABFKW SF set 3 in the pion (NP=2,NG=6,NS=3)    |  *
 *  |  S  ACFGP1   : ACFGP SF set 1 in the photon (NP=3,NG=6,NS=1)  |  *
 *  |  S  ACFGP2   : ACFGP SF set 1 in the photon (mc) (NP=3,NG=6,NS=2)*
 *  |  F  AFCPLU   : to get PD with massive charm for ACFGP         |  *
 *  |  F  AFGETFV  : to interpolate PDs in grid for ACFGP           |  *
 *  |  F  AFPOLIN  : to interpolate PDs in grid for ACFGP           |  *
 *  |  F  AFRATIN  : to interpolate PDs in grid for ACFGP           |  *
 *  |  F  ALPHAS2  : to calculate alpha(s) to second order QCD      |  *
 *  |  F  AUGETFV  : to interpolate logarithmically PDs in grid     |  *
 *  |  S  AURGAM   : to get the PDs in the photon for ACFGP-G       |  *
 *  |  S  AURPI1   : to get the PDs in the pion for ABFKW-P set 1   |  *
 *  |  S  AURPI2   : to get the PDs in the pion for ABFKW-P set 2   |  *
 *  |  S  AURPI3   : to get the PDs in the pion for ABFKW-P set 3   |  *
 *  |  S  AURPR    : to get the PDs in the nucleon for ABFOW        |  *
 *  |  F  BBETA    : DOUBLE PRECISION beta function using DGAMMA    |  *
 *  |  F  BETA     : SINGLE PRECISION beta function using GAMMA     |  *
 *  |  S  BXF3     : BEP SF set for up and down valence quarks      |  *
 *  |  S  BXGLUE   : BEP SF set for gluons                          |  *
 *  |  S  BXQBAR   : BEP SF set for antiquarks                      |  *
 *  |  F  CTQ1PD   : CTEQ SF sets 1M, 1MS, 1ML, 1D and 1L           |  *
 *  |  F  CTQ1PF   : Description of CTEQ sets 1M, 1MS, 1ML, 1D, 1L  |  *
 *  |  F  CTQ2OPD  : CTEQ SF sets 2M, 2MS, 2ML, 2D and 2L           |  *
 *  |  F  CTQ2OPF  : Description of CTEQ sets 2M, 2MS, 2ML, 2D, 2L  |  *
 *  |  S  CTQ2OPS  : CTEQ SF sets 2M, 2MS, 2ML, 2D and 2L           |  *
 *  |  F  CTQ2PPD  : CTEQ SF sets 2pM, 2pMS, 2pML, 2pD and 2pL      |  *
 *  |  F  CTQ2PPF  : Description of CTEQ sets 2pM, 2pMS, 2pML, 2pD, 2pL*
 *  |  S  CTQ2PPS  : CTEQ SF sets 2pM, 2pMS, 2pML, 2pD and 2pL      |  *
 *  |  F  DBFINT   : to interpolate linearly PDs in grid (double precision)
 *  |  S  DFLM1    : interface to DFLM 160 (NP=1,NG=2,NS=7)         |  *
 *  |  S  DFLM2    : interface to DFLM 260 (NP=1,NG=2,NS=8)         |  *
 *  |  S  DFLM3    : interface to DFLM 360 (NP=1,NG=2,NS=9)         |  *
 *  |  S  DFLM4    : interface to DFLM soft valence quarks (NP=1,NG=2,NS=1)
 *  |  S  DFLM5    : interface to DFLM hard valence quarks (NP=1,NG=2,NS=2)
 *  |  S  DFLM6    : interface to DFLM soft gluons (NP=1,NG=2,NS=3) |  *
 *  |  S  DFLM7    : interface to DFLM hard gluons (NP=1,NG=2,NS=4) |  *
 *  |  S  DFLM8    : interface to DFLM LO average fit (NP=1,NG=2,NS=5) *
 *  |  S  DFLM9    : interface to DFLM NLL average fit (NP=1,NG=2,NS=6)*
 *  |  S  DGPHO1   : DG SF set 1 in the photon (NP=3,NG=2,NS=1)     |  *
 *  |_______________________________________________________________|  *
 ******************************* - 26 - ********************************
1***********************************************************************
 *  _________________________________________________________________  *
 *  |  Complete list of subprograms, cont'd                         |  *
 *  |                                                               |  *
 *  |  S  DGPHO2   : DG SF set 2 in the photon (NP=3,NG=2,NS=2)     |  *
 *  |  S  DGPHO3   : DG SF set 3 in the photon (NP=3,NG=2,NS=3)     |  *
 *  |  S  DGPHO4   : DG SF set 4 in the photon (NP=3,NG=2,NS=4)     |  *
 *  |  S  DOPHO1   : DO SF set 1 in the photon (NP=3,NG=1,NS=1)     |  *
 *  |  S  DOPHO2   : DO SF set 2 in the photon (NP=3,NG=1,NS=2)     |  *
 *  |  S  FXAVER   : DFLM SF set of LO central average fit          |  *
 *  |  S  FXG160   : DFLM SF set of NLL with QCDL4 = 160            |  *
 *  |  S  FXG260   : DFLM SF set of NLL with QCDL4 = 260            |  *
 *  |  S  FXG360   : DFLM SF set of NLL with QCDL4 = 360            |  *
 *  |  S  FXGHAR   : DFLM SF set of LO hard gluon                   |  *
 *  |  S  FXGSOF   : DFLM SF set of LO soft gluon                   |  *
 *  |  S  FXNLLA   : DFLM SF set of NLL central average fit         |  *
 *  |  S  FXVHAR   : DFLM SF set of LO hard valence                 |  *
 *  |  S  FXVSOF   : DFLM SF set of LO soft valence                 |  *
 *  |  F  GAMFUN   : to calculate GAMMA funct. in either SP or DP prec.*
 *  |  F  GAMMADO  : to calculate GAMMA funct. in SP (Hastings)     |  *
 *  |  F  GHRDV    : GHR SF set for down quarks                     |  *
 *  |  F  GHRGL    : GHR SF set for gluons                          |  *
 *  |  F  GHRUV    : GHR SF set for up quarks                       |  *
 *  |  F  GHRXI    : GHR SF set for antiquarks                      |  *
 *  |  F  GHRXS    : GHR SF set for strange quarks                  |  *
 *  |  F  GRVFGP   : GRV SF set for gluons in the pion              |  *
 *  |  F  GRVFQBP  : GRV SF set for s, c and b quarks in the pion   |  *
 *  |  F  GRVFV    : GRV SF set for up + down valence quarks        |  *
 *  |  F  GRVFVP   : GRV SF set for valence quarks in the pion      |  *
 *  |  F  GRVFW    : GRV SF set for antiquarks and gluons           |  *
 *  |  F  GRVFWS   : GRV SF set for strange, charm and bottom quarks|  *
 *  |  S  GRVGAH0  : GRV SF set lead. terms of NLL in the photon (NP=3,NG=5,NS=1
 *  |  S  GRVGAHO  : GRV SF set of NLL in the photon (NP=3,NG=5,NS=2)  *
 *  |  S  GRVGALO  : GRV SF set of LO in the photon (NP=3,NG=5,NS=3)|  *
 *  |  F  GRVGF    : GRV SF set for u and d quarks and gluons in the photon
 *  |  F  GRVGFS   : GRV SF set for s, c and b quarks in the photon |  *
 *  |  S  GRVHO    : GRV SF set of NLL (NP=1,NG=5,NS=3)             |  *
 *  |  S  GRVLO    : GRV SF set of LO (NP=1,NG=5,NS=4)              |  *
 *  |  S  GRVPIHO  : GRV SF set of NLL in the pion (NP=2,NG=5,NS=1) |  *
 *  |  S  GRVPILO  : GRV SF set of LO in the pion (NP=2,NG=5,NS=2)  |  *
 *  |  S  GSXCOR   : to get X coordinates for SFGSHO + LO sets 1, 2 |  *
 *  |  S  HMRS1EB  : Description of HMRS sets E, B (1.90-retracted) |  *
 *  |  S  HMRS2EB  : Description of HMRS sets E, B (3.90)           |  *
 *  |  S  HMRS3EB  : Description of HMRS sets E, B (4.90)           |  *
 *  |  S  HMRSGEB  : Description of HMRS sets E, E+, E- and B       |  *
 *  |  S  KMRSEB   : Description of KMRS sets E and B               |  *
 *  |  S  KMRSGEB  : Description of KMRS B0 with L=135,160,200,235  |  *
 *  |  S  MRSDSHP  : Description of MRS sets D0', S0', D-', (H)     |  *
 *  |  S  MRSEB    : Description of MRS sets E and B                |  *
 *  |  S  MRSEBP   : Description of MRS sets E' and B'              |  *
 *  |  S  NEWDO1   : DO SF set 1.1 (NP=1,NG=1,NS=10)                |  *
 *  |  S  PDFSET   : to set all parameters for PDFLIB               |  *
 *  |  S  PDFSTA   : to print summary of param. and error statistics|  *
 *  |  S  PDFVERS  : to get date and time of PDFLIB Library release |  *
 *  |  F  PDXMT    : Description of MT sets 1 to 11                 |  *
 *  |  S  PDZXMT   : MT SF sets 1 to 11                             |  *
 *  |  S  PFTOPDG  : to access PDFLIB with output in PDG format     |  *
 *  |  S  PHLAC1   : interface to LAC-G set 1 in the photon         |  *
 *  |_______________________________________________________________|  *
 ******************************* - 27 - ********************************
1***********************************************************************
 *  _________________________________________________________________  *
 *  |  Complete list of subprograms, cont'd                         |  *
 *  |                                                               |  *
 *  |  S  PHLAC2   : interface to LAC-G set 2 in the photon         |  *
 *  |  S  PHLAC3   : interface to LAC-G set 3 in the photon         |  *
 *  |  S  RXDV     : OR SF set for down valence quarks              |  *
 *  |  S  RXGLUE   : OR SF set for gluons                           |  *
 *  |  S  RXQBAR   : OR SF set for antiquarks                       |  *
 *  |  S  RXUD     : OR SF set for the sum of up + down quarks      |  *
 *  |  S  SEADIS   : BEBC SF set for antiquarks                     |  *
 *  |  S  SFABFOW  : ABFOW SF set (NP=1,NG=6,NS=1)                  |  *
 *  |  S  SFBMA    : BM SF set A (NP=1,NG=7,NS=1)                   |  *
 *  |  S  SFBMB    : BM SF set B (NP=1,NG=7,NS=2)                   |  *
 *  |  S  SFCTQ11  : CTEQ SF set 1L (LO) (NP=1,NG=4,NS=12)          |  *
 *  |  S  SFCTQ12  : CTEQ SF set 1M (NP=1,NG=4,NS=13)               |  *
 *  |  S  SFCTQ13  : CTEQ SF set 1MS (NP=1,NG=4,NS=14)              |  *
 *  |  S  SFCTQ14  : CTEQ SF set 1ML (NP=1,NG=4,NS=15)              |  *
 *  |  S  SFCTQ15  : CTEQ SF set 1D (NP=1,NG=4,NS=16)               |  *
 *  |  S  SFCTQ21  : CTEQ SF set 2L (LO) (NP=1,NG=4,NS=17)          |  *
 *  |  S  SFCTQ22  : CTEQ SF set 2M (NP=1,NG=4,NS=18)               |  *
 *  |  S  SFCTQ23  : CTEQ SF set 2MS (NP=1,NG=4,NS=19)              |  *
 *  |  S  SFCTQ24  : CTEQ SF set 2MF (NP=1,NG=4,NS=20)              |  *
 *  |  S  SFCTQ25  : CTEQ SF set 2ML (NP=1,NG=4,NS=21)              |  *
 *  |  S  SFCTQ26  : CTEQ SF set 2D (NP=1,NG=4,NS=22)               |  *
 *  |  S  SFDISD0  : MRS SF set D0' (DIS) (NP=1,NG=3,NS=33)         |  *
 *  |  S  SFDISDM  : MRS SF set D-' (DIS) (NP=1,NG=3,NS=34)         |  *
 *  |  S  SFDISH   : MRS SF set (H) (NP=1,NG=3,NS=36)               |  *
 *  |  S  SFDISS0  : MRS SF set S0' (DIS) (NP=1,NG=3,NS=32)         |  *
 *  |  S  SFEHLQ1  : EHLQ SF set 1 (NP=1,NG=1,NS=8)                 |  *
 *  |  S  SFEHLQ2  : EHLQ SF set 2 (NP=1,NG=1,NS=9)                 |  *
 *  |  S  SFGRVH   : old GRV SF set of NLL (NP=1,NG=5,NS=1)         |  *
 *  |  S  SFGRVL   : old GRV SF set of LO (NP=1,NG=5,NS=2)          |  *
 *  |  S  SFGSHO   : GS SF set of NLL in the photon (NP=3,NG=4,NS=1)|  *
 *  |  S  SFGSLO1  : GS SF set 1 of LO in the photon (NP=3,NG=4,NS=2)  *
 *  |  S  SFGSLO2  : GS SF set 2 of LO in the photon (NP=3,NG=4,NS=3)  *
 *  |  S  SFKBMR2  : KMRS SF set B- with R2 shadowing (NP=1,NG=3,NS=19)*
 *  |  S  SFKBMR5  : KMRS SF set B- with R5 shadowing (NP=1,NG=3,NS=20)*
 *  |  S  SFLACG1  : LAC SF set 1 in the photon (NP=3,NG=3,NS=1)    |  *
 *  |  S  SFLACG2  : LAC SF set 2 in the photon (NP=3,NG=3,NS=2)    |  *
 *  |  S  SFLACG3  : LAC SF set 3 in the photon (NP=3,NG=3,NS=3)    |  *
 *  |  S  SFMRSD0  : MRS SF set D0 (NP=1,NG=3,NS=27)                |  *
 *  |  S  SFMRSDM  : MRS SF set D- (NP=1,NG=3,NS=28)                |  *
 *  |  S  SFMRSS0  : MRS SF set S0 (NP=1,NG=3,NS=26)                |  *
 *  |  S  SFMSBB0  : MRS SF set B0' (NP=1,NG=3,NS=                  |  *
 *  |  S  SFMSBD0  : MRS SF set D0' (NP=1,NG=3.NS=30)               |  *
 *  |  S  SFMSBDM  : MRS SF set D-' (NP=1,NG=3,NS=31)               |  *
 *  |  S  SFMSBH   : MRS SF set (H) (NP=1,NG=3,NS=35)               |  *
 *  |  S  SFMSBS0  : MRS SF set S0' (NP=1,NG=3,NS=29)               |  *
 *  |  S  SFMT10   : interface to MT set 6 (MS_bar)                 |  *
 *  |  S  SFMT11   : interface to MT set LO                         |  *
 *  |  S  SFMTU1   : interface to MT set S1 (DIS)                   |  *
 *  |  S  SFMTU2   : interface to MT set B1 (DIS)                   |  *
 *  |  S  SFMTU3   : interface to MT set B2 (DIS)                   |  *
 *  |  S  SFMTU4   : interface to MT set E1 (DIS)                   |  *
 *  |  S  SFMTU5   : interface to MT set 6 (u_bar + d_bar)/2 (DIS)  |  *
 *  |  S  SFMTU6   : interface to MT set S1 (MS_bar)                |  *
 *  |  S  SFMTU7   : interface to MT set B1 (MS_bar)                |  *
 *  |_______________________________________________________________|  *
 ******************************* - 28 - ********************************
1***********************************************************************
 *  _________________________________________________________________  *
 *  |  Complete list of subprograms, cont'd                         |  *
 *  |                                                               |  *
 *  |  S  SFMTU8   : interface to MT set B2 (MS_bar)                |  *
 *  |  S  SFMTU9   : interface to MT set E1 (MS_bar)                |  *
 *  |  S  SMRSP31  : SMRS SF set 1 in the pion (NP=2,NG=3,NS=1)     |  *
 *  |  S  SMRSP32  : SMRS SF set 2 in the pion (NP=2,NG=3,NS=2)     |  *
 *  |  S  SMRSP33  : SMRS SF set 3 in the pion (NP=2,NG=3,NS=3)     |  *
 *  |  S  SMRSPI   : Description of SMRS sets 1, 2, 3 in the pion   |  *
 *  |  S  SRCTQ21  : CTEQ SF set 2pL (LO) (NP=1,NG=4,NS=23)         |  *
 *  |  S  SRCTQ22  : CTEQ SF set 2pM (NP=1,NG=4,NS=24)              |  *
 *  |  S  SRCTQ23  : CTEQ SF set 2pMS (NP=1,NG=4,NS=25)             |  *
 *  |  S  SRCTQ24  : CTEQ SF set 2pMF (NP=1,NG=4,NS=26)             |  *
 *  |  S  SRCTQ25  : CTEQ SF set 2pML (NP=1,NG=4,NS=27)             |  *
 *  |  S  SRCTQ26  : CTEQ SF set 2pD (NP=1,NG=4,NS=28)              |  *
 *  |  S  STRBEBC  : BEBC SF set (NP=1,NG=1,NS=2)                   |  *
 *  |  S  STRBEP   : BEP SF set (NP=1,NG=1,NS=4)                    |  *
 *  |  S  STRCBP   : MRS SF set B' (NP=1,NG=3,NS=7)                 |  *
 *  |  S  STRCEP   : MRS SF set E' (NP=1,NG=3,NS=6)                 |  *
 *  |  S  STRDO1   : DO SF set 1 (NP=1,NG=1,NS=6)                   |  *
 *  |  S  STRDO2   : DO SF set 2 (NP=1,NG=1,NS=7)                   |  *
 *  |  S  STRGHR   : GHR SF set (NP=1,NG=1,NS=5)                    |  *
 *  |  S  STRH1B   : HMRS SF set B (NP=1,NG=3,NS=9 - retracted)     |  *
 *  |  S  STRH1E   : HMRS SF set E (NP=1,NG=3,NS=8 - retracted)     |  *
 *  |  S  STRH2B   : HMRS SF set B (NP=1,NG=3,NS=11)                |  *
 *  |  S  STRH2E   : HMRS SF set E (NP=1,NG=3,NS=10)                |  *
 *  |  S  STRH3B   : HMRS SF set B (NP=1,NG=3,NS=15)                |  *
 *  |  S  STRH3B1  : HMRS SF set B - L=100 (NP=1,NG=3,NS=16)        |  *
 *  |  S  STRH3B3  : HMRS SF set B - L=300 (NP=1,NG=3,NS=17)        |  *
 *  |  S  STRH3E   : HMRS SF set E (NP=1,NG=3,NS=14)                |  *
 *  |  S  STRH3EM  : HMRS SF set E- (NP=1,NG=3,NS=13)               |  *
 *  |  S  STRH3EP  : HMRS SF set E+ (NP=1,NG=3,NS=12)               |  *
 *  |  S  STRKB0   : KMRS SF set B0 (NP=1,NG=3,NS=21), updated version *
 *  |  S  STRKB02  : KMRS SF set B0 - L=135 (NP=1,NG=3,NS=22)       |  *
 *  |  S  STRKB03  : KMRS SF set B0 - L=160 (NP=1,NG=3,NS=23)       |  *
 *  |  S  STRKB04  : KMRS SF set B0 - L=200 (NP=1,NG=3,NS=24)       |  *
 *  |  S  STRKB05  : KMRS SF set B0 - L=235 (NP=1,NG=3,NS=25)       |  *
 *  |  S  STRKB0O  : KMRS SF set B0 (old version)                   |  *
 *  |  S  STRKBM   : KMRS SF set B- (NP=1,NG=3,NS=18)               |  *
 *  |  S  STRMRS1  : MRS SF set 1 (NP=1,NG=3,NS=1)                  |  *
 *  |  S  STRMRS2  : MRS SF set 2 (NP=1,NG=3,NS=2)                  |  *
 *  |  S  STRMRS3  : MRS SF set 3 (NP=1,NG=3,NS=3)                  |  *
 *  |  S  STROWP1  : OW SF set 1 in the pion (NP=2,NG=1,NS=1)       |  *
 *  |  S  STROWP2  : OW SF set 2 in the pion (NP=2,NG=1,NS=2)       |  *
 *  |  S  STRPRIV  : Dummy subroutine for a new set of SF           |  *
 *  |  S  STRUCB   : MRS SF set B (NP=1,NG=3,NS=4)                  |  *
 *  |  S  STRUCE   : MRS SF set E (NP=1,NG=3,NS=5)                  |  *
 *  |  S  STRUCOR  : OR SF set (NP=1,NG=1,NS=3)                     |  *
 *  |  S  STRUCTF  : to access PDFLIB with output in PDFLIB format  |  *
 *  |  S  STRUCTM  : to access PDFLIB with output in PDFLIB format  |  *
 *  |     TESTGAM  : main program to test PDFLIB for photons        |  *
 *  |     TESTGAMO : main program to test PDFLIB for photons (OLD)  |  *
 *  |     TESTPDF  : main program to test PDFLIB in PDFLIB format   |  *
 *  |     TESTPDFO : main program to test PDFLIB in PDFLIB format (OLD)*
 *  |     TESTPDG  : main program to test PDFLIB in PDG format      |  *
 *  |     TESTPDGO : main program to test PDFLIB in PDG format (OLD)|  *
 *  |     TMAINUPV : main program to show how to use the package    |  *
 *  |_______________________________________________________________|  *
 ******************************* - 29 - ********************************
1***********************************************************************
 *  _________________________________________________________________  *
 *  |  Complete list of subprograms, cont'd                         |  *
 *  |                                                               |  *
 *  |     TMAINUPVO: main program to show how to use the package (OLD) *
 *  |  F  V3       : BEBC SF set to calculate X * VALDIS            |  *
 *  |  S  VALDIS   : BEBC SF set for up + down valence quarks       |  *
 *  |  S  VLAMBD   : Entry in PDZXMT to get Lambda value for MT sets|  *
 *  |  S  WATE32   : defines 32 point Gaussian quadrature routine   |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |                                                               |  *
 *  |_______________________________________________________________|  *
 ******************************* - 30 - ********************************
1***********************************************************************
 *  _________________________________________________________________  *
 *  |  Complete list of subprograms, cont'd                         |  *
 *  |                                                               |  *
 *  |  C  W50510   : IFLPRT - to set print flag                     |  *
 *  |  C  W50511   : NPTYPE,NGROUP,NSET,MODE,NFL,LO,TMAS -          |  *
 *  |                to select structure function set, to define    |  *
 *  |                number of flavours, flag for leading order     |  *
 *  |                and value for the mass of the top quark        |  *
 *  |  C  W50512   : QCDL4,QCDL5 - contains Lambda_4(QCD) and       |  *
 *  |                Lambda_5(QCD) (in GeV)                         |  *
 *  |  C  W50513   : XMIN,XMAX.Q2MIN,Q2MAX - contains minimum       |  *
 *  |                and maximum values of X and Q**2               |  *
 *  |  C  W50514   : WXMIN,WXMAX,WQ2MIN,WQ2MAX,WTXMIN,WTXMAX,WQ2MIN,|  *
 *  |                WQ2MAX - error counting for fractional and     |  *
 *  |                total number of X and Q**2                     |  *
 *  |  C  W50514W  : PDFWGT - to transmit user filled event weight  |  *
 *  |  C  W50515   : IFLSET,IFLSTA - to transmit flag for call to   |  *
 *  |                subroutine PDFSET or PDFSTA                    |  *
 *  |  C  W50516   : FIRST - to force printing name of structure    |  *
 *  |                function set in case of multiple choices       |  *
 *  |  C  W50517   : N6 - contains logical print unit number        |  *
 *  |  C  W50518   : save of COMMONs W50510, W50511, W50512, W50513 |  *
 *  |  C  W50519   : NEWVER - to flag NEW or OLD version of PDFLIB  |  *
 *  |                format                                         |  *
 *  |  C  W505110  : SFNAME - contains NAME (CHARACTER*8 array with |  *
 *  |                dimension SFNAME(NPTYMX,NGRMAX,NSETMX)) for    |  *
 *  |                each set of PDFs                               |  *
 *  |  C  W505120  : NPGSMX(NPTYMX,NGRMAX),NSETFL(NPTYMX,NGRMAX,NSETMX)*
 *  |                - to provide the actual maximum number of stuc-|  *
 *  |                ture function sets per particle type and author|  *
 *  |                group, NSETFL(I,J,K) = 0 for retracted PDF set |  *
 *  |  C  W505121  : NPTYCR(MODEMX),NGROCR(MODEMX),NSETCR(MODEMX)   |  *
 *  |                - to provide the cross reference to the        |  *
 *  |                NPTYPE, NGROUP and NSET value (NEW format) for |  *
 *  |                a given MODE value (OLD format)                |  *
 *  |  C  W505122  : MODECR(NPTYMX,NGRMAX,NSETMX) - to provide the  |  *
 *  |                cross reference to the MODE value (OLD format) |  *
 *  |                for a given NPTYPE, NGROUP and NSET value      |  *
 *  |                (NEW format)                                   |  *
 *  |  P  W5051P0  : LFLPRT - to define print level (D=0)           |  *
 *  |  P  W5051P1  : LPTYPE,LGROUP,LNSET,IFL,LORD,TM - to define    |  *
 *  |                the defaults for a PDF set given by the        |  *
 *  |                particle type, author group and number of set  |  *
 *  |                within an author group, defaults for number of |  *
 *  |                flavours, leading order flag and mass of the   |  *
 *  |                top quark (D=1, 5, 3, 5, 2, 100.D0)            |  *
 *  |  P  W5051P2  : MODEMX - maximum number of structure function  |  *
 *  |                sets in the OLD PDFLIB format                  |  *
 *  |              : NPTYMX,NGRMAX,NSETMX - maximum number of       |  *
 *  |                structure function sets in the NEW PDFLIB      |  *
 *  |                format for particle type, author group number  |  *
 *  |                and number of sets within an author group      |  *
 *  |  P  W5051P2  : MODEMX - to define maximum number of SF sets   |  *
 *  |                in PDFLIB (OLD format) (D=281)                 |  *
 *  |  P  W5051P7  : L6 - to define logical print unit number (D=6) |  *
 *  |  P  W5051P7  : L6 - to define logical print unit number       |  *
 *  |     W50511C  : Comments only - complete list of available SFs |  *
 *  |_______________________________________________________________|  *
 ******************************* - 31 - ********************************
1***********************************************************************
 *                                                                     *
 *  References                                                         *
 *                                                                     *
 *                                                                     *
 *  [1] A.J.Buras and K.J.F. Gaemers : Nucl. Phys. B132 (1978) 249;    *
 *      BEBC WA 59 Collaboration, K. Varnell et al., Z. Phys. C36 (1987) 1.
 *                                                                     *
 *  [2] J.F.Owens and E. Reya : Phys. Rev. D17 (1978) 3003.            *
 *                                                                     *
 *  [3] R. Baier, J. Engels and B. Petersson : Z. Phys. C2 (1979) 265. *
 *                                                                     *
 *  [4] Gl"uck, E. Hoffman and E. Reya : Z. Phys. C13 (1982) 119.      *
 *                                                                     *
 *  [5] D.W. Duke and J.F.Owens : Phys. Rev. D30 (1984) 49.            *
 *                                                                     *
 *  [6] E. Eichten, I Hinchliffe, K. Lane and C. Quigg :               *
 *      Rev. Mod. Phys. 56 (1984) 579 and Rev. Mod. Phys. 58 (1985) 1065.
 *                                                                     *
 *  [7] M. Diemoz, F. Ferroni, E. Longo and G. Martinelli :            *
 *      Z. Phys. C39 (1988) 21.                                        *
 *                                                                     *
 *  [8] J.F.Owens : Phys. Lett. 266B (1991) 126.                       *
 *                                                                     *
 *  [9] A.D. Martin, R.G. Roberts and W.J. Stirling : Phys.Rev. D37 (1988) 1161;
 *      Phys. Lett. 206B (1988) 327 and Mod. Phys. Lett. A4 (1989) 1135.
 *                                                                     *
 * [10] P.N. Harriman, A.D. Martin, R.G. Roberts and W.J. Stirling :   *
 *      Phys. Rev. D42 (1990) 798 and Phys. Lett. 243B (1990) 421.     *
 *                                                                     *
 * [11] J. Kwiecinski, A.D. Martin, R.G. Roberts and W.J. Stirling :   *
 *      Phys. Rev. D42 (1990) 3645;                                    *
 *      A.D. Martin, R.G. Roberts and W.J. Stirling :                  *
 *      Phys. Rev. D43 (1991) 3648.                                    *
 *                                                                     *
 * [12] A.D. Martin, R.G. Roberts and W.J. Stirling :                  *
 *      Phys. Rev. D47 (1993) 867.                                     *
 *                                                                     *
 * [13] J. Morfin and W.K. Tung : Z. Phys. C52 (1991) 13.              *
 *                                                                     *
 * [14] M. Gl"uck, E. Reya and A. Vogt : Z. Phys. C48 (1990) 471;      *
 *      Z. Phys. C53 (1992) 127 and Phys. Lett. 306B (1993) 391.       *
 *                                                                     *
 * [15] P. Aurenche et al. : Phys. Rev. D39 (1989) 3275.               *
 *                                                                     *
 * [16] A.D. Martin, R.G. Roberts and W.J. Stirling :                  *
 *      Phys. Lett. 306B (1993) 145 and Phys. Lett. 309B (1993) 492.   *
 *                                                                     *
 * [17] A.D. Martin, R.G. Roberts and W.J. Stirling :                  *
 *      RAL Preprint, RAL-92-078 (1992) and RAL Preprint, RAL-93-077 (199
 *                                                                     *
 * [18] CTEQ collaboration, J. Botts et al. : Phys. Lett. 304B (1993) 159.
 *                                                                     *
 * [19] E.L. Berger and R. Meng : Phys. Lett. 304B (1993) 318;         *
 *      CERN-TH 6739/92; ANL-HEP-CP-92-108 and                         *
 *      E.L. Berger, R. Meng and J. Qiu : ANL-HEP-CP-92-79.            *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 ******************************* - 32 - ********************************
1***********************************************************************
 *                                                                     *
 * [20] CTEQ collaboration, J. Botts et al. :                          *
 *      to be published                                                *
 *                                                                     *
 *                                                                     *
 * [21] J.F.Owens : Phys. Rev. D30 (1984) 943.                         *
 *                                                                     *
 * [22] P. Aurenche et al. : Phys. Lett. 233B (1989) 517.              *
 *                                                                     *
 * [23] P.J. Sutton, A.D. Martin, R.G. Roberts and W.J. Stirling :     *
 *      Phys. Rev. D45 (1992) 2349.                                    *
 *                                                                     *
 * [24] M. Gl"uck, E. Reya and A. Vogt : Z. Phys. C53 (1992) 651.      *
 *                                                                     *
 * [25] D.W. Duke and J.F.Owens : Phys. Rev. D26 (1982) 1600.          *
 *                                                                     *
 * [26] M. Drees and K. Grassie : Z. Phys. C28 (1985) 451.             *
 *                                                                     *
 * [27] H. Abramowicz, K. Charchula and A. Levy :                      *
 *      Phys. Lett. 269B (1991) 458.                                   *
 *                                                                     *
 * [28] L.E. Gordon and J.K. Storrow : Z. Phys. C56 (1992) 307.        *
 *                                                                     *
 * [29] M. Gl"uck, E. Reya and A. Vogt : Phys. Rev. D46 (1992) 1973    *
 *      and Phys. Rev. D45 (1992) 3986.                                *
 *                                                                     *
 * [30] P. Aurenche et al. : Z. Phys. C56 (1992) 589.                  *
 *                                                                     *
 * [31] G. Altarelli, R.K. Ellis and G. Martinelli :                   *
 *      Nucl. Phys. B143 (1978) 521 and B157 (1979) 461.               *
 *                                                                     *
 * [32] H. Plothow-Besch, 'PDFLIB: a library of all available parton   *
 *      density functions of the nucleon, the pion and the photon      *
 *      and the corresponding alpha(s) calculations',                  *
 *      Comp. Phys. Comm. 75 (1993) 396-416;                           *
 *      H. Plothow-Besch, 'PDFLIB: Structure Functions and alpha(s)    *
 *      Calculations', User's Manual - Version 1.00, W5051 PDFLIB,     *
 *      1991.03.21, CERN-PPE.                                          *
 *                                                                     *
 * [33] Review of Particle Properties, June 1992,                      *
 *      Phys. Rev. D45 (1992), Part. II, p. III.54.                    *
 *                                                                     *
 * [34] BCDMS Collaboration, A.C. Benvenuti et al. : Phys. Lett. 223B (1989) 485
 *      Phys. Lett. 237B (1990) 592 and Phys. Lett. 237B (1990) 599;   *
 *      EMC Collaboration, M. Arneodo et al. : Nucl. Phys. B333 (1990) 1;
 *      NMC Collaboration, D. Allasia et al. : Phys. Lett. 249B (1990) 366;
 *                         P. Amaudruz et al. : Nucl. Phys. B371 (1992) 3.
 *                                                                     *
 * [35] NMC Collaboration, P. Amaudruz et al. : Phys. Lett. 295B (1992) 159.
 *                                                                     *
 * [36] H1 Collaboration, I. Abt et al. : DESY preprint 93-117 (1993); *
 *      ZEUS Collaboration, M. Derrick et al. : DESY preprint 93-110 (1993).
 *                                                                     *
 * [37] JADE Collaboration, W. Bartel et al. : Z. Phys. C24 (1984) 231.*
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 ******************************* - 33 - ********************************
1***********************************************************************
 *                                                                     *
 * [38] PAW, CERN Program Library Entry Q121, Long Write-up;           *
 *      R. Brun, O. Couet, C. Vandoni and P Zanarini;                  *
 *      see also to References therein.                                *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *  List of Tables                                                     *
 *                                                                     *
 *                                                                     *
 *  1  List of available sets of NUCLEON structure functions in        *
 *     PDFLIB version 6.06 ...................................... 35   *
 *                                                                     *
 *                                                                     *
 *  2  List of available sets of PION structure functions in PDFLIB    *
 *     version 6.06 ............................................. 37   *
 *                                                                     *
 *                                                                     *
 *  3  List of available sets of PHOTON structure functions in PDFLIB  *
 *     version 6.06 ............................................. 38   *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 ******************************* - 34 - ********************************
1***********************************************************************
 *                                                                     *
 *                             Table 1                                 *
 *                                                                     *
 *       Available sets of NUCLEON structure functions in PDFLIB       *
 *                                                                     *
 * ------------------------------------------------------------------- *
 *                                                                     *
 * Nptype                                                              *
 *  Ngroup      Q**2  Name of Set   Scheme        Reference            *
 *     Nset     min                                                    *
 *         QCDL4                                                       *
 * ------------------------------------------------------------------- *
 *       0  200 0.3 ->GRV HO        MS_bar  NLL   Default              *
 * ------------------------------------------------------------------- *
 * ------------------------------------------------------------------- *
 * 1  1  1            PRIVATE to user                                  *
 * 1  1  2   45 0.5   BEBC                  LO    [Nucl.Phys. B132 (1978) 249]
 * 1  1  3  500 1.8   OR                    LO    [Phys.Rev. D17 (1978) 3003]
 * 1  1  4  470  5    BEP                   LO    [Z. Phys. C2 (1979) 265]
 * 1  1  5  400  4    GHR                   LO    [Z. Phys. C13 (1982) 119]
 * 1  1  6  200  4    DO Set 1              LO    [Phys.Rev. D30 (1984) 49]
 * 1  1  7  400  4    DO Set 2              LO                         *
 * 1  1  8  200  5    EHLQ Set 1            LO    [Rev.Mod.Phys. 56 (1984) 579;
 * 1  1  9  290  5    EHLQ Set 2            LO     Rev.Mod.Phys. 58 (1985) 1065]
 * 1  1 10  177  4    DO Set 1.1            LO    [Phys.Lett.266B (1991) 3275]
 * ------------------------------------------------------------------- *
 * 1  2  1  200 10    DFLM soft valence     LO    [Z. Phys. C39 (1988) 21
 * 1  2  2  200 10    DFLM hard valence     LO                         *
 * 1  2  3  200 10    DFLM soft gluon       LO                         *
 * 1  2  4  200 10    DFLM hard gluon       LO                         *
 * 1  2  5  200 10    DFLM central av.      LO                         *
 * 1  2  6  300 10    DFLM central av. DIS  NLL                        *
 * 1  2  7  160 10    DFLM 160         DIS  NLL   [unpublished]        *
 * 1  2  8  260 10  ->DFLM 260         DIS  NLL                        *
 * 1  2  9  360 10    DFLM 360         DIS  NLL                        *
 * ------------------------------------------------------------------- *
 * 1  3  1  107  5    MRS 1         MS_bar  NLL   [Phys.Rev. D37 (1988) 1181]
 * 1  3  2  250  5    MRS 2         MS_bar  NLL                        *
 * 1  3  3  178  5    MRS 3         MS_bar  NLL                        *
 * 1  3  4   91  5    MRS E         MS_bar  NLL   [Phys.Lett. 206B (1988) 327]
 * 1  3  5  228  5    MRS B         MS_bar  NLL                        *
 * 1  3  6   91  5    MRS E'        MS_bar  NLL   [Mod.Phys. Lett.A4 (1989) 1135
 * 1  3  7  228  5    MRS B'        MS_bar  NLL                        *
 * 1  3  8  100  5    HMRS E (1.90) MS_bar  NLL   [retracted]          *
 * 1  3  9  190  5    HMRS B (1.90) MS_bar  NLL                        *
 * 1  3 10  100  5    HMRS E (3.90) MS_bar  NLL                        *
 * 1  3 11  190  5    HMRS B (3.90) MS_bar  NLL                        *
 * 1  3 12  100  5    HMRS E+(4.90) MS_bar  NLL   [Phys.Lett. 243B (1990) 421]
 * 1  3 13  100  5    HMRS E-(4.90) MS_bar  NLL                        *
 * 1  3 14  100  5    HMRS E (4.90) MS_bar  NLL   [Phys.Rev. D4 (1990) 798]
 * 1  3 15  190  5    HMRS B (4.90) MS_bar  NLL                        *
 * 1  3 16  100  5    HMRS B (8.90) MS_bar  NLL   [unpublished]        *
 * 1  3 17  300  5    HMRS B (8.90) MS_bar  NLL                        *
 * 1  3 18  190  5    KMRS B-(7.90) MS_bar  NLL   [Phys.Rev. D42 (1990) 3645]
 * 1  3 19  190  5    KMRS B-R2 Shadow (7.90) MS_bar  NLL              *
 * 1  3 20  190  5    KMRS B-R5 Shadow (7.90) MS_bar  NLL              *
 * 1  3 21  190  5    KMRS B0-190 (7.90) MS_bar  NLL (as updated 4.92) *
 *                                                                     *
 ******************************* - 35 - ********************************
1***********************************************************************
 *                                                                     *
 *                          Table 1 (cont'd)                           *
 *                                                                     *
 * ------------------------------------------------------------------- *
 *                                                                     *
 * Nptype                                                              *
 *  Ngroup      Q**2  Name of Set        Scheme      Reference         *
 *     Nset     min                                                    *
 *         QCDL4                                                       *
 * ------------------------------------------------------------------- *
 * 1  3 22  135  5    MRS B0-135 (10.90) MS_bar  NLL [Phys.Rev. D43 (1991) 3648]
 * 1  3 23  160  5    MRS B0-160 (10.90) MS_bar  NLL                   *
 * 1  3 24  200  5    MRS B0-200 (10.90) MS_bar  NLL                   *
 * 1  3 25  235  5    MRS B0-235 (10.90) MS_bar  NLL                   *
 * 1  3 26  215  5    MRS S0     (4.92)  MS_bar  NLL  [retracted]      *
 * 1  3 27  215  5    MRS D0     (4.92)  MS_bar  NLL                   *
 * 1  3 28  215  5    MRS D-     (4.92)  MS_bar  NLL                   *
 * 1  3 29  230  5    MRS S0'    (11.92) MS_bar  NLL  [Phys.Lett. 306B (1993)145
 * 1  3 30  230  5    MRS D0'    (11.92) MS_bar  NLL                   *
 * 1  3 31  230  5    MRS D-'    (11.92) MS_bar  NLL                   *
 * 1  3 32  230  5    MRS S0'    (2.93)  DIS     NLL  [RAL-92-078 (1992)]
 * 1  3 33  230  5    MRS D0'    (2.93)  DIS     NLL                   *
 * 1  3 34  230  5    MRS D-'    (2.93)  DIS     NLL                   *
 * 1  3 35  230  5    MRS (H)    (11.93) MS_bar  NLL                   *
 * 1  3 36  230  5    MRS (H)    (11.93) DIS     NLL                   *
 * 1  3 37  230  5    MRS (A)     (5.94  MS_bar  NLL  [RAL-94-055 (1994)]
 * 1  3 38  230  5    MRS (A) Fit (5.94) MS_bar  NLL                   *
 * 1  3 39  230  5    MRS (Ap)    (2.95  MS_bar  NLL  [RAL-94-104 (1994)]
 * 1  3 40  231  5    MRS (Ap) Fit(2.95) MS_bar  NLL                   *
 * 1  3 41  255  5  ->MRS (G)     (2.95  MS_bar  NLL  [RAL-95-021 (1995)]
 * 1  3 42  255  5    MRS (G) Fit (2.95) MS_bar  NLL                   *
 * ------------------------------------------------------------------- *
 * 1  4  1  212  4    MT S1              DIS     NLL  [Z. Phys. C52 (1991) 13]
 * 1  4  2  194  4    MT B1              DIS     NLL                   *
 * 1  4  3  191  4    MT B2              DIS     NLL                   *
 * 1  4  4  155  4    MT E1              DIS     NLL                   *
 * 1  4  5  237  4    MT 6  (1/2s)       DIS     NLL                   *
 * 1  4  6  212  4    MT S1              MS_bar  NLL                   *
 * 1  4  7  194  4    MT B1              MS_bar  NLL                   *
 * 1  4  8  191  4    MT B2              MS_bar  NLL                   *
 * 1  4  9  155  4    MT E1              MS_bar  NLL                   *
 * 1  4 10  237  4    MT 6  (1/2s)       MS_bar  NLL                   *
 * 1  4 11  144  4    MT LO                      LO                    *
 * 1  4 12  168  4    CTEQ 1L                    LO   [retracted]      *
 * 1  4 13  231  4    CTEQ 1M            MS_bar  NLL                   *
 * 1  4 14  231  4    CTEQ 1MS           MS_bar  NLL                   *
 * 1  4 15  322  4    CTEQ 1ML           MS_bar  NLL                   *
 * 1  4 16  247  4    CTEQ 1D            DIS     NLL                   *
 * 1  4 17  190  4    CTEQ 2L                    LO   [retracted]      *
 * 1  4 18  213  4    CTEQ 2M            MS_bar  NLL                   *
 * 1  4 19  208  4    CTEQ 2MS           MS_bar  NLL                   *
 * 1  4 20  208  4    CTEQ 2MF           MS_bar  NLL                   *
 * 1  4 21  322  4    CTEQ 2ML           MS_bar  NLL                   *
 * 1  4 22  235  4    CTEQ 2D            DIS     NLL                   *
 * 1  4 23  190  4    CTEQ 2pL                   LO   [to be published]*
 * 1  4 24  213  4    CTEQ 2pM           MS_bar  NLL                   *
 * ------------------------------------------------------------------- *
 *                                                                     *
 ******************************* - 36 - ********************************
1***********************************************************************
 *                                                                     *
 *                          Table 1 (cont'd)                           *
 *                                                                     *
 * ------------------------------------------------------------------- *
 *                                                                     *
 * Nptype                                                              *
 *  Ngroup      Q**2  Name of Set        Scheme      Reference         *
 *     Nset     min                                                    *
 *         QCDL4                                                       *
 * ------------------------------------------------------------------- *
 * 1  4 25  208  4    CTEQ 2pMS          MS_bar  NLL                   *
 * 1  4 26  208  4    CTEQ 2pMF          MS_bar  NLL                   *
 * 1  4 27  322  4    CTEQ 2pML          MS_bar  NLL                   *
 * 1  4 28  235  4    CTEQ 2pD           DIS     NLL                   *
 * 1  4 29  177  4    CTEQ 3L                    LO   [MSU-HEP/41024]  *
 * 1  4 30  239  4  ->CTEQ 3M            MS_bar  NLL                   *
 * 1  4 31  247  4    CTEQ 3D            DIS     NLL                   *
 * ------------------------------------------------------------------- *
 * 1  5  1  160 0.2   old GRV HO         MS_bar  NLL  [retracted]      *
 * 1  5  2  220 0.2   old GRV LO                 LO                    *
 * 1  5  3  200 0.3   GRV HO             MS_bar  NLL  [Z. Phys. C53 (1992) 127]
 * 1  5  4  200 0.25  GRV LO                     LO                    *
 * 1  5  5  200 0.4   GRV 94 LO                  LO   [DO-TH 94/24 (1994)]
 * 1  5  6  200 0.4 ->GRV 94 HO          MS_bar  NLL                   *
 * 1  5  7  200 0.4   GRV 94 HO          DIS     NLL                   *
 * ------------------------------------------------------------------- *
 * 1  6  1  230  2    ABFOW              MS_bar  NLL  [Phys.Rev D39 (1989) 3275]
 * ------------------------------------------------------------------- *
 * 1  7  1  254  5    BM A               MS_bar  NLL  [Phys.Lett. 304B (1993)318
 * 1  7  2  254  5    BM B               MS_bar  NLL                   *
 * ------------------------------------------------------------------- *
 *                                                                     *
 *                                                                     *
 *                             Table 2                                 *
 *                                                                     *
 *         Available sets of PION structure functions in PDFLIB        *
 *                                                                     *
 * ------------------------------------------------------------------- *
 *                                                                     *
 * Nptype                                                              *
 *  Ngroup      Q**2  Name of Set   Scheme        Reference            *
 *     Nset     min                                                    *
 *         QCDL4                                                       *
 * ------------------------------------------------------------------- *
 * 2  1  1  200  4    OW-P Set 1            LO    [Phys.Rev. D30 (1984) 943]
 * 2  1  2  400  4    OW-P Set 2            LO                         *
 * ------------------------------------------------------------------- *
 * 2  3  1  190  5    SMRS-P 1      MS_bar  NLL   [Phys.Rev. D45 (1992) 2349]
 * 2  3  2  190  5  ->SMRS-P 2      MS_bar  NLL                        *
 * 2  3  3  190  5    SMRS-P 3      MS_bar  NLL                        *
 * ------------------------------------------------------------------- *
 * 2  5  1  200 0.3 ->GRV-P HO      MS_bar  NLL   [Z.Phys. C53 (1992) 651]
 * 2  5  2  200 0.25  GRV-P LO              LO                         *
 * ------------------------------------------------------------------- *
 * 2  6  1  231  2  ->ABFKW-P 1     MS_bar  NLL   [Phys.Lett.233B (189) 517]
 * 2  6  2  181  2    ABFKW-P 2     MS_bar  NLL                        *
 * 2  6  3  281  2    ABFKW-P 3     MS_bar  NLL                        *
 * ------------------------------------------------------------------- *
 ******************************* - 37 - ********************************
1***********************************************************************
 *                                                                     *
 *                             Table 3                                 *
 *                                                                     *
 *        Available sets of PHOTON structure functions in PDFLIB       *
 *                                                                     *
 * ------------------------------------------------------------------- *
 *                                                                     *
 * Nptype                                                              *
 *  Ngroup      Q**2  Name of Set   Scheme        Reference            *
 *     Nset     min                                                    *
 *         QCDL4                                                       *
 * ------------------------------------------------------------------- *
 * 3  1  1  380 10    DO-G LO              LO     [Phys.Rev. D26 (1982) 1600]
 * 3  1  2  440 10  ->DO-G HO       MS_bar NLL                         *
 * ------------------------------------------------------------------- *
 * 3  2  1  400  1  ->DG-G Set 1           LO     [Z. Phys. C28 (1985) 451]
 * 3  2  2  400  1    DG-G Set 2           LO                          *
 * 3  2  3  400 20    DG-G Set 3           LO                          *
 * 3  2  4  400 200   DG-G Set 4           LO                          *
 * ------------------------------------------------------------------- *
 * 3  3  1  200  5    LAC-G 1              LO     [Phys.Lett.269B (1991) 458]
 * 3  3  2  200  5  ->LAC-G 2              LO                          *
 * 3  3  3  200  5    LAC-G 3              LO                          *
 * ------------------------------------------------------------------- *
 * 3  4  1  200 5.3 ->GS-G HO       MS_bar NLL    [Z.Phys. C56 (1992) 307]
 * 3  4  2  200 5.3   GS-G LO set 1        LO                          *
 * 3  4  3  200 5.3   GS-G LO set 2        LO                          *
 * ------------------------------------------------------------------- *
 * 3  5  1  200 0.3   GRV-G LHO     DIS*)  NLL    [Phys.Rev. D46 (1992) 1973;
 * 3  5  2  200 0.3 ->GRV-G HO      DIS*)  NLL     Phys.Rev. D45 (1992) 3986]
 * 3  5  3  200 0.25  GRV-G LO             LO                          *
 * ------------------------------------------------------------------- *
 * 3  6  1  200  2    ACFGP-G HO    MS_bar NLL    [Z.Phys. C56 (1992) 589]
 * 3  6  2  200  2  ->ACFGP-G HO-mc MS_bar NLL                         *
 * 3  6  3  200  2  ->AFG-G HO      MS_bar NLL                         *
 * ------------------------------------------------------------------- *
 * 3  8  1  400  4    WHIT-G 1             LO     [KEK preprint 93-160]*
 * 3  8  2  400  4  ->WHIT-G 2             LO                          *
 * 3  8  3  400  4    WHIT-G 3             LO                          *
 * 3  8  4  400  4    WHIT-G 4             LO                          *
 * 3  8  5  400  4    WHIT-G 5             LO                          *
 * 3  8  6  400  4    WHIT-G 6             LO                          *
 * ------------------------------------------------------------------- *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *  *) not standard, please consult references.                        *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 ******************************* - 38 - ********************************
1***********************************************************************
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 *                                                                     *
 ******************************* - 39 - ********************************
