4 SUBROUTINE luzdis(KFL1,KFL2,PR,Z)
7 common/ludat1/mstu(200),paru(200),mstj(200),parj(200)
14 IF(kfla.GE.10) kflh=
mod(kfla/1000,10)
17 IF(mstj(11).EQ.1.OR.(mstj(11).EQ.3.AND.kflh.LE.3))
THEN
19 IF(mstj(91).EQ.1)
fa=parj(43)
20 IF(kflb.GE.10)
fa=
fa+parj(45)
22 IF(mstj(91).EQ.1)
fb=parj(44)*pr
24 IF(kfla.GE.10) fc=fc-parj(45)
25 IF(kflb.GE.10) fc=fc+parj(45)
27 IF(abs(fc-1.).GT.0.01) mc=2
33 IF(fc.GT.
fb) zmax=
fb/fc
34 ELSEIF(abs(fc-
fa).LT.0.01)
THEN
40 IF(zmax.GT.0.99.AND.
fb.GT.100.) zmax=1.-
fa/
fb
52 fint=1.+(1.-1./zdivc)/(fc-1.)
54 ELSEIF(zmax.GT.0.85.AND.
fb.GT.1.)
THEN
57 zdiv=fscb-1./zmax-(fc/
fb)*
log(zmax*0.5*(fscb+fc/
fb))
58 IF(ma.GE.2) zdiv=zdiv+(
fa/
fb)*
log(1.-zmax)
59 zdiv=
min(zmax,
max(0.,zdiv))
67 IF(fint*
rlu(0).LE.1.)
THEN
73 z=1./(zdivc+
z*(1.-zdivc))**(1./(1.-fc))
76 ELSEIF(mmax.EQ.3)
THEN
77 IF(fint*
rlu(0).LE.1.)
THEN
86 IF(
z.LE.
fb/(50.+
fb).OR.
z.GE.1.) goto 100
87 fval=(zmax/
z)**fc*
exp(
fb*(1./zmax-1./
z))
88 IF(ma.GE.2) fval=((1.-
z)/(1.-zmax))**
fa*fval
89 IF(fval.LT.
rlu(0)*fpre) goto 100
93 fc=parj(50+
max(1,kflh))
94 IF(mstj(91).EQ.1) fc=parj(59)
96 IF(fc.GE.0..AND.fc.LE.1.)
THEN
97 IF(fc.GT.
rlu(0))
z=1.-
z**(1./3.)
98 ELSEIF(fc.GT.-1.)
THEN
99 IF(-4.*fc*
z*(1.-
z)**2.LT.
rlu(0)*((1.-
z)**2-fc*
z)**2) goto 110
101 IF(fc.GT.0.)
z=1.-
z**(1./fc)
102 IF(fc.LT.0.)
z=
z**(-1./fc)