LHAPDF/lhapdf5.2.2/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 initQCDparams(nset)
  17       implicit real*8(a-h,o-z)
  18       include 'parmsetup.inc'
  19       real*8 parmQCD(nmxset,0:noemax,2),lam4,lam5
  20       integer nset
  21 c      integer iset,imem
  22 c      common/SET/iset,imem
  23       save
  24       read(1,*)nmem,nrep
  25       if(nrep.eq.0) then
  26         do i=0,nmem
  27           read(1,*) (parmQCD(nset,i,j),j=1,2)
  28         enddo
  29       else
  30         read(1,*) (parmQCD(nset,0,j),j=1,2)
  31         do i=1,nmem
  32           do j=1,2
  33             parmQCD(nset,i,j)=parmQCD(nset,0,j)
  34           enddo
  35         enddo
  36       endif
  37       return
  38
  39       entry GetLam4M(nset,mem,lam4)
  40       lam4 = parmQCD(nset,mem,1)
  41       return
  42 c
  43       entry GetLam5M(nset,mem,lam5)
  44       lam5 = parmQCD(nset,mem,2)
  45       return
  46 c
  47       end
  48
  49
  50