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)
17 IF(njoin.LT.2) goto 120
21 IF(
i.LE.0.OR.
i.GT.
n) goto 120
22 IF(
k(
i,1).LT.1.OR.
k(
i,1).GT.3) goto 120
25 kq=kchg(kc,2)*isign(1,
k(
i,2))
27 IF(ijn.NE.1.AND.ijn.NE.njoin.AND.kq.NE.2) goto 120
28 IF(kq.NE.2) kqsum=kqsum+kq
29 100
IF(ijn.EQ.1) kqs=kq
30 IF(kqsum.NE.0) goto 120
34 IF(kqs.EQ.2) kcs=int(4.5+
rlu(0))
38 IF(ijn.NE.1) ip=ijoin(ijn-1)
39 IF(ijn.EQ.1) ip=ijoin(njoin)
40 IF(ijn.NE.njoin)
in=ijoin(ijn+1)
41 IF(ijn.EQ.njoin)
in=ijoin(1)
44 IF(ijn.EQ.1.AND.kqs.NE.2)
k(
i,9-kcs)=0
45 110
IF(ijn.EQ.njoin.AND.kqs.NE.2)
k(
i,kcs)=0
50 &
'(LUJOIN:) given entries can not be joined by one string')