32
33
34
35#include "implicit_f.inc"
36
37
38
39#include "mvsiz_p.inc"
40
41
42
43
44
45
46 INTEGER NEL
47
49 . offg(*)
50
51 DOUBLE PRECISION
52 . XD(MVSIZ,8),YD(MVSIZ,8),ZD(MVSIZ,8), SAV(NEL,21)
53
54
55
56 INTEGER I
57
58 DO i=1,nel
59 IF(abs(offg(i)) <= one )THEN
60 sav(i,1)=xd(i,1)-xd(i,8)
61 sav(i,2)=yd(i,1)-yd(i,8)
62 sav(i,3)=zd(i,1)-zd(i,8)
63 sav(i,4)=xd(i,2)-xd(i,8)
64 sav(i,5)=yd(i,2)-yd(i,8)
65 sav(i,6)=zd(i,2)-zd(i,8)
66 sav(i,7)=xd(i,3)-xd(i,8)
67 sav(i,8)=yd(i,3)-yd(i,8)
68 sav(i,9)=zd(i,3)-zd(i,8)
69 sav(i,10)=xd(i,4)-xd(i,8)
70 sav(i,11)=yd(i,4)-yd(i,8)
71 sav(i,12)=zd(i,4)-zd(i,8)
72 sav(i,13)=xd(i,5)-xd(i,8)
73 sav(i,14)=yd(i,5)-yd(i,8)
74 sav(i,15)=zd(i,5)-zd(i,8)
75 sav(i,16)=xd(i,6)-xd(i,8)
76 sav(i,17)=yd(i,6)-yd(i,8)
77 sav(i,18)=zd(i,6)-zd(i,8)
78 sav(i,19)=xd(i,7)-xd(i,8)
79 sav(i,20)=yd(i,7)-yd(i,8)
80 sav(i,21)=zd(i,7)-zd(i,8)
81 ENDIF
82 ENDDO
83
84 RETURN