39
40
41
42#include "implicit_f.inc"
43
44
45
46#include "mvsiz_p.inc"
47
48
49
50#include "scr05_c.inc"
51
52
53
54 INTEGER, INTENT(IN) :: NEL
55 INTEGER, INTENT(IN) :: JHBE
56
58 . r11(mvsiz),r12(mvsiz),r13(mvsiz),
59 . r21(mvsiz),r22(mvsiz),r23(mvsiz),
60 . r31(mvsiz),r32(mvsiz),r33(mvsiz),
61 . x(3,*),off(*)
62 INTEGER
63 . NC1(*),NC2(*),NC3(*),NC4(*),NC5(*),NC6(*),NC7(*),NC8(*)
64
65 DOUBLE PRECISION
66 . XDP(3,*)
67
68
69
70 INTEGER I
71
72
74 . rx(mvsiz) , ry(mvsiz) , rz(mvsiz) ,
75 . sx(mvsiz) , sy(mvsiz) , sz(mvsiz) ,
76 . tx(mvsiz) , ty(mvsiz) , tz(mvsiz)
77 DOUBLE PRECISION
78 . (MVSIZ), XD2(MVSIZ), XD3(MVSIZ), XD4(MVSIZ),
79 . XD5(MVSIZ), XD6(MVSIZ), XD7(MVSIZ), XD8(MVSIZ),
80 . YD1(MVSIZ), YD2(MVSIZ), YD3(MVSIZ), YD4(MVSIZ),
81 . YD5(MVSIZ), YD6(MVSIZ), YD7(MVSIZ), YD8(MVSIZ),
82 . ZD1(MVSIZ), ZD2(MVSIZ), ZD3(MVSIZ), ZD4(MVSIZ),
83 . ZD5(MVSIZ), ZD6(MVSIZ), ZD7(MVSIZ), ZD8(MVSIZ)
84
85
86
87
88
89 IF(iresp==1)THEN
90 DO i=1,nel
91 IF (off(i) <= one ) cycle
92 xd1(i)=xdp(1,nc1(i))
93 yd1(i)=xdp(2,nc1(i))
94 zd1(i)=xdp(3,nc1(i))
95 xd2(i)=xdp(1,nc2(i))
96 yd2(i)=xdp(2,nc2(i))
97 zd2(i)=xdp(3,nc2(i))
98 xd3(i)=xdp(1,nc3(i))
99 yd3(i)=xdp(2,nc3(i))
100 zd3(i)=xdp(3,nc3(i))
101 xd4(i)=xdp(1,nc4(i))
102 yd4(i)=xdp(2,nc4(i))
103 zd4(i)=xdp(3,nc4(i))
104 xd5(i)=xdp(1,nc5(i))
105 yd5(i)=xdp(2,nc5(i))
106 zd5(i)=xdp(3,nc5(i))
107 xd6(i)=xdp(1,nc6(i))
108 yd6(i)=xdp(2,nc6(i))
109 zd6(i)=xdp(3,nc6(i))
110 xd7(i)=xdp(1,nc7(i))
111 yd7(i)=xdp(2,nc7(i))
112 zd7(i)=xdp(3,nc7(i))
113 xd8(i)=xdp(1,nc8(i))
114 yd8(i)=xdp(2,nc8(i))
115 zd8(i)=xdp(3,nc8(i))
116 ENDDO
117 ELSE
118 DO i=1,nel
119 IF (off(i) <= one ) cycle
120 xd1(i)=x(1,nc1(i))
121 yd1(i)=x(2,nc1(i))
122 zd1(i)=x(3,nc1(i))
123 xd2(i)=x(1,nc2(i))
124 yd2(i)=x(2,nc2(i))
125 zd2(i)=x(3,nc2(i))
126 xd3(i)=x(1,nc3(i))
127 yd3(i)=x(2,nc3(i))
128 zd3(i)=x(3,nc3(i))
129 xd4(i)=x(1,nc4(i))
130 yd4(i)=x(2,nc4(i))
131 zd4(i)=x(3,nc4(i))
132 xd5(i)=x(1,nc5(i))
133 yd5(i)=x(2,nc5(i))
134 zd5(i)=x(3,nc5(i))
135 xd6(i)=x(1,nc6(i))
136 yd6(i)=x(2,nc6(i))
137 zd6(i)=x(3,nc6(i))
138 xd7(i)=x(1,nc7(i))
139 yd7(i)=x(2,nc7(i))
140 zd7(i)=x(3,nc7(i))
141 xd8(i)=x(1,nc8(i))
142 yd8(i)=x(2,nc8(i))
143 zd8(i)=x(3,nc8(i))
144 ENDDO
145 ENDIF
146
147
148
150 1 off, xd1, xd2, xd3,
151 2 xd4, xd5, xd6, xd7,
152 3 xd8, yd1, yd2, yd3,
153 4 yd4, yd5, yd6, yd7,
154 5 yd8, zd1, zd2, zd3,
155 6 zd4, zd5, zd6, zd7,
156 7 zd8, rx, ry, rz,
157 8 sx, sy, sz, tx,
158 9 ty, tz, nel)
159
160 IF (jhbe == 14 .OR. jhbe == 24) THEN
162 1 off, rx, ry, rz,
163 2 sx, sy, sz, tx,
164 3 ty, tz, r12, r13,
165 4 r11, r22, r23, r21,
166 5 r32, r33, r31, nel)
167 ELSE
169 1 off, rx, ry, rz,
170 2 sx, sy, sz, tx,
171 3 ty, tz, r11, r12,
172 4 r13, r21, r22, r23,
173 5 r31, r32, r33, nel)
174 ENDIF
175
176 RETURN
subroutine sortho12(off, rx, ry, rz, sx, sy, sz, tx, ty, tz, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, nel)
subroutine srepiso12(off, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, rx, ry, rz, sx, sy, sz, tx, ty, tz, nel)
1.15.0