+++ /dev/null
- subroutine GetLam4(mem,lam4)
- implicit none
- integer mem,nset
- real*8 lam4,lam5
- nset = 1
- call GetLam4M(nset,mem,lam4)
- return
-
- entry GetLam5(mem,lam5)
- nset = 1
- call GetLam5M(nset,mem,lam5)
- return
-
- end
-
- subroutine GetXmin(mem,xmin)
- implicit none
- integer mem,nset
- real*8 xmin,xmax,q2min,q2max
- nset = 1
- call GetXminM(nset,mem,xmin)
- return
-
- entry GetXmax(mem,xmax)
- nset = 1
- call GetXmaxM(nset,mem,xmax)
- return
-
- entry GetQ2min(mem,q2min)
- nset = 1
- call GetQ2minM(nset,mem,q2min)
- return
-
- entry GetQ2max(mem,q2max)
- nset = 1
- call GetQ2maxM(nset,mem,q2max)
- return
-
- entry GetMinMax(mem,xmin,xmax,q2min,q2max)
- nset = 1
- call GetMinMaxM(nset,mem,xmin,xmax,q2min,q2max)
- return
-
- end
-
- subroutine initQCDparams(nset)
- implicit real*8(a-h,o-z)
- include 'parmsetup.inc'
- real*8 parmQCD(nmxset,0:noemax,2),lam4,lam5
- integer nset
-c integer iset,imem
-c common/SET/iset,imem
- save
- read(1,*)nmem,nrep
- if(nrep.eq.0) then
- do i=0,nmem
- read(1,*) (parmQCD(nset,i,j),j=1,2)
- enddo
- else
- read(1,*) (parmQCD(nset,0,j),j=1,2)
- do i=1,nmem
- do j=1,2
- parmQCD(nset,i,j)=parmQCD(nset,0,j)
- enddo
- enddo
- endif
- return
-
- entry GetLam4M(nset,mem,lam4)
- lam4 = parmQCD(nset,mem,1)
- return
-c
- entry GetLam5M(nset,mem,lam5)
- lam5 = parmQCD(nset,mem,2)
- return
-c
- end
-
- subroutine initMinMax(nset)
- implicit real*8(a-h,o-z)
- include 'parmsetup.inc'
- real*8 parmXmin(nmxset,0:noemax),xmin
- real*8 parmXmax(nmxset,0:noemax),xmax
- real*8 parmQ2min(nmxset,0:noemax),q2min
- real*8 parmQ2max(nmxset,0:noemax),q2max
- integer nset
- save
- read(1,*)nmem,nrep
- if(nrep.eq.0) then
- do i=0,nmem
- read(1,*) parmXmin(nset,i),
- + parmXmax(nset,i),
- + parmQ2min(nset,i),
- + parmQ2max(nset,i)
- enddo
- else
- read(1,*) parmXmin(nset,0),
- + parmXmax(nset,0),
- + parmQ2min(nset,0),
- + parmQ2max(nset,0)
- do i=1,nmem
- parmXmin(nset,i) = parmXmin(nset,0)
- parmXmax(nset,i) = parmXmax(nset,0)
- parmQ2min(nset,i) = parmQ2min(nset,0)
- parmQ2max(nset,i) = parmQ2max(nset,0)
- enddo
- endif
- return
-
- entry GetXminM(nset,mem,xmin)
- xmin = parmXmin(nset,mem)
- return
-
- entry GetXmaxM(nset,mem,xmax)
- xmax = parmXmax(nset,mem)
- return
-
- entry GetQ2minM(nset,mem,q2min)
- q2min = parmQ2min(nset,mem)
- return
-
- entry GetQ2maxM(nset,mem,q2max)
- q2max = parmQ2max(nset,mem)
- return
-
- entry GetMinMaxM(nset,mem,xmin,xmax,q2min,q2max)
- xmin = parmXmin(nset,mem)
- xmax = parmXmax(nset,mem)
- q2min = parmQ2min(nset,mem)
- q2max = parmQ2max(nset,mem)
- return
-
- end
-
-
-