33 common/ludat1/mstu(200),paru(200),mstj(200),parj(200)
35 common/ludat2/kchg(500,3),pmas(500,4),parf(2000),vckm(4,4)
37 common/pyhipars/mstp(200),parp(200),msti(200),pari(200)
39 dimension xpq(-6:6),xfdflm(9)
40 CHARACTER chdflm(9)*5,header*40
41 DATA chdflm/
'UPVAL',
'DOVAL',
'GLUON',
'QBAR ',
'UBAR ',
'SBAR ',
42 &
'CBAR ',
'BBAR ',
'TBAR '/
43 DATA header/
'Tung evolution package has been invoked'/
50 IF(mstp(51).GE.11.AND.mstp(51).LE.13.AND.mstp(52).LE.1)
THEN
52 q2dflm=
max(10.,
min(1e8,q2))
53 IF(mstp(52).EQ.0) q2dflm=10.
55 IF(mstp(52).EQ.1.AND.j.EQ.9)
THEN
56 q2dflm=q2dflm*(40./pmas(6,1))**2
57 q2dflm=
max(10.,
min(1e8,q2))
65 IF(
x.LT.0.51
e-4.AND.abs(parp(51)-1.).GT.0.01)
THEN
66 cxs=(0.51
e-4/
x)**(parp(51)-1.)
68 110 xfdflm(j)=xfdflm(j)*cxs
71 xpq(1)=xfdflm(2)+xfdflm(5)
72 xpq(2)=xfdflm(1)+xfdflm(5)
90 IF(mstp(52).EQ.4) i1=1
94 IF(mstp(51).GE.11) i2=mstp(51)-3
96 IF(mstp(52).EQ.3) i3=1