7 dimension psc1(5),psc2(5)
16 IF(jp.EQ.0 .OR. jt.EQ.0) go to 25
26 pp(jp,6)=pp(jp,6)+dpp1/2.0
27 pp(jp,7)=pp(jp,7)+dpp2/2.0
28 pp(jp,8)=pp(jp,8)+dpp1/2.0
29 pp(jp,9)=pp(jp,9)+dpp2/2.0
30 pt(jt,6)=
pt(jt,6)+dpt1/2.0
31 pt(jt,7)=
pt(jt,7)+dpt2/2.0
32 pt(jt,8)=
pt(jt,8)+dpt1/2.0
33 pt(jt,9)=
pt(jt,9)+dpt2/2.0
44 25
IF(jp.EQ.0) go to 45
45 pabs=
sqrt(pp(jp,1)**2+pp(jp,2)**2+pp(jp,3)**2)
73 pp(jp,6)=pp(jp,6)+dpp1/2.0
74 pp(jp,7)=pp(jp,7)+dpp2/2.0
75 pp(jp,8)=pp(jp,8)+dpp1/2.0
76 pp(jp,9)=pp(jp,9)+dpp2/2.0
77 pp(i,6)=pp(i,6)+dpt1/2.0
78 pp(i,7)=pp(i,7)+dpt2/2.0
79 pp(i,8)=pp(i,8)+dpt1/2.0
80 pp(i,9)=pp(i,9)+dpt2/2.0
88 45
IF(jt.EQ.0) go to 80
89 50 pabs=
sqrt(
pt(jt,1)**2+
pt(jt,2)**2+
pt(jt,3)**2)
113 dpp1=psc1(1)-
pt(jt,1)
114 dpp2=psc1(2)-
pt(jt,2)
117 pt(jt,6)=
pt(jt,6)+dpp1/2.0
118 pt(jt,7)=
pt(jt,7)+dpp2/2.0
119 pt(jt,8)=
pt(jt,8)+dpp1/2.0
120 pt(jt,9)=
pt(jt,9)+dpp2/2.0
121 pt(i,6)=
pt(i,6)+dpt1/2.0
122 pt(i,7)=
pt(i,7)+dpt2/2.0
123 pt(i,8)=
pt(i,8)+dpt1/2.0
124 pt(i,9)=
pt(i,9)+dpt2/2.0