OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
rdcomm.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!|| rdcomi ../engine/source/output/restart/rdcomm.F
25!||--- called by ------------------------------------------------------
26!|| rdresa ../engine/source/output/restart/rdresa.F
27!||--- calls -----------------------------------------------------
28!|| read_array ../common_source/comm/read_array.F
29!||--- uses -----------------------------------------------------
30!|| ale_mod ../common_source/modules/ale/ale_mod.F
31!|| alefvm_mod ../common_source/modules/ale/alefvm_mod.F
32!|| alemuscl_mod ../common_source/modules/ale/alemuscl_mod.F
33!|| damp_mod ../engine/share/modules/damp_mod.F
34!|| drape_mod ../engine/share/modules/drape_mod.F
35!|| eos_param_mod ../common_source/modules/mat_elem/eos_param_mod.F90
36!|| glob_therm_mod ../common_source/modules/mat_elem/glob_therm_mod.F90
37!|| h3d_mod ../engine/share/modules/h3d_mod.F
38!|| intbufdef_mod ../common_source/modules/interfaces/intbufdef_mod.F90
39!|| loads_mod ../common_source/modules/loads/loads_mod.F90
40!|| mat_elem_mod ../common_source/modules/mat_elem/mat_elem_mod.F90
41!|| multi_fvm_mod ../common_source/modules/ale/multi_fvm_mod.F90
42!|| outmax_mod ../common_source/modules/outmax_mod.F
43!|| output_mod ../common_source/modules/output/output_mod.F90
44!|| parameters_mod ../common_source/modules/interfaces/parameters_mod.F
45!|| pblast_mod ../common_source/modules/loads/pblast_mod.F90
46!|| pinchtype_mod ../common_source/modules/pinchtype_mod.F
47!|| pload_cyl_mod ../common_source/modules/loads/pload_cyl_mod.F
48!|| rbe3_mod ../common_source/modules/constraints/rbe3_mod.F90
49!|| rwall_mod ../common_source/modules/constraints/rwall_mod.F90
50!|| sensor_mod ../common_source/modules/sensor_mod.F90
51!|| skew_mod ../common_source/modules/skew_mod.F90
52!|| state_mod ../common_source/modules/state_mod.F
53!|| th_mod ../engine/share/modules/th_mod.F
54!||====================================================================
55 SUBROUTINE rdcomi(OUTPUT,TABVINT,MULTI_FVM,H3D_DATA,DYNAIN_DATA,
56 . INTER_PARAMETERS,SENSORS,LOADS ,MAT_ELEM,SKEWS,
57 . GLOB_THERM,PBLAST,RBE3 ,RWALL)
58C-----------------------------------------------
59C M o d u l e s
60C-----------------------------------------------
61 USE multi_fvm_mod
62 USE h3d_mod
63 USE pinchtype_mod
64 USE outmax_mod
65 USE pload_cyl_mod
66 USE drape_mod
67 USE state_mod
68 USE intbufdef_mod
69 USE alemuscl_mod , only:alemuscl_param
70 USE alefvm_mod , only:alefvm_param
71 USE ale_mod
72 USE parameters_mod
73 USE sensor_mod
74 USE loads_mod
75 USE mat_elem_mod
76 USE damp_mod
77 USE skew_mod
78 USE glob_therm_mod
79 USE pblast_mod
80 USE rbe3_mod
81 USE th_mod , ONLY : th_has_noda_pext
82 USE output_mod
83 USE eos_param_mod , ONLY : analy_temp
84 USE rwall_mod
85C-----------------------------------------------
86C I m p l i c i t T y p e s
87C-----------------------------------------------
88#include "implicit_f.inc"
89 INTEGER TABVINT(*)
90 TYPE(OUTPUT_), INTENT(INOUT) :: OUTPUT
91 TYPE(multi_fvm_struct) :: MULTI_FVM
92 TYPE(h3d_database) :: H3D_DATA
93 TYPE (DYNAIN_DATABASE),INTENT(INOUT) :: DYNAIN_DATA
94 TYPE (PARAMETERS_) ,INTENT(INOUT) :: INTER_PARAMETERS
95 TYPE (SENSORS_) ,INTENT(INOUT) :: SENSORS
96 TYPE (LOADS_) ,INTENT(INOUT) :: LOADS
97 TYPE (MAT_ELEM_) ,INTENT(INOUT) :: MAT_ELEM
98 TYPE (SKEW_) ,INTENT(INOUT) :: SKEWS
99 TYPE (GLOB_THERM_) ,INTENT(INOUT) :: GLOB_THERM
100 TYPE (PBLAST_) ,INTENT(INOUT) :: PBLAST
101 TYPE (RBE3_) ,INTENT(INOUT) :: RBE3
102 TYPE (RWALL_) ,INTENT(INOUT) :: RWALL
103C-----------------------------------------------
104C C o m m o n B l o c k s
105C-----------------------------------------------
106#include "lagmult.inc"
107#include "r4r8_p.inc"
108#include "com01_c.inc"
109#include "com04_c.inc"
110#include "com09_c.inc"
111#include "com10_c.inc"
112#include "sphcom.inc"
113#include "param_c.inc"
114#include "scr03_c.inc"
115#include "scr05_c.inc"
116#include "scr06_c.inc"
117#include "scr07_c.inc"
118#include "scr14_c.inc"
119#include "scr16_c.inc"
120#include "cong1_c.inc"
121#include "task_c.inc"
122#include "scrnoi_c.inc"
123#include "scrfs_c.inc"
124#include "scrcut_c.inc"
125#include "warn_c.inc"
126#include "couple_c.inc"
127#include "scr18_c.inc"
128#include "scr02_c.inc"
129#include "check.inc"
130#include "stati_c.inc"
131#include "timeri_c.inc"
132#include "scr17_c.inc"
133#include "rad2r_c.inc"
134#include "spmd_c.inc"
135#include "fxbcom.inc"
136#include "eigcom.inc"
137#include "chara_c.inc"
138#include "flowcom.inc"
139#include "impl1_c.inc"
140#include "tabsiz_c.inc"
141#include "remesh_c.inc"
142#include "com_xfem1.inc"
143#include "sysunit.inc"
144#include "intstamp_c.inc"
145#include "sms_c.inc"
146#include "submodel.inc"
147#include "inter18.inc"
148#include "inter22.inc"
149#include "drape_c.inc"
150#include "ige3d_c.inc"
151C----------------------------------------------
152C L o c a l V a r i a b l e s
153C-----------------------------------------------
154 INTEGER LWASPH44C, I, NEXT,IDUM1,ISUSED,SENS_NSTOP,SENS_NSTAT,SENS_NOUTP,SENS_NANIM
155C-----------------------------------------------
156C reading of integer variables.
157C-----------------------------------------------
158 n2d =tabvint(1)
159 ncpri =tabvint(2)
160 !IUSERWI =TABVINT(3)
161 iale =tabvint(4)
162 ngroup =tabvint(5)
163 mat_elem%NGROUP = ngroup
164 ncycle =tabvint(6)
165 irun =tabvint(7)
166 iger =tabvint(8)
167 lbufel =tabvint(9)
168 iroddl =tabvint(10)
169 ieuler =tabvint(11)
170 glob_therm%ITHERM =tabvint(12)
171 ihsh =tabvint(13)
172 itestv =tabvint(14)
173 iturb =tabvint(15)
174 ilag =tabvint(16)
175 isecut =tabvint(17)
176 idamp =tabvint(18)
177 irxdp =tabvint(19)
178 nmult =tabvint(20)
179 integ8 =tabvint(21)
180 isigi =tabvint(22)
181 nspmd =tabvint(23)
182 lenwa =tabvint(24)
183 analy_temp = tabvint(25)
184 th_has_noda_pext = tabvint(26)
185 nnods =tabvint(27)
186 ncnois =tabvint(28)
187 lcne0 =tabvint(29)
188 ipari0 =tabvint(30)
189 imaximp =tabvint(31)
190 nnoiser =tabvint(32)
191 nspgroup=tabvint(33)
192 iresmd =tabvint(34)
193 ifrwv =tabvint(35)
194C-----
195 itestv = iabs(itestv)
196 IF (ir4r8==2) itestv = -itestv
197C ITESTV = IABS(ITESTV)
198C-----
199C -------- Size of integer arrays --------
200 sicode = tabvint(37)
201 ! = TABVINT(38)
202 siskew = tabvint(39)
203 siskwn = tabvint(40)
204 siframe = tabvint(41)
205 sneth = tabvint(42)
206 sibcslag = tabvint(43)
207 sipart = tabvint(44)
208 snom_opt = tabvint(45)
209 snpc = tabvint(46)
210C Empty field = TABVINT(47)
211 sixtg = tabvint(48)
212 sixtg1 = tabvint(49)
213 sixs = tabvint(50)
214 sixq = tabvint(51)
215 sixc = tabvint(52)
216 sixt = tabvint(53)
217 sixp = tabvint(54)
218 sixr = tabvint(55)
219 sitab = tabvint(56)
220 sitabm1 = tabvint(57)
221 sgjbufi = tabvint(58)
222 ale%GLOBAL%SNALE = tabvint(59)
223 ale%GLOBAL%SIELVS = tabvint(60)
224 sifill = tabvint(62)
225 sims = tabvint(63)
226!! to empty+++
227 sisubs = tabvint(64)
228 sigrn = tabvint(65)
229 sisurf = tabvint(66)
230 sislin = tabvint(67)
231 sibufssg = tabvint(68)
232!! to empty---
233 skxx = tabvint(69)
234 sixx = tabvint(70)
235 skxsp = tabvint(71)
236 sixsp = tabvint(72)
237 snod2sp = tabvint(73)
238 sispsym = tabvint(74)
239 sispcond = tabvint(75)
240 sisphio = tabvint(76)
241 slprtsph = tabvint(77)
242 slonfsph = tabvint(78)
243 sibcl = tabvint(79)
244 sibfv = tabvint(80)
245 siecran = tabvint(81)
246 silas = tabvint(82)
247 siebcs = tabvint(83)
248 slaccelm = tabvint(84)
249
250 snnlink = tabvint(86)
251 slnlink = tabvint(87)
252 !SIUSER = TABVINT(88) - unused
253 siparg = tabvint(89)
254 sdd_iad = tabvint(90)
255 sigrv = tabvint(91)
256 slgrav = tabvint(92)
257 siexlnk = tabvint(93)
258 sibvel = tabvint(94)
259 slbvel = tabvint(95)
260 siactiv = tabvint(96)
261 sibmpc = tabvint(97)
262 skinet = tabvint(98)
263 sipari = tabvint(99)
264 sinbuf = 0 !TABVINT(100)
265 snprw = tabvint(101)
266 ! TABVINT(102) obsolete option removed
267 siconx = tabvint(103)
268 smadprt = tabvint(104)
269 smadsh4 = tabvint(105)
270 smadsh3 = tabvint(106)
271 smadsol = tabvint(107)
272 smadnod = tabvint(108)
273 smadfail = tabvint(109)
274 snpby = tabvint(110)
275 snpbyl = tabvint(111)
276 slpby = tabvint(112)
277 slpbyl = tabvint(113)
278 slrivet = tabvint(114)
279 snstrf = tabvint(115)
280 sljoint = tabvint(116)
281 snodpor = tabvint(117)
282 smonvol = tabvint(118)
283 slagbuf = tabvint(119)
284! SITHGRP = TABVINT(120) - obsolete
285 siad_elem = tabvint(121)
286 siad_rby = tabvint(122)
287 siad_rby2 = tabvint(123)
288 siad_i2m = tabvint(124)
289 siadcj = tabvint(125)
290 siad_rbm = tabvint(126)
291 siad_rbm2 = tabvint(127)
292 siad_sec = tabvint(128)
293 siad_cut = tabvint(129)
294 sfr_elem = tabvint(130)
295 sfr_rby = tabvint(131)
296 sfr_wall = tabvint(132)
297 sfr_rby2 = tabvint(133)
298 sfr_i2m = tabvint(134)
299 sfr_mv = tabvint(135)
300 sfr_ll = tabvint(136)
301 sfr_cj = tabvint(137)
302 sfr_rbm = tabvint(138)
303 sfr_rbm2 = tabvint(139)
304 sfr_sec = tabvint(140)
305 sfr_cut = tabvint(141)
306 srg_cut = tabvint(142)
307 sfr_mad = tabvint(143)
308 sfr_i18 = tabvint(144)
309 sweight = tabvint(145)
310 snewfront = tabvint(146)
311 snodglob = tabvint(147)
312 snbrcvois = tabvint(148)
313 snbsdvois = tabvint(149)
314 slnrcvois = tabvint(150) ! NRCVOIS
315 slnsdvois = tabvint(151)
316 snercvois = tabvint(152)
317 snesdvois = tabvint(153)
318 slercvois = tabvint(154)
319 slesdvois = tabvint(155)
320 !TABVINT(156) obsolete option removed
321 !TABVINT(157) obsolete option removed
322 snpsegcom = tabvint(158)
323 slsegcom = tabvint(159)
324 snporgeo = tabvint(160)
325 slnodpor = tabvint(161)
326C-----
327C siz of real arrays
328C
329C TABVINT(162) = 3*(NUMNOD_L+NRCVVOIS_L)
330C TABVINT(163) = 3*(NUMNOD_L+NRCVVOIS_L)
331C TABVINT(164) = 3*(NUMNOD_L+NRCVVOIS_L)
332C TABVINT(165) = 3*NUMNOD_L*IRODDL
333c IF(ISECUT > 0 .OR. IISROT > 0 .OR. IMPOSE_DR /= 0 .OR. IDROT > 0)THEN
334c TABVINT(166) = 3*NUMNOD_L*IRODDL
335c ELSE
336c TABVINT(166) = 0
337c ENDIF
338
339 sx = tabvint(162)
340 sd = tabvint(163)
341 sv = tabvint(164)
342 svr = tabvint(165)
343 sdr = tabvint(166)
344 sthke = tabvint(167)
345 sdampr = tabvint(168)
346 sdamp = tabvint(169)
347 sms = tabvint(170) ! NUMNOD
348 sin = tabvint(171) ! NUMNOD * IRODDL
349 stf = tabvint(172)
350 spm = tabvint(173)
351 sskew = tabvint(174)
352 sxframe = tabvint(175)
353 sgeo = tabvint(176)
354 seani = tabvint(177)
355 sbufmat = tabvint(178)
356 sbufgeo = tabvint(179)
357 sbufsf = tabvint(180)
358 srbmpc = tabvint(181)
359 sgjbufr = tabvint(182)
360 sw = tabvint(183) ! IALE > 0 => 3 * NUMNOD
361 sveul = tabvint(184)
362 sfill = tabvint(185)
363 sdfill = tabvint(186)
364 salph = tabvint(187)
365 swb = tabvint(188) ! 0 or 3 * NUMNOD or 4 * NUMNOD
366 sdsave = tabvint(189) ! 0 or 3 * NUMNOD
367 ! = TABVINT(190)
368 sasave = tabvint(191)
369 ! = TABVINT(192)
370 ! = TABVINT(193)
371 ! = TABVINT(194)
372 ! = TABVINT(195)
373 sspbuf = tabvint(196)
374 svsphio = tabvint(200)
375 ssphveln = tabvint(201)
376 smsnf = tabvint(202)
377 sforc = tabvint(203)
378 svel = tabvint(204)
379 sfsav = tabvint(205)
380 sfzero = tabvint(206) !3*NUMNOD or 12 *(NUMELC+NUMELTG/)
381 sxlas = tabvint(207)
382 srebcs = tabvint(208)
383 saccelm = tabvint(209)
384
385 sfbvel = tabvint(211)
386 sgrav = tabvint(212)
387 sfr_wave = tabvint(213)
388 sparts0 = tabvint(214)
389 !SUSER = TABVINT(215) unused
390 selbuf = tabvint(216)
391 sbufin = 0 !TABVINT(217)
392 srwbuf = tabvint(218)
393 srwsav = tabvint(219)
394 rwall%SZ_RWSAV = srwsav
395 ! TABVINT(220) obsolete option removed
396 srby = tabvint(221)
397 srbyl = tabvint(222)
398 srivet = tabvint(223)
399 ssecbuf = tabvint(224)
400 svolmon = tabvint(225)
401 slambda = tabvint(226)
402 srconx = tabvint(227)
403 sfasolfr = tabvint(228)
404 sdd_r2r = tabvint(229)
405 slprw = tabvint(230)
406 rwall%SZ_LPRW = slprw
407 rwall%nrwlp = nrwlp
408 !TABVINT(231) obsolete option removed
409 siexmad = tabvint(232)
410 sigeo = tabvint(233)
411 sipm = tabvint(234)
412 skippon = tabvint(235)
413 ssegquadfr= tabvint(236)
414C-----1
415 nummat =tabvint(238)
416 mat_elem%NUMMAT = nummat
417 numnod =tabvint(239) ! NUMNOD
418 numskw =tabvint(240)
419 numbcs =tabvint(241)
420 nanaly =tabvint(242)
421 numelq =tabvint(243)
422 numels =tabvint(244)
423 numelc =tabvint(245)
424 numelt =tabvint(246)
425 numgeo =tabvint(247)
426 nfunct =tabvint(248)
427 nconld =tabvint(249)
428 ninvel =tabvint(250)
429 nlaser =tabvint(251)
430 ninter =tabvint(252)
431 nrwall =tabvint(253)
432 rwall%NRWALL = nrwall
433 nrbody =tabvint(254)
434 nodmas =tabvint(255)
435 nfxvel =tabvint(256)
436 nrivet =tabvint(257)
437 numelr =tabvint(258)
438 numelp =tabvint(259)
439 !obsolete option removed : TABVINT(260)=0
440 nsect =tabvint(261)
441 nrbag =tabvint(262)
442 njoint =tabvint(263)
443 ale%GRID%NWALE_RST = tabvint(264)
444 ale%GRID%NWALE = ale%GRID%NWALE_RST
445!------------------------------!
446 numeltg=tabvint(265)
447 nslag =tabvint(266)
448 nfacx =tabvint(267)
449 numpor =tabvint(268)
450 naccelm=tabvint(269)
451 sensors%NSENSOR=tabvint(270)
452 npreten=tabvint(271)
453 nvolu =tabvint(272)
454 npart =tabvint(273)
455 nsurf =tabvint(274)
456 nsubs =tabvint(275)
457 ngrav =tabvint(276)
458 nrby2 =tabvint(277)
459 ngrnod =tabvint(278)
460 ngrbric=tabvint(279)
461 ngrquad=tabvint(280)
462 ngrshel=tabvint(281)
463 ngrsh3n=tabvint(282)
464 ngrbeam=tabvint(283)
465 ngrtrus=tabvint(284)
466 ngrspri=tabvint(285)
467 nlink =tabvint(286)
468 invstr =tabvint(287)
469 nslin =tabvint(288)
470 numelx =tabvint(289)
471 nconx =tabvint(290)
472 isumnx =tabvint(291)
473 nanim1d=tabvint(292)
474 nr2rlnk=tabvint(293)
475C-----
476 output%DATA%ANIM_HAS_NODA_PEXT = max(tabvint(294)/10 , output%DATA%ANIM_HAS_NODA_PEXT) !last digit, unless /ANIM/NODA/PEXT is present
477 output%DATA%H3D_HAS_NODA_PEXT = max(mod(tabvint(294),10), output%DATA%H3D_HAS_NODA_PEXT) !2nd digit, unless /H3D/NODA/PEXT is present
478 ianim =tabvint(295)
479 ioutp =tabvint(296)
480 nhin2 =tabvint(297)
481 ideli7=tabvint(298)
482C-----2
483 nsnod =tabvint(299)
484 nselq =tabvint(300)
485 nsels =tabvint(301)
486 nselc =tabvint(302)
487 nselt =tabvint(303)
488 nselp =tabvint(304)
489 nselr =tabvint(305)
490 nsrby =tabvint(306)
491 nsmat =tabvint(307)
492 nseltg =tabvint(308)
493 nthgrp =tabvint(309)
494 nthgrp0=tabvint(310)
495 nselu =tabvint(311)
496C-----
497 npropm =tabvint(312)
498 nvsiz =tabvint(313)
499 npropg =tabvint(314)
500 nparg =tabvint(315)
501 lveul =tabvint(316)
502 nixfr1 =tabvint(325)
503 nixfr2 =tabvint(326)
504 npari =tabvint(327)
505 lwamp =tabvint(328)
506 lwanmp =tabvint(329)
507 nisx =tabvint(330)
508 ligrn1 =tabvint(331)
509 lisub1 =tabvint(332)
510 nimv =tabvint(333)
511 nnpby =tabvint(334)
512 nrby =tabvint(335)
513 nifv =tabvint(336)
514 nthvki =tabvint(337)
515 npsav =tabvint(338)
516 nrcnx =tabvint(339)
517 nr2r =tabvint(340)
518 numfram=tabvint(379)
519 numsph =tabvint(380)
520 nselsp =tabvint(381)
521 nspcond =tabvint(382)
522 nsphsym =tabvint(384)
523 maxpjet =tabvint(385)
524 nisp =tabvint(386)
525 kvoisph =tabvint(387)
526 nspbuf =tabvint(388)
527 nispcond=tabvint(389)
528 intbag=tabvint(425)
529 ibagsurf=tabvint(426)
530 nrvolu =tabvint(427)
531 IF(nrvolu==0)nrvolu=24
532 numels10 =tabvint(428)
533 numels20 =tabvint(429)
534 idampg=tabvint(430)
535 iclose=tabvint(431)
536c TABVINT(432)
537c TABVINT(433)
538 nibvel=tabvint(434)
539 lwasph=tabvint(435)
540 IF(lwasph==0)lwasph=16*numsph
541C from version 4.4c :
542 lwasph44c=17*numsph
543 IF(lwasph<lwasph44c)lwasph=lwasph44c
544C-----
545 numels16 =tabvint(436)
546 numels8 =tabvint(437)
547C test for compatibility with old restarts
548 IF(numels8==0)numels8 =
549 . numels - (numels10+numels20+numels16)
550 nrivf = tabvint(438)
551 IF (nrivf==0) nrivf = 1
552 nxframe=tabvint(439)
553 nsphio =tabvint(440)
554 nisphio =tabvint(441)
555 lag_ncf =tabvint(442)
556 lag_nkf =tabvint(443)
557 lag_nhf =tabvint(444)
558 lag_ncl =tabvint(445)
559 lag_nkl =tabvint(446)
560 lag_nhl =tabvint(447)
561 nbcslag =tabvint(448)
562 lactiv =tabvint(449)
563 nactiv =tabvint(450)
564 ndamp =tabvint(451)
565C-----
566 nicbag=tabvint(452)
567C versions 4.4X, X <= D :
568 IF(nicbag==0)nicbag=3
569 nrcbag=tabvint(453)
570C versions 4.4X, X <= D :
571 IF(nrcbag==0)nrcbag=3
572C not available for version <= 4.4D :
573 licbag =tabvint(454)
574 lrcbag =tabvint(455)
575 nibjet =tabvint(456)
576 nrbjet =tabvint(457)
577 nibhol =tabvint(458)
578 nrbhol =tabvint(459)
579 libagjet=tabvint(460)
580 lrbagjet=tabvint(461)
581 libaghol=tabvint(462)
582 lrbaghol=tabvint(463)
583 ngjoint =tabvint(464)
584 nummpc =tabvint(465)
585 lvsphio =tabvint(466)
586! SITHGRPA = TABVINT(467) - obsolete
587! SITHGRPB = TABVINT(468)
588! SITHGRPC = TABVINT(469)
589! SITHGRPD = TABVINT(470)
590! SITHGRPE = TABVINT(471)
591! SITHGRPF = TABVINT(472)
592! SITHGRPG = TABVINT(473)
593! SITHGRPH = TABVINT(474)
594! SITHGRPI = TABVINT(475)
595C
596 nthgrp01(1) = tabvint(476)
597 nthgrp01(2) = tabvint(477)
598 nthgrp01(3) = tabvint(478)
599 nthgrp01(4) = tabvint(479)
600 nthgrp01(5) = tabvint(480)
601 nthgrp01(6) = tabvint(481)
602 nthgrp01(7) = tabvint(482)
603 nthgrp01(8) = tabvint(483)
604 nthgrp01(9) = tabvint(484)
605C
606 nthgrp1(1) = tabvint(485)
607 nthgrp1(2) = tabvint(486)
608 nthgrp1(3) = tabvint(487)
609 nthgrp1(4) = tabvint(488)
610 nthgrp1(5) = tabvint(489)
611 nthgrp1(6) = tabvint(490)
612 nthgrp1(7) = tabvint(491)
613 nthgrp1(8) = tabvint(492)
614 nthgrp1(9) = tabvint(493)
615
616 next = 493
617
618 !output data must be read in case of mcheck/=0 (dynamic size)
619 !otherwise LVARINT is increased to skip related range in TABVINT
620 CALL read_array(outp_v ,tabvint, next, mcheck)
621 CALL read_array(outp_ss,tabvint, next, mcheck)
622 CALL read_array(outp_st,tabvint, next, mcheck)
623 CALL read_array(outp_cs,tabvint, next, mcheck)
624 CALL read_array(outp_ct,tabvint, next, mcheck)
625 CALL read_array(outp_ts,tabvint, next, mcheck)
626 CALL read_array(outp_ps,tabvint, next, mcheck)
627 CALL read_array(outp_pt,tabvint, next, mcheck)
628 CALL read_array(outp_rs,tabvint, next, mcheck)
629 CALL read_array(outp_rt,tabvint, next, mcheck)
630 CALL read_array(anim_v ,tabvint, next, mcheck)
631 CALL read_array(anim_t ,tabvint, next, mcheck)
632 CALL read_array(anim_e ,tabvint, next, mcheck)
633 CALL read_array(anim_ce,tabvint, next, mcheck)
634 CALL read_array(anim_se,tabvint, next, mcheck)
635 CALL read_array(anim_fe,tabvint, next, mcheck)
636 CALL read_array(anim_ct,tabvint, next, mcheck)
637 CALL read_array(anim_st,tabvint, next, mcheck)
638 CALL read_array(anim_ft,tabvint, next, mcheck)
639
640 IF(mcheck/=0)THEN
641 next=next+1
642 anim_m = tabvint(next)
643 next=next+1
644 anim_k = tabvint(next)
645 next=next+1
646 anim_u = tabvint(next)
647 next=next+1
648 anim_mat = tabvint(next)
649 next=next+1
650 ncuts = tabvint(next)
651 next=next+1
652 iepsdot = tabvint(next)
653 next=next+1
654 fmt_ani = tabvint(next)
655 next=next+1
656 anim_vers = tabvint(next)
657 next=next+1
658 izip = tabvint(next)
659 ELSE
660 next=next+9
661 ENDIF
662
663 next=next+1
664 nrbylag = tabvint(next)
665 next = next+1
666 nfvlag = tabvint(next)
667 next = next+1
668 nnprw = tabvint(next)
669 rwall%NNPRW = nnprw
670 next = next+1
671 nrwlag = tabvint(next)
672 rwall%NRWLAG = nrwlag
673 next = next+1
674 nhmax = tabvint(next)
675 next = next+1
676 lagmod = tabvint(next)
677 next = next+1
678 lagopt = tabvint(next)
679 next = next+1
680 nsegflu = tabvint(next)
681 next = next+1
682 ale%GLOBAL%NVCONV = tabvint(next)
683C-----
684 DO i=1,ale%GLOBAL%LCONV
685 ale%GLOBAL%CODV(i)=tabvint(next+i)
686 ENDDO
687C-----
688 next = next+ale%GLOBAL%LCONV+1
689 npebc = tabvint(next)
690 next = next+1
691 nebcs = tabvint(next)
692 next = next+1
693 nrbykin =tabvint(next)
694 next = next+1
695 nbcskin = tabvint(next)
696 next = next+1
697 numels8a = tabvint(next)
698 next = next+1
699 numeltg6 = tabvint(next)
700C--------------------------------------
701 IF(mcheck/=0)THEN
702 DO i=1,12
703 idtmin(i)= tabvint(next + i)
704 ENDDO
705 idtgr(11) = tabvint(next + 13 )
706 nodadt = tabvint(next + 14 )
707 nproc = tabvint(next + 15 )
708C No more used TABVINT(NEXT + 16 )
709 th_vers = tabvint(next + 17 )
710 istatg = tabvint(next + 18 )
711 istat = tabvint(next + 19 )
712 nth = tabvint(next + 20 )
713 nanim = tabvint(next + 21 )
714C No more used TABVINT(NEXT + 22 )
715 imon = tabvint(next + 23 )
716 idum1 = tabvint(next + 24 )
717 DO i =1,10
718 debug(i) = tabvint(next + 24 + i )
719 ENDDO
720 itrace = tabvint(next + 35 )
721 ale%GLOBAL%INCOMP = tabvint(next + 36 )
722
723 ale%SUB%IALESUB = tabvint(next + 38 )
724 impl = tabvint(next + 39 )
725 nitmx = tabvint(next + 40 )
726 ale%UPWIND%UPWM = tabvint(next + 41 )
727 iddw = tabvint(next + 43 )
728 ipread = tabvint(next + 44 )
729 nsflsw = tabvint(next + 45 )
730 ntflsw = tabvint(next + 46 )
731 ! IUSERW = TABVINT(NEXT + 47 )
732 ! NUSERW = TABVINT(NEXT + 48 )
733 irad2r = tabvint(next + 49 )
734 nnoise = tabvint(next + 50 )
735 rnoi = tabvint(next + 51 )
736 noisev = tabvint(next + 52 )
737 noisea = tabvint(next + 53 )
738 noisep = tabvint(next + 54 )
739c link
740 nrlink = tabvint(next + 55 )
741 nalelk = tabvint(next + 56 )
742 lllink = tabvint(next + 57 )
743 ENDIF
744
745 next=next+58
746 CALL read_array(anim_n ,tabvint, next, mcheck)
747C
748 next = next+1
749 ioutp_fmt = tabvint(next)
750 next = next+1
751 npropgi = tabvint(next)
752 next = next+1
753 npropmi = tabvint(next)
754 next = next+1
755 nexmad =tabvint(next)
756 next = next+1
757 nmadprt =tabvint(next)
758 next = next+1
759 nmadsh4 =tabvint(next)
760 next = next+1
761 nmadsh3 =tabvint(next)
762 next = next+1
763 nmadsol =tabvint(next)
764 next = next+1
765 nmadnod =tabvint(next)
766 next = next+1
767 ishfram =tabvint(next)
768C--------------------------------------
769 next = next+1
770 nskyrw0 = tabvint(next)
771 next = next+1
772 nskyrbk0 = tabvint(next)
773 next = next+1
774 nbddnod = tabvint(next)
775 next = next+1
776 nbddnrb = tabvint(next)
777 next = next+1
778 multimax = tabvint(next)
779 next = next+1
780 nbddi2m = tabvint(next)
781 next = next+1
782 niskyi2 = tabvint(next)
783 next = next+1
784 lcni2 = tabvint(next)
785 next = next+1
786 i2nsnt = tabvint(next)
787 next = next+1
788 nskymv0 = tabvint(next)
789 next = next+1
790 nskymvc0 = tabvint(next)
791 next = next+1
792 nskyll0 = tabvint(next)
793 next = next+1
794 nbddncj = tabvint(next)
795 next = next+1
796 nbddnrbm = tabvint(next)
797 next = next+1
798 nskyrbm0 = tabvint(next)
799 next = next+1
800 nnodtcut = tabvint(next)
801 next = next+1
802 nnodlcut = tabvint(next)
803 next = next+1
804 nbddrbe3m = tabvint(next)
805 next = next+1
806 nskyrbe30 = tabvint(next)
807C
808 next = next+1
809 nparir = tabvint(next)
810 next = next+1
811 numnodg = tabvint(next)
812 next = next+1
813 numnodm = tabvint(next)
814 next = next+1
815 numelqg = tabvint(next)
816 next = next+1
817 numelcg = tabvint(next)
818 next = next+1
819 numeltgg = tabvint(next)
820 next = next+1
821C empty field
822 next = next+1
823 numelsg = tabvint(next)
824 next = next+1
825 numelpg = tabvint(next)
826 next = next+1
827 numelrg = tabvint(next)
828 next = next+1
829 numeltrg = tabvint(next)
830 next = next+1
831 nb1dg = tabvint(next)
832 next = next+1
833 nskewag = tabvint(next)
834 next = next+1
835 sectiom = tabvint(next)
836 next = next+1
837 numels16g = tabvint(next)
838 next = next+1
839 num16shift = tabvint(next)
840 next = next+1
841 outyy_fmt = tabvint(next)
842 next = next+1
843 nfxbody=tabvint(next)
844 next = next+1
845 lennod=tabvint(next)
846 next = next+1
847 lenmod=tabvint(next)
848 next = next+1
849 lenglm=tabvint(next)
850 next = next+1
851 lencp=tabvint(next)
852 next = next+1
853 lenlm=tabvint(next)
854 next = next+1
855 lenfls=tabvint(next)
856 next = next+1
857 lendls=tabvint(next)
858 next = next+1
859 lenvar=tabvint(next)
860 next = next+1
861 lenrpm=tabvint(next)
862 next = next+1
863 lenmcd=tabvint(next)
864 next = next+1
865 lenelm=tabvint(next)
866 next = next+1
867 lensig=tabvint(next)
868 next = next+1
869 lengrvi=tabvint(next)
870 next = next+1
871 lengrvr=tabvint(next)
872 next = next+1
873 nbipm=tabvint(next)
874 next = next+1
875 nsfxby=tabvint(next)
876 next = next+1
877 neig=tabvint(next)
878 next=next+1
879 neipm=tabvint(next)
880 next=next+1
881 nerpm=tabvint(next)
882 next=next+1
883 leibuf=tabvint(next)
884 next = next+1
885 nintsub = tabvint(next)
886C--------------------------------------
887 next = next+1
888 irootyy = tabvint(next)
889 IF(mcheck/=0)THEN
890 invers = tabvint(next + 1)
891 ncrst = tabvint(next + 2)
892 irltyp = tabvint(next + 3)
893 iretyp = tabvint(next + 4)
894 multirest = tabvint(next + 5)
895 irprev = tabvint(next + 6)
896 ENDIF
897 next = next + 6
898 next = next+1
899 trimat=tabvint(next)
900 next = next+1
901 sphshift = tabvint(next)
902 next = next+1
903 numsphg = tabvint(next)
904 next = next+1
905 numelxg=tabvint(next)
906 next = next+1
907 nithgr=tabvint(next)
908 next = next+1
909 npebcr=tabvint(next)
910 next = next+1
911 nbi18=tabvint(next)
912 next = next+1
913 nskyi18=tabvint(next)
914C common part SPMD CFD
915 next = next+1
916 nsvois=tabvint(next)
917 next = next+1
918 nqvois=tabvint(next)
919 next = next+1
920 ntgvois = tabvint(next)
921 next = next + 1
922 nrcvvois=tabvint(next)
923 next = next+1
924 nsndvois=tabvint(next)
925 next = next+1
926 nervois=tabvint(next)
927 next = next+1
928 nesvois=tabvint(next)
929 next = next+1
930 ilaw11=tabvint(next)
931 next = next+1
932 nsegcom=tabvint(next)
933 next = next+1
934 numporl=tabvint(next)
935 next = next+1
936 ale%GLOBAL%ICAA=tabvint(next)
937 next = next+1
938 IF(ale%GLOBAL%HLLC_LOWMACH == 0)ale%GLOBAL%HLLC_LOWMACH=tabvint(next)
939 next = next+1
940 isphbuc=tabvint(next)
941 next = next+1
942 nrdamp=tabvint(next)
943 IF(nrdamp==0) nrdamp=4
944 next=next+1
945 nsproc=tabvint(next)
946 next=next+1
947 ndsolv=tabvint(next)
948 next=next+1
949 nsbmax=tabvint(next)
950 next=next+1
951 nsvmaxt=tabvint(next)
952 next=next+1
953 nsmod=tabvint(next)
954 next=next+1
955 nflow=tabvint(next)
956 next=next+1
957 liflow=tabvint(next)
958 next=next+1
959 lrflow=tabvint(next)
960 next=next+1
961 niflow=tabvint(next)
962 next=next+1
963 niioflow=tabvint(next)
964 next=next+1
965 nrflow=tabvint(next)
966 next=next+1
967 nrioflow=tabvint(next)
968 next=next+1
969 icondp=tabvint(next)
970 next=next+1
971 impl_s0=tabvint(next)
972 next=next+1
973 nfasolfr=tabvint(next)
974 next=next+1
975 nsegquadfr=tabvint(next)
976 next=next+1
977 sfr_lagf=tabvint(next)
978 next=next+1
979 sllagf=tabvint(next)
980 next=next+1
981 nventtot=tabvint(next)
982 next = next + 1
983 numelc0= tabvint(next)
984 next = next + 1
985 numeltg0= tabvint(next)
986 next = next + 1
987 numnod0= tabvint(next)
988 next = next + 1
989 nadmesh = tabvint(next)
990 next = next + 1
991 levelmax = tabvint(next)
992 next = next + 1
993 iadmrule = tabvint(next)
994 next = next + 1
995 ksh4tree= tabvint(next)
996 next = next + 1
997 ksh3tree= tabvint(next)
998 next = next + 1
999 kpadmesh= tabvint(next)
1000 next = next + 1
1001 kcontact= tabvint(next)
1002 next = next + 1
1003 kipadmesh= tabvint(next)
1004 next = next + 1
1005 sicontact= tabvint(next)
1006 next = next + 1
1007 srcontact= tabvint(next)
1008 next = next + 1
1009 siskwp= tabvint(next)
1010 next = next + 1
1011 snskwp= tabvint(next)
1012 next = next + 1
1013 sisensp= tabvint(next)
1014 next = next + 1
1015 snsensp= tabvint(next)
1016 next = next + 1
1017 siaccp= tabvint(next)
1018 next = next + 1
1019 snaccp= tabvint(next)
1020 next = next + 1
1021 iisrots = tabvint(next)
1022C
1023C .sta files
1024 next = next + 1
1025 mx_stat0 = tabvint(next)
1026 next = next + 1
1027 istatf = tabvint(next)
1028 IF(mcheck/=0)THEN
1029 DO i=1, mx_stat0
1030 next = next + 1
1031 stat_c(i) = tabvint(next)
1032 ENDDO
1033 ELSE
1034 next = next + mx_stat0
1035 END IF
1036C .dyain files
1037 next = next + 1
1038 dynain_data%IDYNAINF = tabvint(next)
1039 IF(mcheck/=0)THEN
1040 DO i=1, dynain_data%MX_DYNAIN
1041 next = next + 1
1042 dynain_data%DYNAIN_C(i) = tabvint(next)
1043 ENDDO
1044 ELSE
1045 next = next + dynain_data%MX_DYNAIN
1046 END IF
1047 next = next + 1
1048 dynain_data%DYNAIN_CHECK = tabvint(next)
1049C heat transfer + int7
1050 next=next+1
1051 glob_therm%INTHEAT = tabvint(next)
1052C xfem
1053 next=next+1
1054 numelcrk=tabvint(next)
1055 next=next+1
1056 nlevset=tabvint(next)
1057C---------------------------------------
1058C heat transfer
1059 next=next+1
1060 glob_therm%ITHERM_FE=tabvint(next)
1061C Static condensation
1062 next=next+1
1063 istatcnd=tabvint(next)
1064C Unit info
1065 next = next + 1
1066 lunit = tabvint(next)
1067 next = next + 1
1068 nunit0 = tabvint(next)
1069 next = next + 1
1070 nunits = tabvint(next)
1071C
1072 next = next + 1
1073 lsh4trim = tabvint(next)
1074 next = next + 1
1075 lsh3trim = tabvint(next)
1076C
1077 next = next + 1
1078 iadmerrt = tabvint(next)
1079C
1080 next = next + 1
1081 impose_dr = tabvint(next)
1082C
1083 next = next + 1
1084 nintstamp = tabvint(next)
1085C
1086 next = next + 1
1087 inter_ithknod = tabvint(next) !defined in interface module (common_source directory)
1088C
1089C rigid material
1090C
1091 next = next + 1
1092 irigid_mat = tabvint(next)
1093C
1094 next = next + 1
1095 nrbym = tabvint(next)
1096C
1097 next = next + 1
1098 ngslnrbym = tabvint(next)
1099C
1100 next = next + 1
1101 nfrbym = tabvint(next)
1102C
1103 next = next + 1
1104 nirbym = tabvint(next)
1105c
1106 next = next + 1
1107 nskyrbmk0 = tabvint(next)
1108C-----------
1109 next = next + 1
1110 glob_therm%NIFT = tabvint(next)
1111 next = next + 1
1112 glob_therm%NICONV = tabvint(next)
1113 next = next + 1
1114 glob_therm%LFACTHER = tabvint(next)
1115 next = next + 1
1116 glob_therm%NUMCONV = tabvint(next)
1117 next = next + 1
1118 glob_therm%NFXTEMP = tabvint(next)
1119C--------------------------------------
1120C stapming : modification of warning messages in engine
1121 next = next + 1
1122 istamping = tabvint(next)
1123C-----RBE3
1124 next = next + 1
1125 nrbe3 = tabvint(next)
1126 rbe3%NRBE3 = nrbe3 ! Number of RBE3 on domain
1127 next = next + 1
1128 sirbe3 = tabvint(next) ! zize of irbe3 = nrbe3 * irbe3_variables (10)
1129 next = next + 1
1130 slrbe3 = tabvint(next)
1131 rbe3%lrbe3_sz = slrbe3 ! Size of lrbe3
1132 next = next + 1
1133 sfrbe3 = tabvint(next)
1134 rbe3%frbe3_sz = sfrbe3 ! Size of frbe3
1135 next = next + 1
1136 siad_rbe3m = tabvint(next)
1137 ! = nspmd+1
1138 next = next + 1
1139 sfr_rbe3m = tabvint(next)
1140 rbe3%mpi%fr_rbe3_sz = sfr_rbe3m
1141C-------------------------------------------------------------------
1142C spmd rigid material
1143 next = next + 1
1144 sfr_rbym =tabvint(next)
1145 next = next + 1
1146 sfr_rbym2 =tabvint(next)
1147C Noise arrays
1148 next = next + 1
1149 sinoise = tabvint(next)
1150 next = next + 1
1151 sfnoise = tabvint(next)
1152 next = next + 1
1153 nnoise_sav = tabvint(next)
1154C--------------------------------------
1155C shell xfem
1156 next = next + 1
1157 iplyxfem = tabvint(next)
1158 next = next + 1
1159 nplymax = tabvint(next)
1160 next = next + 1
1161 iplybcs = tabvint(next)
1162 next = next + 1
1163 nplyxfe = tabvint(next)
1164 next = next + 1
1165 eplyxfe = tabvint(next)
1166C--------------------------------------
1167 next = next + 1
1168 ngrpart = tabvint(next)
1169C--------------------------------------
1170C Madymo array
1171 next = next + 1
1172 smadidx = tabvint(next)
1173C--------------------------------------
1174C Madcl nodes
1175 next = next + 1
1176 madclnods = tabvint(next)
1177C
1178 next = next + 1
1179 idtmins_old = tabvint(next)
1180 next = next + 1
1181 idtgrs_old = tabvint(next)
1182C--------------------------------------
1183 next = next + 1
1184 ireac = tabvint(next)
1185C--------------------------------------
1186 next = next + 1
1187 snom_sect = tabvint(next)
1188C--------------------------------------
1189 next = next + 1
1190 igrelem = tabvint(next)
1191C--------------------------------------
1192 next = next + 1
1193 nthpart = tabvint(next)
1194C--------------------------------------
1195 next = next + 1
1196 irest_mselt = tabvint(next)
1197C--------------------------------------
1198 next = next + 1
1199 nthread0= tabvint(next)
1200C--------------------------------------
1201 next = next + 1
1202 idrot = tabvint(next)
1203 next = next + 1
1204 glob_therm%NIRADIA = tabvint(next)
1205 next = next + 1
1206 glob_therm%NUMRADIA = tabvint(next)
1207C--------------------------------------
1208 next = next+1
1209 glob_therm%NITFLUX = tabvint(next)
1210 next = next+1
1211 glob_therm%NFXFLUX = tabvint(next)
1212C--------------------------------------
1213 next = next + 1
1214 IF(mcheck/=0)nsmspcg = tabvint(next)
1215C--------------------------------------
1216 next = next + 1
1217 ntable = tabvint(next)
1218C----Rbe2----
1219 next = next+1
1220 nbddnrbe2 = tabvint(next)
1221 next = next + 1
1222 nrbe2 = tabvint(next)
1223 next = next + 1
1224 sirbe2 = tabvint(next)
1225 next = next + 1
1226
1227 slrbe2 = tabvint(next)
1228 next = next + 1
1229 siad_rbe2 = tabvint(next)
1230 next = next + 1
1231 sfr_rbe2 = tabvint(next)
1232 next = next + 1
1233 nhrbe2 = tabvint(next)
1234C
1235C --- simplified ale option + lagrange
1236C
1237 next = next + 1
1238 ialelag = tabvint(next)
1239C
1240 next = next + 1
1241 nrbe2g=tabvint(next)
1242C
1243 next = next + 1
1244 nrbe3g=tabvint(next)
1245C
1246C--------------------------------------
1247C ... shell 4-N + shell 3-N (X-FEM)
1248 next=next+1
1249 icrack3d=tabvint(next)
1250C ply xfem (parit on)
1251 next=next+1
1252 lcnepxfem =tabvint(next)
1253 next = next + 1
1254 sithvar = tabvint(next)
1255 next=next+1
1256 nplynodg = tabvint(next)
1257C --- Need NXREF Option for Modif files / SPMD
1258 next = next + 1
1259 nxref = tabvint(next)
1260C
1261 next = next + 1
1262 idtmins_int_old = tabvint(next)
1263C
1264 next = next +1
1265 slcfield = tabvint(next)
1266 next = next + 1
1267 nload = tabvint(next)
1268 next = next + 1
1269 idamp_rdof = tabvint(next)
1270 next = next + 1
1271 lskew = tabvint(next)
1272 next = next + 1
1273 liskn = tabvint(next)
1274 next = next + 1
1275 lfxvelr = tabvint(next)
1276C
1277 next = next + 1
1278 nloadc = tabvint(next)
1279 next = next + 1
1280 slloadp = tabvint(next)
1281 next = next + 1
1282 nloadp = tabvint(next)
1283 next = next + 1
1284 nloadp_hyd = tabvint(next)
1285 next = next + 1
1286 nintloadp = tabvint(next)
1287 next = next + 1
1288 nintloadp21 = tabvint(next)
1289C xfem for layered shell (cracks)
1290 next=next+1
1291 lcnecrkxfem=tabvint(next)
1292 next = next + 1
1293 ncrkxfe = tabvint(next)
1294 next = next + 1
1295 ecrkxfe = tabvint(next)
1296 next = next + 1
1297 nlevmax = tabvint(next)
1298 next = next + 1
1299 ncrknodg = tabvint(next)
1300 next = next + 1
1301 ienrnod = tabvint(next)
1302 next = next + 1
1303 sibufssg_io = tabvint(next)
1304 next = next + 1
1305 nibsph = tabvint(next)
1306 next = next + 1
1307 slinale = tabvint(next)
1308 next = next + 1
1309 nsubmod = tabvint(next)
1310 next = next + 1
1311 ale%GLOBAL%NALENOVS = tabvint(next)
1312C
1313 next = next + 1
1314 lnopt1 = tabvint(next)
1315 next = next + 1
1316 lilset1 = tabvint(next)
1317 next = next + 1
1318 lislin1 = tabvint(next)
1319 next = next + 1
1320 ltitr = tabvint(next)
1321 next = next + 1
1322 lisurf1 = tabvint(next)
1323 next = next + 1
1324 lipart1 = tabvint(next)
1325 next = next + 1
1326 sizfield = tabvint(next)
1327 next = next + 1
1328 sizloadp = tabvint(next)
1329 next = next + 1
1330 i2nsn25 = tabvint(next)
1331 next = next + 1
1332 lvoisph = tabvint(next)
1333 next = next + 1
1334 IF(mcheck /= 0 ) anim_ply = tabvint(next)
1335 next = next + 1
1336 sibordnode = tabvint(next)
1337 next = next + 1
1338 siedgesh = tabvint(next)
1339C multilayer xfem shell
1340 next = next + 1
1341 ecrkxfec = tabvint(next)
1342 next = next + 1
1343 ecrkxfetg = tabvint(next)
1344 next = next + 1
1345 numedges = tabvint(next)
1346c inter22
1347 next = next + 1
1348 int22 = tabvint(next)
1349 next = next + 1
1350 i22_aleul= tabvint(next)
1351 next = next + 1
1352 nsub22 = tabvint(next)
1353 next = next + 1
1354 ioutp22 = tabvint(next)
1355 next = next + 1
1356 i22law37 = tabvint(next)
1357 next = next + 1
1358 i22law51 = tabvint(next)
1359c
1360 next = next + 1
1361 icrash = tabvint(next)
1362C
1363 next = next + 1
1364 siad_edge = tabvint(next)
1365 next = next + 1
1366 sfr_edge = tabvint(next)
1367 next = next + 1
1368 sfr_nbedge= tabvint(next)
1369C
1370 next = next + 1
1371 nbgauge = tabvint(next)
1372 next = next + 1
1373 llaccelm = tabvint(next)
1374 next = next + 1
1375 llgauge = tabvint(next)
1376c
1377 next = next + 1
1378 sigaup= tabvint(next)
1379 next = next + 1
1380 sngaup= tabvint(next)
1381C
1382 next=next+1
1383 sdd_r2r_elem=tabvint(next)
1384 next=next+1
1385 r2r_siu=tabvint(next)
1386C heat trasfert for interface Type21 : main temperature is variable
1387 next = next + 1
1388 ftempvar21 = tabvint(next)
1389C
1390 next = next + 1
1391 numnodxfe=tabvint(next)
1392 next = next + 1
1393 snodglobxfe=tabvint(next)
1394c
1395 next = next + 1
1396 nsphsol= tabvint(next)
1397 next = next + 1
1398 first_sphsol= tabvint(next)
1399C
1400 next = next + 1
1401 iddom=tabvint(next)
1402 next = next + 1
1403 nummat0=tabvint(next)
1404 next = next + 1
1405 numgeo0=tabvint(next)
1406C
1407 next = next + 1
1408 intplyxfem = tabvint(next)
1409C
1410 next = next + 1
1411 i22len = tabvint(next)
1412C
1413 next = next + 1
1414 ncluster = tabvint(next)
1415C
1416 next = next + 1
1417 nxlaymax = tabvint(next)
1418C
1419 next = next+1
1420 srthbuf = tabvint(next)
1421C
1422 next = next + 1
1423 numelig3d = tabvint(next)
1424 next = next + 1
1425 skxig3d = tabvint(next)
1426 next = next + 1
1427 sixig3d = tabvint(next)
1428 next = next + 1
1429 sknot = tabvint(next)
1430 next = next + 1
1431 sknotlocpc = tabvint(next)
1432 next = next + 1
1433 sknotlocel = tabvint(next)
1434 next = next + 1
1435 numfakenodigeo = tabvint(next)
1436 next = next + 1
1437 deg_max = tabvint(next)
1438C
1439 next = next + 1
1440 nxel = tabvint(next)
1441C
1442 next = next + 1
1443 nfilsol = tabvint(next)
1444C
1445 next=next+1
1446 sensors%SFSAV = tabvint(next)
1447C
1448 next=next+1
1449 sensors%STABSEN = tabvint(next)
1450C Stack info
1451 next = next + 1
1452 ipart_stack = tabvint(next)
1453 next = next + 1
1454 npt_stack = tabvint(next)
1455 next = next + 1
1456 ns_stack = tabvint(next)
1457C
1458 next = next + 1
1459 alefvm_param%IEnabled = tabvint(next)
1460 next = next + 1
1461 alefvm_param%ISOLVER = tabvint(next)
1462C MUSCL keyword, compression parameter beta
1463 next = next + 1
1464 alemuscl_param%IALEMUSCL = tabvint(next)
1465 alemuscl_param%IALEMUSCL = alemuscl_param%IALEMUSCL * (1-alemuscl_param%I_MUSCL_OFF) !/ALE/MUSCL/OFF : freform.F
1466 next = next + 1
1467 alemuscl_param%BETA = tabvint(next)
1468 next = next + 1
1469 alemuscl_param%I_LAW = tabvint(next)
1470C ITETRA=3
1471 next=next+1
1472 isfem =tabvint(next)
1473C /LOAD/PBLAST
1474 next = next + 1
1475 nloadp_f = tabvint(next)
1476 next = next + 1
1477 pblast%NLOADP_B = tabvint(next)
1478 next = next + 1
1479 loads%NLOAD_CYL = tabvint(next)
1480C
1481 next = next + 1
1482 IF(mcheck/=0) idtmin(51)= tabvint(next)
1483C /DT/FVMBAG
1484 next = next + 1
1485 IF(mcheck/=0) idtmin(52)= tabvint(next)
1486C
1487 next = next + 1
1488 nbr_gpmp = tabvint(next)
1489C
1490 next = next + 1
1491 taille = tabvint(next)
1492C
1493 next = next + 1
1494 ninter25 = tabvint(next)
1495C
1496 next = next + 1
1497 numnor = tabvint(next)
1498C
1499 next = next + 1
1500 nbddnort = tabvint(next)
1501C
1502 next = next + 1
1503 nbddnor_max = tabvint(next)
1504C
1505 next = next + 1
1506 nbccfr25 = tabvint(next)
1507C
1508 next = next + 1
1509 nbccnor = tabvint(next)
1510C
1511 next = next + 1
1512 nbddedgt = tabvint(next)
1513C
1514 next = next + 1
1515 nbddedg_max = tabvint(next)
1516C
1517 next = next + 1
1518 nsnt25 = tabvint(next)
1519C
1520 next = next + 1
1521 nrtmx25 = tabvint(next)
1522C--------------------------------------
1523C drape option
1524 next = next + 1
1525 ndrape = tabvint(next)
1526C Element activation
1527 next = next + 1
1528 lractiv = tabvint(next)
1529C ID_PLY -> animation file
1530 next = next + 1
1531 numply = tabvint(next)
1532C ID_PLY -> animation file
1533 next = next + 1
1534 numstack = tabvint(next)
1535C
1536 next = next + 1
1537 tabmp_l = tabvint(next)
1538C SOL2SPH flag
1539 next = next + 1
1540 sol2sph_flag = tabvint(next)
1541 next = next + 1
1542 itsol2sph = tabvint(next)
1543C ISFINT flag
1544 next = next + 1
1545 IF(ale%GLOBAL%ISFINT==0)THEN
1546 !no engine keyword (/ALE/SOLVER) affecting ISFINT, then use starter setting
1547 ale%GLOBAL%ISFINT = tabvint(next)
1548 ENDIF
1549C ICNDS10 dyna condensation for S10
1550 next = next + 1
1551 ns10e = tabvint(next)
1552 next = next + 1
1553 lcncnd = tabvint(next)
1554 next = next + 1
1555 siad_cndm = tabvint(next)
1556 next = next + 1
1557 sfr_cndm = tabvint(next)
1558C---- old local sys of shell3n
1559 next = next+1
1560 ish3nfram =tabvint(next)
1561 next = next+1
1562 stab_mat = tabvint(next)
1563C SMS auto switch
1564 next = next + 1
1565 isms_selec = tabvint(next)
1566 next = next + 1
1567 isused = tabvint(next)
1568 multi_fvm%IS_USED = .false.
1569 IF (isused == 1) THEN
1570 multi_fvm%IS_USED = .true.
1571 multi_fvm%BETA = alemuscl_param%BETA
1572 ENDIF
1573 next = next + 1
1574 isused = tabvint(next)
1575 multi_fvm%NBMAT = isused
1576C---- old Isolid17
1577 next = next+1
1578 is17_old =tabvint(next)
1579C IH3D -> H3D count
1580 next = next + 1
1581 h3d_data%IH3D = tabvint(next)
1582C
1583 next = next + 1
1584 numels10g = tabvint(next)
1585C
1586 next = next + 1
1587 numels20g = tabvint(next)
1588!
1589 next = next + 1
1590 ale%GLOBAL%I_DT_NODA_ALE_ON = tabvint(next)
1591c interface adhesion
1592 next = next + 1
1593 interadhesion = tabvint(next)
1594C Friction interface
1595 next = next + 1
1596 ninterfric = tabvint(next)
1597!
1598 next = next + 1
1599 lensubs = tabvint(next)
1600!
1601 next = next + 1
1602 lenigrnod = tabvint(next)
1603!
1604 next = next + 1
1605 lenigrbric = tabvint(next)
1606!
1607 next = next + 1
1608 lenigrquad = tabvint(next)
1609!
1610 next = next + 1
1611 lenigrsh4n = tabvint(next)
1612!
1613 next = next + 1
1614 lenigrsh3n = tabvint(next)
1615!
1616 next = next + 1
1617 lenigrtrus = tabvint(next)
1618!
1619 next = next + 1
1620 lenigrbeam = tabvint(next)
1621!
1622 next = next + 1
1623 lenigrspri = tabvint(next)
1624!
1625 next = next + 1
1626 lenigrpart = tabvint(next)
1627!
1628 next = next + 1
1629 lenisurf = tabvint(next)
1630!
1631 next = next + 1
1632 lenislin = tabvint(next)
1633!
1634 IF(mcheck/=0)THEN
1635 next = next + 1
1636 idts6 = tabvint(next)
1637 next = next + 1
1638 idt1sol = tabvint(next)
1639 next = next + 1
1640 idt1tet10= tabvint(next)
1641 next = next + 1
1642 idt1sh = tabvint(next)
1643
1644 ELSE
1645 next = next + 4
1646 END IF
1647!
1648! H3D -> number of inter skid lines
1649 next = next + 1
1650 nintskidold= tabvint(next)
1651! H3D ANIM-> Is maximum contact force /contact pressure computed
1652 next = next + 1
1653 ifcontmax = tabvint(next)
1654 next = next + 1
1655 ifcontpmax = tabvint(next)
1656 next = next + 1
1657 ifcont2max = tabvint(next)
1658 next = next + 1
1659 ifcontp2max= tabvint(next)
1660 next = next + 1
1661 ifcont2min = tabvint(next)
1662 next = next + 1
1663 ifcontp2min= tabvint(next)
1664 next = next + 1
1665 output%DATA%S_EFRIC = tabvint(next)
1666 next = next + 1
1667 output%DATA%S_EFRICINT = tabvint(next)
1668 next = next + 1
1669 output%DATA%S_EFRICG = tabvint(next)
1670 next = next + 1
1671 output%DATA%S_EFRICINTG = tabvint(next)
1672 next = next + 1
1673 output%DATA%NINEFRIC = tabvint(next)
1674 next = next + 1
1675 output%DATA%NINEFRIC_STAMP = tabvint(next)
1676C NITSCHE METHOD
1677 next=next+1
1678 nitsche = tabvint(next)
1679C Pinching
1680 next = next+1
1681 npinch = tabvint(next)
1682 next = next+1
1683 nploadpinch = tabvint(next)
1684
1685!
1686 next=next+1
1687 IF(ale%UPWIND%I_SUPG_ON_OFF==0) ale%UPWIND%I_SUPG_ON_OFF = tabvint(next)
1688
1689 next = next + 1
1690 nledge = tabvint(next)
1691!
1692 next = next + 1
1693 inispri = tabvint(next)
1694!
1695 next = next + 1
1696 sens_nstop = tabvint(next)
1697 next = next + 1
1698 sens_nstat = tabvint(next)
1699 next = next + 1
1700 sens_noutp = tabvint(next)
1701 next = next + 1
1702 sens_nanim = tabvint(next)
1703 IF (mcheck /= 0) THEN ! rerun after chkpoint restart => use previous run settings
1704 sensors%NSTOP = sens_nstop
1705 sensors%NSTAT = sens_nstat
1706 sensors%NOUTP = sens_noutp
1707 sensors%NANIM = sens_nanim
1708 ALLOCATE (sensors%STOP(sensors%NSTOP))
1709 ALLOCATE (sensors%OUTP(sensors%NOUTP))
1710 ALLOCATE (sensors%STAT(sensors%NSTAT))
1711 ALLOCATE (sensors%ANIM(sensors%NANIM))
1712 sensors%STOP = 0
1713 sensors%OUTP = 0
1714 sensors%STAT = 0
1715 sensors%ANIM = 0
1716 END IF
1717!
1718 next = next + 1
1719 nibcld = tabvint(next)
1720 next = next + 1
1721 lfaccld = tabvint(next)
1722!
1723 next = next + 1
1724 nigrv = tabvint(next)
1725 next = next + 1
1726 lfacgrv = tabvint(next)
1727!
1728C BCS/CYCLIC
1729 next = next + 1
1730 nbcscyc = tabvint(next)
1731 next = next + 1
1732 sibcscyc = tabvint(next)
1733 next = next + 1
1734 slbcscyc = tabvint(next)
1735
1736 next = next + 1
1737 inter18_autoparam = tabvint(next)
1738 next = next + 1
1739 IF(tabvint(next) == 1)inter18_is_variable_gap_defined = .true.
1740
1741 IF ((iale > 0 .AND. ale%GRID%NWALE_RST == 6) .OR.
1742 . inter18_autoparam == 1 .OR. inter18_is_variable_gap_defined) THEN
1743 next = next + 1
1744 ne_slercvois = tabvint(next)
1745 next = next + 1
1746 ne_slesdvois = tabvint(next)
1747 next = next + 1
1748 ne_nsvois = tabvint(next)
1749 next = next + 1
1750 ne_nqvois = tabvint(next)
1751 next = next + 1
1752 ne_ntgvois = tabvint(next)
1753 ENDIF
1754
1755 next = next + 1
1756 pdel = tabvint(next)
1757
1758C-----Seatbelts - slipring / retractors
1759 next = next + 1
1760 nslipring = tabvint(next)
1761 next = next + 1
1762 nretractor = tabvint(next)
1763 next = next + 1
1764 nslipring_g = tabvint(next)
1765 next = next + 1
1766 nretractor_g = tabvint(next)
1767 next = next + 1
1768 n_anchor_remote = tabvint(next)
1769 next = next + 1
1770 n_anchor_remote_send = tabvint(next)
1771 next = next + 1
1772 n_seatbelt = tabvint(next)
1773 next = next + 1
1774 n_seatbelt_2d = tabvint(next)
1775C
1776 next = next + 1
1777 lastanimcycle = tabvint(next)
1778
1779 next = next + 1
1780 lasth3dcycle = tabvint(next)
1781
1782 next = next + 1
1783 lastdyncycle = tabvint(next)
1784
1785 next = next + 1
1786 laststatcycle = tabvint(next)
1787C
1788C---- Starter input version
1789 next = next + 1
1790 st_invers = tabvint(next)
1791C
1792C---- Number of RBODY before merge
1793 next = next + 1
1794 nrbody0 = tabvint(next)
1795C---- MAX (/H3D)
1796 next = next + 1
1797 lmax_dis = tabvint(next)
1798 next = next + 1
1799 lmax_vel = tabvint(next)
1800 next = next + 1
1801 lmax_nsig = tabvint(next)
1802 next = next + 1
1803 lmax_nstra = tabvint(next)
1804C---- Time history strain flag
1805 next = next + 1
1806 th_strain = tabvint(next)
1807C draped shell element
1808 next = next + 1
1809 numelc_drape= tabvint(next)
1810 next = next + 1
1811 numeltg_drape= tabvint(next)
1812C
1813 next = next + 1
1814 inter_parameters%ISTIF_DT = tabvint(next)
1815C
1816 next = next + 1
1817 inter_parameters%INTCAREA = tabvint(next)
1818C
1819C /LOAD/CLOAD
1820 next = next + 1
1821 loads%NLOAD_CLOAD = tabvint(next)
1822C /LOAD/PLOAD
1823 next = next + 1
1824 loads%NLOAD_PLOAD = tabvint(next)
1825C-----/DAMP/VREL
1826 next = next + 1
1827 ndamp_vrel = tabvint(next)
1828!
1829 next = next + 1
1830 inter_parameters%INT25_EROSION_SOLID = tabvint(next)
1831C /INIVEL
1832 next = next + 1
1833 loads%NINIVELT = tabvint(next)
1834! ALE%REZON
1835 next = next + 1
1836 ale%REZON%NUM_NUVAR_MAT = tabvint(next)
1837 next = next + 1
1838 ale%REZON%NUM_NUVAR_EOS = tabvint(next)
1839C
1840 ! Initialize variables from Restart
1841 ! ----------------------------------
1842 ! Skews
1843 skews%SKEW_NUMBER=numskw ! Copy the Skew Number in SKEW_ type
1844 skews%N_SKEW_VAR= lskew ! Copy the number of variables per skew
1845 skews%TOTAL_SKEW_NUMBER = numskw + 1 + min(1,nspcond)*numsph + nsubmod ! skew array contains deck skews + 1 skew per SPH + 1 Skew per Submodel
1846C--------------------------------------
1847 RETURN
1848 END
1849!||====================================================================
1850!|| rdcomr ../engine/source/output/restart/rdcomm.F
1851!||--- called by ------------------------------------------------------
1852!|| rdresb ../engine/source/output/restart/rdresb.F
1853!||--- uses -----------------------------------------------------
1854!|| ale_mod ../common_source/modules/ale/ale_mod.F
1855!|| glob_therm_mod ../common_source/modules/mat_elem/glob_therm_mod.F90
1856!|| multimat_param_mod ../common_source/modules/multimat_param_mod.F90
1857!|| output_mod ../common_source/modules/output/output_mod.f90
1858!|| parameters_mod ../common_source/modules/interfaces/parameters_mod.F
1859!|| sph_mod ../engine/share/modules/sph_mod.F
1860!|| state_mod ../common_source/modules/state_mod.F
1861!||====================================================================
1862 SUBROUTINE rdcomr(LMXVREA,TABVREA,DYNAIN_DATA,INTER_PARAMETERS,
1863 . OUTPUT,GLOB_THERM)
1864C-----------------------------------------------
1865C M o d u l e s
1866C-----------------------------------------------
1867 USE output_mod
1868 USE state_mod
1869 USE ale_mod
1870 USE parameters_mod
1871 USE output_mod
1872 USE multimat_param_mod , ONLY : m51_iflg6
1873 USE glob_therm_mod
1874 USE sph_mod, ONLY : xi_res,yi_res,zi_res
1875C-----------------------------------------------
1876C I m p l i c i t T y p e s
1877C-----------------------------------------------
1878#include "implicit_f.inc"
1879C-----------------------------------------------
1880C A r g u m e n t s
1881C-----------------------------------------------
1882 INTEGER LMXVREA
1883 my_real TABVREA(LMXVREA)
1884 TYPE (DYNAIN_DATABASE), INTENT(INOUT) :: DYNAIN_DATA
1885 TYPE (PARAMETERS_) ,INTENT(INOUT):: INTER_PARAMETERS
1886 TYPE(output_), INTENT(INOUT) :: OUTPUT !< output structure
1887 TYPE (glob_therm_) , INTENT(INOUT) :: GLOB_THERM
1888C-----------------------------------------------
1889C C o m m o n B l o c k s
1890C-----------------------------------------------
1891#include "com06_c.inc"
1892#include "com08_c.inc"
1893#include "cong2_c.inc"
1894#include "couple_c.inc"
1895#include "lagmult.inc"
1896#include "remesh_c.inc"
1897#include "scr06_c.inc"
1898#include "scr07_c.inc"
1899#include "scr11_c.inc"
1900#include "scr18_c.inc"
1901#include "scr19_c.inc"
1902#include "scrnoi_c.inc"
1903#include "sms_c.inc"
1904#include "sphcom.inc"
1905#include "sysunit.inc"
1906#include "statr_c.inc"
1907#include "inter22.inc"
1908#include "itet2_c.inc"
1909C-----------------------------------------------
1910#include "ddspmd_c.inc"
1911C-----------------------------------------------
1912C L o c a l V a r i a b l e s
1913C-----------------------------------------------
1914 INTEGER I,RVARINT,J
1915C-----------------------------------------------
1916C-----
1917 tt =tabvrea(1)
1918 dt1 =tabvrea(2)
1919 dt2 =tabvrea(3)
1920 dt12 =tabvrea(4)
1921 dt2old =tabvrea(5)
1922 tstop =tabvrea(6)
1923 output%DATA%INOUT%DM_IN = tabvrea(7)
1924 output%TH%THIS =tabvrea(8)
1925C-----
1926 output%DATA%INOUT%DM_OUT = tabvrea(9)
1927 output%TH%DTHIS =tabvrea(10)
1928 dtfac =tabvrea(11)
1929 dtmin =tabvrea(12)
1930 volmin =tabvrea(13)
1931 output%TH%WFEXT =tabvrea(14)
1932 reint =tabvrea(15)
1933 output%DATA%INOUT%DE_IN = tabvrea(16)
1934 output%DATA%INOUT%DE_OUT = tabvrea(17)
1935 output%DTANIM =tabvrea(18)
1936 output%TANIM =tabvrea(19)
1937 ureint =tabvrea(20)
1938 econtv =tabvrea(21)
1939 ehour =tabvrea(22)
1940 dtoutp =tabvrea(23)
1941 toutp =tabvrea(24)
1942 t1s =tabvrea(25)
1943 dt2s =tabvrea(26)
1944 epor =tabvrea(27)
1945 ! unused USREINT=TABVREA(28)
1946 r2rfx1 =tabvrea(29)
1947 r2rfx2 =tabvrea(30)
1948 dampa =tabvrea(31)
1949 dampb =tabvrea(32)
1950 dampw =tabvrea(33)
1951C-----
1952 ale%GRID%ALPHA = tabvrea(34)
1953 ale%GRID%GAMMA = tabvrea(35)
1954 ale%GRID%VGX = tabvrea(36)
1955 ale%GRID%VGY = tabvrea(37)
1956 ale%GRID%VGZ = tabvrea(38)
1957C-----
1958 spasort=tabvrea(39)
1959 spatrue=tabvrea(40)
1960 DO i=1,9
1961 output%TH%DTHIS1(i) = tabvrea(40+i)
1962 output%TH%THIS1(i) = tabvrea(49+i)
1963 ENDDO
1964 mass0 = tabvrea(59)
1965 entot0 = tabvrea(60)
1966 lagm_tol =tabvrea(61)
1967 lag_alph =tabvrea(62)
1968 lag_alphs=tabvrea(63)
1969
1970C-----
1971C read/write next real variable at TABVREA(LVR44A+1) if LVARREA > LVR44A, else, set its value to zero.
1972C (whatever the real common to store the variable).
1973C--------------------------------------
1974 IF(mcheck/=0)THEN
1975 DO i=1,12
1976 dtfac1(i) = tabvrea(63 + i)
1977 dtmin1(i) = tabvrea(63 + 12 +i)
1978 ENDDO
1979 beta = tabvrea(88)
1980 betate = tabvrea(89)
1981 period = tabvrea(90)
1982 dmtmxs = tabvrea(91)
1983 dmnmxs = tabvrea(92)
1984 dmnmxs = tabvrea(93)
1985 dtin = tabvrea(94)
1986 dtmx = tabvrea(95)
1987 madul = tabvrea(96)
1988 madut = tabvrea(97)
1989 madum = tabvrea(98)
1990 ale%SUB%DTFSUB = tabvrea(99)
1991 ale%SUB%DTMSUB = tabvrea(100)
1992 eps = tabvrea(101)
1993 eps2 = tabvrea(102)
1994 ale%UPWIND%CUPWM = tabvrea(103)
1995c noise
1996 tnoise = tabvrea(104)
1997 dtnoise= tabvrea(105)
1998 ENDIF
1999
2000 dwmad = tabvrea(106)
2001 fac_mass = tabvrea(107)
2002 fac_length = tabvrea(108)
2003 fac_time = tabvrea(109)
2004 t1sh = tabvrea(110)
2005 ale%UPWIND%UPWSM = tabvrea(111)
2006
2007 rvarint = 111
2008
2009 rvarint = rvarint + 1
2010 tpsref = tabvrea(rvarint)
2011C THERMAL ANALYSIS
2012 rvarint = rvarint + 1
2013 glob_therm%HEAT_STORED = tabvrea(rvarint)
2014 rvarint = rvarint + 1
2015 glob_therm%HEAT_MECA = tabvrea(rvarint)
2016 rvarint = rvarint + 1
2017 glob_therm%HEAT_CONV = tabvrea(rvarint)
2018 rvarint = rvarint + 1
2019 glob_therm%HEAT_RADIA = tabvrea(rvarint)
2020 rvarint = rvarint + 1
2021 glob_therm%HEAT_FFLUX = tabvrea(rvarint)
2022C SOL1TNL array
2023 DO i=0, maxlaw
2024 rvarint = rvarint + 1
2025 sol1tnl(i,1) = tabvrea(rvarint)
2026 rvarint = rvarint + 1
2027 sol1tnl(i,2) = tabvrea(rvarint)
2028 ENDDO
2029C SOL8TNL array
2030 DO i=0, maxlaw
2031 rvarint = rvarint + 1
2032 sol8tnl(i,1) = tabvrea(rvarint)
2033 rvarint = rvarint + 1
2034 sol8tnl(i,2) = tabvrea(rvarint)
2035 ENDDO
2036
2037C TET4TNL array
2038 DO i=0, maxlaw
2039 rvarint = rvarint + 1
2040 tet4tnl(i,1) = tabvrea(rvarint)
2041 ENDDO
2042C TET10TNL array
2043 DO i=0, maxlaw
2044 rvarint = rvarint + 1
2045 tet10tnl(i,1) = tabvrea(rvarint)
2046 ENDDO
2047
2048C SHTNL array
2049 DO i=0, maxlaw
2050 DO j=0,3
2051 rvarint = rvarint + 1
2052 shtnl(i,j,1) = tabvrea(rvarint)
2053 rvarint = rvarint + 1
2054 shtnl(i,j,2) = tabvrea(rvarint)
2055 rvarint = rvarint + 1
2056 shtnl(i,j,3) = tabvrea(rvarint)
2057 ENDDO
2058 ENDDO
2059C TRITNL array
2060 DO i=0, maxlaw
2061 DO j=0,3
2062 rvarint = rvarint + 1
2063 tritnl(i,j,1) = tabvrea(rvarint)
2064 rvarint = rvarint + 1
2065 tritnl(i,j,2) = tabvrea(rvarint)
2066 rvarint = rvarint + 1
2067 tritnl(i,j,3) = tabvrea(rvarint)
2068 ENDDO
2069 ENDDO
2070C SOLTELT array
2071 DO i=1,8
2072 rvarint = rvarint + 1
2073 soltelt(i) = tabvrea(rvarint)
2074 ENDDO
2075
2076 rvarint = rvarint + 1
2077 tet4telt(1) = tabvrea(rvarint)
2078 rvarint = rvarint + 1
2079 tet10telt(1) = tabvrea(rvarint)
2080
2081C SHTELT array
2082 DO i=1,10
2083 rvarint = rvarint + 1
2084 shtelt(i) = tabvrea(rvarint)
2085 ENDDO
2086C TRITELT array
2087 DO i=1,10
2088 rvarint = rvarint + 1
2089 tritelt(i) = tabvrea(rvarint)
2090 ENDDO
2091C
2092C .sta files
2093 rvarint = rvarint + 1
2094 dtstat = tabvrea(rvarint)
2095 rvarint = rvarint + 1
2096 tstat = tabvrea(rvarint)
2097C
2098C .dynain files
2099 rvarint = rvarint + 1
2100 dynain_data%DTDYNAIN = tabvrea(rvarint)
2101 rvarint = rvarint + 1
2102 dynain_data%TDYNAIN = tabvrea(rvarint)
2103C
2104 rvarint = rvarint + 1
2105 tadmesh = tabvrea(rvarint)
2106 rvarint = rvarint + 1
2107 dtadmesh= tabvrea(rvarint)
2108C
2109 rvarint = rvarint + 1
2110 dtfacs_old = tabvrea(rvarint)
2111 rvarint = rvarint + 1
2112 dtmins_old = tabvrea(rvarint)
2113C
2114 rvarint = rvarint + 1
2115 IF(mcheck/=0) tol_sms = tabvrea(rvarint)
2116C abf files
2117 DO i=1,10
2118 rvarint = rvarint + 1
2119 tabfis(i) = tabvrea(rvarint)
2120 rvarint = rvarint + 1
2121 dtabf(i) = tabvrea(rvarint)
2122 rvarint = rvarint + 1
2123 dtabfwr(i) = tabvrea(rvarint)
2124 rvarint = rvarint + 1
2125 tabfwr(i) = tabvrea(rvarint)
2126 ENDDO
2127C
2128 rvarint = rvarint + 1
2129 ratio22 = tabvrea(rvarint)
2130C
2131 rvarint = rvarint + 1
2132 m51_iflg6 = nint(tabvrea(rvarint))
2133C
2134 rvarint = rvarint + 1
2135 ms_nd = tabvrea(rvarint)
2136 rvarint = rvarint + 1
2137 dmsi2 = tabvrea(rvarint)
2138 rvarint = rvarint + 1
2139 msi20 = tabvrea(rvarint)
2140C
2141 rvarint = rvarint + 1
2142 IF(mcheck/=0) dtfac1(51)= tabvrea(rvarint)
2143 rvarint = rvarint + 1
2144 IF(mcheck/=0) dtmin1(51)= tabvrea(rvarint)
2145C /DT/FVMBAG
2146 rvarint = rvarint + 1
2147 IF(mcheck/=0) dtfac1(52)= tabvrea(rvarint)
2148 rvarint = rvarint + 1
2149 IF(mcheck/=0) dtmin1(52)= tabvrea(rvarint)
2150C
2151 rvarint = rvarint + 1
2152 IF(mcheck/=0) dtfac1(102)= tabvrea(rvarint)
2153 rvarint = rvarint + 1
2154 IF(mcheck/=0) dtmin1(102)= tabvrea(rvarint)
2155C
2156 rvarint = rvarint + 1
2157 IF(mcheck/=0) dtfacs_int_old= tabvrea(rvarint)
2158 rvarint = rvarint + 1
2159 IF(mcheck/=0) dtmins_int_old= tabvrea(rvarint)
2160C
2161C /ADYREL :BETATE_0, EINT_0,ENCIN_0,PIMAX ,PCMAX,PINT ,PCIN ,FIL_
2162 rvarint = rvarint + 1
2163 betate_0 = tabvrea(rvarint)
2164 rvarint = rvarint + 1
2165 eint_0 = tabvrea(rvarint)
2166 rvarint = rvarint + 1
2167 encin_0 = tabvrea(rvarint)
2168 rvarint = rvarint + 1
2169 pimax = tabvrea(rvarint)
2170 rvarint = rvarint + 1
2171 pcmax = tabvrea(rvarint)
2172 rvarint = rvarint + 1
2173 pint = tabvrea(rvarint)
2174 rvarint = rvarint + 1
2175 pcin = tabvrea(rvarint)
2176 rvarint = rvarint + 1
2177 fil_ie(1:4)=tabvrea(rvarint:rvarint+3)
2178 rvarint = rvarint + 4
2179 fil_ke(1:4)=tabvrea(rvarint:rvarint+3)
2180 rvarint = rvarint + 4
2181 nfirst = tabvrea(rvarint)
2182C
2183 rvarint = rvarint + 1
2184 IF(mcheck/=0) ale%GRID%VGY0= tabvrea(rvarint)
2185 rvarint = rvarint + 1
2186 IF(mcheck/=0) ale%GRID%VGZ0= tabvrea(rvarint)
2187C
2188 rvarint = rvarint + 1
2189 econtd= tabvrea(rvarint)
2190 rvarint = rvarint + 1
2191 econt_cumu= tabvrea(rvarint)
2192
2193C INITIAL MASS IN STARTER
2194 rvarint = rvarint + 1
2195 mass0_start = tabvrea(rvarint)
2196 rvarint = rvarint + 1
2197 percent_addmass_old = tabvrea(rvarint)
2198C ANIM FILES FINAL TIME
2199 rvarint = rvarint + 1
2200 output%TANIM_STOP = tabvrea(rvarint)
2201C INITIAL COORDS OF SPH RESERVE
2202 rvarint = rvarint + 1
2203 xi_res = tabvrea(rvarint)
2204 rvarint = rvarint + 1
2205 yi_res = tabvrea(rvarint)
2206 rvarint = rvarint + 1
2207 zi_res = tabvrea(rvarint)
2208C
2209 rvarint = rvarint + 1
2210 inter_parameters%DT_STIFINT = tabvrea(rvarint)
2211C KINETIC ENERGY
2212 rvarint = rvarint + 1
2213 encin = tabvrea(rvarint)
2214 rvarint = rvarint + 1
2215 enrot = tabvrea(rvarint)
2216
2217 RETURN
2218 END
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
ale_mod
Definition ale_mod.F:150
ale_mod::ale
type(ale_) ale
Definition ale_mod.F:253
alefvm_mod
Definition alefvm_mod.F:88
alefvm_mod::alefvm_param
type(alefvm_param_), target alefvm_param
Definition alefvm_mod.F:121
alemuscl_mod
Definition alemuscl_mod.F:55
alemuscl_mod::alemuscl_param
type(alemuscl_param_) alemuscl_param
Definition alemuscl_mod.F:76
damp_mod
Definition damp_mod.F:33
damp_mod::ndamp_vrel
integer ndamp_vrel
Definition damp_mod.F:44
drape_mod
Definition drape_mod.F:54
drape_mod::numeltg_drape
integer numeltg_drape
Definition drape_mod.F:92
drape_mod::numelc_drape
integer numelc_drape
Definition drape_mod.F:92
h3d_mod
Definition h3d_mod.F:311
outmax_mod
Definition outmax_mod.F:52
outmax_mod::ifcontp2max
integer ifcontp2max
Definition outmax_mod.F:69
outmax_mod::lmax_vel
integer lmax_vel
Definition outmax_mod.F:61
outmax_mod::ifcontmax
integer ifcontmax
Definition outmax_mod.F:69
outmax_mod::ifcont2max
integer ifcont2max
Definition outmax_mod.F:69
outmax_mod::ifcontp2min
integer ifcontp2min
Definition outmax_mod.F:69
outmax_mod::lmax_nstra
integer lmax_nstra
Definition outmax_mod.F:63
outmax_mod::ifcontpmax
integer ifcontpmax
Definition outmax_mod.F:69
outmax_mod::lmax_dis
integer lmax_dis
Definition outmax_mod.F:60
outmax_mod::lmax_nsig
integer lmax_nsig
Definition outmax_mod.F:62
outmax_mod::ifcont2min
integer ifcont2min
Definition outmax_mod.F:69
parameters_mod
Definition parameters_mod.F:49
pinchtype_mod
Definition pinchtype_mod.F:69
pinchtype_mod::nploadpinch
integer nploadpinch
Definition pinchtype_mod.F:74
pinchtype_mod::npinch
integer npinch
Definition pinchtype_mod.F:73
pload_cyl_mod
Definition pload_cyl_mod.F:35
sph_mod
Definition sph_mod.F:31
state_mod
Definition state_mod.F:54
th_mod
Definition th_mod.F:42
th_mod::th_has_noda_pext
integer th_has_noda_pext
Definition th_mod.F:121
rdcomr
subroutine rdcomr(lmxvrea, tabvrea, dynain_data, inter_parameters, output, glob_therm)
Definition rdcomm.F:1864
rdcomi
subroutine rdcomi(output, tabvint, multi_fvm, h3d_data, dynain_data, inter_parameters, sensors, loads, mat_elem, skews, glob_therm, pblast, rbe3, rwall)
Definition rdcomm.F:58
read_array
subroutine read_array(array, tabvint, lvarint, m_check)
Definition read_array.F:29
h3d_mod::h3d_database
Definition h3d_mod.F:438