setctq5.F

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

      Subroutine setctq5 (Iset)
      Implicit Double Precision (a-h,o-z)
      include "pepadm.inc"
      parameter(isetmax=9)
      Character flnm(isetmax)*100, tablefile*40
      Data (flnm(i), i=1,isetmax)
     > / './pdf/cteq5m.tbl','./pdf/cteq5d.tbl','./pdf/cteq5l.tbl', 
     >   './pdf/cteq5hj.tbl',
     >   './pdf/cteq5hq.tbl','./pdf/cteq5f3.tbl','./pdf/cteq5f4.tbl',
     >   './pdf/cteq5m1.tbl','./pdf/ctq5hq1.tbl'  /
      Data tablefile / 'test.tbl' /
      Data isetold, isetmin, isettest / -987, 1, 911 /

C             If data file not initialized, do so.
      If(iset.ne.isetold) then
        iu= nextun5()
         If (iset .eq. isettest) then
            print* ,'Opening ', tablefile
 21         Open(iu, file=tablefile, status='OLD', err=101)
            goto 22
 101        print*, tablefile, ' cannot be opened '
            print*, 'Please input the .tbl file:'
            Read (*,'(A)') tablefile
            goto 21
 22         Continue
         ElseIf (iset.lt.isetmin .or. iset.gt.isetmax) Then
           print *, 'Invalid Iset number in SetCtq5 :', iset
           stop
         Else
c            Tablefile=Flnm(Iset)
c            Open(IU, File=Tablefile, Status='OLD', Err=100)
            print*, 'We are here yeah'
            print*, cunpol
            Open(iplst(2), file=cunpol, status='OLD', err=100)
            iu=iplst(2)
        Endif
         Call readtbl5(iu)
         Close (iu)
        isetold=iset
      Endif
      Return

 100  print *, ' Data file ', tablefile, ' cannot be opened '
     >//'in SetCtq5!!'
      stop
C                             ********************
      End