LHAPDF 5.2.2 source code.
[u/mrichter/AliRoot.git] / LHAPDF / lhapdf5.2.2 / evolution.f
1       subroutine GetOrderPDF(order)
2       implicit none
3       integer nset,order
4       real*8 Q,Q2
5       nset = 1
6       call GetOrderPDFM(nset,order)
7       return
8 c
9       entry GetRenFac(Q)
10       nset = 1
11       call GetRenFacM(nset,Q)
12       return
13 c
14       entry GetQ2fit(Q2)
15       nset = 1
16       call GetQ2fitM(nset,Q2)
17       return
18       end
19 c
20       subroutine InitEvolve(nset)
21       implicit none
22       character*16 s2
23       real*8 x,Q,Q2,Q2fit,muR
24       integer nset,order,EvlOrd
25 c
26       save EvlOrd,Q2fit,muR
27 *
28       EvlOrd=-1
29       read(1,*) s2,Q2fit,muR
30       if (index(s2,'lo').eq.1) EvlOrd=0
31       if (index(s2,'nlo').eq.1) EvlOrd=1
32       if (index(s2,'nnlo').eq.1) EvlOrd=2
33       if (EvlOrd.lt.0) then
34          write(*,*) 'File description error:'
35          write(*,*) 'Unknown PDF evolution order ',s2
36          stop
37       endif
38       if (muR.ne.1.0d0) then
39          write(*,*) '***********************************************'
40          write(*,*) '* Note than the renormalization scale is      *'
41          write(*,*) '* unequal to the factorization scale for this *'
42          write(*,*) '* particular PDF set.                         *'
43          write(*,*) '* See manual for proper use.                  *'
44          write(*,*) '***********************************************'
45       endif
46 c      print *,'calling readevolve', nset
47       call readevolve(nset)
48       return
49 *
50       entry InitEvolveCode(nset)
51 c      print *,'calling initevolution', nset,Evlord,Q2fit
52       call initevolution(nset,EvlOrd,Q2fit)
53       return
54 *
55       entry GetOrderPDFM(nset,order)
56       order=EvlOrd
57       return
58 *
59       entry GetRenFacM(nset,Q)
60       Q=muR
61       return
62 *
63       entry GetQ2fitM(nset,Q2)
64       Q2=Q2fit
65       return
66 *
67       end