60
61
62
63
64
65
66 USE elbufdef_mod
70 USE matparam_def_mod
72 use glob_therm_mod
73 use element_mod , only : nixs
74
75
76
77#include "implicit_f.inc"
78
79
80
81#include "mvsiz_p.inc"
82
83
84
85#include "com04_c.inc"
86#include "param_c.inc"
87#include "scr03_c.inc"
88#include "scr12_c.inc"
89#include "scr17_c.inc"
90#include "scry_c.inc"
91#include "vect01_c.inc"
92
93
94
95 INTEGER IXS(NIXS,*), IPARG(*),IPARTS(*),
96 . NEL, IPART(LIPART1,*),PERTURB(NPERTURB),
97 . IPM(NPROPMI,*), PTSOL(*), NSIGI, IUSER, NSIGS, NPF(*)
98 INTEGER IGEO(NPROPGI,*),STRSGLOB(*),STRAGLOB(*),ORTHOGLOB(*),
99 . FAIL_INI(*)
101 . mas(*), pm(npropm,*), x(*), geo(npropg,*),
102 . veul(lveul,*), dtelem(*),sigi(nsigs,*),skew(lskew,*),stifn(*),
103 . partsav(20,*), v(*), mss(8,*),
104 . sigsp(nsigi,*),msnf(*), mssf(8,*), wma(*),rnoise(nperturb,*),
105 . volnod(*), bvolnod(*), vns(8,*), bns(8,*),bufmat(*),
106 . mcp(*), mcps(8,*),temp(*), tf(*),xrefs(8,3,*), mssa(*)
107 TYPE(ELBUF_STRUCT_), TARGET :: ELBUF_STR
108 INTEGER,INTENT(IN) :: ILOADP(SIZLOADP,*)
109 my_real,
INTENT(IN) :: facload(lfacload,*)
110 TYPE(DETONATORS_STRUCT_)::DETONATORS
111 TYPE(t_ale_connectivity), INTENT(INOUT) :: ALE_CONNECTIVITY
112 TYPE (MATPARAM_STRUCT_) ,DIMENSION(NUMMAT) ,INTENT(INOUT) :: MAT_PARAM
113 type (glob_therm_) ,intent(in) :: glob_therm
114
115
116
117 INTEGER I,NF1,IBID,IGTYP,IP,IR,IS,IT,IL,NLAY,NPTR,NPTS,NPTT,NCC,
118 . JHBE,IREP,MPT,NUVAR,NUVARR,IDEF,NREFSTA,
119 . IPTHK, IPPOS,IG,IM,MTN0,ICSTR,IPID1,L_PLA,L_SIGB
120 INTEGER PID(MVSIZ), NGL(MVSIZ),MAT(MVSIZ), MAT0(MVSIZ),
121 . IX1(MVSIZ), IX2(MVSIZ), IX3(MVSIZ), IX4(MVSIZ),
122 . IX5(MVSIZ), IX6(MVSIZ), IX7(MVSIZ), IX8(MVSIZ)
124 . x1(mvsiz),x2(mvsiz),x3(mvsiz),x4(mvsiz),x5(mvsiz),x6(mvsiz),
125 . x7(mvsiz),x8(mvsiz),y1(mvsiz),y2(mvsiz),y3(mvsiz),y4(mvsiz),
126 . y5(mvsiz),y6(mvsiz),y7(mvsiz),y8(mvsiz),z1(mvsiz),z2(mvsiz),
127 . z3(mvsiz),z4(mvsiz),z5(mvsiz),z6(mvsiz),z7(mvsiz),z8(mvsiz),
128 . rx(mvsiz) ,ry(mvsiz) ,rz(mvsiz) ,sx(mvsiz) ,
129 . sy(mvsiz) ,sz(mvsiz) ,tx(mvsiz) ,ty(mvsiz) ,tz(mvsiz) ,
130 . e1x(mvsiz),e1y(mvsiz),e1z(mvsiz),e2x(mvsiz),
131 . e2y(mvsiz),e2z(mvsiz),e3x(mvsiz),e3y(mvsiz),e3z(mvsiz),
132 . f1x(mvsiz) ,f1y(mvsiz) ,f1z(mvsiz) ,
133 . f2x(mvsiz) ,f2y(mvsiz) ,f2z(mvsiz) ,
134 . ajc1(mvsiz) , ajc2(mvsiz) , ajc3(mvsiz) ,
135 . ajc4(mvsiz) , ajc5(mvsiz) , ajc6(mvsiz) ,
136 . ajc7(mvsiz) , ajc8(mvsiz) , ajc9(mvsiz) ,
137 . hx(4,mvsiz) , hy(4,mvsiz), hz(4,mvsiz),gama(6,mvsiz),
138 . smax(mvsiz) , volu(mvsiz), dtx(mvsiz), deltax(mvsiz),
139 . rhocp(mvsiz),temp0(mvsiz), aire(mvsiz),llsh(mvsiz)
141 . bid(mvsiz), fv, sti, wi
143 . angle(mvsiz),dtx0(mvsiz),wt,zr,zs,zt,zz
145 DOUBLE PRECISION
146 . XD1(MVSIZ), XD2(MVSIZ), XD3(), XD4(MVSIZ),
147 . XD5(MVSIZ), XD6(MVSIZ), XD7(MVSIZ), XD8(MVSIZ),
148 . YD1(MVSIZ), YD2(MVSIZ), YD3(MVSIZ), YD4(MVSIZ),
149 . YD5(MVSIZ), YD6(MVSIZ), (MVSIZ), YD8(MVSIZ),
150 . ZD1(MVSIZ), ZD2(MVSIZ), ZD3(MVSIZ), (MVSIZ),
151 . ZD5(MVSIZ), ZD6(MVSIZ), ZD7(MVSIZ), (MVSIZ)
152 INTEGER NLYMAX, IPMAT,IPANG
153 CHARACTER(LEN=NCHARTITLE)::TITR
154 parameter(nlymax = 200,ipmat = 100,ipang = 200)
155
156 TYPE() ,POINTER :: LBUF
157 TYPE(G_BUFEL_) ,POINTER :: GBUF
158 TYPE(BUF_MAT_) ,POINTER ::
159
161 . w_gauss(9,9),a_gauss(9,9)
162 DATA w_gauss /
163 1 2. ,0. ,0. ,
164 1 0. ,0. ,0. ,
165 1 0. ,0. ,0. ,
166 2 1. ,1. ,0. ,
167 2 0. ,0. ,0. ,
168 2 0. ,0. ,0. ,
169 3 0.555555555555556,0.888888888888889,0.555555555555556,
170 3 0. ,0. ,0. ,
171 3 0. ,0. ,0. ,
172 4 0.347854845137454,0.652145154862546,0.652145154862546,
173 4 0.347854845137454,0. ,0. ,
174 4 0. ,0. ,0. ,
175 5 0.236926885056189,0.478628670499366,0.568888888888889,
176 5 0.478628670499366,0.236926885056189,0. ,
177 5 0. ,0. ,0. ,
178 6 0.171324492379170,0.360761573048139,0.467913934572691,
179 6 0.467913934572691,0.360761573048139,0.171324492379170,
180 6 0. ,0. ,0. ,
181 7 0.129484966168870,0.279705391489277,0.381830050505119,
182 7 0.417959183673469,0.381830050505119,0.279705391489277,
183 7 0.129484966168870,0. ,0. ,
184 8 0.101228536290376,0.222381034453374,0.313706645877887,
185 8 0.362683783378362,0.362683783378362,0.313706645877887,
186 8 0.222381034453374,0.101228536290376,0. ,
187 9 0.081274388361574,0.180648160694857,0.260610696402935,
188 9 0.312347077040003,0.330239355001260,0.312347077040003,
189 9 0.260610696402935,0.180648160694857,0.081274388361574/
190 DATA a_gauss /
191 1 0. ,0. ,0. ,
192 1 0. ,0. ,0. ,
193 1 0. ,0. ,0. ,
194 2 -.577350269189626,0.577350269189626,0. ,
195 2 0. ,0. ,0. ,
196 2 0. ,0. ,0. ,
197 3 -.774596669241483,0. ,0.774596669241483,
198 3 0. ,0. ,0. ,
199 3 0. ,0. ,0. ,
200 4 -.861136311594053,-.339981043584856,0.339981043584856,
201 4 0.861136311594053,0. ,0. ,
202 4 0. ,0. ,0. ,
203 5 -.906179845938664,-.538469310105683,0. ,
204 5 0.538469310105683,0.906179845938664,0. ,
205 5 0. ,0. ,0. ,
206 6 -.932469514203152,-.661209386466265,-.238619186083197,
207 6 0.238619186083197,0.661209386466265,0.932469514203152,
208 6 0. ,0. ,0. ,
209 7 -.949107912342759,-.741531185599394,-.405845151377397,
210 7 0. ,0.405845151377397,0.741531185599394,
211 7 0.949107912342759,0. ,0. ,
212 8 -.960289856497536,-.796666477413627,-.525532409916329,
213 8 -.183434642495650,0.183434642495650,0.525532409916329,
214 8 0.796666477413627,0.960289856497536,0. ,
215 9 -.968160239507626,-.836031107326636,-.613371432700590,
216 9 -.324253423403809,0. ,0.324253423403809,
217 9 0.613371432700590,0.836031107326636,0.968160239507626/
218
219
220
221 gbuf => elbuf_str%GBUF
222 mbuf => elbuf_str%BUFLY(1)%MAT(1,1,1)
223 lbuf => elbuf_str%BUFLY(1)%LBUF(1,1,1)
224
225 bid(1:mvsiz) = zero
226 nrefsta = nxref
227 nxref = 0
228 DO i=1,nel
229 deltax(i)=ep30
230 ENDDO
231 jhbe = iparg(23)
232 irep = iparg(35)
233 igtyp = iparg(38)
234 IF (jcvt==1.AND.isorth/=0) jcvt=2
235
236 nf1=nft+1
237 IF (igtyp /= 22) isorth = 0
238 icstr=iparg(17)
239
240 DO i=1,nel
241 rhocp(i) = pm(69,ixs(1,nft+i))
242 temp0(i) = pm(79,ixs(1,nft+i))
243 ENDDO
244
245 IF (jcvt == 0) THEN
246 CALL scoor3(x ,bid(1),ixs(1,nf1) ,geo ,mat ,pid ,ngl ,
247 . ix1 ,ix2 ,ix3 ,ix4 ,ix5 ,ix6 ,ix7 ,ix8 ,
248 . x1 ,x2 ,x3 ,x4 ,x5 ,x6 ,x7 ,x8 ,
249 . y1 ,y2 ,y3 ,y4 ,y5 ,y6 ,y7 ,y8 ,
250 . z1 ,z2 ,z3 ,z4 ,z5 ,z6 ,z7 ,z8 ,
251 . rx ,ry ,rz ,sx ,sy ,sz ,tx ,ty ,tz ,
252 . e1x ,e1y ,e1z ,e2x ,e2y ,e2z ,e3x ,e3y ,e3z ,
253 . f1x ,f1y ,f1z ,f2x ,f2y ,f2z ,temp0, temp,glob_therm%NINTEMP,
254 . xd1 ,xd2 ,xd3 ,xd4 ,xd5 ,xd6 ,xd7 ,xd8 ,
255 . yd1 ,yd2 ,yd3 ,yd4 ,yd5 ,yd6 ,yd7 ,yd8 ,
256 . zd1 ,zd2 ,zd3 ,zd4 ,zd5 ,zd6 ,zd7 ,zd8 )
257 ELSE
258 CALL srcoor3(x,bid(1),ixs(1,nf1) ,geo ,mat ,pid ,ngl ,jhbe ,
259 . ix1 ,ix2 ,ix3 ,ix4 ,ix5 ,ix6 ,ix7 ,ix8 ,
260 . x1 ,x2 ,x3 ,x4 ,x5 ,x6 ,x7 ,x8 ,
261 . y1 ,y2 ,y3 ,y4 ,y5 ,y6 ,y7 ,y8 ,
262 . z1 ,z2 ,z3 ,z4 ,z5 ,z6 ,z7 ,z8 ,
263 . rx ,ry ,rz ,sx ,sy ,sz ,tx ,ty ,tz ,
264 . e1x ,e1y ,e1z ,e2x ,e2y ,e2z ,e3x ,e3y ,e3z ,
265 . f1x ,f1y ,f1z ,f2x ,f2y ,f2z ,temp0,temp,glob_therm%NINTEMP,
266 . xd1 ,xd2 ,xd3 ,xd4 ,xd5 ,xd6 ,xd7 ,xd8 ,
267 . yd1 ,yd2 ,yd3 ,yd4 ,yd5 ,yd6 ,yd7 ,yd8 ,
268 . zd1 ,zd2 ,zd3 ,zd4 ,zd5 ,zd6 ,zd7 ,zd8 )
269 ENDIF
270
271 SELECT CASE (igtyp)
272
273 CASE(21)
274 DO i=1,nel
275 angle(i) = geo(1,pid(i))
276 END DO
277 CALL scmorth3(pid ,geo ,igeo ,skew ,irep ,gbuf%GAMA ,
278 . rx ,ry ,rz ,sx ,sy ,sz ,tx ,ty ,tz ,
279 . e1x ,e1y ,e1z ,e2x ,e2y ,e2z ,e3x ,e3y ,e3z ,
280 . ngl ,angle,nsigi,sigsp,nsigs,sigi ,ixs(1,nf1) ,1 ,
281 . orthoglob ,ptsol,nel)
282
283 CASE(22)
284 DO i=1,nel
285 angle(i) = geo(1,pid(i))
286 END DO
287 CALL scmorth3(pid ,geo ,igeo ,skew ,irep ,gbuf%GAMA ,
288 . rx ,ry ,rz ,sx ,sy ,sz ,tx ,ty ,tz ,
289 . e1x ,e1y ,e1z ,e2x ,e2y ,e2z ,e3x ,e3y ,e3z ,
290 . ngl ,angle ,nsigi,sigsp,nsigs,sigi ,ixs(1,nf1),1 ,
291 . orthoglob ,ptsol,nel)
292 ipthk = ipang+nlymax
293 ippos = ipthk+nlymax
294 ig = pid(1)
295 mtn0 = mtn
296 DO i=1,nel
297 mat0(i) = mat(i)
298 dtx0(i) = ep20
299 ENDDO
300 END SELECT
301
303 . ajc1,ajc2,ajc3,
304 . ajc4,ajc5,ajc6,
305 . ajc7,ajc8,ajc9,smax, volu, ngl,
306 . xd1 ,xd2 ,xd3 ,xd4 ,xd5 ,xd6 ,xd7 ,xd8 ,
307 . yd1 ,yd2 ,yd3 ,yd4 ,yd5 ,yd6 ,yd7 ,yd8 ,
308 . zd1 ,zd2 ,zd3 ,zd4 ,zd5 ,zd6 ,zd7 ,zd8 )
309 IF (idttsh > 0) THEN
311 . x1, x2, x3, x4, x5, x6, x7, x8,
312 . y1, y2, y3, y4, y5, y6, y7, y8,
313 . z1, z2, z3, z4, z5, z6, z7, z8,icstr,idt1sol)
314 END IF
315
316
317
318
319 IF (jthe == 0 .and. glob_therm%NINTEMP > 0) THEN
320 DO i=1,nel
321 tempel(i) = one_over_8 *(temp(ixs(2,i)) + temp(ixs(3,i))
322 . + temp(ixs(4,i)) + temp(ixs(5,i))
323 . + temp(ixs(6,i)) + temp(ixs(7,i))
324 . + temp(ixs(8,i)) + temp(ixs(9,i)))
325 ENDDO
326 ELSE
327 tempel(1:nel) = temp0(1:nel)
328 END IF
329
330 ip=0
331 CALL matini(pm ,ixs ,nixs ,x ,
332 . geo ,ale_connectivity ,detonators ,iparg ,
333 . sigi ,nel ,skew ,igeo ,
334 . ipart ,iparts ,
335 . mat ,ipm ,nsigs ,numsol ,ptsol ,
336 . ip ,ngl ,npf ,tf ,bufmat ,
337 . gbuf ,lbuf ,mbuf ,elbuf_str ,iloadp ,
338 . facload, deltax ,tempel ,mat_param )
339
340 CALL sczero3(gbuf%RHO,gbuf%SIG,gbuf%EINT,nel)
341
342
343 IF(jthe /=0)
CALL atheri(mat,pm,gbuf%TEMP)
344
345 nlay = elbuf_str%NLAY
346 nptr = elbuf_str%NPTR
347 npts = elbuf_str%NPTS
348 nptt = elbuf_str%NPTT
349 it = 1
350
351
352 DO ir=1,nptr
353 DO is=1,npts
354 DO il=1,nlay
355
356 lbuf => elbuf_str%BUFLY(il)%LBUF(ir,is,it)
357 mbuf => elbuf_str%BUFLY(il)%MAT(ir,is,it)
358 l_pla = elbuf_str%BUFLY(il)%L_PLA
359 l_sigb= elbuf_str%BUFLY(il)%L_SIGB
360
361 IF (igtyp == 22) THEN
362 wt = geo(ipthk+il,ig)
363 zz = geo(ippos+il,ig)
364 im =igeo(ipmat+il,ig)
365 mtn=nint(pm(19,im))
366 DO i=1,nel
367 mat(i)=im
368 angle(i) = geo(ipang+il,pid(i))
369 ENDDO
370 ELSE
371 zz = a_gauss(il,nlay)
372 wt = w_gauss(il,nlay)
373 ENDIF
374
375 IF (icstr == 10) THEN
376 zr = a_gauss(ir,nptr)
377 zs = a_gauss(is,npts)
378 zt = zz
379 ELSEIF (icstr == 100) THEN
380 zr = a_gauss(ir,nptr)
381 zs = zz
382 zt = a_gauss(is,npts)
383 ELSEIF (icstr == 1) THEN
384 zr = zz
385 zs = a_gauss(ir,nptr)
386 zt = a_gauss(is,npts)
387 ENDIF
388 ip = ir + ( (is-1) + (il-1)*npts )*nptr
389 wi = w_gauss(ir,nptr)*w_gauss(is,npts)*wt
390
391 CALL s8zderi3(lbuf%VOL,veul(1,nf1),geo,
392 . zr,zs,zt,wi,
393 . hx, hy, hz,
394 . ajc1,ajc2,ajc3,
395 . ajc4,ajc5,ajc6,
396 . ajc7,ajc8,ajc9,smax, deltax, ngl,lbuf%VOL0DP)
397 IF (idttsh > 0) THEN
398 DO i=1,nel
399 IF (gbuf%IDT_TSH(i)>0)
400 . deltax(i)=
max(llsh(i),deltax(i))
401 ENDDO
402 END IF
403 IF (igtyp == 22)
404 .
CALL scmorth3(pid ,geo ,igeo ,skew ,irep ,lbuf%GAMA,
405 . rx ,ry ,rz ,sx ,sy ,sz ,tx ,ty ,tz ,
406 . e1x ,e1y ,e1z ,e2x ,e2y ,e2z ,e3x ,e3y ,e3z ,
407 . ngl ,angle ,nsigi,sigsp,nsigs,sigi ,ixs(1,nf1),il ,
408 . orthoglob, ptsol,nel)
409
410 IF (jthe == 0 .and. glob_therm%NINTEMP > 0) THEN
411 DO i=1,nel
412 tempel(i) = one_over_8 *(temp(ixs(2,i)) + temp(ixs(3,i))
413 . + temp(ixs(4,i)) + temp(ixs(5,i))
414 . + temp(ixs(6,i)) + temp(ixs(7,i))
415 . + temp(ixs(8,i)) + temp(ixs(9,i)))
416 ENDDO
417 ELSE
418 tempel(1:nel) = temp0(1:nel)
419 END IF
420
421 CALL matini(pm ,ixs ,nixs ,x ,
422 . geo ,ale_connectivity ,detonators,iparg ,
423 . sigi ,nel ,skew ,igeo ,
424 . ipart ,iparts ,
425 . mat ,ipm ,nsigs ,numsol ,ptsol ,
426 . ip ,ngl ,npf ,tf ,bufmat ,
427 . gbuf ,lbuf ,mbuf ,elbuf_str ,iloadp ,
428 . facload,deltax ,tempel ,mat_param )
429
430 idef =0
431 IF (mtn >= 28) THEN
432 nuvar = ipm(8,ixs(1,nft+1))
433 idef =1
434 ELSE
435 nuvar = 0
436 IF (mtn == 14 .OR. mtn == 12 .OR. mtn == 24) THEN
437 idef =1
438 ELSEIF (istrain == 1 .AND.
439 . (mtn == 1 .OR. mtn == 2 .OR. mtn == 3 .OR.
440 . mtn == 4 .OR. mtn == 6 .OR. mtn == 10 .OR.
441 . mtn == 21 .OR. mtn == 22 .OR. mtn == 23 .OR.
442 . mtn == 49)) THEN
443 idef =1
444 ENDIF
445 ENDIF
446
448 . lbuf%SIG ,pm ,lbuf%VOL ,sigsp ,
449 . sigi ,lbuf%EINT,lbuf%RHO ,mbuf%VAR ,lbuf%STRA,
450 . ixs ,nixs ,nsigi ,ip ,nuvar ,
451 . nel ,iuser ,idef ,nsigs ,strsglob ,
452 . straglob ,jhbe ,igtyp ,x ,gbuf%GAMA,
453 . mat ,lbuf%PLA ,l_pla ,ptsol ,lbuf%SIGB,
454 . l_sigb ,ipm ,bufmat ,lbuf%VOL0DP)
455
456 IF (igtyp == 22) THEN
457 aire(:) = zero
458 CALL dtmain(geo ,pm ,ipm ,pid ,mat ,fv ,
459 . lbuf%EINT ,lbuf%TEMP ,lbuf%DELTAX ,lbuf%RK ,lbuf%RE ,bufmat, deltax, aire,
460 . volu, dtx,igeo,igtyp)
461
463 . lbuf%RHO,lbuf%VOL,lbuf%OFF,lbuf%SIG,lbuf%EINT,dtx,
464 . gbuf%RHO,gbuf%VOL,gbuf%OFF,gbuf%SIG,gbuf%EINT,dtx0,
465 . nel )
466 ELSE
468 . lbuf%RHO,lbuf%VOL,lbuf%OFF,lbuf%SIG,lbuf%EINT,dtx,
469 . gbuf%RHO,gbuf%VOL,gbuf%OFF,gbuf%SIG,gbuf%EINT,dtx0,
470 . nel )
471 ENDIF
472
473 ENDDO
474 ENDDO
475 ENDDO
476
477 IF (igtyp == 22) THEN
478 mtn=mtn0
479 DO i=1,nel
480 mat(i)=mat0(i)
481 ENDDO
482 ENDIF
483
484
486 . gbuf%RHO,mas,partsav,x,v,
487 . iparts(nf1),mss(1,nf1),volu ,
488 . msnf ,mssf(1,nf1) ,bid(1) ,
489 . bid(1) ,bid(1) ,wma ,rhocp ,mcp ,
490 . mcps(1,nf1) ,mssa ,bid(1) ,bid(1),gbuf%FILL,
491 . ix1, ix2, ix3, ix4, ix5, ix6, ix7, ix8)
492
493
494
495 IF (i7stifs /= 0) THEN
496 ncc=8
497 CALL sbulk3(volu ,ix1 ,ncc,mat,pm ,
498 . volnod,bvolnod,vns(1,nf1),bns(1,nf1),bid(1),
499 . bid(1) ,gbuf%FILL)
500 ENDIF
501
502
503 CALL failini(elbuf_str,nptr,npts,nptt,nlay,
504 . ipm,sigsp,nsigi,fail_ini ,
505 . sigi,nsigs,ixs,nixs,ptsol,rnoise,perturb,mat_param)
506
507
508 aire(:) = zero
509 CALL dtmain(geo ,pm ,ipm ,pid ,mat ,fv ,
510 . lbuf%EINT ,lbuf%TEMP ,lbuf%DELTAX ,lbuf%RK ,lbuf%RE ,bufmat, deltax, aire,
511 . volu, dtx,igeo,igtyp)
512
513 IF (igtyp == 22) THEN
514 DO i=1,nel
515 dtx(i)=dtx0(i)
516 ENDDO
517 ENDIF
518 DO i=1,nel
519 IF (ixs(10,i+nft)/=0.AND.invers>14) THEN
520 IF(igtyp/=0.AND.igtyp/=6.AND.igtyp/=14.AND.igtyp/=15
521 . .AND.igtyp/=20.AND.igtyp/=21.AND.igtyp/=22)THEN
522 ipid1=ixs(nixs-1,i+nft)
523 CALL fretitl2(titr,igeo(npropgi-ltitr+1,ipid1),ltitr)
525 . msgtype=msgerror,
526 . anmode=aninfo_blind_1,
527 . i1=igeo(1,ipid1),
528 . c1=titr,
529 . i2=igtyp)
530 ENDIF
531 ENDIF
532 dtelem(nft+i)=dtx(i)
533
534 sti = fourth * gbuf%FILL(i) * gbuf%RHO(i) * volu(i)
535 . /
max(em20,dtx(i)*dtx(i))
536 stifn(ixs(2,i+nft))=stifn(ixs(2,i+nft))+sti
537 stifn(ixs(3,i+nft))=stifn(ixs(3,i+nft))+sti
538 stifn(ixs(4,i+nft))=stifn(ixs(4,i+nft))+sti
539 stifn(ixs(5,i+nft))=stifn(ixs(5,i+nft))+sti
540 stifn(ixs(6,i+nft))=stifn(ixs(6,i+nft))+sti
541 stifn(ixs(7,i+nft))=stifn(ixs(7,i+nft))+sti
542 stifn(ixs(8,i+nft))=stifn(ixs(8,i+nft))+sti
543 stifn(ixs(9,i+nft))=stifn(ixs(9,i+nft))+sti
544 ENDDO
545
546 nxref = nrefsta
547
548 RETURN
subroutine atheri(mat, pm, temp)
subroutine dtmain(geo, pm, ipm, pid, mat, fv, eint, temp, deltax, rk, re, bufmat, ddeltax, aire, vol, dtx, igeo, igtyp)
subroutine failini(elbuf_str, nptr, npts, nptt, nlay, ipm, sigsp, nsigi, fail_ini, sigi, nsigs, ix, nix, pt, rnoise, perturb, mat_param)
subroutine matini(pm, ix, nix, x, geo, ale_connectivity, detonators, iparg, sigi, nel, skew, igeo, ipart, ipartel, mat, ipm, nsig, nums, pt, ipt, ngl, npf, tf, bufmat, gbuf, lbuf, mbuf, elbuf_str, iloadp, facload, ddeltax, tempel, mat_param)
integer, parameter nchartitle
subroutine sigin20b(sig, pm, vol, sigsp, sigi, eint, rho, uvar, eps, ix, nix, nsigi, ipt, nuvar, nel, iuser, idef, nsigs, strsglob, straglob, jhbe, igtyp, x, bufgama, mat, epsp, l_pla, pt, sigb, l_sigb, ipm, bufmat, voldp)
subroutine sbulk3(volu, nc, nnc, mat, pm, volnod, bvolnod, vns, bns, vnsx, bnsx, fill)
subroutine sczero3(rhog, sigg, eintg, nel)
subroutine svalue0(rho, vol, off, sig, eint, dtx, rhog, volg, offg, sigg, eintg, dtxg, nel)
subroutine scmorth3(pid, geo, igeo, skew, irep, gama, rx, ry, rz, sx, sy, sz, tx, ty, tz, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, ngl, angle, nsigi, sigsp, nsigs, sigi, ixs, ilay, orthoglob, pt, nel)
subroutine smass3(rho, ms, partsav, x, v, ipart, mss, volu, msnf, mssf, in, vr, ins, wma, rhocp, mcp, mcps, mssa, rhof, frac, fill, nc1, nc2, nc3, nc4, nc5, nc6, nc7, nc8)
subroutine s8zderi3(vol, veul, geo, ksi, eta, zeta, wi, hx, hy, hz, cj1, cj2, cj3, cj4, cj5, cj6, cj7, cj8, cj9, smax, deltax, ngl, voldp)
subroutine s8zderic3(vol, hx, hy, hz, jac1, jac2, jac3, jac4, jac5, jac6, jac7, jac8, jac9, smax, det, ngl, xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8)
subroutine scoor3(x, xrefs, ixs, geo, mxt, ngeo, ngl, ix1, ix2, ix3, ix4, ix5, ix6, ix7, ix8, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, rx, ry, rz, sx, sy, sz, tx, ty, tz, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, f1x, f1y, f1z, f2x, f2y, f2z, temp0, temp, nintemp, xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8)
subroutine srcoor3(x, xrefs, ixs, geo, mxt, ngeo, ngl, jhbe, ix1, ix2, ix3, ix4, ix5, ix6, ix7, ix8, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, rx, ry, rz, sx, sy, sz, tx, ty, tz, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, f1x, f1y, f1z, f2x, f2y, f2z, temp0, temp, nintemp, xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8)
subroutine sdlensh14(nel, llsh, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, ics, idt1sol)
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
1.15.0