7 common/ludat1/mstu(200),paru(200),mstj(200),parj(200)
9 common/ludat2/kchg(500,3),pmas(500,4),parf(2000),vckm(4,4)
13 IF(mstu(111).LE.0)
THEN
23 IF(mstu(115).GE.2) q2eff=
max(q2,paru(114))
26 100
IF(nf.GT.
max(2,mstu(113)))
THEN
27 q2thr=paru(113)*pmas(nf,1)**2
28 IF(q2eff.LT.q2thr)
THEN
30 alam2=alam2*(q2thr/alam2)**(2./(33.-2.*nf))
34 110
IF(nf.LT.
min(8,mstu(114)))
THEN
35 q2thr=paru(113)*pmas(nf+1,1)**2
36 IF(q2eff.GT.q2thr)
THEN
38 alam2=alam2*(alam2/q2thr)**(2./(33.-2.*nf))
42 IF(mstu(115).EQ.1) q2eff=q2eff+alam2
48 IF(mstu(111).EQ.1)
THEN
52 ulalps=paru(2)/(b0*algq)*(1.-b1*
log(algq)/(b0**2*algq))