57
58
59
60#include "implicit_f.inc"
61
62
63
64
65
66
67 INTEGER, INTENT(IN) :: NEL
68
70 . px1(*), px2(*), px3(*), px4(*),
71 . px5(*), px6(*), px7(*), px8(*),
72 . py1(*), py2(*), py3(*), py4(*),
73 . py5(*), py6(*), py7(*), py8(*),
74 . pz1(*), pz2(*), pz3(*), pz4(*),
75 . pz5(*), pz6(*), pz7(*), pz8(*),
76 . pxc1(*), pxc2(*), pxc3(*), pxc4(*),
77 . pyc1(*), pyc2(*), pyc3(*), pyc4(*),
78 . pzc1(*), pzc2(*), pzc3(*), pzc4(*),
79 . bxy1(*),bxy2(*),bxy3(*),bxy4(*),
80 . bxy5(*),bxy6(*),bxy7(*),bxy8(*),
81 . byx1(*),byx2(*),byx3(*),byx4(*),
82 . byx5(*),byx6(*),byx7(*),byx8(*),
83 . bxz1(*),bxz2(*),bxz3(*),bxz4(*),
84 . bxz5(*),bxz6(*),bxz7(*),bxz8(*),
85 . bzx1(*),bzx2(*),bzx3(*),bzx4(*),
86 . bzx5(*),bzx6(*),bzx7(*),bzx8(*),
87 . byz1(*),byz2(*),byz3(*),byz4(*),
88 . byz5(*),byz6(*),byz7(*),byz8(*),
89 . bzy1(*),bzy2(*),bzy3(*),bzy4(*),
90 . bzy5(*),bzy6(*),bzy7(*),bzy8(*),
91 . bxx1(*),bxx2(*),bxx3(*),bxx4(*),
92 . bxx5(*),bxx6(*),bxx7(*),bxx8(*),
93 . byy1(*),byy2(*),byy3(*),byy4(*),
94 . byy5(*),byy6(*),byy7(*),byy8(*),
95 . bzz1(*),bzz2(*),bzz3(*),bzz4(*),
96 . bzz5(*),bzz6(*),bzz7(*),bzz8(*),nu(*)
97
98
99
100 INTEGER I, J
101
102
104 . bxhi,byhi,bzhi,temp
105
106
107
108
109 temp=half
110 DO i=1,nel
111
112 bxhi=temp*(px1(i)-pxc1(i))
113 byhi=temp*(py1(i)-pyc1(i))
114 bzhi=temp*(pz1(i)-pzc1(i))
115
116
117
118 bxy1(i)=-bxhi
119 bxz1(i)=-bxhi
120 byx1(i)=-byhi
121 byz1(i)=-byhi
122 bzx1(i)=-bzhi
123 bzy1(i)=-bzhi
124 ENDDO
125 DO i=1,nel
126
127 bxhi=temp*(px2(i)-pxc2(i))
128 byhi=temp*(py2(i)-pyc2(i))
129 bzhi=temp*(pz2(i)-pzc2(i))
130
131
132
133 bxy2(i)=-bxhi
134 bxz2(i)=-bxhi
135 byx2(i)=-byhi
136 byz2(i)=-byhi
137 bzx2(i)=-bzhi
138 bzy2(i)=-bzhi
139 ENDDO
140 DO i=1,nel
141
142 bxhi=temp*(px3(i)-pxc3(i))
143 byhi=temp*(py3(i)-pyc3(i))
144 bzhi=temp*(pz3(i)-pzc3(i))
145
146
147
148 bxy3(i)=-bxhi
149 bxz3(i)=-bxhi
150 byx3(i)=-byhi
151 byz3(i)=-byhi
152 bzx3(i)=-bzhi
153 bzy3(i)=-bzhi
154 ENDDO
155 DO i=1,nel
156
157 bxhi=temp*(px4(i)-pxc4(i))
158 byhi=temp*(py4(i)-pyc4(i))
159 bzhi=temp*(pz4(i)-pzc4(i))
160
161
162
163 bxy4(i)=-bxhi
164 bxz4(i)=-bxhi
165 byx4(i)=-byhi
166 byz4(i)=-byhi
167 bzx4(i)=-bzhi
168 bzy4(i)=-bzhi
169 ENDDO
170 DO i=1,nel
171
172 bxhi=temp*(px5(i)+pxc3(i))
173 byhi=temp*(py5(i)+pyc3(i))
174 bzhi=temp*(pz5(i)+pzc3(i))
175
176
177
178 bxy5(i)=-bxhi
179 bxz5(i)=-bxhi
180 byx5(i)=-byhi
181 byz5(i)=-byhi
182 bzx5(i)=-bzhi
183 bzy5(i)=-bzhi
184 ENDDO
185 DO i=1,nel
186
187 bxhi=temp*(px6(i)+pxc4(i))
188 byhi=temp*(py6(i)+pyc4(i))
189 bzhi=temp*(pz6(i)+pzc4(i))
190
191
192
193 bxy6(i)=-bxhi
194 bxz6(i)=-bxhi
195 byx6(i)=-byhi
196 byz6(i)=-byhi
197 bzx6(i)=-bzhi
198 bzy6(i)=-bzhi
199 ENDDO
200 DO i=1,nel
201
202 bxhi=temp*(px7(i)+pxc1(i))
203 byhi=temp*(py7(i)+pyc1(i))
204 bzhi=temp*(pz7(i)+pzc1(i))
205
206
207
208 bxy7(i)=-bxhi
209 bxz7(i)=-bxhi
210 byx7(i)=-byhi
211 byz7(i)=-byhi
212 bzx7(i)=-bzhi
213 bzy7(i)=-bzhi
214 ENDDO
215 DO i=1,nel
216
217 bxhi=temp*(px8(i)+pxc2(i))
218 byhi=temp*(py8(i)+pyc2(i))
219 bzhi=temp*(pz8(i)+pzc2(i))
220
221
222
223 bxy8(i)=-bxhi
224 bxz8(i)=-bxhi
225 byx8(i)=-byhi
226 byz8(i)=-byhi
227 bzx8(i)=-bzhi
228 bzy8(i)=-bzhi
229 ENDDO
230 DO i=1,nel
231 bxx1(i) = px1(i)
232 byy1(i) = py1(i)
233 bzz1(i) = pz1(i)
234 bxx2(i) = px2(i)
235 byy2(i) = py2(i)
236 bzz2(i) = pz2(i)
237 bxx3(i) = px3(i)
238 byy3(i) = py3(i)
239 bzz3(i) = pz3(i)
240 bxx4(i) = px4(i)
241 byy4(i) = py4(i)
242 bzz4(i) = pz4(i)
243 bxx5(i) = px5(i)
244 byy5(i) = py5(i)
245 bzz5(i) = pz5(i)
246 bxx6(i) = px6(i)
247 byy6(i) = py6(i)
248 bzz6(i) = pz6(i)
249 bxx7(i) = px7(i)
250 byy7(i) = py7(i)
251 bzz7(i) = pz7(i)
252 bxx8(i) = px8(i)
253 byy8(i) = py8(i)
254 bzz8(i) = pz8(i)
255 ENDDO
256
257 RETURN