OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
contrl.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!|| contrl ../starter/source/starter/contrl.F
25!||--- called by ------------------------------------------------------
26!|| starter0 ../starter/source/starter/starter0.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../starter/source/output/message/message.F
29!|| arret ../starter/source/system/arret.F
30!|| ascii_encoding_mu_letter ../starter/source/starter/ascii_encoding_mu_letter.F90
31!|| contrbe2 ../starter/source/constraints/general/rbe2/hm_read_rbe2.F
32!|| contrbe3 ../starter/source/constraints/general/rbe3/contrbe3.F
33!|| find_dt1brick_engine ../starter/source/starter/contrl.F
34!|| hm_elem_count ../starter/source/devtools/hm_reader/hm_elem_count.F
35!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.F
36!|| hm_get_string ../starter/source/devtools/hm_reader/hm_get_string.F
37!|| hm_option_count ../starter/source/devtools/hm_reader/hm_option_count.F
38!|| hm_option_next ../starter/source/devtools/hm_reader/hm_option_next.f
39!|| hm_option_read_key ../starter/source/devtools/hm_reader/hm_option_read_key.F
40!|| hm_option_start ../starter/source/devtools/hm_reader/hm_option_start.F
41!|| hm_prelce16s ../starter/source/elements/reader/hm_read_solid.F
42!|| hm_preread_node ../starter/source/elements/reader/hm_preread_node.F
43!|| hm_preread_part ../starter/source/model/assembling/hm_read_part.F
44!|| hm_preread_properties ../starter/source/properties/hm_preread_properties.f
45!|| hm_read_ale_grid ../starter/source/general_controls/ale_grid/hm_read_ale_grid.F
46!|| hm_read_ale_muscl ../starter/source/general_controls/ale_cfd/hm_read_ale_muscl.F
47!|| hm_read_ale_solver ../starter/source/general_controls/ale_cfd/hm_read_ale_solver.F
48!|| hm_read_analy ../starter/source/general_controls/computation/hm_read_analy.F
49!|| hm_read_caa ../starter/source/general_controls/computation/hm_read_caa.F
50!|| hm_read_definter ../starter/source/general_controls/default_values/hm_read_definter.F
51!|| hm_read_defshell ../starter/source/general_controls/default_values/hm_read_defshell.F
52!|| hm_read_defsolid ../starter/source/general_controls/default_values/hm_read_defsolid.F
53!|| hm_read_implicit ../starter/source/general_controls/computation/hm_read_implicit.F
54!|| hm_read_inista ../starter/source/initial_conditions/inista/hm_read_inista.F
55!|| hm_read_ioflag ../starter/source/general_controls/inputoutput/hm_read_ioflag.f
56!|| hm_read_lagmul ../starter/source/tools/lagmul/hm_read_lagmul.F
57!|| hm_read_refsta ../starter/source/loads/reference_state/refsta/hm_read_refsta.F
58!|| hm_read_sms ../starter/source/general_controls/computation/hm_read_sms.F
59!|| hm_read_sphglo ../starter/source/general_controls/computation/hm_read_sphglo.F
60!|| hm_read_spmd ../starter/source/general_controls/computation/hm_read_spmd.F
61!|| hm_read_unit ../starter/source/general_controls/computation/hm_read_unit.F
62!|| hm_read_upwind ../starter/source/general_controls/computation/hm_read_upwind.f
63!|| init_def_elem ../starter/source/modules/defaults_mod.F90
64!|| init_def_zero ../starter/source/modules/defaults_mod.F90
65!|| istr ../starter/source/tools/univ/istr.F
66!|| nbadigemesh ../starter/source/elements/ige3d/nbadigemesh.F
67!|| nbadmesh ../starter/source/model/remesh/nbadmesh.F
68!|| nbsph ../starter/source/elements/sph/nbsph.F
69!||--- uses -----------------------------------------------------
70!|| defaults_mod ../starter/source/modules/defaults_mod.F90
71!|| detonators_mod ../starter/share/modules1/detonators_mod.F
72!|| format_mod ../starter/share/modules1/format_mod.F90
73!|| grp_size_mod ../starter/share/modules1/grp_size_mod.f
74!|| hm_option_read_mod ../starter/share/modules1/hm_option_read_mod.F
75!|| inivol_def_mod ../starter/share/modules1/inivol_mod.F
76!|| message_mod ../starter/share/message_module/message_mod.F
77!|| refsta_mod ../starter/share/modules1/refsta_mod.F
78!|| restmod ../starter/share/modules1/restart_mod.F
79!|| submodel_mod ../starter/share/modules1/submodel_mod.F
80!||====================================================================
81 SUBROUTINE contrl(MULTI_FVM,LSUBMODEL,IS_DYNA,DETONATORS,USER_WINDOWS,MAT_ELEM,
82 . NAMES_AND_TITLES,LIPART1,DEFAULTS,GLOB_THERM,PBLAST,OUTPUT)
83C----------------------------------------------------------
84C M o d u l e s
85C-----------------------------------------------
86 USE message_mod
87 USE multi_fvm_mod
88 USE submodel_mod
89 USE pinchtype_mod
90 USE hm_option_read_mod
91 USE check_mod
92 USE inoutfile_mod
93 USE setdef_mod
94 USE detonators_mod
95 USE inivol_def_mod , ONLY : num_inivol
96 USE refsta_mod
97 USE ale_ebcs_mod
98 USE restmod
99 USE grp_size_mod
100 USE anim_mod
101 USE alemuscl_mod , only:alemuscl_param
102 USE alefvm_mod , only:alefvm_param
103 USE user_windows_mod
104 USE ale_mod
105 USE mat_elem_mod
106 USE names_and_titles_mod
107 USE bcs_mod , only : bcs
108 USE defaults_mod
109 USE format_mod
110 use glob_therm_mod
111 USE pblast_mod
112 USE output_mod , ONLY : output_
113C-----------------------------------------------
114C I m p l i c i t T y p e s
115C-----------------------------------------------
116#include "implicit_f.inc"
117C-----------------------------------------------
118C G l o b a l P a r a m e t e r s
119C-----------------------------------------------
120#include "mvsiz_p.inc"
121C-----------------------------------------------
122C C o m m o n B l o c k s
123C-----------------------------------------------
124#include "units_c.inc"
125#include "warn_c.inc"
126#include "com01_c.inc"
127#include "com04_c.inc"
128#include "com06_c.inc"
129#include "com08_c.inc"
130#include "com09_c.inc"
131#include "com10_c.inc"
132#include "com_xfem1.inc"
133#include "intstamp_c.inc"
134#include "random_c.inc"
135#include "scr03_c.inc"
136#include "scr05_c.inc"
137#include "scr06_c.inc"
138#include "scr10_c.inc"
139#include "scr12_c.inc"
140#include "scr15_c.inc"
141#include "scr16_c.inc"
142#include "scr22_c.inc"
143#include "scr23_c.inc"
144#include "titr_c.inc"
145#include "param_c.inc"
146#include "sphcom.inc"
147#include "lagmult.inc"
148#include "fxbcom.inc"
149#include "scr14_c.inc"
150#include "remesh_c.inc"
151#include "sysunit.inc"
152#include "commandline.inc"
153#include "r2r_c.inc"
154#include "userlib.inc"
155#include "spmd_c.inc"
156#include "drape_c.inc"
157#include "inter22.inc"
158#include "boltpr_c.inc"
159#include "ige3d_c.inc"
160#include "com_engcards_c.inc"
161#include "sms_c.inc"
162C-----------------------------------------------
163C D u m m y A r g u m e n t s
164C-----------------------------------------------
165 TYPE(multi_fvm_struct), INTENT(IN) :: MULTI_FVM
166 TYPE(submodel_data) LSUBMODEL(NSUBMOD)
167 INTEGER,INTENT(IN)::IS_DYNA
168 TYPE(detonators_struct_) :: DETONATORS
169 TYPE(user_windows_), INTENT(INOUT) :: USER_WINDOWS
170 TYPE(mat_elem_), INTENT(INOUT) :: MAT_ELEM
171 TYPE(names_and_titles_),INTENT(INOUT) :: NAMES_AND_TITLES !< NAMES_AND_TITLES host the input deck names and titles for outputs
172 INTEGER,INTENT(IN) :: LIPART1 !< Number of variables of IPART
173 TYPE(defaults_), INTENT(INOUT) :: DEFAULTS !< Defaults mod
174 type (glob_therm_) ,intent(inout) :: glob_therm
175 TYPE(pblast_), INTENT(INOUT) :: PBLAST !< PBLAST load type
176 TYPE(output_),INTENT(INOUT) :: OUTPUT
177C-----------------------------------------------
178C L o c a l V a r i a b l e s
179C-----------------------------------------------
180 INTEGER I, J, N, ISFIL, IPCT, IBID,INSH,
181 . ihbe_dd,iparith,jale,jeul,
182 . iun,nshfram,jupw,imat,ig,ish3n_dd,iframe_dds,
183 . ipla_dds,npts_dd,uid,ihbe,ismstr,iplas,icpre,icstr,npt,
184 . isrot,l1,npts, nimpdisp,nimpvel,nimpacc,icr,
185 . stat,ierror,id,idx,idy,idz,nctrl,bid,nsh3nfram,ns17_old,
186 . ignore_spmd,ignore_threads,got_variable,krot,old_rsb,
187 . nrafx,nrafy,nrafz,nfunct0,ntable0,ntable1,nsensor,
188 . irfe,irfl, hm_ninter_def,nb_ams,numnusr,nperturb_hm,icr3,
189 . npyfun
190 INTEGER IARCHS(8)
191 INTEGER IS_BEGIN,SCHAR
192 INTEGER IHBE_DS,ISST_DS,IPLA_DS,IFRAME_DS,ITET4_D,ITET10_D,ICPRE_D,IMAS_DS,ICONTROL_D,
193 . ihbe_d,ipla_d,istr_d,ithk_d,ishea_d,isst_d,
194 . ish3n_d, istra_d,npts_d,idril_d,ioffset_d,def_inter(100)
195 my_real dtini, dtx ,rbid
196 CHARACTER (LEN=NCHARLINE) :: CART
197 CHARACTER (LEN=NCHARLINE) :: XRFILE ! NCHARLINE as #define is set to 500 in Starter
198 CHARACTER (LEN=NCHARLINE) :: KEY
199 CHARACTER (LEN=NCHARLINE) :: KEY2
200 CHARACTER (LEN=NCHARLINE) :: TMPLINE
201 CHARACTER (LEN=NCHARLINE) :: LINE
202 CHARACTER (LEN=NCHARTITLE) :: TITR
203 CHARACTER (LEN=255) :: STR_NBTHREADS
204 CHARACTER MESS*40, ERRMSG*40
205 CHARACTER*3 :: LABEL_DEF,LABEL_ROT
206C-----------------------------------------------
207C OpenMP specific
208#if defined(_OPENMP)
209 INTEGER OMP_GET_THREAD_NUM, OMP_GET_NUM_THREADS,NTHREAD1
210 EXTERNAL omp_get_thread_num, omp_get_num_threads
211 CHARACTER (LEN=255) :: STR
212#endif
213 INTEGER NTHREAD_S
214 my_real :: dt_input
215C
216 INTEGER , DIMENSION(:), ALLOCATABLE :: NPLY,NSUB,NISUB
217
218 INTEGER :: NITER,IFORM
219
220 INTEGER :: LEN_TMP_NAME
221 CHARACTER(len=4096) :: TMP_NAME
222 LOGICAL :: IS_AVAILABLE
223 INTEGER :: NUMTITLE, NGRTRIA, IDUMMY, NANIM_VERS
224 INTEGER :: NALEMUSCL
225 INTEGER :: NB_INISHE,NB_INISH3,NB_INIBRI,NB_INIQUAD,
226 . nb_inibeam,nb_initruss,nb_inisprig,nb_inisphcel
227 INTEGER :: LEN_LINE
228 CHARACTER*20 UNIT_NAME
229 INTEGER IS_U_STRING
230C-----------------------------------------------
231C E x t e r n a l F u n c t i o n
232C-----------------------------------------------
233 INTEGER NBCS_LAGM
234 EXTERNAL nbcs_lagm
235C-----------------------------------------------
236C data
237C-----------------------------------------------
238 DATA s0file/' '/
239 DATA xrfile/' '/
240 DATA iun/1/
241
242C-----------------------------------------------
243C S o u r c e L i n e s
244C-----------------------------------------------
245
246C=======================================================================
247C Read title
248C=======================================================================
249 CALL hm_option_count('/TITLE', numtitle)
250 IF (numtitle > 0) THEN
251 CALL hm_option_start('/TITLE')
252 CALL hm_option_read_key(lsubmodel, option_titr = line)
253 CALL hm_get_string('MY_TITLE', line, ncharline, is_available)
254 ELSE
255 line = ' '
256 ENDIF
257! Standard output
258 WRITE(istdo,'(A)') line(1:len_trim(line))
259
260 ! Store the input deck title in Structure
261 len_line= min(len_trim(line),ltitle) ! Truncate to LTITLE
262 names_and_titles%TITLE(1:len_line)=line(1:len_line)
263
264 imot=0
265 reel=zep66
266C=======================================================================
267C UNIT SYSTEM READER
268C First option to be read in order to convert other options.
269C=======================================================================
270 WRITE(istdo,'(A)')titre(12)
271 CALL hm_option_count('/UNIT',nunit0)
272 unitab%NUNIT0 = nunit0
273c
274 CALL hm_option_count('/BEGIN',is_begin)
275
276 schar = 20
277 IF (is_begin /= 0) THEN
278 CALL hm_option_start('/BEGIN')
279 CALL hm_option_next()
280
281 CALL hm_get_intv('INVERS', invers, is_available, lsubmodel)
282 invers_src = invers
283 invers_init = invers
284
285 CALL hm_get_string('length_inputunit_code',keyli,schar,is_available)
286 CALL hm_get_string('mass_inputunit_code',keymi,schar,is_available)
287 CALL hm_get_string('time_inputunit_code',keyti,schar,is_available)
288 CALL hm_get_string('length_workunit_code',keyl,schar,is_available)
289 CALL hm_get_string('mass_workunit_code',keym,schar,is_available)
290 CALL hm_get_string('time_workunit_code',keyt,schar,is_available)
291
292 !convert into ascci format specific encoding of greek letter \mu
293 CALL ascii_encoding_mu_letter(keyli, keymi, keyti, keyl, keym, keyt)
294
295 IF(invers <= 90)THEN
296 ! /UNIT/LENGTH/[value]
297 ! /UNIT/MASS/[value]
298 ! /UNIT/TIME/[value]
299 CALL hm_option_start('/UNIT')
300 DO n=1,nunit0
301 CALL hm_option_read_key(lsubmodel,option_id = id)
302 unit_name = ''
303 CALL hm_get_string('UNIT_NAME',unit_name,2*ncharfield,is_available)
304c
305 IF(unit_name == 'LENGTH') THEN
306 CALL hm_get_string('TYPE_UNIT',keyl,2*ncharfield,is_available)
307 ENDIF
308 IF(unit_name == 'MASS') THEN
309 CALL hm_get_string('TYPE_UNIT',keym,2*ncharfield,is_available)
310 ENDIF
311 IF(unit_name == 'TIME') THEN
312 CALL hm_get_string('TYPE_UNIT',keyt,2*ncharfield,is_available)
313 ENDIF
314 ENDDO
315 ENDIF !(INVERS <= 90)
316
317 ENDIF
318 CALL hm_read_unit(unitab,lsubmodel)
319C=======================================================================
320C READING CONTROL CARDS
321C=======================================================================
322 WRITE(istdo,'(A)')titre(10)
323C-------------------------------------------------------------------
324C READING OF /SPMD
325C-------------------------------------------------------------------
326 CALL hm_read_spmd(lsubmodel)
327C
328 ignore_spmd=0
329 IF (got_ncpu ==1) THEN
330 IF (nspmd/=0) THEN
331 ignore_spmd=1
332 END IF
333 nspmd = ncpu
334 ENDIF
335C
336C thread number
337C
338 ignore_threads=0
339 got_variable=0
340#if defined(_OPENMP)
341 str = ' '
342 CALL getenv('OMP_NUM_THREADS',str)
343 nthread1=0
344C nthread1 : thread number determined by environment variable
345 READ(str,'(I10)',err=999)nthread1
346 IF(nthread1>0)THEN
347 nthread = nthread1
348 ENDIF
349 IF (got_nth ==1) THEN
350 IF (nthread/=0) THEN
351 ignore_threads=1
352 END IF
353 nthread = nth
354 ELSE
355 IF(nthread1>0)THEN
356 got_variable=1
357 END IF
358 ENDIF
359C case no -nt, no OMP_NUM_THREADS, no /SPMD
360 nthread = max(nthread,1)
361 nthread_s = nthread
362 nthread_r2r = nthread
363 CALL omp_set_num_threads(nthread_s)
364
365c else (open mp not defined)
366#elif 1
367 IF (got_nth ==1) THEN
368 IF (nthread/=0) THEN
369 ignore_threads=1
370 END IF
371 nthread = nth
372 ENDIF
373 nthread_s = 1
374#endif
375C-------------------------------------------------------------------
376C /ARCH option not defined, IARCH= IBUILTIN & Grpsiz = ARCHINFO(IBUILTIN)
377 iarch = ibuiltin
378 nvsiz = archinfo(ibuiltin,1)
379 ivector = archinfo(ibuiltin,2)
380
381 ! ---------------------------
382 ! -grp_size hidden option
383 IF(grp_size_bool) THEN
384 nvsiz = grp_size
385 ENDIF
386 ! ---------------------------
387C-------------------------------------------------------------------
388C READING OF /IOFLAG
389C-------------------------------------------------------------------
390 CALL hm_read_ioflag(lsubmodel)
391C-------------------------------------------------------------------
392C READING OF /ANALY
393C-------------------------------------------------------------------
394 ipari0= 1
395 integ8= 0
396 CALL hm_read_analy(nanaly,iparith,ipari0,lsubmodel)
397 n2d = nanaly
398C-------------------------------------------------------------------
399C READING OF /IMPLICIT
400C-------------------------------------------------------------------
401 CALL hm_read_implicit(lsubmodel)
402C-------------------------------------------------------------------
403C READING OF /AMS
404C-------------------------------------------------------------------
405 CALL hm_option_count('/ams', NB_AMS)
406 CALL HM_READ_SMS(LSUBMODEL)
407C-------------------------------------------------------------------
408C READING OF /CAA
409C-------------------------------------------------------------------
410 CALL HM_READ_CAA(LSUBMODEL)
411C-------------------------------------------------------------------
412C READING OF /RANDOM
413C-------------------------------------------------------------------
414 CALL HM_OPTION_COUNT('/random',NRAND)
415C-------------------------------------------------------------------
416C READING OF /LAGMUL
417C-------------------------------------------------------------------
418 CALL HM_READ_LAGMUL(LSUBMODEL)
419C-------------------------------------------------------------------
420 CALL HM_OPTION_COUNT('/PRIVATE/metadata/fatxml',IPRIVATE)
421 IF(IPRIVATE > 0 ) IPRIVATE = 1
422C-------------------------------------------------------------------
423C READING OF /STAMPING
424C-------------------------------------------------------------------
425 CALL HM_OPTION_COUNT('/stamping',ISTAMPING)
426 IF(ISTAMPING > 0 ) ISTAMPING = 1
427C-------------------------------------------------------------------
428 ICRASH = 0
429C-------------------------------------------------------------------
430C READING DEFAULT VALUES FOR SHELLS, SOLIDS, INTERFACES ...
431C-------------------------------------------------------------------
432!--- remove first /DEF/SOLID in src21_c.inc
433 CALL init_def_zero(DEFAULTS)
434!---- obsolet flag values fixed here
435 IPLA_DS =2 ! obsolet
436 ISTR_D=1 ! istrain
437 ISHEA_D=0 ! Old hidden flag in /DEF_SHELL
438 INER_9_12 = ZERO
439 INSH = 0 ! Old hidden flag in /DEF_SHELL
440 NPTS_D = -1 ! Old hidden flag in /DEF_SHELL
441C-------------------------------------------------------------------
442C READING OF SHELL PROPERTIES DEFAULT VALUES
443C-------------------------------------------------------------------
444 CALL HM_READ_DEFSHELL(LSUBMODEL,DEFAULTS%SHELL)
445C-------------------------------------------------------------------
446C READING OF SOLID PROPERTIES DEFAULT VALUES
447C-------------------------------------------------------------------
448 CALL HM_READ_DEFSOLID(LSUBMODEL,DEFAULTS%SOLID)
449C----------------------------------------------
450 IS17OLD = 1 ! To be cleaned
451C-----
452 IF(INSH==1)INER_9_12 = TWELVE
453 IF(INSH==2)INER_9_12 = SIXTEEN
454 IF(INSH==3)INER_9_12 = FOURTY8
455 IF(INSH==4)INER_9_12 = THIRTY
456 IF(INSH==5)INER_9_12 = NINE
457! default values in def_solid/def_shell
458 CALL init_def_elem(N2D,IIMPLICIT,DEFAULTS)
459C-------------------------------------------------------------------
460C READING OF INTERFACES DEFAULT VALUES
461C-------------------------------------------------------------------
462 HM_NINTER_DEF = 0
463 DEF_INTER(1:100) = 0
464 CALL HM_OPTION_COUNT('/default/inter',HM_NINTER_DEF)
465 CALL HM_READ_DEFINTER(HM_NINTER_DEF,DEF_INTER,LSUBMODEL)
466 DEFAULTS%interface%DEF_INTER(1:100) = DEF_INTER(1:100)
467C-------------------------------------------------------------------
468 CALL HM_OPTION_COUNT('/intthick/v5',IINTTHICK)
469 IF (IINTTHICK > 1) THEN
470 CALL ANCMSG(MSGID=725,MSGTYPE=MSGERROR,ANMODE=ANINFO,I1=IINTTHICK)
471 ENDIF
472C-------------------------------------------------------------------
473C READGING OF /SHFRAM
474C-------------------------------------------------------------------
475 ISHFRAM = 0
476 CALL HM_OPTION_COUNT('/shfra/v4',NSHFRAM)
477 IF (NSHFRAM == 1) THEN
478 ISHFRAM = 2
479 ELSEIF (NSHFRAM > 1) THEN
480 CALL ANCMSG(MSGID=546,MSGTYPE=MSGERROR,ANMODE=ANINFO,I1=NSHFRAM)
481 ENDIF
482C-------------------------------------------------------------------
483C READGING OF /SH_3NFR
484C-------------------------------------------------------------------
485 ISH3NFRAM = 0
486C=======================================================================
487C Initialisations vs DOMDEC (Must stay after reading the control cards, because of output to STDO)
488C=======================================================================
489 IF(DECNEQ==0) DECNEQ = 100
490 IF(DECTYP < 0) THEN
491 DECTYP = - DECTYP
492 EDGE_FILTERING = 0
493 ELSE
494 EDGE_FILTERING = 1
495 ENDIF
496 OLD_RSB=0
497 IF(DECTYP == 2)THEN
498 DECTYP=0
499 OLD_RSB=1
500 ENDIF
501 IF(DECTYP==0)THEN
502.AND. IF(NB_AMS==0IIMPLICIT==0) THEN
503 DECTYP = 3
504 ELSE
505 DECTYP = 5
506 END IF
507 END IF
508
509 DDNOD_SMS=0
510 IF(DECTYP==7)THEN
511 IF(NB_AMS/=0)THEN
512 DECTYP = 5
513 DDNOD_SMS=1
514 ELSEIF(IIMPLICIT/=0) THEN
515 DECTYP = 5
516 ELSE
517 DECTYP = 3
518 END IF
519 END IF
520 IF(NSPMD < 1) NSPMD=1
521 IF(NSPMD > PARASIZ) NSPMD=PARASIZ
522 IF(NTHREAD < 1) NTHREAD=1
523C maximum number of SMP threads equal to NTHMAX
524 IF(NTHREAD > NTHMAX) NTHREAD=NTHMAX
525C
526 IF (GOT_INSPIRE_ALM == 1)THEN
527 IF (NTHREAD_S==1)THEN
528 WRITE(ISTDO,'(a,i4,a)')' .. solver running on ',NTHREAD_S,' thread'
529 ELSE
530 WRITE(ISTDO,'(a,i4,a)')' .. solver running on ',NTHREAD_S,' threads'
531 ENDIF
532 ELSE
533 IF (NTHREAD_S==1)THEN
534 WRITE(ISTDO,'(a,i4,a)')' .. starter running on ',NTHREAD_S,' thread'
535 ELSE
536 WRITE(ISTDO,'(a,i4,a)')' .. starter running on ',NTHREAD_S,' threads'
537 ENDIF
538 ENDIF
539
540C=======================================================================
541C OPTIONS NUMBERING
542C=======================================================================
543C WRITE(ISTDO,'(/,A)')' .. OPTIONS COUNTING'
544C-------------------------------------------------------------------
545C User-Defined Nodes & Cnodes numbering
546C Check nodes within some tolerance and possibly merge nodes
547C-------------------------------------------------------------------
548 CALL CPP_NODES_COUNT(NUMNUSR,NUMCNOD)
549C------
550C Pre-read Nodes & Cnodes and compute NUMNOD taking into account that some nodes may be merged.
551 CALL HM_PREREAD_NODE(UNITAB,LSUBMODEL,NUMNUSR,IS_DYNA) ! NUMNOD is computed here
552C------
553C numnod=0 : ask user to provide relevant input file.
554 IF(NUMNOD==0)THEN
555 CALL ANCMSG(MSGID=3,MSGTYPE=MSGERROR,ANMODE=ANINFO)
556 CALL ARRET(2)
557 ENDIF
558C-------------------------------------------------------------------
559 ISUMNX = 0
560 NANIM1D = 0
561 MAXNX = 0
562 NANIM2D = 0
563 NANIM3D = 0
564 NTHREAD_R2R = 1
565 IBID = 0
566 RBID = ZERO
567 NPINCH = 0
568 NSUBDOM = 0
569C---------------------------------------------------
570 CALL HM_OPTION_COUNT('/subdomain',NSUBDOM)
571 CALL HM_OPTION_COUNT('/extern/link',NR2RLNK)
572C---------------------------------------------------
573 LENMOD=0
574 LENGLM=0
575 LENCP=0
576 LENLM=0
577 LENFLS=0
578 LENDLS=0
579 LENVAR=0
580 LENRPM=0
581 LENMCD=0
582 LENELM=0
583 LENSIG=0
584 LENGRVI=0
585 LENGRVR=0
586C--- Xfem -------------------------------------------------------
587C---
588C WARNING: NUMNOD, NUMELC, NUMELTG ==> will be changed (by xfem)
589C---
590 NLEVMAX = 0
591 CALL HM_OPTION_COUNT('/inicrack',NINICRACK)
592C---------------------------------------------------
593C COUNTING ELEMENTS
594C---------------------------------------------------
595 NUMELX = 0
596 NUMBRICK = 0
597 NUMTETRA4 = 0
598 NUMPENTA6 = 0
599 NUMELS10 = 0
600 NUMELS20 = 0
601 CALL HM_ELEM_COUNT('xelem',NUMELX,IS_DYNA)
602 CALL HM_ELEM_COUNT('brick',NUMBRICK,IS_DYNA)
603 CALL HM_ELEM_COUNT('tetra4',NUMTETRA4,IS_DYNA)
604 CALL HM_ELEM_COUNT('penta6',NUMPENTA6,IS_DYNA)
605 NUMELS8 = NUMBRICK+NUMTETRA4+NUMPENTA6
606C
607 CALL HM_ELEM_COUNT('tetra10',NUMELS10,IS_DYNA)
608 IF (IS_DYNA == 0) CALL HM_ELEM_COUNT('brick20',NUMELS20,IS_DYNA)
609 CALL HM_ELEM_COUNT('shel16',NUMELS16,IS_DYNA)
610C
611C
612 NUMELS = NUMELS8+NUMELS10+NUMELS20+NUMELS16
613C
614C Look for /DT1/BRICK & /DT1TET10 /DTTSH in RADIOSS Engine input deck
615C
616 IF(NUMELS>0)THEN
617 CALL FIND_DT1BRICK_ENGINE()
618 END IF
619C-----
620 CALL HM_ELEM_COUNT('shell',NUMELC,IS_DYNA)
621.AND. IF(NUMELC > 0 NANALY /= 0)THEN
622 CALL ANCMSG(MSGID=285,MSGTYPE=MSGERROR,ANMODE=ANINFO)
623 NUMELC = 0
624 ENDIF
625C-----
626 CALL HM_ELEM_COUNT('sh3n',NUMELTG,IS_DYNA)
627.AND. IF(NUMELTG > 0 NANALY /= 0)THEN
628 CALL ANCMSG(MSGID=287,MSGTYPE=MSGERROR,ANMODE=ANINFO)
629 NUMELTG = 0
630 ENDIF
631C-----
632 CALL HM_ELEM_COUNT('truss' ,NUMELT,IS_DYNA)
633 CALL HM_ELEM_COUNT('beam' ,NUMELP,IS_DYNA)
634 CALL HM_ELEM_COUNT('spring',NUMELR,IS_DYNA)
635C-----
636 CALL HM_ELEM_COUNT('rivet',NRIVET,IS_DYNA)
637C-----
638 NUMELQ = 0
639 CALL HM_ELEM_COUNT('quad',numelq,is_dyna)
640C-----
641 CALL hm_elem_count('TRIA',numeltria,is_dyna)
642 numeltg = numeltg + numeltria !2d shell + 3d tria (only one type depending on N2D flag)
643C-----
644C---------------------------------------------------
645 CALL hm_option_count('PART',npart)
646 IF(npart==0)THEN
647 CALL ancmsg(msgid=1114,
648 . msgtype=msgwarning,
649 . anmode=aninfo)
650 ENDIF
651 CALL hm_option_count('SUBSET',nsubs)
652C add 1 : for global subset
653 nsubs = nsubs+1
654C---------------------------------------------------
655 CALL hm_option_count('/THPART',nthpart)
656C---------------------------------------------------
657C LOOKING FOR /ADMESH and COUNTING ADDITIONAL NODES & ELEMENTS
658C---------------------------------------------------
659 CALL hm_option_count('/ADMESH/GLOBAL',nadmeshg)
660 CALL hm_option_count('/ADMESH/SET',nadmeshset)
661 CALL hm_option_count('/ADMESH/STATE',nadmeshstat)
662 nadmesh = nadmeshg + nadmeshset + nadmeshstat
663 numnod0 = numnod
664 numelc0 = numelc
665 numeltg0 = numeltg
666 istatcnd = 0
667 iadmerrt = 0
668 IF(nadmesh/=0)THEN
669C
670C NUMNOD, NUMELC, NUMELTG will be changed
671 CALL nbadmesh(lsubmodel,numnusr,unitab)
672 END IF
673C-----
674 IF(nspmd > 1)THEN
675C Tests option non disponible en SPMD (ie NSPMD > 1)
676 IF(nadmesh/=0)THEN
677 CALL ancmsg(msgid=704,
678 . msgtype=msgerror,
679 . anmode=aninfo)
680 END IF
681 END IF
682C---------------------------------------------------
683C LOOKING FOR IGE REFINEMENT and COUNTING ADDITIONAL NODES & ELEMENTS
684C---------------------------------------------------
685 numnodige0 = numnod
686 numelig3d0 = numelig3d
687 IF(nrafmax/=0)THEN
688C
689C NUMNOD, NUMELIG3D will be changed, SIXIG3D and KXIG3D too
690 CALL nbadigemesh(lsubmodel,numnusr)
691 END IF
692C---------------------------------------------------
693C MATERIALS (+1 materiau fictif ressort)
694C---------------------------------------------------
695 hm_nummat = 0
696 CALL hm_option_count('MATERIAL',hm_nummat)
697 nummat = hm_nummat + 1
698 mat_elem%NUMMAT = nummat
699C---------------------------------------------------
700C PROPERTIES
701C---------------------------------------------------
702 hm_numgeo = 0
703 CALL hm_option_count('PROPERTY',hm_numgeo)
704 numgeo = hm_numgeo
705C-----
706 CALL hm_option_count('/PLY', numply)
707 CALL hm_option_count('/DRAPE',ndrape)
708 CALL hm_option_count('/STACK',numstack)
709C---------------------------------------------------
710C TOOLS
711C---------------------------------------------------
712 CALL hm_option_count('/ACCEL', naccelm)
713C-----
714 CALL hm_option_count('/GAUGE', nbgauge)
715C-----
716 CALL hm_option_count('/ACTIV', nactiv)
717C-----
718 CALL hm_option_count('/ADMAS',nodmas)
719C-----
720 CALL hm_option_count('/CLUSTER',ncluster)
721C-----
722 CALL hm_option_count('/PYTHON_FUNCT' , npyfun)
723
724 CALL hm_option_count('/FUNCT' , nfunct0)
725
726 CALL hm_option_count('/TABLE/0', ntable0)
727 CALL hm_option_count('/TABLE/1', ntable1)
728 nfunct = nfunct0 + ntable0 + ntable1 + npyfun
729 ntable = nfunct
730C-----
731 CALL hm_option_count('/FUNC_2D', nfunc2d)
732C-----
733 iperturb = 0
734 CALL hm_option_count('/PERTURB',nperturb)
735 IF (nperturb > 0) iperturb = 1
736C-----
737 CALL hm_option_count('/SENSOR', nsensor)
738C-----
739 CALL hm_option_count('/SKEW',numskw)
740 CALL hm_option_count('/FRAME',numfram)
741C-----
742 CALL hm_option_count('TRANSFORM',ntransf)
743C---------------------------------------------------
744C LOADS
745C---------------------------------------------------
746 CALL hm_option_count('/CLOAD',nconld)
747 CALL hm_option_count('/PLOAD',npreld)
748 CALL hm_option_count('/LOAD/CENTRI',nloadc)
749 CALL hm_option_count('/LOAD/PFLUID',nloadp_f)
750 CALL hm_option_count('/LOAD/PBLAST',pblast%NLOADP_B)
751 CALL hm_option_count('/LOAD/PRESSURE',nloadp_hyd)
752C---------------------------------------------------
753C GRAVITY
754C---------------------------------------------------
755 CALL hm_option_count('/GRAV',ngrav)
756C---------------------------------------------------
757C BOUNDARY CONDITIONS
758C---------------------------------------------------
759 CALL hm_option_count('/BCS/LAGMUL',nbcslag)
760 CALL hm_option_count('/BCS/CYCLIC',nbcscyc)
761 CALL hm_option_count('/BCS',numbcs)
762 CALL hm_option_count('/NBCS',numbcsn)
763 nbcskin = numbcs - nbcslag
764 CALL hm_option_count('/BCS/WALL',bcs%NUM_WALL)
765C---------------------------------------------------
766C KINEMATIC CONDITIONS
767C---------------------------------------------------
768 CALL hm_option_count('/IMPDISP',nimpdisp)
769 CALL hm_option_count('/impvel' ,NIMPVEL)
770 CALL HM_OPTION_COUNT('/impacc' ,NIMPACC)
771 NFXVEL = NIMPDISP + NIMPVEL + NIMPACC
772C-----
773 CALL HM_OPTION_COUNT('/rbody',NRBODY)
774 NRBODY0 = NRBODY
775 CALL HM_OPTION_COUNT('/rbody/lagmul',NRBYLAG)
776 NRBYKIN = NRBODY - NRBYLAG
777 CALL HM_OPTION_COUNT('/fxbody',NFXBODY)
778C-----
779 CALL HM_OPTION_COUNT('/merge/rbody',NRBMERGE)
780C-----
781 CALL HM_OPTION_COUNT('/merge/node',NB_MERGE_NODE)
782C-----
783 CALL HM_OPTION_COUNT('/rbe2', NRBE2)
784 CALL HM_OPTION_COUNT('/rbe3', NRBE3)
785C-----
786 CALL HM_OPTION_COUNT('/rlink', NLINK)
787C-----
788 CALL HM_OPTION_COUNT('/rwall',NRWALL)
789C-----
790 CALL HM_OPTION_COUNT('/mpc',NUMMPC)
791 CALL HM_OPTION_COUNT('/cyl_joint', NJOINT)
792 CALL HM_OPTION_COUNT('/gjoint', NGJOINT)
793C---------------------------------------------------
794C SEATBELT TOOLS
795C---------------------------------------------------
796 CALL HM_OPTION_COUNT('/slipring', NSLIPRING)
797 CALL HM_OPTION_COUNT('/retractor', NRETRACTOR)
798C---------------------------------------------------
799C MONITORED VOLUMES
800C---------------------------------------------------
801 CALL HM_OPTION_COUNT('/monvol', NMONVOL)
802C---------------------------------------------------
803C INTERFACES
804C---------------------------------------------------
805 NINTER = 0
806 HM_NINTER= 0
807 CALL HM_OPTION_COUNT('/inter',HM_NINTER)
808C-----
809 CALL HM_OPTION_COUNT('/inter/type22',INT22)
810 IF(INT22>0)NSUB22=2
811 IF(INT22>0)ALEFVM_Param%IEnabled=1 !AUTOMATICALLY ENABLING ALE FVM SCHEME FOR FSI INTER22
812C-----
813 CALL HM_OPTION_COUNT('/inter/sub',NINTSUB)
814 NINTER = HM_NINTER - NINTSUB
815C-----
816 CALL HM_OPTION_COUNT('/friction',NINTERFRIC)
817 CALL HM_OPTION_COUNT('/fric_orient',NFRIC_ORIENT)
818C---------------------------------------------------
819 CALL HM_OPTION_COUNT('/damp',NDAMP)
820C---------------------------------------------------
821 CALL HM_OPTION_COUNT('/preload',NPRELOAD) !Bolt preloading
822C---------------------------------------------------
823 CALL HM_OPTION_COUNT('/sect',NSECT)
824C---------------------------------------------------
825C---------------------------------------------------
826C Box, Groups, Lines, Surfaces, Sets
827C---------------------------------------------------
828 CALL HM_OPTION_COUNT('/box' ,NBBOX)
829 CALL HM_OPTION_COUNT('/surf',NSURF)
830 CALL HM_OPTION_COUNT('/line' ,NSLIN)
831 CALL HM_OPTION_COUNT('/grnod' ,NGRNOD )
832 CALL HM_OPTION_COUNT('/grbric',NGRBRIC)
833 CALL HM_OPTION_COUNT('/grquad',NGRQUAD)
834 CALL HM_OPTION_COUNT('/grpart',NGRPART)
835 CALL HM_OPTION_COUNT('/grshel',NGRSHEL)
836 CALL HM_OPTION_COUNT('/grsh3n',NGRSH3N)
837 CALL HM_OPTION_COUNT('/grtria',NGRTRIA)
838 NGRSH3N = NGRSH3N + NGRTRIA ! 3D or 2D (same buffer)
839 CALL HM_OPTION_COUNT('/grtrus',NGRTRUS)
840 CALL HM_OPTION_COUNT('/grbeam',NGRBEAM)
841 CALL HM_OPTION_COUNT('/grspri',NGRSPRI)
842 NSETS = 0
843 CALL HM_OPTION_COUNT('/set',NSETS)
844 NGPE = NGRNOD + NGRBRIC + NGRQUAD + NGRSHEL + NGRSH3N + NGRTRUS + NGRBEAM + NGRSPRI + NGRPART
845C---------------------------------------------------
846C Initial conditions
847C---------------------------------------------------
848 CALL HM_OPTION_COUNT('/inivel',HM_NINVEL)
849 NINVEL = HM_NINVEL
850C-----
851 CALL HM_OPTION_COUNT('/refsta',IREFSTA)
852 IS_REFSTA = .FALSE.
853 IF(IREFSTA > 0)IS_REFSTA = .TRUE.
854C-----
855 CALL HM_OPTION_COUNT('/xref',NXREF)
856C-----
857 CALL HM_OPTION_COUNT('/eref',NEREF)
858C
859C---- READING OF /INISTA
860 CALL HM_READ_INISTA(S0FILE, ISIGI, IOUTP_FMT, IROOTYY_R, LSUBMODEL)
861C
862 IRFE=IRFORM/5
863 IRFL=IRFORM-5*IRFE
864 IRFORM=5*IRFE+IRFL
865 IF (IRFORM /= 12) THEN
866 CALL ANCMSG(MSGID=636,MSGTYPE=MSGWARNING,ANMODE=ANINFO_BLIND_1)
867 IRFORM = 12
868 ENDIF
869C
870C---- READING OF /REFSTA
871 CALL HM_READ_REFSTA(LSUBMODEL, XRFILE)
872C-----
873 CALL HM_OPTION_COUNT('/inigrav',NINIGRAV)
874C-----
875 CALL HM_OPTION_COUNT('/inimap1d', NINIMAP1D)
876 CALL HM_OPTION_COUNT('/inimap2d', NINIMAP2D)
877C-----
878 CALL HM_OPTION_COUNT('/inivol' ,NUM_INIVOL)
879C---------------------------------------------------
880C Thermal FE options
881C---------------------------------------------------
882 CALL HM_OPTION_COUNT('/initemp' ,GLOB_THERM%NINTEMP)
883 CALL HM_OPTION_COUNT('/imptemp' ,GLOB_THERM%NIMTEMP)
884 CALL HM_OPTION_COUNT('/impflux' ,GLOB_THERM%NIMPFLUX)
885 CALL HM_OPTION_COUNT('/convec' ,GLOB_THERM%NCONVEC)
886 CALL HM_OPTION_COUNT('/radiation',GLOB_THERM%NRADIA)
887C---------------------------------------------------
888C SPH
889C---------------------------------------------------
890 CALL HM_OPTION_COUNT('/sphcel',NUMSPH)
891
892 NSELSP = 0
893 CALL HM_OPTION_COUNT('/sphbcs',NSPCOND)
894 NSPHSYM= 0
895 MAXPJET= 0
896C
897c call NBSPH only once NSPMD is known in order to compute NSPHRES=NSPHRES*NSPMD to prepare NSPHRES by proc
898c for inlets treatment in SPMD
899c NSPHRES is the global number of SPH reserve
900C NBSPH computes NSPHIO & NSPHRES :
901 CALL NBSPH(LSUBMODEL)
902C
903C add NSPHRES number of SPH reserve for inlet
904 NUMNOD=NUMNOD+NSPHRES
905 NUMSPH=NUMSPH+NSPHRES
906C---------------------------------------------------
907C MADYMO
908C---------------------------------------------------
909 NEXMAD =0
910 NCONX =0
911 CALL HM_OPTION_COUNT('/madymo/exfem', NEXMAD)
912 CALL HM_OPTION_COUNT('/madymo/link', NCONX)
913 IF (NCONX > 0) CALL ANCMSG(MSGID=2023, MSGTYPE=MSGERROR, ANMODE=ANINFO)
914C---------------------------------------------------
915 NEIG = 0
916 CALL HM_OPTION_COUNT('/eig',NEIG)
917 IF (NEIG>0) IPARI0 = 0
918C---------------------------------------------------
919 NFLOW = 0
920 CALL HM_OPTION_COUNT('/bem/flow', NFLOW)
921 IDUMMY = 0
922 CALL HM_OPTION_COUNT('/bem/daa', IDUMMY)
923 NFLOW = NFLOW + IDUMMY
924C---------------------------------------------------
925C DFS, ALE, EULER
926C---------------------------------------------------
927 FLG_FSI = 0
928 ALE%GLOBAL%IS_BOUNDARY_MATERIAL = .FALSE.
929C-----
930 CALL HM_OPTION_COUNT('/ebcs',NEBCS)
931C---------------------------------------------------
932! Counting /DFS
933C---------------------------------------------------
934 CALL HM_OPTION_COUNT('/dfs/detpoin',DETONATORS%N_DET_POINT)
935 CALL HM_OPTION_COUNT('/dfs/detline',DETONATORS%N_DET_LINE)
936 CALL HM_OPTION_COUNT('/dfs/wav_sha',DETONATORS%N_DET_WAVE_SHAPER)
937 CALL HM_OPTION_COUNT('/dfs/detplan',DETONATORS%N_DET_PLANE)
938 CALL HM_OPTION_COUNT('/dfs/detcord',DETONATORS%N_DET_CORD)
939 DETONATORS%N_DET=DETONATORS%N_DET_POINT
940 . +DETONATORS%N_DET_LINE
941 . +DETONATORS%N_DET_WAVE_SHAPER
942 . +DETONATORS%N_DET_PLANE
943 . +DETONATORS%N_DET_CORD
944C---------------------------------------------------
945! /DFS/LASER
946C---------------------------------------------------
947 NLASER = 0
948 CALL HM_OPTION_COUNT('/dfs/laser',NLASER)
949C---------------------------------------------------
950 JEUL=0
951 CALL HM_OPTION_COUNT('/euler/mat', JEUL)
952 JEUL=MIN(JEUL,1)
953 LVEUL=32
954 IF(INTEG8==1)LVEUL=52
955C---------------------------------------------------
956 JALE=0
957 CALL HM_OPTION_COUNT('/ale/mat', JALE)
958 JALE=MIN(JALE,1)
959C---------------------------------------------------
960 NALEBCS = 0
961 CALL HM_OPTION_COUNT('/ale/bcs', NALEBCS)
962 IF(ALE%GLOBAL%ICAA == 1) ALE%GRID%NWALE = -1
963 DT_INPUT = ZERO
964 !********************************
965 ! /ALE/GRID: ALE grid formulation
966 !********************************
967 CALL HM_READ_ALE_GRID(DT_INPUT, ALE%GRID%ALPHA, ALE%GRID%GAMMA, ALE%GRID%VGX, ALE%GRID%VGY, ALE%GRID%VGZ,
968 . VOLMIN, LSUBMODEL, UNITAB)
969
970 !*******************************************************************
971 ! /ALE/MUSCL: activation of second order methods for LAW51 or LAW151
972 !*******************************************************************
973 CALL HM_READ_ALE_MUSCL(LSUBMODEL, UNITAB)
974
975 !************************
976 ! /ALE/LINK/VEL: counting
977 !************************
978 CALL HM_OPTION_COUNT('/ale/link/vel', NALELK)
979
980 !*****************
981 ! /ALE/SOLVER/FINT
982 !*****************
983 CALL HM_READ_ALE_SOLVER(LSUBMODEL, UNITAB, ALE%GLOBAL%ICAA, ALE%GLOBAL%ISFINT)
984
985 !********
986 ! /UPWIND
987 !********
988 CALL HM_READ_UPWIND(JUPW, ALE%UPWIND%UPWMG, ALE%UPWIND%UPWOG, ALE%UPWIND%UPWSM, LSUBMODEL, UNITAB)
989C---------------------------------------------------
990C OPTIONS DE L'ENGINE LUES DANS LE STARTER
991C { ENGINE OPTIONS READ IN THE STARTER }
992C---------------------------------------------------
993 CALL HM_OPTION_COUNT('/run',NRUN_ENG)
994 CALL HM_OPTION_COUNT('/anim',NANIM_ENG)
995 CALL HM_OPTION_COUNT('/tfile',NTFILE_ENG)
996 CALL HM_OPTION_COUNT('/rfile',NRFILE_ENG)
997 CALL HM_OPTION_COUNT('/dt',NDT_ENG)
998 CALL HM_OPTION_COUNT('/stop',NSTOP_ENG)
999 CALL HM_OPTION_COUNT('/vers',NVERS_ENG)
1000 NGINE=NRUN_ENG+NANIM_ENG+NTFILE_ENG+NRFILE_ENG+NDT_ENG+NSTOP_ENG+NVERS_ENG
1001C---------------------------------------------------
1002 NDSOLV=0
1003C----------------------------------------------------------
1004C Anim version 5 non encore supporte
1005 ANIM_VERS = 44
1006 NANIM_VERS = 0
1007 CALL HM_OPTION_COUNT('/anim/vers',NANIM_VERS)
1008 IF (NANIM_VERS > 0) THEN
1009 CALL HM_OPTION_START('/anim/vers')
1010 CALL HM_OPTION_NEXT()
1011 CALL HM_GET_INTV('anim_vers', ANIM_VERS, IS_AVAILABLE, LSUBMODEL)
1012 ENDIF
1013C=======================================================================
1014C Some initializations ... (some flags / values need a 1st pre-reading of some options)
1015C=======================================================================
1016 IFORM8 = 2
1017 DTINI = ZERO
1018 DTFAC = ZERO
1019 DTHIS = ZERO
1020 DO I = 1,10
1021 DTABF(I) = EP30
1022 DTABFWR(I) = EP30
1023 ENDDO
1024 DTMIN = ZERO
1025 IHSH = 0
1026 DO I = 1,9
1027 DTHIS1(I)= 0
1028 ENDDO
1029C----------------------------------------------------------
1030C solids with rotation
1031C----------------------------------------------------------
1032 KROT = 0
1033 CALL HM_OPTION_COUNT('/inivel/node',I)
1034 IF(I > 0)KROT=1
1035 IRODDL=MIN(1,NUMELC+NUMELP+NRBODY+NUMELR+NUMELTG+NGJOINT+NUMMPC+NFXBODY+NUMELX+KROT)
1036 IRODDL0 = 0
1037 IISROT = 0
1038 !for SIN initialization and IN allocation in lectur.F
1039 CALL HM_OPTION_START('material')
1040 DO I = 1, HM_NUMMAT
1041 MAT_NUMBER = I
1042 CALL HM_OPTION_READ_KEY(LSUBMODEL,KEYWORD2 = KEY)
1043 IF ((KEY(1:5) == 'law13.AND.'
1044 . KEY(1:6) /= 'law131.AND.' KEY(1:6) /= 'law132.AND.'
1045 . KEY(1:6) /= 'law133.AND.' KEY(1:6) /= 'law134.AND.'
1046 . KEY(1:6) /= 'law135.AND.' KEY(1:6) /= 'law136.AND.'
1047 . KEY(1:6) /= 'law137.AND.' KEY(1:6) /= 'law138.AND.'
1048 . KEY(1:6) /= 'law139' )
1049.OR. . KEY(1:5) == 'rigid') IRODDL0 = 1
1050 IF (KEY(1:5) == 'law68.OR.' KEY(1:5) == 'cosse') IISROT = 1
1051 ENDDO
1052C--------------------------------------------
1053C PRE-READ OF Prop IDS for triangles SH3N6 & solid rotations
1054C PRE-READ OF Prop IDS , PARTS & SOLIDSfor sol to SPH
1055C--------------------------------------------
1056 ALLOCATE(IGEO(NPROPGI*NUMGEO),STAT=stat)
1057 IF (STAT /= 0) CALL ANCMSG(MSGID=268,ANMODE=ANINFO,MSGTYPE=MSGERROR,C1='igeo')
1058 IGEO=0
1059 NUMELTG6 = 0
1060 NSPHSOL = 0
1061 ALLOCATE(NSUB (NUMGEO + NUMSTACK) ,STAT=stat)
1062 ALLOCATE(NISUB(NUMGEO + NUMSTACK) ,STAT=stat)
1063 ALLOCATE(NPLY (NUMGEO + NUMSTACK) ,STAT=stat)
1064c
1065 CALL HM_PREREAD_PROPERTIES(IGEO,NSPHSOL,NPLY,NSUB,NISUB,LSUBMODEL,DEFAULTS)
1066
1067.AND. IF(NSPHSOL/=0NUMELS8/=0)THEN
1068 ALLOCATE(IPART(LIPART1*NPART),STAT=stat)
1069 IF (STAT /= 0) CALL ANCMSG(MSGID=268,ANMODE=ANINFO,MSGTYPE=MSGERROR,C1='ipart')
1070 IPART=0
1071 CALL HM_PREREAD_PART(IPART,IGEO,LSUBMODEL)
1072C
1073 ALLOCATE(IXS(NIXS*NUMELS8),STAT=stat)
1074 IF (STAT /= 0) CALL ANCMSG(MSGID=268,ANMODE=ANINFO,MSGTYPE=MSGERROR,C1='ipart')
1075 CALL HM_PRELCE16S(IPART,IGEO,IXS,NSPHSOL,LSUBMODEL,IS_DYNA)
1076C
1077 IREST_MSELT=1
1078C
1079 ENDIF
1080 NUMNOD=NUMNOD+NSPHSOL
1081 NUMSPH=NUMSPH+NSPHSOL
1082C-----
1083 IF(NUMSPH/=0)THEN
1084 CALL HM_READ_SPHGLO(LSUBMODEL)
1085 ENDIF
1086C----------------------------------------------------------
1087 CALL CONTRBE2(ICR,LSUBMODEL)
1088 CALL CONTRBE3(ICR3,LSUBMODEL)
1089.OR..OR. IF(IISROT==1ICR>0ICR3>0)IRODDL = 1
1090C----------------------------------------------------------
1091 CALL HM_OPTION_COUNT('/userwi',USER_WINDOWS%HAS_USER_WINDOW)
1092 NCPRI=1
1093 IF(DTINI==ZERO)DTINI=EP06
1094 IF(VOLMIN==ZERO)VOLMIN=-EP20
1095 DT2OLD=DTINI/ONEP1
1096 TT=ZERO
1097 DT1=ZERO
1098 DT2=ZERO
1099C
1100 T1S= ZERO
1101 DT2S=ZERO
1102C
1103 DTX=ZERO
1104 TSTOP=ZERO
1105 DTANIM = ZERO
1106 TANIM=ZERO
1107 TANIM_STOP = EP20
1108 DTOUTP = ZERO
1109 TOUTP=ZERO
1110 THIS=ZERO
1111 TABFIS=ZERO
1112 TABFWR=ZERO
1113 ECONTV = ZERO
1114 OUTPUT%TH%WFEXT = ZERO
1115 REINT = ZERO
1116 UREINT = ZERO
1117 ECONTD = ZERO
1118 ECONT_CUMU = ZERO
1119
1120 IRUN=0
1121 IGER=0
1122 IANIM=0
1123 IH3D=0
1124C
1125 DO I = 1,9
1126 THIS1(I)= 0
1127 ENDDO
1128C
1129 N2D = NANALY
1130.AND. IF(DTFAC==ZERON2D/=0)DTFAC=0.67 !0.67 => 0.670000016689301 !ZEP67 => 0.670000000000000
1131.AND. IF(DTFAC==ZERON2D==0)DTFAC=0.90 !0.90 => 0.899999976158142 !ZEP9 => 0.900000000000000
1132C=======================================================================
1133C Opening of some input files (.sty, etc...)
1134C=======================================================================
1135.OR. IF(ISIGI==1ISIGI==2) THEN
1136C-- Lecture fichier S00 supprimee
1137.OR..OR. ELSEIF(ISIGI==3ISIGI==4ISIGI==5) THEN
1138C-- Lecture fichier Y00/Ynn
1139 IF(S0FILE==' ') THEN
1140 ISIGI=-ISIGI
1141 ELSE
1142 J = 0
1143 IF(IROOTYY_R==2)THEN
1144 DO I=1,ncharline
1145 IF(S0FILE(I:I)/=' ')J = J + 1
1146 ENDDO
1147 N = J - 3
1148 TMP_NAME=INFILE_NAME(1:INFILE_NAME_LEN)//S0FILE(1:len_trim(S0FILE))
1149 LEN_TMP_NAME = INFILE_NAME_LEN+len_trim(S0FILE)
1150 IF(S0FILE(N:N)=='y')THEN
1151 OPEN(UNIT=IIN4,FILE=TMP_NAME(1:LEN_TMP_NAME),ACCESS='sequential',FORM='formatted',STATUS='old')
1152
1153 TMP_NAME=INFILE_NAME(1:INFILE_NAME_LEN)//S0FILE(1:N)//'000'
1154 LEN_TMP_NAME = INFILE_NAME_LEN+N+3
1155 OPEN(UNIT=IIN5,FILE=TMP_NAME(1:LEN_TMP_NAME),ACCESS='sequential',FORM='formatted',STATUS='old')
1156 ELSE
1157 OPEN(UNIT=IIN4,FILE=TMP_NAME(1:LEN_TMP_NAME),ERR=100,ACCESS='sequential',FORM='formatted',STATUS='old')
1158 CALL ANCMSG(MSGID=169,MSGTYPE=MSGERROR,ANMODE=ANINFO,C1=S0FILE)
1159100 CALL ANCMSG(MSGID=2062,MSGTYPE=MSGERROR,ANMODE=ANINFO,C1=S0FILE)
1160 ENDIF
1161 ELSE
1162 J = 0
1163 I = 1
1164 DO WHILE(S0FILE(I:I)/=' ')
1165 J = J + 1
1166 I = I + 1
1167 ENDDO
1168 N = J-3
1169 TMP_NAME=INFILE_NAME(1:INFILE_NAME_LEN)//S0FILE(1:len_trim(S0FILE))
1170 LEN_TMP_NAME = INFILE_NAME_LEN+len_trim(S0FILE)
1171 IF(S0FILE(N:N+4)=='.sty')THEN
1172 OPEN(UNIT=IIN4,FILE=TMP_NAME(1:LEN_TMP_NAME),ACCESS='sequential',FORM='formatted',STATUS='old')
1173
1174 TMP_NAME=INFILE_NAME(1:INFILE_NAME_LEN)//S0FILE(1:N-6)//'_0000.sty'
1175 LEN_TMP_NAME = INFILE_NAME_LEN+N-6+9
1176 OPEN(UNIT=IIN5,FILE=TMP_NAME(1:LEN_TMP_NAME),ACCESS='sequential',FORM='formatted',STATUS='old')
1177 ELSE
1178 OPEN(UNIT=IIN4,FILE=TMP_NAME(1:LEN_TMP_NAME),ERR=200,ACCESS='sequential',FORM='formatted',STATUS='old')
1179 CALL ANCMSG(MSGID=169,MSGTYPE=MSGERROR,ANMODE=ANINFO,C1=S0FILE)
1180200 CALL ANCMSG(MSGID=2062,MSGTYPE=MSGERROR,ANMODE=ANINFO,C1=S0FILE)
1181 ENDIF
1182 ENDIF
1183 ENDIF
1184 ENDIF
1185C---------------------------------------------------
1186C COUNT /INI CARDS READ BY HM_READER
1187C---------------------------------------------------
1188 CALL HM_OPTION_COUNT('/inishe', NB_INISHE)
1189 CALL HM_OPTION_COUNT('/inish3', NB_INISH3)
1190 CALL HM_OPTION_COUNT('/inibri', NB_INIBRI)
1191 CALL HM_OPTION_COUNT('/iniqua', NB_INIQUAD)
1192 CALL HM_OPTION_COUNT('/inibeam', NB_INIBEAM)
1193 CALL HM_OPTION_COUNT('/initruss', NB_INITRUSS)
1194 CALL HM_OPTION_COUNT('/inispri', NB_INISPRIG)
1195 CALL HM_OPTION_COUNT('/inisphcel', NB_INISPHCEL)
1196
1197.AND. IF(ISIGI==0
1198 . (NB_INISHE+NB_INISH3
1199 . +NB_INIBRI+NB_INIQUAD
1200 . +NB_INIBEAM+NB_INITRUSS
1201 . +NB_INISPRIG+NB_INISPHCEL)/=0) ISIGI=-3
1202C=======================================================================
1203C Writing to 0000.out FILE
1204C=======================================================================
1205 WRITE (IOUT,'(a)') TRIM(NAMES_AND_TITLES%TITLE)
1206 WRITE (IOUT,'(//a/a/,(a,1pg20.13))')TITRE(110),TITRE(111)
1207 WRITE (IOUT,'(a,i10)') TITRE(117)(1:57),IPRI,
1208 . TITRE(120)(1:57),INVERS_SRC
1209 WRITE (IOUT,'(a,i10)')
1210 . ' nsubs: number of subsets. . . . . . . . . . . . . .',
1211 . NSUBS,
1212 . ' npart: number of parts. . . . . . . . . . . . . . .',
1213 . NPART
1214 WRITE (IOUT,'(a,i10)') TITRE(125)(1:57),NUMMAT-1,
1215 . TITRE(135)(1:57),NUMGEO,
1216 . TITRE(126)(1:57),NUMNOD,TITRE(127)(1:57),NUMSKW,
1217 . TITRE(128)(1:57),NUMBCS
1218 IF (NALEBCS /= 0)
1219 . WRITE (IOUT,'(a,i10)') TITRE(189)(1:57), NALEBCS
1220 WRITE (IOUT,'(a,i10)')
1221 . ' numfram: number of reference frames . . . . . . . . .',
1222 . numfram
1223 WRITE (iout,'(A,I10)')
1224 . titre(131)(1:57),numelq ,
1225 . titre(191)(1:57),numeltria,
1226 . titre(132)(1:57),numels ,
1227 . titre(133)(1:57),numelc ,
1228 . titre(134)(1:57),numelt ,
1229 . titre(136)(1:57),numelp ,
1230 . titre(137)(1:57),numelr ,
1231 . titre(180)(1:57),numeltg-numeltria,
1232 . titre(182)(1:57),numelx ,
1233 . titre(186)(1:57),numelig3d
1234 WRITE (iout,'(A,I10)')
1235 .' NUMSPH : NUMBER OF SMOOTH PARTICLES (SPH CELLS) . . .',
1236 . numsph
1237 WRITE (iout,'(A,I10)')
1238 .' NSPHBCS: NUMBER OF SPH SYMMETRY CONDITIONS. . . . . .',
1239 . nspcond
1240 WRITE (iout,'(A,I10)')
1241 .' NSPHIO : NUMBER OF SPH INLET/OUTLET CONDITIONS. . . .',
1242 . nsphio
1243 IF(nsphres/=0)THEN
1244 WRITE (iout,'(A,/,A,I10)')
1245 .' NSPHRES:NUMBER OF PARTICLES FROM SPH RESERVES AMONG NUMSPH,',
1246 .' NUMBER OF NODES FROM SPH RESERVES AMONG NUMNOD. . .',
1247 . nsphres
1248 ENDIF
1249 IF(nsphsol/=0)THEN
1250 WRITE (iout,'(A,/,A,I10)')
1251 .' NSPHSOL:NUMBER OF PARTICLES & NODES CREATED FROM SOLIDS .',
1252 .' (AMONG NUMSPH & NUMNOD) . . . . . . . . . . . . .',
1253 . nsphsol
1254 ENDIF
1255 WRITE (iout,'(/(A,I10))') titre(138)(1:57),nfunct,
1256 .' NGRAV: NUMBER OF GRAVITY LOADS . . . . . . . . . . .',
1257 .ngrav,
1258 .' NFUNC2D: NUMBER OF USER 2D FUNCTIONS . . . . . . . .',
1259 .nfunc2d,
1260 .' NINIGRV:NUMBER OF INITIAL GRAVITY LOADS . . . . . . .',
1261 .ninigrav,
1262 .' NCONLD: NUMBER OF CONCENTRATED LOADS. . . . . . . . .',
1263 .nconld,
1264 .' NINVEL: NUMBER OF INITIAL VELOCITIES. . . . . . . . .',
1265 .ninvel,
1266 .' NPRELD: NUMBER OF PRESSURE LOADS. . . . . . . . . . .',
1267 .npreld,
1268 .' NINIMAP1D: NUMBER OF INITIAL 1D MAPPING.. . . . . . .',
1269 .ninimap1d,
1270 .' NINIMAP2D: NUMBER OF INITIAL 2D MAPPING.. . . . . . .',
1271 .ninimap2d
1272 IF(npreload>0) THEN
1273 WRITE (iout,'(A,I10)')
1274 . ' NPRELOAD: NUMBER OF BOLT PRELOADINGS. . . . . . . . .',
1275 . npreload
1276 ENDIF
1277 IF(detonators%N_DET > 0)THEN
1278 WRITE (iout,'(/(A,I10))') titre(141)(1:57),detonators%N_DET_POINT,
1279 . titre(171)(1:57),detonators%N_DET_LINE,
1280 . titre(172)(1:57),detonators%N_DET_WAVE_SHAPER,
1281 . titre(187)(1:57),detonators%N_DET_CORD,
1282 . titre(188)(1:57),detonators%N_DET_PLANE
1283 ENDIF
1284 IF(nlaser > 0) WRITE (iout,'((A,I10))') titre(178)(1:57),nlaser
1285 WRITE (iout,'(4X,A,I10)')
1286 . 'NUMBER OF ACCELEROMETERS. . . . . . . . . . . . . . .',
1287 . naccelm,
1288 . 'NUMBER OF SENSORS . . . . . . . . . . . . . . . . . .',
1289 . nsensor,
1290 . 'NUMBER OF GAUGES. . . . . . . . . . . . . . . . . . .',
1291 . nbgauge
1292 WRITE (iout,'(A,I10)') titre(146)(1:57),ninter,
1293 . titre(147)(1:57),nrwall,titre(148)(1:57),nrbody,
1294 .' NFXBODY: NUMBER OF FLEXIBLE BODIES. . . . . . . . . .',
1295 .nfxbody,
1296 . titre(149)(1:57),nconx,
1297 . titre(150)(1:57),nodmas,titre(183)(1:57),nimpdisp,
1298 . titre(184)(1:57),nimpvel,titre(185)(1:57),nimpacc,
1299 . titre(152)(1:57),nrivet,titre(153)(1:57),nsect,
1300 . titre(155)(1:57),njoint
1301C
1302 IF(ninterfric > 0) WRITE (iout,'(A,I10)') titre(190)(1:57),ninterfric
1303 IF(nalelk>0)WRITE(iout,'(A,I10)')
1304 .' NALELK: NUMBER OF ALE LINKS. . . . . . . . . . . . . ',
1305 .nalelk
1306 WRITE (iout, 5051) nactiv
1307 WRITE (iout, 5052) ndamp
1308 WRITE (iout, 5053) ngjoint
1309 WRITE (iout, 5054) nummpc
1310 WRITE (iout, 5050) nr2rlnk
1311 WRITE (iout, 5055) nsubdom
1312 WRITE (iout,5000)nvolu+nmonvol
1313 WRITE(iout,'(A,I10)')
1314 .' neig: number of eigen and static modes problems . .',
1315 .NEIG
1316 WRITE(IOUT,'(a,i10)')
1317 .' nbem: number of bem solved problems . . . . . . . .',
1318 .NFLOW
1319 WRITE(IOUT,'(a,i10)')
1320 .' nrbe2: number of rbe2 rigid elements . . . . . . . .',
1321 .NRBE2
1322 WRITE(IOUT,'(a,i10)')
1323 .' nrbe3: number of rbe3 constraint elements . . . . .',
1324 .NRBE3
1325 WRITE (IOUT,'(2a/a,i10)')' initial stress file =',TRIM(S0FILE),
1326 . ' flag isigi. . .',
1327 . ISIGI
1328C
1329 IF (IREFSTA/=0) THEN
1330 WRITE (IOUT,'(2a/a,i10/a,i10)')
1331 . ' reference metric file =',XRFILE,
1332 . ' flag irefsta. .',
1333 . IREFSTA,
1334 . ' number of steps , nitrs. .',
1335 . NITRS
1336 ENDIF
1337
1338 IF(JALE+JEUL/=0)THEN
1339 WRITE(IOUT,'(a,i10)')' isfint: ale/euler momentum integration formulation. .',ALE%GLOBAL%ISFINT
1340 ENDIF
1341 IF(ALE%GLOBAL%ICAA==1)THEN
1342 WRITE(IOUT,5380)
1343 END IF
1344 IF(JALE/=0)THEN
1345 WRITE (IOUT,5100)ALE%GRID%NWALE
1346 SELECT CASE (ALE%GRID%NWALE)
1347 CASE(0);WRITE (IOUT,5199) ALE%GRID%ALPHA,ALE%GRID%GAMMA,ALE%GRID%VGX,ALE%GRID%VGY,ALE%GRID%VGZ,VOLMIN
1348 CASE(1);WRITE (IOUT,5200) ALE%GRID%ALPHA,VOLMIN
1349 CASE(2);WRITE (IOUT,5300) DT_INPUT, ALE%GRID%ALPHA,ALE%GRID%GAMMA,ALE%GRID%VGX,ALE%GRID%VGY,VOLMIN
1350 CASE(3);WRITE (IOUT,5350)
1351 CASE(4);WRITE (IOUT,5351) ALE%GRID%ALPHA,ALE%GRID%GAMMA,ALE%GRID%VGX,ALE%GRID%VGY
1352 CASE(5);WRITE (IOUT,5353) ALE%GRID%ALPHA,NINT(ALE%GRID%VGX)
1353 CASE(6);WRITE (IOUT,5354)
1354 CASE(7);
1355 LABEL_DEF = ' no'
1356 IF(INT(ALE%GRID%VGX) == 1)LABEL_DEF = 'yes'
1357 LABEL_ROT = ' no'
1358 IF(INT(ALE%GRID%VGY) == 1)LABEL_ROT = 'yes'
1359 WRITE (IOUT,5355) LABEL_DEF,LABEL_ROT,ALE%GRID%ALPHA,ALE%GRID%GAMMA
1360 END SELECT
1361 ENDIF
1362
1363 IF(ALEFVM_Param%ISOLVER /= 0)ALEMUSCL_Param%IALEMUSCL=0 !muscl not compatible with FVM solver for int22
1364 CALL HM_OPTION_COUNT('/ale/muscl', NALEMUSCL)
1365 IF(ALE%UPWIND%UPWSM/=ONE)THEN
1366.OR. IF(NALEMUSCL>0 ALEMUSCL_Param%IALEMUSCL>0)THEN
1367 NALEMUSCL=0
1368 ALEMUSCL_Param%IALEMUSCL=0
1369 !ignore muscl & use upwsm3 (backward compatibility)
1370 !double check engine file with /ALE/MUSCL/OFF
1371 ENDIF
1372 ENDIF
1373
1374 IF(JALE+JEUL/=0)THEN
1375.AND. IF(ALEMUSCL_Param%IALEMUSCL == 1 JUPW/=0 ) THEN !eta3 + muscl(law51) : conflict
1376 CALL ANCMSG(MSGID=1564, MSGTYPE=MSGWARNING, ANMODE=ANINFO)
1377.OR. ELSEIF(ALEMUSCL_Param%IALEMUSCL == 0 JUPW/=0)THEN
1378 WRITE (IOUT,5360)ALE%UPWIND%UPWMG,ALE%UPWIND%UPWOG !,UPWSM
1379 END IF
1380.OR. IF (ALEMUSCL_Param%IALEMUSCL == 1 ALEMUSCL_Param%IALEMUSCL==2) THEN
1381 WRITE (IOUT, 5504) ALEMUSCL_Param%BETA,ALEMUSCL_Param%IALEMUSCL-1
1382 ENDIF
1383 ENDIF
1384
1385 WRITE (IOUT,'(//,a,i10)') TITRE(130)(1:57),NANALY
1386
1387! for print out
1388 IHBE_DS= DEFAULTS%SOLID%ISOLID
1389 ISST_DS= DEFAULTS%SOLID%ISMSTR
1390 ICPRE_D= DEFAULTS%SOLID%ICPRE
1391 ITET4_D= DEFAULTS%SOLID%ITETRA4
1392 ITET10_D= DEFAULTS%SOLID%ITETRA10
1393 IFRAME_DS = DEFAULTS%SOLID%IFRAME
1394 IMAS_DS= DEFAULTS%SOLID%IMAS
1395!
1396 IOFFSET_D= DEFAULTS%SHELL%ioffset !< offset support contact
1397 IHBE_D = DEFAULTS%SHELL%ishell
1398 ISH3N_D= DEFAULTS%SHELL%ish3n
1399 ISST_D = DEFAULTS%SHELL%ismstr
1400 IPLA_D = DEFAULTS%SHELL%iplas
1401 ITHK_D = DEFAULTS%SHELL%ithick
1402 IDRIL_D= DEFAULTS%SHELL%idrill
1403
1404 ICONTROL_D= DEFAULTS%SOLID%ICONTROL
1405 IF(N2D==0)THEN
1406.AND. IF(NUMELQ/=0. OR. (NUMELTG/=0N2D/=0) )THEN
1407 CALL ANCMSG(MSGID = 286,MSGTYPE = MSGERROR,ANMODE = ANINFO)
1408 ENDIF
1409 ENDIF
1410C 012
1411 WRITE(IOUT,5500)INTEG8,IPARITH,USER_WINDOWS%HAS_USER_WINDOW,IHBE_DS,
1412 . ITET4_D,ITET10_D,ISST_DS,ICPRE_D,IHBE_D,ISST_D,ITHK_D,
1413 . IPLA_D,ISTR_D,ISHEA_D,INSH,ISH3N_D, NPTS_D, IFRAME_DS,
1414 . IOFFSET_D,
1415 . ICONTROL_D
1416 WRITE(IOUT,5501)ISHFRAM
1417 IF(ISH3NFRAM>0)WRITE(IOUT,5901)ISH3NFRAM
1418.AND. IF(JALE/=0 ALEFVM_Param%IEnabled/=0)WRITE(IOUT,5502)ALEFVM_Param%ISOLVER
1419.AND. IF(JALE/=0 ALE%GLOBAL%ISFINT/=0) WRITE(IOUT,5503)ALE%GLOBAL%ISFINT
1420 WRITE(IOUT,5700)LAGM_TOL
1421 WRITE(IOUT,5800)IMAS_DS
1422.AND. IF ((NUMSPH>0)(NSPHSOL==0)) WRITE (IOUT,5600) SPASORT,LVOISPH,KVOISPH
1423.AND. IF ((NUMSPH>0)(NSPHSOL>0)) THEN
1424 IF (ITSOL2SPH==1) WRITE (IOUT,5610) SPASORT,LVOISPH,KVOISPH,ITSOL2SPH
1425 IF (ITSOL2SPH==2) WRITE (IOUT,5620) SPASORT,LVOISPH,KVOISPH,ITSOL2SPH
1426 ENDIF
1427C-------------------------------------------------------------------
1428 WRITE(IOUT,1000)
1429 IF(IEXPM==1) THEN
1430 ELSE IF(IMOT/=0) THEN
1431 WRITE(IOUT,*)'memory request not efficient on this computer'
1432 ENDIF
1433 IF(ICRAY<2) THEN
1434C
1435 IPCT=NINT(REEL*100)
1436c WRITE(IOUT,3000) LMA,LAM,IPCT
1437 ELSE IF(ICRAY==2) THEN
1438C MasPar Special
1439c WRITE(IOUT,3002) LMA,LAM
1440 ENDIF
1441
1442 IF(NTHREAD_S>1)THEN
1443 STR_NBTHREADS = 'threads'
1444 ELSE
1445 STR_NBTHREADS = 'thread'
1446 ENDIF
1447 IF(GOT_INSPIRE_ALM == 1)THEN
1448 WRITE(IOUT,4001) NTHREAD_S,STR_NBTHREADS,NSPMD,NTHREAD,IARCH,ARCHN(IARCH),NVSIZ
1449 ELSE
1450 WRITE(IOUT,4000) NTHREAD_S,STR_NBTHREADS,NSPMD,NTHREAD,IARCH,ARCHN(IARCH),NVSIZ
1451 ENDIF
1452C
1453.OR. IF(IGNORE_SPMD==1IGNORE_THREADS==1)THEN
1454 WRITE(IOUT,'(a)') ' '
1455 WRITE(IOUT,'(a)') ' info : number of spmd domain and threads defined in command line'
1456 WRITE(IOUT,'(a)') ' /spmd card PARAMETER are ignored'
1457 ENDIF
1458 IF(GOT_VARIABLE==1)THEN
1459 WRITE(IOUT,'(a)') ' '
1460 WRITE(IOUT,*) ' info : omp_num_threads set, default setting ignored'
1461 WRITE(IOUT,*)' '
1462 ENDIF
1463
1464 IF(OLD_RSB==1)THEN
1465 WRITE(IOUT,'(a)') ' '
1466 WRITE(IOUT,'(a)') ' info : /spmd option, dectyp=2 '
1467 WRITE(IOUT,'(a)') ' rsb domain decomposition deprecated, changing to default VALUE'
1468 ENDIF
1469
1470
1471 IF (USERL_AVAIL==1)THEN
1472 WRITE(IOUT,4500)
1473 IF(GOT_INSPIRE_ALM == 1)THEN
1474 WRITE(IOUT,4601) DLIBFILE(1:DLIBFILE_SIZE),DLIBTKVERS
1475 ELSE
1476 WRITE(IOUT,4600) DLIBFILE(1:DLIBFILE_SIZE),DLIBTKVERS
1477 ENDIF
1478 ENDIF
1479C
1480C--------------------------------------------------------------------
1481C DEALLOCATE
1482C--------------------------------------------------------------------
1483 IF(ALLOCATED(IXS)) DEALLOCATE(IXS)
1484 IF(ALLOCATED(IGEO)) DEALLOCATE(IGEO)
1485 IF(ALLOCATED(IPART))DEALLOCATE(IPART)
1486C--------------------------------------------------------------------
1487 1000 FORMAT(//
1488 & 4X,'speed parameters '/
1489 & 4X,'---------------- '/)
1490c 2000 FORMAT(
1491c & 4X,'MEMORY REQUESTED BY USER (KWORDS). . . . . . . . . .',I10)
1492c 3000 FORMAT(
1493c & 4X,'MEMORY AVAILABLE FOR INTEGERS . . . . . . . . . . . .',I10/
1494c & 4X,'MEMORY AVAILABLE FOR REALS. . . . . . . . . . . . . .',I10/
1495c & 4X,'PERCENTAGE OF MEMORY FOR REALS . . . . . . . . . . .',I10)
1496c 3002 FORMAT(
1497c & 4X,'MEMORY AVAILABLE FOR INTEGERS . . . . . . . . . . . .',I10/
1498c & 4X,'MEMORY AVAILABLE FOR REALS. . . . . . . . . . . . . .',I10)
1499 4000 FORMAT(
1500 & 4X,'starter running on. . . . . . . . . . . . . . . . . .',I10,
1501 . ' ',A20/
1502 & 4X,'number of spmd domains. . . . . . . . . . . . . . . .',I10/
1503 & 4X,'number of threads per domain. . . . . . . . . . . . .',I10/
1504 & 4X,'architecture optimization . . . . . . . . . . . . . .',I10,
1505 . ', ',A20/
1506 & 4X,'SIZE of element buffer. . . . . . . . . . . . . . . .',I10)
1507 4001 FORMAT(
1508 & 4X,'solver running on . . . . . . . . . . . . . . . . . .',I10,
1509 . ' ',A20/
1510 & 4X,'number of spmd domains. . . . . . . . . . . . . . . .',I10/
1511 & 4X,'number of threads per domain. . . . . . . . . . . . .',I10/
1512 & 4X,'architecture optimization . . . . . . . . . . . . . .',I10,
1513 . ', ',A20/
1514 & 4X,'SIZE of element buffer. . . . . . . . . . . . . . . .',I10)
1515
1516 4500 FORMAT(//
1517 & 4X,'EXTERNAL library for users code INTERFACE '/
1518 & 4X,'----------------------------------------- '/)
1519 4600 FORMAT(
1520 & 4X,'library name . . . . . . . . . . . . . . . . . . . . ',A/
1521 & 4X,'radioss users code INTERFACE version . . . . . . . .',I10)
1522 4601 FORMAT(
1523 & 4X,'library name . . . . . . . . . . . . . . . . . . . . ',A/
1524 & 4X,'solver users code INTERFACE version . . . . . . . . .',I10)
1525
1526
1527 5000 FORMAT(
1528 & ' nvolu: number of monitored volumes . . . . . . . . .',I10)
1529 5050 FORMAT(
1530 & ' nr2rlnk: number of EXTERNAL coupling links . . . . . ',I10)
1531 5051 FORMAT(
1532 & ' nactiv: number of element deactivation groups . . . .',I10)
1533 5052 FORMAT(
1534 & ' ndamp: number of rayleigh damping groups . . . . . .',I10)
1535 5053 FORMAT(
1536 & ' ngjoint: number of gear TYPE joints . . . . . . . . .',I10)
1537 5054 FORMAT(
1538 & ' nummpc: number of multi-point constraints . . . . . .',I10)
1539 5055 FORMAT(
1540 & ' nsubdom: number of subdomains . . . . . .. . . . . . ',I10)
1541 5100 FORMAT(/
1542 & 4X,'nwale : ale grid velocity formulation . . . . . . . .',I10)
1543 5199 FORMAT(//
1544 & 4X,'ale grid smoothing formulation'/
1545 & 4X,'------------------------------'/
1546 & 5X,'donea grid velocity computation method '//
1547 & 5X,'alpha : donea coefficient. . . . . . . . . . ',1PG20.13/
1548 & 5X,'gamma : grid velocity limitation factor. . . ',1PG20.13/
1549 & 5X,'fscalex : x-grid velocity scale factor . . . . ',1PG20.13/
1550 & 5X,'fscaley : y-grid velocity scale factor . . . . ',1PG20.13/
1551 & 5X,'fscalez : z-grid velocity scale factor . . . . ',1PG20.13/
1552 & 5X,'volmin : minimum volume for element deletion. ',1PG20.13)
1553 5200 FORMAT(//
1554 & 4X,'ale grid smoothing formulation'/
1555 & 4X,'------------------------------'/
1556 & 5X,'altair average displacement grid formulation '//
1557 & 5X,'umax : maximum absolute grid velocity . . . . ',1PG20.13/
1558 & 5X,'vmin : minimum volume for element deletion. . ',1PG20.13)
1559 5300 FORMAT(//
1560 & 4X,'ale grid smoothing formulation'/
1561 & 4X,'------------------------------'/
1562 &5X,'altair spring method for grid velocity computation '//
1563 &5X,'dt0 : typical time step . . . . . . . . . . . . ',1PG20.13/
1564 &5X,'dt0* : effective time step. . . . . . . . . . . . ',1PG20.13/
1565 &5X,'gamma : non linearity factor . . . . . . . . . . . ',1PG20.13/
1566 &5X,'eta : damping coefficient . . . . . . . . . . . ',1PG20.13/
1567 &5X,'nu : shear factor . . . . . . . . . . . . . . . ',1PG20.13/
1568 &5X,'volmin: minimum volume for element deletion. . . . ',1PG20.13)
1569 5350 FORMAT(//
1570 & 4X,'ale grid smoothing formulation'/
1571 & 4X,'------------------------------'/
1572 & 5X,'grid velocity is not computed(quasi euler) ')
1573 5351 FORMAT(//
1574 & 4X,'ale grid smoothing formulation'/
1575 & 4X,'------------------------------'/
1576 & 5X,'altair standard method for grid velocity computation '//
1577 & 5X,'alpha : stability factor . . . . . . . . . . . . ',1PG20.13/
1578 & 5X,'gamma : non linearity factor . . . . . . . . . . ',1PG20.13/
1579 & 5X,'beta : damping coefficient. . . . . . . . . . . ',1PG20.13/
1580 & 5X,'lc : characteristic length. . . . . . . . . . ',1PG20.13)
1581 5353 FORMAT(//
1582 & 4X,'ale grid smoothing formulation'/
1583 & 4X,'------------------------------'/
1584 & 5X,'laplacian smoothing '//
1585 & 5X,'lambda:. . . . . . . . . . . . . . . . . . . . . ',1PG20.13/
1586 & 5X,'niter :. . . . . . . . . . . . . . . . . . . . . ',I10)
1587 5354 FORMAT(//
1588 & 4X,'ale grid smoothing formulation'/
1589 & 4X,'------------------------------'/
1590 & 5X,'volume smoothing ')
1591 5355 FORMAT(//
1592 & 4X,'ale grid smoothing formulation'/
1593 & 4X,'------------------------------'/
1594 & 5X,'flow-tracking(mass weighted averaged velocity)'//
1595 & 5X,'deformation enabled : . . . . . . . . . . . . . ',A3/
1596 & 5X,'rotation enabled : . . . . . . . . . . . . . . . ',A3/
1597 & 5X,'deformation scale factor : . . . . . . . . . . . ',1PG20.13/
1598 & 5X,'rotation scale factor: . . . . . . . . . . . . . ',1PG20.13/)
1599
1600
1601 5360 FORMAT(//
1602 & 4X,'ale upwind parameters'/
1603 & 4X,'---------------------'/
1604 & 5X,'upwind for momentum transport . . . . .=',1PG20.13/,
1605 & 5X,'upwind for other transport. . . . . . .=',1PG20.13/)
1606 5380 FORMAT(/
1607 & 4X,'caa fluid formulation activated')
1608 5400 FORMAT(
1609 & 4X,'alpha m + beta k damping'/,
1610 & 4X,'------------------------'/,
1611 . 4X,'node group id(=0 all nodes). . . . . . ',I5/,
1612 & 5X,'alpha . . . . . . . . . . . . . . . . .=',1PG20.13/,
1613 & 5X,'beta. . . . . . . . . . . . . . . . . .=',1PG20.13)
1614 5500 FORMAT(//4X,'analysis options'/
1615 & 4X,'----------------'//
1616 & 4X,'integ8 : 8 gauss point condensed integration . . . . .',I10/
1617 & 4X,'iparith: parallel arithmetic flag(2 off, 1 0n). . . .',I10/
1618 & 4X,'iuserw : general user window flag. . . . . . . . . . .',I10/
1619 & 4X,'isolid : default brick formulation flag. . . . . . . .',I10/
1620 & 4X,'itet4 : default tetra4 formulation flag . . . . . . .',I10/
1621 & 4X,'itet10 : default tetra10 formulation flag. . . . . . .',I10/
1622 & 4X,'ismstr : default brick small strain flag . . . . . . .',I10/
1623 & 4X,'icpre : default solid constant pressure flag. . . . .',I10/
1624 & 4X,'ishell : default shell formulation flag. . . . . . . .',I10/
1625 & 4X,'isst : default shell small strain flag . . . . . . .',I10/
1626 & 4X,'ithk : default shell thickness variation flag. . . .',I10/
1627 & 4X,'ipla : default plane stress plasticity flag. . . . .',I10/
1628 & 4X,'istr : default shell strain computation flag . . . .',I10/
1629 & 4X,'ishea : default shell shear computation flag. . . . .',I10/
1630 & 4X,'insh : shell inertia flag. . . . . . . . . . . . . .',I10/
1631 & 4X,'ish3n : default 3 node shell formulation flag . . . .',I10/
1632 & 4X,'npts : shell properties default number of ',/
1633 & 4X,' integration points or layers. . . . . . . . .',I10/
1634 & 4X,'iframe : default solid frame formulation flag . . . . ',I10/
1635 & 4X,'ioffset: default shell offset flag. . . . . . . . . . ',I10/
1636 & 4X,'icontrol:default solid distortion control flag. . . . ',i10)
1637 5501 FORMAT(
1638 & 4x,'ISHFRAM: LOCAL SHELL FRAME DEFINITION. . . . . . . . .',i10)
1639 5502 FORMAT(
1640 & 4x,'IALEFVM: FVM FOR ALE MOMENTUM EQUATION. . . . . . . .',i10)
1641 5503 FORMAT(
1642 & 4x,'ISFINT : INTERNAL FORCES FORMULATION . . . . . . . . .',i10)
1643 5504 FORMAT(//
1644 & 4x,'MUSCL (MONOTONIC UPSTREAM-CENTERED SCHEME FOR CONSERVATION LAWS)'/
1645 & 4x,'----------------------------------------------------------------'/
1646 & 5x,'COMPRESSION COEFFICIENT (BETA). . . . . . . . . . . : ',1pg20.13/
1647 & 5x,'FORMULATION FLAG (IFLAG). . . . . . . . . . . . . . : ',i10)
1648 5700 FORMAT(//
1649 & 4x,'lagrange multiplier options'/
1650 & 4X,'---------------------------'/
1651 & 4X,'lagm_tol:convergence criterion. . . . . . . . . . . .',
1652 & 1PG20.13)
1653 RETURN
1654 5600 FORMAT(//4X,'sph global parameters'/
1655 & 4x,'---------------------'//
1656 & 4x,'ALPHA SORT :SORTING SECURITY COEFFICIENT . . . . . .',
1657 & 1pg20.13/,
1658 & 4x,'LNEIGH :MAXIMUM NUMBER OF COMPUTED NEIGHBOURS. .',i10/,
1659 & 4x,'NNEIGH :MAXIMUM NUMBER OF STORED NEIGHBOURS. . .',i10)
1660 5610 FORMAT(//4x,'SPH GLOBAL PARAMETERS'/
1661 & 4x,'---------------------'//
1662 & 4x,'ALPHA SORT :SORTING SECURITY COEFFICIENT . . . . . .',
1663 & 1pg20.13/,
1664 & 4x,'LNEIGH :MAXIMUM NUMBER OF COMPUTED NEIGHBOURS. .',i10/,
1665 & 4x,'NNEIGH :MAXIMUM NUMBER OF STORED NEIGHBOURS. . .',i10/,
1666 & 4x,'ITSOL2SPH :SOL2SPH PARTICLES ACTIVATION TYPE. . . .',i10/,
1667 & 4x,'(SOL2SPH PARTICLES ACTIVATION BASED ON PARTS)')
1668 5620 FORMAT(//4x,'SPH GLOBAL PARAMETERS'/
1669 & 4x,'---------------------'//
1670 & 4x,'ALPHA SORT :SORTING SECURITY COEFFICIENT . . . . . .',
1671 & 1pg20.13/,
1672 & 4x,'LNEIGH :MAXIMUM NUMBER OF COMPUTED NEIGHBOURS. .',i10/,
1673 & 4x,'NNEIGH :MAXIMUM NUMBER OF STORED NEIGHBOURS. . .',i10/,
1674 & 4x,'ITSOL2SPH :SOL2SPH PARTICLES ACTIVATION TYPE. . . .',i10/,
1675 & 4x,'(SOL2SPH PARTICLES ACTIVATION BASED ON SUBSETS)')
1676 5800 FORMAT(
1677 & //4x,'NODAL MASS DISTRIBUTION FLAG . . . . . . . . .',i10)
1678 5901 FORMAT(
1679 & 4x,'ISH3NFRAM: OLD LOCAL SH3N FRAME ACTIVATION . . . . . .',i10)
1680 999 CALL freerr(1)
1681 RETURN
1682 END
1683!||====================================================================
1684!|| find_dt1brick_engine ../starter/source/starter/contrl.F
1685!||--- called by ------------------------------------------------------
1686!|| contrl ../starter/source/starter/contrl.F
1687!||--- calls -----------------------------------------------------
1688!|| ancmsg ../starter/source/output/message/message.F
1689!||--- uses -----------------------------------------------------
1690!|| message_mod ../starter/share/message_module/message_mod.F
1691!||====================================================================
1692 SUBROUTINE find_dt1brick_engine()
1693C-----------------------------------------------
1694C M o d u l e s
1695C-----------------------------------------------
1696 USE message_mod
1697 USE inoutfile_mod
1698C-----------------------------------------------
1699C I m p l i c i t T y p e s
1700C-----------------------------------------------
1701#include "implicit_f.inc"
1702C-----------------------------------------------
1703C C o m m o n B l o c k s
1704C-----------------------------------------------
1705#include "com04_c.inc"
1706#include "scr15_c.inc"
1707#include "scr17_c.inc"
1708C-----------------------------------------------
1709C L o c a l V a r i a b l e s
1710C-----------------------------------------------
1711 INTEGER IO_ERR1, NITER
1712 CHARACTER FILNAM*109, KEYA*80, KEYA2*80
1713
1714 INTEGER :: LEN_TMP_NAME
1715 CHARACTER(len=4096) :: TMP_NAME
1716C-----------------------------------------------
1717 io_err1=0
1718 idttsh=0
1719C-----------------------------------------------
1720C Lecture des donnees dans le fichier engine
1721C { Reading data from the engine file }
1722C-----------------------------------------------
1723 filnam=rootnam(1:rootlen)//'_0001.rad'
1724
1725 tmp_name=infile_name(1:infile_name_len)//filnam(1:rootlen+9)
1726 len_tmp_name = infile_name_len+rootlen+9
1727 OPEN(unit=71,file=tmp_name(1:len_tmp_name),
1728 . access='SEQUENTIAL',status='OLD',iostat=io_err1)
1729C
1730 IF (io_err1/=0) THEN
1731 filnam=rootnam(1:rootlen)//'D01'
1732 tmp_name=infile_name(1:infile_name_len)//filnam(1:rootlen+3)
1733 len_tmp_name = infile_name_len+rootlen+3
1734 OPEN(unit=71,file=tmp_name(1:len_tmp_name),
1735 . access='SEQUENTIAL',status='OLD',iostat=io_err1)
1736 ENDIF
1737
1738 IF (io_err1==0) THEN
1739C
1740 10 READ(71,'(A)',END=20) keya
1741C
1742 IF(keya(1:10)=='/DT1/BRICK') THEN
1743 idt1sol=1
1744 ELSEIF(keya(1:9)=='/DT1TET10') THEN
1745 niter=0
1746 READ(keya(11:20),'(I10)',err=30) niter
1747 30 CONTINUE
1748 IF(niter == 0)niter=1
1749 idt1tet10=niter+1
1750 ELSEIF(keya(1:6)=='/DTTSH') THEN
1751 idttsh=1
1752 ENDIF
1753C
1754 GOTO 10
1755C
1756 20 CONTINUE
1757
1758 CLOSE(71)
1759C
1760 ELSE
1761C
1762 IF(numtetra4/=0) CALL ancmsg(msgid=1589,
1763 . msgtype=msgwarning,
1764 . anmode=aninfo_blind_2,
1765 . c1=rootnam(1:rootlen)//'_0001.rad',
1766 . c2=rootnam(1:rootlen)//'d01')
1767C
1768 IF(NUMELS10/=0) CALL ANCMSG(MSGID=1606,
1769 . MSGTYPE=MSGWARNING,
1770 . ANMODE=ANINFO_BLIND_2,
1771 . C1=ROOTNAM(1:ROOTLEN)//'_0001.rad',
1772 . C2=ROOTNAM(1:ROOTLEN)//'d01')
1773C
1774 ENDIF
1775C-------------------------------------------
1776 RETURN
1777 END
1778!||====================================================================
1779!|| ini_h3dtmax_engine ../starter/source/starter/contrl.F
1780!||--- called by ------------------------------------------------------
1781!|| lectur ../starter/source/starter/lectur.F
1782!||--- calls -----------------------------------------------------
1783!|| read_h3dtmax_key ../starter/source/starter/contrl.F
1784!||--- uses -----------------------------------------------------
1785!|| message_mod ../starter/share/message_module/message_mod.F
1786!||====================================================================
1787 SUBROUTINE INI_H3DTMAX_ENGINE(IPARG,IPART,IPARTS,IPARTC,IPARTG,IDDLEVEL)
1788C-----------------------------------------------
1789C M o d u l e s
1790C-----------------------------------------------
1791 USE MESSAGE_MOD
1792 USE INOUTFILE_MOD
1793 USE OUTMAX_MOD
1794 USE NAMES_AND_TITLES_MOD , ONLY : NCHARLINE
1795C-----------------------------------------------
1796C I m p l i c i t T y p e s
1797C-----------------------------------------------
1798#include "implicit_f.inc"
1799C-----------------------------------------------
1800C C o m m o n B l o c k s
1801C-----------------------------------------------
1802#include "com01_c.inc"
1803#include "com04_c.inc"
1804#include "scr15_c.inc"
1805#include "scr17_c.inc"
1806#include "param_c.inc"
1807C-----------------------------------------------
1808C D u m m y A r g u m e n t s
1809C-----------------------------------------------
1810 INTEGER ,INTENT(IN) :: IDDLEVEL
1811 INTEGER, DIMENSION(NPARG,NGROUP) ,INTENT(IN):: IPARG
1812 INTEGER ,DIMENSION(LIPART1,NPART),INTENT(IN):: IPART
1813 INTEGER ,DIMENSION(NUMELS),INTENT(IN):: IPARTS
1814 INTEGER ,DIMENSION(NUMELC),INTENT(IN):: IPARTC
1815 INTEGER ,DIMENSION(NUMELTG),INTENT(IN):: IPARTG
1816C-----------------------------------------------
1817C L o c a l V a r i a b l e s
1818C-----------------------------------------------
1819 INTEGER I, J, K, N ,NELC , NELTG , IP , NPRT , IPRT
1820 INTEGER IH3D,NG,ITY,NFT,IKEY,K_LEN
1821
1822 INTEGER :: TMAX_IPART(NPART),NKPART(NPART+1,NKEYMAX+1)
1823 CHARACTER(LEN=NCHARLINE) :: CARTE
1824C-----------------------------------------------
1825 IF (ALLOCATED(IKEYMAX)) DEALLOCATE(IKEYMAX)
1826 IF (ALLOCATED(IPART_OK)) DEALLOCATE(IPART_OK)
1827 ALLOCATE(IKEYMAX(NKEYMAX),IPART_OK(NGROUP,NKEYMAX))
1828 IKEYMAX =0
1829 IPART_OK = 0
1830 LMAX_DIS = 0
1831 LMAX_VEL = 0
1832 LMAX_NSIG = 0
1833 LMAX_NSTRA = 0
1834 NKPART(1:NPART+1,1:NKEYMAX+1) = 0
1835C-----------------------------------------------
1836C Lecture des donnees dans le fichier engine
1837C { Reading data from the engine file }
1838C-----------------------------------------------
1839 IH3D=1
1840 K_LEN = 19
1841 CARTE(1:K_LEN) = '/h3d/elem/vonm/tmax'
1842 CALL READ_H3DTMAX_KEY(CARTE,K_LEN,IKEY,NKPART(1,IH3D),IPART)
1843 IF (IKEY >0) THEN
1844 IKEYMAX(IH3D) = 1
1845 ELSE
1846 K_LEN = 20
1847 CARTE(1:K_LEN) = '/h3d/solid/vonm/tmax'
1848 CALL READ_H3DTMAX_KEY(CARTE,K_LEN,IKEY,NKPART(1,IH3D),IPART)
1849 IF (IKEY >0) IKEYMAX(IH3D) = 1
1850 CARTE(1:K_LEN) = '/h3d/shell/vonm/tmax'
1851 CALL READ_H3DTMAX_KEY(CARTE,K_LEN,IKEY,NKPART(1,IH3D),IPART)
1852 IF (IKEY >0) IKEYMAX(IH3D) = 1
1853 END IF
1854 IH3D=2
1855 K_LEN = 20
1856 CARTE(1:K_LEN) = '/h3d/elem/sigeq/tmax'
1857 CALL READ_H3DTMAX_KEY(CARTE,K_LEN,IKEY,NKPART(1,IH3D),IPART)
1858 IF (IKEY >0) THEN
1859 IKEYMAX(IH3D) = 1
1860 ELSE
1861 K_LEN = 21
1862 CARTE(1:K_LEN) = '/h3d/solid/sigeq/tmax'
1863 CALL READ_H3DTMAX_KEY(CARTE,K_LEN,IKEY,NKPART(1,IH3D),IPART)
1864 IF (IKEY >0) IKEYMAX(IH3D) = 1
1865 CARTE(1:K_LEN) = '/h3d/shell/sigeq/tmax'
1866 CALL READ_H3DTMAX_KEY(CARTE,K_LEN,IKEY,NKPART(1,IH3D),IPART)
1867 IF (IKEY >0) IKEYMAX(IH3D) = 1
1868 END IF
1869 IH3D=3
1870 K_LEN = 19
1871 CARTE(1:K_LEN) = '/h3d/elem/ener/tmax'
1872 CALL READ_H3DTMAX_KEY(CARTE,K_LEN,IKEY,NKPART(1,IH3D),IPART)
1873 IF (IKEY >0) THEN
1874 IKEYMAX(IH3D) = 1
1875 ELSE
1876 K_LEN = 20
1877 CARTE(1:K_LEN) = '/h3d/solid/ener/tmax'
1878 CALL READ_H3DTMAX_KEY(CARTE,K_LEN,IKEY,NKPART(1,IH3D),IPART)
1879 IF (IKEY >0) IKEYMAX(IH3D) = 1
1880 CARTE(1:K_LEN) = '/h3d/shell/ener/tmax'
1881 CALL READ_H3DTMAX_KEY(CARTE,K_LEN,IKEY,NKPART(1,IH3D),IPART)
1882 IF (IKEY >0) IKEYMAX(IH3D) = 1
1883 END IF
1884 IH3D=4
1885 K_LEN = 19
1886 CARTE(1:K_LEN) = '/h3d/elem/dama/tmax'
1887 CALL READ_H3DTMAX_KEY(CARTE,K_LEN,IKEY,NKPART(1,IH3D),IPART)
1888 IF (IKEY >0) THEN
1889 IKEYMAX(IH3D) = 1
1890 ELSE
1891 K_LEN = 20
1892 CARTE(1:K_LEN) = '/h3d/solid/dama/tmax'
1893 CALL READ_H3DTMAX_KEY(CARTE,K_LEN,IKEY,NKPART(1,IH3D),IPART)
1894 IF (IKEY >0) IKEYMAX(IH3D) = 1
1895 CARTE(1:K_LEN) = '/h3d/shell/dama/tmax'
1896 CALL READ_H3DTMAX_KEY(CARTE,K_LEN,IKEY,NKPART(1,IH3D),IPART)
1897 IF (IKEY >0) IKEYMAX(IH3D) = 1
1898 END IF
1899 IH3D=5
1900 K_LEN = 26
1901 CARTE(1:K_LEN) = '/h3d/elem/tens/stress/tmax'
1902 CALL READ_H3DTMAX_KEY(CARTE,K_LEN,IKEY,NKPART(1,IH3D),IPART)
1903 IF (IKEY >0) THEN
1904 IKEYMAX(IH3D) = 1
1905 ELSE
1906 K_LEN = 27
1907 CARTE(1:K_LEN) = '/h3d/solid/tens/stress/tmax'
1908 CALL READ_H3DTMAX_KEY(CARTE,K_LEN,IKEY,NKPART(1,IH3D),IPART)
1909 IF (IKEY >0) IKEYMAX(IH3D) = 1
1910 CARTE(1:K_LEN) = '/h3d/shell/tens/stress/tmax'
1911 CALL READ_H3DTMAX_KEY(CARTE,K_LEN,IKEY,NKPART(1,IH3D),IPART)
1912 IF (IKEY >0) IKEYMAX(IH3D) = 1
1913 END IF
1914 IH3D=6
1915 K_LEN = 26
1916 CARTE(1:K_LEN) = '/h3d/elem/tens/strain/tmax'
1917 CALL READ_H3DTMAX_KEY(CARTE,K_LEN,IKEY,NKPART(1,IH3D),IPART)
1918 IF (IKEY >0) THEN
1919 IKEYMAX(IH3D) = 1
1920 ELSE
1921 K_LEN = 27
1922 CARTE(1:K_LEN) = '/h3d/solid/tens/strain/tmax'
1923 CALL READ_H3DTMAX_KEY(CARTE,K_LEN,IKEY,NKPART(1,IH3D),IPART)
1924 IF (IKEY >0) IKEYMAX(IH3D) = 1
1925 CARTE(1:K_LEN) = '/h3d/shell/tens/strain/tmax'
1926 CALL READ_H3DTMAX_KEY(CARTE,K_LEN,IKEY,NKPART(1,IH3D),IPART)
1927 IF (IKEY >0) IKEYMAX(IH3D) = 1
1928 END IF
1929 IH3D=7 ! NKEYMAX+1
1930 K_LEN = 9
1931 CARTE(1:K_LEN) = '/h3d/part'
1932 CALL READ_H3DTMAX_KEY(CARTE,K_LEN,IKEY,NKPART(1,IH3D),IPART)
1933 IF (IKEY == 0 ) NKPART(1:NPART,IH3D)=1
1934 IH3D=8
1935 K_LEN =18
1936 CARTE(1:K_LEN) = '/h3d/noda/dis/tmax'
1937 CALL READ_H3DTMAX_KEY(CARTE,K_LEN,IKEY,TMAX_IPART,IPART)
1938 IF (IKEY > 0 ) LMAX_DIS = 3
1939 IH3D=9
1940 K_LEN =18
1941 CARTE(1:K_LEN) = '/h3d/noda/vel/tmax'
1942 CALL read_h3dtmax_key(carte,k_len,ikey,tmax_ipart,ipart)
1943 IF (ikey > 0 ) lmax_vel = 3
1944 ih3d=10
1945 k_len =18
1946 carte(1:k_len) = '/H3D/NODA/GPS/TMAX'
1947 CALL read_h3dtmax_key(carte,k_len,ikey,tmax_ipart,ipart)
1948 IF (ikey > 0 ) lmax_nsig = 6
1949 ih3d=11
1950 k_len =23
1951 carte(1:k_len) = '/H3D/NODA/GPSTRAIN/TMAX'
1952 CALL read_h3dtmax_key(carte,k_len,ikey,tmax_ipart,ipart)
1953 IF (ikey > 0 ) lmax_nstra = 6
1954C
1955 ih3d = 0
1956 DO i=1,nkeymax
1957 ih3d = ih3d + ikeymax(i)
1958 END DO
1959 IF (ih3d>0) THEN
1960 DO i=1,nkeymax
1961 nkpart(1:npart,i)=nkpart(1:npart,nkeymax+1)*nkpart(1:npart,i)
1962 END DO
1963C
1964 DO ng=1,ngroup
1965 nft=iparg(3,ng)+1
1966 ity=iparg(5,ng)
1967C
1968 iprt = 0
1969 SELECT CASE (ity)
1970 CASE(1)
1971 iprt = iparts(nft)
1972 CASE(3)
1973 iprt = ipartc(nft)
1974 CASE(7)
1975 iprt = ipartg(nft)
1976 END SELECT
1977 IF(iprt>0) ipart_ok(ng,1:nkeymax) = nkpart(iprt,1:nkeymax)
1978 END DO
1979 END IF
1980C-------------------------------------------
1981 RETURN
1982 END
1983!||====================================================================
1984!|| read_h3dtmax_key ../starter/source/starter/contrl.f
1985!||--- called by ------------------------------------------------------
1986!|| ini_h3dtmax_engine ../starter/source/starter/contrl.F
1987!||--- uses -----------------------------------------------------
1988!|| message_mod ../starter/share/message_module/message_mod.F
1989!||====================================================================
1990 SUBROUTINE read_h3dtmax_key(KEY_TM,KEY_LEN,IFUND,NTM_PART,IPART)
1991C-----------------------------------------------
1992C M o d u l e s
1993C-----------------------------------------------
1994 USE message_mod
1995 USE inoutfile_mod
1996 USE names_and_titles_mod , ONLY : ncharline
1997C-----------------------------------------------
1998C I m p l i c i t T y p e s
1999C-----------------------------------------------
2000#include "implicit_f.inc"
2001C-----------------------------------------------
2002C C o m m o n B l o c k s
2003C-----------------------------------------------
2004#include "com01_c.inc"
2005#include "com04_c.inc"
2006#include "scr15_c.inc"
2007#include "scr17_c.inc"
2008#include "param_c.inc"
2009C-----------------------------------------------
2010C D u m m y A r g u m e n t s
2011C-----------------------------------------------
2012 INTEGER, INTENT(IN):: KEY_LEN
2013 CHARACTER(len=KEY_LEN) :: KEY_TM
2014 INTEGER ,DIMENSION(LIPART1,NPART),INTENT(IN):: IPART
2015 INTEGER, INTENT(OUT):: IFUND
2016 INTEGER ,DIMENSION(NPART),INTENT(INOUT):: NTM_PART
2017C-----------------------------------------------
2018C L o c a l V a r i a b l e s
2019C-----------------------------------------------
2020 INTEGER I, J, K, N ,NELC , NELTG , IP , NPRT , IPRT
2021 INTEGER IH3D,NG,ITY,NFT,IO_ERR1,LEN_C,IC
2022 CHARACTER FILNAM*109, KEYA*80, KEYA2*80
2023 CHARACTER(LEN=NCHARLINE) :: CARTE
2024
2025 INTEGER :: LEN_TMP_NAME,TMAX_IPART(NPART)
2026 CHARACTER(len=4096) :: TMP_NAME
2027C-----------------------------------------------
2028 ifund=0
2029 io_err1 = 0
2030 nprt = 0
2031C-----------------------------------------------
2032C Lecture des donnees dans le fichier engine
2033C-----------------------------------------------
2034 filnam=rootnam(1:rootlen)//'_0001.rad'
2035
2036 tmp_name=infile_name(1:infile_name_len)//filnam(1:rootlen+9)
2037 len_tmp_name = infile_name_len+rootlen+9
2038 OPEN(unit=71,file=tmp_name(1:len_tmp_name),
2039 . access='SEQUENTIAL',status='OLD',iostat=io_err1)
2040
2041 DO WHILE (io_err1 == 0 .AND. ifund==0)
2042C
2043 READ(71,fmt='(A)',END=20) keya
2044 IF(keya(1:key_len) == key_tm ) THEN
2045C---- read ipart_id if there are
2046 ifund = 1
2047 ic = 1
2048 nprt = npart
2049 DO WHILE(ic == 1)
2050 READ(71,fmt='(A)',err=20,END=20) carte
2051 IF(carte(1:1) == '#' .OR. carte(1:1) == '$') cycle
2052 ic = 0
2053 len_c = len_trim(carte)
2054 IF(carte(len_c:len_c)==char(13)) len_c = len_c - 1
2055 IF (carte(1:1) == '/' .OR. len_c==0) THEN
2056 ELSE
2057 j=1
2058 nprt = 0
2059 DO WHILE(carte(1:1) /= '/' .AND. carte(1:1) /= '#' .AND.
2060 . carte(1:1) /= '$' .AND. len_c /= 0)
2061 DO WHILE (j <= len_c)
2062 IF(carte(j:j) /= ' ') THEN
2063 k=j
2064 DO WHILE(carte(k:k) /= ' ' .AND. k < len_c)
2065 k=k+1
2066 ENDDO
2067 nprt = nprt + 1
2068 READ(carte(j:k),'(I10)') iprt
2069 tmax_ipart(nprt) = iprt
2070 j = k
2071 ENDIF
2072 j = j +1
2073 END DO
2074 READ(71,fmt='(a)',ERR=20,END=20) CARTE
2075 LEN_C = LEN_TRIM(CARTE)
2076 IF(CARTE(LEN_C:LEN_C)==CHAR(13)) LEN_C = LEN_C - 1
2077 END DO
2078 IF (NPRT==0) THEN
2079 NPRT = NPART
2080 NTM_PART(1:NPRT) = 1
2081 END IF
2082 END IF
2083 END DO !(IC == 1) remove comments
2084 END IF
2085C
2086 END DO
2087C
2088 20 CONTINUE
2089
2090 CLOSE(71)
2091C
2092.AND..OR. IF (IFUND>0 (NPRT == NPART NPRT ==0)) THEN
2093 NTM_PART(1:NPRT) = 1
2094 ELSEIF (IFUND>0) THEN
2095 DO I=1,NPRT
2096 IP=0
2097 IPRT = TMAX_IPART(I)
2098 DO J=1,NPART
2099 IF(IPART(4,J)==IPRT)IP=J
2100 END DO
2101 IF (IP > 0) NTM_PART(IP)=1
2102 IP = 0
2103 DO J=1,NPART
2104 IP = IP +NTM_PART(J)
2105 END DO
2106 IF (IP==0) NTM_PART(1:NPART) = 1
2107 ENDDO
2108 END IF
2109C-------------------------------------------
2110 RETURN
2111 END
read_h3dtmax_key
subroutine read_h3dtmax_key(key_tm, key_len, ifund, ntm_part, ipart)
Definition contrl.F:1991
find_dt1brick_engine
subroutine find_dt1brick_engine()
Definition contrl.F:1693
contrl
subroutine contrl(multi_fvm, lsubmodel, is_dyna, detonators, user_windows, mat_elem, names_and_titles, lipart1, defaults, glob_therm, pblast, output)
Definition contrl.F:83
convec
subroutine convec(ibcv, fconv, npc, tf, x, temp, nsensor, sensor_tab, fthe, iad, fthesky, python, glob_therm)
Definition convec.F:38
my_real
#define my_real
Definition cppsort.cpp:32
damping
subroutine damping(nodft, nodlt, v, vr, a, ar, damp, ms, in, igrnod, dim, itask, weight, tagslv_rby, wfext)
Definition damping.F:882
detcord
subroutine detcord(detonator_cord, x, xc, yc, zc, vdet, vdet2, alt, bt, tb, has_detonator, iopt)
Definition detcord.F:34
eig
subroutine eig(k_diag, k_lt, iadk, jdik, ms, in, nddl, ndof, nnzl, x, d, v, a, bufel, ixs, ixq, ixc, ixt, ixp, ixr, ixtg, pm, geo, cont, icut, skew, xcut, fint, itab, fext, fopt, anin, lpby, npby, nstrf, rwbuf, nprw, tani, elbuf_tab, matparam_tab, dd_iad, fr_iad, dd_front, cluster, weight, eani, ipart, rby, nom_opt, igrsurf, bufsf, idata, rdata, bufmat, bufgeo, kxx, ixx, kxsp, ixsp, nod2sp, spbuf, ixs10, ixs20, ixs16, vr, monvol, volmon, ipm, igeo, iparg, eigipm, eigibuf, eigrpm, ldiag, ljdik, ljdik2, ikc, maxncv, thke, nms, nint2, iint2, ipari, intbuf_tab, nodglob, iad_elem, fr_elem, fr_sec, fr_rby2, iad_rby2, fr_wall, inloc, iddl, partsav, fncont, ftcont, temp, err_thk_sh4, err_thk_sh3, irbe2, irbe3, lrbe2, lrbe3, fr_rbe2, fr_rbe3m, iad_rbe2, weight_md, fcluster, mcluster, xfem_tab, w, nv46, nercvois, nesdvois, lercvois, lesdvois, crkedge, indx_crk, xedge4n, xedge3n, stack, sph2sol, stifn, stifr, drape_q4, drape_t3, h3d_data, subset, igrnod, fcont_max, fncontp2, ftcontp2, ale_connectivity, glob_therm)
Definition eig.F:73
alpha
#define alpha
Definition eval.h:35
hm_elem_count
subroutine hm_elem_count(elem_type, hm_elem_number, is_dyna)
Definition hm_elem_count.F:33
hm_get_intv
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
Definition hm_get_intv.F:523
hm_get_string
subroutine hm_get_string(name, sval, size, is_available)
Definition hm_get_string.F:76
hm_option_count
subroutine hm_option_count(entity_type, hm_option_number)
Definition hm_option_count.F:105
hm_option_next
subroutine hm_option_next()
Definition hm_option_next.F:39
hm_option_start
subroutine hm_option_start(entity_type)
Definition hm_option_start.F:251
hm_preread_properties
subroutine hm_preread_properties(igeo, nsphsol, nply, nsub, nisub, lsubmodel, defaults)
Definition hm_preread_properties.F:40
hm_read_analy
subroutine hm_read_analy(nanaly, iparith, ipari0, lsubmodel)
Definition hm_read_analy.F:37
hm_read_implicit
subroutine hm_read_implicit(lsubmodel)
Definition hm_read_implicit.F:36
hm_read_ioflag
subroutine hm_read_ioflag(lsubmodel)
Definition hm_read_ioflag.F:38
hm_read_spmd
subroutine hm_read_spmd(lsubmodel)
Definition hm_read_spmd.F:38
hm_read_unit
subroutine hm_read_unit(unitab, lsubmodel)
Definition hm_read_unit.F:45
hm_read_upwind
subroutine hm_read_upwind(jupw, eta1, eta2, eta3, lsubmodel, unitab)
Definition hm_read_upwind.F:39
inivel
subroutine inivel(v, vr, svr, itabm1)
Definition inivel.F:35
hm_option_read_mod::hm_option_read_key
Definition hm_option_read_mod.F:410
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
merge
subroutine merge(x, itab, itabm1, cmerge, imerge, imerge2, iadmerge2, nmerge_tot)
Definition merge.F:36
id
initmumps id
Definition multiplerhs_example.m:4
for
for(i8=*sizetab-1;i8 >=0;i8--)
Definition mumps_common.c:91
ale_ebcs_mod
Definition ale_ebcs_mod.F:40
ale_mod
Definition ale_mod.F:149
ale_mod::ale
type(ale_) ale
Definition ale_mod.F:249
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:54
alemuscl_mod::alemuscl_param
type(alemuscl_param_) alemuscl_param
Definition alemuscl_mod.F:75
anim_mod
Definition anim_mod.F:96
check_mod
Definition check_mod.F:34
detonators_mod
Definition detonators_mod.F:79
grp_size_mod
Definition grp_size_mod.F:29
grp_size_mod::grp_size_bool
logical grp_size_bool
Definition grp_size_mod.F:39
grp_size_mod::grp_size
integer grp_size
Definition grp_size_mod.F:40
hm_option_read_mod
Definition hm_option_read_mod.F:408
inivol_def_mod
Definition inivol_mod.F:45
inivol_def_mod::inivol
type(inivol_struct_), dimension(:), allocatable inivol
Definition inivol_mod.F:84
inivol_def_mod::num_inivol
integer num_inivol
Definition inivol_mod.F:85
inoutfile_mod
Definition inoutfile_mod.F:96
inoutfile_mod::infile_name_len
integer infile_name_len
Definition inoutfile_mod.F:114
inoutfile_mod::infile_name
character(len=infile_char_len) infile_name
Definition inoutfile_mod.F:116
message_mod
Definition message_mod.F:1249
names_and_titles_mod
Definition names_and_titles_mod.F:997
names_and_titles_mod::ltitle
integer, parameter ltitle
Definition names_and_titles_mod.F:1007
names_and_titles_mod::ncharfield
integer, parameter ncharfield
Definition names_and_titles_mod.F:1005
names_and_titles_mod::ncharline
integer, parameter ncharline
Definition names_and_titles_mod.F:1001
pinchtype_mod
Definition pinchtype_mod.F:69
refsta_mod
Definition refsta_mod.F:32
restmod::iframe
integer, dimension(:), pointer iframe
Definition restart_mod.F:108
restmod::unitab
type(unit_type_) unitab
Definition restart_mod.F:211
restmod::ipart
integer, dimension(:), allocatable, target ipart
Definition restart_mod.F:60
restmod::monvol
integer, dimension(:), allocatable monvol
Definition restart_mod.F:60
restmod::igeo
integer, dimension(:), allocatable igeo
Definition restart_mod.F:83
restmod::iparg
integer, dimension(:), allocatable iparg
Definition restart_mod.F:60
setdef_mod
Definition setdef_mod.F:107
submodel_mod
Definition submodel_mod.F:712
user_windows_mod
Definition user_windows_mod.F:48
nbadigemesh
subroutine nbadigemesh(lsubmodel, numnusr)
Definition nbadigemesh.F:45
nbadmesh
subroutine nbadmesh(lsubmodel, numnusr, unitab)
Definition nbadmesh.F:43
radiation
subroutine radiation(ibcr, fradia, npc, tf, x, temp, nsensor, sensor_tab, fthe, iad, fthesky, python, glob_therm)
Definition radiation.F:38
radioss
program radioss
Definition radioss.F:34
ancmsg
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)
Definition message.F:889
freerr
subroutine freerr(it)
Definition freform.F:506
istr
integer function istr(str)
Definition istr.F:31
starter
program starter
Definition starter.F:39
static
subroutine static(v, vr, a, ar, ms, in, igrnod, weight_md, wfext)
Definition static.F:33
detonators_mod::detonators_struct_
Definition detonators_mod.F:154
names_and_titles_mod::names_and_titles_
Definition names_and_titles_mod.F:1016
submodel_mod::submodel_data
Definition submodel_mod.F:723
user_windows_mod::user_windows_
Definition user_windows_mod.F:58
tmax_ipart
subroutine tmax_ipart(iparg, ipart, iparts, ipartc, ipartg, h3d_data)
Definition tmax_ipart.F:34
upwind
subroutine upwind(rho, vis, vdx, vdy, vdz, r, s, t, deltax, gam, nel)
Definition upwind.F:35