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