OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s20cumg3.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!|| s20cumg3 ../engine/source/elements/solid/solide20/s20cumg3.f
25!||--- called by ------------------------------------------------------
26!|| s20ke3 ../engine/source/elements/solid/solide20/s20ke3.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 s20cumg3(
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, K110, K111, K112,
39 6 K113, K114, K115, K116,
40 7 K117, K118, K119, K120,
41 8 K22, K23, K24, K25,
42 9 K26, K27, K28, K29,
43 A K210, K211, K212, K213,
44 B K214, K215, K216, K217,
45 C K218, K219, K220, K33,
46 D K34, K35, K36, K37,
47 E K38, K39, K310, K311,
48 F K312, K313, K314, K315,
49 G K316, K317, K318, K319,
50 H K320, K44, K45, K46,
51 I K47, K48, K49, K410,
52 J K411, K412, K413, K414,
53 K K415, K416, K417, K418,
54 L K419, K420, K55, K56,
55 M K57, K58, K59, K510,
56 N K511, K512, K513, K514,
57 O K515, K516, K517, K518,
58 P K519, K520, K66, K67,
59 Q K68, K69, K610, K611,
60 R K612, K613, K614, K615,
61 S K616, K617, K618, K619,
62 T K620, K77, K78, K79,
63 U K710, K711, K712, K713,
64 V K714, K715, K716, K717,
65 W K718, K719, K720, K88,
66 X K89, K810, K811, K812,
67 Y K813, K814, K815, K816,
68 Z K817, K818, K819, K820,
69 1 K99, K910, K911, K912,
70 2 K913, K914, K915, K916,
71 3 K917, K918, K919, K920,
72 4 K1010, K1011, K1012, K1013,
73 5 K1014, K1015, K1016, K1017,
74 6 K1018, K1019, K1020, K1111,
75 7 K1112, K1113, K1114, K1115,
76 8 K1116, K1117, K1118, K1119,
77 9 K1120, K1212, K1213, K1214,
78 A K1215, K1216, K1217, K1218,
79 B K1219, K1220, K1313, K1314,
80 C K1315, K1316, K1317, K1318,
81 D K1319, K1320, K1414, K1415,
82 E K1416, K1417, K1418, K1419,
83 F K1420, K1515, K1516, K1517,
84 G K1518, K1519, K1520, K1616,
85 H K1617, K1618, K1619, K1620,
86 I K1717, K1718, K1719, K1720,
87 J K1818, K1819, K1820, K1919,
88 K K1920, K2020, NEL)
89C-----------------------------------------------
90C I m p l i c i t T y p e s
91C-----------------------------------------------
92#include "implicit_f.inc"
93C-----------------------------------------------
94C G l o b a l P a r a m e t e r s
95C-----------------------------------------------
96#include "mvsiz_p.inc"
97C-----------------------------------------------
98C C o m m o n B l o c k s
99C-----------------------------------------------
100C-----------------------------------------------
101C D u m m y A r g u m e n t s
102C-----------------------------------------------
103 INTEGER, INTENT(IN) :: NEL
104 INTEGER IKSUP
105C REAL
106 my_real
107 . PX(MVSIZ,20),PY(MVSIZ,20),PZ(MVSIZ,20),
108 . DD(3,3,*),GG(*),DG(9,*) ,G33(9,*)
109 my_real
110 . K11(9,*) ,K12(9,*) ,K13(9,*) ,K14(9,*) ,K15(9,*) ,
111 . K16(9,*) ,K17(9,*) ,K18(9,*) ,K19(9,*) ,K110(9,*) ,
112 . K111(9,*) ,K112(9,*) ,K113(9,*) ,K114(9,*) ,K115(9,*) ,
113 . K116(9,*) ,K117(9,*) ,K118(9,*) ,K119(9,*) ,K120(9,*) ,
114 . K22(9,*) ,K23(9,*) ,K24(9,*) ,K25(9,*) ,K26(9,*) ,
115 . K27(9,*) ,K28(9,*) ,K29(9,*) ,K210(9,*) ,K211(9,*) ,
116 . K212(9,*) ,K213(9,*) ,K214(9,*) ,K215(9,*) ,K216(9,*) ,
117 . K217(9,*) ,K218(9,*) ,K219(9,*) ,K220(9,*) ,K33(9,*) ,
118 . K34(9,*) ,K35(9,*) ,K36(9,*) ,K37(9,*) ,K38(9,*) ,
119 . K39(9,*) ,K310(9,*) ,K311(9,*) ,K312(9,*) ,K313(9,*) ,
120 . K314(9,*) ,K315(9,*) ,K316(9,*) ,K317(9,*) ,K318(9,*) ,
121 . K319(9,*) ,K320(9,*) ,K44(9,*) ,K45(9,*) ,K46(9,*) ,
122 . K47(9,*) ,K48(9,*) ,K49(9,*) ,K410(9,*) ,K411(9,*) ,
123 . K412(9,*) ,K413(9,*) ,K414(9,*) ,K415(9,*) ,K416(9,*) ,
124 . K417(9,*) ,K418(9,*) ,K419(9,*) ,K420(9,*) ,K55(9,*) ,
125 . K56(9,*) ,K57(9,*) ,K58(9,*) ,K59(9,*) ,K510(9,*) ,
126 . K511(9,*) ,K512(9,*) ,K513(9,*) ,K514(9,*) ,K515(9,*) ,
127 . K516(9,*) ,K517(9,*) ,K518(9,*) ,K519(9,*) ,K520(9,*) ,
128 . K66(9,*) ,K67(9,*) ,K68(9,*) ,K69(9,*) ,K610(9,*) ,
129 . K611(9,*) ,K612(9,*) ,K613(9,*) ,K614(9,*) ,K615(9,*) ,
130 . K616(9,*) ,K617(9,*) ,K618(9,*) ,K619(9,*) ,K620(9,*) ,
131 . K77(9,*) ,K78(9,*) ,K79(9,*) ,K710(9,*) ,K711(9,*) ,
132 . K712(9,*) ,K713(9,*) ,K714(9,*) ,K715(9,*) ,K716(9,*) ,
133 . K717(9,*) ,K718(9,*) ,K719(9,*) ,K720(9,*) ,K88(9,*) ,
134 . K89(9,*) ,K810(9,*) ,K811(9,*) ,K812(9,*) ,K813(9,*) ,
135 . K814(9,*) ,K815(9,*) ,K816(9,*) ,K817(9,*) ,K818(9,*) ,
136 . K819(9,*) ,K820(9,*) ,K99(9,*) ,K910(9,*) ,K911(9,*) ,
137 . K912(9,*) ,K913(9,*) ,K914(9,*) ,K915(9,*) ,K916(9,*) ,
138 . K917(9,*) ,K918(9,*) ,K919(9,*) ,K920(9,*) ,K1010(9,*) ,
139 . K1011(9,*) ,K1012(9,*) ,K1013(9,*) ,K1014(9,*) ,K1015(9,*) ,
140 . K1016(9,*) ,K1017(9,*) ,K1018(9,*) ,K1019(9,*) ,K1020(9,*) ,
141 . K1111(9,*) ,K1112(9,*) ,K1113(9,*) ,K1114(9,*) ,K1115(9,*) ,
142 . K1116(9,*) ,K1117(9,*) ,K1118(9,*) ,K1119(9,*) ,K1120(9,*) ,
143 . K1212(9,*) ,K1213(9,*) ,K1214(9,*) ,K1215(9,*) ,K1216(9,*) ,
144 . K1217(9,*) ,K1218(9,*) ,K1219(9,*) ,K1220(9,*) ,K1313(9,*) ,
145 . K1314(9,*) ,K1315(9,*) ,K1316(9,*) ,K1317(9,*) ,K1318(9,*) ,
146 . K1319(9,*) ,K1320(9,*) ,K1414(9,*) ,K1415(9,*) ,K1416(9,*) ,
147 . K1417(9,*) ,K1418(9,*) ,K1419(9,*) ,K1420(9,*) ,K1515(9,*) ,
148 . K1516(9,*) ,K1517(9,*) ,K1518(9,*) ,K1519(9,*) ,K1520(9,*) ,
149 . K1616(9,*) ,K1617(9,*) ,K1618(9,*) ,K1619(9,*) ,K1620(9,*) ,
150 . K1717(9,*) ,K1718(9,*) ,K1719(9,*) ,K1720(9,*) ,K1818(9,*) ,
151 . K1819(9,*) ,K1820(9,*) ,K1919(9,*) ,K1920(9,*) ,K2020(9,*)
152C-----------------------------------------------
153C L o c a l V a r i a b l e s
154C-----------------------------------------------
155 INTEGER I, J,IS
156C-----------------------------------------------
157C -------symm(diag)------
158 is =1
159C 12
160 CALL s8zked3(1,nel,
161 . px(1,1),py(1,1),pz(1,1),px(1,1),py(1,1),pz(1,1),
162 . dd, k11,is)
163 CALL s8zked3(1,nel,
164 . px(1,2),py(1,2),pz(1,2),px(1,2),py(1,2),pz(1,2),
165 . dd, k22,is)
166 CALL s8zked3(1,nel,
167 . px(1,3),py(1,3),pz(1,3),px(1,3),py(1,3),pz(1,3),
168 . dd, k33,is)
169 CALL s8zked3(1,nel,
170 . px(1,4),py(1,4),pz(1,4),px(1,4),py(1,4),pz(1,4),
171 . dd, k44,is)
172 CALL s8zked3(1,nel,
173 . px(1,5),py(1,5),pz(1,5),px(1,5),py(1,5),pz(1,5),
174 . dd, k55,is)
175 CALL s8zked3(1,nel,
176 . px(1,6),py(1,6),pz(1,6),px(1,6),py(1,6),pz(1,6),
177 . dd, k66,is)
178 CALL s8zked3(1,nel,
179 . px(1,7),py(1,7),pz(1,7),px(1,7),py(1,7),pz(1,7),
180 . dd, k77,is)
181 CALL s8zked3(1,nel,
182 . px(1,8),py(1,8),pz(1,8),px(1,8),py(1,8),pz(1,8),
183 . dd, k88,is)
184 CALL s8zked3(1,nel,
185 . px(1,9),py(1,9),pz(1,9),px(1,9),py(1,9),pz(1,9),
186 . dd, k99,is)
187 CALL s8zked3(1,nel,
188 . px(1,10),py(1,10),pz(1,10),px(1,10),py(1,10),pz(1,10),
189 . dd, k1010,is)
190 CALL s8zked3(1,nel,
191 . px(1,11),py(1,11),pz(1,11),px(1,11),py(1,11),pz(1,11),
192 . dd, k1111,is)
193 CALL s8zked3(1,nel,
194 . px(1,12),py(1,12),pz(1,12),px(1,12),py(1,12),pz(1,12),
195 . dd, k1212,is)
196 CALL s8zked3(1,nel,
197 . px(1,13),py(1,13),pz(1,13),px(1,13),py(1,13),pz(1,13),
198 . dd, k1313,is)
199 CALL s8zked3(1,nel,
200 . px(1,14),py(1,14),pz(1,14),px(1,14),py(1,14),pz(1,14),
201 . dd, k1414,is)
202 CALL s8zked3(1,nel,
203 . px(1,15),py(1,15),pz(1,15),px(1,15),py(1,15),pz(1,15),
204 . dd, k1515,is)
205 CALL s8zked3(1,nel,
206 . px(1,16),py(1,16),pz(1,16),px(1,16),py(1,16),pz(1,16),
207 . dd, k1616,is)
208 CALL s8zked3(1,nel,
209 . px(1,17),py(1,17),pz(1,17),px(1,17),py(1,17),pz(1,17),
210 . dd, k1717,is)
211 CALL s8zked3(1,nel,
212 . px(1,18),py(1,18),pz(1,18),px(1,18),py(1,18),pz(1,18),
213 . dd, k1818,is)
214 CALL s8zked3(1,nel,
215 . px(1,19),py(1,19),pz(1,19),px(1,19),py(1,19),pz(1,19),
216 . dd, k1919,is)
217 CALL s8zked3(1,nel,
218 . px(1,20),py(1,20),pz(1,20),px(1,20),py(1,20),pz(1,20),
219 . dd, k2020,is)
220C
221 IF (iksup==0) THEN
222 CALL s8zkeg3(1,nel,
223 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
224 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
225 . gg,k11,is)
226 CALL s8zkeg3(1,nel,
227 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
228 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
229 . gg,k22,is)
230 CALL s8zkeg3(1,nel,
231 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
232 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
233 . gg,k33,is)
234 CALL s8zkeg3(1,nel,
235 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
236 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
237 . gg,k44,is)
238 CALL s8zkeg3(1,nel,
239 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
240 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
241 . gg,k55,is)
242 CALL s8zkeg3(1,nel,
243 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
244 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
245 . gg,k66,is)
246 CALL s8zkeg3(1,nel,
247 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
248 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
249 . gg,k77,is)
250 CALL s8zkeg3(1,nel,
251 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
252 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
253 . gg,k88,is)
254 CALL s8zkeg3(1,nel,
255 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
256 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
257 . gg,k99,is)
258 CALL s8zkeg3(1,nel,
259 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
260 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
261 . gg,k1010,is)
262 CALL s8zkeg3(1,nel,
263 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
264 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
265 . gg,k1111,is)
266 CALL s8zkeg3(1,nel,
267 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
268 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
269 . gg,k1212,is)
270 CALL s8zkeg3(1,nel,
271 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
272 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
273 . gg,k1313,is)
274 CALL s8zkeg3(1,nel,
275 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
276 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
277 . gg,k1414,is)
278 CALL s8zkeg3(1,nel,
279 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
280 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
281 . gg,k1515,is)
282 CALL s8zkeg3(1,nel,
283 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
284 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
285 . gg,k1616,is)
286 CALL s8zkeg3(1,nel,
287 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
288 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
289 . gg,k1717,is)
290 CALL s8zkeg3(1,nel,
291 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
292 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
293 . gg,k1818,is)
294 CALL s8zkeg3(1,nel,
295 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
296 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
297 . gg,k1919,is)
298 CALL s8zkeg3(1,nel,
299 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
300 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
301 . gg,k2020,is)
302
303 ENDIF
304C -------non symm------
305 is =0
306 CALL s8zked3(1,nel,
307 . px(1,1),py(1,1),pz(1,1),px(1,2),py(1,2),pz(1,2),
308 . dd, k12,is)
309 CALL s8zked3(1,nel,
310 . px(1,1),py(1,1),pz(1,1),px(1,3),py(1,3),pz(1,3),
311 . dd, k13,is)
312 CALL s8zked3(1,nel,
313 . px(1,1),py(1,1),pz(1,1),px(1,4),py(1,4),pz(1,4),
314 . dd, k14,is)
315 CALL s8zked3(1,nel,
316 . px(1,1),py(1,1),pz(1,1),px(1,5),py(1,5),pz(1,5),
317 . dd, k15,is)
318 CALL s8zked3(1,nel,
319 . px(1,1),py(1,1),pz(1,1),px(1,6),py(1,6),pz(1,6),
320 . dd, k16,is)
321 CALL s8zked3(1,nel,
322 . px(1,1),py(1,1),pz(1,1),px(1,7),py(1,7),pz(1,7),
323 . dd, k17,is)
324 CALL s8zked3(1,nel,
325 . px(1,1),py(1,1),pz(1,1),px(1,8),py(1,8),pz(1,8),
326 . dd, k18,is)
327 CALL s8zked3(1,nel,
328 . px(1,1),py(1,1),pz(1,1),px(1,9),py(1,9),pz(1,9),
329 . dd, k19,is)
330 CALL s8zked3(1,nel,
331 . px(1,1),py(1,1),pz(1,1),px(1,10),py(1,10),pz(1,10),
332 . dd, k110,is)
333 CALL s8zked3(1,nel,
334 . px(1,1),py(1,1),pz(1,1),px(1,11),py(1,11),pz(1,11),
335 . dd, k111,is)
336 CALL s8zked3(1,nel,
337 . px(1,1),py(1,1),pz(1,1),px(1,12),py(1,12),pz(1,12),
338 . dd, k112,is)
339 CALL s8zked3(1,nel,
340 . px(1,1),py(1,1),pz(1,1),px(1,13),py(1,13),pz(1,13),
341 . dd, k113,is)
342 CALL s8zked3(1,nel,
343 . px(1,1),py(1,1),pz(1,1),px(1,14),py(1,14),pz(1,14),
344 . dd, k114,is)
345 CALL s8zked3(1,nel,
346 . px(1,1),py(1,1),pz(1,1),px(1,15),py(1,15),pz(1,15),
347 . dd, k115,is)
348 CALL s8zked3(1,nel,
349 . px(1,1),py(1,1),pz(1,1),px(1,16),py(1,16),pz(1,16),
350 . dd, k116,is)
351 CALL s8zked3(1,nel,
352 . px(1,1),py(1,1),pz(1,1),px(1,17),py(1,17),pz(1,17),
353 . dd, k117,is)
354 CALL s8zked3(1,nel,
355 . px(1,1),py(1,1),pz(1,1),px(1,18),py(1,18),pz(1,18),
356 . dd, k118,is)
357 CALL s8zked3(1,nel,
358 . px(1,1),py(1,1),pz(1,1),px(1,19),py(1,19),pz(1,19),
359 . dd, k119,is)
360 CALL s8zked3(1,nel,
361 . px(1,1),py(1,1),pz(1,1),px(1,20),py(1,20),pz(1,20),
362 . dd, k120,is)
363
364
365
366 CALL s8zked3(1,nel,
367 . px(1,2),py(1,2),pz(1,2),px(1,3),py(1,3),pz(1,3),
368 . dd, k23,is)
369 CALL s8zked3(1,nel,
370 . px(1,2),py(1,2),pz(1,2),px(1,4),py(1,4),pz(1,4),
371 . dd, k24,is)
372 CALL s8zked3(1,nel,
373 . px(1,2),py(1,2),pz(1,2),px(1,5),py(1,5),pz(1,5),
374 . dd, k25,is)
375 CALL s8zked3(1,nel,
376 . px(1,2),py(1,2),pz(1,2),px(1,6),py(1,6),pz(1,6),
377 . dd, k26,is)
378 CALL s8zked3(1,nel,
379 . px(1,2),py(1,2),pz(1,2),px(1,7),py(1,7),pz(1,7),
380 . dd, k27,is)
381 CALL s8zked3(1,nel,
382 . px(1,2),py(1,2),pz(1,2),px(1,8),py(1,8),pz(1,8),
383 . dd, k28,is)
384 CALL s8zked3(1,nel,
385 . px(1,2),py(1,2),pz(1,2),px(1,9),py(1,9),pz(1,9),
386 . dd, k29,is)
387 CALL s8zked3(1,nel,
388 . px(1,2),py(1,2),pz(1,2),px(1,10),py(1,10),pz(1,10),
389 . dd, k210,is)
390 CALL s8zked3(1,nel,
391 . px(1,2),py(1,2),pz(1,2),px(1,11),py(1,11),pz(1,11),
392 . dd, k211,is)
393 CALL s8zked3(1,nel,
394 . px(1,2),py(1,2),pz(1,2),px(1,12),py(1,12),pz(1,12),
395 . dd, k212,is)
396 CALL s8zked3(1,nel,
397 . px(1,2),py(1,2),pz(1,2),px(1,13),py(1,13),pz(1,13),
398 . dd, k213,is)
399 CALL s8zked3(1,nel,
400 . px(1,2),py(1,2),pz(1,2),px(1,14),py(1,14),pz(1,14),
401 . dd, k214,is)
402 CALL s8zked3(1,nel,
403 . px(1,2),py(1,2),pz(1,2),px(1,15),py(1,15),pz(1,15),
404 . dd, k215,is)
405 CALL s8zked3(1,nel,
406 . px(1,2),py(1,2),pz(1,2),px(1,16),py(1,16),pz(1,16),
407 . dd, k216,is)
408 CALL s8zked3(1,nel,
409 . px(1,2),py(1,2),pz(1,2),px(1,17),py(1,17),pz(1,17),
410 . dd, k217,is)
411 CALL s8zked3(1,nel,
412 . px(1,2),py(1,2),pz(1,2),px(1,18),py(1,18),pz(1,18),
413 . dd, k218,is)
414 CALL s8zked3(1,nel,
415 . px(1,2),py(1,2),pz(1,2),px(1,19),py(1,19),pz(1,19),
416 . dd, k219,is)
417 CALL s8zked3(1,nel,
418 . px(1,2),py(1,2),pz(1,2),px(1,20),py(1,20),pz(1,20),
419 . dd, k220,is)
420
421
422
423 CALL s8zked3(1,nel,
424 . px(1,3),py(1,3),pz(1,3),px(1,4),py(1,4),pz(1,4),
425 . dd, k34,is)
426 CALL s8zked3(1,nel,
427 . px(1,3),py(1,3),pz(1,3),px(1,5),py(1,5),pz(1,5),
428 . dd, k35,is)
429 CALL s8zked3(1,nel,
430 . px(1,3),py(1,3),pz(1,3),px(1,6),py(1,6),pz(1,6),
431 . dd, k36,is)
432 CALL s8zked3(1,nel,
433 . px(1,3),py(1,3),pz(1,3),px(1,7),py(1,7),pz(1,7),
434 . dd, k37,is)
435 CALL s8zked3(1,nel,
436 . px(1,3),py(1,3),pz(1,3),px(1,8),py(1,8),pz(1,8),
437 . dd, k38,is)
438 CALL s8zked3(1,nel,
439 . px(1,3),py(1,3),pz(1,3),px(1,9),py(1,9),pz(1,9),
440 . dd, k39,is)
441 CALL s8zked3(1,nel,
442 . px(1,3),py(1,3),pz(1,3),px(1,10),py(1,10),pz(1,10),
443 . dd, k310,is)
444 CALL s8zked3(1,nel,
445 . px(1,3),py(1,3),pz(1,3),px(1,11),py(1,11),pz(1,11),
446 . dd, k311,is)
447 CALL s8zked3(1,nel,
448 . px(1,3),py(1,3),pz(1,3),px(1,12),py(1,12),pz(1,12),
449 . dd, k312,is)
450 CALL s8zked3(1,nel,
451 . px(1,3),py(1,3),pz(1,3),px(1,13),py(1,13),pz(1,13),
452 . dd, k313,is)
453 CALL s8zked3(1,nel,
454 . px(1,3),py(1,3),pz(1,3),px(1,14),py(1,14),pz(1,14),
455 . dd, k314,is)
456 CALL s8zked3(1,nel,
457 . px(1,3),py(1,3),pz(1,3),px(1,15),py(1,15),pz(1,15),
458 . dd, k315,is)
459 CALL s8zked3(1,nel,
460 . px(1,3),py(1,3),pz(1,3),px(1,16),py(1,16),pz(1,16),
461 . dd, k316,is)
462 CALL s8zked3(1,nel,
463 . px(1,3),py(1,3),pz(1,3),px(1,17),py(1,17),pz(1,17),
464 . dd, k317,is)
465 CALL s8zked3(1,nel,
466 . px(1,3),py(1,3),pz(1,3),px(1,18),py(1,18),pz(1,18),
467 . dd, k318,is)
468 CALL s8zked3(1,nel,
469 . px(1,3),py(1,3),pz(1,3),px(1,19),py(1,19),pz(1,19),
470 . dd, k319,is)
471 CALL s8zked3(1,nel,
472 . px(1,3),py(1,3),pz(1,3),px(1,20),py(1,20),pz(1,20),
473 . dd, k320,is)
474
475
476 CALL s8zked3(1,nel,
477 . px(1,4),py(1,4),pz(1,4),px(1,5),py(1,5),pz(1,5),
478 . dd, k45,is)
479 CALL s8zked3(1,nel,
480 . px(1,4),py(1,4),pz(1,4),px(1,6),py(1,6),pz(1,6),
481 . dd, k46,is)
482 CALL s8zked3(1,nel,
483 . px(1,4),py(1,4),pz(1,4),px(1,7),py(1,7),pz(1,7),
484 . dd, k47,is)
485 CALL s8zked3(1,nel,
486 . px(1,4),py(1,4),pz(1,4),px(1,8),py(1,8),pz(1,8),
487 . dd, k48,is)
488 CALL s8zked3(1,nel,
489 . px(1,4),py(1,4),pz(1,4),px(1,9),py(1,9),pz(1,9),
490 . dd, k49,is)
491 CALL s8zked3(1,nel,
492 . px(1,4),py(1,4),pz(1,4),px(1,10),py(1,10),pz(1,10),
493 . dd, k410,is)
494 CALL s8zked3(1,nel,
495 . px(1,4),py(1,4),pz(1,4),px(1,11),py(1,11),pz(1,11),
496 . dd, k411,is)
497 CALL s8zked3(1,nel,
498 . px(1,4),py(1,4),pz(1,4),px(1,12),py(1,12),pz(1,12),
499 . dd, k412,is)
500 CALL s8zked3(1,nel,
501 . px(1,4),py(1,4),pz(1,4),px(1,13),py(1,13),pz(1,13),
502 . dd, k413,is)
503 CALL s8zked3(1,nel,
504 . px(1,4),py(1,4),pz(1,4),px(1,14),py(1,14),pz(1,14),
505 . dd, k414,is)
506 CALL s8zked3(1,nel,
507 . px(1,4),py(1,4),pz(1,4),px(1,15),py(1,15),pz(1,15),
508 . dd, k415,is)
509 CALL s8zked3(1,nel,
510 . px(1,4),py(1,4),pz(1,4),px(1,16),py(1,16),pz(1,16),
511 . dd, k416,is)
512 CALL s8zked3(1,nel,
513 . px(1,4),py(1,4),pz(1,4),px(1,17),py(1,17),pz(1,17),
514 . dd, k417,is)
515 CALL s8zked3(1,nel,
516 . px(1,4),py(1,4),pz(1,4),px(1,18),py(1,18),pz(1,18),
517 . dd, k418,is)
518 CALL s8zked3(1,nel,
519 . px(1,4),py(1,4),pz(1,4),px(1,19),py(1,19),pz(1,19),
520 . dd, k419,is)
521 CALL s8zked3(1,nel,
522 . px(1,4),py(1,4),pz(1,4),px(1,20),py(1,20),pz(1,20),
523 . dd, k420,is)
524
525
526
527 CALL s8zked3(1,nel,
528 . px(1,5),py(1,5),pz(1,5),px(1,6),py(1,6),pz(1,6),
529 . dd, k56,is)
530 CALL s8zked3(1,nel,
531 . px(1,5),py(1,5),pz(1,5),px(1,7),py(1,7),pz(1,7),
532 . dd, k57,is)
533 CALL s8zked3(1,nel,
534 . px(1,5),py(1,5),pz(1,5),px(1,8),py(1,8),pz(1,8),
535 . dd, k58,is)
536 CALL s8zked3(1,nel,
537 . px(1,5),py(1,5),pz(1,5),px(1,9),py(1,9),pz(1,9),
538 . dd, k59,is)
539 CALL s8zked3(1,nel,
540 . px(1,5),py(1,5),pz(1,5),px(1,10),py(1,10),pz(1,10),
541 . dd, k510,is)
542 CALL s8zked3(1,nel,
543 . px(1,5),py(1,5),pz(1,5),px(1,11),py(1,11),pz(1,11),
544 . dd, k511,is)
545 CALL s8zked3(1,nel,
546 . px(1,5),py(1,5),pz(1,5),px(1,12),py(1,12),pz(1,12),
547 . dd, k512,is)
548 CALL s8zked3(1,nel,
549 . px(1,5),py(1,5),pz(1,5),px(1,13),py(1,13),pz(1,13),
550 . dd, k513,is)
551 CALL s8zked3(1,nel,
552 . px(1,5),py(1,5),pz(1,5),px(1,14),py(1,14),pz(1,14),
553 . dd, k514,is)
554 CALL s8zked3(1,nel,
555 . px(1,5),py(1,5),pz(1,5),px(1,15),py(1,15),pz(1,15),
556 . dd, k515,is)
557 CALL s8zked3(1,nel,
558 . px(1,5),py(1,5),pz(1,5),px(1,16),py(1,16),pz(1,16),
559 . dd, k516,is)
560 CALL s8zked3(1,nel,
561 . px(1,5),py(1,5),pz(1,5),px(1,17),py(1,17),pz(1,17),
562 . dd, k517,is)
563 CALL s8zked3(1,nel,
564 . px(1,5),py(1,5),pz(1,5),px(1,18),py(1,18),pz(1,18),
565 . dd, k518,is)
566 CALL s8zked3(1,nel,
567 . px(1,5),py(1,5),pz(1,5),px(1,19),py(1,19),pz(1,19),
568 . dd, k519,is)
569 CALL s8zked3(1,nel,
570 . px(1,5),py(1,5),pz(1,5),px(1,20),py(1,20),pz(1,20),
571 . dd, k520,is)
572
573
574
575 CALL s8zked3(1,nel,
576 . px(1,6),py(1,6),pz(1,6),px(1,7),py(1,7),pz(1,7),
577 . dd, k67,is)
578 CALL s8zked3(1,nel,
579 . px(1,6),py(1,6),pz(1,6),px(1,8),py(1,8),pz(1,8),
580 . dd, k68,is)
581 CALL s8zked3(1,nel,
582 . px(1,6),py(1,6),pz(1,6),px(1,9),py(1,9),pz(1,9),
583 . dd, k69,is)
584 CALL s8zked3(1,nel,
585 . px(1,6),py(1,6),pz(1,6),px(1,10),py(1,10),pz(1,10),
586 . dd, k610,is)
587 CALL s8zked3(1,nel,
588 . px(1,6),py(1,6),pz(1,6),px(1,11),py(1,11),pz(1,11),
589 . dd, k611,is)
590 CALL s8zked3(1,nel,
591 . px(1,6),py(1,6),pz(1,6),px(1,12),py(1,12),pz(1,12),
592 . dd, k612,is)
593 CALL s8zked3(1,nel,
594 . px(1,6),py(1,6),pz(1,6),px(1,13),py(1,13),pz(1,13),
595 . dd, k613,is)
596 CALL s8zked3(1,nel,
597 . px(1,6),py(1,6),pz(1,6),px(1,14),py(1,14),pz(1,14),
598 . dd, k614,is)
599 CALL s8zked3(1,nel,
600 . px(1,6),py(1,6),pz(1,6),px(1,15),py(1,15),pz(1,15),
601 . dd, k615,is)
602 CALL s8zked3(1,nel,
603 . px(1,6),py(1,6),pz(1,6),px(1,16),py(1,16),pz(1,16),
604 . dd, k616,is)
605 CALL s8zked3(1,nel,
606 . px(1,6),py(1,6),pz(1,6),px(1,17),py(1,17),pz(1,17),
607 . dd, k617,is)
608 CALL s8zked3(1,nel,
609 . px(1,6),py(1,6),pz(1,6),px(1,18),py(1,18),pz(1,18),
610 . dd, k618,is)
611 CALL s8zked3(1,nel,
612 . px(1,6),py(1,6),pz(1,6),px(1,19),py(1,19),pz(1,19),
613 . dd, k619,is)
614 CALL s8zked3(1,nel,
615 . px(1,6),py(1,6),pz(1,6),px(1,20),py(1,20),pz(1,20),
616 . dd, k620,is)
617
618
619
620 CALL s8zked3(1,nel,
621 . px(1,7),py(1,7),pz(1,7),px(1,8),py(1,8),pz(1,8),
622 . dd, k78,is)
623 CALL s8zked3(1,nel,
624 . px(1,7),py(1,7),pz(1,7),px(1,9),py(1,9),pz(1,9),
625 . dd, k79,is)
626 CALL s8zked3(1,nel,
627 . px(1,7),py(1,7),pz(1,7),px(1,10),py(1,10),pz(1,10),
628 . dd, k710,is)
629 CALL s8zked3(1,nel,
630 . px(1,7),py(1,7),pz(1,7),px(1,11),py(1,11),pz(1,11),
631 . dd, k711,is)
632 CALL s8zked3(1,nel,
633 . px(1,7),py(1,7),pz(1,7),px(1,12),py(1,12),pz(1,12),
634 . dd, k712,is)
635 CALL s8zked3(1,nel,
636 . px(1,7),py(1,7),pz(1,7),px(1,13),py(1,13),pz(1,13),
637 . dd, k713,is)
638 CALL s8zked3(1,nel,
639 . px(1,7),py(1,7),pz(1,7),px(1,14),py(1,14),pz(1,14),
640 . dd, k714,is)
641 CALL s8zked3(1,nel,
642 . px(1,7),py(1,7),pz(1,7),px(1,15),py(1,15),pz(1,15),
643 . dd, k715,is)
644 CALL s8zked3(1,nel,
645 . px(1,7),py(1,7),pz(1,7),px(1,16),py(1,16),pz(1,16),
646 . dd, k716,is)
647 CALL s8zked3(1,nel,
648 . px(1,7),py(1,7),pz(1,7),px(1,17),py(1,17),pz(1,17),
649 . dd, k717,is)
650 CALL s8zked3(1,nel,
651 . px(1,7),py(1,7),pz(1,7),px(1,18),py(1,18),pz(1,18),
652 . dd, k718,is)
653 CALL s8zked3(1,nel,
654 . px(1,7),py(1,7),pz(1,7),px(1,19),py(1,19),pz(1,19),
655 . dd, k719,is)
656 CALL s8zked3(1,nel,
657 . px(1,7),py(1,7),pz(1,7),px(1,20),py(1,20),pz(1,20),
658 . dd, k720,is)
659
660
661 CALL s8zked3(1,nel,
662 . px(1,8),py(1,8),pz(1,8),px(1,9),py(1,9),pz(1,9),
663 . dd, k89,is)
664 CALL s8zked3(1,nel,
665 . px(1,8),py(1,8),pz(1,8),px(1,10),py(1,10),pz(1,10),
666 . dd, k810,is)
667 CALL s8zked3(1,nel,
668 . px(1,8),py(1,8),pz(1,8),px(1,11),py(1,11),pz(1,11),
669 . dd, k811,is)
670 CALL s8zked3(1,nel,
671 . px(1,8),py(1,8),pz(1,8),px(1,12),py(1,12),pz(1,12),
672 . dd, k812,is)
673 CALL s8zked3(1,nel,
674 . px(1,8),py(1,8),pz(1,8),px(1,13),py(1,13),pz(1,13),
675 . dd, k813,is)
676 CALL s8zked3(1,nel,
677 . px(1,8),py(1,8),pz(1,8),px(1,14),py(1,14),pz(1,14),
678 . dd, k814,is)
679 CALL s8zked3(1,nel,
680 . px(1,8),py(1,8),pz(1,8),px(1,15),py(1,15),pz(1,15),
681 . dd, k815,is)
682 CALL s8zked3(1,nel,
683 . px(1,8),py(1,8),pz(1,8),px(1,16),py(1,16),pz(1,16),
684 . dd, k816,is)
685 CALL s8zked3(1,nel,
686 . px(1,8),py(1,8),pz(1,8),px(1,17),py(1,17),pz(1,17),
687 . dd, k817,is)
688 CALL s8zked3(1,nel,
689 . px(1,8),py(1,8),pz(1,8),px(1,18),py(1,18),pz(1,18),
690 . dd, k818,is)
691 CALL s8zked3(1,nel,
692 . px(1,8),py(1,8),pz(1,8),px(1,19),py(1,19),pz(1,19),
693 . dd, k819,is)
694 CALL s8zked3(1,nel,
695 . px(1,8),py(1,8),pz(1,8),px(1,20),py(1,20),pz(1,20),
696 . dd, k820,is)
697
698
699
700 CALL s8zked3(1,nel,
701 . px(1,9),py(1,9),pz(1,9),px(1,10),py(1,10),pz(1,10),
702 . dd, k910,is)
703 CALL s8zked3(1,nel,
704 . px(1,9),py(1,9),pz(1,9),px(1,11),py(1,11),pz(1,11),
705 . dd, k911,is)
706 CALL s8zked3(1,nel,
707 . px(1,9),py(1,9),pz(1,9),px(1,12),py(1,12),pz(1,12),
708 . dd, k912,is)
709 CALL s8zked3(1,nel,
710 . px(1,9),py(1,9),pz(1,9),px(1,13),py(1,13),pz(1,13),
711 . dd, k913,is)
712 CALL s8zked3(1,nel,
713 . px(1,9),py(1,9),pz(1,9),px(1,14),py(1,14),pz(1,14),
714 . dd, k914,is)
715 CALL s8zked3(1,nel,
716 . px(1,9),py(1,9),pz(1,9),px(1,15),py(1,15),pz(1,15),
717 . dd, k915,is)
718 CALL s8zked3(1,nel,
719 . px(1,9),py(1,9),pz(1,9),px(1,16),py(1,16),pz(1,16),
720 . dd, k916,is)
721 CALL s8zked3(1,nel,
722 . px(1,9),py(1,9),pz(1,9),px(1,17),py(1,17),pz(1,17),
723 . dd, k917,is)
724 CALL s8zked3(1,nel,
725 . px(1,9),py(1,9),pz(1,9),px(1,18),py(1,18),pz(1,18),
726 . dd, k918,is)
727 CALL s8zked3(1,nel,
728 . px(1,9),py(1,9),pz(1,9),px(1,19),py(1,19),pz(1,19),
729 . dd, k919,is)
730 CALL s8zked3(1,nel,
731 . px(1,9),py(1,9),pz(1,9),px(1,20),py(1,20),pz(1,20),
732 . dd, k920,is)
733
734
735
736 CALL s8zked3(1,nel,
737 . px(1,10),py(1,10),pz(1,10),px(1,11),py(1,11),pz(1,11),
738 . dd, k1011,is)
739 CALL s8zked3(1,nel,
740 . px(1,10),py(1,10),pz(1,10),px(1,12),py(1,12),pz(1,12),
741 . dd, k1012,is)
742 CALL s8zked3(1,nel,
743 . px(1,10),py(1,10),pz(1,10),px(1,13),py(1,13),pz(1,13),
744 . dd, k1013,is)
745 CALL s8zked3(1,nel,
746 . px(1,10),py(1,10),pz(1,10),px(1,14),py(1,14),pz(1,14),
747 . dd, k1014,is)
748 CALL s8zked3(1,nel,
749 . px(1,10),py(1,10),pz(1,10),px(1,15),py(1,15),pz(1,15),
750 . dd, k1015,is)
751 CALL s8zked3(1,nel,
752 . px(1,10),py(1,10),pz(1,10),px(1,16),py(1,16),pz(1,16),
753 . dd, k1016,is)
754 CALL s8zked3(1,nel,
755 . px(1,10),py(1,10),pz(1,10),px(1,17),py(1,17),pz(1,17),
756 . dd, k1017,is)
757 CALL s8zked3(1,nel,
758 . px(1,10),py(1,10),pz(1,10),px(1,18),py(1,18),pz(1,18),
759 . dd, k1018,is)
760 CALL s8zked3(1,nel,
761 . px(1,10),py(1,10),pz(1,10),px(1,19),py(1,19),pz(1,19),
762 . dd, k1019,is)
763 CALL s8zked3(1,nel,
764 . px(1,10),py(1,10),pz(1,10),px(1,20),py(1,20),pz(1,20),
765 . dd, k1020,is)
766
767
768
769 CALL s8zked3(1,nel,
770 . px(1,11),py(1,11),pz(1,11),px(1,12),py(1,12),pz(1,12),
771 . dd, k1112,is)
772 CALL s8zked3(1,nel,
773 . px(1,11),py(1,11),pz(1,11),px(1,13),py(1,13),pz(1,13),
774 . dd, k1113,is)
775 CALL s8zked3(1,nel,
776 . px(1,11),py(1,11),pz(1,11),px(1,14),py(1,14),pz(1,14),
777 . dd, k1114,is)
778 CALL s8zked3(1,nel,
779 . px(1,11),py(1,11),pz(1,11),px(1,15),py(1,15),pz(1,15),
780 . dd, k1115,is)
781 CALL s8zked3(1,nel,
782 . px(1,11),py(1,11),pz(1,11),px(1,16),py(1,16),pz(1,16),
783 . dd, k1116,is)
784 CALL s8zked3(1,nel,
785 . px(1,11),py(1,11),pz(1,11),px(1,17),py(1,17),pz(1,17),
786 . dd, k1117,is)
787 CALL s8zked3(1,nel,
788 . px(1,11),py(1,11),pz(1,11),px(1,18),py(1,18),pz(1,18),
789 . dd, k1118,is)
790 CALL s8zked3(1,nel,
791 . px(1,11),py(1,11),pz(1,11),px(1,19),py(1,19),pz(1,19),
792 . dd, k1119,is)
793 CALL s8zked3(1,nel,
794 . px(1,11),py(1,11),pz(1,11),px(1,20),py(1,20),pz(1,20),
795 . dd, k1120,is)
796
797
798 CALL s8zked3(1,nel,
799 . px(1,12),py(1,12),pz(1,12),px(1,13),py(1,13),pz(1,13),
800 . dd, k1213,is)
801 CALL s8zked3(1,nel,
802 . px(1,12),py(1,12),pz(1,12),px(1,14),py(1,14),pz(1,14),
803 . dd, k1214,is)
804 CALL s8zked3(1,nel,
805 . px(1,12),py(1,12),pz(1,12),px(1,15),py(1,15),pz(1,15),
806 . dd, k1215,is)
807 CALL s8zked3(1,nel,
808 . px(1,12),py(1,12),pz(1,12),px(1,16),py(1,16),pz(1,16),
809 . dd, k1216,is)
810 CALL s8zked3(1,nel,
811 . px(1,12),py(1,12),pz(1,12),px(1,17),py(1,17),pz(1,17),
812 . dd, k1217,is)
813 CALL s8zked3(1,nel,
814 . px(1,12),py(1,12),pz(1,12),px(1,18),py(1,18),pz(1,18),
815 . dd, k1218,is)
816 CALL s8zked3(1,nel,
817 . px(1,12),py(1,12),pz(1,12),px(1,19),py(1,19),pz(1,19),
818 . dd, k1219,is)
819 CALL s8zked3(1,nel,
820 . px(1,12),py(1,12),pz(1,12),px(1,20),py(1,20),pz(1,20),
821 . dd, k1220,is)
822
823
824 CALL s8zked3(1,nel,
825 . px(1,13),py(1,13),pz(1,13),px(1,14),py(1,14),pz(1,14),
826 . dd, k1314,is)
827 CALL s8zked3(1,nel,
828 . px(1,13),py(1,13),pz(1,13),px(1,15),py(1,15),pz(1,15),
829 . dd, k1315,is)
830 CALL s8zked3(1,nel,
831 . px(1,13),py(1,13),pz(1,13),px(1,16),py(1,16),pz(1,16),
832 . dd, k1316,is)
833 CALL s8zked3(1,nel,
834 . px(1,13),py(1,13),pz(1,13),px(1,17),py(1,17),pz(1,17),
835 . dd, k1317,is)
836 CALL s8zked3(1,nel,
837 . px(1,13),py(1,13),pz(1,13),px(1,18),py(1,18),pz(1,18),
838 . dd, k1318,is)
839 CALL s8zked3(1,nel,
840 . px(1,13),py(1,13),pz(1,13),px(1,19),py(1,19),pz(1,19),
841 . dd, k1319,is)
842 CALL s8zked3(1,nel,
843 . px(1,13),py(1,13),pz(1,13),px(1,20),py(1,20),pz(1,20),
844 . dd, k1320,is)
845
846
847
848 CALL s8zked3(1,nel,
849 . px(1,14),py(1,14),pz(1,14),px(1,15),py(1,15),pz(1,15),
850 . dd, k1415,is)
851 CALL s8zked3(1,nel,
852 . px(1,14),py(1,14),pz(1,14),px(1,16),py(1,16),pz(1,16),
853 . dd, k1416,is)
854 CALL s8zked3(1,nel,
855 . px(1,14),py(1,14),pz(1,14),px(1,17),py(1,17),pz(1,17),
856 . dd, k1417,is)
857 CALL s8zked3(1,nel,
858 . px(1,14),py(1,14),pz(1,14),px(1,18),py(1,18),pz(1,18),
859 . dd, k1418,is)
860 CALL s8zked3(1,nel,
861 . px(1,14),py(1,14),pz(1,14),px(1,19),py(1,19),pz(1,19),
862 . dd, k1419,is)
863 CALL s8zked3(1,nel,
864 . px(1,14),py(1,14),pz(1,14),px(1,20),py(1,20),pz(1,20),
865 . dd, k1420,is)
866
867
868
869 CALL s8zked3(1,nel,
870 . px(1,15),py(1,15),pz(1,15),px(1,16),py(1,16),pz(1,16),
871 . dd, k1516,is)
872 CALL s8zked3(1,nel,
873 . px(1,15),py(1,15),pz(1,15),px(1,17),py(1,17),pz(1,17),
874 . dd, k1517,is)
875 CALL s8zked3(1,nel,
876 . px(1,15),py(1,15),pz(1,15),px(1,18),py(1,18),pz(1,18),
877 . dd, k1518,is)
878 CALL s8zked3(1,nel,
879 . px(1,15),py(1,15),pz(1,15),px(1,19),py(1,19),pz(1,19),
880 . dd, k1519,is)
881 CALL s8zked3(1,nel,
882 . px(1,15),py(1,15),pz(1,15),px(1,20),py(1,20),pz(1,20),
883 . dd, k1520,is)
884
885
886
887 CALL s8zked3(1,nel,
888 . px(1,16),py(1,16),pz(1,16),px(1,17),py(1,17),pz(1,17),
889 . dd, k1617,is)
890 CALL s8zked3(1,nel,
891 . px(1,16),py(1,16),pz(1,16),px(1,18),py(1,18),pz(1,18),
892 . dd, k1618,is)
893 CALL s8zked3(1,nel,
894 . px(1,16),py(1,16),pz(1,16),px(1,19),py(1,19),pz(1,19),
895 . dd, k1619,is)
896 CALL s8zked3(1,nel,
897 . px(1,16),py(1,16),pz(1,16),px(1,20),py(1,20),pz(1,20),
898 . dd, k1620,is)
899
900
901
902 CALL s8zked3(1,nel,
903 . px(1,17),py(1,17),pz(1,17),px(1,18),py(1,18),pz(1,18),
904 . dd, k1718,is)
905 CALL s8zked3(1,nel,
906 . px(1,17),py(1,17),pz(1,17),px(1,19),py(1,19),pz(1,19),
907 . dd, k1719,is)
908 CALL s8zked3(1,nel,
909 . px(1,17),py(1,17),pz(1,17),px(1,20),py(1,20),pz(1,20),
910 . dd, k1720,is)
911
912
913
914 CALL s8zked3(1,nel,
915 . px(1,18),py(1,18),pz(1,18),px(1,19),py(1,19),pz(1,19),
916 . dd, k1819,is)
917 CALL s8zked3(1,nel,
918 . px(1,18),py(1,18),pz(1,18),px(1,20),py(1,20),pz(1,20),
919 . dd, k1820,is)
920
921
922 CALL s8zked3(1,nel,
923 . px(1,19),py(1,19),pz(1,19),px(1,20),py(1,20),pz(1,20),
924 . dd, k1920,is)
925
926
927
928C
929 IF (iksup==0) THEN
930 CALL s8zkeg3(1,nel,
931 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
932 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
933 . gg,k12,is)
934 CALL s8zkeg3(1,nel,
935 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
936 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
937 . gg,k13,is)
938 CALL s8zkeg3(1,nel,
939 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
940 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
941 . gg,k14,is)
942 CALL s8zkeg3(1,nel,
943 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
944 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
945 . gg,k15,is)
946 CALL s8zkeg3(1,nel,
947 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
948 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
949 . gg,k16,is)
950 CALL s8zkeg3(1,nel,
951 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
952 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
953 . gg,k17,is)
954 CALL s8zkeg3(1,nel,
955 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
956 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
957 . gg,k18,is)
958 CALL s8zkeg3(1,nel,
959 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
960 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
961 . gg,k19,is)
962 CALL s8zkeg3(1,nel,
963 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
964 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
965 . gg,k110,is)
966 CALL s8zkeg3(1,nel,
967 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
968 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
969 . gg,k111,is)
970 CALL s8zkeg3(1,nel,
971 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
972 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
973 . gg,k112,is)
974 CALL s8zkeg3(1,nel,
975 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
976 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
977 . gg,k113,is)
978 CALL s8zkeg3(1,nel,
979 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
980 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
981 . gg,k114,is)
982 CALL s8zkeg3(1,nel,
983 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
984 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
985 . gg,k115,is)
986 CALL s8zkeg3(1,nel,
987 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
988 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
989 . gg,k116,is)
990 CALL s8zkeg3(1,nel,
991 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
992 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
993 . gg,k117,is)
994 CALL s8zkeg3(1,nel,
995 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
996 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
997 . gg,k118,is)
998 CALL s8zkeg3(1,nel,
999 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
1000 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1001 . gg,k119,is)
1002 CALL s8zkeg3(1,nel,
1003 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
1004 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1005 . gg,k120,is)
1006
1007
1008
1009
1010 CALL s8zkeg3(1,nel,
1011 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1012 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1013 . gg,k23,is)
1014 CALL s8zkeg3(1,nel,
1015 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1016 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1017 . gg,k24,is)
1018 CALL s8zkeg3(1,nel,
1019 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1020 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1021 . gg,k25,is)
1022 CALL s8zkeg3(1,nel,
1023 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1024 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1025 . gg,k26,is)
1026 CALL s8zkeg3(1,nel,
1027 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1028 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1029 . gg,k27,is)
1030 CALL s8zkeg3(1,nel,
1031 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1032 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1033 . gg,k28,is)
1034 CALL s8zkeg3(1,nel,
1035 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1036 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1037 . gg,k29,is)
1038 CALL s8zkeg3(1,nel,
1039 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1040 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1041 . gg,k210,is)
1042 CALL s8zkeg3(1,nel,
1043 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1044 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1045 . gg,k211,is)
1046 CALL s8zkeg3(1,nel,
1047 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1048 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
1049 . gg,k212,is)
1050 CALL s8zkeg3(1,nel,
1051 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1052 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
1053 . gg,k213,is)
1054 CALL s8zkeg3(1,nel,
1055 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1056 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
1057 . gg,k214,is)
1058 CALL s8zkeg3(1,nel,
1059 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1060 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
1061 . gg,k215,is)
1062 CALL s8zkeg3(1,nel,
1063 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1064 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
1065 . gg,k216,is)
1066 CALL s8zkeg3(1,nel,
1067 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1068 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
1069 . gg,k217,is)
1070 CALL s8zkeg3(1,nel,
1071 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1072 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
1073 . gg,k218,is)
1074 CALL s8zkeg3(1,nel,
1075 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1076 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1077 . gg,k219,is)
1078 CALL s8zkeg3(1,nel,
1079 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1080 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1081 . gg,k220,is)
1082
1083
1084
1085
1086 CALL s8zkeg3(1,nel,
1087 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1088 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1089 . gg,k34,is)
1090 CALL s8zkeg3(1,nel,
1091 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1092 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1093 . gg,k35,is)
1094 CALL s8zkeg3(1,nel,
1095 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1096 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1097 . gg,k36,is)
1098 CALL s8zkeg3(1,nel,
1099 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1100 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1101 . gg,k37,is)
1102 CALL s8zkeg3(1,nel,
1103 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1104 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1105 . gg,k38,is)
1106 CALL s8zkeg3(1,nel,
1107 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1108 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1109 . gg,k39,is)
1110 CALL s8zkeg3(1,nel,
1111 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1112 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1113 . gg,k310,is)
1114 CALL s8zkeg3(1,nel,
1115 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1116 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1117 . gg,k311,is)
1118 CALL s8zkeg3(1,nel,
1119 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1120 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
1121 . gg,k312,is)
1122 CALL s8zkeg3(1,nel,
1123 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1124 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
1125 . gg,k313,is)
1126 CALL s8zkeg3(1,nel,
1127 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1128 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
1129 . gg,k314,is)
1130 CALL s8zkeg3(1,nel,
1131 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1132 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
1133 . gg,k315,is)
1134 CALL s8zkeg3(1,nel,
1135 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1136 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
1137 . gg,k316,is)
1138 CALL s8zkeg3(1,nel,
1139 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1140 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
1141 . gg,k317,is)
1142 CALL s8zkeg3(1,nel,
1143 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1144 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
1145 . gg,k318,is)
1146 CALL s8zkeg3(1,nel,
1147 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1148 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1149 . gg,k319,is)
1150 CALL s8zkeg3(1,nel,
1151 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1152 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1153 . gg,k320,is)
1154
1155
1156
1157
1158 CALL s8zkeg3(1,nel,
1159 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1160 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1161 . gg,k45,is)
1162 CALL s8zkeg3(1,nel,
1163 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1164 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1165 . gg,k46,is)
1166 CALL s8zkeg3(1,nel,
1167 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1168 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1169 . gg,k47,is)
1170 CALL s8zkeg3(1,nel,
1171 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1172 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1173 . gg,k48,is)
1174 CALL s8zkeg3(1,nel,
1175 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1176 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1177 . gg,k49,is)
1178 CALL s8zkeg3(1,nel,
1179 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1180 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1181 . gg,k410,is)
1182 CALL s8zkeg3(1,nel,
1183 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1184 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1185 . gg,k411,is)
1186 CALL s8zkeg3(1,nel,
1187 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1188 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
1189 . gg,k412,is)
1190 CALL s8zkeg3(1,nel,
1191 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1192 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
1193 . gg,k413,is)
1194 CALL s8zkeg3(1,nel,
1195 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1196 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
1197 . gg,k414,is)
1198 CALL s8zkeg3(1,nel,
1199 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1200 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
1201 . gg,k415,is)
1202 CALL s8zkeg3(1,nel,
1203 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1204 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
1205 . gg,k416,is)
1206 CALL s8zkeg3(1,nel,
1207 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1208 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
1209 . gg,k417,is)
1210 CALL s8zkeg3(1,nel,
1211 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1212 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
1213 . gg,k418,is)
1214 CALL s8zkeg3(1,nel,
1215 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1216 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1217 . gg,k419,is)
1218 CALL s8zkeg3(1,nel,
1219 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1220 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1221 . gg,k420,is)
1222
1223
1224
1225
1226 CALL s8zkeg3(1,nel,
1227 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1228 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1229 . gg,k56,is)
1230 CALL s8zkeg3(1,nel,
1231 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1232 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1233 . gg,k57,is)
1234 CALL s8zkeg3(1,nel,
1235 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1236 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1237 . gg,k58,is)
1238 CALL s8zkeg3(1,nel,
1239 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1240 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1241 . gg,k59,is)
1242 CALL s8zkeg3(1,nel,
1243 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1244 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1245 . gg,k510,is)
1246 CALL s8zkeg3(1,nel,
1247 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1248 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1249 . gg,k511,is)
1250 CALL s8zkeg3(1,nel,
1251 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1252 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
1253 . gg,k512,is)
1254 CALL s8zkeg3(1,nel,
1255 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1256 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
1257 . gg,k513,is)
1258 CALL s8zkeg3(1,nel,
1259 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1260 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
1261 . gg,k514,is)
1262 CALL s8zkeg3(1,nel,
1263 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1264 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
1265 . gg,k515,is)
1266 CALL s8zkeg3(1,nel,
1267 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1268 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
1269 . gg,k516,is)
1270 CALL s8zkeg3(1,nel,
1271 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1272 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
1273 . gg,k517,is)
1274 CALL s8zkeg3(1,nel,
1275 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1276 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
1277 . gg,k518,is)
1278 CALL s8zkeg3(1,nel,
1279 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1280 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1281 . gg,k519,is)
1282 CALL s8zkeg3(1,nel,
1283 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1284 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1285 . gg,k520,is)
1286
1287
1288
1289 CALL s8zkeg3(1,nel,
1290 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1291 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1292 . gg,k67,is)
1293 CALL s8zkeg3(1,nel,
1294 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1295 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1296 . gg,k68,is)
1297 CALL s8zkeg3(1,nel,
1298 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1299 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1300 . gg,k69,is)
1301 CALL s8zkeg3(1,nel,
1302 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1303 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1304 . gg,k610,is)
1305 CALL s8zkeg3(1,nel,
1306 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1307 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1308 . gg,k611,is)
1309 CALL s8zkeg3(1,nel,
1310 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1311 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
1312 . gg,k612,is)
1313 CALL s8zkeg3(1,nel,
1314 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1315 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
1316 . gg,k613,is)
1317 CALL s8zkeg3(1,nel,
1318 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1319 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
1320 . gg,k614,is)
1321 CALL s8zkeg3(1,nel,
1322 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1323 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
1324 . gg,k615,is)
1325 CALL s8zkeg3(1,nel,
1326 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1327 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
1328 . gg,k616,is)
1329 CALL s8zkeg3(1,nel,
1330 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1331 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
1332 . gg,k617,is)
1333 CALL s8zkeg3(1,nel,
1334 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1335 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
1336 . gg,k618,is)
1337 CALL s8zkeg3(1,nel,
1338 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1339 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1340 . gg,k619,is)
1341 CALL s8zkeg3(1,nel,
1342 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1343 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1344 . gg,k620,is)
1345
1346
1347
1348 CALL s8zkeg3(1,nel,
1349 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1350 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1351 . gg,k78,is)
1352 CALL s8zkeg3(1,nel,
1353 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1354 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1355 . gg,k79,is)
1356 CALL s8zkeg3(1,nel,
1357 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1358 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1359 . gg,k710,is)
1360 CALL s8zkeg3(1,nel,
1361 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1362 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1363 . gg,k711,is)
1364 CALL s8zkeg3(1,nel,
1365 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1366 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
1367 . gg,k712,is)
1368 CALL s8zkeg3(1,nel,
1369 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1370 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
1371 . gg,k713,is)
1372 CALL s8zkeg3(1,nel,
1373 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1374 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
1375 . gg,k714,is)
1376 CALL s8zkeg3(1,nel,
1377 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1378 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
1379 . gg,k715,is)
1380 CALL s8zkeg3(1,nel,
1381 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1382 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
1383 . gg,k716,is)
1384 CALL s8zkeg3(1,nel,
1385 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1386 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
1387 . gg,k717,is)
1388 CALL s8zkeg3(1,nel,
1389 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1390 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
1391 . gg,k718,is)
1392 CALL s8zkeg3(1,nel,
1393 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1394 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1395 . gg,k719,is)
1396 CALL s8zkeg3(1,nel,
1397 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1398 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1399 . gg,k720,is)
1400
1401
1402
1403 CALL s8zkeg3(1,nel,
1404 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1405 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1406 . gg,k89,is)
1407 CALL s8zkeg3(1,nel,
1408 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1409 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1410 . gg,k810,is)
1411 CALL s8zkeg3(1,nel,
1412 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1413 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1414 . gg,k811,is)
1415 CALL s8zkeg3(1,nel,
1416 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1417 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
1418 . gg,k812,is)
1419 CALL s8zkeg3(1,nel,
1420 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1421 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
1422 . gg,k813,is)
1423 CALL s8zkeg3(1,nel,
1424 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1425 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
1426 . gg,k814,is)
1427 CALL s8zkeg3(1,nel,
1428 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1429 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
1430 . gg,k815,is)
1431 CALL s8zkeg3(1,nel,
1432 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1433 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
1434 . gg,k816,is)
1435 CALL s8zkeg3(1,nel,
1436 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1437 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
1438 . gg,k817,is)
1439 CALL s8zkeg3(1,nel,
1440 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1441 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
1442 . gg,k818,is)
1443 CALL s8zkeg3(1,nel,
1444 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1445 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1446 . gg,k819,is)
1447 CALL s8zkeg3(1,nel,
1448 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1449 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1450 . gg,k820,is)
1451
1452
1453
1454 CALL s8zkeg3(1,nel,
1455 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1456 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1457 . gg,k910,is)
1458 CALL s8zkeg3(1,nel,
1459 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1460 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1461 . gg,k911,is)
1462 CALL s8zkeg3(1,nel,
1463 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1464 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
1465 . gg,k912,is)
1466 CALL s8zkeg3(1,nel,
1467 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1468 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
1469 . gg,k913,is)
1470 CALL s8zkeg3(1,nel,
1471 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1472 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
1473 . gg,k914,is)
1474 CALL s8zkeg3(1,nel,
1475 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1476 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
1477 . gg,k915,is)
1478 CALL s8zkeg3(1,nel,
1479 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1480 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
1481 . gg,k916,is)
1482 CALL s8zkeg3(1,nel,
1483 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1484 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
1485 . gg,k917,is)
1486 CALL s8zkeg3(1,nel,
1487 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1488 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
1489 . gg,k918,is)
1490 CALL s8zkeg3(1,nel,
1491 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1492 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1493 . gg,k919,is)
1494 CALL s8zkeg3(1,nel,
1495 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1496 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1497 . gg,k920,is)
1498
1499
1500
1501 CALL s8zkeg3(1,nel,
1502 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1503 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1504 . gg,k1011,is)
1505 CALL s8zkeg3(1,nel,
1506 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1507 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
1508 . gg,k1012,is)
1509 CALL s8zkeg3(1,nel,
1510 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1511 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
1512 . gg,k1013,is)
1513 CALL s8zkeg3(1,nel,
1514 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1515 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
1516 . gg,k1014,is)
1517 CALL s8zkeg3(1,nel,
1518 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1519 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
1520 . gg,k1015,is)
1521 CALL s8zkeg3(1,nel,
1522 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1523 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
1524 . gg,k1016,is)
1525 CALL s8zkeg3(1,nel,
1526 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1527 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
1528 . gg,k1017,is)
1529 CALL s8zkeg3(1,nel,
1530 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1531 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
1532 . gg,k1018,is)
1533 CALL s8zkeg3(1,nel,
1534 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1535 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1536 . gg,k1019,is)
1537 CALL s8zkeg3(1,nel,
1538 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1539 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1540 . gg,k1020,is)
1541
1542
1543 CALL s8zkeg3(1,nel,
1544 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1545 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
1546 . gg,k1112,is)
1547 CALL s8zkeg3(1,nel,
1548 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1549 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
1550 . gg,k1113,is)
1551 CALL s8zkeg3(1,nel,
1552 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1553 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
1554 . gg,k1114,is)
1555 CALL s8zkeg3(1,nel,
1556 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1557 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
1558 . gg,k1115,is)
1559 CALL s8zkeg3(1,nel,
1560 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1561 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
1562 . gg,k1116,is)
1563 CALL s8zkeg3(1,nel,
1564 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1565 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
1566 . gg,k1117,is)
1567 CALL s8zkeg3(1,nel,
1568 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1569 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
1570 . gg,k1118,is)
1571 CALL s8zkeg3(1,nel,
1572 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1573 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1574 . gg,k1119,is)
1575 CALL s8zkeg3(1,nel,
1576 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1577 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1578 . gg,k1120,is)
1579
1580
1581 CALL s8zkeg3(1,nel,
1582 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
1583 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
1584 . gg,k1213,is)
1585 CALL s8zkeg3(1,nel,
1586 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
1587 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
1588 . gg,k1214,is)
1589 CALL s8zkeg3(1,nel,
1590 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
1591 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
1592 . gg,k1215,is)
1593 CALL s8zkeg3(1,nel,
1594 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
1595 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
1596 . gg,k1216,is)
1597 CALL s8zkeg3(1,nel,
1598 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
1599 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
1600 . gg,k1217,is)
1601 CALL s8zkeg3(1,nel,
1602 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
1603 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
1604 . gg,k1218,is)
1605 CALL s8zkeg3(1,nel,
1606 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
1607 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1608 . gg,k1219,is)
1609 CALL s8zkeg3(1,nel,
1610 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
1611 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1612 . gg,k1220,is)
1613
1614
1615 CALL s8zkeg3(1,nel,
1616 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
1617 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
1618 . gg,k1314,is)
1619 CALL s8zkeg3(1,nel,
1620 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
1621 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
1622 . gg,k1315,is)
1623 CALL s8zkeg3(1,nel,
1624 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
1625 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
1626 . gg,k1316,is)
1627 CALL s8zkeg3(1,nel,
1628 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
1629 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
1630 . gg,k1317,is)
1631 CALL s8zkeg3(1,nel,
1632 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
1633 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
1634 . gg,k1318,is)
1635 CALL s8zkeg3(1,nel,
1636 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
1637 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1638 . gg,k1319,is)
1639 CALL s8zkeg3(1,nel,
1640 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
1641 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1642 . gg,k1320,is)
1643
1644
1645
1646 CALL s8zkeg3(1,nel,
1647 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
1648 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
1649 . gg,k1415,is)
1650 CALL s8zkeg3(1,nel,
1651 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
1652 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
1653 . gg,k1416,is)
1654 CALL s8zkeg3(1,nel,
1655 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
1656 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
1657 . gg,k1417,is)
1658 CALL s8zkeg3(1,nel,
1659 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
1660 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
1661 . gg,k1418,is)
1662 CALL s8zkeg3(1,nel,
1663 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
1664 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1665 . gg,k1419,is)
1666 CALL s8zkeg3(1,nel,
1667 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
1668 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1669 . gg,k1420,is)
1670
1671
1672 CALL s8zkeg3(1,nel,
1673 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
1674 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
1675 . gg,k1516,is)
1676 CALL s8zkeg3(1,nel,
1677 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
1678 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
1679 . gg,k1517,is)
1680 CALL s8zkeg3(1,nel,
1681 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
1682 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
1683 . gg,k1518,is)
1684 CALL s8zkeg3(1,nel,
1685 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
1686 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1687 . gg,k1519,is)
1688 CALL s8zkeg3(1,nel,
1689 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
1690 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1691 . gg,k1520,is)
1692
1693
1694 CALL s8zkeg3(1,nel,
1695 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
1696 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
1697 . gg,k1617,is)
1698 CALL s8zkeg3(1,nel,
1699 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
1700 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
1701 . gg,k1618,is)
1702 CALL s8zkeg3(1,nel,
1703 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
1704 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1705 . gg,k1619,is)
1706 CALL s8zkeg3(1,nel,
1707 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
1708 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1709 . gg,k1620,is)
1710
1711
1712 CALL s8zkeg3(1,nel,
1713 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
1714 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
1715 . gg,k1718,is)
1716 CALL s8zkeg3(1,nel,
1717 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
1718 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1719 . gg,k1719,is)
1720 CALL s8zkeg3(1,nel,
1721 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
1722 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1723 . gg,k1720,is)
1724
1725
1726 CALL s8zkeg3(1,nel,
1727 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
1728 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1729 . gg,k1819,is)
1730 CALL s8zkeg3(1,nel,
1731 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
1732 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1733 . gg,k1820,is)
1734
1735
1736 CALL s8zkeg3(1,nel,
1737 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1738 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1739 . gg,k1920,is)
1740
1741
1742 ELSE
1743 CALL s8zkeg33(1,nel,
1744 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
1745 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1746 . g33,k12,is)
1747 CALL s8zkeg33(1,nel,
1748 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
1749 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1750 . g33,k13,is)
1751 CALL s8zkeg33(1,nel,
1752 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
1753 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1754 . g33,k14,is)
1755 CALL s8zkeg33(1,nel,
1756 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
1757 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1758 . g33,k15,is)
1759 CALL s8zkeg33(1,nel,
1760 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
1761 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1762 . g33,k16,is)
1763 CALL s8zkeg33(1,nel,
1764 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
1765 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1766 . g33,k17,is)
1767 CALL s8zkeg33(1,nel,
1768 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
1769 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1770 . g33,k18,is)
1771 CALL s8zkeg33(1,nel,
1772 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
1773 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1774 . g33,k19,is)
1775 CALL s8zkeg33(1,nel,
1776 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
1777 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1778 . g33,k110,is)
1779 CALL s8zkeg33(1,nel,
1780 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
1781 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1782 . g33,k111,is)
1783 CALL s8zkeg33(1,nel,
1784 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
1785 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
1786 . g33,k112,is)
1787 CALL s8zkeg33(1,nel,
1788 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
1789 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
1790 . g33,k113,is)
1791 CALL s8zkeg33(1,nel,
1792 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
1793 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
1794 . g33,k114,is)
1795 CALL s8zkeg33(1,nel,
1796 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
1797 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
1798 . g33,k115,is)
1799 CALL s8zkeg33(1,nel,
1800 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
1801 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
1802 . g33,k116,is)
1803 CALL s8zkeg33(1,nel,
1804 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
1805 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
1806 . g33,k117,is)
1807 CALL s8zkeg33(1,nel,
1808 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
1809 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
1810 . g33,k118,is)
1811 CALL s8zkeg33(1,nel,
1812 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
1813 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1814 . g33,k119,is)
1815 CALL s8zkeg33(1,nel,
1816 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
1817 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1818 . g33,k120,is)
1819
1820
1821
1822 CALL s8zkeg33(1,nel,
1823 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1824 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1825 . g33,k23,is)
1826 CALL s8zkeg33(1,nel,
1827 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1828 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1829 . g33,k24,is)
1830 CALL s8zkeg33(1,nel,
1831 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1832 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1833 . g33,k25,is)
1834 CALL s8zkeg33(1,nel,
1835 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1836 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1837 . g33,k26,is)
1838 CALL s8zkeg33(1,nel,
1839 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1840 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1841 . g33,k27,is)
1842 CALL s8zkeg33(1,nel,
1843 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1844 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1845 . g33,k28,is)
1846 CALL s8zkeg33(1,nel,
1847 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1848 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1849 . g33,k29,is)
1850 CALL s8zkeg33(1,nel,
1851 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1852 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1853 . g33,k210,is)
1854 CALL s8zkeg33(1,nel,
1855 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1856 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1857 . g33,k211,is)
1858 CALL s8zkeg33(1,nel,
1859 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1860 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
1861 . g33,k212,is)
1862 CALL s8zkeg33(1,nel,
1863 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1864 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
1865 . g33,k213,is)
1866 CALL s8zkeg33(1,nel,
1867 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1868 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
1869 . g33,k214,is)
1870 CALL s8zkeg33(1,nel,
1871 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1872 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
1873 . g33,k215,is)
1874 CALL s8zkeg33(1,nel,
1875 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1876 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
1877 . g33,k216,is)
1878 CALL s8zkeg33(1,nel,
1879 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1880 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
1881 . g33,k217,is)
1882 CALL s8zkeg33(1,nel,
1883 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1884 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
1885 . g33,k218,is)
1886 CALL s8zkeg33(1,nel,
1887 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1888 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1889 . g33,k219,is)
1890 CALL s8zkeg33(1,nel,
1891 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
1892 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1893 . g33,k220,is)
1894
1895
1896
1897 CALL s8zkeg33(1,nel,
1898 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1899 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1900 . g33,k34,is)
1901 CALL s8zkeg33(1,nel,
1902 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1903 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1904 . g33,k35,is)
1905 CALL s8zkeg33(1,nel,
1906 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1907 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1908 . g33,k36,is)
1909 CALL s8zkeg33(1,nel,
1910 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1911 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1912 . g33,k37,is)
1913 CALL s8zkeg33(1,nel,
1914 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1915 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1916 . g33,k38,is)
1917 CALL s8zkeg33(1,nel,
1918 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1919 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1920 . g33,k39,is)
1921 CALL s8zkeg33(1,nel,
1922 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1923 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1924 . g33,k310,is)
1925 CALL s8zkeg33(1,nel,
1926 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1927 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1928 . g33,k311,is)
1929 CALL s8zkeg33(1,nel,
1930 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1931 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
1932 . g33,k312,is)
1933 CALL s8zkeg33(1,nel,
1934 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1935 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
1936 . g33,k313,is)
1937 CALL s8zkeg33(1,nel,
1938 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1939 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
1940 . g33,k314,is)
1941 CALL s8zkeg33(1,nel,
1942 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1943 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
1944 . g33,k315,is)
1945 CALL s8zkeg33(1,nel,
1946 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1947 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
1948 . g33,k316,is)
1949 CALL s8zkeg33(1,nel,
1950 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1951 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
1952 . g33,k317,is)
1953 CALL s8zkeg33(1,nel,
1954 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1955 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
1956 . g33,k318,is)
1957 CALL s8zkeg33(1,nel,
1958 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1959 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
1960 . g33,k319,is)
1961 CALL s8zkeg33(1,nel,
1962 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
1963 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
1964 . g33,k320,is)
1965
1966
1967
1968
1969 CALL s8zkeg33(1,nel,
1970 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1971 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
1972 . g33,k45,is)
1973 CALL s8zkeg33(1,nel,
1974 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1975 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
1976 . g33,k46,is)
1977 CALL s8zkeg33(1,nel,
1978 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1979 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
1980 . g33,k47,is)
1981 CALL s8zkeg33(1,nel,
1982 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1983 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
1984 . g33,k48,is)
1985 CALL s8zkeg33(1,nel,
1986 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1987 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
1988 . g33,k49,is)
1989 CALL s8zkeg33(1,nel,
1990 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1991 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
1992 . g33,k410,is)
1993 CALL s8zkeg33(1,nel,
1994 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1995 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
1996 . g33,k411,is)
1997 CALL s8zkeg33(1,nel,
1998 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
1999 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
2000 . g33,k412,is)
2001 CALL s8zkeg33(1,nel,
2002 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
2003 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
2004 . g33,k413,is)
2005 CALL s8zkeg33(1,nel,
2006 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
2007 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
2008 . g33,k414,is)
2009 CALL s8zkeg33(1,nel,
2010 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
2011 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
2012 . g33,k415,is)
2013 CALL s8zkeg33(1,nel,
2014 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
2015 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
2016 . g33,k416,is)
2017 CALL s8zkeg33(1,nel,
2018 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
2019 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
2020 . g33,k417,is)
2021 CALL s8zkeg33(1,nel,
2022 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
2023 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
2024 . g33,k418,is)
2025 CALL s8zkeg33(1,nel,
2026 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
2027 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
2028 . g33,k419,is)
2029 CALL s8zkeg33(1,nel,
2030 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
2031 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
2032 . g33,k420,is)
2033
2034
2035
2036 CALL s8zkeg33(1,nel,
2037 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
2038 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
2039 . g33,k56,is)
2040 CALL s8zkeg33(1,nel,
2041 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
2042 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
2043 . g33,k57,is)
2044 CALL s8zkeg33(1,nel,
2045 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
2046 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
2047 . g33,k58,is)
2048 CALL s8zkeg33(1,nel,
2049 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
2050 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
2051 . g33,k59,is)
2052 CALL s8zkeg33(1,nel,
2053 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
2054 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
2055 . g33,k510,is)
2056 CALL s8zkeg33(1,nel,
2057 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
2058 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
2059 . g33,k511,is)
2060 CALL s8zkeg33(1,nel,
2061 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
2062 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
2063 . g33,k512,is)
2064 CALL s8zkeg33(1,nel,
2065 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
2066 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
2067 . g33,k513,is)
2068 CALL s8zkeg33(1,nel,
2069 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
2070 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
2071 . g33,k514,is)
2072 CALL s8zkeg33(1,nel,
2073 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
2074 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
2075 . g33,k515,is)
2076 CALL s8zkeg33(1,nel,
2077 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
2078 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
2079 . g33,k516,is)
2080 CALL s8zkeg33(1,nel,
2081 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
2082 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
2083 . g33,k517,is)
2084 CALL s8zkeg33(1,nel,
2085 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
2086 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
2087 . g33,k518,is)
2088 CALL s8zkeg33(1,nel,
2089 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
2090 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
2091 . g33,k519,is)
2092 CALL s8zkeg33(1,nel,
2093 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
2094 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
2095 . g33,k520,is)
2096
2097
2098
2099 CALL s8zkeg33(1,nel,
2100 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
2101 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
2102 . g33,k67,is)
2103 CALL s8zkeg33(1,nel,
2104 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
2105 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
2106 . g33,k68,is)
2107 CALL s8zkeg33(1,nel,
2108 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
2109 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
2110 . g33,k69,is)
2111 CALL s8zkeg33(1,nel,
2112 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
2113 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
2114 . g33,k610,is)
2115 CALL s8zkeg33(1,nel,
2116 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
2117 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
2118 . g33,k611,is)
2119 CALL s8zkeg33(1,nel,
2120 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
2121 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
2122 . g33,k612,is)
2123 CALL s8zkeg33(1,nel,
2124 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
2125 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
2126 . g33,k613,is)
2127 CALL s8zkeg33(1,nel,
2128 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
2129 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
2130 . g33,k614,is)
2131 CALL s8zkeg33(1,nel,
2132 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
2133 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
2134 . g33,k615,is)
2135 CALL s8zkeg33(1,nel,
2136 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
2137 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
2138 . g33,k616,is)
2139 CALL s8zkeg33(1,nel,
2140 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
2141 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
2142 . g33,k617,is)
2143 CALL s8zkeg33(1,nel,
2144 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
2145 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
2146 . g33,k618,is)
2147 CALL s8zkeg33(1,nel,
2148 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
2149 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
2150 . g33,k619,is)
2151 CALL s8zkeg33(1,nel,
2152 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
2153 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
2154 . g33,k620,is)
2155
2156
2157
2158 CALL s8zkeg33(1,nel,
2159 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
2160 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
2161 . g33,k78,is)
2162 CALL s8zkeg33(1,nel,
2163 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
2164 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
2165 . g33,k79,is)
2166 CALL s8zkeg33(1,nel,
2167 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
2168 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
2169 . g33,k710,is)
2170 CALL s8zkeg33(1,nel,
2171 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
2172 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
2173 . g33,k711,is)
2174 CALL s8zkeg33(1,nel,
2175 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
2176 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
2177 . g33,k712,is)
2178 CALL s8zkeg33(1,nel,
2179 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
2180 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
2181 . g33,k713,is)
2182 CALL s8zkeg33(1,nel,
2183 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
2184 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
2185 . g33,k714,is)
2186 CALL s8zkeg33(1,nel,
2187 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
2188 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
2189 . g33,k715,is)
2190 CALL s8zkeg33(1,nel,
2191 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
2192 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
2193 . g33,k716,is)
2194 CALL s8zkeg33(1,nel,
2195 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
2196 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
2197 . g33,k717,is)
2198 CALL s8zkeg33(1,nel,
2199 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
2200 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
2201 . g33,k718,is)
2202 CALL s8zkeg33(1,nel,
2203 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
2204 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
2205 . g33,k719,is)
2206 CALL s8zkeg33(1,nel,
2207 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
2208 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
2209 . g33,k720,is)
2210
2211
2212
2213 CALL s8zkeg33(1,nel,
2214 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
2215 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
2216 . g33,k89,is)
2217 CALL s8zkeg33(1,nel,
2218 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
2219 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
2220 . g33,k810,is)
2221 CALL s8zkeg33(1,nel,
2222 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
2223 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
2224 . g33,k811,is)
2225 CALL s8zkeg33(1,nel,
2226 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
2227 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
2228 . g33,k812,is)
2229 CALL s8zkeg33(1,nel,
2230 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
2231 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
2232 . g33,k813,is)
2233 CALL s8zkeg33(1,nel,
2234 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
2235 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
2236 . g33,k814,is)
2237 CALL s8zkeg33(1,nel,
2238 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
2239 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
2240 . g33,k815,is)
2241 CALL s8zkeg33(1,nel,
2242 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
2243 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
2244 . g33,k816,is)
2245 CALL s8zkeg33(1,nel,
2246 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
2247 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
2248 . g33,k817,is)
2249 CALL s8zkeg33(1,nel,
2250 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
2251 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
2252 . g33,k818,is)
2253 CALL s8zkeg33(1,nel,
2254 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
2255 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
2256 . g33,k819,is)
2257 CALL s8zkeg33(1,nel,
2258 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
2259 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
2260 . g33,k820,is)
2261
2262
2263
2264 CALL s8zkeg33(1,nel,
2265 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
2266 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
2267 . g33,k910,is)
2268 CALL s8zkeg33(1,nel,
2269 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
2270 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
2271 . g33,k911,is)
2272 CALL s8zkeg33(1,nel,
2273 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
2274 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
2275 . g33,k912,is)
2276 CALL s8zkeg33(1,nel,
2277 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
2278 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
2279 . g33,k913,is)
2280 CALL s8zkeg33(1,nel,
2281 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
2282 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
2283 . g33,k914,is)
2284 CALL s8zkeg33(1,nel,
2285 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
2286 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
2287 . g33,k915,is)
2288 CALL s8zkeg33(1,nel,
2289 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
2290 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
2291 . g33,k916,is)
2292 CALL s8zkeg33(1,nel,
2293 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
2294 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
2295 . g33,k917,is)
2296 CALL s8zkeg33(1,nel,
2297 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
2298 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
2299 . g33,k918,is)
2300 CALL s8zkeg33(1,nel,
2301 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
2302 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
2303 . g33,k919,is)
2304 CALL s8zkeg33(1,nel,
2305 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
2306 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
2307 . g33,k920,is)
2308
2309
2310 CALL s8zkeg33(1,nel,
2311 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
2312 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
2313 . g33,k1011,is)
2314 CALL s8zkeg33(1,nel,
2315 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
2316 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
2317 . g33,k1012,is)
2318 CALL s8zkeg33(1,nel,
2319 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
2320 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
2321 . g33,k1013,is)
2322 CALL s8zkeg33(1,nel,
2323 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
2324 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
2325 . g33,k1014,is)
2326 CALL s8zkeg33(1,nel,
2327 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
2328 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
2329 . g33,k1015,is)
2330 CALL s8zkeg33(1,nel,
2331 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
2332 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
2333 . g33,k1016,is)
2334 CALL s8zkeg33(1,nel,
2335 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
2336 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
2337 . g33,k1017,is)
2338 CALL s8zkeg33(1,nel,
2339 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
2340 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
2341 . g33,k1018,is)
2342 CALL s8zkeg33(1,nel,
2343 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
2344 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
2345 . g33,k1019,is)
2346 CALL s8zkeg33(1,nel,
2347 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
2348 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
2349 . g33,k1020,is)
2350
2351
2352
2353 CALL s8zkeg33(1,nel,
2354 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
2355 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
2356 . g33,k1112,is)
2357 CALL s8zkeg33(1,nel,
2358 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
2359 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
2360 . g33,k1113,is)
2361 CALL s8zkeg33(1,nel,
2362 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
2363 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
2364 . g33,k1114,is)
2365 CALL s8zkeg33(1,nel,
2366 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
2367 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
2368 . g33,k1115,is)
2369 CALL s8zkeg33(1,nel,
2370 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
2371 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
2372 . g33,k1116,is)
2373 CALL s8zkeg33(1,nel,
2374 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
2375 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
2376 . g33,k1117,is)
2377 CALL s8zkeg33(1,nel,
2378 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
2379 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
2380 . g33,k1118,is)
2381 CALL s8zkeg33(1,nel,
2382 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
2383 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
2384 . g33,k1119,is)
2385 CALL s8zkeg33(1,nel,
2386 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
2387 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
2388 . g33,k1120,is)
2389
2390
2391
2392 CALL s8zkeg33(1,nel,
2393 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
2394 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
2395 . g33,k1213,is)
2396 CALL s8zkeg33(1,nel,
2397 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
2398 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
2399 . g33,k1214,is)
2400 CALL s8zkeg33(1,nel,
2401 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
2402 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
2403 . g33,k1215,is)
2404 CALL s8zkeg33(1,nel,
2405 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
2406 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
2407 . g33,k1216,is)
2408 CALL s8zkeg33(1,nel,
2409 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
2410 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
2411 . g33,k1217,is)
2412 CALL s8zkeg33(1,nel,
2413 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
2414 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
2415 . g33,k1218,is)
2416 CALL s8zkeg33(1,nel,
2417 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
2418 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
2419 . g33,k1219,is)
2420 CALL s8zkeg33(1,nel,
2421 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
2422 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
2423 . g33,k1220,is)
2424
2425
2426
2427 CALL s8zkeg33(1,nel,
2428 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
2429 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
2430 . g33,k1314,is)
2431 CALL s8zkeg33(1,nel,
2432 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
2433 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
2434 . g33,k1315,is)
2435 CALL s8zkeg33(1,nel,
2436 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
2437 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
2438 . g33,k1316,is)
2439 CALL s8zkeg33(1,nel,
2440 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
2441 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
2442 . g33,k1317,is)
2443 CALL s8zkeg33(1,nel,
2444 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
2445 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
2446 . g33,k1318,is)
2447 CALL s8zkeg33(1,nel,
2448 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
2449 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
2450 . g33,k1319,is)
2451 CALL s8zkeg33(1,nel,
2452 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
2453 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
2454 . g33,k1320,is)
2455
2456
2457 CALL s8zkeg33(1,nel,
2458 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
2459 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
2460 . g33,k1415,is)
2461 CALL s8zkeg33(1,nel,
2462 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
2463 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
2464 . g33,k1416,is)
2465 CALL s8zkeg33(1,nel,
2466 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
2467 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
2468 . g33,k1417,is)
2469 CALL s8zkeg33(1,nel,
2470 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
2471 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
2472 . g33,k1418,is)
2473 CALL s8zkeg33(1,nel,
2474 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
2475 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
2476 . g33,k1419,is)
2477 CALL s8zkeg33(1,nel,
2478 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
2479 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
2480 . g33,k1420,is)
2481
2482
2483 CALL s8zkeg33(1,nel,
2484 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
2485 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
2486 . g33,k1516,is)
2487 CALL s8zkeg33(1,nel,
2488 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
2489 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
2490 . g33,k1517,is)
2491 CALL s8zkeg33(1,nel,
2492 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
2493 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
2494 . g33,k1518,is)
2495 CALL s8zkeg33(1,nel,
2496 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
2497 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
2498 . g33,k1519,is)
2499 CALL s8zkeg33(1,nel,
2500 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
2501 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
2502 . g33,k1520,is)
2503
2504
2505
2506 CALL s8zkeg33(1,nel,
2507 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
2508 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
2509 . g33,k1617,is)
2510 CALL s8zkeg33(1,nel,
2511 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
2512 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
2513 . g33,k1618,is)
2514 CALL s8zkeg33(1,nel,
2515 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
2516 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
2517 . g33,k1619,is)
2518 CALL s8zkeg33(1,nel,
2519 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
2520 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
2521 . g33,k1620,is)
2522
2523
2524 CALL s8zkeg33(1,nel,
2525 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
2526 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
2527 . g33,k1718,is)
2528 CALL s8zkeg33(1,nel,
2529 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
2530 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
2531 . g33,k1719,is)
2532 CALL s8zkeg33(1,nel,
2533 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
2534 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
2535 . g33,k1720,is)
2536
2537
2538 CALL s8zkeg33(1,nel,
2539 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
2540 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
2541 . g33,k1819,is)
2542 CALL s8zkeg33(1,nel,
2543 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
2544 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
2545 . g33,k1820,is)
2546
2547
2548 CALL s8zkeg33(1,nel,
2549 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
2550 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
2551 . g33,k1920,is)
2552
2553
2554C
2555
2556 CALL s8zkebg1(1,nel,
2557 . px(1,1),py(1,1),pz(1,1),px(1,2),py(1,2),pz(1,2),
2558 . dg, k12,is)
2559 CALL s8zkebg1(1,nel,
2560 . px(1,1),py(1,1),pz(1,1),px(1,3),py(1,3),pz(1,3),
2561 . dg, k13,is)
2562 CALL s8zkebg1(1,nel,
2563 . px(1,1),py(1,1),pz(1,1),px(1,4),py(1,4),pz(1,4),
2564 . dg, k14,is)
2565 CALL s8zkebg1(1,nel,
2566 . px(1,1),py(1,1),pz(1,1),px(1,5),py(1,5),pz(1,5),
2567 . dg, k15,is)
2568 CALL s8zkebg1(1,nel,
2569 . px(1,1),py(1,1),pz(1,1),px(1,6),py(1,6),pz(1,6),
2570 . dg, k16,is)
2571 CALL s8zkebg1(1,nel,
2572 . px(1,1),py(1,1),pz(1,1),px(1,7),py(1,7),pz(1,7),
2573 . dg, k17,is)
2574 CALL s8zkebg1(1,nel,
2575 . px(1,1),py(1,1),pz(1,1),px(1,8),py(1,8),pz(1,8),
2576 . dg, k18,is)
2577 CALL s8zkebg1(1,nel,
2578 . px(1,1),py(1,1),pz(1,1),px(1,9),py(1,9),pz(1,9),
2579 . dg, k19,is)
2580 CALL s8zkebg1(1,nel,
2581 . px(1,1),py(1,1),pz(1,1),px(1,10),py(1,10),pz(1,10),
2582 . dg, k110,is)
2583 CALL s8zkebg1(1,nel,
2584 . px(1,1),py(1,1),pz(1,1),px(1,11),py(1,11),pz(1,11),
2585 . dg, k111,is)
2586 CALL s8zkebg1(1,nel,
2587 . px(1,1),py(1,1),pz(1,1),px(1,12),py(1,12),pz(1,12),
2588 . dg, k112,is)
2589 CALL s8zkebg1(1,nel,
2590 . px(1,1),py(1,1),pz(1,1),px(1,13),py(1,13),pz(1,13),
2591 . dg, k113,is)
2592 CALL s8zkebg1(1,nel,
2593 . px(1,1),py(1,1),pz(1,1),px(1,14),py(1,14),pz(1,14),
2594 . dg, k114,is)
2595 CALL s8zkebg1(1,nel,
2596 . px(1,1),py(1,1),pz(1,1),px(1,15),py(1,15),pz(1,15),
2597 . dg, k115,is)
2598 CALL s8zkebg1(1,nel,
2599 . px(1,1),py(1,1),pz(1,1),px(1,16),py(1,16),pz(1,16),
2600 . dg, k116,is)
2601 CALL s8zkebg1(1,nel,
2602 . px(1,1),py(1,1),pz(1,1),px(1,17),py(1,17),pz(1,17),
2603 . dg, k117,is)
2604 CALL s8zkebg1(1,nel,
2605 . px(1,1),py(1,1),pz(1,1),px(1,18),py(1,18),pz(1,18),
2606 . dg, k118,is)
2607 CALL s8zkebg1(1,nel,
2608 . px(1,1),py(1,1),pz(1,1),px(1,19),py(1,19),pz(1,19),
2609 . dg, k119,is)
2610 CALL s8zkebg1(1,nel,
2611 . px(1,1),py(1,1),pz(1,1),px(1,20),py(1,20),pz(1,20),
2612 . dg, k120,is)
2613
2614
2615
2616 CALL s8zkebg1(1,nel,
2617 . px(1,2),py(1,2),pz(1,2),px(1,3),py(1,3),pz(1,3),
2618 . dg, k23,is)
2619 CALL s8zkebg1(1,nel,
2620 . px(1,2),py(1,2),pz(1,2),px(1,4),py(1,4),pz(1,4),
2621 . dg, k24,is)
2622 CALL s8zkebg1(1,nel,
2623 . px(1,2),py(1,2),pz(1,2),px(1,5),py(1,5),pz(1,5),
2624 . dg, k25,is)
2625 CALL s8zkebg1(1,nel,
2626 . px(1,2),py(1,2),pz(1,2),px(1,6),py(1,6),pz(1,6),
2627 . dg, k26,is)
2628 CALL s8zkebg1(1,nel,
2629 . px(1,2),py(1,2),pz(1,2),px(1,7),py(1,7),pz(1,7),
2630 . dg, k27,is)
2631 CALL s8zkebg1(1,nel,
2632 . px(1,2),py(1,2),pz(1,2),px(1,8),py(1,8),pz(1,8),
2633 . dg, k28,is)
2634 CALL s8zkebg1(1,nel,
2635 . px(1,2),py(1,2),pz(1,2),px(1,9),py(1,9),pz(1,9),
2636 . dg, k29,is)
2637 CALL s8zkebg1(1,nel,
2638 . px(1,2),py(1,2),pz(1,2),px(1,10),py(1,10),pz(1,10),
2639 . dg, k210,is)
2640 CALL s8zkebg1(1,nel,
2641 . px(1,2),py(1,2),pz(1,2),px(1,11),py(1,11),pz(1,11),
2642 . dg, k211,is)
2643 CALL s8zkebg1(1,nel,
2644 . px(1,2),py(1,2),pz(1,2),px(1,12),py(1,12),pz(1,12),
2645 . dg, k212,is)
2646 CALL s8zkebg1(1,nel,
2647 . px(1,2),py(1,2),pz(1,2),px(1,13),py(1,13),pz(1,13),
2648 . dg, k213,is)
2649 CALL s8zkebg1(1,nel,
2650 . px(1,2),py(1,2),pz(1,2),px(1,14),py(1,14),pz(1,14),
2651 . dg, k214,is)
2652 CALL s8zkebg1(1,nel,
2653 . px(1,2),py(1,2),pz(1,2),px(1,15),py(1,15),pz(1,15),
2654 . dg, k215,is)
2655 CALL s8zkebg1(1,nel,
2656 . px(1,2),py(1,2),pz(1,2),px(1,16),py(1,16),pz(1,16),
2657 . dg, k216,is)
2658 CALL s8zkebg1(1,nel,
2659 . px(1,2),py(1,2),pz(1,2),px(1,17),py(1,17),pz(1,17),
2660 . dg, k217,is)
2661 CALL s8zkebg1(1,nel,
2662 . px(1,2),py(1,2),pz(1,2),px(1,18),py(1,18),pz(1,18),
2663 . dg, k218,is)
2664 CALL s8zkebg1(1,nel,
2665 . px(1,2),py(1,2),pz(1,2),px(1,19),py(1,19),pz(1,19),
2666 . dg, k219,is)
2667 CALL s8zkebg1(1,nel,
2668 . px(1,2),py(1,2),pz(1,2),px(1,20),py(1,20),pz(1,20),
2669 . dg, k220,is)
2670
2671
2672
2673 CALL s8zkebg1(1,nel,
2674 . px(1,3),py(1,3),pz(1,3),px(1,4),py(1,4),pz(1,4),
2675 . dg, k34,is)
2676 CALL s8zkebg1(1,nel,
2677 . px(1,3),py(1,3),pz(1,3),px(1,5),py(1,5),pz(1,5),
2678 . dg, k35,is)
2679 CALL s8zkebg1(1,nel,
2680 . px(1,3),py(1,3),pz(1,3),px(1,6),py(1,6),pz(1,6),
2681 . dg, k36,is)
2682 CALL s8zkebg1(1,nel,
2683 . px(1,3),py(1,3),pz(1,3),px(1,7),py(1,7),pz(1,7),
2684 . dg, k37,is)
2685 CALL s8zkebg1(1,nel,
2686 . px(1,3),py(1,3),pz(1,3),px(1,8),py(1,8),pz(1,8),
2687 . dg, k38,is)
2688 CALL s8zkebg1(1,nel,
2689 . px(1,3),py(1,3),pz(1,3),px(1,9),py(1,9),pz(1,9),
2690 . dg, k39,is)
2691 CALL s8zkebg1(1,nel,
2692 . px(1,3),py(1,3),pz(1,3),px(1,10),py(1,10),pz(1,10),
2693 . dg, k310,is)
2694 CALL s8zkebg1(1,nel,
2695 . px(1,3),py(1,3),pz(1,3),px(1,11),py(1,11),pz(1,11),
2696 . dg, k311,is)
2697 CALL s8zkebg1(1,nel,
2698 . px(1,3),py(1,3),pz(1,3),px(1,12),py(1,12),pz(1,12),
2699 . dg, k312,is)
2700 CALL s8zkebg1(1,nel,
2701 . px(1,3),py(1,3),pz(1,3),px(1,13),py(1,13),pz(1,13),
2702 . dg, k313,is)
2703 CALL s8zkebg1(1,nel,
2704 . px(1,3),py(1,3),pz(1,3),px(1,14),py(1,14),pz(1,14),
2705 . dg, k314,is)
2706 CALL s8zkebg1(1,nel,
2707 . px(1,3),py(1,3),pz(1,3),px(1,15),py(1,15),pz(1,15),
2708 . dg, k315,is)
2709 CALL s8zkebg1(1,nel,
2710 . px(1,3),py(1,3),pz(1,3),px(1,16),py(1,16),pz(1,16),
2711 . dg, k316,is)
2712 CALL s8zkebg1(1,nel,
2713 . px(1,3),py(1,3),pz(1,3),px(1,17),py(1,17),pz(1,17),
2714 . dg, k317,is)
2715 CALL s8zkebg1(1,nel,
2716 . px(1,3),py(1,3),pz(1,3),px(1,18),py(1,18),pz(1,18),
2717 . dg, k318,is)
2718 CALL s8zkebg1(1,nel,
2719 . px(1,3),py(1,3),pz(1,3),px(1,19),py(1,19),pz(1,19),
2720 . dg, k319,is)
2721 CALL s8zkebg1(1,nel,
2722 . px(1,3),py(1,3),pz(1,3),px(1,20),py(1,20),pz(1,20),
2723 . dg, k320,is)
2724
2725
2726
2727
2728 CALL s8zkebg1(1,nel,
2729 . px(1,4),py(1,4),pz(1,4),px(1,5),py(1,5),pz(1,5),
2730 . dg, k45,is)
2731 CALL s8zkebg1(1,nel,
2732 . px(1,4),py(1,4),pz(1,4),px(1,6),py(1,6),pz(1,6),
2733 . dg, k46,is)
2734 CALL s8zkebg1(1,nel,
2735 . px(1,4),py(1,4),pz(1,4),px(1,7),py(1,7),pz(1,7),
2736 . dg, k47,is)
2737 CALL s8zkebg1(1,nel,
2738 . px(1,4),py(1,4),pz(1,4),px(1,8),py(1,8),pz(1,8),
2739 . dg, k48,is)
2740 CALL s8zkebg1(1,nel,
2741 . px(1,4),py(1,4),pz(1,4),px(1,9),py(1,9),pz(1,9),
2742 . dg, k49,is)
2743 CALL s8zkebg1(1,nel,
2744 . px(1,4),py(1,4),pz(1,4),px(1,10),py(1,10),pz(1,10),
2745 . dg, k410,is)
2746 CALL s8zkebg1(1,nel,
2747 . px(1,4),py(1,4),pz(1,4),px(1,11),py(1,11),pz(1,11),
2748 . dg, k411,is)
2749 CALL s8zkebg1(1,nel,
2750 . px(1,4),py(1,4),pz(1,4),px(1,12),py(1,12),pz(1,12),
2751 . dg, k412,is)
2752 CALL s8zkebg1(1,nel,
2753 . px(1,4),py(1,4),pz(1,4),px(1,13),py(1,13),pz(1,13),
2754 . dg, k413,is)
2755 CALL s8zkebg1(1,nel,
2756 . px(1,4),py(1,4),pz(1,4),px(1,14),py(1,14),pz(1,14),
2757 . dg, k414,is)
2758 CALL s8zkebg1(1,nel,
2759 . px(1,4),py(1,4),pz(1,4),px(1,15),py(1,15),pz(1,15),
2760 . dg, k415,is)
2761 CALL s8zkebg1(1,nel,
2762 . px(1,4),py(1,4),pz(1,4),px(1,16),py(1,16),pz(1,16),
2763 . dg, k416,is)
2764 CALL s8zkebg1(1,nel,
2765 . px(1,4),py(1,4),pz(1,4),px(1,17),py(1,17),pz(1,17),
2766 . dg, k417,is)
2767 CALL s8zkebg1(1,nel,
2768 . px(1,4),py(1,4),pz(1,4),px(1,18),py(1,18),pz(1,18),
2769 . dg, k418,is)
2770 CALL s8zkebg1(1,nel,
2771 . px(1,4),py(1,4),pz(1,4),px(1,19),py(1,19),pz(1,19),
2772 . dg, k419,is)
2773 CALL s8zkebg1(1,nel,
2774 . px(1,4),py(1,4),pz(1,4),px(1,20),py(1,20),pz(1,20),
2775 . dg, k420,is)
2776
2777
2778
2779 CALL s8zkebg1(1,nel,
2780 . px(1,5),py(1,5),pz(1,5),px(1,6),py(1,6),pz(1,6),
2781 . dg, k56,is)
2782 CALL s8zkebg1(1,nel,
2783 . px(1,5),py(1,5),pz(1,5),px(1,7),py(1,7),pz(1,7),
2784 . dg, k57,is)
2785 CALL s8zkebg1(1,nel,
2786 . px(1,5),py(1,5),pz(1,5),px(1,8),py(1,8),pz(1,8),
2787 . dg, k58,is)
2788 CALL s8zkebg1(1,nel,
2789 . px(1,5),py(1,5),pz(1,5),px(1,9),py(1,9),pz(1,9),
2790 . dg, k59,is)
2791 CALL s8zkebg1(1,nel,
2792 . px(1,5),py(1,5),pz(1,5),px(1,10),py(1,10),pz(1,10),
2793 . dg, k510,is)
2794 CALL s8zkebg1(1,nel,
2795 . px(1,5),py(1,5),pz(1,5),px(1,11),py(1,11),pz(1,11),
2796 . dg, k511,is)
2797 CALL s8zkebg1(1,nel,
2798 . px(1,5),py(1,5),pz(1,5),px(1,12),py(1,12),pz(1,12),
2799 . dg, k512,is)
2800 CALL s8zkebg1(1,nel,
2801 . px(1,5),py(1,5),pz(1,5),px(1,13),py(1,13),pz(1,13),
2802 . dg, k513,is)
2803 CALL s8zkebg1(1,nel,
2804 . px(1,5),py(1,5),pz(1,5),px(1,14),py(1,14),pz(1,14),
2805 . dg, k514,is)
2806 CALL s8zkebg1(1,nel,
2807 . px(1,5),py(1,5),pz(1,5),px(1,15),py(1,15),pz(1,15),
2808 . dg, k515,is)
2809 CALL s8zkebg1(1,nel,
2810 . px(1,5),py(1,5),pz(1,5),px(1,16),py(1,16),pz(1,16),
2811 . dg, k516,is)
2812 CALL s8zkebg1(1,nel,
2813 . px(1,5),py(1,5),pz(1,5),px(1,17),py(1,17),pz(1,17),
2814 . dg, k517,is)
2815 CALL s8zkebg1(1,nel,
2816 . px(1,5),py(1,5),pz(1,5),px(1,18),py(1,18),pz(1,18),
2817 . dg, k518,is)
2818 CALL s8zkebg1(1,nel,
2819 . px(1,5),py(1,5),pz(1,5),px(1,19),py(1,19),pz(1,19),
2820 . dg, k519,is)
2821 CALL s8zkebg1(1,nel,
2822 . px(1,5),py(1,5),pz(1,5),px(1,20),py(1,20),pz(1,20),
2823 . dg, k520,is)
2824
2825
2826
2827 CALL s8zkebg1(1,nel,
2828 . px(1,6),py(1,6),pz(1,6),px(1,7),py(1,7),pz(1,7),
2829 . dg, k67,is)
2830 CALL s8zkebg1(1,nel,
2831 . px(1,6),py(1,6),pz(1,6),px(1,8),py(1,8),pz(1,8),
2832 . dg, k68,is)
2833 CALL s8zkebg1(1,nel,
2834 . px(1,6),py(1,6),pz(1,6),px(1,9),py(1,9),pz(1,9),
2835 . dg, k69,is)
2836 CALL s8zkebg1(1,nel,
2837 . px(1,6),py(1,6),pz(1,6),px(1,10),py(1,10),pz(1,10),
2838 . dg, k610,is)
2839 CALL s8zkebg1(1,nel,
2840 . px(1,6),py(1,6),pz(1,6),px(1,11),py(1,11),pz(1,11),
2841 . dg, k611,is)
2842 CALL s8zkebg1(1,nel,
2843 . px(1,6),py(1,6),pz(1,6),px(1,12),py(1,12),pz(1,12),
2844 . dg, k612,is)
2845 CALL s8zkebg1(1,nel,
2846 . px(1,6),py(1,6),pz(1,6),px(1,13),py(1,13),pz(1,13),
2847 . dg, k613,is)
2848 CALL s8zkebg1(1,nel,
2849 . px(1,6),py(1,6),pz(1,6),px(1,14),py(1,14),pz(1,14),
2850 . dg, k614,is)
2851 CALL s8zkebg1(1,nel,
2852 . px(1,6),py(1,6),pz(1,6),px(1,15),py(1,15),pz(1,15),
2853 . dg, k615,is)
2854 CALL s8zkebg1(1,nel,
2855 . px(1,6),py(1,6),pz(1,6),px(1,16),py(1,16),pz(1,16),
2856 . dg, k616,is)
2857 CALL s8zkebg1(1,nel,
2858 . px(1,6),py(1,6),pz(1,6),px(1,17),py(1,17),pz(1,17),
2859 . dg, k617,is)
2860 CALL s8zkebg1(1,nel,
2861 . px(1,6),py(1,6),pz(1,6),px(1,18),py(1,18),pz(1,18),
2862 . dg, k618,is)
2863 CALL s8zkebg1(1,nel,
2864 . px(1,6),py(1,6),pz(1,6),px(1,19),py(1,19),pz(1,19),
2865 . dg, k619,is)
2866 CALL s8zkebg1(1,nel,
2867 . px(1,6),py(1,6),pz(1,6),px(1,20),py(1,20),pz(1,20),
2868 . dg, k620,is)
2869
2870
2871
2872 CALL s8zkebg1(1,nel,
2873 . px(1,7),py(1,7),pz(1,7),px(1,8),py(1,8),pz(1,8),
2874 . dg, k78,is)
2875 CALL s8zkebg1(1,nel,
2876 . px(1,7),py(1,7),pz(1,7),px(1,9),py(1,9),pz(1,9),
2877 . dg, k79,is)
2878 CALL s8zkebg1(1,nel,
2879 . px(1,7),py(1,7),pz(1,7),px(1,10),py(1,10),pz(1,10),
2880 . dg, k710,is)
2881 CALL s8zkebg1(1,nel,
2882 . px(1,7),py(1,7),pz(1,7),px(1,11),py(1,11),pz(1,11),
2883 . dg, k711,is)
2884 CALL s8zkebg1(1,nel,
2885 . px(1,7),py(1,7),pz(1,7),px(1,12),py(1,12),pz(1,12),
2886 . dg, k712,is)
2887 CALL s8zkebg1(1,nel,
2888 . px(1,7),py(1,7),pz(1,7),px(1,13),py(1,13),pz(1,13),
2889 . dg, k713,is)
2890 CALL s8zkebg1(1,nel,
2891 . px(1,7),py(1,7),pz(1,7),px(1,14),py(1,14),pz(1,14),
2892 . dg, k714,is)
2893 CALL s8zkebg1(1,nel,
2894 . px(1,7),py(1,7),pz(1,7),px(1,15),py(1,15),pz(1,15),
2895 . dg, k715,is)
2896 CALL s8zkebg1(1,nel,
2897 . px(1,7),py(1,7),pz(1,7),px(1,16),py(1,16),pz(1,16),
2898 . dg, k716,is)
2899 CALL s8zkebg1(1,nel,
2900 . px(1,7),py(1,7),pz(1,7),px(1,17),py(1,17),pz(1,17),
2901 . dg, k717,is)
2902 CALL s8zkebg1(1,nel,
2903 . px(1,7),py(1,7),pz(1,7),px(1,18),py(1,18),pz(1,18),
2904 . dg, k718,is)
2905 CALL s8zkebg1(1,nel,
2906 . px(1,7),py(1,7),pz(1,7),px(1,19),py(1,19),pz(1,19),
2907 . dg, k719,is)
2908 CALL s8zkebg1(1,nel,
2909 . px(1,7),py(1,7),pz(1,7),px(1,20),py(1,20),pz(1,20),
2910 . dg, k720,is)
2911
2912
2913
2914 CALL s8zkebg1(1,nel,
2915 . px(1,8),py(1,8),pz(1,8),px(1,9),py(1,9),pz(1,9),
2916 . dg, k89,is)
2917 CALL s8zkebg1(1,nel,
2918 . px(1,8),py(1,8),pz(1,8),px(1,10),py(1,10),pz(1,10),
2919 . dg, k810,is)
2920 CALL s8zkebg1(1,nel,
2921 . px(1,8),py(1,8),pz(1,8),px(1,11),py(1,11),pz(1,11),
2922 . dg, k811,is)
2923 CALL s8zkebg1(1,nel,
2924 . px(1,8),py(1,8),pz(1,8),px(1,12),py(1,12),pz(1,12),
2925 . dg, k812,is)
2926 CALL s8zkebg1(1,nel,
2927 . px(1,8),py(1,8),pz(1,8),px(1,13),py(1,13),pz(1,13),
2928 . dg, k813,is)
2929 CALL s8zkebg1(1,nel,
2930 . px(1,8),py(1,8),pz(1,8),px(1,14),py(1,14),pz(1,14),
2931 . dg, k814,is)
2932 CALL s8zkebg1(1,nel,
2933 . px(1,8),py(1,8),pz(1,8),px(1,15),py(1,15),pz(1,15),
2934 . dg, k815,is)
2935 CALL s8zkebg1(1,nel,
2936 . px(1,8),py(1,8),pz(1,8),px(1,16),py(1,16),pz(1,16),
2937 . dg, k816,is)
2938 CALL s8zkebg1(1,nel,
2939 . px(1,8),py(1,8),pz(1,8),px(1,17),py(1,17),pz(1,17),
2940 . dg, k817,is)
2941 CALL s8zkebg1(1,nel,
2942 . px(1,8),py(1,8),pz(1,8),px(1,18),py(1,18),pz(1,18),
2943 . dg, k818,is)
2944 CALL s8zkebg1(1,nel,
2945 . px(1,8),py(1,8),pz(1,8),px(1,19),py(1,19),pz(1,19),
2946 . dg, k819,is)
2947 CALL s8zkebg1(1,nel,
2948 . px(1,8),py(1,8),pz(1,8),px(1,20),py(1,20),pz(1,20),
2949 . dg, k820,is)
2950
2951
2952
2953 CALL s8zkebg1(1,nel,
2954 . px(1,9),py(1,9),pz(1,9),px(1,10),py(1,10),pz(1,10),
2955 . dg, k910,is)
2956 CALL s8zkebg1(1,nel,
2957 . px(1,9),py(1,9),pz(1,9),px(1,11),py(1,11),pz(1,11),
2958 . dg, k911,is)
2959 CALL s8zkebg1(1,nel,
2960 . px(1,9),py(1,9),pz(1,9),px(1,12),py(1,12),pz(1,12),
2961 . dg, k912,is)
2962 CALL s8zkebg1(1,nel,
2963 . px(1,9),py(1,9),pz(1,9),px(1,13),py(1,13),pz(1,13),
2964 . dg, k913,is)
2965 CALL s8zkebg1(1,nel,
2966 . px(1,9),py(1,9),pz(1,9),px(1,14),py(1,14),pz(1,14),
2967 . dg, k914,is)
2968 CALL s8zkebg1(1,nel,
2969 . px(1,9),py(1,9),pz(1,9),px(1,15),py(1,15),pz(1,15),
2970 . dg, k915,is)
2971 CALL s8zkebg1(1,nel,
2972 . px(1,9),py(1,9),pz(1,9),px(1,16),py(1,16),pz(1,16),
2973 . dg, k916,is)
2974 CALL s8zkebg1(1,nel,
2975 . px(1,9),py(1,9),pz(1,9),px(1,17),py(1,17),pz(1,17),
2976 . dg, k917,is)
2977 CALL s8zkebg1(1,nel,
2978 . px(1,9),py(1,9),pz(1,9),px(1,18),py(1,18),pz(1,18),
2979 . dg, k918,is)
2980 CALL s8zkebg1(1,nel,
2981 . px(1,9),py(1,9),pz(1,9),px(1,19),py(1,19),pz(1,19),
2982 . dg, k919,is)
2983 CALL s8zkebg1(1,nel,
2984 . px(1,9),py(1,9),pz(1,9),px(1,20),py(1,20),pz(1,20),
2985 . dg, k920,is)
2986
2987
2988 CALL s8zkebg1(1,nel,
2989 . px(1,10),py(1,10),pz(1,10),px(1,11),py(1,11),pz(1,11),
2990 . dg, k1011,is)
2991 CALL s8zkebg1(1,nel,
2992 . px(1,10),py(1,10),pz(1,10),px(1,12),py(1,12),pz(1,12),
2993 . dg, k1012,is)
2994 CALL s8zkebg1(1,nel,
2995 . px(1,10),py(1,10),pz(1,10),px(1,13),py(1,13),pz(1,13),
2996 . dg, k1013,is)
2997 CALL s8zkebg1(1,nel,
2998 . px(1,10),py(1,10),pz(1,10),px(1,14),py(1,14),pz(1,14),
2999 . dg, k1014,is)
3000 CALL s8zkebg1(1,nel,
3001 . px(1,10),py(1,10),pz(1,10),px(1,15),py(1,15),pz(1,15),
3002 . dg, k1015,is)
3003 CALL s8zkebg1(1,nel,
3004 . px(1,10),py(1,10),pz(1,10),px(1,16),py(1,16),pz(1,16),
3005 . dg, k1016,is)
3006 CALL s8zkebg1(1,nel,
3007 . px(1,10),py(1,10),pz(1,10),px(1,17),py(1,17),pz(1,17),
3008 . dg, k1017,is)
3009 CALL s8zkebg1(1,nel,
3010 . px(1,10),py(1,10),pz(1,10),px(1,18),py(1,18),pz(1,18),
3011 . dg, k1018,is)
3012 CALL s8zkebg1(1,nel,
3013 . px(1,10),py(1,10),pz(1,10),px(1,19),py(1,19),pz(1,19),
3014 . dg, k1019,is)
3015 CALL s8zkebg1(1,nel,
3016 . px(1,10),py(1,10),pz(1,10),px(1,20),py(1,20),pz(1,20),
3017 . dg, k1020,is)
3018
3019
3020 CALL s8zkebg1(1,nel,
3021 . px(1,11),py(1,11),pz(1,11),px(1,12),py(1,12),pz(1,12),
3022 . dg, k1112,is)
3023 CALL s8zkebg1(1,nel,
3024 . px(1,11),py(1,11),pz(1,11),px(1,13),py(1,13),pz(1,13),
3025 . dg, k1113,is)
3026 CALL s8zkebg1(1,nel,
3027 . px(1,11),py(1,11),pz(1,11),px(1,14),py(1,14),pz(1,14),
3028 . dg, k1114,is)
3029 CALL s8zkebg1(1,nel,
3030 . px(1,11),py(1,11),pz(1,11),px(1,15),py(1,15),pz(1,15),
3031 . dg, k1115,is)
3032 CALL s8zkebg1(1,nel,
3033 . px(1,11),py(1,11),pz(1,11),px(1,16),py(1,16),pz(1,16),
3034 . dg, k1116,is)
3035 CALL s8zkebg1(1,nel,
3036 . px(1,11),py(1,11),pz(1,11),px(1,17),py(1,17),pz(1,17),
3037 . dg, k1117,is)
3038 CALL s8zkebg1(1,nel,
3039 . px(1,11),py(1,11),pz(1,11),px(1,18),py(1,18),pz(1,18),
3040 . dg, k1118,is)
3041 CALL s8zkebg1(1,nel,
3042 . px(1,11),py(1,11),pz(1,11),px(1,19),py(1,19),pz(1,19),
3043 . dg, k1119,is)
3044 CALL s8zkebg1(1,nel,
3045 . px(1,11),py(1,11),pz(1,11),px(1,20),py(1,20),pz(1,20),
3046 . dg, k1120,is)
3047
3048
3049 CALL s8zkebg1(1,nel,
3050 . px(1,12),py(1,12),pz(1,12),px(1,13),py(1,13),pz(1,13),
3051 . dg, k1213,is)
3052 CALL s8zkebg1(1,nel,
3053 . px(1,12),py(1,12),pz(1,12),px(1,14),py(1,14),pz(1,14),
3054 . dg, k1214,is)
3055 CALL s8zkebg1(1,nel,
3056 . px(1,12),py(1,12),pz(1,12),px(1,15),py(1,15),pz(1,15),
3057 . dg, k1215,is)
3058 CALL s8zkebg1(1,nel,
3059 . px(1,12),py(1,12),pz(1,12),px(1,16),py(1,16),pz(1,16),
3060 . dg, k1216,is)
3061 CALL s8zkebg1(1,nel,
3062 . px(1,12),py(1,12),pz(1,12),px(1,17),py(1,17),pz(1,17),
3063 . dg, k1217,is)
3064 CALL s8zkebg1(1,nel,
3065 . px(1,12),py(1,12),pz(1,12),px(1,18),py(1,18),pz(1,18),
3066 . dg, k1218,is)
3067 CALL s8zkebg1(1,nel,
3068 . px(1,12),py(1,12),pz(1,12),px(1,19),py(1,19),pz(1,19),
3069 . dg, k1219,is)
3070 CALL s8zkebg1(1,nel,
3071 . px(1,12),py(1,12),pz(1,12),px(1,20),py(1,20),pz(1,20),
3072 . dg, k1220,is)
3073
3074
3075 CALL s8zkebg1(1,nel,
3076 . px(1,13),py(1,13),pz(1,13),px(1,14),py(1,14),pz(1,14),
3077 . dg, k1314,is)
3078 CALL s8zkebg1(1,nel,
3079 . px(1,13),py(1,13),pz(1,13),px(1,15),py(1,15),pz(1,15),
3080 . dg, k1315,is)
3081 CALL s8zkebg1(1,nel,
3082 . px(1,13),py(1,13),pz(1,13),px(1,16),py(1,16),pz(1,16),
3083 . dg, k1316,is)
3084 CALL s8zkebg1(1,nel,
3085 . px(1,13),py(1,13),pz(1,13),px(1,17),py(1,17),pz(1,17),
3086 . dg, k1317,is)
3087 CALL s8zkebg1(1,nel,
3088 . px(1,13),py(1,13),pz(1,13),px(1,18),py(1,18),pz(1,18),
3089 . dg, k1318,is)
3090 CALL s8zkebg1(1,nel,
3091 . px(1,13),py(1,13),pz(1,13),px(1,19),py(1,19),pz(1,19),
3092 . dg, k1319,is)
3093 CALL s8zkebg1(1,nel,
3094 . px(1,13),py(1,13),pz(1,13),px(1,20),py(1,20),pz(1,20),
3095 . dg, k1320,is)
3096
3097
3098
3099 CALL s8zkebg1(1,nel,
3100 . px(1,14),py(1,14),pz(1,14),px(1,15),py(1,15),pz(1,15),
3101 . dg, k1415,is)
3102 CALL s8zkebg1(1,nel,
3103 . px(1,14),py(1,14),pz(1,14),px(1,16),py(1,16),pz(1,16),
3104 . dg, k1416,is)
3105 CALL s8zkebg1(1,nel,
3106 . px(1,14),py(1,14),pz(1,14),px(1,17),py(1,17),pz(1,17),
3107 . dg, k1417,is)
3108 CALL s8zkebg1(1,nel,
3109 . px(1,14),py(1,14),pz(1,14),px(1,18),py(1,18),pz(1,18),
3110 . dg, k1418,is)
3111 CALL s8zkebg1(1,nel,
3112 . px(1,14),py(1,14),pz(1,14),px(1,19),py(1,19),pz(1,19),
3113 . dg, k1419,is)
3114 CALL s8zkebg1(1,nel,
3115 . px(1,14),py(1,14),pz(1,14),px(1,20),py(1,20),pz(1,20),
3116 . dg, k1420,is)
3117
3118
3119 CALL s8zkebg1(1,nel,
3120 . px(1,15),py(1,15),pz(1,15),px(1,16),py(1,16),pz(1,16),
3121 . dg, k1516,is)
3122 CALL s8zkebg1(1,nel,
3123 . px(1,15),py(1,15),pz(1,15),px(1,17),py(1,17),pz(1,17),
3124 . dg, k1517,is)
3125 CALL s8zkebg1(1,nel,
3126 . px(1,15),py(1,15),pz(1,15),px(1,18),py(1,18),pz(1,18),
3127 . dg, k1518,is)
3128 CALL s8zkebg1(1,nel,
3129 . px(1,15),py(1,15),pz(1,15),px(1,19),py(1,19),pz(1,19),
3130 . dg, k1519,is)
3131 CALL s8zkebg1(1,nel,
3132 . px(1,15),py(1,15),pz(1,15),px(1,20),py(1,20),pz(1,20),
3133 . dg, k1520,is)
3134
3135
3136 CALL s8zkebg1(1,nel,
3137 . px(1,16),py(1,16),pz(1,16),px(1,17),py(1,17),pz(1,17),
3138 . dg, k1617,is)
3139 CALL s8zkebg1(1,nel,
3140 . px(1,16),py(1,16),pz(1,16),px(1,18),py(1,18),pz(1,18),
3141 . dg, k1618,is)
3142 CALL s8zkebg1(1,nel,
3143 . px(1,16),py(1,16),pz(1,16),px(1,19),py(1,19),pz(1,19),
3144 . dg, k1619,is)
3145 CALL s8zkebg1(1,nel,
3146 . px(1,16),py(1,16),pz(1,16),px(1,20),py(1,20),pz(1,20),
3147 . dg, k1620,is)
3148
3149
3150 CALL s8zkebg1(1,nel,
3151 . px(1,17),py(1,17),pz(1,17),px(1,18),py(1,18),pz(1,18),
3152 . dg, k1718,is)
3153 CALL s8zkebg1(1,nel,
3154 . px(1,17),py(1,17),pz(1,17),px(1,19),py(1,19),pz(1,19),
3155 . dg, k1719,is)
3156 CALL s8zkebg1(1,nel,
3157 . px(1,17),py(1,17),pz(1,17),px(1,20),py(1,20),pz(1,20),
3158 . dg, k1720,is)
3159
3160
3161 CALL s8zkebg1(1,nel,
3162 . px(1,18),py(1,18),pz(1,18),px(1,19),py(1,19),pz(1,19),
3163 . dg, k1819,is)
3164 CALL s8zkebg1(1,nel,
3165 . px(1,18),py(1,18),pz(1,18),px(1,20),py(1,20),pz(1,20),
3166 . dg, k1820,is)
3167
3168
3169 CALL s8zkebg1(1,nel,
3170 . px(1,19),py(1,19),pz(1,19),px(1,20),py(1,20),pz(1,20),
3171 . dg, k1920,is)
3172
3173
3174C
3175 is =1
3176 CALL s8zkebg1(1,nel,
3177 . px(1,1),py(1,1),pz(1,1),px(1,1),py(1,1),pz(1,1),
3178 . dg, k11,is)
3179 CALL s8zkebg1(1,nel,
3180 . px(1,2),py(1,2),pz(1,2),px(1,2),py(1,2),pz(1,2),
3181 . dg, k22,is)
3182 CALL s8zkebg1(1,nel,
3183 . px(1,3),py(1,3),pz(1,3),px(1,3),py(1,3),pz(1,3),
3184 . dg, k33,is)
3185 CALL s8zkebg1(1,nel,
3186 . px(1,4),py(1,4),pz(1,4),px(1,4),py(1,4),pz(1,4),
3187 . dg, k44,is)
3188 CALL s8zkebg1(1,nel,
3189 . px(1,5),py(1,5),pz(1,5),px(1,5),py(1,5),pz(1,5),
3190 . dg, k55,is)
3191 CALL s8zkebg1(1,nel,
3192 . px(1,6),py(1,6),pz(1,6),px(1,6),py(1,6),pz(1,6),
3193 . dg, k66,is)
3194 CALL s8zkebg1(1,nel,
3195 . px(1,7),py(1,7),pz(1,7),px(1,7),py(1,7),pz(1,7),
3196 . dg, k77,is)
3197 CALL s8zkebg1(1,nel,
3198 . px(1,8),py(1,8),pz(1,8),px(1,8),py(1,8),pz(1,8),
3199 . dg, k88,is)
3200 CALL s8zkebg1(1,nel,
3201 . px(1,9),py(1,9),pz(1,9),px(1,9),py(1,9),pz(1,9),
3202 . dg, k99,is)
3203 CALL s8zkebg1(1,nel,
3204 . px(1,10),py(1,10),pz(1,10),px(1,10),py(1,10),pz(1,10),
3205 . dg, k1010,is)
3206 CALL s8zkebg1(1,nel,
3207 . px(1,11),py(1,11),pz(1,11),px(1,11),py(1,11),pz(1,11),
3208 . dg, k1111,is)
3209 CALL s8zkebg1(1,nel,
3210 . px(1,12),py(1,12),pz(1,12),px(1,12),py(1,12),pz(1,12),
3211 . dg, k1212,is)
3212 CALL s8zkebg1(1,nel,
3213 . px(1,13),py(1,13),pz(1,13),px(1,13),py(1,13),pz(1,13),
3214 . dg, k1313,is)
3215 CALL s8zkebg1(1,nel,
3216 . px(1,14),py(1,14),pz(1,14),px(1,14),py(1,14),pz(1,14),
3217 . dg, k1414,is)
3218 CALL s8zkebg1(1,nel,
3219 . px(1,15),py(1,15),pz(1,15),px(1,15),py(1,15),pz(1,15),
3220 . dg, k1515,is)
3221 CALL s8zkebg1(1,nel,
3222 . px(1,16),py(1,16),pz(1,16),px(1,16),py(1,16),pz(1,16),
3223 . dg, k1616,is)
3224 CALL s8zkebg1(1,nel,
3225 . px(1,17),py(1,17),pz(1,17),px(1,17),py(1,17),pz(1,17),
3226 . dg, k1717,is)
3227 CALL s8zkebg1(1,nel,
3228 . px(1,18),py(1,18),pz(1,18),px(1,18),py(1,18),pz(1,18),
3229 . dg, k1818,is)
3230 CALL s8zkebg1(1,nel,
3231 . px(1,19),py(1,19),pz(1,19),px(1,19),py(1,19),pz(1,19),
3232 . dg, k1919,is)
3233 CALL s8zkebg1(1,nel,
3234 . px(1,20),py(1,20),pz(1,20),px(1,20),py(1,20),pz(1,20),
3235 . dg, k2020,is)
3236
3237
3238C
3239 CALL s8zkeg33(1,nel,
3240 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
3241 . py(1,1),px(1,1),pz(1,1),px(1,1),pz(1,1),py(1,1),
3242 . g33,k11,is)
3243 CALL s8zkeg33(1,nel,
3244 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
3245 . py(1,2),px(1,2),pz(1,2),px(1,2),pz(1,2),py(1,2),
3246 . g33,k22,is)
3247 CALL s8zkeg33(1,nel,
3248 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
3249 . py(1,3),px(1,3),pz(1,3),px(1,3),pz(1,3),py(1,3),
3250 . g33,k33,is)
3251 CALL s8zkeg33(1,nel,
3252 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
3253 . py(1,4),px(1,4),pz(1,4),px(1,4),pz(1,4),py(1,4),
3254 . g33,k44,is)
3255 CALL s8zkeg33(1,nel,
3256 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
3257 . py(1,5),px(1,5),pz(1,5),px(1,5),pz(1,5),py(1,5),
3258 . g33,k55,is)
3259 CALL s8zkeg33(1,nel,
3260 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
3261 . py(1,6),px(1,6),pz(1,6),px(1,6),pz(1,6),py(1,6),
3262 . g33,k66,is)
3263 CALL s8zkeg33(1,nel,
3264 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
3265 . py(1,7),px(1,7),pz(1,7),px(1,7),pz(1,7),py(1,7),
3266 . g33,k77,is)
3267 CALL s8zkeg33(1,nel,
3268 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
3269 . py(1,8),px(1,8),pz(1,8),px(1,8),pz(1,8),py(1,8),
3270 . g33,k88,is)
3271 CALL s8zkeg33(1,nel,
3272 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
3273 . py(1,9),px(1,9),pz(1,9),px(1,9),pz(1,9),py(1,9),
3274 . g33,k99,is)
3275 CALL s8zkeg33(1,nel,
3276 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
3277 . py(1,10),px(1,10),pz(1,10),px(1,10),pz(1,10),py(1,10),
3278 . g33,k1010,is)
3279 CALL s8zkeg33(1,nel,
3280 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
3281 . py(1,11),px(1,11),pz(1,11),px(1,11),pz(1,11),py(1,11),
3282 . g33,k1111,is)
3283 CALL s8zkeg33(1,nel,
3284 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
3285 . py(1,12),px(1,12),pz(1,12),px(1,12),pz(1,12),py(1,12),
3286 . g33,k1212,is)
3287 CALL s8zkeg33(1,nel,
3288 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
3289 . py(1,13),px(1,13),pz(1,13),px(1,13),pz(1,13),py(1,13),
3290 . g33,k1313,is)
3291 CALL s8zkeg33(1,nel,
3292 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
3293 . py(1,14),px(1,14),pz(1,14),px(1,14),pz(1,14),py(1,14),
3294 . g33,k1414,is)
3295 CALL s8zkeg33(1,nel,
3296 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
3297 . py(1,15),px(1,15),pz(1,15),px(1,15),pz(1,15),py(1,15),
3298 . g33,k1515,is)
3299 CALL s8zkeg33(1,nel,
3300 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
3301 . py(1,16),px(1,16),pz(1,16),px(1,16),pz(1,16),py(1,16),
3302 . g33,k1616,is)
3303 CALL s8zkeg33(1,nel,
3304 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
3305 . py(1,17),px(1,17),pz(1,17),px(1,17),pz(1,17),py(1,17),
3306 . g33,k1717,is)
3307 CALL s8zkeg33(1,nel,
3308 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
3309 . py(1,18),px(1,18),pz(1,18),px(1,18),pz(1,18),py(1,18),
3310 . g33,k1818,is)
3311 CALL s8zkeg33(1,nel,
3312 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
3313 . py(1,19),px(1,19),pz(1,19),px(1,19),pz(1,19),py(1,19),
3314 . g33,k1919,is)
3315 CALL s8zkeg33(1,nel,
3316 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
3317 . py(1,20),px(1,20),pz(1,20),px(1,20),pz(1,20),py(1,20),
3318 . g33,k2020,is)
3319
3320
3321 ENDIF
3322
3323C
3324 RETURN
3325 END
s20cumg3
subroutine s20cumg3(px, py, pz, dd, gg, dg, g33, iksup, 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 s20cumg3.F:89
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