LHAPDF/lhapdf5.3.1/QCDparams.f

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