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)
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)