setctq4.F

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file setctq4.F

      Subroutine setctq4 (Iset)
      Implicit Double Precision (a-h,o-z)
      INTEGER imxpdf
      parameter(imxpdf=40)
      COMMON /pepadm/cpdfnam(2,imxpdf),ipdfnam(2,imxpdf),
     &       iplst(10),cunpol,cpol
      CHARACTER*256 cpdfnam,cunpol,cpol
      INTEGER iplst,ipdfnam
      SAVE /pepadm/

**************************************************************
*      
*   IPLST(1) =  0 (default) : number of PDF warnings
*   IPLST(2) = 11 (default) : unit -1- for pdf files
*   IPLST(3) = 12 (default) : unit -2- for pdf files      
*
**************************************************************      

      parameter(isetmax=14)
      Character flnm(isetmax)*12, tablefile*40
      Data (flnm(i), i=1,isetmax)
     > / 'cteq4m.tbl', 'cteq4d.tbl', 'cteq4l.tbl'
     > , 'cteq4a1.tbl', 'cteq4a2.tbl', 'cteq4m.tbl', 'cteq4a4.tbl'
     > , 'cteq4a5.tbl', 'cteq4hj.tbl', 'cteq4lq.tbl'
     > , 'cteq4hq.tbl', 'cteq4hq1.tbl', 'cteq4f3.tbl', 'cteq4f4.tbl' /
      Data tablefile / 'test.tbl' /
      Data isetold, isetmin, isettest / -987, 1, 911 /

C             If data file not initialized, do so.
C    Rewrite file name to fit Pepsi distribution format
C
      If(iset.ne.isetold) then
        iu= nextun()
         If (iset .eq. isettest) then
            print* ,'Opening ', tablefile
            Open(iu, file=tablefile, status='OLD', err=101)
         ElseIf (iset.lt.isetmin .or. iset.gt.isetmax) Then
            print *, 'Invalid Iset number in SetCtq4 :', iset
            stop
         Else
*            Tablefile=Flnm(Iset)
*            Tablefile='./pdf/' 
            Open(iplst(2), file=cunpol, status='OLD', err=100)
            iu=iplst(2)
         Endif
         Call readtbl(iu)
         Close (iu)
         isetold=iset
      Endif
      Return

 100  print *, ' Data file ', tablefile, ' cannot be opened '
     > 
      stop
 101  print*, tablefile, ' cannot be opened '
      print*, 'Please input the .tbl file:'
      Read (*,'(A)') tablefile

C                             ********************
      End