]>
Commit | Line | Data |
---|---|---|
4e9e3152 | 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 |