LHAPDF/lhapdf-5.9.1/src/QCDparams.f

   1 ! -*- F90 -*-
   2
   3
   4 subroutine GetLam4(mem,lam4)
   5   implicit none
   6   integer mem,nset
   7   double precision lam4,lam5
   8   nset = 1
   9   call GetLam4M(nset,mem,lam4)
  10   return
  11
  12   entry GetLam5(mem,lam5)
  13   nset = 1
  14   call GetLam5M(nset,mem,lam5)
  15   return
  16 end subroutine GetLam4
  17
  18
  19 subroutine GetXmin(mem,xmin)
  20   implicit none
  21   integer mem,nset
  22   double precision xmin,xmax,q2min,q2max
  23   nset = 1
  24   call GetXminM(nset,mem,xmin)
  25   return
  26
  27   entry GetXmax(mem,xmax)
  28   nset = 1
  29   call GetXmaxM(nset,mem,xmax)
  30   return
  31
  32   entry GetQ2min(mem,q2min)
  33   nset = 1
  34   call GetQ2minM(nset,mem,q2min)
  35   return
  36
  37   entry GetQ2max(mem,q2max)
  38   nset = 1
  39   call GetQ2maxM(nset,mem,q2max)
  40   return
  41
  42   entry GetMinMax(mem,xmin,xmax,q2min,q2max)
  43   nset = 1
  44   call GetMinMaxM(nset,mem,xmin,xmax,q2min,q2max)
  45   return
  46 end subroutine GetXmin
  47
  48
  49 subroutine initQCDparams(nset)
  50   implicit double precision(a-h,o-z)
  51   include 'parmsetup.inc'
  52   double precision parmQCD(nmxset,0:noemax,2),lam4,lam5
  53   integer nset
  54   ! integer iset,imem
  55   ! common/SET/iset,imem
  56   save
  57   read(1,*)nmem,nrep
  58   if(nrep.eq.0) then
  59      do i=0,nmem
  60         read(1,*) (parmQCD(nset,i,j),j=1,2)
  61      enddo
  62   else
  63      read(1,*) (parmQCD(nset,0,j),j=1,2)
  64      do i=1,nmem
  65         do j=1,2
  66            parmQCD(nset,i,j)=parmQCD(nset,0,j)
  67         enddo
  68      enddo
  69   endif
  70   return
  71
  72   entry GetLam4M(nset,mem,lam4)
  73   lam4 = parmQCD(nset,mem,1)
  74   return
  75
  76   entry GetLam5M(nset,mem,lam5)
  77   lam5 = parmQCD(nset,mem,2)
  78   return
  79 end subroutine initQCDparams
  80
  81
  82 subroutine initMinMax(nset)
  83   implicit double precision(a-h,o-z)
  84   include 'parmsetup.inc'
  85   double precision parmXmin(nmxset,0:noemax),xmin
  86   double precision parmXmax(nmxset,0:noemax),xmax
  87   double precision parmQ2min(nmxset,0:noemax),q2min
  88   double precision parmQ2max(nmxset,0:noemax),q2max
  89   integer nset
  90   save
  91   read(1,*)nmem,nrep
  92   if(nrep.eq.0) then
  93      do i=0,nmem
  94         read(1,*) &
  95              parmXmin(nset,i),&
  96              parmXmax(nset,i),&
  97              parmQ2min(nset,i),&
  98              parmQ2max(nset,i)
  99      enddo
 100   else
 101      read(1,*) &
 102           parmXmin(nset,0),&
 103           parmXmax(nset,0),&
 104           parmQ2min(nset,0),&
 105           parmQ2max(nset,0)
 106      do i=1,nmem
 107         parmXmin(nset,i) = parmXmin(nset,0)
 108         parmXmax(nset,i) = parmXmax(nset,0)
 109         parmQ2min(nset,i) = parmQ2min(nset,0)
 110         parmQ2max(nset,i) = parmQ2max(nset,0)
 111      enddo
 112   endif
 113   return
 114
 115   entry GetXminM(nset,mem,xmin)
 116   xmin = parmXmin(nset,mem)
 117   return
 118
 119   entry GetXmaxM(nset,mem,xmax)
 120   xmax = parmXmax(nset,mem)
 121   return
 122
 123   entry GetQ2minM(nset,mem,q2min)
 124   q2min = parmQ2min(nset,mem)
 125   return
 126
 127   entry GetQ2maxM(nset,mem,q2max)
 128   q2max = parmQ2max(nset,mem)
 129   return
 130
 131   entry GetMinMaxM(nset,mem,xmin,xmax,q2min,q2max)
 132   xmin = parmXmin(nset,mem)
 133   xmax = parmXmax(nset,mem)
 134   q2min = parmQ2min(nset,mem)
 135   q2max = parmQ2max(nset,mem)
 136   return
 137 end subroutine initMinMax