13 kfl(1)=mod(kfa/1000,10)
14 kfl(2)=mod(kfa/100,10)
21 kfl(2)=kfl(2)*(-1)**kfl(2)
22 kfl(3)=-kfl(3)*(-1)**iabs(kfl(2))
23 IF(kflr.EQ.kfl(2))
THEN
25 ELSEIF(kflr.EQ.kfl(3))
THEN
27 ELSEIF(iabs(kflr).EQ.21.AND.
rlu(0).GT.0.5)
THEN
30 ELSEIF(iabs(kflr).EQ.21)
THEN
33 ELSEIF(kflr*kfl(2).GT.0)
THEN
34 CALL
lukfdi(-kflr,kfl(2),kfdump,kflch)
37 CALL
lukfdi(-kflr,kfl(3),kfdump,kflch)
45 100
IF(kflr.EQ.kfl(
j)) nagr=nagr+1
47 ragr=0.00001+(nagr-0.00002)*
rlu(0)
50 IF(kflr.EQ.kfl(
j)) ragr=ragr-1.
51 110
IF(iagr.EQ.0.AND.ragr.LE.0.) iagr=
j
53 iagr=1.00001+2.99998*
rlu(0)
57 IF(iagr.EQ.1.AND.kfl(3).GT.kfl(2)) id1=3
60 IF(mod(kfa,10).EQ.2.AND.kfl(1).EQ.kfl(2))
THEN
61 IF(iagr.NE.3.AND.
rlu(0).GT.0.25) ksp=1
62 ELSEIF(mod(kfa,10).EQ.2.AND.kfl(2).GE.kfl(3))
THEN
63 IF(iagr.NE.1.AND.
rlu(0).GT.0.25) ksp=1
64 ELSEIF(mod(kfa,10).EQ.2)
THEN
66 IF(iagr.NE.1.AND.
rlu(0).GT.0.75) ksp=1
68 kflsp=1000*kfl(id1)+100*kfl(id2)+ksp
71 ELSEIF(nagr.EQ.0.AND.kflr.GT.0)
THEN
72 CALL
lukfdi(-kflr,kfl(iagr),kfdump,kflch)
73 ELSEIF(nagr.EQ.0)
THEN
74 CALL
lukfdi(10000+kflsp,-kflr,kfdump,kflch)