LHAPDF 5.2.2 source code.
[u/mrichter/AliRoot.git] / 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