8 common/lujets/
n,
k(9000,5),
p(9000,5),
v(9000,5)
10 common/ludat1/mstu(200),paru(200),mstj(200),parj(200)
12 common/ludat2/kchg(500,3),pmas(500,4),parf(2000),vckm(4,4)
14 dimension ns(2),pts(2),pls(2)
17 IF((medit.GE.0.AND.medit.LE.3).OR.medit.EQ.5)
THEN
19 IF(mstu(2).GT.0) imax=mstu(2)
21 DO 110
i=
max(1,mstu(1)),imax
22 IF(
k(
i,1).EQ.0.OR.
k(
i,1).GT.20) goto 110
24 IF(
k(
i,1).GT.10) goto 110
25 ELSEIF(medit.EQ.2)
THEN
26 IF(
k(
i,1).GT.10) goto 110
28 IF(kc.EQ.0.OR.kc.EQ.12.OR.kc.EQ.14.OR.kc.EQ.16.OR.kc.EQ.18)
30 ELSEIF(medit.EQ.3)
THEN
31 IF(
k(
i,1).GT.10) goto 110
34 IF(kchg(kc,2).EQ.0.AND.
luchge(
k(
i,2)).EQ.0) goto 110
35 ELSEIF(medit.EQ.5)
THEN
36 IF(
k(
i,1).EQ.13.OR.
k(
i,1).EQ.14) goto 110
39 IF(
k(
i,1).GE.11.AND.kchg(kc,2).EQ.0) goto 110
53 ELSEIF(medit.GE.11.AND.medit.LE.15)
THEN
56 k(
i,3)=mod(
k(
i,3),mstu(5))
57 IF(medit.EQ.11.AND.
k(
i,1).LT.0) goto 120
58 IF(medit.EQ.12.AND.
k(
i,1).EQ.0) goto 120
59 IF(medit.EQ.13.AND.(
k(
i,1).EQ.11.OR.
k(
i,1).EQ.12.OR.
60 &
k(
i,1).EQ.15).AND.
k(
i,2).NE.94) goto 120
61 IF(medit.EQ.14.AND.(
k(
i,1).EQ.13.OR.
k(
i,1).EQ.14.OR.
62 &
k(
i,2).EQ.94)) goto 120
63 IF(medit.EQ.15.AND.
k(
i,1).GE.21) goto 120
70 IF(
k(
i,1).LE.0.OR.
k(
i,1).GT.20.OR.
k(
i,3)/mstu(5).EQ.0) goto 140
72 130 im=mod(
k(
id,3),mstu(5))
73 IF(medit.EQ.13.AND.im.GT.0.AND.im.LE.
n)
THEN
74 IF((
k(im,1).EQ.11.OR.
k(im,1).EQ.12.OR.
k(im,1).EQ.15).AND.
79 ELSEIF(medit.EQ.14.AND.im.GT.0.AND.im.LE.
n)
THEN
80 IF(
k(im,1).EQ.13.OR.
k(im,1).EQ.14.OR.
k(im,2).EQ.94)
THEN
85 k(
i,3)=mstu(5)*(
k(
i,3)/mstu(5))
86 IF(im.NE.0)
k(
i,3)=
k(
i,3)+
k(im,3)/mstu(5)
87 IF(
k(
i,1).NE.3.AND.
k(
i,1).NE.13.AND.
k(
i,1).NE.14)
THEN
88 IF(
k(
i,4).GT.0.AND.
k(
i,4).LE.mstu(4))
k(
i,4)=
90 IF(
k(
i,5).GT.0.AND.
k(
i,5).LE.mstu(4))
k(
i,5)=
93 kcm=mod(
k(
i,4)/mstu(5),mstu(5))
94 IF(kcm.GT.0.AND.kcm.LE.mstu(4)) kcm=
k(kcm,3)/mstu(5)
95 kcd=mod(
k(
i,4),mstu(5))
96 IF(kcd.GT.0.AND.kcd.LE.mstu(4)) kcd=
k(kcd,3)/mstu(5)
97 k(
i,4)=mstu(5)**2*(
k(
i,4)/mstu(5)**2)+mstu(5)*kcm+kcd
98 kcm=mod(
k(
i,5)/mstu(5),mstu(5))
99 IF(kcm.GT.0.AND.kcm.LE.mstu(4)) kcm=
k(kcm,3)/mstu(5)
100 kcd=mod(
k(
i,5),mstu(5))
101 IF(kcd.GT.0.AND.kcd.LE.mstu(4)) kcd=
k(kcd,3)/mstu(5)
102 k(
i,5)=mstu(5)**2*(
k(
i,5)/mstu(5)**2)+mstu(5)*kcm+kcd
109 IF(
k(
i,3)/mstu(5).EQ.0) goto 160
115 k(
i1,3)=mod(
k(
i1,3),mstu(5))
120 ELSEIF(medit.EQ.21)
THEN
121 IF(2*
n.GE.mstu(4))
THEN
122 CALL
luerrm(11,
'(LUEDIT:) no more memory left in LUJETS')
133 ELSEIF(medit.EQ.22)
THEN
142 ELSEIF(medit.EQ.23)
THEN
147 IF(
k(kh,1).GT.20) kh=0
151 190
IF(
k(
i,1).GT.10.AND.
k(
i,1).LE.20)
k(
i,1)=
k(
i,1)-10
155 ELSEIF(medit.EQ.31.OR.medit.EQ.32)
THEN
156 CALL ludbrb(1,
n+mstu(3),0.,-
ulangl(
p(mstu(61),1),
157 &
p(mstu(61),2)),0d0,0d0,0d0)
158 CALL ludbrb(1,
n+mstu(3),-
ulangl(
p(mstu(61),3),
159 &
p(mstu(61),1)),0.,0d0,0d0,0d0)
160 CALL ludbrb(1,
n+mstu(3),0.,-
ulangl(
p(mstu(61)+1,1),
161 &
p(mstu(61)+1,2)),0d0,0d0,0d0)
162 IF(medit.EQ.31)
RETURN
170 IF(
k(
i,1).LE.0.OR.
k(
i,1).GT.10) goto 220
171 IF(mstu(41).GE.2)
THEN
173 IF(kc.EQ.0.OR.kc.EQ.12.OR.kc.EQ.14.OR.kc.EQ.16.OR.
175 IF(mstu(41).GE.3.AND.kchg(kc,2).EQ.0.AND.
luchge(
k(
i,2)).EQ.0)
180 pts(is)=pts(is)+sqrt(
p(
i,1)**2+
p(
i,2)**2)
182 IF(ns(1)*pts(2)**2.LT.ns(2)*pts(1)**2)
183 & CALL ludbrb(1,
n+mstu(3),paru(1),0.,0d0,0d0,0d0)
187 IF(
p(
i,3).GE.0.) goto 230
188 IF(
k(
i,1).LE.0.OR.
k(
i,1).GT.10) goto 230
189 IF(mstu(41).GE.2)
THEN
191 IF(kc.EQ.0.OR.kc.EQ.12.OR.kc.EQ.14.OR.kc.EQ.16.OR.
193 IF(mstu(41).GE.3.AND.kchg(kc,2).EQ.0.AND.
luchge(
k(
i,2)).EQ.0)
197 pls(is)=pls(is)-
p(
i,3)
199 IF(pls(2).GT.pls(1)) CALL ludbrb(1,
n+mstu(3),0.,paru(1),