OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s10cumg3.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!|| s10cumg3 ../engine/source/elements/solid/solide10/s10cumg3.F
25!||--- called by ------------------------------------------------------
26!|| s10ke3 ../engine/source/elements/solid/solide10/s10ke3.F
27!||--- calls -----------------------------------------------------
28!|| s8zkebg1 ../engine/source/elements/solid/solide8z/s8zkebg1.F
29!|| s8zked3 ../engine/source/elements/solid/solide8z/s8zked3.F
30!|| s8zkeg3 ../engine/source/elements/solid/solide8z/s8zkeg3.F
31!|| s8zkeg33 ../engine/source/elements/solid/solide8z/s8zkeg33.F
32!||====================================================================
33 SUBROUTINE s10cumg3(
34 1 PX, PY, PZ, DD,
35 2 GG, DG, G33, IKSUP,
36 3 K11, K12, K13, K14,
37 4 K15, K16, K17, K18,
38 5 K19, K10, K22, K23,
39 6 K24, K25, K26, K27,
40 7 K28, K29, K20, K33,
41 8 K34, K35, K36, K37,
42 9 K38, K39, K30, K44,
43 A K45, K46, K47, K48,
44 B K49, K40, K55, K56,
45 C K57, K58, K59, K50,
46 D K66, K67, K68, K69,
47 E K60, K77, K78, K79,
48 F K70, K88, K89, K80,
49 G K99, K90, K00, NEL)
50C-----------------------------------------------
51C I m p l i c i t T y p e s
52C-----------------------------------------------
53#include "implicit_f.inc"
54C-----------------------------------------------
55C G l o b a l P a r a m e t e r s
56C-----------------------------------------------
57#include "mvsiz_p.inc"
58C-----------------------------------------------
59C C o m m o n B l o c k s
60C-----------------------------------------------
61C-----------------------------------------------
62C D u m m y A r g u m e n t s
63C-----------------------------------------------
64 INTEGER, INTENT(IN) :: NEL
65 INTEGER IKSUP
66 my_real
67 . PX(MVSIZ,10),PY(MVSIZ,10),PZ(MVSIZ,10),
68 . DD(3,3,*),GG(*),DG(9,*) ,G33(9,*)
69 my_real
70 . K11(9,*) ,K12(9,*) ,K13(9,*) ,K14(9,*) ,K15(9,*) ,
71 . K16(9,*) ,K17(9,*) ,K18(9,*) ,K19(9,*) ,K10(9,*) ,
72 . K22(9,*) ,K23(9,*) ,K24(9,*) ,K25(9,*) ,K26(9,*) ,
73 . K27(9,*) ,K28(9,*) ,K29(9,*) ,K20(9,*) ,K33(9,*) ,
74 . K34(9,*) ,K35(9,*) ,K36(9,*) ,K37(9,*) ,K38(9,*) ,
75 . K39(9,*) ,K30(9,*) ,K44(9,*) ,K45(9,*) ,K46(9,*) ,
76 . K47(9,*) ,K48(9,*) ,K49(9,*) ,K40(9,*) ,K55(9,*) ,
77 . K56(9,*) ,K57(9,*) ,K58(9,*) ,K59(9,*) ,K50(9,*) ,
78 . K66(9,*) ,K67(9,*) ,K68(9,*) ,K69(9,*) ,K60(9,*) ,
79 . K77(9,*) ,K78(9,*) ,K79(9,*) ,K70(9,*) ,K88(9,*) ,
80 . k89(9,*) ,k80(9,*) ,k99(9,*) ,k90(9,*) ,k00(9,*)
81C-----------------------------------------------
82C L o c a l V a r i a b l e s
83C-----------------------------------------------
84 INTEGER I, J,IS
85C-----------------------------------------------
86C -------symm(diag)------
87 is =1
88C 12
89 CALL s8zked3(1,nel,
90 . px(1,1),py(1,1),pz(1,1),px(1,1),py(1,1),pz(1,1),
91 . dd, k11,is)
92 CALL s8zked3(1,nel,
93 . px(1,2),py(1,2),pz(1,2),px(1,2),py(1,2),pz(1,2),
94 . dd, k22,is)
95 CALL s8zked3(1,nel,
96 . px(1,3),py(1,3),pz(1,3),px(1,3),py(1,3),pz(1,3),
97 . dd, k33,is)
98 CALL s8zked3(1,nel,
99 . px(1,4),py(1,4),pz(1,4),px(1,4),py(1,4),pz(1,4),
100 . dd, k44,is)
101 CALL s8zked3(1,nel,
102 . px(1,5),py(1,5),pz(1,5),px(1,5),py(1,5),pz(1,5),
103 . dd, k55,is)
104 CALL s8zked3(1,nel,
105 . px(1,6),py(1,6),pz(1,6),px(1,6),py(1,6),pz(1,6),
106 . dd, k66,is)
107 CALL s8zked3(1,nel,
108 . px(1,7),py(1,7),pz(1,7),px(1,7),py(1,7),pz(1,7),
109 . dd, k77,is)
110 CALL s8zked3(1,nel,
111 . px(1,8),py(1,8),pz(1,8),px(1,8),py(1,8),pz(1,8),
112 . dd, k88,is)
113 CALL s8zked3(1,nel,
114 . px(1,9),py(1,9),pz(1,9),px(1,9),py(1,9),pz(1,9),
115 . dd, k99,is)
116 CALL s8zked3(1,nel,
117 . px(1,10),py(1,10),pz(1,10),px(1,10),py(1,10),pz(1,10),
118 . dd, k00,is)
119C
120 IF (iksup==0) THEN
121 CALL s8zkeg3(1,nel,
122 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
123 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
124 . gg,k11,is)
125 CALL s8zkeg3(1,nel,
126 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
127 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
128 . gg,k22,is)
129 CALL s8zkeg3(1,nel,
130 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
131 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
132 . gg,k33,is)
133 CALL s8zkeg3(1,nel,
134 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
135 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
136 . gg,k44,is)
137 CALL s8zkeg3(1,nel,
138 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
139 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
140 . gg,k55,is)
141 CALL s8zkeg3(1,nel,
142 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
143 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
144 . gg,k66,is)
145 CALL s8zkeg3(1,nel,
146 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
147 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
148 . gg,k77,is)
149 CALL s8zkeg3(1,nel,
150 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
151 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
152 . gg,k88,is)
153 CALL s8zkeg3(1,nel,
154 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
155 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
156 . gg,k99,is)
157 CALL s8zkeg3(1,nel,
158 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
159 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
160 . gg,k00,is)
161 ENDIF
162C -------non symm------
163 is =0
164 CALL s8zked3(1,nel,
165 . px(1,1),py(1,1),pz(1,1),px(1,2),py(1,2),pz(1,2),
166 . dd, k12,is)
167 CALL s8zked3(1,nel,
168 . px(1,1),py(1,1),pz(1,1),px(1,3),py(1,3),pz(1,3),
169 . dd, k13,is)
170 CALL s8zked3(1,nel,
171 . px(1,1),py(1,1),pz(1,1),px(1,4),py(1,4),pz(1,4),
172 . dd, k14,is)
173 CALL s8zked3(1,nel,
174 . px(1,1),py(1,1),pz(1,1),px(1,5),py(1,5),pz(1,5),
175 . dd, k15,is)
176 CALL s8zked3(1,nel,
177 . px(1,1),py(1,1),pz(1,1),px(1,6),py(1,6),pz(1,6),
178 . dd, k16,is)
179 CALL s8zked3(1,nel,
180 . px(1,1),py(1,1),pz(1,1),px(1,7),py(1,7),pz(1,7),
181 . dd, k17,is)
182 CALL s8zked3(1,nel,
183 . px(1,1),py(1,1),pz(1,1),px(1,8),py(1,8),pz(1,8),
184 . dd, k18,is)
185 CALL s8zked3(1,nel,
186 . px(1,1),py(1,1),pz(1,1),px(1,9),py(1,9),pz(1,9),
187 . dd, k19,is)
188 CALL s8zked3(1,nel,
189 . px(1,1),py(1,1),pz(1,1),px(1,10),py(1,10),pz(1,10),
190 . dd, k10,is)
191 CALL s8zked3(1,nel,
192 . px(1,2),py(1,2),pz(1,2),px(1,3),py(1,3),pz(1,3),
193 . dd, k23,is)
194 CALL s8zked3(1,nel,
195 . px(1,2),py(1,2),pz(1,2),px(1,4),py(1,4),pz(1,4),
196 . dd, k24,is)
197 CALL s8zked3(1,nel,
198 . px(1,2),py(1,2),pz(1,2),px(1,5),py(1,5),pz(1,5),
199 . dd, k25,is)
200 CALL s8zked3(1,nel,
201 . px(1,2),py(1,2),pz(1,2),px(1,6),py(1,6),pz(1,6),
202 . dd, k26,is)
203 CALL s8zked3(1,nel,
204 . px(1,2),py(1,2),pz(1,2),px(1,7),py(1,7),pz(1,7),
205 . dd, k27,is)
206 CALL s8zked3(1,nel,
207 . px(1,2),py(1,2),pz(1,2),px(1,8),py(1,8),pz(1,8),
208 . dd, k28,is)
209 CALL s8zked3(1,nel,
210 . px(1,2),py(1,2),pz(1,2),px(1,9),py(1,9),pz(1,9),
211 . dd, k29,is)
212 CALL s8zked3(1,nel,
213 . px(1,2),py(1,2),pz(1,2),px(1,10),py(1,10),pz(1,10),
214 . dd, k20,is)
215 CALL s8zked3(1,nel,
216 . px(1,3),py(1,3),pz(1,3),px(1,4),py(1,4),pz(1,4),
217 . dd, k34,is)
218 CALL s8zked3(1,nel,
219 . px(1,3),py(1,3),pz(1,3),px(1,5),py(1,5),pz(1,5),
220 . dd, k35,is)
221 CALL s8zked3(1,nel,
222 . px(1,3),py(1,3),pz(1,3),px(1,6),py(1,6),pz(1,6),
223 . dd, k36,is)
224 CALL s8zked3(1,nel,
225 . px(1,3),py(1,3),pz(1,3),px(1,7),py(1,7),pz(1,7),
226 . dd, k37,is)
227 CALL s8zked3(1,nel,
228 . px(1,3),py(1,3),pz(1,3),px(1,8),py(1,8),pz(1,8),
229 . dd, k38,is)
230 CALL s8zked3(1,nel,
231 . px(1,3),py(1,3),pz(1,3),px(1,9),py(1,9),pz(1,9),
232 . dd, k39,is)
233 CALL s8zked3(1,nel,
234 . px(1,3),py(1,3),pz(1,3),px(1,10),py(1,10),pz(1,10),
235 . dd, k30,is)
236 CALL s8zked3(1,nel,
237 . px(1,4),py(1,4),pz(1,4),px(1,5),py(1,5),pz(1,5),
238 . dd, k45,is)
239 CALL s8zked3(1,nel,
240 . px(1,4),py(1,4),pz(1,4),px(1,6),py(1,6),pz(1,6),
241 . dd, k46,is)
242 CALL s8zked3(1,nel,
243 . px(1,4),py(1,4),pz(1,4),px(1,7),py(1,7),pz(1,7),
244 . dd, k47,is)
245 CALL s8zked3(1,nel,
246 . px(1,4),py(1,4),pz(1,4),px(1,8),py(1,8),pz(1,8),
247 . dd, k48,is)
248 CALL s8zked3(1,nel,
249 . px(1,4),py(1,4),pz(1,4),px(1,9),py(1,9),pz(1,9),
250 . dd, k49,is)
251 CALL s8zked3(1,nel,
252 . px(1,4),py(1,4),pz(1,4),px(1,10),py(1,10),pz(1,10),
253 . dd, k40,is)
254 CALL s8zked3(1,nel,
255 . px(1,5),py(1,5),pz(1,5),px(1,6),py(1,6),pz(1,6),
256 . dd, k56,is)
257 CALL s8zked3(1,nel,
258 . px(1,5),py(1,5),pz(1,5),px(1,7),py(1,7),pz(1,7),
259 . dd, k57,is)
260 CALL s8zked3(1,nel,
261 . px(1,5),py(1,5),pz(1,5),px(1,8),py(1,8),pz(1,8),
262 . dd, k58,is)
263 CALL s8zked3(1,nel,
264 . px(1,5),py(1,5),pz(1,5),px(1,9),py(1,9),pz(1,9),
265 . dd, k59,is)
266 CALL s8zked3(1,nel,
267 . px(1,5),py(1,5),pz(1,5),px(1,10),py(1,10),pz(1,10),
268 . dd, k50,is)
269 CALL s8zked3(1,nel,
270 . px(1,6),py(1,6),pz(1,6),px(1,7),py(1,7),pz(1,7),
271 . dd, k67,is)
272 CALL s8zked3(1,nel,
273 . px(1,6),py(1,6),pz(1,6),px(1,8),py(1,8),pz(1,8),
274 . dd, k68,is)
275 CALL s8zked3(1,nel,
276 . px(1,6),py(1,6),pz(1,6),px(1,9),py(1,9),pz(1,9),
277 . dd, k69,is)
278 CALL s8zked3(1,nel,
279 . px(1,6),py(1,6),pz(1,6),px(1,10),py(1,10),pz(1,10),
280 . dd, k60,is)
281 CALL s8zked3(1,nel,
282 . px(1,7),py(1,7),pz(1,7),px(1,8),py(1,8),pz(1,8),
283 . dd, k78,is)
284 CALL s8zked3(1,nel,
285 . px(1,7),py(1,7),pz(1,7),px(1,9),py(1,9),pz(1,9),
286 . dd, k79,is)
287 CALL s8zked3(1,nel,
288 . px(1,7),py(1,7),pz(1,7),px(1,10),py(1,10),pz(1,10),
289 . dd, k70,is)
290 CALL s8zked3(1,nel,
291 . px(1,8),py(1,8),pz(1,8),px(1,9),py(1,9),pz(1,9),
292 . dd, k89,is)
293 CALL s8zked3(1,nel,
294 . px(1,8),py(1,8),pz(1,8),px(1,10),py(1,10),pz(1,10),
295 . dd, k80,is)
296 CALL s8zked3(1,nel,
297 . px(1,9),py(1,9),pz(1,9),px(1,10),py(1,10),pz(1,10),
298 . dd, k90,is)
299C
300 IF (iksup==0) THEN
301 CALL s8zkeg3(1,nel,
302 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
303 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
304 . gg,k12,is)
305 CALL s8zkeg3(1,nel,
306 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
307 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
308 . gg,k13,is)
309 CALL s8zkeg3(1,nel,
310 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
311 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
312 . gg,k14,is)
313 CALL s8zkeg3(1,nel,
314 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
315 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
316 . gg,k15,is)
317 CALL s8zkeg3(1,nel,
318 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
319 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
320 . gg,k16,is)
321 CALL s8zkeg3(1,nel,
322 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
323 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
324 . gg,k17,is)
325 CALL s8zkeg3(1,nel,
326 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
327 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
328 . gg,k18,is)
329 CALL s8zkeg3(1,nel,
330 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
331 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
332 . gg,k19,is)
333 CALL s8zkeg3(1,nel,
334 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
335 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
336 . gg,k10,is)
337 CALL s8zkeg3(1,nel,
338 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
339 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
340 . gg,k23,is)
341 CALL s8zkeg3(1,nel,
342 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
343 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
344 . gg,k24,is)
345 CALL s8zkeg3(1,nel,
346 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
347 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
348 . gg,k25,is)
349 CALL s8zkeg3(1,nel,
350 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
351 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
352 . gg,k26,is)
353 CALL s8zkeg3(1,nel,
354 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
355 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
356 . gg,k27,is)
357 CALL s8zkeg3(1,nel,
358 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
359 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
360 . gg,k28,is)
361 CALL s8zkeg3(1,nel,
362 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
363 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
364 . gg,k29,is)
365 CALL s8zkeg3(1,nel,
366 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
367 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
368 . gg,k20,is)
369 CALL s8zkeg3(1,nel,
370 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
371 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
372 . gg,k34,is)
373 CALL s8zkeg3(1,nel,
374 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
375 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
376 . gg,k35,is)
377 CALL s8zkeg3(1,nel,
378 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
379 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
380 . gg,k36,is)
381 CALL s8zkeg3(1,nel,
382 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
383 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
384 . gg,k37,is)
385 CALL s8zkeg3(1,nel,
386 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
387 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
388 . gg,k38,is)
389 CALL s8zkeg3(1,nel,
390 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
391 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
392 . gg,k39,is)
393 CALL s8zkeg3(1,nel,
394 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
395 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
396 . gg,k30,is)
397 CALL s8zkeg3(1,nel,
398 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
399 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
400 . gg,k45,is)
401 CALL s8zkeg3(1,nel,
402 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
403 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
404 . gg,k46,is)
405 CALL s8zkeg3(1,nel,
406 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
407 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
408 . gg,k47,is)
409 CALL s8zkeg3(1,nel,
410 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
411 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
412 . gg,k48,is)
413 CALL s8zkeg3(1,nel,
414 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
415 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
416 . gg,k49,is)
417 CALL s8zkeg3(1,nel,
418 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
419 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
420 . gg,k40,is)
421 CALL s8zkeg3(1,nel,
422 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
423 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
424 . gg,k56,is)
425 CALL s8zkeg3(1,nel,
426 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
427 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
428 . gg,k57,is)
429 CALL s8zkeg3(1,nel,
430 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
431 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
432 . gg,k58,is)
433 CALL s8zkeg3(1,nel,
434 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
435 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
436 . gg,k59,is)
437 CALL s8zkeg3(1,nel,
438 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
439 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
440 . gg,k50,is)
441 CALL s8zkeg3(1,nel,
442 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
443 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
444 . gg,k67,is)
445 CALL s8zkeg3(1,nel,
446 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
447 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
448 . gg,k68,is)
449 CALL s8zkeg3(1,nel,
450 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
451 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
452 . gg,k69,is)
453 CALL s8zkeg3(1,nel,
454 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
455 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
456 . gg,k60,is)
457 CALL s8zkeg3(1,nel,
458 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
459 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
460 . gg,k78,is)
461 CALL s8zkeg3(1,nel,
462 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
463 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
464 . gg,k79,is)
465 CALL s8zkeg3(1,nel,
466 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
467 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
468 . gg,k70,is)
469 CALL s8zkeg3(1,nel,
470 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
471 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
472 . gg,k89,is)
473 CALL s8zkeg3(1,nel,
474 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
475 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
476 . gg,k80,is)
477 CALL s8zkeg3(1,nel,
478 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
479 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
480 . gg,k90,is)
481 ELSE
482 CALL s8zkeg33(1,nel,
483 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
484 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
485 . g33,k12,is)
486 CALL s8zkeg33(1,nel,
487 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
488 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
489 . g33,k13,is)
490 CALL s8zkeg33(1,nel,
491 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
492 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
493 . g33,k14,is)
494 CALL s8zkeg33(1,nel,
495 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
496 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
497 . g33,k15,is)
498 CALL s8zkeg33(1,nel,
499 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
500 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
501 . g33,k16,is)
502 CALL s8zkeg33(1,nel,
503 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
504 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
505 . g33,k17,is)
506 CALL s8zkeg33(1,nel,
507 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
508 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
509 . g33,k18,is)
510 CALL s8zkeg33(1,nel,
511 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
512 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
513 . g33,k19,is)
514 CALL s8zkeg33(1,nel,
515 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
516 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
517 . g33,k10,is)
518 CALL s8zkeg33(1,nel,
519 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
520 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
521 . g33,k23,is)
522 CALL s8zkeg33(1,nel,
523 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
524 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
525 . g33,k24,is)
526 CALL s8zkeg33(1,nel,
527 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
528 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
529 . g33,k25,is)
530 CALL s8zkeg33(1,nel,
531 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
532 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
533 . g33,k26,is)
534 CALL s8zkeg33(1,nel,
535 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
536 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
537 . g33,k27,is)
538 CALL s8zkeg33(1,nel,
539 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
540 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
541 . g33,k28,is)
542 CALL s8zkeg33(1,nel,
543 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
544 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
545 . g33,k29,is)
546 CALL s8zkeg33(1,nel,
547 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
548 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
549 . g33,k20,is)
550 CALL s8zkeg33(1,nel,
551 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
552 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
553 . g33,k34,is)
554 CALL s8zkeg33(1,nel,
555 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
556 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
557 . g33,k35,is)
558 CALL s8zkeg33(1,nel,
559 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
560 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
561 . g33,k36,is)
562 CALL s8zkeg33(1,nel,
563 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
564 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
565 . g33,k37,is)
566 CALL s8zkeg33(1,nel,
567 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
568 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
569 . g33,k38,is)
570 CALL s8zkeg33(1,nel,
571 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
572 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
573 . g33,k39,is)
574 CALL s8zkeg33(1,nel,
575 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
576 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
577 . g33,k30,is)
578 CALL s8zkeg33(1,nel,
579 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
580 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
581 . g33,k45,is)
582 CALL s8zkeg33(1,nel,
583 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
584 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
585 . g33,k46,is)
586 CALL s8zkeg33(1,nel,
587 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
588 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
589 . g33,k47,is)
590 CALL s8zkeg33(1,nel,
591 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
592 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
593 . g33,k48,is)
594 CALL s8zkeg33(1,nel,
595 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
596 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
597 . g33,k49,is)
598 CALL s8zkeg33(1,nel,
599 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
600 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
601 . g33,k40,is)
602 CALL s8zkeg33(1,nel,
603 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
604 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
605 . g33,k56,is)
606 CALL s8zkeg33(1,nel,
607 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
608 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
609 . g33,k57,is)
610 CALL s8zkeg33(1,nel,
611 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
612 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
613 . g33,k58,is)
614 CALL s8zkeg33(1,nel,
615 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
616 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
617 . g33,k59,is)
618 CALL s8zkeg33(1,nel,
619 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
620 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
621 . g33,k50,is)
622 CALL s8zkeg33(1,nel,
623 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
624 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
625 . g33,k67,is)
626 CALL s8zkeg33(1,nel,
627 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
628 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
629 . g33,k68,is)
630 CALL s8zkeg33(1,nel,
631 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
632 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
633 . g33,k69,is)
634 CALL s8zkeg33(1,nel,
635 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
636 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
637 . g33,k60,is)
638 CALL s8zkeg33(1,nel,
639 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
640 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
641 . g33,k78,is)
642 CALL s8zkeg33(1,nel,
643 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
644 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
645 . g33,k79,is)
646 CALL s8zkeg33(1,nel,
647 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
648 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
649 . g33,k70,is)
650 CALL s8zkeg33(1,nel,
651 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
652 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
653 . g33,k89,is)
654 CALL s8zkeg33(1,nel,
655 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
656 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
657 . g33,k80,is)
658 CALL s8zkeg33(1,nel,
659 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
660 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
661 . g33,k90,is)
662C
663 CALL s8zkebg1(1,nel,
664 . px(1,1),py(1,1),pz(1,1),px(1,2),py(1,2),pz(1,2),
665 . dg, k12,is)
666 CALL s8zkebg1(1,nel,
667 . px(1,1),py(1,1),pz(1,1),px(1,3),py(1,3),pz(1,3),
668 . dg, k13,is)
669 CALL s8zkebg1(1,nel,
670 . px(1,1),py(1,1),pz(1,1),px(1,4),py(1,4),pz(1,4),
671 . dg, k14,is)
672 CALL s8zkebg1(1,nel,
673 . px(1,1),py(1,1),pz(1,1),px(1,5),py(1,5),pz(1,5),
674 . dg, k15,is)
675 CALL s8zkebg1(1,nel,
676 . px(1,1),py(1,1),pz(1,1),px(1,6),py(1,6),pz(1,6),
677 . dg, k16,is)
678 CALL s8zkebg1(1,nel,
679 . px(1,1),py(1,1),pz(1,1),px(1,7),py(1,7),pz(1,7),
680 . dg, k17,is)
681 CALL s8zkebg1(1,nel,
682 . px(1,1),py(1,1),pz(1,1),px(1,8),py(1,8),pz(1,8),
683 . dg, k18,is)
684 CALL s8zkebg1(1,nel,
685 . px(1,1),py(1,1),pz(1,1),px(1,9),py(1,9),pz(1,9),
686 . dg, k19,is)
687 CALL s8zkebg1(1,nel,
688 . px(1,1),py(1,1),pz(1,1),px(1,10),py(1,10),pz(1,10),
689 . dg, k10,is)
690 CALL s8zkebg1(1,nel,
691 . px(1,2),py(1,2),pz(1,2),px(1,3),py(1,3),pz(1,3),
692 . dg, k23,is)
693 CALL s8zkebg1(1,nel,
694 . px(1,2),py(1,2),pz(1,2),px(1,4),py(1,4),pz(1,4),
695 . dg, k24,is)
696 CALL s8zkebg1(1,nel,
697 . px(1,2),py(1,2),pz(1,2),px(1,5),py(1,5),pz(1,5),
698 . dg, k25,is)
699 CALL s8zkebg1(1,nel,
700 . px(1,2),py(1,2),pz(1,2),px(1,6),py(1,6),pz(1,6),
701 . dg, k26,is)
702 CALL s8zkebg1(1,nel,
703 . px(1,2),py(1,2),pz(1,2),px(1,7),py(1,7),pz(1,7),
704 . dg, k27,is)
705 CALL s8zkebg1(1,nel,
706 . px(1,2),py(1,2),pz(1,2),px(1,8),py(1,8),pz(1,8),
707 . dg, k28,is)
708 CALL s8zkebg1(1,nel,
709 . px(1,2),py(1,2),pz(1,2),px(1,9),py(1,9),pz(1,9),
710 . dg, k29,is)
711 CALL s8zkebg1(1,nel,
712 . px(1,2),py(1,2),pz(1,2),px(1,10),py(1,10),pz(1,10),
713 . dg, k20,is)
714 CALL s8zkebg1(1,nel,
715 . px(1,3),py(1,3),pz(1,3),px(1,4),py(1,4),pz(1,4),
716 . dg, k34,is)
717 CALL s8zkebg1(1,nel,
718 . px(1,3),py(1,3),pz(1,3),px(1,5),py(1,5),pz(1,5),
719 . dg, k35,is)
720 CALL s8zkebg1(1,nel,
721 . px(1,3),py(1,3),pz(1,3),px(1,6),py(1,6),pz(1,6),
722 . dg, k36,is)
723 CALL s8zkebg1(1,nel,
724 . px(1,3),py(1,3),pz(1,3),px(1,7),py(1,7),pz(1,7),
725 . dg, k37,is)
726 CALL s8zkebg1(1,nel,
727 . px(1,3),py(1,3),pz(1,3),px(1,8),py(1,8),pz(1,8),
728 . dg, k38,is)
729 CALL s8zkebg1(1,nel,
730 . px(1,3),py(1,3),pz(1,3),px(1,9),py(1,9),pz(1,9),
731 . dg, k39,is)
732 CALL s8zkebg1(1,nel,
733 . px(1,3),py(1,3),pz(1,3),px(1,10),py(1,10),pz(1,10),
734 . dg, k30,is)
735 CALL s8zkebg1(1,nel,
736 . px(1,4),py(1,4),pz(1,4),px(1,5),py(1,5),pz(1,5),
737 . dg, k45,is)
738 CALL s8zkebg1(1,nel,
739 . px(1,4),py(1,4),pz(1,4),px(1,6),py(1,6),pz(1,6),
740 . dg, k46,is)
741 CALL s8zkebg1(1,nel,
742 . px(1,4),py(1,4),pz(1,4),px(1,7),py(1,7),pz(1,7),
743 . dg, k47,is)
744 CALL s8zkebg1(1,nel,
745 . px(1,4),py(1,4),pz(1,4),px(1,8),py(1,8),pz(1,8),
746 . dg, k48,is)
747 CALL s8zkebg1(1,nel,
748 . px(1,4),py(1,4),pz(1,4),px(1,9),py(1,9),pz(1,9),
749 . dg, k49,is)
750 CALL s8zkebg1(1,nel,
751 . px(1,4),py(1,4),pz(1,4),px(1,10),py(1,10),pz(1,10),
752 . dg, k40,is)
753 CALL s8zkebg1(1,nel,
754 . px(1,5),py(1,5),pz(1,5),px(1,6),py(1,6),pz(1,6),
755 . dg, k56,is)
756 CALL s8zkebg1(1,nel,
757 . px(1,5),py(1,5),pz(1,5),px(1,7),py(1,7),pz(1,7),
758 . dg, k57,is)
759 CALL s8zkebg1(1,nel,
760 . px(1,5),py(1,5),pz(1,5),px(1,8),py(1,8),pz(1,8),
761 . dg, k58,is)
762 CALL s8zkebg1(1,nel,
763 . px(1,5),py(1,5),pz(1,5),px(1,9),py(1,9),pz(1,9),
764 . dg, k59,is)
765 CALL s8zkebg1(1,nel,
766 . px(1,5),py(1,5),pz(1,5),px(1,10),py(1,10),pz(1,10),
767 . dg, k50,is)
768 CALL s8zkebg1(1,nel,
769 . px(1,6),py(1,6),pz(1,6),px(1,7),py(1,7),pz(1,7),
770 . dg, k67,is)
771 CALL s8zkebg1(1,nel,
772 . px(1,6),py(1,6),pz(1,6),px(1,8),py(1,8),pz(1,8),
773 . dg, k68,is)
774 CALL s8zkebg1(1,nel,
775 . px(1,6),py(1,6),pz(1,6),px(1,9),py(1,9),pz(1,9),
776 . dg, k69,is)
777 CALL s8zkebg1(1,nel,
778 . px(1,6),py(1,6),pz(1,6),px(1,10),py(1,10),pz(1,10),
779 . dg, k60,is)
780 CALL s8zkebg1(1,nel,
781 . px(1,7),py(1,7),pz(1,7),px(1,8),py(1,8),pz(1,8),
782 . dg, k78,is)
783 CALL s8zkebg1(1,nel,
784 . px(1,7),py(1,7),pz(1,7),px(1,9),py(1,9),pz(1,9),
785 . dg, k79,is)
786 CALL s8zkebg1(1,nel,
787 . px(1,7),py(1,7),pz(1,7),px(1,10),py(1,10),pz(1,10),
788 . dg, k70,is)
789 CALL s8zkebg1(1,nel,
790 . px(1,8),py(1,8),pz(1,8),px(1,9),py(1,9),pz(1,9),
791 . dg, k89,is)
792 CALL s8zkebg1(1,nel,
793 . px(1,8),py(1,8),pz(1,8),px(1,10),py(1,10),pz(1,10),
794 . dg, k80,is)
795 CALL s8zkebg1(1,nel,
796 . px(1,9),py(1,9),pz(1,9),px(1,10),py(1,10),pz(1,10),
797 . dg, k90,is)
798C
799 is =1
800 CALL s8zkebg1(1,nel,
801 . px(1,1),py(1,1),pz(1,1),px(1,1),py(1,1),pz(1,1),
802 . dg, k11,is)
803 CALL s8zkebg1(1,nel,
804 . px(1,2),py(1,2),pz(1,2),px(1,2),py(1,2),pz(1,2),
805 . dg, k22,is)
806 CALL s8zkebg1(1,nel,
807 . px(1,3),py(1,3),pz(1,3),px(1,3),py(1,3),pz(1,3),
808 . dg, k33,is)
809 CALL s8zkebg1(1,nel,
810 . px(1,4),py(1,4),pz(1,4),px(1,4),py(1,4),pz(1,4),
811 . dg, k44,is)
812 CALL s8zkebg1(1,nel,
813 . px(1,5),py(1,5),pz(1,5),px(1,5),py(1,5),pz(1,5),
814 . dg, k55,is)
815 CALL s8zkebg1(1,nel,
816 . px(1,6),py(1,6),pz(1,6),px(1,6),py(1,6),pz(1,6),
817 . dg, k66,is)
818 CALL s8zkebg1(1,nel,
819 . px(1,7),py(1,7),pz(1,7),px(1,7),py(1,7),pz(1,7),
820 . dg, k77,is)
821 CALL s8zkebg1(1,nel,
822 . px(1,8),py(1,8),pz(1,8),px(1,8),py(1,8),pz(1,8),
823 . dg, k88,is)
824 CALL s8zkebg1(1,nel,
825 . px(1,9),py(1,9),pz(1,9),px(1,9),py(1,9),pz(1,9),
826 . dg, k99,is)
827 CALL s8zkebg1(1,nel,
828 . px(1,10),py(1,10),pz(1,10),px(1,10),py(1,10),pz(1,10),
829 . dg, k00,is)
830C
831 CALL s8zkeg33(1,nel,
832 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
833 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
834 . g33,k11,is)
835 CALL s8zkeg33(1,nel,
836 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
837 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
838 . g33,k22,is)
839 CALL s8zkeg33(1,nel,
840 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
841 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
842 . g33,k33,is)
843 CALL s8zkeg33(1,nel,
844 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
845 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
846 . g33,k44,is)
847 CALL s8zkeg33(1,nel,
848 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
849 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
850 . g33,k55,is)
851 CALL s8zkeg33(1,nel,
852 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
853 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
854 . g33,k66,is)
855 CALL s8zkeg33(1,nel,
856 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
857 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
858 . g33,k77,is)
859 CALL s8zkeg33(1,nel,
860 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
861 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
862 . g33,k88,is)
863 CALL s8zkeg33(1,nel,
864 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
865 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
866 . g33,k99,is)
867 CALL s8zkeg33(1,nel,
868 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
869 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
870 . g33,k00,is)
871 ENDIF
872C
873 RETURN
874 END
s10cumg3
subroutine s10cumg3(px, py, pz, dd, gg, dg, g33, iksup, k11, k12, k13, k14, k15, k16, k17, k18, k19, k10, k22, k23, k24, k25, k26, k27, k28, k29, k20, k33, k34, k35, k36, k37, k38, k39, k30, k44, k45, k46, k47, k48, k49, k40, k55, k56, k57, k58, k59, k50, k66, k67, k68, k69, k60, k77, k78, k79, k70, k88, k89, k80, k99, k90, k00, nel)
Definition s10cumg3.F:50
s8zkebg1
subroutine s8zkebg1(lft, llt, pxi, pyi, pzi, pxj, pyj, pzj, dg, kij, is)
Definition s8zkebg1.F:32
s8zked3
subroutine s8zked3(lft, llt, pxi, pyi, pzi, pxj, pyj, pzj, dd, kij, is)
Definition s8zked3.F:34
s8zkeg33
subroutine s8zkeg33(lft, llt, pxyi, pyxi, pxzi, pzxi, pyzi, pzyi, pxyj, pyxj, pxzj, pzxj, pyzj, pzyj, g33, kij, is)
Definition s8zkeg33.F:34
s8zkeg3
subroutine s8zkeg3(lft, llt, pxyi, pyxi, pxzi, pzxi, pyzi, pzyi, pxyj, pyxj, pxzj, pzxj, pyzj, pzyj, g, kij, is)
Definition s8zkeg3.F:34