59 1 IPARI ,X ,A ,IGROUPS ,ALE_CONNECTIVITY,
60 2 ICODT ,FSAV ,V ,MS ,DT2T ,
61 3 NELTST ,ITYPTST ,ITAB ,STIFN ,TF ,
62 4 FSKYI ,ISKY ,VR ,FCONT ,SECFCUM ,
63 5 JTASK ,NISKYFI ,KINET ,NEWFRONT ,NSTRF ,
64 6 ICONTACT ,VISCN ,XCELL ,
65 8 NUM_IMP ,NS_IMP ,NE_IMP ,IND_IMP ,NT_IMP ,
66 9 FR_I18 ,IGRBRIC ,EMINX ,
67 A IXS ,IXS16 ,IXS20 ,FNCONT ,FTCONT ,
68 B IAD_ELEM ,FR_ELEM ,RCONTACT ,ACONTACT ,PCONTACT ,
69 C TEMP ,FTHE ,FTHESKYI ,IPARG ,NSENSOR ,
70 D PM ,INTSTAMP ,WEIGHT ,NISKYFIE ,IRLEN20 ,
71 E ISLEN20 ,IRLEN20T ,ISLEN20T ,IRLEN20E ,ISLEN20E ,
72 F MSKYI_SMS ,ISKYI_SMS ,NODNX_SMS ,INT18ADD ,FCONTG ,
73 G FNCONTG ,FTCONTG ,NODGLOB ,MS0 ,NPC ,
74 H WA ,SENSOR_TAB ,QFRICINT ,NCONT ,INDEXCONT ,
75 I TAGCONT ,INOD_PXFEM ,MS_PLY ,WAGAP ,ELBUF_TAB ,
76 J CONDN ,CONDNSKYI ,NV46 ,
77 K SFBSAV6 ,FBSAV6 ,NODADT_THERM,THEACCFACT,
78 L ISENSINT ,NISUBMAX ,NB25_CANDT,NB25_IMPCT ,
79 M NB25_DST1 ,NB25_DST2 ,IXIG3D ,KXIG3D ,WIGE ,
80 N KNOT ,IGEO ,MULTI_FVM ,H3D_DATA ,
81 P PSKIDS ,T2MAIN_SMS ,FORNEQS ,KNOTLOCPC,KNOTLOCEL,
82 Q APINCH ,STIFPINCH ,T2FAC_SMS ,TAGNCONT ,
83 R KLOADPINTER,LOADPINTER ,LOADP_HYD_INTER,DGAPLOADINT,S_LOADPINTER,
84 S INTERFACES,XCELL_REMOTE)
106#include "implicit_f.inc"
107#include "comlock.inc"
111#include "com01_c.inc"
112#include "com04_c.inc"
113#include "com08_c.inc"
115#include "param_c.inc"
117#include "units_c.inc"
118#include "parit_c.inc"
119#include "impl1_c.inc"
120#include "timeri_c.inc"
121#include "tabsiz_c.inc"
125 TYPE(output_),
intent(inout) :: output
126 TYPE(timer_),
INTENT(INOUT) :: TIMERS
127 INTEGER ,
INTENT(IN) :: NSENSOR
128 INTEGER NELTST,ITYPTST,NEWFRONT(*),NSTRF(*),
129 . IRLEN20, ISLEN20, IRLEN20T, ISLEN20T,
130 . IRLEN20E, ISLEN20E,INT18ADD(*),NCONT,NV46,SFBSAV6,IGROUPS(NUMELS)
131 INTEGER IPARI(NPARI,*), ICODT(*),ITAB(*),
132 . ISKY(*), JTASK, NISKYFI(*),ICONTACT(*),
133 . IPARG(NPARG,*),INOD_PXFEM(*),TAGNCONT(NLOADP_HYD_INTER,*)
135 INTEGER,
TARGET :: FR_I18(NSPMD+2,*)
136 INTEGER,
TARGET :: KINET(*)
137 INTEGER IAD_ELEM(2,*),FR_ELEM(*), NISKYFIE(*)
138 INTEGER NUM_IMP(*),NS_IMP(*),NE_IMP(*),IND_IMP(*),NT_IMP
139 INTEGER IXS(*) ,IXS16(*) ,IXS20(*)
140 INTEGER WEIGHT(*), ISKYI_SMS(*), NODNX_SMS(*), NODGLOB(*), NPC(*)
141 INTEGER INDEXCONT(*),TAGCONT(*),
143 . NB25_CANDT(PARASIZ),NB25_IMPCT(PARASIZ),
144 . NB25_DST1(PARASIZ) ,NB25_DST2(PARASIZ)
145 INTEGER KXIG3D(NIXIG3D,*),IXIG3D(*),IGEO(NPROPGI,*),T2MAIN_SMS(6,*)
146 INTEGER ,
INTENT(IN) :: S_LOADPINTER
147 INTEGER ,
INTENT(IN) :: KLOADPINTER(NINTER+1),LOADPINTER(S_LOADPINTER),
148 . LOADP_HYD_INTER(NLOADP_HYD)
149 INTEGER ,
INTENT(IN) :: NODADT_THERM
150 my_real ,
intent(in) :: theaccfact
151 my_real ,
INTENT(IN) :: dgaploadint(s_loadpinter)
156 . stifn(*), tf(*),fskyi(lskyi,nfskyi),vr(3,*),fcont(3,*),
157 . secfcum(7,numnod,nsect),
158 . viscn(*), fncont(3,*),ftcont
159 . ftheskyi(lskyi),pm(npropm,*),fthe(*),
160 . acontact(*), pcontact(*), mskyi_sms(*),
161 . fcontg(3,*),fncontg(3,*),ftcontg(3,*),wa(*),
162 . qfricint(*),ms_ply(*),wagap(*),condn(*),condnskyi(lskyi),
163 . wige(*),knot(*),pskids(*),forneqs(3,*),knotlocpc(*),knotlocel(*),
164 . apinch(3,*),stifpinch(*),t2fac_sms(*)
165 my_real,
TARGET :: x(3*numnod),v(3*numnod),ms(numnod),xcell(3,sxcell)
167 TYPE (ELBUF_STRUCT_),
DIMENSION(NGROUP)
168DOUBLE PRECISION FBSAV6(12,6,SFBSAV6)
169 INTEGER ISENSINT(NISUBMAX+1,NINTER)
170 TYPE(MULTI_FVM_STRUCT),
INTENT(INOUT),
TARGET :: MULTI_FVM
173 TYPE (SENSOR_STR_) ,
DIMENSION(NSENSOR) :: SENSOR_TAB
174 TYPE (INTERFACES_) ,
INTENT(INOUT):: INTERFACES
175 TYPE(
array_type),
DIMENSION(NINTER),
INTENT(in) :: XCELL_REMOTE
177 TYPE (GROUP_) ,
DIMENSION(NGRBRIC) :: IGRBRIC
181 INTEGER N, LINDMAX,IAD,IDUM,NRTMDIM, NI18,
182 . IAD17, IGN, IGE, NME, NMES, ISTAMP,NN1,NN2,ISENS
184 INTEGER NB_JLT(PARASIZ),NB_JLT_NEW(PARASIZ),NB_STOK_N(PARASIZ),NTH
186 SAVE NB_JLT,NB_JLT_NEW,NB_STOK_N
188 . PCT,PCT1,PCT2, , BID
191 my_real,
DIMENSION(:),
ALLOCATABLE,
TARGET :: xe,ve,mse
192 my_real,
DIMENSION(:),
POINTER :: ptr_x,ptr_v,ptr_ms
193 INTEGER,
DIMENSION(:),
POINTER :: PTR_KINET
203 IF (debug(3)>=1.AND.ncycle==0)
THEN
205 nb_jlt_new(jtask) = 0
212 IF( multi_fvm%IS_INT18_LAW151 )
THEN
213 CALL int18_law151_nsv_shift(
'+',jtask-1,nthread,multi_fvm,ipari,interfaces%INTBUF_TAB,npari,ninter,numnod)
227 lindmax = ipari(18,n)*ipari(23,n)
229 IF(int18add(ninter+1) >= 0)
THEN
230 int18add(n+1) = int18add(n)
235 IF(nty == 7.OR.nty == 11.OR.nty == 21.OR.nty == 24.OR.nty == 25)
236 . isens = ipari(64,n)
242 interefric = h3d_data%N_CSE_FRIC_INTER (n)
246 IF (ipari(71,n)>0) nth = ipari(71,n)
249 IF((nty == 7.AND.tt>=ts).OR.nty == 17)
THEN
251 IF(ipari(22,n)==7)
THEN
254 IF(int18add(ninter+1) >= 0.and.jtask == 1
255 . .and.ipari(34,n)==-2)
THEN
256 int18add(n+1) = int18add(n) + lindmax
258 ELSEIF(ipari(44,n)/=0)
THEN
264 nn1 = ninter+nrwall+nrbody+nsect+njoint
265 . +nvolu+nrbag+nfxbody+1
266 nn2 = ninter+nrwall+nrbody+nsect+njoint+1
268 IF(interfaces%INTBUF_TAB(n)%S_NIGE/=0)
THEN
269 ALLOCATE(xe(3*(numnod+interfaces%INTBUF_TAB(n)%S_NIGE)))
270 xe(1:3*numnod) = x(1:3*numnod)
271 xe(3*numnod+1:3*(numnod+interfaces%INTBUF_TAB(n)%S_NIGE)) =
272 * interfaces%INTBUF_TAB(n)%XIGE(1:3*interfaces%INTBUF_TAB(n)%S_NIGE)
273 ALLOCATE(ve(3*(numnod+interfaces%INTBUF_TAB(n)%S_NIGE)))
274 ve(1:3*numnod) = v(1:3*numnod)
275 ve(3*numnod+1:3*(numnod+interfaces%INTBUF_TAB(n)%S_NIGE)) =
276 * interfaces%INTBUF_TAB(n)%VIGE(1:3*interfaces%INTBUF_TAB(n)%S_NIGE)
277 ALLOCATE(mse(numnod+interfaces%INTBUF_TAB(n)%S_NIGE))
278 mse(1:numnod) = ms(1:numnod)
279 mse(numnod+1:numnod+interfaces%INTBUF_TAB(n)%S_NIGE) =
280 * interfaces%INTBUF_TAB(n)%MASSIGE(1:interfaces%INTBUF_TAB(n)%S_NIGE)
284 ptr_kinet => kinet(1:numnod)
285 ELSEIF (multi_fvm%IS_USED .AND. ipari(22,n) == 7)
THEN
286 ptr_x => multi_fvm%X_APPEND
287 ptr_v => multi_fvm%V_APPEND
288 ptr_ms => multi_fvm%MASS_APPEND
289 ptr_kinet => multi_fvm%KINET_APPEND
295 ptr_kinet => kinet(1:numnod)
300 1 ipari ,ptr_x ,a ,ale_connectivity,xcell ,
301 2 icodt ,fsav(1,nth) ,ptr_v ,ptr_ms ,dt2t ,
302 3 neltst ,ityptst ,itab ,stifn ,fskyi ,
303 4 isky ,fcont ,n ,lindmax ,ptr_kinet ,
304 5 jtask ,nb_jlt(jtask ) ,nb_jlt_new(jtask),nb_stok_n(jtask),elbuf_tab ,
305 6 niskyfi(n) ,newfront(n) ,nstrf ,secfcum ,igroups ,
306 7 icontact ,viscn ,idum ,
307 9 idum ,idum ,idum ,fsav(1,nn1) ,nrtmdim ,
309 b eminx ,ixs ,ixs16 ,ixs20 ,fncont ,
310 c ftcont ,iad_elem ,fr_elem ,rcontact ,acontact ,
311 d pcontact ,temp ,fthe ,ftheskyi ,
312 e pm ,iparg ,iad17 ,mskyi_sms ,iskyi_sms ,
313 f nodnx_sms ,ms0 ,qfricint ,npc ,tf ,
314 g condn ,condnskyi ,interfaces%INTBUF_TAB(n),nodadt_therm ,theaccfact ,
315 h fbsav6 ,isensint(1,n) ,sfbsav6 ,ixig3d ,
316 i kxig3d ,wige ,knot ,igeo ,multi_fvm ,
317 j h3d_data ,interfaces%INTBUF_FRIC_TAB,knotlocpc ,knotlocel ,jtask ,
318 k tagncont ,kloadpinter ,loadpinter ,loadp_hyd_inter ,dgaploadint ,
319 l s_loadpinter ,interefric ,xcell_remote(n)%SIZE_MY_REAL_ARRAY_1D,xcell_remote(n)%MY_REAL_ARRAY_1D)
322 1 ipari ,ptr_x ,a ,ale_connectivity,xcell ,
323 2 icodt ,fsav(1,nth) ,ptr_v ,ptr_ms ,dt2t ,
324 3 neltst ,ityptst ,itab ,stifn ,fskyi ,
325 4 isky ,fcont ,n ,lindmax ,ptr_kinet ,
326 5 jtask ,nb_jlt(jtask) ,nb_jlt_new(jtask) ,nb_stok_n(jtask),elbuf_tab ,
327 6 niskyfi(n) ,newfront(n) ,nstrf ,secfcum ,igroups ,
328 7 icontact ,viscn ,num_imp(n) ,
329 9 ns_imp(iad),ne_imp(iad) ,ind_imp(iad) ,fsav(1,nn1) ,nrtmdim ,
331 b eminx ,ixs ,ixs16 ,ixs20 ,fncont ,
332 c ftcont ,iad_elem ,fr_elem ,rcontact ,acontact ,
333 d pcontact ,temp ,fthe ,ftheskyi ,
334 e pm ,iparg ,iad17 ,mskyi_sms ,iskyi_sms ,
335 f nodnx_sms ,ms0 ,qfricint ,npc ,tf ,
336 g condn ,condnskyi ,interfaces%INTBUF_TAB(n),nodadt_therm,theaccfact ,
337 h fbsav6 ,isensint(1,n) ,sfbsav6 ,ixig3d ,
338 i kxig3d ,wige ,knot ,igeo ,multi_fvm,
339 j h3d_data ,interfaces%INTBUF_FRIC_TAB,knotlocpc ,knotlocel ,jtask ,
340 k tagncont ,kloadpinter ,loadpinter ,loadp_hyd_inter ,dgaploadint ,
341 l s_loadpinter,interefric ,xcell_remote(n)%SIZE_MY_REAL_ARRAY_1D,xcell_remote(n)%MY_REAL_ARRAY_1D)
344 IF(nty == 17 .AND. ipari(33,n) == 0)
THEN
347 nmes =igrbric(ign)%NENTITY
348 nme =igrbric(ige)%NENTITY
349 iad17 = iad17+6*(nme+nmes)
351 IF(interfaces%INTBUF_TAB(n)%S_NIGE/=0)
THEN
352 DEALLOCATE(xe,ve,mse)
361 2 icodt ,fsav(1,n) ,v ,ms ,dt2t ,
362 3 neltst ,ityptst ,itab ,stifn ,fskyi ,
363 4 isky ,fcont ,lindmax ,
365 6 niskyfi(n),nstrf , secfcum ,viscn ,n ,
366 . fsav(1,ninter+nrwall+nrbody+nsect+njoint+nvolu+nrbag+nfxbody) ,
367 9 idum ,idum ,idum ,idum ,fncont ,
368 a ftcont ,mskyi_sms ,iskyi_sms ,nodnx_sms ,
369 b icontact ,interfaces%INTBUF_TAB(n),fbsav6,isensint(1,n),sfbsav6 ,
370 c h3d_data ,nodadt_therm)
374 2 icodt ,fsav(1,n) ,v ,ms ,dt2t ,
375 3 neltst ,ityptst ,itab ,stifn ,fskyi ,
376 4 isky ,fcont ,lindmax ,
377 5 jtask,nb_jlt(jtask),nb_jlt_new(jtask),nb_stok_n(jtask),
378 6 niskyfi(n),nstrf , secfcum ,viscn ,n ,
379 . fsav(1,ninter+nrwall+nrbody+nsect+njoint+nvolu+nrbag+nfxbody) ,
380 9 num_imp(n),ns_imp(iad) ,ne_imp(iad) ,ind_imp(iad),fncont ,
381 a ftcont ,mskyi_sms ,iskyi_sms ,nodnx_sms ,
382 b icontact ,interfaces%INTBUF_TAB(n),fbsav6,isensint(1,n),sfbsav6 ,
383 c h3d_data ,nodadt_therm)
388 ELSEIF(nty==11.AND.tt>=ts)
THEN
391 nn1 = ninter+nrwall+nrbody+nsect+njoint
392 . +nvolu+nrbag+nfxbody+1
397 2 icodt ,fsav(1,nth) ,v ,ms ,dt2t ,
398 3 neltst ,ityptst ,itab ,stifn ,fskyi ,
399 4 isky ,fcont ,n ,lindmax ,jtask ,
400 5 nb_jlt(jtask),nb_jlt_new(jtask),nb_stok_n(jtask) ,niskyfi(n) ,
401 6 newfront(n),nstrf,secfcum, viscn ,idum ,idum ,
402 7 idum ,mskyi_sms ,iskyi_sms ,nodnx_sms,icontact ,
403 8 interfaces%INTBUF_TAB(n),pm ,temp , fthe ,ftheskyi,
404 9 npc , tf ,condn ,condnskyi ,fbsav6 ,
405 a isensint(1,n),sfbsav6 ,fsav(1,nn1) ,h3d_data ,interfaces%INTBUF_FRIC_TAB,
406 b jtask ,tagncont ,kloadpinter ,loadpinter,loadp_hyd_inter,
407 c dgaploadint,s_loadpinter ,nodadt_therm)
411 2 icodt ,fsav(1,nth) ,v ,ms ,dt2t ,
412 3 neltst ,ityptst ,itab ,stifn ,fskyi ,
413 4 isky ,fcont ,n ,lindmax ,jtask ,
414 5 nb_jlt(jtask),nb_jlt_new(jtask),nb_stok_n(jtask) ,niskyfi(n) ,
415 6 newfront(n),nstrf,secfcum, viscn ,num_imp(n),ns_imp(iad),
416 7 ne_imp(iad),mskyi_sms ,iskyi_sms ,nodnx_sms,icontact ,
417 8 interfaces%INTBUF_TAB(n),pm ,temp , fthe ,ftheskyi ,
418 9 npc , tf ,condn ,condnskyi ,fbsav6 ,
419 a isensint(1,n),sfbsav6 ,fsav(1,nn1) ,h3d_data ,interfaces%INTBUF_FRIC_TAB,
420 b jtask ,tagncont ,kloadpinter ,loadpinter,loadp_hyd_inter,
421 c dgaploadint,s_loadpinter ,nodadt_therm)
424 ELSEIF(nty == 20)
THEN
426 IF(ipari(22,n)==7)
THEN
429 ELSEIF(ipari(44,n)/=0)
THEN
435 nn1 = ninter+nrwall+nrbody+nsect+njoint
436 . +nvolu+nrbag+nfxbody+1
437 nn2 = ninter+nrwall+nrbody+nsect+njoint+1
439 IF(impl_s == 0) iad=1
443 2 icodt ,fsav(1,n) ,v ,ms ,dt2t ,
444 3 neltst ,ityptst ,itab ,stifn ,fskyi ,
445 4 isky ,fcont ,n ,lindmax ,kinet ,
446 5 jtask ,nb_jlt(jtask),nb_jlt_new(jtask) ,nb_stok_n(jtask),
447 6 niskyfi(n),newfront(n) , nstrf ,secfcum ,icontact ,
448 7 viscn ,num_imp(n) ,
449 9 ns_imp(iad),ne_imp(iad) ,ind_imp(iad) ,fsav(1,nn1),nrtmdim ,
451 b eminx ,ixs ,ixs16 ,ixs20 ,fncont ,
452 c ftcont ,iad_elem ,fr_elem ,rcontact ,acontact ,
453 d pcontact ,temp ,fthe ,ftheskyi ,
454 e pm ,iparg ,iad17 ,weight ,niskyfie(n) ,
455 f irlen20 ,islen20 ,irlen20t ,islen20t ,irlen20e ,
456 g islen20e ,mskyi_sms ,iskyi_sms ,nodnx_sms ,npc ,
457 h tf ,interfaces%INTBUF_TAB(n),fbsav6,isensint(1,n) ,sfbsav6 ,
458 i h3d_data ,theaccfact)
460 IF(impl_s == 1) iad=iad+num_imp(n)
463 ELSEIF(nty == 21.AND.tt>=ts)
THEN
465 IF(ipari(44,n)/=0)
THEN
473 1 ipari ,interfaces%INTBUF_TAB(n) ,x ,a ,
474 2 icodt ,fsav(1,n) ,v ,ms ,itab ,
475 3 stifn ,fskyi ,isky ,fcontg ,n ,
477 + jtask ,nb_jlt(jtask),nb_jlt_new(jtask),nb_stok_n(jtask),
478 5 nstrf ,secfcum ,icontact ,viscn ,idum ,
479 6 idum ,idum ,idum ,nrtmdim ,fncontg ,
480 7 ftcontg ,rcontact ,acontact ,pcontact ,intstamp(istamp),
481 8 weight ,temp ,fthe ,ftheskyi ,mskyi_sms ,
482 9 iskyi_sms ,nodnx_sms ,nodglob ,npc ,tf ,
483 a qfricint ,ncont ,indexcont ,tagcont ,condn ,
484 b condnskyi ,dt2t ,neltst ,ityptst ,kinet ,
485 c fbsav6 ,isensint(1,n),sfbsav6 ,niskyfi(n) ,h3d_data ,
486 d pskids ,tagncont ,kloadpinter,loadpinter ,loadp_hyd_inter,
487 e dgaploadint,s_loadpinter,interefric ,nodadt_therm,theaccfact ,
488 f interfaces%INTBUF_FRIC_TAB)
492 ELSEIF(nty == 22)
THEN
496 nn1 = ninter+nrwall+nrbody+nsect+njoint
497 . +nvolu+nrbag+nfxbody+1
498 nn2 = ninter+nrwall+nrbody+nsect+njoint+1
500 1 ipari ,interfaces%INTBUF_TAB(n) ,x ,a ,
501 2 icodt ,fsav(1,n) ,v ,ms ,dt2t ,
502 3 neltst ,ityptst ,itab ,stifn ,fskyi
503 4 isky ,fcont ,n ,lindmax ,kinet ,
504 5 jtask ,nb_jlt(jtask) ,nb_jlt_new(jtask) ,nb_stok_n(jtask) ,
505 6 niskyfi(n) ,newfront(n) ,nstrf ,secfcum ,icontact ,
507 9 idum ,idum ,idum ,fsav(1,nn1) ,nrtmdim ,
508 a fsav(1,nn2) ,igrbric ,
509 b eminx ,ixs ,ixs16 ,ixs20 ,fncont ,
510 c ftcont ,iad_elem ,fr_elem ,rcontact ,acontact ,
511 d pcontact ,temp ,fthe ,ftheskyi ,
512 e pm ,iparg ,iad17 ,mskyi_sms ,iskyi_sms ,
513 f nodnx_sms ,ms0 ,elbuf_tab ,nv46 ,fbsav6 ,
514 g isensint(1,n),sfbsav6 ,h3d_data )
517 ELSEIF(nty == 23)
THEN
519 IF(ipari(44,n)/=0)
THEN
526 1 ipari ,interfaces%INTBUF_TAB(n) ,x ,a ,
527 2 icodt ,fsav(1,n) ,v ,ms ,itab ,
528 3 stifn ,fskyi ,isky ,fcont ,n ,
530 + jtask ,nb_jlt(jtask),nb_jlt_new(jtask),nb_stok_n(jtask),
531 5 nstrf ,secfcum ,icontact ,viscn ,idum ,
532 6 idum ,idum ,idum ,nrtmdim ,fncont ,
533 7 ftcont ,rcontact ,acontact ,pcontact ,kinet ,
534 8 weight ,mskyi_sms ,iskyi_sms ,nodnx_sms ,nodglob ,
535 9 npc ,tf ,niskyfi(n) ,newfront(n),wa ,
536 a fbsav6,isensint(1,n) ,sfbsav6 ,dt2t ,h3d_data )
538 ELSEIF(nty == 24.AND.tt>=ts)
THEN
540 IF(ipari(44,n)/=0)
THEN
546 nn1 = ninter+nrwall+nrbody+nsect+njoint
547 . +nvolu+nrbag+nfxbody+1
548 nn2 = ninter+nrwall+nrbody+nsect+njoint+1
551 1 ipari ,interfaces%INTBUF_TAB(n),x ,a ,
552 2 icodt ,fsav(1,n) ,v ,ms ,dt2t ,
553 3 neltst ,ityptst ,itab ,stifn ,fskyi ,
554 4 isky ,fcont ,n ,lindmax ,kinet ,
555 5 jtask ,nb_jlt(jtask),nb_jlt_new(jtask),nb_stok_n(jtask),
556 6 niskyfi(n),newfront(n) , nstrf ,secfcum ,icontact ,
558 9 idum ,idum ,idum ,fsav(1,nn1) ,nrtmdim ,
560 b eminx ,ixs ,ixs16 ,ixs20 ,fncont ,
561 c ftcont ,iad_elem ,fr_elem ,rcontact ,acontact ,
562 d pcontact ,temp ,fthe ,ftheskyi ,
563 e pm ,iparg ,iad17 ,mskyi_sms ,iskyi_sms ,
564 f nodnx_sms ,ms0 ,inod_pxfem ,ms_ply ,wagap ,
565 g fbsav6 ,isensint(1,n),
566 h sfbsav6 ,h3d_data ,interfaces%INTBUF_FRIC_TAB ,t2main_sms,forneqs ,
567 i t2fac_sms ,npc ,tf ,tagncont ,kloadpinter ,
568 j loadpinter ,loadp_hyd_inter,dgaploadint,s_loadpinter,interefric,
569 k interfaces,nisubmax)
573 1 ipari ,interfaces%INTBUF_TAB(n),x ,a ,
574 2 icodt ,fsav(1,n) ,v ,ms ,dt2t ,
575 3 neltst ,ityptst ,itab ,stifn ,fskyi ,
576 4 isky ,fcont ,n ,lindmax ,kinet ,
577 5 jtask ,nb_jlt(jtask),nb_jlt_new(jtask) ,nb_stok_n(jtask),
578 6 niskyfi(n),newfront(n) , nstrf ,secfcum ,icontact ,
579 7 viscn ,num_imp(n) ,
580 9 ns_imp(iad),ne_imp(iad) ,ind_imp(iad) ,fsav(1,nn1),nrtmdim ,
582 b eminx ,ixs ,ixs16 ,ixs20 ,fncont ,
583 c ftcont ,iad_elem ,fr_elem ,rcontact ,acontact ,
584 d pcontact ,temp ,fthe ,ftheskyi ,
585 e pm ,iparg ,iad17 ,mskyi_sms ,iskyi_sms ,
586 f nodnx_sms ,ms0 ,inod_pxfem ,ms_ply ,wagap ,
587 g fbsav6 ,isensint(1,n),
588 h sfbsav6 ,h3d_data ,interfaces%INTBUF_FRIC_TAB ,t2main_sms,forneqs ,
589 i t2fac_sms ,npc ,tf ,tagncont ,kloadpinter ,
590 j loadpinter ,loadp_hyd_inter,dgaploadint,s_loadpinter,interefric,
591 k interfaces,nisubmax)
595 ELSEIF(nty == 25.AND.tt>=ts)
THEN
597 IF(ipari(44,n)/=0)
THEN
603 nn1 = ninter+nrwall+nrbody+nsect+njoint
604 . +nvolu+nrbag+nfxbody
605 nn2 = ninter+nrwall+nrbody+nsect+njoint+1
607 lindmax = interfaces%INTBUF_TAB(n)%I_STOK(2)
608 IF(ipari(58,n)/=0) lindmax=
max(lindmax,
609 . interfaces%INTBUF_TAB(n)%I_STOK_E(1),interfaces%INTBUF_TAB(n
611 1 ipari ,interfaces%INTBUF_TAB(n),x ,a ,
612 2 icodt ,fsav(1,n) ,v ,ms ,dt2t ,
613 3 neltst ,ityptst ,itab ,stifn ,fskyi ,
614 4 isky ,fcont ,n ,lindmax ,kinet ,
615 5 jtask ,nb25_impct(jtask),
616 6 niskyfi(n),newfront(n) , nstrf ,secfcum ,icontact ,
618 9 idum ,idum ,idum ,fsav(1,nn1) ,nrtmdim ,
620 b eminx ,ixs ,ixs16 ,ixs20 ,fncont ,
621 c ftcont ,iad_elem ,fr_elem ,rcontact ,acontact ,
622 d pcontact ,temp ,fthe ,ftheskyi ,
623 e pm ,iparg ,iad17 ,mskyi_sms ,iskyi_sms ,
624 f nodnx_sms ,ms0 ,inod_pxfem ,ms_ply ,wagap ,
625 g fbsav6 ,isensint(1,n),nodadt_therm,theaccfact ,
626 h sfbsav6 ,h3d_data ,interfaces%INTBUF_FRIC_TAB,niskyfie(n),
627 i apinch ,stifpinch ,npc ,tf , condn ,
628 j condnskyi ,qfricint ,tagncont ,kloadpinter,loadpinter ,
629 k loadp_hyd_inter,dgaploadint
631 IF(impl_s == 1) iad=iad+num_imp(n)
646 IF( multi_fvm%IS_INT18_LAW151 )
THEN
648 CALL int18_law151_nsv_shift(
'-',jtask-1,nthread,multi_fvm,ipari,interfaces%INTBUF_TAB,npari,ninter,numnod)
651 IF (debug(3)>=1)
THEN
652 IF(mod(ncycle+1,debug(3))==0)
THEN
653 IF (nb_jlt(jtask)==0)
THEN
657 pct = hundred*nb_jlt_new(jtask)/nb_jlt(jtask)
658 pct1 = hundred - hundred*nb_stok_n(jtask)/nb_jlt(jtask)
662 .
'(A,I6,A,I4,A,I4,A,I10,A,I10,2X,F5.2,A,A,I10,4X,F5.2,A)')
663 .
' NCYCLE = ',ncycle,
664 .
' NSPMD = ',ispmd+1,
666 .
' CANDIDATS = ',nb_jlt(jtask),
667 .
' OPT CAND = ',nb_stok_n(jtask),pct1,
'%',
668 .
' IMPACTS = ',nb_jlt_new(jtask),pct,
'%'
669#include "lockoff.inc"
670 nb_jlt_new(jtask) = 0
677 IF(ninter25 /= 0)
THEN
678 IF(mod(ncycle+1,debug(3))==0)
THEN
679 IF(nb25_candt(1)==0)
THEN
682 pct2 = hundred * nb25_dst2(jtask) / nb25_candt(1)
686 .
'(A,I6,A,I4,A,I4,A,I10,A,I10,2X,F6.2,A,A,I10,A,I10)')
687 .
' NCYCLE = ',ncycle,
688 .
' NSPMD = ',ispmd+1,
690 .
' CANDIDATS = ',nb25_candt(jtask),
691 .
' OPT CAND VS NEW IMPACTS = ',nb25_dst2(jtask),pct2,
'%',
692 .
' CAND VS OLD IMPACTS = ',nb25_dst1(jtask),
693 .
' IMPACTS = ',nb25_impct(jtask)
694#include "lockoff.inc"
698 nb25_impct(jtask) = 0
699 nb25_candt(jtask) = 0