55
56
57
58 USE elbufdef_mod
59 use element_mod , only : nixs
60
61
62
63#include "implicit_f.inc"
64
65
66
67#include "mvsiz_p.inc"
68
69
70
71#include "com04_c.inc"
72#include "param_c.inc"
73
74
75
76 INTEGER, INTENT(IN) :: IREP
77 INTEGER, INTENT(IN) :: NFT
78 INTEGER, INTENT(IN) :: MTN
79 INTEGER, INTENT(IN) :: ISMSTR
80 INTEGER, INTENT(IN) :: JHBE
81 INTEGER, INTENT(IN) :: JCVT
82 INTEGER, INTENT(IN) :: IGTYP
83 INTEGER, INTENT(IN) :: ISORTH
84 INTEGER NEL ,ICP, ICSIG,IKGEO,MPT
85 INTEGER NDOF(*) ,IADK(*) ,JDIK(*),
86 . IXS(NIXS,*),IPM(NPROPMI,*),IGEO(NPROPGI,*),ETAG(*),IDDL(*)
87
89 . pm(npropm,*), geo(npropg,*), x(*),
90 . k11(9,mvsiz) , k12(9,mvsiz) ,k13(9,mvsiz) ,k14(9,mvsiz) ,
91 . k15(9,mvsiz) , k16(9,mvsiz) ,k17(9,mvsiz) ,k18(9,mvsiz) ,
92 . k22(9,mvsiz) ,k23(9,mvsiz) ,k24(9,mvsiz) ,k25(9,mvsiz) ,
93 . k26(9,mvsiz) ,k27(9,mvsiz) ,k28(9,mvsiz) ,k33(9,mvsiz) ,
94 . k34(9,mvsiz) ,k35(9,mvsiz) ,k36(9,mvsiz) ,k37(9,mvsiz) ,
95 . k38(9,mvsiz) ,k44(9,mvsiz) ,k45(9,mvsiz) ,k46(9,mvsiz) ,
96 . k47(9,mvsiz) ,k48(9,mvsiz) ,k55(9,mvsiz) ,k56(9,mvsiz) ,
97 . k57(9,mvsiz) ,k58(9,mvsiz) ,k66(9,mvsiz) ,k67(9,mvsiz) ,
98 . k68(9,mvsiz) ,k77(9,mvsiz) ,k78(9,mvsiz) ,k88(9,mvsiz) ,
99 . offg(mvsiz) ,bufmat(*) ,k_diag(*) ,k_lt(*)
100 TYPE (ELBUF_STRUCT_), TARGET :: ELBUF_STR
101
102
103
104 INTEGER NF1, I, IKORTH, IADBUF, IPREDU,
105 . IR, IS, IT,IL,IP,ICPG,NPTR,NPTS,NPTT,NLAY,NPPT,
106 . IAD0,PID
107 INTEGER MXT(MVSIZ),NGL(MVSIZ),NGEO(MVSIZ),IDEG(MVSIZ)
109 . voln(mvsiz), deltax(mvsiz),
110 . ajc1(mvsiz) , ajc2(mvsiz) , ajc3(mvsiz) ,
111 . ajc4(mvsiz) , ajc5(mvsiz) , ajc6(mvsiz) ,
112 . ajc7(mvsiz) , ajc8(mvsiz) , ajc9(mvsiz) ,
113 . aj1(mvsiz) , aj2(mvsiz) , aj3(mvsiz) ,
114 . aj4(mvsiz) , aj5(mvsiz) , aj6(mvsiz)
115
116 INTEGER NC1(MVSIZ), NC2(MVSIZ), NC3(MVSIZ), NC4(MVSIZ),
117 . NC5(MVSIZ), NC6(MVSIZ), NC7(MVSIZ), NC8(MVSIZ)
119 . off(mvsiz) ,
120 . x1(mvsiz), x2(mvsiz), x3(mvsiz), x4(mvsiz),
121 . x5(mvsiz), x6(mvsiz), x7(mvsiz), x8(mvsiz),
122 . y1(mvsiz), y2(mvsiz), y3(mvsiz), y4(mvsiz),
123 . y5(mvsiz), y6(mvsiz), y7(mvsiz), y8(mvsiz),
124 . z1(mvsiz), z2(mvsiz), z3(mvsiz), z4(mvsiz),
125 . z5(mvsiz), z6(mvsiz), z7(mvsiz), z8(mvsiz),
126 . px1(mvsiz),px2(mvsiz),px3(mvsiz),px4(mvsiz),
127 . px5(mvsiz),px6(mvsiz),px7(mvsiz),px8(mvsiz),
128 . py1(mvsiz),py2(mvsiz),py3(mvsiz),py4(mvsiz),
129 . py5(mvsiz),py6(mvsiz),py7(mvsiz),py8(mvsiz),
130 . pz1(mvsiz),pz2(mvsiz),pz3(mvsiz),pz4(mvsiz),
131 . pz5(mvsiz),pz6(mvsiz),pz7(mvsiz),pz8(mvsiz),
132 . pxy1(mvsiz),pxy2(mvsiz),pxy3(mvsiz),pxy4(mvsiz),
133 . pxy5(mvsiz),pxy6(mvsiz),pxy7(mvsiz),pxy8(mvsiz),
134 . pyx1(mvsiz),pyx2(mvsiz),pyx3(mvsiz),pyx4(mvsiz),
135 . pyx5(mvsiz),pyx6(mvsiz),pyx7(mvsiz),pyx8(mvsiz),
136 . pxz1(mvsiz),pxz2(mvsiz),pxz3(mvsiz),pxz4(mvsiz),
137 . pxz5(mvsiz),pxz6(mvsiz),pxz7(mvsiz),pxz8(mvsiz),
138 . pzx1(mvsiz),pzx2(mvsiz),pzx3(mvsiz),pzx4(mvsiz),
139 . pzx5(mvsiz),pzx6(mvsiz),pzx7(mvsiz),pzx8(mvsiz),
140 . pyz1(mvsiz),pyz2(mvsiz),pyz3(mvsiz),pyz4(mvsiz),
141 . pyz5(mvsiz),pyz6(mvsiz),pyz7(mvsiz),pyz8(mvsiz),
142 . pzy1(mvsiz),pzy2(mvsiz),pzy3(mvsiz),pzy4(mvsiz),
143 . pzy5(mvsiz),pzy6(mvsiz),pzy7(mvsiz),pzy8(mvsiz),
144 . pxc1(mvsiz),pxc2(mvsiz),pxc3(mvsiz),pxc4(mvsiz),
145 . pyc1(mvsiz),pyc2(mvsiz),pyc3(mvsiz),pyc4(mvsiz),
146 . pzc1(mvsiz),pzc2(mvsiz),pzc3(mvsiz),pzc4(mvsiz),
147 . bxy1(mvsiz),bxy2(mvsiz),bxy3(mvsiz),bxy4(mvsiz),
148 . bxy5(mvsiz),bxy6(mvsiz),bxy7(mvsiz),bxy8(mvsiz),
149 . byx1(mvsiz),byx2(mvsiz),byx3(mvsiz),byx4(mvsiz),
150 . byx5(mvsiz),byx6(mvsiz),byx7(mvsiz),byx8(mvsiz),
151 . bxz1(mvsiz),bxz2(mvsiz),bxz3(mvsiz),bxz4(mvsiz),
152 . bxz5(mvsiz),bxz6(mvsiz),bxz7(mvsiz),bxz8(mvsiz),
153 . bzx1(mvsiz),bzx2(mvsiz),bzx3(mvsiz),bzx4(mvsiz),
154 . bzx5(mvsiz),bzx6(mvsiz),bzx7(mvsiz),bzx8(mvsiz),
155 . byz1(mvsiz),byz2(mvsiz),byz3(mvsiz),byz4(mvsiz),
156 . byz5(mvsiz),byz6(mvsiz),byz7(mvsiz),byz8(mvsiz),
157 . bzy1(mvsiz),bzy2(mvsiz),bzy3(mvsiz),bzy4(mvsiz),
158 . bzy5(mvsiz),bzy6(mvsiz),bzy7(mvsiz),bzy8(mvsiz),
159 . px1h1(mvsiz),px2h1(mvsiz),px3h1(mvsiz),px4h1(mvsiz),
160 . px1h2(mvsiz),px2h2(mvsiz),px3h2(mvsiz),px4h2(mvsiz),
161 . px1h3(mvsiz),px2h3(mvsiz),px3h3(mvsiz),px4h3(mvsiz),
162 . px1h4(mvsiz),px2h4(mvsiz),px3h4(mvsiz),px4h4(mvsiz),
163 . hx(mvsiz,4), hy(mvsiz,4), hz(mvsiz,4)
165 . r11(mvsiz),r12(mvsiz),r13(mvsiz),
166 . r21(mvsiz),r22(mvsiz),r23(mvsiz),
167 . r31(mvsiz),r32(mvsiz),r33(mvsiz)
169 . bid,wi,smax(mvsiz),volg(mvsiz),nu(mvsiz),
170 . nu1(mvsiz),gama(mvsiz,6)
171 TYPE(G_BUFEL_) ,POINTER :: GBUF
172 TYPE(L_BUFEL_) ,POINTER :: LBUF
174 . w_gauss(9,9),a_gauss(9,9)
175 DATA w_gauss /
176 1 2. ,0. ,0. ,
177 1 0. ,0. ,0. ,
178 1 0. ,0. ,0. ,
179 2 1. ,1. ,0. ,
180 2 0. ,0. ,0. ,
181 2 0. ,0. ,0. ,
182 3 0.555555555555556,0.888888888888889,0.555555555555556,
183 3 0. ,0. ,0. ,
184 3 0. ,0. ,0. ,
185 4 0.347854845137454,0.652145154862546,0.652145154862546,
186 4 0.347854845137454,0. ,0. ,
187 4 0. ,0. ,0. ,
188 5 0.236926885056189,0.478628670499366,0.568888888888889,
189 5 0.478628670499366,0.236926885056189,0. ,
190 5 0. ,0. ,0. ,
191 6 0.171324492379170,0.360761573048139,0.467913934572691,
192 6 0.467913934572691,0.360761573048139,0.171324492379170,
193 6 0. ,0. ,0. ,
194 7 0.129484966168870,0.279705391489277,0.381830050505119,
195 7 0.417959183673469,0.381830050505119,0.279705391489277,
196 7 0.129484966168870,0. ,0. ,
197 8 0.101228536290376,0.222381034453374,0.313706645877887,
198 8 0.362683783378362,0.362683783378362,0.313706645877887,
199 8 0.222381034453374,0.101228536290376,0. ,
200 9 0.081274388361574,0.180648160694857,0.260610696402935,
201 9 0.312347077040003,0.330239355001260,0.312347077040003,
202 9 0.260610696402935,0.180648160694857,0.081274388361574/
203 DATA a_gauss /
204 1 0. ,0. ,0. ,
205 1 0. ,0. ,0. ,
206 1 0. ,0. ,0. ,
207 2 -.577350269189626,0.577350269189626,0. ,
208 2 0. ,0. ,0. ,
209 2 0. ,0. ,0. ,
210 3 -.774596669241483,0. ,0.774596669241483,
211 3 0. ,0. ,0. ,
212 3 0. ,0. ,0. ,
213 4 -.861136311594053,-.339981043584856,0.339981043584856,
214 4 0.861136311594053,0. ,0. ,
215 4 0. ,0. ,0. ,
216 5 -.906179845938664,-.538469310105683,0. ,
217 5 0.538469310105683,0.906179845938664,0. ,
218 5 0. ,0. ,0. ,
219 6 -.932469514203152,-.661209386466265,-.238619186083197,
220 6 0.238619186083197,0.661209386466265,0.932469514203152,
221 6 0. ,0. ,0. ,
222 7 -.949107912342759,-.741531185599394,-.405845151377397,
223 7 0. ,0.405845151377397,0.741531185599394,
224 7 0.949107912342759,0. ,0. ,
225 8 -.960289856497536,-.796666477413627,-.525532409916329,
226 8 -.183434642495650,0.183434642495650,0.525532409916329,
227 8 0.796666477413627,0.960289856497536,0. ,
228 9 -.968160239507626,-.836031107326636,-.613371432700590,
229 9 -.324253423403809,0. ,0.324253423403809,
230 9 0.613371432700590,0.836031107326636,0.968160239507626/
231
233 . hh(2,mvsiz),fac(mvsiz),c1,e0(mvsiz),dd(9,mvsiz),gg(mvsiz),
234 . dm(9,mvsiz),gm(9,mvsiz),dgm(9,mvsiz),dg(9,mvsiz),
235 . g33(9,mvsiz),hh1(2,mvsiz)
236
237
238
239
240 gbuf => elbuf_str%GBUF
241 nlay = elbuf_str%NLAY
242 IF (mpt == 222) THEN
243 nptr = 2
244 npts = 2
245 nptt = 2
246 ELSE
247 nptr = elbuf_str%NPTR
248 npts = elbuf_str%NPTS
249 nptt = elbuf_str%NPTT
250 ENDIF
251 iad0 = 1
252 IF (isorth > 0) iad0 = 1 + 6*nel
253 IF (igtyp == 21.OR.igtyp == 22) THEN
254 ikorth=2
255 ELSEIF (isorth>0) THEN
256 ikorth=1
257 ELSE
258 ikorth=0
259 ENDIF
260
261 nf1=nft+1
262
264 1 x, ixs(1,nf1),x1, x2,
265 2 x3, x4, x5, x6,
266 3 x7, x8, y1, y2,
267 4 y3, y4, y5, y6,
268 5 y7, y8, z1, z2,
269 6 z3, z4, z5, z6,
270 7 z7, z8, gbuf%OFF, offg,
271 8 gbuf%SMSTR,r11, r12, r13,
272 9 r21, r22, r23, r31,
273 a r32, r33, nc1, nc2,
274 b nc3, nc4, nc5, nc6,
275 c nc7, nc8, ngl, mxt,
276 d ngeo, k11, k12, k13,
277 e k14, k15, k16, k17,
278 f k18, k22, k23, k24,
279 g k25, k26, k27, k28,
280 h k33, k34, k35, k36,
281 i k37, k38, k44, k45,
282 j k46, k47, k48, k55,
283 k k56, k57, k58, k66,
284 l k67, k68, k77, k78,
285 m k88, jhbe, gbuf%GAMA, gama,
286 n nel, ismstr, jhbe, jcvt,
287 o irep, igtyp, isorth)
288
289 pid=ngeo(1)
290
291
292 DO i=1,nel
293 nu(i)=
min(half,pm(21,mxt(i)))
294 c1 =pm(32,mxt(i))
295 e0(i) =three*(one-two*nu(i))*c1
296
297 deltax(i)=ep30
298 ideg(i)=0
299 ENDDO
300 IF (jhbe/=14.AND.jhbe<20) THEN
301 DO i=1,nel
302 nu(i)=zero
303 ENDDO
304 ELSE
305 ENDIF
306 icpg = icp
307 IF (jhbe/=14.AND.jhbe/=17.AND.jhbe/=24) THEN
308 icpg = 1
309 ENDIF
310
312 1 offg, volg, ngl, x1,
313 2 x2, x3, x4, x5,
314 3 x6, x7, x8, y1,
315 4 y2, y3, y4, y5,
316 5 y6, y7, y8, z1,
317 6 z2, z3, z4, z5,
318 7 z6, z7, z8, pxc1,
319 8 pxc2, pxc3, pxc4, pyc1,
320 9 pyc2, pyc3, pyc4, pzc1,
321 a pzc2, pzc3, pzc4, px1h1,
322 b px1h2, px1h3, px1h4, px2h1,
323 c px2h2, px2h3, px2h4, px3h1,
324 d px3h2, px3h3, px3h4, px4h1,
325 e px4h2, px4h3, px4h4, hx,
326 f hy, hz, ajc1, ajc2,
327 g ajc3, ajc4, ajc5, ajc6,
328 h ajc7, ajc8, ajc9, smax,
329 i nel)
330 IF (icpg==2) THEN
331 CALL s8zsigp3(1 ,nel ,gbuf%SIG,e0 ,gbuf%PLA,
332 2 fac ,gbuf%G_PLA,nel )
333 DO i=1,nel
334 nu1(i)=nu(i)+(half-nu(i))*fac(i)
335 ENDDO
336 ELSEIF (jhbe==17.AND.icpg==1.AND.ismstr==10) THEN
337 DO i=1,nel
338 fac(i)=two*nu(i)
339 nu1(i)=nu(i)
340 ENDDO
341 ELSE
342 DO i=1,nel
343 nu1(i)=nu(i)
344 ENDDO
345 ENDIF
346 IF (mtn>=28) THEN
347 iadbuf = ipm(7,mxt(1))
348 ELSE
349 iadbuf = 1
350 ENDIF
351 nppt = npts*nptt*nptr*nlay
352 CALL mmats(1 ,nel ,pm ,mxt ,hh ,
353 . mtn ,ikorth ,ipm ,igeo ,gama ,
354 . bufmat(iadbuf) ,dm ,dgm ,gm ,
355 . jhbe ,gbuf%SIG ,bid ,nppt ,nel )
357 1 hh, hh1, fac, icpg,
358 2 ipredu, nel, mtn, ismstr,
359 3 jhbe)
360
361
362
363 il = 1
364 DO ir=1,nptr
365 DO is=1,npts
366 DO it=1,nptt
367 IF (jhbe == 14.OR.jhbe == 12.OR.jhbe == 17) THEN
368 lbuf => elbuf_str%BUFLY(il)%LBUF(ir,is,it)
369 ELSE
370 lbuf => elbuf_str%BUFLY(il)%LBUF(1,1,1)
371 ENDIF
372
373 ip = ir + ( (is-1) + (it-1)*npts )*nptr
374 wi = w_gauss(ir,nptr)*w_gauss(is,npts)*w_gauss(it,nptt)
375
377 1 offg, off, voln, ngl,
378 2 a_gauss(ir,nptr),a_gauss(is,npts),a_gauss(it,nptt),wi,
379 3 pxc1, pxc2, pxc3, pxc4,
380 4 pyc1, pyc2, pyc3, pyc4,
381 5 pzc1, pzc2, pzc3, pzc4,
382 6 px1h1, px1h2, px1h3, px1h4,
383 7 px2h1, px2h2, px2h3, px2h4,
384 8 px3h1, px3h2, px3h3, px3h4,
385 9 px4h1, px4h2, px4h3, px4h4,
386 a hx, hy, hz, px1,
387 b px2, px3, px4, px5,
388 c px6, px7, px8, py1,
389 d py2, py3, py4, py5,
390 e py6, py7, py8, pz1,
391 f pz2, pz3, pz4, pz5,
392 g pz6, pz7, pz8, pxy1,
393 h pxy2, pxy3, pxy4, pxy5,
394 i pxy6, pxy7, pxy8, pyx1,
395 j pyx2, pyx3, pyx4, pyx5,
396 k pyx6, pyx7, pyx8, pxz1,
397 l pxz2, pxz3, pxz4, pxz5,
398 m pxz6, pxz7, pxz8, pzx1,
399 n pzx2, pzx3, pzx4, pzx5,
400 o pzx6, pzx7, pzx8, pyz1,
401 p pyz2, pyz3, pyz4, pyz5,
402 q pyz6, pyz7, pyz8, pzy1,
403 r pzy2, pzy3, pzy4, pzy5,
404 s pzy6, pzy7, pzy8, bxy1,
405 t bxy2, bxy3, bxy4, bxy5,
406 u bxy6, bxy7, bxy8, byx1,
407 v byx2, byx3, byx4, byx5,
408 w byx6, byx7, byx8, bxz1,
409 x bxz2, bxz3, bxz4, bxz5,
410 y bxz6, bxz7, bxz8, bzx1,
411 z bzx2, bzx3, bzx4, bzx5,
412 1 bzx6, bzx7, bzx8, byz1,
413 2 byz2, byz3, byz4, byz5,
414 3 byz6, byz7, byz8, bzy1,
415 4 bzy2, bzy3, bzy4, bzy5,
416 5 bzy6, bzy7, bzy8, ajc1,
417 6 ajc2, ajc3, ajc4, ajc5,
418 7 ajc6, ajc7, ajc8, ajc9,
419 8 aj1, aj2, aj3, aj4,
420 9 aj5, aj6, smax, deltax,
421 a nu1, icpg, ideg, nel,
422 b mtn, ismstr, jhbe)
424 1 pm, mxt, hh1, voln,
425 2 icsig, dd, gg, dg,
426 3 g33, dm, gm, dgm,
427 4 ikorth, lbuf%SIG,ir, is,
428 5 it, nel, jhbe, mtn)
429
431 1 px1, px2, px3, px4,
432 2 py1, py2, py3, py4,
433 3 pz1, pz2, pz3, pz4,
434 4 px5, px6, px7, px8,
435 5 py5, py6, py7, py8,
436 6 pz5, pz6, pz7, pz8,
437 7 pxy1, pxy2, pxy3, pxy4,
438 8 pxy5, pxy6, pxy7, pxy8,
439 9 pyx1, pyx2, pyx3, pyx4,
440 a pyx5, pyx6, pyx7, pyx8,
441 b pxz1, pxz2, pxz3, pxz4,
442 c pxz5, pxz6, pxz7, pxz8,
443 d pzx1, pzx2, pzx3, pzx4,
444 e pzx5, pzx6, pzx7, pzx8,
445 f pyz1, pyz2, pyz3, pyz4,
446 g pyz5, pyz6, pyz7, pyz8,
447 h pzy1, pzy2, pzy3, pzy4,
448 i pzy5, pzy6, pzy7, pzy8,
449 j bxy1, bxy2, bxy3, bxy4,
450 k bxy5, bxy6, bxy7, bxy8,
451 l byx1, byx2, byx3, byx4,
452 m byx5, byx6, byx7, byx8,
453 n bxz1, bxz2, bxz3, bxz4,
454 o bxz5, bxz6, bxz7, bxz8,
455 p bzx1, bzx2, bzx3, bzx4,
456 q bzx5, bzx6, bzx7, bzx8,
457 r byz1, byz2, byz3, byz4,
458 s byz5, byz6, byz7, byz8,
459 t bzy1, bzy2, bzy3, bzy4,
460 u bzy5, bzy6, bzy7, bzy8,
461 v k11, k12, k13, k14,
462 w k15, k16, k17, k18,
463 x k22, k23, k24, k25,
464 y k26, k27, k28, k33,
465 z k34, k35, k36, k37,
466 1 k38, k44, k45, k46,
467 2 k47, k48, k55, k56,
468 3 k57, k58, k66, k67,
469 4 k68, k77, k78, k88,
470 5 dd, gg, icpg, offg,
471 6 ikorth, dg, g33, nel)
472 IF (ikgeo<0.AND.jhbe==14) THEN
474 1 lbuf%SIG,voln, px1, px2,
475 2 px3, px4, py1, py2,
476 3 py3, py4, pz1, pz2,
477 4 pz3, pz4, px5, px6,
478 5 px7, px8, py5, py6,
479 6 py7, py8, pz5, pz6,
480 7 pz7, pz8, k11, k12,
481 8 k13, k14, k15, k16,
482 9 k17, k18, k22, k23,
483 a k24, k25, k26, k27,
484 b k28, k33, k34, k35,
485 c k36, k37, k38, k44,
486 d k45, k46, k47, k48,
487 e k55, k56, k57, k58,
488 f k66, k67, k68, k77,
489 g k78, k88, nel)
490 ENDIF
491 ENDDO
492 ENDDO
493 ENDDO
494
495 IF (ipredu > 0) THEN
497 1 pxc1, pxc2, pxc3, pxc4,
498 2 pyc1, pyc2, pyc3, pyc4,
499 3 pzc1, pzc2, pzc3, pzc4,
500 4 k11, k12, k13, k14,
501 5 k15, k16, k17, k18,
502 6 k22, k23, k24, k25,
503 7 k26, k27, k28, k33,
504 8 k34, k35, k36, k37,
505 9 k38, k44, k45, k46,
506 a k47, k48, k55, k56,
507 b k57, k58, k66, k67,
508 c k68, k77, k78, k88,
509 d hh, volg, fac, icpg,
510 e nel)
511 END IF
512
513 IF (ikgeo>0) THEN
515 1 gbuf%SIG,volg, pxc1, pxc2,
516 2 pxc3, pxc4, pyc1, pyc2,
517 3 pyc3, pyc4, pzc1, pzc2,
518 4 pzc3, pzc4, k11, k12,
519 5 k13, k14, k15, k16,
520 6 k17, k18, k22, k23,
521 7 k24, k25, k26, k27,
522 8 k28, k33, k34, k35,
523 9 k36, k37, k38, k44,
524 a k45, k46, k47, k48,
525 b k55, k56, k57, k58,
526 c k66, k67, k68, k77,
527 d k78, k88, nel)
528 ENDIF
529
530
531
533 1 r11, r21, r31, r12,
534 2 r22, r32, r13, r23,
535 3 r33, k11, k12, k13,
536 4 k14, k15, k16, k17,
537 5 k18, k22, k23, k24,
538 6 k25, k26, k27, k28,
539 7 k33, k34, k35, k36,
540 8 k37, k38, k44, k45,
541 9 k46, k47, k48, k55,
542 a k56, k57, k58, k66,
543 b k67, k68, k77, k78,
544 c k88, nel)
545
547 1 1, nel, ixs(1,nf1), etag, offg)
549 1 ixs(1,nf1),nel ,iddl ,ndof ,k_diag,
550 2 k_lt ,iadk ,jdik ,k11 ,k12 ,
551 3 k13 ,k14 ,k15 ,k16 ,k17 ,
552 4 k18 ,k22 ,k23 ,k24 ,k25 ,
553 5 k26 ,k27 ,k28 ,k33 ,k34 ,
554 6 k35 ,k36 ,k37 ,k38 ,k44 ,
555 7 k45 ,k46 ,k47 ,k48 ,k55 ,
556 8 k56 ,k57 ,k58 ,k66 ,k67 ,
557 9 k68 ,k77 ,k78 ,k88 ,offg )
558
559 RETURN
subroutine assem_s8(ixs, nel, iddl, ndof, k_diag, k_lt, iadk, jdik, k11, k12, k13, k14, k15, k16, k17, k18, k22, k23, k24, k25, k26, k27, k28, k33, k34, k35, k36, k37, k38, k44, k45, k46, k47, k48, k55, k56, k57, k58, k66, k67, k68, k77, k78, k88, off)
subroutine mmat_h1(hh, hh1, fac, icp, ipredu, nel, mtn, ismstr, jhbe)
subroutine mmats(jft, jlt, pm, mat, hh, mtn, iorth, ipm, igeo, gama, uparam, cc, cg, g33, jhbe, sig, eps, nppt, nel)
subroutine mmstifs(pm, mat, hh, vol, icsig, dd, gg, dg, g33, dm, gm, dgm, iorth, sig, ir, is, it, nel, jhbe, mtn)
subroutine s8eoff(jft, jlt, ixs, etag, off)
subroutine s8zcumg3(r11, r21, r31, r12, r22, r32, r13, r23, r33, k11, k12, k13, k14, k15, k16, k17, k18, k22, k23, k24, k25, k26, k27, k28, k33, k34, k35, k36, k37, k38, k44, k45, k46, k47, k48, k55, k56, k57, k58, k66, k67, k68, k77, k78, k88, nel)
subroutine s8zderict3(off, det, ngl, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, px1h1, px1h2, px1h3, px1h4, px2h1, px2h2, px2h3, px2h4, px3h1, px3h2, px3h3, px3h4, px4h1, px4h2, px4h3, px4h4, hx, hy, hz, jac1, jac2, jac3, jac4, jac5, jac6, jac7, jac8, jac9, smax, nel)
subroutine s8zderii3(offg, off, vol, ngl, ksi, eta, zeta, wi, pxc1, pxc2, pxc3, pxc4, pyc1, pyc2, pyc3, pyc4, pzc1, pzc2, pzc3, pzc4, px1h1, px1h2, px1h3, px1h4, px2h1, px2h2, px2h3, px2h4, px3h1, px3h2, px3h3, px3h4, px4h1, px4h2, px4h3, px4h4, hx, hy, hz, px1, px2, px3, px4, px5, px6, px7, px8, py1, py2, py3, py4, py5, py6, py7, py8, pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8, pxy1, pxy2, pxy3, pxy4, pxy5, pxy6, pxy7, pxy8, pyx1, pyx2, pyx3, pyx4, pyx5, pyx6, pyx7, pyx8, pxz1, pxz2, pxz3, pxz4, pxz5, pxz6, pxz7, pxz8, pzx1, pzx2, pzx3, pzx4, pzx5, pzx6, pzx7, pzx8, pyz1, pyz2, pyz3, pyz4, pyz5, pyz6, pyz7, pyz8, pzy1, pzy2, pzy3, pzy4, pzy5, pzy6, pzy7, pzy8, bxy1, bxy2, bxy3, bxy4, bxy5, bxy6, bxy7, bxy8, byx1, byx2, byx3, byx4, byx5, byx6, byx7, byx8, bxz1, bxz2, bxz3, bxz4, bxz5, bxz6, bxz7, bxz8, bzx1, bzx2, bzx3, bzx4, bzx5, bzx6, bzx7, bzx8, byz1, byz2, byz3, byz4, byz5, byz6, byz7, byz8, bzy1, bzy2, bzy3, bzy4, bzy5, bzy6, bzy7, bzy8, cj1, cj2, cj3, cj4, cj5, cj6, cj7, cj8, cj9, jac4, jac5, jac6, jac7, jac8, jac9, smax, deltax, nu, icp, ideg, nel, mtn, ismstr, jhbe)
subroutine s8zkgei3(sig, vol, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, px5, px6, px7, px8, py5, py6, py7, py8, pz5, pz6, pz7, pz8, k11, k12, k13, k14, k15, k16, k17, k18, k22, k23, k24, k25, k26, k27, k28, k33, k34, k35, k36, k37, k38, k44, k45, k46, k47, k48, k55, k56, k57, k58, k66, k67, k68, k77, k78, k88, nel)
subroutine s8zkgeo3(sig, vol, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, k11, k12, k13, k14, k15, k16, k17, k18, k22, k23, k24, k25, k26, k27, k28, k33, k34, k35, k36, k37, k38, k44, k45, k46, k47, k48, k55, k56, k57, k58, k66, k67, k68, k77, k78, k88, nel)
subroutine s8zlke3(px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, px5, px6, px7, px8, py5, py6, py7, py8, pz5, pz6, pz7, pz8, pxy1, pxy2, pxy3, pxy4, pxy5, pxy6, pxy7, pxy8, pyx1, pyx2, pyx3, pyx4, pyx5, pyx6, pyx7, pyx8, pxz1, pxz2, pxz3, pxz4, pxz5, pxz6, pxz7, pxz8, pzx1, pzx2, pzx3, pzx4, pzx5, pzx6, pzx7, pzx8, pyz1, pyz2, pyz3, pyz4, pyz5, pyz6, pyz7, pyz8, pzy1, pzy2, pzy3, pzy4, pzy5, pzy6, pzy7, pzy8, bxy1, bxy2, bxy3, bxy4, bxy5, bxy6, bxy7, bxy8, byx1, byx2, byx3, byx4, byx5, byx6, byx7, byx8, bxz1, bxz2, bxz3, bxz4, bxz5, bxz6, bxz7, bxz8, bzx1, bzx2, bzx3, bzx4, bzx5, bzx6, bzx7, bzx8, byz1, byz2, byz3, byz4, byz5, byz6, byz7, byz8, bzy1, bzy2, bzy3, bzy4, bzy5, bzy6, bzy7, bzy8, k11, k12, k13, k14, k15, k16, k17, k18, k22, k23, k24, k25, k26, k27, k28, k33, k34, k35, k36, k37, k38, k44, k45, k46, k47, k48, k55, k56, k57, k58, k66, k67, k68, k77, k78, k88, dd, gg, icp, off, iksup, dg, g33, nel)
subroutine s8zlkep3(px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, k11, k12, k13, k14, k15, k16, k17, k18, k22, k23, k24, k25, k26, k27, k28, k33, k34, k35, k36, k37, k38, k44, k45, k46, k47, k48, k55, k56, k57, k58, k66, k67, k68, k77, k78, k88, hh, vol, fac, icp, nel)
subroutine s8zsigp3(lft, llt, sig, e0, defp, fac, g_pla, nel)
subroutine srcoork(x, ixs, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, offg, off, sav, r11, r12, r13, r21, r22, r23, r31, r32, r33, nc1, nc2, nc3, nc4, nc5, nc6, nc7, nc8, ngl, mxt, ngeo, k11, k12, k13, k14, k15, k16, k17, k18, k22, k23, k24, k25, k26, k27, k28, k33, k34, k35, k36, k37, k38, k44, k45, k46, k47, k48, k55, k56, k57, k58, k66, k67, k68, k77, k78, k88, khbe, gama0, gama, nel, ismstr, jhbe, jcvt, irep, igtyp, isorth)
1.15.0