setctq6.F

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

      Subroutine setctq6 (Iset)
      Implicit Double Precision (a-h,o-z)
      include "pepadm.inc"
      parameter(isetmax0=3)
      Character flnm(isetmax0)*20, nn*3, tablefile*40
      Data (flnm(i), i=1,isetmax0)
     & /'./pdf/cteq6l.tbl','./pdf/cteq6d.tbl','./pdf/cteq6m.tbl'/
      Data isetold, isetmin0, isetmin1, isetmax1 /-987,1,101,140/

C             If data file not initialized, do so.
      If(iset.ne.isetold) then
         iu= nextun6()
         If (iset.ge.isetmin0 .and. iset.le.isetmax0) Then
C            Tablefile=Flnm(Iset)//'.tbl'
            tablefile=flnm(iset)
         Elseif (iset.ge.isetmin1 .and. iset.le.isetmax1) Then
            write(nn,'(I3)') iset
            tablefile=flnm(1)//nn//'.tbl'
         Else
            print *, 'Invalid Iset number in SetCtq6 :', iset
            stop
         Endif
         Open(iu, file=tablefile, status='OLD', err=100)
 21      Call readtbl6(iu)
         Close (iu)
         isetold=iset
      Endif
      Return

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