LHAPDF/lhapdf5.3.1/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