33
34
35
36#include "implicit_f.inc"
37
38
39
40 INTEGER JFT, JLT,NEL
41
43 . sig(nel,6), g1x(*),g1y(*),g1z(*),
44 . g2x(*),g2y(*),g2z(*),g3x(*),g3y(*),g3z(*),off(*)
45
46
47
48 INTEGER I
49
51 . sx,sy,sz,sign(nel,6)
52
53 DO i=jft,jlt
54 IF (off(i) <=one) cycle
55 sx = sig(i,1)*g1x(i)+sig(i,4)*g1y(i)+sig(i,6)*g1z(i)
56 sy = sig(i,4)*g1x(i)+sig(i,2)*g1y(i)+sig(i,5)*g1z(i)
57 sz = sig(i,6)*g1x(i)+sig(i,5)*g1y(i)+sig(i,3)*g1z(i)
58 sign(i,1) = sx*g1x(i)+sy*g1y(i)+sz*g1z(i)
59 sign(i,4) = sx*g2x(i)+sy*g2y(i)+sz*g2z(i)
60 sign(i,6) = sx*g3x(i)+sy*g3y(i)+sz*g3z(i)
61 sx = sig(i,1)*g2x(i)+sig(i,4)*g2y(i)+sig(i,6)*g2z(i)
62 sy = sig(i,4)*g2x(i)+sig(i,2)*g2y(i)+sig(i,5)*g2z(i)
63 sz = sig(i,6)*g2x(i)+sig(i,5)*g2y(i)+sig(i,3)*g2z(i)
64 sign(i,2) = sx*g2x(i)+sy*g2y(i)+sz*g2z(i)
65 sign(i,5) = sx*g3x(i)+sy*g3y(i)+sz*g3z(i)
66 sx = sig(i,1)*g3x(i)+sig(i,4)*g3y(i)+sig(i,6)*g3z(i)
67 sy = sig(i,4)*g3x(i)+sig(i,2)*g3y(i)+sig(i,5)*g3z(i)
68 sz = sig(i,6)*g3x(i)+sig(i,5)*g3y(i)+sig(i,3)*g3z(i)
69 sign(i,3) = sx*g3x(i)+sy*g3y(i)+sz*g3z(i)
70 sig(i,1) = sign(i,1)
71 sig(i,2) = sign(i,2)
72 sig(i,3) = sign(i,3)
73 sig(i,4) = sign(i,4)
74 sig(i,5) = sign(i,5)
75 sig(i,6) = sign(i,6)
76 ENDDO
77
78 RETURN