readtbl.F

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

      Subroutine readtbl (Nu)
      Implicit Double Precision (a-h,o-z)
      Character line*80
      parameter(mxx = 105, mxq = 25, mxf = 6)
      parameter(mxpqx = (mxf *2 +2) * mxq * mxx)
      Common 
     > / ctqpar1 / al, xv(0:mxx), ql(0:mxq), upd(mxpqx)
     > / ctqpar2 / nx, nt, nfmx
     > / xqrange / qini, qmax, xmin
     > / qcdtable /  alambda, nfl, iorder
     > / masstbl / amass(6)
      
      Read  (nu, '(A)') line     
      Read  (nu, '(A)') line
      Read  (nu, *) dr, fl, al, (amass(i),i=1,6)
      iorder = nint(dr)
      nfl = nint(fl)
      alambda = al

      Read  (nu, '(A)') line 
      Read  (nu, *) nx,  nt, nfmx

      Read  (nu, '(A)') line
      Read  (nu, *) qini, qmax, (ql(i), i =0, nt)

      Read  (nu, '(A)') line
      Read  (nu, *) xmin, (xv(i), i =0, nx)

      Do 11 iq = 0, nt
         ql(iq) = log(ql(iq) /al)
   11 Continue
C
C                  Since quark = anti-quark for nfl>2 at this stage, 
C                  we Read  out only the non-redundent data points
C     No of flavors = NfMx (sea) + 1 (gluon) + 2 (valence) 

      nblk = (nx+1) * (nt+1)
      npts =  nblk  * (nfmx+3)
      Read  (nu, '(A)') line
      Read  (nu, *, iostat=iret) (upd(i), i=1,npts)

      Return
C                        ****************************
      End