readtbl6.F

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

      Subroutine readtbl6 (Nu)
      Implicit Double Precision (a-h,o-z)
      Character line*80
      parameter(mxx = 96, mxq = 20, mxf = 5)
      parameter(mxpqx = (mxf + 3) * mxq * mxx)
      Common
     > / ctqpar1 / al, xv(0:mxx), tv(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, (tv(i), i =0, nt)

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

      Do 11 iq = 0, nt
         tv(iq) = log(log(tv(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