OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s20kgeo3.F
Go to the documentation of this file.
1Copyright> OpenRadioss
2Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3Copyright>
4Copyright> This program is free software: you can redistribute it and/or modify
5Copyright> it under the terms of the GNU Affero General Public License as published by
6Copyright> the Free Software Foundation, either version 3 of the License, or
7Copyright> (at your option) any later version.
8Copyright>
9Copyright> This program is distributed in the hope that it will be useful,
10Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12Copyright> GNU Affero General Public License for more details.
13Copyright>
14Copyright> You should have received a copy of the GNU Affero General Public License
15Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16Copyright>
17Copyright>
18Copyright> Commercial Alternative: Altair Radioss Software
19Copyright>
20Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21Copyright> software under a commercial license. Contact Altair to discuss further if the
22Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23
24!||====================================================================
25!|| s20kgeo3 ../engine/source/elements/solid/solide20/s20kgeo3.F
26!||--- called by ------------------------------------------------------
27!|| s20ke3 ../engine/source/elements/solid/solide20/s20ke3.F
28!||--- calls -----------------------------------------------------
29!|| s4keh3 ../engine/source/elements/solid/solide4/s4keh3.F
30!||====================================================================
31 SUBROUTINE s20kgeo3(
32 1 SIG, VOL, PX, PY,
33 2 PZ, K11, K12, K13,
34 3 K14, K15, K16, K17,
35 4 K18, K19, K110, K111,
36 5 K112, K113, K114, K115,
37 6 K116, K117, K118, K119,
38 7 K120, K22, K23, K24,
39 8 K25, K26, K27, K28,
40 9 K29, K210, K211, K212,
41 A K213, K214, K215, K216,
42 B K217, K218, K219, K220,
43 C K33, K34, K35, K36,
44 D K37, K38, K39, K310,
45 E K311, K312, K313, K314,
46 F K315, K316, K317, K318,
47 G K319, K320, K44, K45,
48 H K46, K47, K48, K49,
49 I K410, K411, K412, K413,
50 J K414, K415, K416, K417,
51 K K418, K419, K420, K55,
52 L K56, K57, K58, K59,
53 M K510, K511, K512, K513,
54 N K514, K515, K516, K517,
55 O K518, K519, K520, K66,
56 P K67, K68, K69, K610,
57 Q K611, K612, K613, K614,
58 R K615, K616, K617, K618,
59 S K619, K620, K77, K78,
60 T K79, K710, K711, K712,
61 U K713, K714, K715, K716,
62 V K717, K718, K719, K720,
63 W K88, K89, K810, K811,
64 X K812, K813, K814, K815,
65 Y K816, K817, K818, K819,
66 Z K820, K99, K910, K911,
67 1 K912, K913, K914, K915,
68 2 K916, K917, K918, K919,
69 3 K920, K1010, K1011, K1012,
70 4 K1013, K1014, K1015, K1016,
71 5 K1017, K1018, K1019, K1020,
72 6 K1111, K1112, K1113, K1114,
73 7 K1115, K1116, K1117, K1118,
74 8 K1119, K1120, K1212, K1213,
75 9 K1214, K1215, K1216, K1217,
76 A K1218, K1219, K1220, K1313,
77 B K1314, K1315, K1316, K1317,
78 C K1318, K1319, K1320, K1414,
79 D K1415, K1416, K1417, K1418,
80 E K1419, K1420, K1515, K1516,
81 F K1517, K1518, K1519, K1520,
82 G K1616, K1617, K1618, K1619,
83 H K1620, K1717, K1718, K1719,
84 I K1720, K1818, K1819, K1820,
85 J K1919, K1920, K2020, NEL)
86C-----------------------------------------------
87C I m p l i c i t T y p e s
88C-----------------------------------------------
89#include "implicit_f.inc"
90C-----------------------------------------------
91C G l o b a l P a r a m e t e r s
92C-----------------------------------------------
93#include "mvsiz_p.inc"
94C-----------------------------------------------
95C C o m m o n B l o c k s
96C-----------------------------------------------
97C-----------------------------------------------
98C D u m m y A r g u m e n t s
99C-----------------------------------------------
100 INTEGER NEL
101C REAL
102 my_real
103 . PX(MVSIZ,20),PY(MVSIZ,20),PZ(MVSIZ,20),
104 . VOL(*),SIG(NEL,6)
105 my_real
106 . K11(9,*) ,K12(9,*) ,K13(9,*) ,K14(9,*) ,K15(9,*) ,
107 . K16(9,*) ,K17(9,*) ,K18(9,*) ,K19(9,*) ,K110(9,*) ,
108 . K111(9,*) ,K112(9,*) ,K113(9,*) ,K114(9,*) ,K115(9,*) ,
109 . K116(9,*) ,K117(9,*) ,K118(9,*) ,K119(9,*) ,K120(9,*) ,
110 . K22(9,*) ,K23(9,*) ,K24(9,*) ,K25(9,*) ,K26(9,*) ,
111 . K27(9,*) ,K28(9,*) ,K29(9,*) ,K210(9,*) ,K211(9,*) ,
112 . K212(9,*) ,K213(9,*) ,K214(9,*) ,K215(9,*) ,K216(9,*) ,
113 . K217(9,*) ,K218(9,*) ,K219(9,*) ,K220(9,*) ,K33(9,*) ,
114 . K34(9,*) ,K35(9,*) ,K36(9,*) ,K37(9,*) ,K38(9,*) ,
115 . K39(9,*) ,K310(9,*) ,K311(9,*) ,K312(9,*) ,K313(9,*) ,
116 . K314(9,*) ,K315(9,*) ,K316(9,*) ,K317(9,*) ,K318(9,*) ,
117 . K319(9,*) ,K320(9,*) ,K44(9,*) ,K45(9,*) ,K46(9,*) ,
118 . K47(9,*) ,K48(9,*) ,K49(9,*) ,K410(9,*) ,K411(9,*) ,
119 . K412(9,*) ,K413(9,*) ,K414(9,*) ,K415(9,*) ,K416(9,*) ,
120 . K417(9,*) ,K418(9,*) ,K419(9,*) ,K420(9,*) ,K55(9,*) ,
121 . K56(9,*) ,K57(9,*) ,K58(9,*) ,K59(9,*) ,K510(9,*) ,
122 . K511(9,*) ,K512(9,*) ,K513(9,*) ,K514(9,*) ,K515(9,*) ,
123 . K516(9,*) ,K517(9,*) ,K518(9,*) ,K519(9,*) ,K520(9,*) ,
124 . K66(9,*) ,K67(9,*) ,K68(9,*) ,K69(9,*) ,K610(9,*) ,
125 . K611(9,*) ,K612(9,*) ,K613(9,*) ,K614(9,*) ,K615(9,*) ,
126 . K616(9,*) ,K617(9,*) ,K618(9,*) ,K619(9,*) ,K620(9,*) ,
127 . K77(9,*) ,K78(9,*) ,K79(9,*) ,K710(9,*) ,K711(9,*) ,
128 . K712(9,*) ,K713(9,*) ,K714(9,*) ,K715(9,*) ,K716(9,*) ,
129 . K717(9,*) ,K718(9,*) ,K719(9,*) ,K720(9,*) ,K88(9,*) ,
130 . K89(9,*) ,K810(9,*) ,K811(9,*) ,K812(9,*) ,K813(9,*) ,
131 . K814(9,*) ,K815(9,*) ,K816(9,*) ,K817(9,*) ,K818(9,*) ,
132 . K819(9,*) ,K820(9,*) ,K99(9,*) ,K910(9,*) ,K911(9,*) ,
133 . K912(9,*) ,K913(9,*) ,K914(9,*) ,K915(9,*) ,K916(9,*) ,
134 . K917(9,*) ,K918(9,*) ,K919(9,*) ,K920(9,*) ,K1010(9,*) ,
135 . K1011(9,*) ,K1012(9,*) ,K1013(9,*) ,K1014(9,*) ,K1015(9,*) ,
136 . K1016(9,*) ,K1017(9,*) ,K1018(9,*) ,K1019(9,*) ,K1020(9,*) ,
137 . K1111(9,*) ,K1112(9,*) ,K1113(9,*) ,K1114(9,*) ,K1115(9,*) ,
138 . K1116(9,*) ,K1117(9,*) ,K1118(9,*) ,K1119(9,*) ,K1120(9,*) ,
139 . K1212(9,*) ,K1213(9,*) ,K1214(9,*) ,K1215(9,*) ,K1216(9,*) ,
140 . K1217(9,*) ,K1218(9,*) ,K1219(9,*) ,K1220(9,*) ,K1313(9,*) ,
141 . K1314(9,*) ,K1315(9,*) ,K1316(9,*) ,K1317(9,*) ,K1318(9,*) ,
142 . K1319(9,*) ,K1320(9,*) ,K1414(9,*) ,K1415(9,*) ,K1416(9,*) ,
143 . K1417(9,*) ,K1418(9,*) ,K1419(9,*) ,K1420(9,*) ,K1515(9,*) ,
144 . K1516(9,*) ,K1517(9,*) ,K1518(9,*) ,K1519(9,*) ,K1520(9,*) ,
145 . K1616(9,*) ,K1617(9,*) ,K1618(9,*) ,K1619(9,*) ,K1620(9,*) ,
146 . K1717(9,*) ,K1718(9,*) ,K1719(9,*) ,K1720(9,*) ,K1818(9,*) ,
147 . K1819(9,*) ,K1820(9,*) ,K1919(9,*) ,K1920(9,*) ,K2020(9,*)
148C-----------------------------------------------
149C L o c a l V a r i a b l e s
150C-----------------------------------------------
151 INTEGER I,J
152 my_real
153 . FXX(MVSIZ),FYY(MVSIZ),FZZ(MVSIZ),FXY(MVSIZ),FXZ(MVSIZ),
154 . fyz(mvsiz)
155C-----------------------------------------------
156 DO i=1,nel
157 fxx(i)=sig(i,1)*vol(i)
158 fyy(i)=sig(i,2)*vol(i)
159 fzz(i)=sig(i,3)*vol(i)
160 fxy(i)=sig(i,4)*vol(i)
161 fyz(i)=sig(i,5)*vol(i)
162 fxz(i)=sig(i,6)*vol(i)
163 ENDDO
164C -----------------------------------------------
165 CALL s4keh3(1,nel,
166 . px(1,1),py(1,1),pz(1,1),px(1,1),py(1,1),pz(1,1),
167 . fxx, fyy ,fzz ,fxy, fyz, fxz,k11 )
168 CALL s4keh3(1,nel,
169 . px(1,2),py(1,2),pz(1,2),px(1,2),py(1,2),pz(1,2),
170 . fxx, fyy ,fzz ,fxy, fyz, fxz,k22 )
171 CALL s4keh3(1,nel,
172 . px(1,3),py(1,3),pz(1,3),px(1,3),py(1,3),pz(1,3),
173 . fxx, fyy ,fzz ,fxy, fyz, fxz,k33 )
174 CALL s4keh3(1,nel,
175 . px(1,4),py(1,4),pz(1,4),px(1,4),py(1,4),pz(1,4),
176 . fxx, fyy ,fzz ,fxy, fyz, fxz,k44 )
177 CALL s4keh3(1,nel,
178 . px(1,5),py(1,5),pz(1,5),px(1,5),py(1,5),pz(1,5),
179 . fxx, fyy ,fzz ,fxy, fyz, fxz,k55 )
180 CALL s4keh3(1,nel,
181 . px(1,6),py(1,6),pz(1,6),px(1,6),py(1,6),pz(1,6),
182 . fxx, fyy ,fzz ,fxy, fyz, fxz,k66 )
183 CALL s4keh3(1,nel,
184 . px(1,7),py(1,7),pz(1,7),px(1,7),py(1,7),pz(1,7),
185 . fxx, fyy ,fzz ,fxy, fyz, fxz,k77 )
186 CALL s4keh3(1,nel,
187 . px(1,8),py(1,8),pz(1,8),px(1,8),py(1,8),pz(1,8),
188 . fxx, fyy ,fzz ,fxy, fyz, fxz,k88 )
189 CALL s4keh3(1,nel,
190 . px(1,9),py(1,9),pz(1,9),px(1,9),py(1,9),pz(1,9),
191 . fxx, fyy ,fzz ,fxy, fyz, fxz,k99 )
192 CALL s4keh3(1,nel,
193 . px(1,10),py(1,10),pz(1,10),px(1,10),py(1,10),pz(1,10),
194 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1010 )
195 CALL s4keh3(1,nel,
196 . px(1,11),py(1,11),pz(1,11),px(1,11),py(1,11),pz(1,11),
197 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1111 )
198 CALL s4keh3(1,nel,
199 . px(1,12),py(1,12),pz(1,12),px(1,12),py(1,12),pz(1,12),
200 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1212 )
201 CALL s4keh3(1,nel,
202 . px(1,13),py(1,13),pz(1,13),px(1,13),py(1,13),pz(1,13),
203 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1313 )
204 CALL s4keh3(1,nel,
205 . px(1,14),py(1,14),pz(1,14),px(1,14),py(1,14),pz(1,14),
206 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1414 )
207 CALL s4keh3(1,nel,
208 . px(1,15),py(1,15),pz(1,15),px(1,15),py(1,15),pz(1,15),
209 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1515 )
210 CALL s4keh3(1,nel,
211 . px(1,16),py(1,16),pz(1,16),px(1,16),py(1,16),pz(1,16),
212 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1616 )
213 CALL s4keh3(1,nel,
214 . px(1,17),py(1,17),pz(1,17),px(1,17),py(1,17),pz(1,17),
215 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1717 )
216 CALL s4keh3(1,nel,
217 . px(1,18),py(1,18),pz(1,18),px(1,18),py(1,18),pz(1,18),
218 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1818 )
219 CALL s4keh3(1,nel,
220 . px(1,19),py(1,19),pz(1,19),px(1,19),py(1,19),pz(1,19),
221 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1919 )
222 CALL s4keh3(1,nel,
223 . px(1,20),py(1,20),pz(1,20),px(1,20),py(1,20),pz(1,20),
224 . fxx, fyy ,fzz ,fxy, fyz, fxz,k2020 )
225
226
227C
228 CALL s4keh3(1,nel,
229 . px(1,1),py(1,1),pz(1,1),px(1,2),py(1,2),pz(1,2),
230 . fxx, fyy ,fzz ,fxy, fyz, fxz,k12 )
231 CALL s4keh3(1,nel,
232 . px(1,1),py(1,1),pz(1,1),px(1,3),py(1,3),pz(1,3),
233 . fxx, fyy ,fzz ,fxy, fyz, fxz,k13 )
234 CALL s4keh3(1,nel,
235 . px(1,1),py(1,1),pz(1,1),px(1,4),py(1,4),pz(1,4),
236 . fxx, fyy ,fzz ,fxy, fyz, fxz,k14 )
237 CALL s4keh3(1,nel,
238 . px(1,1),py(1,1),pz(1,1),px(1,5),py(1,5),pz(1,5),
239 . fxx, fyy ,fzz ,fxy, fyz, fxz,k15 )
240 CALL s4keh3(1,nel,
241 . px(1,1),py(1,1),pz(1,1),px(1,6),py(1,6),pz(1,6),
242 . fxx, fyy ,fzz ,fxy, fyz, fxz,k16 )
243 CALL s4keh3(1,nel,
244 . px(1,1),py(1,1),pz(1,1),px(1,7),py(1,7),pz(1,7),
245 . fxx, fyy ,fzz ,fxy, fyz, fxz,k17 )
246 CALL s4keh3(1,nel,
247 . px(1,1),py(1,1),pz(1,1),px(1,8),py(1,8),pz(1,8),
248 . fxx, fyy ,fzz ,fxy, fyz, fxz,k18 )
249 CALL s4keh3(1,nel,
250 . px(1,1),py(1,1),pz(1,1),px(1,9),py(1,9),pz(1,9),
251 . fxx, fyy ,fzz ,fxy, fyz, fxz,k19 )
252 CALL s4keh3(1,nel,
253 . px(1,1),py(1,1),pz(1,1),px(1,10),py(1,10),pz(1,10),
254 . fxx, fyy ,fzz ,fxy, fyz, fxz,k110 )
255 CALL s4keh3(1,nel,
256 . px(1,1),py(1,1),pz(1,1),px(1,11),py(1,11),pz(1,11),
257 . fxx, fyy ,fzz ,fxy, fyz, fxz,k111 )
258 CALL s4keh3(1,nel,
259 . px(1,1),py(1,1),pz(1,1),px(1,12),py(1,12),pz(1,12),
260 . fxx, fyy ,fzz ,fxy, fyz, fxz,k112 )
261 CALL s4keh3(1,nel,
262 . px(1,1),py(1,1),pz(1,1),px(1,13),py(1,13),pz(1,13),
263 . fxx, fyy ,fzz ,fxy, fyz, fxz,k113 )
264 CALL s4keh3(1,nel,
265 . px(1,1),py(1,1),pz(1,1),px(1,14),py(1,14),pz(1,14),
266 . fxx, fyy ,fzz ,fxy, fyz, fxz,k114 )
267 CALL s4keh3(1,nel,
268 . px(1,1),py(1,1),pz(1,1),px(1,15),py(1,15),pz(1,15),
269 . fxx, fyy ,fzz ,fxy, fyz, fxz,k115 )
270 CALL s4keh3(1,nel,
271 . px(1,1),py(1,1),pz(1,1),px(1,16),py(1,16),pz(1,16),
272 . fxx, fyy ,fzz ,fxy, fyz, fxz,k116 )
273 CALL s4keh3(1,nel,
274 . px(1,1),py(1,1),pz(1,1),px(1,17),py(1,17),pz(1,17),
275 . fxx, fyy ,fzz ,fxy, fyz, fxz,k117 )
276 CALL s4keh3(1,nel,
277 . px(1,1),py(1,1),pz(1,1),px(1,18),py(1,18),pz(1,18),
278 . fxx, fyy ,fzz ,fxy, fyz, fxz,k118 )
279 CALL s4keh3(1,nel,
280 . px(1,1),py(1,1),pz(1,1),px(1,19),py(1,19),pz(1,19),
281 . fxx, fyy ,fzz ,fxy, fyz, fxz,k119 )
282 CALL s4keh3(1,nel,
283 . px(1,1),py(1,1),pz(1,1),px(1,20),py(1,20),pz(1,20),
284 . fxx, fyy ,fzz ,fxy, fyz, fxz,k120 )
285
286
287 CALL s4keh3(1,nel,
288 . px(1,2),py(1,2),pz(1,2),px(1,3),py(1,3),pz(1,3),
289 . fxx, fyy ,fzz ,fxy, fyz, fxz,k23 )
290 CALL s4keh3(1,nel,
291 . px(1,2),py(1,2),pz(1,2),px(1,4),py(1,4),pz(1,4),
292 . fxx, fyy ,fzz ,fxy, fyz, fxz,k24 )
293 CALL s4keh3(1,nel,
294 . px(1,2),py(1,2),pz(1,2),px(1,5),py(1,5),pz(1,5),
295 . fxx, fyy ,fzz ,fxy, fyz, fxz,k25 )
296 CALL s4keh3(1,nel,
297 . px(1,2),py(1,2),pz(1,2),px(1,6),py(1,6),pz(1,6),
298 . fxx, fyy ,fzz ,fxy, fyz, fxz,k26 )
299 CALL s4keh3(1,nel,
300 . px(1,2),py(1,2),pz(1,2),px(1,7),py(1,7),pz(1,7),
301 . fxx, fyy ,fzz ,fxy, fyz, fxz,k27 )
302 CALL s4keh3(1,nel,
303 . px(1,2),py(1,2),pz(1,2),px(1,8),py(1,8),pz(1,8),
304 . fxx, fyy ,fzz ,fxy, fyz, fxz,k28 )
305 CALL s4keh3(1,nel,
306 . px(1,2),py(1,2),pz(1,2),px(1,9),py(1,9),pz(1,9),
307 . fxx, fyy ,fzz ,fxy, fyz, fxz,k29 )
308 CALL s4keh3(1,nel,
309 . px(1,2),py(1,2),pz(1,2),px(1,10),py(1,10),pz(1,10),
310 . fxx, fyy ,fzz ,fxy, fyz, fxz,k210 )
311 CALL s4keh3(1,nel,
312 . px(1,2),py(1,2),pz(1,2),px(1,11),py(1,11),pz(1,11),
313 . fxx, fyy ,fzz ,fxy, fyz, fxz,k211 )
314 CALL s4keh3(1,nel,
315 . px(1,2),py(1,2),pz(1,2),px(1,12),py(1,12),pz(1,12),
316 . fxx, fyy ,fzz ,fxy, fyz, fxz,k212 )
317 CALL s4keh3(1,nel,
318 . px(1,2),py(1,2),pz(1,2),px(1,13),py(1,13),pz(1,13),
319 . fxx, fyy ,fzz ,fxy, fyz, fxz,k213 )
320 CALL s4keh3(1,nel,
321 . px(1,2),py(1,2),pz(1,2),px(1,14),py(1,14),pz(1,14),
322 . fxx, fyy ,fzz ,fxy, fyz, fxz,k214 )
323 CALL s4keh3(1,nel,
324 . px(1,2),py(1,2),pz(1,2),px(1,15),py(1,15),pz(1,15),
325 . fxx, fyy ,fzz ,fxy, fyz, fxz,k215 )
326 CALL s4keh3(1,nel,
327 . px(1,2),py(1,2),pz(1,2),px(1,16),py(1,16),pz(1,16),
328 . fxx, fyy ,fzz ,fxy, fyz, fxz,k216 )
329 CALL s4keh3(1,nel,
330 . px(1,2),py(1,2),pz(1,2),px(1,17),py(1,17),pz(1,17),
331 . fxx, fyy ,fzz ,fxy, fyz, fxz,k217 )
332 CALL s4keh3(1,nel,
333 . px(1,2),py(1,2),pz(1,2),px(1,18),py(1,18),pz(1,18),
334 . fxx, fyy ,fzz ,fxy, fyz, fxz,k218 )
335 CALL s4keh3(1,nel,
336 . px(1,2),py(1,2),pz(1,2),px(1,19),py(1,19),pz(1,19),
337 . fxx, fyy ,fzz ,fxy, fyz, fxz,k219 )
338 CALL s4keh3(1,nel,
339 . px(1,2),py(1,2),pz(1,2),px(1,20),py(1,20),pz(1,20),
340 . fxx, fyy ,fzz ,fxy, fyz, fxz,k220 )
341
342
343 CALL s4keh3(1,nel,
344 . px(1,3),py(1,3),pz(1,3),px(1,4),py(1,4),pz(1,4),
345 . fxx, fyy ,fzz ,fxy, fyz, fxz,k34 )
346 CALL s4keh3(1,nel,
347 . px(1,3),py(1,3),pz(1,3),px(1,5),py(1,5),pz(1,5),
348 . fxx, fyy ,fzz ,fxy, fyz, fxz,k35 )
349 CALL s4keh3(1,nel,
350 . px(1,3),py(1,3),pz(1,3),px(1,6),py(1,6),pz(1,6),
351 . fxx, fyy ,fzz ,fxy, fyz, fxz,k36 )
352 CALL s4keh3(1,nel,
353 . px(1,3),py(1,3),pz(1,3),px(1,7),py(1,7),pz(1,7),
354 . fxx, fyy ,fzz ,fxy, fyz, fxz,k37 )
355 CALL s4keh3(1,nel,
356 . px(1,3),py(1,3),pz(1,3),px(1,8),py(1,8),pz(1,8),
357 . fxx, fyy ,fzz ,fxy, fyz, fxz,k38 )
358 CALL s4keh3(1,nel,
359 . px(1,3),py(1,3),pz(1,3),px(1,9),py(1,9),pz(1,9),
360 . fxx, fyy ,fzz ,fxy, fyz, fxz,k39 )
361 CALL s4keh3(1,nel,
362 . px(1,3),py(1,3),pz(1,3),px(1,10),py(1,10),pz(1,10),
363 . fxx, fyy ,fzz ,fxy, fyz, fxz,k310 )
364 CALL s4keh3(1,nel,
365 . px(1,3),py(1,3),pz(1,3),px(1,11),py(1,11),pz(1,11),
366 . fxx, fyy ,fzz ,fxy, fyz, fxz,k311 )
367 CALL s4keh3(1,nel,
368 . px(1,3),py(1,3),pz(1,3),px(1,12),py(1,12),pz(1,12),
369 . fxx, fyy ,fzz ,fxy, fyz, fxz,k312 )
370 CALL s4keh3(1,nel,
371 . px(1,3),py(1,3),pz(1,3),px(1,13),py(1,13),pz(1,13),
372 . fxx, fyy ,fzz ,fxy, fyz, fxz,k313 )
373 CALL s4keh3(1,nel,
374 . px(1,3),py(1,3),pz(1,3),px(1,14),py(1,14),pz(1,14),
375 . fxx, fyy ,fzz ,fxy, fyz, fxz,k314 )
376 CALL s4keh3(1,nel,
377 . px(1,3),py(1,3),pz(1,3),px(1,15),py(1,15),pz(1,15),
378 . fxx, fyy ,fzz ,fxy, fyz, fxz,k315 )
379 CALL s4keh3(1,nel,
380 . px(1,3),py(1,3),pz(1,3),px(1,16),py(1,16),pz(1,16),
381 . fxx, fyy ,fzz ,fxy, fyz, fxz,k316 )
382 CALL s4keh3(1,nel,
383 . px(1,3),py(1,3),pz(1,3),px(1,17),py(1,17),pz(1,17),
384 . fxx, fyy ,fzz ,fxy, fyz, fxz,k317 )
385 CALL s4keh3(1,nel,
386 . px(1,3),py(1,3),pz(1,3),px(1,18),py(1,18),pz(1,18),
387 . fxx, fyy ,fzz ,fxy, fyz, fxz,k318 )
388 CALL s4keh3(1,nel,
389 . px(1,3),py(1,3),pz(1,3),px(1,19),py(1,19),pz(1,19),
390 . fxx, fyy ,fzz ,fxy, fyz, fxz,k319 )
391 CALL s4keh3(1,nel,
392 . px(1,3),py(1,3),pz(1,3),px(1,20),py(1,20),pz(1,20),
393 . fxx, fyy ,fzz ,fxy, fyz, fxz,k320 )
394
395
396 CALL s4keh3(1,nel,
397 . px(1,4),py(1,4),pz(1,4),px(1,5),py(1,5),pz(1,5),
398 . fxx, fyy ,fzz ,fxy, fyz, fxz,k45 )
399 CALL s4keh3(1,nel,
400 . px(1,4),py(1,4),pz(1,4),px(1,6),py(1,6),pz(1,6),
401 . fxx, fyy ,fzz ,fxy, fyz, fxz,k46 )
402 CALL s4keh3(1,nel,
403 . px(1,4),py(1,4),pz(1,4),px(1,7),py(1,7),pz(1,7),
404 . fxx, fyy ,fzz ,fxy, fyz, fxz,k47 )
405 CALL s4keh3(1,nel,
406 . px(1,4),py(1,4),pz(1,4),px(1,8),py(1,8),pz(1,8),
407 . fxx, fyy ,fzz ,fxy, fyz, fxz,k48 )
408 CALL s4keh3(1,nel,
409 . px(1,4),py(1,4),pz(1,4),px(1,9),py(1,9),pz(1,9),
410 . fxx, fyy ,fzz ,fxy, fyz, fxz,k49 )
411 CALL s4keh3(1,nel,
412 . px(1,4),py(1,4),pz(1,4),px(1,10),py(1,10),pz(1,10),
413 . fxx, fyy ,fzz ,fxy, fyz, fxz,k410 )
414 CALL s4keh3(1,nel,
415 . px(1,4),py(1,4),pz(1,4),px(1,11),py(1,11),pz(1,11),
416 . fxx, fyy ,fzz ,fxy, fyz, fxz,k411 )
417 CALL s4keh3(1,nel,
418 . px(1,4),py(1,4),pz(1,4),px(1,12),py(1,12),pz(1,12),
419 . fxx, fyy ,fzz ,fxy, fyz, fxz,k412 )
420 CALL s4keh3(1,nel,
421 . px(1,4),py(1,4),pz(1,4),px(1,13),py(1,13),pz(1,13),
422 . fxx, fyy ,fzz ,fxy, fyz, fxz,k413 )
423 CALL s4keh3(1,nel,
424 . px(1,4),py(1,4),pz(1,4),px(1,14),py(1,14),pz(1,14),
425 . fxx, fyy ,fzz ,fxy, fyz, fxz,k414 )
426 CALL s4keh3(1,nel,
427 . px(1,4),py(1,4),pz(1,4),px(1,15),py(1,15),pz(1,15),
428 . fxx, fyy ,fzz ,fxy, fyz, fxz,k415 )
429 CALL s4keh3(1,nel,
430 . px(1,4),py(1,4),pz(1,4),px(1,16),py(1,16),pz(1,16),
431 . fxx, fyy ,fzz ,fxy, fyz, fxz,k416 )
432 CALL s4keh3(1,nel,
433 . px(1,4),py(1,4),pz(1,4),px(1,17),py(1,17),pz(1,17),
434 . fxx, fyy ,fzz ,fxy, fyz, fxz,k417 )
435 CALL s4keh3(1,nel,
436 . px(1,4),py(1,4),pz(1,4),px(1,18),py(1,18),pz(1,18),
437 . fxx, fyy ,fzz ,fxy, fyz, fxz,k418 )
438 CALL s4keh3(1,nel,
439 . px(1,4),py(1,4),pz(1,4),px(1,19),py(1,19),pz(1,19),
440 . fxx, fyy ,fzz ,fxy, fyz, fxz,k419 )
441 CALL s4keh3(1,nel,
442 . px(1,4),py(1,4),pz(1,4),px(1,20),py(1,20),pz(1,20),
443 . fxx, fyy ,fzz ,fxy, fyz, fxz,k420 )
444
445
446 CALL s4keh3(1,nel,
447 . px(1,5),py(1,5),pz(1,5),px(1,6),py(1,6),pz(1,6),
448 . fxx, fyy ,fzz ,fxy, fyz, fxz,k56 )
449 CALL s4keh3(1,nel,
450 . px(1,5),py(1,5),pz(1,5),px(1,7),py(1,7),pz(1,7),
451 . fxx, fyy ,fzz ,fxy, fyz, fxz,k57 )
452 CALL s4keh3(1,nel,
453 . px(1,5),py(1,5),pz(1,5),px(1,8),py(1,8),pz(1,8),
454 . fxx, fyy ,fzz ,fxy, fyz, fxz,k58 )
455 CALL s4keh3(1,nel,
456 . px(1,5),py(1,5),pz(1,5),px(1,9),py(1,9),pz(1,9),
457 . fxx, fyy ,fzz ,fxy, fyz, fxz,k59 )
458 CALL s4keh3(1,nel,
459 . px(1,5),py(1,5),pz(1,5),px(1,10),py(1,10),pz(1,10),
460 . fxx, fyy ,fzz ,fxy, fyz, fxz,k510 )
461 CALL s4keh3(1,nel,
462 . px(1,5),py(1,5),pz(1,5),px(1,11),py(1,11),pz(1,11),
463 . fxx, fyy ,fzz ,fxy, fyz, fxz,k511 )
464 CALL s4keh3(1,nel,
465 . px(1,5),py(1,5),pz(1,5),px(1,12),py(1,12),pz(1,12),
466 . fxx, fyy ,fzz ,fxy, fyz, fxz,k512 )
467 CALL s4keh3(1,nel,
468 . px(1,5),py(1,5),pz(1,5),px(1,13),py(1,13),pz(1,13),
469 . fxx, fyy ,fzz ,fxy, fyz, fxz,k513 )
470 CALL s4keh3(1,nel,
471 . px(1,5),py(1,5),pz(1,5),px(1,14),py(1,14),pz(1,14),
472 . fxx, fyy ,fzz ,fxy, fyz, fxz,k514 )
473 CALL s4keh3(1,nel,
474 . px(1,5),py(1,5),pz(1,5),px(1,15),py(1,15),pz(1,15),
475 . fxx, fyy ,fzz ,fxy, fyz, fxz,k515 )
476 CALL s4keh3(1,nel,
477 . px(1,5),py(1,5),pz(1,5),px(1,16),py(1,16),pz(1,16),
478 . fxx, fyy ,fzz ,fxy, fyz, fxz,k516 )
479 CALL s4keh3(1,nel,
480 . px(1,5),py(1,5),pz(1,5),px(1,17),py(1,17),pz(1,17),
481 . fxx, fyy ,fzz ,fxy, fyz, fxz,k517 )
482 CALL s4keh3(1,nel,
483 . px(1,5),py(1,5),pz(1,5),px(1,18),py(1,18),pz(1,18),
484 . fxx, fyy ,fzz ,fxy, fyz, fxz,k518 )
485 CALL s4keh3(1,nel,
486 . px(1,5),py(1,5),pz(1,5),px(1,19),py(1,19),pz(1,19),
487 . fxx, fyy ,fzz ,fxy, fyz, fxz,k519 )
488 CALL s4keh3(1,nel,
489 . px(1,5),py(1,5),pz(1,5),px(1,20),py(1,20),pz(1,20),
490 . fxx, fyy ,fzz ,fxy, fyz, fxz,k520 )
491
492
493 CALL s4keh3(1,nel,
494 . px(1,6),py(1,6),pz(1,6),px(1,7),py(1,7),pz(1,7),
495 . fxx, fyy ,fzz ,fxy, fyz, fxz,k67 )
496 CALL s4keh3(1,nel,
497 . px(1,6),py(1,6),pz(1,6),px(1,8),py(1,8),pz(1,8),
498 . fxx, fyy ,fzz ,fxy, fyz, fxz,k68 )
499 CALL s4keh3(1,nel,
500 . px(1,6),py(1,6),pz(1,6),px(1,9),py(1,9),pz(1,9),
501 . fxx, fyy ,fzz ,fxy, fyz, fxz,k69 )
502 CALL s4keh3(1,nel,
503 . px(1,6),py(1,6),pz(1,6),px(1,10),py(1,10),pz(1,10),
504 . fxx, fyy ,fzz ,fxy, fyz, fxz,k610 )
505 CALL s4keh3(1,nel,
506 . px(1,6),py(1,6),pz(1,6),px(1,11),py(1,11),pz(1,11),
507 . fxx, fyy ,fzz ,fxy, fyz, fxz,k611 )
508 CALL s4keh3(1,nel,
509 . px(1,6),py(1,6),pz(1,6),px(1,12),py(1,12),pz(1,12),
510 . fxx, fyy ,fzz ,fxy, fyz, fxz,k612 )
511 CALL s4keh3(1,nel,
512 . px(1,6),py(1,6),pz(1,6),px(1,13),py(1,13),pz(1,13),
513 . fxx, fyy ,fzz ,fxy, fyz, fxz,k613 )
514 CALL s4keh3(1,nel,
515 . px(1,6),py(1,6),pz(1,6),px(1,14),py(1,14),pz(1,14),
516 . fxx, fyy ,fzz ,fxy, fyz, fxz,k614 )
517 CALL s4keh3(1,nel,
518 . px(1,6),py(1,6),pz(1,6),px(1,15),py(1,15),pz(1,15),
519 . fxx, fyy ,fzz ,fxy, fyz, fxz,k615 )
520 CALL s4keh3(1,nel,
521 . px(1,6),py(1,6),pz(1,6),px(1,16),py(1,16),pz(1,16),
522 . fxx, fyy ,fzz ,fxy, fyz, fxz,k616 )
523 CALL s4keh3(1,nel,
524 . px(1,6),py(1,6),pz(1,6),px(1,17),py(1,17),pz(1,17),
525 . fxx, fyy ,fzz ,fxy, fyz, fxz,k617 )
526 CALL s4keh3(1,nel,
527 . px(1,6),py(1,6),pz(1,6),px(1,18),py(1,18),pz(1,18),
528 . fxx, fyy ,fzz ,fxy, fyz, fxz,k618 )
529 CALL s4keh3(1,nel,
530 . px(1,6),py(1,6),pz(1,6),px(1,19),py(1,19),pz(1,19),
531 . fxx, fyy ,fzz ,fxy, fyz, fxz,k619 )
532 CALL s4keh3(1,nel,
533 . px(1,6),py(1,6),pz(1,6),px(1,20),py(1,20),pz(1,20),
534 . fxx, fyy ,fzz ,fxy, fyz, fxz,k620 )
535
536
537 CALL s4keh3(1,nel,
538 . px(1,7),py(1,7),pz(1,7),px(1,8),py(1,8),pz(1,8),
539 . fxx, fyy ,fzz ,fxy, fyz, fxz,k78 )
540 CALL s4keh3(1,nel,
541 . px(1,7),py(1,7),pz(1,7),px(1,9),py(1,9),pz(1,9),
542 . fxx, fyy ,fzz ,fxy, fyz, fxz,k79 )
543 CALL s4keh3(1,nel,
544 . px(1,7),py(1,7),pz(1,7),px(1,10),py(1,10),pz(1,10),
545 . fxx, fyy ,fzz ,fxy, fyz, fxz,k710 )
546 CALL s4keh3(1,nel,
547 . px(1,7),py(1,7),pz(1,7),px(1,11),py(1,11),pz(1,11),
548 . fxx, fyy ,fzz ,fxy, fyz, fxz,k711 )
549 CALL s4keh3(1,nel,
550 . px(1,7),py(1,7),pz(1,7),px(1,12),py(1,12),pz(1,12),
551 . fxx, fyy ,fzz ,fxy, fyz, fxz,k712 )
552 CALL s4keh3(1,nel,
553 . px(1,7),py(1,7),pz(1,7),px(1,13),py(1,13),pz(1,13),
554 . fxx, fyy ,fzz ,fxy, fyz, fxz,k713 )
555 CALL s4keh3(1,nel,
556 . px(1,7),py(1,7),pz(1,7),px(1,14),py(1,14),pz(1,14),
557 . fxx, fyy ,fzz ,fxy, fyz, fxz,k714 )
558 CALL s4keh3(1,nel,
559 . px(1,7),py(1,7),pz(1,7),px(1,15),py(1,15),pz(1,15),
560 . fxx, fyy ,fzz ,fxy, fyz, fxz,k715 )
561 CALL s4keh3(1,nel,
562 . px(1,7),py(1,7),pz(1,7),px(1,16),py(1,16),pz(1,16),
563 . fxx, fyy ,fzz ,fxy, fyz, fxz,k716 )
564 CALL s4keh3(1,nel,
565 . px(1,7),py(1,7),pz(1,7),px(1,17),py(1,17),pz(1,17),
566 . fxx, fyy ,fzz ,fxy, fyz, fxz,k717 )
567 CALL s4keh3(1,nel,
568 . px(1,7),py(1,7),pz(1,7),px(1,18),py(1,18),pz(1,18),
569 . fxx, fyy ,fzz ,fxy, fyz, fxz,k718 )
570 CALL s4keh3(1,nel,
571 . px(1,7),py(1,7),pz(1,7),px(1,19),py(1,19),pz(1,19),
572 . fxx, fyy ,fzz ,fxy, fyz, fxz,k719 )
573 CALL s4keh3(1,nel,
574 . px(1,7),py(1,7),pz(1,7),px(1,20),py(1,20),pz(1,20),
575 . fxx, fyy ,fzz ,fxy, fyz, fxz,k720 )
576
577
578 CALL s4keh3(1,nel,
579 . px(1,8),py(1,8),pz(1,8),px(1,9),py(1,9),pz(1,9),
580 . fxx, fyy ,fzz ,fxy, fyz, fxz,k89 )
581 CALL s4keh3(1,nel,
582 . px(1,8),py(1,8),pz(1,8),px(1,10),py(1,10),pz(1,10),
583 . fxx, fyy ,fzz ,fxy, fyz, fxz,k810 )
584 CALL s4keh3(1,nel,
585 . px(1,8),py(1,8),pz(1,8),px(1,11),py(1,11),pz(1,11),
586 . fxx, fyy ,fzz ,fxy, fyz, fxz,k811 )
587 CALL s4keh3(1,nel,
588 . px(1,8),py(1,8),pz(1,8),px(1,12),py(1,12),pz(1,12),
589 . fxx, fyy ,fzz ,fxy, fyz, fxz,k812 )
590 CALL s4keh3(1,nel,
591 . px(1,8),py(1,8),pz(1,8),px(1,13),py(1,13),pz(1,13),
592 . fxx, fyy ,fzz ,fxy, fyz, fxz,k813 )
593 CALL s4keh3(1,nel,
594 . px(1,8),py(1,8),pz(1,8),px(1,14),py(1,14),pz(1,14),
595 . fxx, fyy ,fzz ,fxy, fyz, fxz,k814 )
596 CALL s4keh3(1,nel,
597 . px(1,8),py(1,8),pz(1,8),px(1,15),py(1,15),pz(1,15),
598 . fxx, fyy ,fzz ,fxy, fyz, fxz,k815 )
599 CALL s4keh3(1,nel,
600 . px(1,8),py(1,8),pz(1,8),px(1,16),py(1,16),pz(1,16),
601 . fxx, fyy ,fzz ,fxy, fyz, fxz,k816 )
602 CALL s4keh3(1,nel,
603 . px(1,8),py(1,8),pz(1,8),px(1,17),py(1,17),pz(1,17),
604 . fxx, fyy ,fzz ,fxy, fyz, fxz,k817 )
605 CALL s4keh3(1,nel,
606 . px(1,8),py(1,8),pz(1,8),px(1,18),py(1,18),pz(1,18),
607 . fxx, fyy ,fzz ,fxy, fyz, fxz,k818 )
608 CALL s4keh3(1,nel,
609 . px(1,8),py(1,8),pz(1,8),px(1,19),py(1,19),pz(1,19),
610 . fxx, fyy ,fzz ,fxy, fyz, fxz,k819 )
611 CALL s4keh3(1,nel,
612 . px(1,8),py(1,8),pz(1,8),px(1,20),py(1,20),pz(1,20),
613 . fxx, fyy ,fzz ,fxy, fyz, fxz,k820 )
614
615
616 CALL s4keh3(1,nel,
617 . px(1,9),py(1,9),pz(1,9),px(1,10),py(1,10),pz(1,10),
618 . fxx, fyy ,fzz ,fxy, fyz, fxz,k910 )
619 CALL s4keh3(1,nel,
620 . px(1,9),py(1,9),pz(1,9),px(1,11),py(1,11),pz(1,11),
621 . fxx, fyy ,fzz ,fxy, fyz, fxz,k911 )
622 CALL s4keh3(1,nel,
623 . px(1,9),py(1,9),pz(1,9),px(1,12),py(1,12),pz(1,12),
624 . fxx, fyy ,fzz ,fxy, fyz, fxz,k912 )
625 CALL s4keh3(1,nel,
626 . px(1,9),py(1,9),pz(1,9),px(1,13),py(1,13),pz(1,13),
627 . fxx, fyy ,fzz ,fxy, fyz, fxz,k913 )
628 CALL s4keh3(1,nel,
629 . px(1,9),py(1,9),pz(1,9),px(1,14),py(1,14),pz(1,14),
630 . fxx, fyy ,fzz ,fxy, fyz, fxz,k914 )
631 CALL s4keh3(1,nel,
632 . px(1,9),py(1,9),pz(1,9),px(1,15),py(1,15),pz(1,15),
633 . fxx, fyy ,fzz ,fxy, fyz, fxz,k915 )
634 CALL s4keh3(1,nel,
635 . px(1,9),py(1,9),pz(1,9),px(1,16),py(1,16),pz(1,16),
636 . fxx, fyy ,fzz ,fxy, fyz, fxz,k916 )
637 CALL s4keh3(1,nel,
638 . px(1,9),py(1,9),pz(1,9),px(1,17),py(1,17),pz(1,17),
639 . fxx, fyy ,fzz ,fxy, fyz, fxz,k917 )
640 CALL s4keh3(1,nel,
641 . px(1,9),py(1,9),pz(1,9),px(1,18),py(1,18),pz(1,18),
642 . fxx, fyy ,fzz ,fxy, fyz, fxz,k918 )
643 CALL s4keh3(1,nel,
644 . px(1,9),py(1,9),pz(1,9),px(1,19),py(1,19),pz(1,19),
645 . fxx, fyy ,fzz ,fxy, fyz, fxz,k919 )
646 CALL s4keh3(1,nel,
647 . px(1,9),py(1,9),pz(1,9),px(1,20),py(1,20),pz(1,20),
648 . fxx, fyy ,fzz ,fxy, fyz, fxz,k920 )
649
650
651 CALL s4keh3(1,nel,
652 . px(1,10),py(1,10),pz(1,10),px(1,11),py(1,11),pz(1,11),
653 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1011 )
654 CALL s4keh3(1,nel,
655 . px(1,10),py(1,10),pz(1,10),px(1,12),py(1,12),pz(1,12),
656 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1012 )
657 CALL s4keh3(1,nel,
658 . px(1,10),py(1,10),pz(1,10),px(1,13),py(1,13),pz(1,13),
659 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1013 )
660 CALL s4keh3(1,nel,
661 . px(1,10),py(1,10),pz(1,10),px(1,14),py(1,14),pz(1,14),
662 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1014 )
663 CALL s4keh3(1,nel,
664 . px(1,10),py(1,10),pz(1,10),px(1,15),py(1,15),pz(1,15),
665 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1015 )
666 CALL s4keh3(1,nel,
667 . px(1,10),py(1,10),pz(1,10),px(1,16),py(1,16),pz(1,16),
668 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1016 )
669 CALL s4keh3(1,nel,
670 . px(1,10),py(1,10),pz(1,10),px(1,17),py(1,17),pz(1,17),
671 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1017 )
672 CALL s4keh3(1,nel,
673 . px(1,10),py(1,10),pz(1,10),px(1,18),py(1,18),pz(1,18),
674 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1018 )
675 CALL s4keh3(1,nel,
676 . px(1,10),py(1,10),pz(1,10),px(1,19),py(1,19),pz(1,19),
677 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1019 )
678 CALL s4keh3(1,nel,
679 . px(1,10),py(1,10),pz(1,10),px(1,20),py(1,20),pz(1,20),
680 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1020 )
681
682
683 CALL s4keh3(1,nel,
684 . px(1,11),py(1,11),pz(1,11),px(1,12),py(1,12),pz(1,12),
685 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1112 )
686 CALL s4keh3(1,nel,
687 . px(1,11),py(1,11),pz(1,11),px(1,13),py(1,13),pz(1,13),
688 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1113 )
689 CALL s4keh3(1,nel,
690 . px(1,11),py(1,11),pz(1,11),px(1,14),py(1,14),pz(1,14),
691 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1114 )
692 CALL s4keh3(1,nel,
693 . px(1,11),py(1,11),pz(1,11),px(1,15),py(1,15),pz(1,15),
694 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1115 )
695 CALL s4keh3(1,nel,
696 . px(1,11),py(1,11),pz(1,11),px(1,16),py(1,16),pz(1,16),
697 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1116 )
698 CALL s4keh3(1,nel,
699 . px(1,11),py(1,11),pz(1,11),px(1,17),py(1,17),pz(1,17),
700 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1117 )
701 CALL s4keh3(1,nel,
702 . px(1,11),py(1,11),pz(1,11),px(1,18),py(1,18),pz(1,18),
703 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1118 )
704 CALL s4keh3(1,nel,
705 . px(1,11),py(1,11),pz(1,11),px(1,19),py(1,19),pz(1,19),
706 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1119 )
707 CALL s4keh3(1,nel,
708 . px(1,11),py(1,11),pz(1,11),px(1,20),py(1,20),pz(1,20),
709 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1120 )
710
711
712 CALL s4keh3(1,nel,
713 . px(1,12),py(1,12),pz(1,12),px(1,13),py(1,13),pz(1,13),
714 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1213 )
715 CALL s4keh3(1,nel,
716 . px(1,12),py(1,12),pz(1,12),px(1,14),py(1,14),pz(1,14),
717 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1214 )
718 CALL s4keh3(1,nel,
719 . px(1,12),py(1,12),pz(1,12),px(1,15),py(1,15),pz(1,15),
720 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1215 )
721 CALL s4keh3(1,nel,
722 . px(1,12),py(1,12),pz(1,12),px(1,16),py(1,16),pz(1,16),
723 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1216 )
724 CALL s4keh3(1,nel,
725 . px(1,12),py(1,12),pz(1,12),px(1,17),py(1,17),pz(1,17),
726 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1217 )
727 CALL s4keh3(1,nel,
728 . px(1,12),py(1,12),pz(1,12),px(1,18),py(1,18),pz(1,18),
729 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1218 )
730 CALL s4keh3(1,nel,
731 . px(1,12),py(1,12),pz(1,12),px(1,19),py(1,19),pz(1,19),
732 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1219 )
733 CALL s4keh3(1,nel,
734 . px(1,12),py(1,12),pz(1,12),px(1,20),py(1,20),pz(1,20),
735 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1220 )
736
737
738 CALL s4keh3(1,nel,
739 . px(1,13),py(1,13),pz(1,13),px(1,14),py(1,14),pz(1,14),
740 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1314 )
741 CALL s4keh3(1,nel,
742 . px(1,13),py(1,13),pz(1,13),px(1,15),py(1,15),pz(1,15),
743 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1315 )
744 CALL s4keh3(1,nel,
745 . px(1,13),py(1,13),pz(1,13),px(1,16),py(1,16),pz(1,16),
746 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1316 )
747 CALL s4keh3(1,nel,
748 . px(1,13),py(1,13),pz(1,13),px(1,17),py(1,17),pz(1,17),
749 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1317 )
750 CALL s4keh3(1,nel,
751 . px(1,13),py(1,13),pz(1,13),px(1,18),py(1,18),pz(1,18),
752 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1318 )
753 CALL s4keh3(1,nel,
754 . px(1,13),py(1,13),pz(1,13),px(1,19),py(1,19),pz(1,19),
755 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1319 )
756 CALL s4keh3(1,nel,
757 . px(1,13),py(1,13),pz(1,13),px(1,20),py(1,20),pz(1,20),
758 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1320 )
759
760
761
762 CALL s4keh3(1,nel,
763 . px(1,14),py(1,14),pz(1,14),px(1,15),py(1,15),pz(1,15),
764 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1415 )
765 CALL s4keh3(1,nel,
766 . px(1,14),py(1,14),pz(1,14),px(1,16),py(1,16),pz(1,16),
767 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1416 )
768 CALL s4keh3(1,nel,
769 . px(1,14),py(1,14),pz(1,14),px(1,17),py(1,17),pz(1,17),
770 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1417 )
771 CALL s4keh3(1,nel,
772 . px(1,14),py(1,14),pz(1,14),px(1,18),py(1,18),pz(1,18),
773 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1418 )
774 CALL s4keh3(1,nel,
775 . px(1,14),py(1,14),pz(1,14),px(1,19),py(1,19),pz(1,19),
776 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1419 )
777 CALL s4keh3(1,nel,
778 . px(1,14),py(1,14),pz(1,14),px(1,20),py(1,20),pz(1,20),
779 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1420 )
780
781
782 CALL s4keh3(1,nel,
783 . px(1,15),py(1,15),pz(1,15),px(1,16),py(1,16),pz(1,16),
784 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1516 )
785 CALL s4keh3(1,nel,
786 . px(1,15),py(1,15),pz(1,15),px(1,17),py(1,17),pz(1,17),
787 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1517 )
788 CALL s4keh3(1,nel,
789 . px(1,15),py(1,15),pz(1,15),px(1,18),py(1,18),pz(1,18),
790 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1518 )
791 CALL s4keh3(1,nel,
792 . px(1,15),py(1,15),pz(1,15),px(1,19),py(1,19),pz(1,19),
793 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1519 )
794 CALL s4keh3(1,nel,
795 . px(1,15),py(1,15),pz(1,15),px(1,20),py(1,20),pz(1,20),
796 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1520 )
797
798
799
800 CALL s4keh3(1,nel,
801 . px(1,16),py(1,16),pz(1,16),px(1,17),py(1,17),pz(1,17),
802 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1617 )
803 CALL s4keh3(1,nel,
804 . px(1,16),py(1,16),pz(1,16),px(1,18),py(1,18),pz(1,18),
805 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1618 )
806 CALL s4keh3(1,nel,
807 . px(1,16),py(1,16),pz(1,16),px(1,19),py(1,19),pz(1,19),
808 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1619 )
809 CALL s4keh3(1,nel,
810 . px(1,16),py(1,16),pz(1,16),px(1,20),py(1,20),pz(1,20),
811 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1620 )
812
813 CALL s4keh3(1,nel,
814 . px(1,17),py(1,17),pz(1,17),px(1,18),py(1,18),pz(1,18),
815 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1718 )
816 CALL s4keh3(1,nel,
817 . px(1,17),py(1,17),pz(1,17),px(1,19),py(1,19),pz(1,19),
818 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1719 )
819 CALL s4keh3(1,nel,
820 . px(1,17),py(1,17),pz(1,17),px(1,20),py(1,20),pz(1,20),
821 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1720 )
822
823
824 CALL s4keh3(1,nel,
825 . px(1,18),py(1,18),pz(1,18),px(1,19),py(1,19),pz(1,19),
826 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1819 )
827 CALL s4keh3(1,nel,
828 . px(1,18),py(1,18),pz(1,18),px(1,20),py(1,20),pz(1,20),
829 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1820 )
830
831
832 CALL s4keh3(1,nel,
833 . px(1,19),py(1,19),pz(1,19),px(1,20),py(1,20),pz(1,20),
834 . fxx, fyy ,fzz ,fxy, fyz, fxz,k1920 )
835
836C
837 RETURN
838 END
s20kgeo3
subroutine s20kgeo3(sig, vol, px, py, pz, k11, k12, k13, k14, k15, k16, k17, k18, k19, k110, k111, k112, k113, k114, k115, k116, k117, k118, k119, k120, k22, k23, k24, k25, k26, k27, k28, k29, k210, k211, k212, k213, k214, k215, k216, k217, k218, k219, k220, k33, k34, k35, k36, k37, k38, k39, k310, k311, k312, k313, k314, k315, k316, k317, k318, k319, k320, k44, k45, k46, k47, k48, k49, k410, k411, k412, k413, k414, k415, k416, k417, k418, k419, k420, k55, k56, k57, k58, k59, k510, k511, k512, k513, k514, k515, k516, k517, k518, k519, k520, k66, k67, k68, k69, k610, k611, k612, k613, k614, k615, k616, k617, k618, k619, k620, k77, k78, k79, k710, k711, k712, k713, k714, k715, k716, k717, k718, k719, k720, k88, k89, k810, k811, k812, k813, k814, k815, k816, k817, k818, k819, k820, k99, k910, k911, k912, k913, k914, k915, k916, k917, k918, k919, k920, k1010, k1011, k1012, k1013, k1014, k1015, k1016, k1017, k1018, k1019, k1020, k1111, k1112, k1113, k1114, k1115, k1116, k1117, k1118, k1119, k1120, k1212, k1213, k1214, k1215, k1216, k1217, k1218, k1219, k1220, k1313, k1314, k1315, k1316, k1317, k1318, k1319, k1320, k1414, k1415, k1416, k1417, k1418, k1419, k1420, k1515, k1516, k1517, k1518, k1519, k1520, k1616, k1617, k1618, k1619, k1620, k1717, k1718, k1719, k1720, k1818, k1819, k1820, k1919, k1920, k2020, nel)
Definition s20kgeo3.F:86
s4keh3
subroutine s4keh3(lft, llt, pxi, pyi, pzi, pxj, pyj, pzj, fxx, fyy, fzz, fxy, fyz, fxz, kij)
Definition s4keh3.F:35