Version update.
[u/mrichter/AliRoot.git] / 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