OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
lectur.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!|| lectur ../engine/source/input/lectur.F
25!||--- called by ------------------------------------------------------
26!|| radioss2 ../engine/source/engine/radioss2.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../engine/source/output/message/message.F
29!|| anim_build_index_all ../engine/source/output/anim/reader/anim_build_index_all.F
30!|| arret ../engine/source/system/arret.F
31!|| chkipari ../engine/source/interfaces/interf/chkstfn3.F
32!|| finter ../engine/source/tools/curve/finter.F
33!|| fr_rlale ../engine/source/mpi/kinematic_conditions/fr_rlink1.F
34!|| fr_rlink1 ../engine/source/mpi/kinematic_conditions/fr_rlink1.F
35!|| lcbcsf ../engine/source/constraints/general/bcs/lcbcsf.f
36!|| leccut ../engine/source/tools/sect/leccut.F
37!|| lecdamp ../engine/source/input/lecdamp.F
38!|| leceig ../engine/source/input/leceig.F
39!|| lecflsw ../engine/source/fluid/lecflsw.F
40!|| lecfun ../engine/source/tools/curve/lecfun.F
41!|| lecfvbag ../engine/source/input/lecfvbag.F
42!|| lecfvbag1 ../engine/source/input/lecfvbag1.F
43!|| lecfxinp ../engine/source/input/lecfxinp.F
44!|| lech3d ../engine/source/output/h3d/h3d_build_fortran/lech3d.F
45!|| lecimpl ../engine/source/input/lectur.F
46!|| lecinv ../engine/source/input/lecinv.F
47!|| lecnoise ../engine/source/general_controls/computation/lecnoise.F
48!|| lecstat ../engine/source/input/lecstat.F
49!|| ngr2usr ../engine/source/input/freform.F
50!|| prelech3d ../engine/source/output/h3d/h3d_build_fortran/prelech3d.F90
51!|| rbyonf ../engine/source/constraints/general/rbody/rbyonf.F
52!|| read_dynain ../engine/source/output/dynain/read_dynain.F
53!|| rlink0 ../engine/source/constraints/general/rlink/rlink0.F
54!|| spmd_chkw ../engine/source/mpi/generic/spmd_chkw.F
55!|| spmd_glob_isum9 ../engine/source/mpi/interfaces/spmd_th.F
56!|| spmd_ibcast ../engine/source/mpi/generic/spmd_ibcast.F
57!|| spmd_wiout ../engine/source/mpi/generic/spmd_wiout.F
58!||--- uses -----------------------------------------------------
59!|| ale_mod ../common_source/modules/ale/ale_mod.F
60!|| alemuscl_mod ../common_source/modules/ale/alemuscl_mod.F
61!|| check_mod ../common_source/modules/check_mod.F
62!|| dt_mod ../engine/source/modules/dt_mod.F
63!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
64!|| element_mod ../common_source/modules/elements/element_mod.F90
65!|| fvbag_mod ../engine/share/modules/fvbag_mod.f
66!|| fxb_mod ../engine/share/modules/fxb_mod.F
67!|| glob_therm_mod ../common_source/modules/mat_elem/glob_therm_mod.F90
68!|| groupdef_mod ../common_source/modules/groupdef_mod.F
69!|| h3d_mod ../engine/share/modules/h3d_mod.F
70!|| intbufdef_mod ../common_source/modules/interfaces/intbufdef_mod.F90
71!|| intstamp_glob_mod ../engine/share/modules/intstamp_glob_mod.F
72!|| loads_mod ../common_source/modules/loads/loads_mod.F90
73!|| mat_elem_mod ../common_source/modules/mat_elem/mat_elem_mod.F90
74!|| message_mod ../engine/share/message_module/message_mod.F
75!|| multi_fvm_mod ../common_source/modules/ale/multi_fvm_mod.F90
76!|| my_alloc_mod ../common_source/tools/memory/my_alloc.F90
77!|| names_and_titles_mod ../common_source/modules/names_and_titles_mod.F
78!|| output_mod ../common_source/modules/output/output_mod.F90
79!|| pblast_mod ../common_source/modules/loads/pblast_mod.F90
80!|| prelech3d_mod ../engine/source/output/h3d/h3d_build_fortran/prelech3d.F90
81!|| sensor_mod ../common_source/modules/sensor_mod.F90
82!|| stack_mod ../engine/share/modules/stack_mod.F
83!|| state_mod ../common_source/modules/state_mod.F
84!|| table_mod ../engine/share/modules/table_mod.F
85!||====================================================================
86 SUBROUTINE lectur(ICODE ,ISKEW ,ISKWN ,IXTG ,IXS ,
87 2 IXQ ,IXC ,IXT ,IXP ,IXR ,
88 3 ITAB ,ITABM1,NPC ,IPARG ,IGRV ,
89 4 IBGR ,IPARI ,
90 5 NPBY ,LPBY ,NNLINK ,LLINK ,LINALE,
91 6 NEFLSW ,NNFLSW,ICUT ,INOISE ,
92 7 X ,V ,VR ,MS ,IN ,
93 8 SKEW ,PLD ,RBY ,
94 9 WA ,CRFLSW,XCUT ,DAMPR ,
95 A IGRNOD ,KXSP ,WEIGHT ,FR_RBY2,FR_RL ,
96 B PARTSAV,IPART ,PM ,
97 C MONVOL ,VOLMON,IPART_STATE,GEO,TABLE ,
98 D IFRAME ,XFRAME,ELBUF_STR,IGEO ,INTBUF_TAB,
99 E IPM ,H3D_DATA, MULTI_FVM,
100 F IGRPART ,TAG_SKINS6,ICFIELD,LCFIELD,TAGSLV_RBY,
101 G MDS_LABEL,MDS_OUTPUT_TABLE,MDS_NMAT,MAX_DEPVAR,
102 H MDS_NDEPSVAR,STACK,IBCL,ILOADP,LLOADP,SENSORS,
103 I DYNAIN_DATA,DT ,LOADS ,OUTPUT,NAMES_AND_TITLES,
104 J MAT_PARAM,GLOB_THERM,PBLAST)
105C-----------------------------------------------
106C M o d u l e s
107C-----------------------------------------------
108 USE table_mod
109 USE message_mod
110 USE mat_elem_mod
111 USE intbufdef_mod
112 USE stack_mod
113 USE h3d_mod
114 USE multi_fvm_mod
115 USE groupdef_mod
116 USE intstamp_glob_mod
117 USE fvbag_mod
118 USE check_mod
119 USE sensor_mod
120 USE output_mod
121 USE state_mod
122 USE alemuscl_mod , only:alemuscl_param
123 USE ale_mod
124 USE dt_mod
125 USE loads_mod
126 USE fxb_mod
127 USE my_alloc_mod
128 USE names_and_titles_mod
129 USE elbufdef_mod
130 USE glob_therm_mod
131 USE pblast_mod
132 USE prelech3d_mod
133 use element_mod , only : nixs,nixq,nixc,nixt,nixp,nixr,nixtg
134C-----------------------------------------------
135C I m p l i c i t T y p e s
136C-----------------------------------------------
137#include "implicit_f.inc"
138C-----------------------------------------------
139C C o m m o n B l o c k s
140C-----------------------------------------------
141#include "com01_c.inc"
142#include "com04_c.inc"
143#include "com06_c.inc"
144#include "com08_c.inc"
145#include "com09_c.inc"
146#include "param_c.inc"
147#include "warn_c.inc"
148#include "scr02_c.inc"
149#include "scr03_c.inc"
150#include "scr05_c.inc"
151#include "scr06_c.inc"
152#include "scr07_c.inc"
153#include "scr14_c.inc"
154#include "scr16_c.inc"
155#include "scr17_c.inc"
156#include "scr18_c.inc"
157#include "cong1_c.inc"
158#include "cong2_c.inc"
159#include "scrfs_c.inc"
160#include "stati_c.inc"
161#include "statr_c.inc"
162#include "units_c.inc"
163#include "scrcut_c.inc"
164#include "scrnoi_c.inc"
165#include "parit_c.inc"
166#include "chara_c.inc"
167#include "task_c.inc"
168#include "sphcom.inc"
169#include "impl1_c.inc"
170#include "tabsiz_c.inc"
171#include "remesh_c.inc"
172#include "sms_c.inc"
173#include "rad2r_c.inc"
174#include "inter22.inc"
175#include "userlib.inc"
176#include "spmd_c.inc"
177#include "intstamp_c.inc"
178#include "couple_c.inc"
179C-----------------------------------------------
180C D u m m y A r g u m e n t s
181C-----------------------------------------------
182 INTEGER IPARG(NPARG,NGROUP), IPARI(NPARI,*), IXS(NIXS,*),
183 . IXQ(NIXQ,*), IXC(NIXC,*), IXT(NIXT,*), IXP(NIXP,*),
184 . IXR(NIXR,*), IXTG(NIXTG,*),ITAB(*), ITABM1(*),
185 . ISKWN(LISKN,*), NPBY(NNPBY,*),NNLINK(*) ,LLINK(*) ,LINALE(*),
186 . ICODE(*) ,ISKEW(*),NPC(*),NEFLSW(*),NNFLSW(*),ICUT(*),
187 . INOISE(*),IGRV(NIGRV,*),IBGR(*),
188 . LPBY(*),KXSP(NISP,*),WEIGHT(*),FR_RBY2(*),
189 . FR_RL(NSPMD+2,*),
190 . IPART(*), MONVOL(*), IPART_STATE(*),IFRAME(LISKN,*),
191 . IGEO(NPROPGI,*),IPM(NPROPMI,*),TAG_SKINS6(*),
192 . ICFIELD(SIZFIELD,*), LCFIELD(*), TAGSLV_RBY(*),
193 . MDS_OUTPUT_TABLE(MDS_NMAT*MAX_DEPVAR),MDS_NMAT,MAX_DEPVAR,
194 . MDS_NDEPSVAR(*),IBCL(*),ILOADP(*),LLOADP(*)
195 my_real
196 . SKEW(LSKEW,*),RBY(*),MS(*),IN(*),
197 . X(*), V(*), VR(*),WA(*),PLD(*),CRFLSW(*),XCUT(*),
198 . DAMPR(NRDAMP,*), PARTSAV(*),PM(NPROPM,*),
199 . VOLMON(*),GEO(NPROPG,*),XFRAME(NXFRAME,*)
200 TYPE(ttable) TABLE(*)
201 TYPE(ELBUF_STRUCT_), TARGET, DIMENSION(NGROUP) :: ELBUF_STR
202 TYPE(INTBUF_STRUCT_) INTBUF_TAB(*)
203 TYPE (H3D_DATABASE) :: H3D_DATA
204 TYPe (MULTI_FVM_STRUCT) :: MULTI_FVM
205 TYPE(GROUP_) ,DIMENSION(NGRNOD) :: IGRNOD
206 TYPE(GROUP_) ,DIMENSION(NGRPART) :: IGRPART
207 CHARACTER MDS_LABEL(1024,MDS_NMAT)
208 TYPE (STACK_PLY) :: STACK
209 TYPE (SENSORS_) ,INTENT(INOUT) ,TARGET :: SENSORS
210 TYPE (DYNAIN_DATABASE), INTENT(INOUT) :: DYNAIN_DATA
211 TYPE (DT_), INTENT(INOUT) :: DT
212 TYPE (LOADS_), INTENT(IN) :: LOADS
213 TYPE(output_), INTENT(INOUT) :: OUTPUT !< output structure
214 TYPE(NAMES_AND_TITLES_),INTENT(INOUT) :: NAMES_AND_TITLES !< NAMES_AND_TITLES host the input deck names and titles for outputs
215 TYPE (MATPARAM_STRUCT_) ,DIMENSION(NUMMAT) ,INTENT(INOUT) :: MAT_PARAM
216 type (glob_therm_) ,intent(inout) :: glob_therm
217 type (pblast_) ,intent(inout) :: pblast
218C-----------------------------------------------
219C L o c a l V a r i a b l e s
220C-----------------------------------------------
221 INTEGER NLEC(16), I, NSLIOF, NELSOF, NELQOF, NDAMPN,
222 . nelcof, neltof, nelpof, nelrof, nintch, nubcsn, itfor0, irfe0,
223 . irfl0, k, nelof, nn, nbc, il, ii, j, klg, kug, nblk, ng,
224 . ity, nel, nft, iad, igof, k1, k2, n, i1, i2, i3, ir1, ir2, ir3,
225 . isk, ic, icr, m1, m2, im, noint, nsearch, jpri, npts,
226 . nrbyon, nrbyof, neltgof,niniv,nsliofn,nsliofs,kk,nty,
227 . nsn,nrts,nrtm,nmn,l,ncpri1,nrlink0,nalelk0,
228 . nsphof,iwiout,nlecsph(10),nfxinp,ncrst1,neigoff,neoff,
229 . naleof,neulerof,nthermof,nfvmesh,mlw,ll,isk1,isk2,
230 . nfvmodi,iok, alestrl(0:8), nalelink0, gr_id,
231 . nalelk_starter, uid, nalelk_removed, jale,userl_count,
232 . idel2,is,
233 . stext1,multirest1,nlpri1,ierror,
234 . ninterskid,ivolu,ifv,ifv_type,istatg_p,
235 . k3,k4,k5,k6,k7,k8,k9,ninefricg,ns,ni,npartof,ift,ilt,k10,k11,jj,dtfct,ok
236
237 INTEGER :: NBPARTALEON, NBPARTALEOFF, PARTID, ISON, ION, IOFF, NBALEON_PART, NBALEOFF_PART
238 INTEGER, DIMENSION(:), ALLOCATABLE :: ALEOFF_PARTIDS_TMP, ALEOFF_PARTIDS,
239 . aleon_partids_tmp, aleon_partids
240 INTEGER, DIMENSION(:), ALLOCATABLE :: ALELIN_ON_OFF
241 my_real
242 . TFI, DTF, DTM, ALP, GAM,
243 . CV1, CV2, CV3, VOLM, TSTART, TFIN,DT_INPUT, DT_CRIT
244 CHARACTER*9 Label1,Label2,Label3,Label4,Label5,Label6,Label7,Label8
245 CHARACTER*16 ALEform(0:8),Cale1,Cale2,CHAR1,CHAR2,Label9
246 CHARACTER*3 :: LABEL_DEF,LABEL_ROT
247 CHARACTER(LEN=4), DIMENSION(20) :: TITLE
248 my_real eta1,eta2,MIN_ASPECT,MIN_DEFV,DYDX
249 my_real, DIMENSION(:), POINTER :: OFFG
250 TYPE (SENSOR_STR_) ,DIMENSION(:) ,POINTER :: SENSOR_TAB
251 LOGICAL IS_ALREADY_PRINTED
252 REAL(KIND=8) :: DTH
253 REAL(KIND=8), DIMENSION(9) :: DTH1
254 INTEGER, DIMENSION(:), ALLOCATABLE :: IPARTOF
255C-----------------------------------------------
256C E x t e r n a l F u n c t i o n s
257C-----------------------------------------------
258 INTEGER,EXTERNAL :: NGR2USR
259 my_real,EXTERNAL :: finter
260C-----------------------------------------------
261 ! setting NOINT in IPARI(13)
262 ierr = 0
263 dt_input= zero
264 DO i=1,ninter
265 ipari(13,i)=10
266 IF(ipari(15,i) == 0) ipari(15,i)=i
267 ENDDO
268 sensor_tab => sensors%SENSOR_TAB(1:sensors%NSENSOR)
269
270 CALL chkipari(ipari)
271
272 READ (iin,'(20A4)') title
273 READ (iin,'(4F16.0,I8,I10,2I8)')tfi,dth,dtf,dtm,ncpri1,ncrst1,multirest1,nlpri1
274 READ (iin,'(5F16.0)')(dth1(i),i=1,5)
275 READ (iin,'(4F16.0)')(dth1(i+5),i=1,4)
276 READ (iin,'(6F10.0)')alp,gam,cv1,cv2,cv3,volm
277 READ (iin,'(10I8)') nsliof,nelsof,nelqof,nelcof,neltof,nelpof,nelrof,neltgof,nsliofn,nsliofs
278 READ (iin,'(I8)') nsphof
279 READ (iin,'(4I8)') naleof,neulerof,nthermof,npartof
280 READ (iin,'(2I8)') nrlink0,nalelk0
281 READ (iin,'(2I8)') nalelink0
282C--------------
283C ALE ON / OFF
284C--------------
285 READ (iin, '(I8, I8)') nbpartaleon, nbpartaleoff
286 IF(mcheck == 0)THEN
287 nrlink=nrlink0
288 nalelk=nalelk0
289 nalelink=nalelink0
290 irprev = 0
291 ENDIF
292 nrbyon = 0
293 nrbyof = 0
294 nubcsn = 0
295 niniv = 0
296 READ (iin,'(2I8)') nubcsn,niniv
297 READ (iin,'(3I8)') nintch,nrbyon,nrbyof
298 READ (iin,'(I8)') ndampn
299 READ(iin,'(I8)') nfxinp
300C-----
301 READ(iin,'(2I8)') neigoff, neoff
302 READ(iin,'(2I8)') nfvmesh, nfvmodi
303C-----
304 IF(ncpri1 == 0) ncpri1=1
305 IF(mcheck == 0)ncpri = ncpri1
306C-----
307 IF(mcheck == 0)nlpri = nlpri1
308C
309 IF(ncrst1 == 0.AND.mcheck == 0) ncrst=10000000
310 IF(tfi /= 0.0) tstop=tfi
311C
312 IF(toutp0 /= zero) toutp = toutp0
313 IF(dtoutp0 > zero) dtoutp= dtoutp0
314 IF(dtoutp<=zero) toutp = ep30
315 IF (toutp < tt-dt2.AND.dtoutp > zero)toutp = toutp
316 . + int((tt-dt2-toutp)/dtoutp)*dtoutp
317 IF (toutp < tt-dt2)toutp = toutp+dtoutp
318C
319C .sta files
320 IF(tstat0 /= zero) tstat = tstat0
321 IF(dtstat0 > zero) dtstat= dtstat0
322 IF(dtstat<=zero) tstat = ep30
323 IF (tstat < tt-dt2.AND.dtstat > zero)tstat = tstat
324 . + int((tt-dt2-tstat)/dtstat)*dtstat
325 IF (tstat < tt-dt2)tstat = tstat+dtstat
326C .dynain files
327 IF(dynain_data%TDYNAIN0 /= zero) dynain_data%TDYNAIN = dynain_data%TDYNAIN0
328 IF(dynain_data%DTDYNAIN0 > zero) dynain_data%DTDYNAIN= dynain_data%DTDYNAIN0
329 IF(dynain_data%DTDYNAIN<=zero) dynain_data%TDYNAIN = ep30
330 IF (dynain_data%TDYNAIN < tt-dt2.AND.dynain_data%DTDYNAIN > zero)dynain_data%TDYNAIN = dynain_data%TDYNAIN
331 . + int((tt-dt2-dynain_data%TDYNAIN)/dynain_data%DTDYNAIN)*dynain_data%DTDYNAIN
332 IF (dynain_data%TDYNAIN< tt-dt2)dynain_data%TDYNAIN = dynain_data%TDYNAIN+dynain_data%DTDYNAIN
333C abf files
334 DO i=1,10
335 IF(dtabf0(i) > zero) dtabf(i)= dtabf0(i)
336 IF(dtabfwr0(i) > zero) dtabfwr(i)= dtabfwr0(i)
337 ENDDO
338 IF (abfile(1) /= 0) tabfis(1) = tt
339C
340 IF(dth /= zero.AND.mcheck == 0) output%TH%DTHIS=dth
341 DO i= 1, 9
342 IF(dth1(i) /= zero.AND.mcheck == 0) output%TH%DTHIS1(i)=dth1(i)
343 ENDDO
344
345C
346 IF(mcheck == 0)THEN
347 ale%GRID%VGY0=zero
348 ale%GRID%VGZ0=zero
349 IF(dtf /= zero) dtfac=dtf
350 IF(dtm /= zero) dtmin=dtm
351 IF(alp /= zero) ale%GRID%ALPHA=alp
352 IF(gam /= zero) ale%GRID%GAMMA=gam
353 IF(cv1 /= zero) ale%GRID%VGX =cv1
354 IF(cv2 /= zero) ale%GRID%VGY =cv2
355 IF(cv3 /= zero) ale%GRID%VGZ =cv3
356 IF(ale%GRID%NWALE_ENGINE == 2)THEN
357 IF(alp /= zero) THEN
358 ale%GRID%ALPHA = alp/(-ale%GRID%VGX+sqrt(ale%GRID%VGX**2+ one))
359 dt_input = alp
360 ELSE
361 dt_input = ale%GRID%ALPHA*(-ale%GRID%VGX+sqrt(ale%GRID%VGX**2+ one))
362 ENDIF
363 ENDIF
364 IF(ale%GRID%NWALE_ENGINE == 1.AND.alp == zero) ale%GRID%ALPHA=ep30
365 IF(volm /= zero)volmin=volm
366 IF(int22>0)THEN
367 dtfac22 = one
368 IF(dtfac>half)dtfac22 = half / dtfac
369 ENDIF
370 ENDIF
371C
372 IF(mcheck==0)THEN
373 IF(idtmins/=0)THEN
374 IF(dtmins == zero)dtmins = dtmin
375 IF(dtfacs == zero)dtfacs = dtfac
376 IF(tol_sms == zero) tol_sms = em03
377 IF(nsmspcg==0)nsmspcg=1000
378 idtgrs =idtgrs_old
379 IF(ispmd==0)THEN
380C
381 IF(idtmins==2.AND.irest_mselt==0)THEN
382 CALL ancmsg(msgid=120,anmode=aninfo_blind)
383 CALL arret(2)
384 END IF
385C
386C IF((IDTMINS==2.OR.IDTMINS_INT/=0).AND.NODADT/=0)THEN
387C CALL ANCMSG(MSGID=209,ANMODE=ANINFO_BLIND)
388C CALL ARRET(2)
389C END IF
390 END IF
391 ENDIF
392 IF(idtmins_int/=0)THEN
393 IF(dtmins_int == zero)dtmins_int = dtmin
394 IF(dtfacs_int == zero)dtfacs_int = dtfac
395 IF(tol_sms == zero) tol_sms = em03
396 IF(nsmspcg==0)nsmspcg=1000
397 ENDIF
398 ELSE
399 idtmins=idtmins_old
400 dtmins =dtmins_old
401 dtfacs =dtfacs_old
402 idtgrs =idtgrs_old
403 idtmins_int=idtmins_int_old
404 dtmins_int =dtmins_int_old
405 dtfacs_int =dtfacs_int_old
406 END IF
407C
408 dtfacx = one
409 DO i=1,51
410 IF(dtmin1(i) == zero)dtmin1(i) = dtmin
411 IF(dtfac1(i) == zero)dtfac1(i) = dtfac
412 dtfacx = min(dtfac1(i), dtfacx)
413 ENDDO
414
415 IF(idtmin(52) == 0)idtmin(52) = 1
416 IF(dtfac1(52) == zero)dtfac1(52) = zep9
417 IF(dtmin1(52) == zero)dtmin1(52) = em20
418
419 i=102
420 IF(ale%GLOBAL%IDT_ALE==-1)THEN
421 IF(dtmin1(i) == zero)dtmin1(i) = dtmin
422 IF(dtfac1(i) == zero)dtfac1(i) = dtfac
423 ELSE
424 IF(dtmin1(i) == zero)dtmin1(i) = zero
425 IF(dtfac1(i) == zero)dtfac1(i) = half
426 ENDIF
427C
428 IF(nodadt == 0)THEN
429 IF(idtmin(1) == 0)idtmin(1) = 1
430 IF(idtmin(2) == 0)idtmin(2) = 1
431 IF(idtmin(3) == 0)idtmin(3) = 2
432 IF(idtmin(7) == 0)idtmin(7) = 2
433 ENDIF
434C
435 IF(ispmd == 0) WRITE(iout,'(//1X,20A4//)') title
436 itfor0=itform + 1
437 irfe0 =(irform/5)
438 irfl0 =(irform-5*irfe0) + 1
439 irfe0 =irfe0 + 1
440
441C --------------------------------------------------------------------
442c CALL RCHECKMASS(IXR ,GEO ,PM ,MSR ,INR ,
443c . MS ,IN ,ITAB ,IGEO ,IPM ,
444c . UPARAM ,IPART ,IPARTR )
445C --------------------------------------------------------------------
446C User Libraries output & checks
447C --------------------------------------------------------------------
448 IF (nspmd > 1)THEN
449 userl_count = userl_avail
450 CALL spmd_glob_isum9(userl_count,1)
451 IF (ispmd==0)THEN
452 IF (userl_count /= 0 .AND. userl_count /= nspmd)THEN
453 CALL ancmsg(msgid=254,anmode=aninfo,
454 . c1=dlibfile(1:dlibfile_size))
455 CALL arret(2)
456 ENDIF
457 ENDIF
458 ENDIF
459C When Dynamical user libraries are used, add a print in 0001.out file
460 IF (userl_avail==1)THEN
461 IF(ispmd==0)THEN
462 WRITE(iout,4500)
463 WRITE(iout,4600) dlibfile(1:dlibfile_size),dlibtkvers
464 ENDIF
465 ENDIF
466C --------------------------------------------------------------------
467
468 IF (irad2r==1) THEN
469 IF(ispmd==0) WRITE(iout,1099)
470 IF (itfor0==1 ) THEN
471 CALL ancmsg(msgid=242,anmode=aninfo_blind)
472 CALL arret(2)
473 ENDIF
474 ELSE
475 IF ((r2r_siu==1).AND.(ispmd==0)) THEN
476 CALL ancmsg(msgid=239,anmode=aninfo_blind)
477 CALL arret(2)
478 ENDIF
479 ENDIF
480C
481 min_aspect = dt%BRICK_CST_COL_MIN
482 min_defv = dt%BRICK_CST_DEFV_MIN
483 IF(ispmd == 0.AND.mcheck == 0)WRITE(iout,1100)tstop,output%TH%DTHIS,dtfac,dtmin
484 IF(nodadt == 0)THEN
485 IF ((min_aspect+min_defv)>zero) THEN
486 IF(ispmd == 0.AND.mcheck == 0)THEN
487 WRITE(iout,4700) dtfac1(1),dtmin1(1),idtmin(1),
488 + min_aspect,min_defv,
489 + dtfac1(2),dtmin1(2),idtmin(2),
490 + dtfac1(3),dtmin1(3),idtmin(3),
491 + dtfac1(4),dtmin1(4),idtmin(4),
492 + dtfac1(5),dtmin1(5),idtmin(5),
493 + dtfac1(6),dtmin1(6),idtmin(6),
494 + dtfac1(9),dtmin1(9),idtmin(9),
495 + dtfac1(10),dtmin1(10),idtmin(10)
496 !WRITE(IOUT,4710) DTFAC1(102),DTMIN1(102)
497 IF(idt1sh /=0) WRITE(iout,4720)
498 IF(idt1sol/=0) WRITE(iout,4730)
499 IF(idt1tet10/=0) WRITE(iout,4740)
500 endif!(ISPMD == 0.AND.MCHECK == 0)
501 ELSE
502 IF(ispmd == 0.AND.mcheck == 0)THEN
503 WRITE(iout,1105) dtfac1(1),dtmin1(1),idtmin(1),
504 + dtfac1(2),dtmin1(2),idtmin(2),
505 + dtfac1(3),dtmin1(3),idtmin(3),
506 + dtfac1(4),dtmin1(4),idtmin(4),
507 + dtfac1(5),dtmin1(5),idtmin(5),
508 + dtfac1(6),dtmin1(6),idtmin(6),
509 + dtfac1(9),dtmin1(9),idtmin(9),
510 + dtfac1(10),dtmin1(10),idtmin(10)
511 !WRITE(IOUT,4710) DTFAC1(102),DTMIN1(102)
512 IF(idt1sh /=0) WRITE(iout,4720)
513 IF(idt1sol/=0) WRITE(iout,4730)
514 IF(idt1tet10/=0) WRITE(iout,4740)
515 endif!(ISPMD == 0.AND.MCHECK == 0)
516 END IF !(MIN_ASPECT+MIN_DEFV)>ZERO)
517 IF (dt%IDEL_BRICK>zero) THEN
518 IF(ispmd == 0.AND.mcheck == 0)THEN
519 WRITE(iout,5020)dt%BRICK_DEL_COL_MIN,dt%BRICK_DEL_DEFV_MIN,dt%BRICK_DEL_ASP_MAX,dt%BRICK_DEL_DEFV_MAX
520 ENDIF
521 ENDIF
522 IF(idtmin(11) == 3 .OR. idtmin(11) == 8) THEN
523 IF(ispmd == 0.AND.mcheck == 0)WRITE(iout,1107)dtfac1(11),dtmin1(11),idtmin(11)
524 ENDIF
525 IF(idtmins_int/=0)THEN
526 IF((ispmd == 0).AND.mcheck == 0)THEN
527 WRITE(iout,1209) dtfacs_int,dtmins_int
528 END IF
529 END IF
530 IF(idtmins/=0)THEN
531 IF((ispmd == 0).AND.mcheck == 0.AND.idtgrs<=0)THEN
532 WRITE(iout,1109) dtfacs,dtmins,tol_sms,nsmspcg,m_vs_sms,ncprisms,-idtgrs
533 ELSEIF((ispmd == 0).AND.mcheck == 0)THEN
534 WRITE(iout,1109) dtfacs,dtmins,tol_sms,nsmspcg,m_vs_sms,ncprisms,igrpart(idtgrs)%ID
535 END IF
536 END IF
537 ELSE
538 IF(ispmd == 0.AND.mcheck == 0)THEN
539 IF(idtmin(11)==0)THEN
540 WRITE(iout,1116) dtfac1(11),zero
541 ELSEIF(idtmin(11)==1.OR.idtmin(11)==3.OR.idtmin(11) == 8)THEN
542 IF ( percent_addmass > zero) THEN
543 WRITE(iout,1206) dtfac1(11),dtmin1(11),percent_addmass,idtmin(11)
544 ELSE
545 WRITE(iout,1106) dtfac1(11),dtmin1(11),idtmin(11)
546 ENDIF
547 END IF
548 END IF
549 IF ((min_aspect+min_defv)>zero) THEN
550 IF(ispmd == 0.AND.mcheck == 0)THEN
551 WRITE(iout,4700) dtfac1(1),dtmin1(1),idtmin(1),
552 + min_aspect,min_defv,
553 + dtfac1(2),dtmin1(2),idtmin(2),
554 + dtfac1(3),dtmin1(3),idtmin(3),
555 + dtfac1(4),dtmin1(4),idtmin(4),
556 + dtfac1(5),dtmin1(5),idtmin(5),
557 + dtfac1(6),dtmin1(6),idtmin(6),
558 + dtfac1(9),dtmin1(9),idtmin(9),
559 + dtfac1(10),dtmin1(10),idtmin(10)
560 !WRITE(IOUT,4710) DTFAC1(102),DTMIN1(102)
561 IF(idt1sh /=0) WRITE(iout,4720)
562 IF(idt1sol/=0) WRITE(iout,4730)
563 IF(idt1tet10/=0) WRITE(iout,4740)
564 ENDIF
565 ELSE
566 IF(ispmd == 0.AND.mcheck == 0)THEN
567 WRITE(iout,1105) dtfac1(1),dtmin1(1),idtmin(1),
568 + dtfac1(2),dtmin1(2),idtmin(2),
569 + dtfac1(3),dtmin1(3),idtmin(3),
570 + dtfac1(4),dtmin1(4),idtmin(4),
571 + dtfac1(5),dtmin1(5),idtmin(5),
572 + dtfac1(6),dtmin1(6),idtmin(6),
573 + dtfac1(9),dtmin1(9),idtmin(9),
574 + dtfac1(10),dtmin1(10),idtmin(10)
575 !WRITE(IOUT,4710) DTFAC1(102),DTMIN1(102)
576 IF(idt1sh /=0) WRITE(iout,4720)
577 IF(idt1sol/=0) WRITE(iout,4730)
578 IF(idt1tet10/=0) WRITE(iout,4740)
579 ENDIF
580 END IF !((MIN_ASPECT+MIN_DEFV)>ZERO) THEN
581 IF (dt%IDEL_BRICK>zero) THEN
582 IF(ispmd == 0.AND.mcheck == 0)THEN
583 WRITE(iout,5020)dt%BRICK_DEL_COL_MIN,dt%BRICK_DEL_DEFV_MIN,dt%BRICK_DEL_ASP_MAX,dt%BRICK_DEL_DEFV_MAX
584 ENDIF
585 END IF
586 IF(idtmins_int/=0)THEN
587 IF((ispmd == 0).AND.mcheck == 0)THEN
588 WRITE(iout,1209) dtfacs_int,dtmins_int
589 END IF
590 END IF
591C
592 IF(idtmins/=0)THEN
593C
594 IF (isms_selec < 3) THEN
595 IF((ispmd == 0).AND.mcheck == 0.AND.idtgrs<=0) THEN
596 WRITE(iout,1108) dtfacs,dtmins,tol_sms,nsmspcg,ncprisms,-idtgrs
597 ELSEIF((ispmd == 0).AND.mcheck == 0)THEN
598 WRITE(iout,1108) dtfacs,dtmins,tol_sms,nsmspcg,ncprisms,igrpart(idtgrs)%ID
599 END IF
600C
601 ELSE
602C- /DT /CST_AMS>- Automatic election
603 dt_crit = dtmins / max(em20,dtfac1(11))
604 IF((ispmd == 0).AND.mcheck == 0.AND.idtgrs<=0)THEN
605 WRITE(iout,2109) dtfacs,dtmins,tol_sms,nsmspcg,m_vs_sms,ncprisms,dt_crit,-idtgrs
606 ELSEIF((ispmd == 0).AND.mcheck == 0)THEN
607 WRITE(iout,2109) dtfacs,dtmins,tol_sms, nsmspcg,m_vs_sms,ncprisms,dt_crit,igrpart(idtgrs)%ID
608 END IF
609C
610 END IF
611C
612 ENDIF
613C
614 ENDIF
615C
616
617 IF(nodadt==0.AND.(istatcnd/=0.AND.impl_s==0))THEN
618 CALL ancmsg(msgid=121,anmode=aninfo_blind)
619 CALL arret(2)
620 END IF
621 IF(ispmd == 0.AND.mcheck == 0.AND.kdtint /= 0)THEN
622 WRITE(iout,'(A)')' NEW (HIDDEN) TIME STEP COMPUTATION',' ON INTERFACE TYPE 7,11 AND 19 IS ON'
623 ENDIF
624 IF(ispmd == 0.AND.mcheck == 0.AND.kdtsmstr == 0)THEN
625 WRITE(iout,'(A)')' BACK TO VERSION 4 COMPUTATION OF NODAL TIME STEP',' CASE OF SMALL STRAIN FOR SOLIDS.'
626 ENDIF
627 IF(codvers>=44.AND.numsph /= 0.AND.mcheck == 0)THEN
628 IF(idtmin(51) == 3)THEN
629 IF(ispmd == 0)THEN
630 WRITE(istdo,*)' ** WARNING SMALL STRAIN FORMULATION FOR SPH'
631 WRITE(iout,*)' ** WARNING SMALL STRAIN FORMULATION IS NOT AVAILABLE FOR SPH,',' OPTION /DT/SPHCEL/CST WILL BE OMITTED.'
632 END IF
633 idtmin(51)=0
634 ENDIF
635 IF(nodadt == 1)THEN
636 IF(ispmd == 0) WRITE(iout,*)'IMPROVED TIME STEP (NODAL) COMPUTATION ON SPH PARTICLES :'
637 dtfac1(51)=dtfac1(11)
638 ENDIF
639 IF(ispmd == 0) WRITE(iout,1151)dtfac1(51),dtmin1(51),idtmin(51)
640 ENDIF
641C
642 IF(ispmd == 0.AND.mcheck == 0) THEN
643 IF(num_opt_dt_fvmbag_0+num_opt_dt_fvmbag_1>0)THEN
644 !if at least on option /DT/FVMBAG/0 or /DT/FVMBAG/1 is defined then output parameter for the retained one (last read one)
645 ! (multiple definition is not expected)
646 IF(num_opt_dt_fvmbag_1>0)THEN
647 WRITE(iout,1156)
648 ELSE
649 WRITE(iout,1155)
650 ENDIF
651 WRITE(iout,1152)dtfac1(52),dtmin1(52),idtmin(52)
652 ENDIF
653 ENDIF
654 IF(ispmd == 0.AND.mcheck == 0) THEN
655 k1 = 1
656 is_already_printed = .false.
657 DO ivolu = 1, nvolu
658 ifv = monvol(k1 - 1 + 45)
659 ifv_type = monvol(k1 - 1 + 2)
660 IF (ifv_type == 6 .OR. ifv_type == 8 .OR. ifv_type == 11) THEN
661 IF(.NOT.is_already_printed)THEN
662 WRITE(iout,1157)
663 is_already_printed=.true.
664 ENDIF
665 ENDIF
666 IF(ifv_type == 6) WRITE(iout,1147)monvol(k1)
667 IF(ifv_type == 8) WRITE(iout,1148)monvol(k1)
668 IF(ifv_type == 11)WRITE(iout,1149)monvol(k1)
669 IF (ifv_type == 6 .OR. ifv_type == 8 .OR. ifv_type == 11) THEN
670 WRITE(iout,1153)fvdata_old(ifv)%CFL_COEF,fvdata_old(ifv)%DTMIN,fvdata_old(ifv)%ID_DT_OPTION
671 IF(fvdata_old(ifv)%ID_DT_OPTION ==2)THEN
672 WRITE(iout,1154)fvdata_old(ifv)%L_TYPE,fvdata_old(ifv)%LAMBDA,fvdata_old(ifv)%DTOLD
673 ENDIF
674 ENDIF
675 k1=k1+nimv
676 ENDDO
677 ENDIF
678C
679 IF(ispmd == 0.AND.mcheck == 0.AND.idt1sol /= 0) WRITE(iout,'(A)')' OPTIMIZED TIME STEP COMPUTATION FOR HEPH.'
680 IF(ispmd == 0.AND.mcheck == 0.AND.idttsh /= 0) WRITE(iout,'(A)')' OPTIMIZED (SHELL) TIME STEP FOR THICK SHELLS.'
681 IF(ispmd == 0.AND.mcheck == 0)WRITE(iout,1110)ncpri,ncrst,multirest,invers,itfor0
682 IF(ispmd == 0.AND.restart_file == 0) WRITE(iout,*) check_message(2)(1:len_trim(check_message(2)) )
683 IF (impl_s == 1.AND.ikg>=5) THEN
684 ikg=ikg-5
685 IF(ispmd == 0) THEN
686 WRITE(iout,*) ' ** WARNING ** : PARITH/ON IS NOT ','COMPATIBLE WITH IMPLICIT OPTION '
687 WRITE(iout,*) ' ** RESETTING ** : PARITH/OFF '
688 ENDIF
689 ENDIF
690C---------------------
691C FE THERMAL ANALYSIS
692C---------------------
693 IF(ispmd == 0 .AND. imassi > 0 ) THEN
694 WRITE(iout,'(A)')' RESET INITIAL MASS FOR THIS RUN '
695 ENDIF
696C---------------------
697C INITIAL MASS
698C---------------------
699 IF (ispmd == 0 .AND. glob_therm%ITHERM_FE > 0 ) THEN
700 WRITE(iout,'(A)')' FINITE ELEMENT THERMAL ANALYSIS '
701 IF (glob_therm%IDT_THERM == 1) THEN
702 WRITE(iout,'(A)')' THERMAL ANALYSIS ONLY'
703 WRITE(iout,'(A,1X,G14.7)')' THERMAL TIME STEP SCALE FACTOR . . . . . . . . .', glob_therm%DTFACTHERM
704 ENDIF
705 IF (glob_therm%THEACCFACT > one ) THEN
706 WRITE(iout,'(A,1X,G14.7/)')' FACTOR TO SPEED-UP THERMAL ANALYSIS. . . . . . .', glob_therm%THEACCFACT
707 ENDIF
708 ENDIF
709C
710C-----------------------
711C The correct indexes in ANIM data structures wrt integration points or layers
712C of all actual elements will be rebuilt ::
713C-----------------------
714 CALL anim_build_index_all(ispmd ,mcheck ,sensors ,igeo ,geo )
715C-------------------------------------------
716C STOP/SENSOR
717C-------------------------------------------
718 IF (mcheck == 0) THEN
719 ALLOCATE(sensors%STOP(sensors%NSTOP))
720 sensors%STOP(:) = 0
721 END IF
722 IF (mcheck == 0) THEN
723 IF (sensors%NSTOP > 0) THEN
724 DO k=1,sensors%NSTOP
725 ierr = 1
726 IF (sensors%STOP_TMP(k) > 0) THEN
727 DO i=1,sensors%NSENSOR
728 IF (sensor_tab(i)%SENS_ID == sensors%STOP_TMP(k)) THEN
729 sensors%STOP(k) = i
730 ierr = 0
731 EXIT
732 ENDIF
733 ENDDO
734 ENDIF
735 IF (ierr == 1) THEN
736 CALL ancmsg(msgid=233, anmode=aninfo,i1=sensors%STOP_TMP(k))
737 ENDIF
738 ENDDO
739 ENDIF
740 END IF !(ISPMD==0.AND.MCHECK==0)
741C-------------------------------------------
742C STATE/SENSOR
743C-------------------------------------------
744 IF (mcheck == 0) THEN
745 ALLOCATE(sensors%STAT(sensors%NSTAT))
746 sensors%STAT(:) = 0
747 END IF
748 IF (ispmd == 0 .AND. mcheck == 0) THEN
749 mstatt = 0
750 IF (sensors%NSTAT > 0) THEN
751 DO k=1,sensors%NSTAT
752 ierr = 1
753 IF(sensors%STAT(k) /= 0) THEN
754 DO i=1,sensors%NSENSOR
755 IF (sensors%SENSOR_TAB(i)%SENS_ID == sensors%STAT(k)) THEN
756 sensors%STAT(k) = i
757 ierr = 0
758 EXIT
759 ENDIF
760 ENDDO
761 ENDIF
762 IF (ierr == 1) THEN
763 CALL ancmsg(msgid=235, anmode=aninfo,i1=sensors%STAT(k))
764 ELSE
765 mstat(k) = 0
766 ENDIF
767 ENDDO
768 ENDIF
769 END IF !(ISPMD==0.AND.MCHECK==0)
770C-------------------------------------------
771C OUTP/SENSOR
772C-------------------------------------------
773 IF (mcheck == 0) THEN
774 ALLOCATE(sensors%OUTP(sensors%NOUTP))
775 sensors%OUTP(:) = 0
776 END IF
777 IF (ispmd == 0 .AND. mcheck == 0) THEN
778 IF (sensors%NOUTP > 0) THEN
779 DO k=1,sensors%NOUTP
780 ierr = 1
781 IF (sensors%OUTP_TMP(k) > 0) THEN
782 DO i=1,sensors%NSENSOR
783 IF (sensors%SENSOR_TAB(i)%SENS_ID == sensors%OUTP_TMP(k)) THEN
784 sensors%OUTP(k) = i
785 ierr = 0
786 EXIT
787 ENDIF
788 ENDDO
789 ENDIF
790 IF (ierr == 1) THEN
791 CALL ancmsg(msgid=236, anmode=aninfo,i1 = sensors%OUTP_TMP(k))
792 ENDIF
793 ENDDO
794 ENDIF
795 END IF !(ISPMD==0.AND.MCHECK==0)
796
797C-------------------------------------------
798C STATE/SENSOR
799C-------------------------------------------
800 IF(ispmd == 0.AND.mcheck == 0) THEN
801 WRITE(iout,1150)dtin,dtmx
802 IF(impl /= 0)
803 . WRITE(iout,1160)eps,eps2,nitmx
804 istatg_p = iabs(istatg)
805 IF(istat == 1)THEN
806 WRITE(iout,1171)istatg_p,beta,period
807 ELSEIF(istat == 2)THEN
808 IF (tst_start>zero.OR.tst_stop>zero) THEN
809 IF (tst_stop==zero) tst_stop = tstop
810 WRITE(iout,5001)istatg_p,tst_start,tst_stop
811 ELSE
812 WRITE(iout,1172)istatg_p
813 END IF
814 ELSEIF(istat == 3)THEN
815 IF (tst_start>zero.OR.tst_stop>zero) THEN
816 IF (tst_stop==zero) tst_stop = tstop
817 WRITE(iout,5011)istatg_p,tst_start,tst_stop
818 ELSE
819 WRITE(iout,5010)istatg_p
820 END IF
821 ENDIF
822C-----------------------
823 IF(iparit == 0)THEN
824 WRITE(iout,1180)
825 ELSEIF(iparit == 1)THEN
826 WRITE(iout,1181)
827 ELSE
828 WRITE(iout,1182)iparit-1
829 ENDIF
830C-----------------------
831 WRITE(iout,1300) nsliof,npartof,nelsof,nelqof,nelcof,neltof,nelpof,nelrof,neltgof,nsphof
832 WRITE(iout,1400) nrlink
833 WRITE(iout,1500) nubcsn
834 IF(iale /= 0.OR.ieuler /= 0)WRITE(iout,1450) nalelink+nalelk !format v12 + format v5
835c
836 ELSEIF(ispmd == 0 .AND. mcheck /= 0) THEN
837 ENDIF !(ISPMD==0.AND.MCHECK==0)
838C-----------------------------------
839C ALE/EULER : UPWIND
840C----------------------------------------
841 ! ALE%UPWIND%UPW_UPDATE == 1 : ENGINE /UPWIND CARD DETECTED
842 !ALE%UPWIND%UPW_UPDATE == 2 : /UPWIND CARD IS CHANGING AT LEAST ONE PARAMETER
843 !if /UPWIND is defined then set new parameters
844 IF(ispmd == 0 .AND. mcheck == 0)THEN
845 IF(iale /= 0 .OR. ieuler /= 0)THEN
846 IF(ale%UPWIND%UPW_UPDATE /= 0 )THEN
847 DO k=1,nummat-1
848 !ILAW=PM(19,K)
849 jale=int(pm(72,k)) ! JALE from PROP (IGEO(62,IPID)) is not checked since this option /UPWIND is obsolete
850 IF(jale /= 0)THEN
851 IF(pm(15,k) /= ale%UPWIND%UPWMG2 .OR. pm(16,k) /= ale%UPWIND%UPWOG2)THEN
852 ale%UPWIND%UPW_UPDATE = 2
853 pm(15,k) = ale%UPWIND%UPWMG2
854 pm(16,k) = ale%UPWIND%UPWOG2
855 ENDIF
856 ENDIF
857 ENDDO
858 IF(ale%UPWIND%UPWSM /= ale%UPWIND%UPWSM2)THEN
859 ale%UPWIND%UPW_UPDATE = 2
860 ale%UPWIND%UPWSM = ale%UPWIND%UPWSM2
861 ENDIF
862 !if /UPWIND is not defined then catch parameter from previous run
863 ELSEIF(ale%UPWIND%UPW_UPDATE == 0)THEN
864 ale%UPWIND%UPWMG2=one
865 ale%UPWIND%UPWOG2=one
866 ale%UPWIND%UPWSM2=one
867 DO k=1,nummat-1
868 !ILAW=PM(19,K)
869 jale=int(pm(72,k))
870 IF(jale /= 0)THEN
871 IF(pm(15,k) /= ale%UPWIND%UPWMG2 .OR. pm(16,k) /= ale%UPWIND%UPWOG2 .OR. ale%UPWIND%UPWSM /= ale%UPWIND%UPWSM2)THEN
872 IF(pm(15,k) /= zero)ale%UPWIND%UPWMG2 = pm(15,k)
873 IF(pm(16,k) /= zero)ale%UPWIND%UPWOG2 = pm(16,k)
874 IF(ale%UPWIND%UPWSM /= zero)ale%UPWIND%UPWSM2 = ale%UPWIND%UPWSM
875 EXIT
876 ENDIF
877 ENDIF
878 ENDDO
879 ENDIF
880 ENDIF
881 ENDIF
882C-----------------------------------
883C ALE/EULER : OUTPUT
884C----------------------------------------
885 IF(iale /= 0 .OR. ieuler /= 0)THEN
886 !-----------------------
887 ! QUASI-INCOMPRESSIBLE /INCMP
888 !-----------------------
889 IF(ale%GLOBAL%INCOMP == 1.AND.mcheck == 0) WRITE(iout,1196)
890 !-----------------------
891 ! GRID FORMULATION
892 !-----------------------
893 !---ALE GRID FORMULATION CHECK---!
894 IF(ale%GRID%NWALE_ENGINE /= ale%GRID%NWALE_RST .AND. ale%GRID%NWALE_ENGINE /= -1) THEN !Forbid ALE grid formulation switch
895 IF(ale%GRID%NWALE_ENGINE /= 3)THEN !expect if new one is /ALE/ZERO
896 aleform=(/'DONEA ','DISP ','SPRING ', 'ZERO ',
897 . 'STANDARD ','LAPLACIAN ','VOLUME ', 'FLOW-TRACKING ',
898 . 'LAGRANGE ' /)
899 alestrl=(/5,4,6,4,8,9,6,13,8/)
900 cale1=aleform(ale%GRID%NWALE_ENGINE)
901 cale2=aleform(ale%GRID%NWALE_RST)
902 CALL ancmsg(msgid=229,anmode=aninfo,
903 . c1=cale1(1:alestrl(ale%GRID%NWALE_ENGINE)),c2=cale2(1:alestrl(ale%GRID%NWALE_RST)))
904 CALL arret(2)
905 ELSE
906 ale%GRID%NWALE_ENGINE = 3
907 END IF
908 END IF
909
910 !--Labels for Staggered Scheme
911 IF(ale%UPWIND%UPWM == 2)THEN
912 label1='TG '
913 ELSEIF(ale%UPWIND%UPWM == 3)THEN
914 label1='SUPG '
915 ELSE
916 label1='UPWIND '
917 ENDIF
918 label2='UPWIND '
919 label3='UPWIND '
920 IF(alemuscl_param%IALEMUSCL == 0)THEN
921 label4='UPWIND '
922 ELSE
923 label4='MUSCL '
924 ENDIF
925 !--Labels for Colocated Scheme
926 IF(alemuscl_param%IALEMUSCL == 0)THEN
927 label5='1ST-ORDER'
928 label6='1ST-ORDER'
929 label7='1ST-ORDER'
930 label8='1ST-ORDER'
931 ELSEIF(alemuscl_param%IALEMUSCL == 1)THEN
932 label5='2ND-ORDER'
933 label6='2ND-ORDER'
934 label7='2ND-ORDER'
935 label8='2ND-ORDER'
936 ELSEIF(alemuscl_param%IALEMUSCL == 2)THEN
937 label5='1ST-ORDER'
938 label6='1ST-ORDER'
939 label7='1ST-ORDER'
940 label8='2ND-ORDER'
941 ENDIF
942 multi_fvm%LOWMACH_OPT = .false.
943 IF (ale%GLOBAL%HLLC_LOWMACH == 1 .AND. multi_fvm%IS_USED) THEN
944 multi_fvm%LOWMACH_OPT = .true.
945 ENDIF
946 aleform=(/'DONEA ','DISP ','SPRING ', 'ZERO ',
947 . 'STANDARD ','LAPLACIAN ','VOLUME ', 'FLOW-TRACKING ',
948 . 'LAGRANGE '/)
949 !effective value which will be retained during numerical solving
950 IF(ale%GRID%NWALE_ENGINE /= -1)THEN
951 ale%GRID%NWALE = ale%GRID%NWALE_ENGINE
952 ELSE
953 ale%GRID%NWALE = ale%GRID%NWALE_RST
954 ENDIF
955 label9(:)=' '
956 label9(1:16)=aleform(ale%GRID%NWALE)
957 eta1=ale%UPWIND%UPWMG2
958 eta2=ale%UPWIND%UPWOG2
959 IF(ale%UPWIND%UPWM>1)eta1=ale%UPWIND%CUPWM
960 !OUTPUT NUMERICAL SCHEME AND ITS PARAMETERS
961 IF(ispmd == 0 .AND. mcheck == 0)THEN
962 WRITE(iout,1001)
963 WRITE(iout,1002)label1,eta1,label2,eta2,label3,eta2,label4
964 IF(multi_fvm%IS_USED)THEN
965 WRITE(iout,1003)label5,label6,label7,label8
966 IF(multi_fvm%LOWMACH_OPT)WRITE(iout,1004)
967 IF(alemuscl_param%IALEMUSCL/=0)WRITE(iout,1005)alemuscl_param%BETA
968 ENDIF
969 !OUTPUT COURANT NUMBER
970 WRITE(iout,1006) dtfac1(102),dtmin1(102)
971 !OUTPUT GRID SMOOTHING FORMULATION AND ITS PARAMETERS
972 WRITE(iout,1007)label9(1:len_trim(label9))
973 SELECT CASE (ale%GRID%NWALE)
974 CASE(0) !DONEA
975 WRITE(iout,1008)
976 WRITE(iout,1200) ale%GRID%ALPHA,ale%GRID%GAMMA,ale%GRID%VGX,ale%GRID%VGY,ale%GRID%VGZ,volmin
977 CASE(1) !DISP
978 WRITE(iout,1008)
979 WRITE(iout,1220) ale%GRID%ALPHA,volmin
980 CASE(2) !SPRING
981 WRITE(iout,1008)
982 WRITE(iout,1250) dt_input,ale%GRID%ALPHA,ale%GRID%GAMMA,ale%GRID%VGX,ale%GRID%VGY,volmin
983 CASE(3) !ZERO
984 !no parameters
985 CASE(4) !STANDARD
986 WRITE(iout,1008)
987 WRITE(iout,1254) ale%GRID%ALPHA,ale%GRID%GAMMA,ale%GRID%VGX,ale%GRID%VGY
988 CASE(5)! LAPLACIAN
989 WRITE(iout,1008)
990 WRITE(iout,1254) ale%GRID%ALPHA,ale%GRID%GAMMA,ale%GRID%VGX,ale%GRID%VGY
991 CASE(6) !VOLUME
992 !no parameters
993 CASE(7) !FLOW-TRACKING
994 WRITE(iout,1008)
995 label_def=' NO'
996 label_rot=' NO'
997 IF(int(ale%GRID%VGX) == 1)label_def='YES'
998 IF(int(ale%GRID%VGY) == 1)label_rot='YES'
999 WRITE(iout,1257) label_def,label_rot,ale%GRID%ALPHA,ale%GRID%GAMMA
1000 CASE(8) !LAGRANGE
1001 !no parameters
1002 END SELECT
1003 ENDIF
1004 ENDIF !(IALE/=0.OR.IEULER/=0)
1005 !---CHECK END-------------------!
1006!-----------------------
1007! FEM MOMENTUM INTEGRATION
1008!-----------------------
1009 IF(ispmd == 0.AND.mcheck == 0)THEN
1010 IF(iale+ieuler /= 0 .AND. ale%GLOBAL%ISFINT /= 3)THEN
1011 SELECT CASE(ale%GLOBAL%ISFINT)
1012 !CASE(3)
1013 ! WRITE(IOUT,1197) !default
1014 CASE(2)
1015 WRITE(iout,1198)
1016 CASE(1)
1017 WRITE(iout,1199)
1018 END SELECT
1019 ENDIF
1020 ENDIF
1021
1022C-------------------------------
1023C ELIMINATION PARTIE ALE EULER THERMIQUE
1024C-------------------------------
1025 IF(iale /= 0.AND.naleof == 1)THEN
1026 IF(ispmd == 0.AND.mcheck == 0)THEN
1027 WRITE(iout,*)'ALE FORMULATION SWITCHED OFF '
1028 END IF
1029 iale = 0
1030 END IF
1031 IF(ieuler /= 0.AND.neulerof == 1)THEN
1032 IF(ispmd == 0.AND.mcheck == 0)THEN
1033 WRITE(iout,*)'EULER FORMULATION SWITCHED OFF '
1034 END IF
1035 ieuler = 0
1036 END IF
1037 IF (glob_therm%ITHERM /= 0.AND.nthermof == 1)THEN
1038 IF(ispmd == 0.AND.mcheck == 0)THEN
1039 WRITE(iout,*)'THERMIC FORMULATION SWITCHED OFF '
1040 END IF
1041 glob_therm%ITHERM = 0
1042 END IF
1043C-------------------------------
1044 IF(idel7 /= 0.AND.mcheck == 0.AND.n2d == 0)THEN
1045 ideli7 = idel7 - 1
1046 IF(ispmd == 0)THEN
1047 WRITE(iout,1550) ideli7
1048 CALL ancmsg(msgid=122,anmode=aninfo_blind)
1049 CALL arret(2)
1050 ENDIF
1051 ELSEIF(n2d /= 0.AND.mcheck == 0)THEN
1052 IF(ispmd == 0)WRITE(iout,1550) idel7-1
1053 END IF
1054C-------------------------------
1055C INTERFACE CLEANING
1056C-------------------------------
1057 IF(nsliof /= 0) THEN
1058C
1059 nn=(nsliof+9)/10
1060 DO il=1,nn
1061 READ (iin,'(10I10)') (nlec(i),i=1,10)
1062 DO i=1,10
1063 IF(nlec(i) == 0)cycle
1064 IF(ispmd == 0.AND.mcheck == 0)WRITE(iout,2000) nlec(i)
1065 DO k=1,ninter
1066 IF(ipari(15,k) == nlec(i))THEN
1067 IF(ipari(7,k)==2)THEN
1068 idel2= ipari(17,k)
1069 IF(idel2 /= 0)THEN
1070C SMASS & SINER have been saved only if IDEL2 /= 0
1071 nsn = ipari(5,k)
1072 DO n=1,nsn
1073 is =intbuf_tab(k)%NSV(n)
1074 IF(is > 0)THEN
1075 ms(is)=intbuf_tab(k)%SMAS(n)
1076 in(is)=intbuf_tab(k)%SINER(n)
1077 intbuf_tab(k)%NSV(n)=-is
1078 END IF
1079 END DO
1080 END IF
1081 END IF
1082C----- change T_stop for int25 otherwise issue w/ int25
1083 IF(ipari(7,k)==25)THEN
1084 intbuf_tab(k)%VARIABLES(11) = zero
1085 ELSE
1086 ipari(7,k)=0
1087 ENDIF
1088 END IF
1089 ENDDO
1090 ENDDO !next I
1091 ENDDO
1092 ENDIF
1093C-----------------------------------
1094C CLEANING INTERFACE NODES
1095C----------------------------------------
1096 IF(nsliofn /= 0) THEN
1097 READ (iin,'(2I8)') i,nn
1098 kk=0
1099 DO k=1,ninter
1100 IF(ipari(15,k) == i)kk=k
1101 ENDDO
1102 IF(kk == 0)THEN
1103 IF(ispmd == 0) THEN
1104 CALL ancmsg(msgid=123,anmode=aninfo_blind)
1105 ENDIF
1106 CALL arret(2)
1107 ENDIF
1108 nrts =ipari(3,kk)
1109 nrtm =ipari(3,kk)
1110 nsn =ipari(5,kk)
1111 nmn =ipari(6,kk)
1112 nty =ipari(7,kk)
1113 IF (ipari(71,kk)/=0) nty = 19 ! Interface type 19
1114 IF(nty == 3)THEN
1115 nsn =nsn + nmn
1116 ELSEIF(nty == 4 .OR. nty == 5)THEN
1117 ELSEIF(nty == 7 .OR. nty == 10 .OR. nty == 22)THEN
1118 ELSE
1119 IF(ispmd == 0) THEN
1120 CALL ancmsg(msgid=124,anmode=aninfo_blind)
1121 ENDIF
1122 CALL arret(2)
1123 ENDIF
1124 DO il=1,nn
1125 READ (iin,'(10I10)') (nlec(i),i=1,10)
1126 DO i=1,10
1127 IF(nlec(i) /= 0)THEN
1128 DO k=1,nsn
1129 IF(itab(intbuf_tab(kk)%NSV(k)) == nlec(i))THEN !NSV(K)
1130 intbuf_tab(kk)%STFNS(k)=zero !STFNS(K)
1131 ENDIF
1132 ENDDO
1133 ENDIF
1134 ENDDO
1135 ENDDO
1136 ENDIF
1137C----------------------------------------------
1138C CLEANING ELEMS
1139C----------------------------------------------
1140 nelof=nelsof+nelqof+nelcof+neltof+nelpof+nelrof+neltgof+nsphof
1141C-----------------------
1142C 1. 3D ELEMS
1143C-----------------------
1144 IF(nelsof > 0)THEN
1145 IF(ispmd == 0.AND.mcheck == 0)WRITE (iout,1810)
1146 nn=(nelsof+4)/5
1147 nbc=5
1148 DO il=1,nn
1149 READ (iin,'(10I10)') (nlec(ii),ii=1,10)
1150 DO j=1,nbc
1151 klg=nlec(2*j-1)
1152 kug=nlec(2*j )
1153 IF(klg<=0) GO TO 120
1154 nblk=(il-1)*nbc+j
1155 DO l=1,numels
1156 IF(ixs(nixs,l)>=klg.AND.ixs(nixs,l)<=kug) THEN
1157 DO k=klg,kug
1158 IF(ixs(nixs,l) == k) THEN
1159 ixs(1,l)=-iabs(ixs(1,l))
1160 GOTO 111
1161 ENDIF
1162 ENDDO
1163 ENDIF
1164111 CONTINUE
1165 ENDDO
1166 ENDDO
1167 ENDDO
1168 120 CONTINUE
1169 k=0
1170 DO j=1,numels
1171 IF(ixs(1,j) < 0)THEN
1172 k=k+1
1173 nlec(k) = ixs(nixs,j)
1174 IF(k == 10)THEN
1175 IF(mcheck == 0)WRITE (iout,'(5I10)') (nlec(ii),ii=1,k)
1176 k=0
1177 ENDIF
1178 ENDIF
1179 ENDDO
1180 IF(k > 0) THEN
1181 IF(mcheck == 0)WRITE (iout,'(5I10)') (nlec(ii),ii=1,k)
1182 ENDIF
1183 IF(nspmd > 1) THEN
1184 ! required treatments to retrieve deleted elems in correct order
1185 iwiout = 0
1186 IF (ispmd /= 0) CALL spmd_chkw(iwiout,iout)
1187 CALL spmd_glob_isum9(iwiout,1)
1188 CALL spmd_ibcast(iwiout,iwiout,1,1,0,2)
1189 IF (iwiout > 0) THEN
1190 CALL spmd_wiout(iout,iwiout)
1191 iwiout = 0
1192 ENDIF
1193 ENDIF
1194 endif!(NELSOF > 0)
1195C-----------------------
1196C 2. 2D-ELEMS
1197C-----------------------
1198 IF(nelqof > 0)THEN
1199 IF(ispmd == 0.AND.mcheck == 0)WRITE (iout,1820)
1200 IF (invers < 18) THEN
1201 nn=(nelqof+7)/8
1202 nbc=8
1203 ELSE
1204 nn=(nelqof+4)/5
1205 nbc=5
1206 ENDIF
1207 DO i=1,nn
1208 READ (iin,'(10I10)') (nlec(ii),ii=1,10)
1209 DO j=1,nbc
1210 klg=nlec(2*j-1)
1211 kug=nlec(2*j )
1212 IF(klg<=0) GO TO 220
1213 nblk=(i-1)*nbc+j
1214 DO l=1,numelq
1215 IF(ixq(nixq,l)>=klg.AND.ixq(nixq,l)<=kug) THEN
1216 DO k=klg,kug
1217 IF(ixq(nixq,l) == k) THEN
1218 ixq(1,l)=-iabs(ixq(1,l))
1219 GOTO 211
1220 ENDIF
1221 ENDDO
1222 ENDIF
1223 211 CONTINUE
1224 ENDDO
1225 ENDDO
1226 ENDDO
1227 220 CONTINUE
1228 k=0
1229 DO j=1,numelq
1230 IF(ixq(1,j) < 0)THEN
1231 k=k+1
1232 nlec(k) = ixq(nixq,j)
1233 IF(k == 10)THEN
1234 IF(mcheck == 0)WRITE (iout,'(5I10)') (nlec(ii),ii=1,k)
1235 k=0
1236 ENDIF
1237 ENDIF
1238 ENDDO
1239 IF(k > 0) THEN
1240 IF(mcheck == 0)WRITE (iout,'(5I10)') (nlec(ii),ii=1,k)
1241 ENDIF
1242 IF(nspmd > 1) THEN
1243 ! required treatments to retrieve deleted elems in correct order
1244 iwiout = 0
1245 IF (ispmd /= 0) CALL spmd_chkw(iwiout,iout)
1246 CALL spmd_glob_isum9(iwiout,1)
1247 CALL spmd_ibcast(iwiout,iwiout,1,1,0,2)
1248 IF (iwiout > 0) THEN
1249 CALL spmd_wiout(iout,iwiout)
1250 iwiout = 0
1251 ENDIF
1252 ENDIF
1253 END IF
1254C-----------------------
1255C 3. SHELL ELEMS
1256C-----------------------
1257 IF(nelcof > 0)THEN
1258 IF(ispmd == 0.AND.mcheck == 0)WRITE (iout,1830)
1259 nn=(nelcof+4)/5
1260 nbc=5
1261 DO i=1,nn
1262 READ (iin,'(10I10)') (nlec(ii),ii=1,10)
1263 DO j=1,nbc
1264 klg=nlec(2*j-1)
1265 kug=nlec(2*j )
1266 IF(klg<=0) GO TO 320
1267 nblk=(i-1)*nbc+j
1268 DO l=1,numelc
1269 IF(ixc(nixc,l)>=klg.AND.ixc(nixc,l)<=kug) THEN
1270 DO k=klg,kug
1271 IF(ixc(nixc,l) == k) THEN
1272 ixc(1,l)=-iabs(ixc(1,l))
1273 GOTO 311
1274 ENDIF
1275 ENDDO
1276 ENDIF
1277 311 CONTINUE
1278 ENDDO
1279 ENDDO
1280 ENDDO
1281 320 CONTINUE
1282 k=0
1283 DO j=1,numelc
1284 IF(ixc(1,j) < 0)THEN
1285 k=k+1
1286 nlec(k) = ixc(nixc,j)
1287 IF(k == 10)THEN
1288 IF(mcheck == 0)WRITE (iout,'(5I10)') (nlec(ii),ii=1,k)
1289 k=0
1290 ENDIF
1291 ENDIF
1292 ENDDO
1293 IF(k > 0) THEN
1294 IF(mcheck == 0)WRITE (iout,'(5I10)') (nlec(ii),ii=1,k)
1295 ENDIF
1296 IF(nspmd > 1) THEN
1297 ! required treatments to retrieve deleted elems in correct order
1298 iwiout = 0
1299 IF (ispmd /= 0) CALL spmd_chkw(iwiout,iout)
1300 CALL spmd_glob_isum9(iwiout,1)
1301 CALL spmd_ibcast(iwiout,iwiout,1,1,0,2)
1302 IF (iwiout > 0) THEN
1303 CALL spmd_wiout(iout,iwiout)
1304 iwiout = 0
1305 ENDIF
1306 ENDIF
1307 END IF
1308C-----------------------
1309C 4. ROD ELEMS
1310C-----------------------
1311 IF(neltof > 0)THEN
1312 IF(ispmd == 0.AND.mcheck == 0)WRITE (iout,1840)
1313 nn=(neltof+4)/5
1314 nbc=5
1315 DO i=1,nn
1316 READ (iin,'(10I10)') (nlec(ii),ii=1,10)
1317 DO j=1,nbc
1318 klg=nlec(2*j-1)
1319 kug=nlec(2*j )
1320 IF(klg<=0) GO TO 420
1321 nblk=(i-1)*nbc+j
1322 DO l=1,numelt
1323 IF(ixt(nixt,l)>=klg.AND.ixt(nixt,l)<=kug) THEN
1324 DO k=klg,kug
1325 IF(ixt(nixt,l) == k) THEN
1326 ixt(1,l)=-iabs(ixt(1,l))
1327 GOTO 411
1328 ENDIF
1329 ENDDO
1330 ENDIF
1331 411 CONTINUE
1332 ENDDO
1333 ENDDO
1334 ENDDO
1335 420 CONTINUE
1336 k=0
1337 DO j=1,numelt
1338 IF(ixt(1,j) < 0)THEN
1339 k=k+1
1340 nlec(k) = ixt(nixt,j)
1341 IF(k == 10)THEN
1342 IF(mcheck == 0)WRITE (iout,'(5I10)') (nlec(ii),ii=1,k)
1343 k=0
1344 ENDIF
1345 ENDIF
1346 ENDDO
1347 IF(k > 0) THEN
1348 IF(mcheck == 0) WRITE (iout,'(5I10)') (nlec(ii),ii=1,k)
1349 ENDIF
1350 IF(nspmd > 1) THEN
1351 ! required treatments to retrieve deleted elems in correct order
1352 iwiout = 0
1353 IF (ispmd /= 0) CALL spmd_chkw(iwiout,iout)
1354 CALL spmd_glob_isum9(iwiout,1)
1355 CALL spmd_ibcast(iwiout,iwiout,1,1,0,2)
1356 IF (iwiout > 0) THEN
1357 CALL spmd_wiout(iout,iwiout)
1358 iwiout = 0
1359 ENDIF
1360 ENDIF
1361 END IF
1362C-----------------------
1363C 5. BEAM ELEMS
1364C-----------------------
1365 IF(nelpof > 0)THEN
1366 IF(ispmd == 0.AND.mcheck == 0)WRITE (iout,1850)
1367 nn=(nelpof+4)/5
1368 nbc=5
1369 DO i=1,nn
1370 READ (iin,'(10I10)') (nlec(ii),ii=1,10)
1371 DO j=1,nbc
1372 klg=nlec(2*j-1)
1373 kug=nlec(2*j )
1374 IF(klg<=0) GO TO 520
1375 nblk=(i-1)*nbc+j
1376 DO l=1,numelp
1377 IF(ixp(nixp,l)>=klg.AND.ixp(nixp,l)<=kug) THEN
1378 DO k=klg,kug
1379 IF(ixp(nixp,l) == k) THEN
1380 ixp(1,l)=-iabs(ixp(1,l))
1381 GOTO 511
1382 ENDIF
1383 ENDDO
1384 ENDIF
1385 511 CONTINUE
1386 ENDDO
1387 ENDDO
1388 ENDDO
1389 520 CONTINUE
1390 k=0
1391 DO j=1,numelp
1392 IF(ixp(1,j) < 0)THEN
1393 k=k+1
1394 nlec(k) = ixp(nixp,j)
1395 IF(k == 10)THEN
1396 IF(mcheck == 0)WRITE (iout,'(5I10)') (nlec(ii),ii=1,k)
1397 k=0
1398 ENDIF
1399 ENDIF
1400 ENDDO
1401 IF(k > 0) THEN
1402 IF(mcheck == 0)WRITE (iout,'(5I10)') (nlec(ii),ii=1,k)
1403 ENDIF
1404 IF(nspmd > 1) THEN
1405 ! required treatments to retrieve deleted elems in correct order
1406 iwiout = 0
1407 IF (ispmd /= 0) CALL spmd_chkw(iwiout,iout)
1408 CALL spmd_glob_isum9(iwiout,1)
1409 CALL spmd_ibcast(iwiout,iwiout,1,1,0,2)
1410 IF (iwiout > 0) THEN
1411 CALL spmd_wiout(iout,iwiout)
1412 iwiout = 0
1413 ENDIF
1414 ENDIF
1415 END IF
1416C-----------------------
1417C 6. SPRING ELEMS
1418C-----------------------
1419 IF(nelrof > 0)THEN
1420 IF(ispmd == 0.AND.mcheck == 0)WRITE (iout,1860)
1421 nn=(nelrof+4)/5
1422 nbc=5
1423 DO i=1,nn
1424 READ (iin,'(10I10)') (nlec(ii),ii=1,10)
1425 DO j=1,nbc
1426 klg=nlec(2*j-1)
1427 kug=nlec(2*j )
1428 IF(klg<=0) GO TO 620
1429 nblk=(i-1)*nbc+j
1430 DO l=1,numelr
1431 IF(ixr(nixr,l)>=klg.AND.ixr(nixr,l)<=kug) THEN
1432 DO k=klg,kug
1433 IF(ixr(nixr,l) == k) THEN
1434 ixr(1,l)=-iabs(ixr(1,l))
1435 GOTO 611
1436 ENDIF
1437 ENDDO
1438 ENDIF
1439 611 CONTINUE
1440 ENDDO
1441 ENDDO
1442 ENDDO
1443 620 CONTINUE
1444 k=0
1445 DO j=1,numelr
1446 IF(ixr(1,j) < 0)THEN
1447 k=k+1
1448 nlec(k) = ixr(nixr,j)
1449 IF(k == 10)THEN
1450 IF(mcheck == 0)WRITE (iout,'(5I10)') (nlec(ii),ii=1,k)
1451 k=0
1452 ENDIF
1453 ENDIF
1454 ENDDO
1455 IF(k > 0) THEN
1456 IF(mcheck == 0) WRITE (iout,'(5I10)') (nlec(ii),ii=1,k)
1457 ENDIF
1458 IF(nspmd > 1) THEN
1459 ! required treatments to retrieve deleted elems in correct order
1460 iwiout = 0
1461 IF (ispmd /= 0) CALL spmd_chkw(iwiout,iout)
1462 CALL spmd_glob_isum9(iwiout,1)
1463 CALL spmd_ibcast(iwiout,iwiout,1,1,0,2)
1464 IF (iwiout > 0) THEN
1465 CALL spmd_wiout(iout,iwiout)
1466 iwiout = 0
1467 ENDIF
1468 ENDIF
1469 ENDIF
1470C-----------------------
1471C 7. 3N-SHELL-ELEMS
1472C-----------------------
1473 IF(neltgof > 0)THEN
1474 IF(ispmd == 0.AND.mcheck == 0)WRITE (iout,1870)
1475 nn=(neltgof+4)/5
1476 nbc=5
1477 DO i=1,nn
1478 READ (iin,'(10I10)') (nlec(ii),ii=1,10)
1479 DO j=1,nbc
1480 klg=nlec(2*j-1)
1481 kug=nlec(2*j )
1482 IF(klg<=0) GO TO 640
1483 nblk=(i-1)*nbc+j
1484 DO l=1,numeltg
1485 IF(ixtg(nixtg,l)>=klg.AND.ixtg(nixtg,l)<=kug) THEN
1486 DO k=klg,kug
1487 IF(ixtg(nixtg,l) == k) THEN
1488 ixtg(1,l)=-iabs(ixtg(1,l))
1489 GOTO 631
1490 ENDIF
1491 ENDDO
1492 ENDIF
1493 631 CONTINUE
1494 ENDDO
1495 ENDDO
1496 ENDDO
1497 640 CONTINUE
1498 k=0
1499 DO j=1,numeltg
1500 IF(ixtg(1,j) < 0)THEN
1501 k=k+1
1502 nlec(k) = ixtg(nixtg,j)
1503 IF(k == 10)THEN
1504 IF(mcheck == 0)WRITE (iout,'(5I10)') (nlec(ii),ii=1,k)
1505 k=0
1506 ENDIF
1507 ENDIF
1508 ENDDO
1509 IF(k > 0) THEN
1510 IF(mcheck == 0) WRITE (iout,'(5I10)') (nlec(ii),ii=1,k)
1511 ENDIF
1512 IF(nspmd > 1) THEN
1513 ! required treatments to retrieve deleted elems in correct order
1514 iwiout = 0
1515 IF (ispmd /= 0) CALL spmd_chkw(iwiout,iout)
1516 CALL spmd_glob_isum9(iwiout,1)
1517 CALL spmd_ibcast(iwiout,iwiout,1,1,0,2)
1518 IF (iwiout > 0) THEN
1519 CALL spmd_wiout(iout,iwiout)
1520 iwiout = 0
1521 ENDIF
1522 ENDIF
1523C
1524 END IF
1525C-----------------------
1526C 8. SPH PARTICLES
1527C-----------------------
1528 IF(nsphof > 0)THEN
1529 IF(ispmd == 0.AND.mcheck == 0)WRITE (iout,1880)
1530 nn=(nsphof+4)/5
1531 nbc=5
1532 k=0
1533 DO i=1,nn
1534 READ (iin,'(10I10)') (nlec(ii),ii=1,10)
1535 DO j=1,nbc
1536 klg=nlec(2*j-1)
1537 kug=nlec(2*j )
1538 IF(klg<=0) GO TO 732
1539 nblk=(i-1)*nbc+j
1540 DO l=1,numsph
1541 IF(kxsp(nisp,l)>=klg.AND.kxsp(nisp,l)<=kug) THEN
1542C number of the related group :
1543 ng =mod(abs(kxsp(2,l)),ngroup+1)
1544 IF(ng /= 0)THEN
1545 nft=iparg(3,ng)
1546 iad=iparg(4,ng)
1547 ii=l-nft
1548! ELBUF_STR(NG)%GBUF%OFF(II) = ZERO
1549 kxsp(2,l) = 0
1550 k=k+1
1551 nlecsph(k) = kxsp(nisp,l)
1552 IF(k == 10)THEN
1553 IF(mcheck == 0) WRITE (iout,'(5I10)') (nlecsph(ii),ii=1,k)
1554 k=0
1555 END IF
1556 END IF
1557 ENDIF
1558 ENDDO
1559 ENDDO
1560 ENDDO
1561 732 CONTINUE
1562 IF(k > 0) THEN
1563 IF(mcheck == 0)WRITE (iout,'(5I10)') (nlecsph(ii),ii=1,k)
1564 END IF
1565 IF(nspmd > 1) THEN
1566 ! required treatments to retrieve deleted elems in correct order
1567 iwiout = 0
1568 IF (ispmd /= 0) CALL spmd_chkw(iwiout,iout)
1569 CALL spmd_glob_isum9(iwiout,1)
1570 CALL spmd_ibcast(iwiout,iwiout,1,1,0,2)
1571 IF (iwiout > 0) THEN
1572 CALL spmd_wiout(iout,iwiout)
1573 iwiout = 0
1574 ENDIF
1575 ENDIF
1576C
1577 END IF
1578!-----------------------
1579! 9. /DEL/PART
1580!-----------------------
1581 IF(npartof > 0)THEN
1582 IF(ispmd == 0.AND.mcheck == 0)WRITE (iout,1890)
1583 nn=(npartof+4)/5
1584 nbc=5
1585 ALLOCATE(ipartof(npart))
1586 ipartof = 0
1587 DO i=1,nn
1588 READ (iin,'(10I10)') (nlec(ii),ii=1,10)
1589 DO j=1,nbc
1590 klg=nlec(2*j-1)
1591 kug=nlec(2*j )
1592 IF(klg<=0) GO TO 832
1593 nblk=(i-1)*nbc+j
1594 DO l = 1, npart
1595 DO k=klg,kug
1596 IF(ipart(4 + (l - 1) * lipart1) == k) THEN
1597 ipartof(l)=1
1598 cycle
1599 ENDIF
1600 ENDDO
1601 ENDDO
1602 ENDDO
1603 ENDDO
1604 832 CONTINUE
1605 k=0
1606 DO l=1,npart
1607 IF(ipartof(l)== 1)THEN
1608 k=k+1
1609 nlec(k) = ipart(4 + (l - 1) * lipart1)
1610 IF(k == 10)THEN
1611 IF(mcheck == 0)WRITE (iout,'(5I10)') (nlec(ii),ii=1,k)
1612 k=0
1613 ENDIF
1614 ENDIF
1615 ENDDO
1616 IF(k > 0) THEN
1617 IF(mcheck == 0) WRITE (iout,'(5I10)') (nlec(ii),ii=1,k)
1618 ENDIF
1619 IF(nspmd > 1) THEN
1620 ! required treatments to retrieve deleted elems in correct order
1621 iwiout = 0
1622 IF (ispmd /= 0) CALL spmd_chkw(iwiout,iout)
1623 CALL spmd_glob_isum9(iwiout,1)
1624 CALL spmd_ibcast(iwiout,iwiout,1,1,0,2)
1625 IF (iwiout > 0) THEN
1626 CALL spmd_wiout(iout,iwiout)
1627 iwiout = 0
1628 ENDIF
1629 ENDIF
1630!
1631 END IF
1632
1633
1634C-----------------------
1635C RESET OFF ARRAY
1636C-----------------------
1637 IF (nelof > 0)THEN
1638 DO ng=1,ngroup
1639 offg => elbuf_str(ng)%GBUF%OFF
1640 mlw=iparg(1,ng)
1641 ity=iparg(5,ng)
1642 nel=iparg(2,ng)
1643 nft=iparg(3,ng)
1644 iad=iparg(4,ng)
1645 igof=0
1646C-----------------------
1647C 1. SOLID ELEMS
1648C-----------------------
1649 IF(ity == 1.AND.nelsof /= 0)THEN
1650 igof=1
1651 DO i=1,nel
1652 ii=i+nft
1653 IF(ixs(1,ii) < 0)THEN
1654 ixs(1,ii)=-ixs(1,ii)
1655 IF (mlw /= 0) THEN
1656 offg(i)=zero
1657 ELSE ! loi0, no off
1658 CALL ancmsg(msgid=238,anmode=aninfo_blind,i1=ixs(nixs,ii),c1='BRICK',c2='BRICK')
1659 ENDIF
1660 ELSE
1661 igof=0
1662 ENDIF
1663 ENDDO
1664C-----------------------
1665C 2. 2D ELEMS
1666C-----------------------
1667 ELSEIF(ity == 2.AND.nelqof /= 0)THEN
1668 igof=1
1669 DO i=1,nel
1670 ii=i+nft
1671 IF(ixq(1,ii) < 0)THEN
1672 ixq(1,ii)=-ixq(1,ii)
1673 offg(i) = zero
1674 ELSE
1675 igof=0
1676 ENDIF
1677 ENDDO
1678C-----------------------
1679C 3. SHELL ELEMS
1680C-----------------------
1681 ELSEIF(ity == 3.AND.nelcof /= 0)THEN
1682 igof=1
1683 DO i=1,nel
1684 ii=i+nft
1685 IF(ixc(1,ii) < 0)THEN
1686 ixc(1,ii)=-ixc(1,ii)
1687 IF (mlw /= 0) THEN
1688 offg(i) = zero
1689 ELSE ! loi0, no off
1690 CALL ancmsg(msgid=238,anmode=aninfo_blind,i1=ixc(nixc,ii), c1='SHELL',c2='SHELL')
1691 ENDIF
1692 ELSE
1693 igof=0
1694 ENDIF
1695 ENDDO
1696C-----------------------
1697C 4. ROD ELEMS
1698C-----------------------
1699 ELSEIF(ity == 4.AND.neltof /= 0)THEN
1700C IGOF=1
1701C removal of the group incompatible with gap truss option
1702 igof=0
1703 DO i=1,nel
1704 ii=i+nft
1705 IF(ixt(1,ii) < 0)THEN
1706 ixt(1,ii)=-ixt(1,ii)
1707 offg(i) = zero
1708C ELSE
1709C IGOF=0
1710 ENDIF
1711 ENDDO
1712C-----------------------
1713C 5. BEAM ELEMS
1714C-----------------------
1715 ELSEIF(ity == 5.AND.nelpof /= 0)THEN
1716 igof=1
1717 DO i=1,nel
1718 ii=i+nft
1719 IF(ixp(1,ii) < 0)THEN
1720 ixp(1,ii)=-ixp(1,ii)
1721 offg(i) = zero
1722 ELSE
1723 igof=0
1724 ENDIF
1725 ENDDO
1726C-----------------------
1727C 6. ELEMENTS RESSORTS
1728C-----------------------
1729 ELSEIF(ity == 6.AND.nelrof /= 0)THEN
1730 igof=0
1731 DO i=1,nel
1732 ii=i+nft
1733 IF(ixr(1,ii) < 0)THEN
1734 ixr(1,ii)=-ixr(1,ii)
1735 offg(i) = zero
1736 ENDIF
1737 ENDDO
1738C-----------------------
1739C 7. 3N-SHELL-ELEMS
1740C-----------------------
1741 ELSEIF(ity == 7.AND.neltgof /= 0)THEN
1742 igof=1
1743 DO i=1,nel
1744 ii=i+nft
1745 IF(ixtg(1,ii) < 0)THEN
1746 ixtg(1,ii)=-ixtg(1,ii)
1747 IF(mlw /= 0) THEN
1748 offg(i) = zero
1749 ELSE ! loi0, no off
1750 CALL ancmsg(msgid=238,anmode=aninfo_blind,i1=ixtg(nixtg,ii), c1='SH_3N',c2='SH_3N')
1751 ENDIF
1752 ELSE
1753 igof=0
1754 ENDIF
1755 ENDDO
1756C-----------------------
1757C 8. SPH PARTICLES
1758C-----------------------
1759 ELSEIF(ity == 51.AND.nsphof /= 0)THEN
1760 igof=1
1761 DO i=1,nel
1762 ii=i+nft
1763 IF(kxsp(2,ii) == 0)THEN
1764 offg(i) = zero
1765 ELSE
1766 igof=0
1767 ENDIF
1768 ENDDO
1769 ENDIF
1770C----------------------------------------
1771C CHECK FOR GROUP CLEANING
1772C----------------------------------------
1773 iparg(8,ng)=max0(iparg(8,ng),igof)
1774 ENDDO ! next NG
1775 END IF
1776! /DEL/PART
1777 IF(npartof > 0)THEN
1778! per group
1779 k1=1+lipart1*(npart+nthpart)+2*9*(npart+nthpart)
1780 k2=k1+numels
1781 k3=k2+numelq
1782 k4=k3+numelc
1783 k5=k4+numelt
1784 k6=k5+numelp
1785 k7=k6+numelr
1786 k8=k7
1787 k9=k8+numeltg
1788 k10= k9 + numelx
1789 k11= k10+ numsph
1790 DO ng = 1, ngroup
1791 offg => elbuf_str(ng)%GBUF%OFF
1792 mlw=iparg(1,ng)
1793 ity = iparg(5, ng)
1794 nel = iparg(2, ng)
1795 nft = iparg(3, ng)
1796 IF (nel == 0) cycle
1797 igof=0
1798 SELECT CASE (ity)
1799 CASE(1)
1800! 1. SOLID ELEMS
1801 igof=1
1802 DO i=1,nel
1803 ii=i+nft-1
1804 partid = ipart(k1 + ii)
1805 IF (ipartof(partid)==1) THEN
1806 IF (mlw/=0) offg(i) = zero
1807 ELSE
1808 igof=0
1809 ENDIF
1810 ENDDO
1811 CASE(2)
1812! 2. 2D ELEMS
1813 igof=1
1814 DO i=1,nel
1815 ii=i+nft-1
1816 partid = ipart(k2 + ii)
1817 IF (ipartof(partid)==1) THEN
1818 offg(i) = zero
1819 ELSE
1820 igof=0
1821 ENDIF
1822 ENDDO
1823 CASE(3)
1824! 3. SHELL ELEMS
1825 igof=1
1826 DO i=1,nel
1827 ii=i+nft-1
1828 partid = ipart(k3 + ii)
1829 IF (ipartof(partid)==1) THEN
1830 IF (mlw/=0) offg(i) = zero
1831 ELSE
1832 igof=0
1833 ENDIF
1834 ENDDO
1835 CASE(4)
1836! 4. ROD ELEMS
1837 igof=0
1838 DO i=1,nel
1839 ii=i+nft-1
1840 partid = ipart(k4 + ii)
1841 IF (ipartof(partid)==1) offg(i) = zero
1842 ENDDO
1843 CASE(5)
1844! 5. BEAM ELEMS
1845 igof=1
1846 DO i=1,nel
1847 ii=i+nft-1
1848 partid = ipart(k5 + ii)
1849 IF (ipartof(partid)==1) THEN
1850 offg(i) = zero
1851 ELSE
1852 igof=0
1853 ENDIF
1854 ENDDO
1855 CASE(6)
1856! 6. ELEMENTS RESSORTS
1857 igof=0
1858 DO i=1,nel
1859 ii=i+nft-1
1860 partid = ipart(k6 + ii)
1861 IF (ipartof(partid)==1) offg(i) = zero
1862 ENDDO
1863 CASE(7)
1864! 7. 3N-SHELL-ELEMS
1865 igof=1
1866 DO i=1,nel
1867 ii=i+nft-1
1868 partid = ipart(k7 + ii)
1869 IF (ipartof(partid)==1) THEN
1870 IF (mlw/=0) offg(i) = zero
1871 ELSE
1872 igof=0
1873 ENDIF
1874 ENDDO
1875 CASE(51)
1876! 8. SPH PARTICLES
1877 igof=1
1878 DO i=1,nel
1879 ii=i+nft-1
1880 partid = ipart(k10 + ii)
1881 IF (ipartof(partid)==1) THEN
1882 offg(i) = zero
1883 kxsp(2,nft+i) = 0
1884 ELSE
1885 igof=0
1886 ENDIF
1887 ENDDO
1888 END SELECT
1889 iparg(8,ng)=max0(iparg(8,ng),igof)
1890 ENDDO
1891 IF (ALLOCATED(ipartof)) DEALLOCATE(ipartof)
1892 END IF
1893C-------------------------------------------
1894C ADDING RIGID LINK BETWEEN NODES
1895C-------------------------------------------
1896 IF(nrlink > 0 .AND. mcheck == 0) THEN
1897 k1 =1
1898 k2 =1
1899 DO k=1,nrlink
1900 READ (iin,'(I8,1X,3I1,1X,3I1,I10)')n,i1,i2,i3,ir1,ir2,ir3,isk
1901 READ (iin,'(10I10)') (llink(k2+i-1),i=1,n)
1902 ic=i3+2*i2+4*i1
1903 icr=ir3+2*ir2+4*ir1
1904 IF(iroddl == 0)icr=0
1905 IF(ispmd == 0.AND.mcheck == 0)WRITE(iout,2100) k,i1,i2,i3,ir1,ir2,ir3,isk,n
1906 IF(ispmd == 0.AND.mcheck == 0)WRITE(iout,'(10I10)') (llink(k2+i-1),i=1,n)
1907 IF ( isk /= 0) THEN
1908 isk1 = 0
1909 isk2 = 0
1910 DO ll=0,numskw
1911 IF(isk == iskwn(4,ll+1)) THEN
1912 isk1 = ll
1913 isk2 = 1
1914 ENDIF
1915 ENDDO
1916 IF ( isk2 == 0) THEN
1917 CALL ancmsg(msgid=125,anmode=aninfo)
1918 CALL arret(2)
1919 ENDIF
1920 isk = isk1
1921 ENDIF
1922 ! shared verification & defining specific data structure
1923 CALL fr_rlink1(llink(k2),itabm1,fr_rl(1,k),n)
1924 CALL rlink0(
1925 1 v ,vr ,ms ,in ,nnlink(k1),
1926 2 nnlink(k1+1),nnlink(k1+2),llink(k2),n ,ic ,
1927 3 icr ,nnlink(k1+3),isk ,skew(1,isk+1),iskwn ,
1928 4 fr_rl(1,k) ,weight )
1929 k1=k1+4
1930 k2=k2+n
1931 IF(ic+icr == 0)THEN
1932 IF(ispmd == 0)CALL ancmsg(msgid=126,anmode=aninfo)
1933 CALL arret(2)
1934 ENDIF
1935 ENDDO
1936 ENDIF
1937C-------------------------------------------
1938C ALE - LINK ON GRID VELOCITY
1939C format v5 ('/VEL/*')
1940C-------------------------------------------
1941C
1942C |ALE LINK 1 |2 |NALELK
1943C +----+----+----+----+----+----+----+----+----+...-+----+----+----+--...+----+----+--...
1944C LINALE |uID | M1 | M2 | N | IC | IM |id1 |id2 |id3 |... |idn | ... |
1945C +----+----+----+----+----+----+----+----+----+...-+----+----+----+--...+----+----+--...
1946C (1:SLINALE) 1 2 3 4 5 6 6+1 6+N
1947C +LLINAL |
1948C (N=1 if grnod_id is used)
1949C
1950C First subarrays are used to define ale links from starter (1:SLINALE). Engine links are in (SLINALE+1:SLINALE+LLINAL)
1951C Only Starter part is written in restart files.
1952C
1953C M1 : MAIN node 1
1954C M2 : MAIN node 2
1955C N : number of nodes
1956C IC : dir XYZ
1957C IM : formulation option -1,0,1
1958C id*: node_id list or grnod_id
1959C
1960C N<0 means that ALE LINK is tagged in FRALNK() as defined from a grnod_id.
1961
1962C
1963C +--RADIOSS2
1964C +--LECINP
1965C | +--FREFORM (direct access writing IUSC1, counting number of /VEL/ALE cards)
1966C | +--FRALNK(src/freef/fralnk.F) (LINALE size calculation:SLINALE, reading/checking GRNOD option)
1967C +--LECTUR (storing data in LINALE)
1968C | +--FR_RLALE(priv/spmd/fr_rlink1.F) (check and tag available nodes on current domain)
1969C +--RESOL
1970C +--ALEWDX
1971C +--ALELIN(priv/ale/alelin.F) (ale link subroutine, updated to be used with grnod_id)
1972C
1973C
1974C
1975
1976C SLINALE : starter size
1977C LLINAL : engine size
1978C NALELK : engine /VEL/ALE card (obsolete)
1979C NALELINK: engine /ALELINK/VEL
1980C after reading : SLINALE <- SLINALE+LLINAL
1981C NALELK <- NALELK+NALELINK
1982
1983 nalelk_starter = 0
1984 IF(mcheck==0)nalelk_starter = slinale/7 !SLINALE : starter linale size
1985 nalelk_removed=0
1986 k=6
1987
1988 !setting MAIN node user id to internal id
1989 DO j=1,nalelk_starter
1990 IF(irun == 1 .AND. mcheck == 0)
1991 . CALL fr_rlale(linale(k-4),linale(k-3),linale(k+1),itabm1,-1)
1992 k=k+1+6
1993 ENDDO
1994
1995 IF(nalelk /= 0.AND.mcheck == 0)THEN
1996 k=slinale+6 !LINALE(1:SLINALE):starter cards ;LINALE(SLINALE+1:SLINALE+LLINE):engine cards
1997 !K=5 !en attendant de copier LIALE starter dans (1:SLINALE)
1998 DO j=1,nalelk
1999 READ (iin,'(3I10,5X,3I1,I10)')m1,m2,n,i1,i2,i3,im
2000 !---------------------------------!
2001 ! ALE LINK DEFINED FROM NODES !
2002 !---------------------------------!
2003 IF(n>0)THEN
2004 READ (iin,'(10I10)')(linale(k+i),i=1,n)
2005 ic=i3+2*i2+4*i1
2006 IF(ic == 0)THEN
2007 ic=7
2008 i1=1
2009 i2=1
2010 i3=1
2011 ENDIF
2012 IF(ispmd == 0.AND.mcheck == 0)WRITE(iout,2200)m1,m2,i1,i2,i3,im
2013 linale(k-4)=m1
2014 linale(k-3)=m2
2015 linale(k-2)=n
2016 linale(k-1)=ic
2017 linale(k)=im
2018 IF(ispmd == 0.AND.mcheck == 0)THEN
2019 WRITE(iout,'(10I10)')(linale(k+i),i=1,n)
2020 WRITE(iout,2201)
2021 ENDIF
2022c Verification shared and tag in - nodes not present
2023 CALL fr_rlale(linale(k-4),linale(k-3),linale(k+1),itabm1,n)
2024 k=k+n+6
2025 !---------------------------------!
2026 ! ALE LINK DEFINED FROM GRNOD !
2027 !---------------------------------!
2028 ELSE !N<=0
2029 !M1 < 0
2030 !GR_ID < 0
2031 READ (iin,'(I10)', err=998,iostat=ierror)gr_id
2032
2033 0998 IF(ierror/=0)THEN
2034 WRITE(istdo,*) ' ** ERROR IN ALE LINK: INVALID GRNOD_ID'
2035 gr_id = 0
2036 ENDIF
2037 linale(k+1)=ngr2usr(gr_id,igrnod,ngrnod)
2038
2039 IF(linale(k+1)==0)CALL arret(2) !node group not found
2040
2041 IF(im==0.AND.igrnod(linale(k+1))%SORTED /= 1)THEN
2042 !option 0 needs to be defined with a /GRNOD/NODENS
2043 WRITE(istdo,*)' ** ERROR IN ALE LINK:'
2044 WRITE(istdo,*) ' UNSORTABLE NODE GROUP REQUIRED WITH OPTION 0 (/GRNOD/NODENS)'
2045 WRITE(iout ,*)' ** ERROR IN ALE LINK:'
2046 WRITE(iout,*) ' UNSORTABLE NODE GROUP REQUIRED WITH OPTION 0 (/GRNOD/NODENS)'
2047 CALL arret(2)
2048 ENDIF
2049
2050 ic=i3+2*i2+4*i1
2051 IF(ic == 0)THEN
2052 ic=7
2053 i1=1
2054 i2=1
2055 i3=1
2056 ENDIF
2057 IF(ispmd == 0.AND.mcheck == 0)THEN
2058 WRITE(iout,2200)m1,m2,i1,i2,i3,im !N=1 with GRNOD option
2059 WRITE(iout,2202)gr_id
2060 ENDIF
2061 linale(k-4)=m1
2062 linale(k-3)=m2
2063 linale(k-2)=n
2064 linale(k-1)=ic
2065 linale(k)=im
2066 !shared verification & untag non present nodes
2067 CALL fr_rlale(linale(k-4),linale(k-3),linale(k+1),itabm1,n)
2068 k=k+1+6
2069 ENDIF
2070 ENDDO !next J
2071 slinale = slinale + llinal
2072 ENDIF !IF(NALELK /= 0.AND.MCHECK == 0)THEN
2073
2074C-------------------------------------------
2075C ALE - LINK ON GRID VELOCITY
2076C format >= v12 ('/ALE/LINK/*')
2077C-------------------------------------------
2078 IF(nalelink /= 0.AND.mcheck == 0)THEN
2079 IF(nalelk==0)k=slinale+6
2080 DO j=1,nalelink
2081 READ (iin,'(3I10,5X,3I1,I10)')m1,m2,n,i1,i2,i3,im
2082 IF(m1>0.AND.m2>0)THEN
2083 !---------------------------------!
2084 ! ALE LINK DEFINED FROM NODES !
2085 !---------------------------------!
2086 IF(n>0)THEN
2087 READ (iin,'(10I10)')(linale(k+i),i=1,n)
2088 ic=i3+2*i2+4*i1
2089 IF(ic == 0)THEN
2090 ic=7
2091 i1=1
2092 i2=1
2093 i3=1
2094 ENDIF
2095 IF(ispmd == 0.AND.mcheck == 0)WRITE(iout,2200)m1,m2,i1,i2,i3,im
2096 linale(k-4)=m1
2097 linale(k-3)=m2
2098 linale(k-2)=n
2099 linale(k-1)=ic
2100 linale(k)=im
2101 IF(ispmd == 0.AND.mcheck == 0)THEN
2102 WRITE(iout,'(10I10)')(linale(k+i),i=1,n)
2103 WRITE(iout,2201)
2104 ENDIF
2105 ! Verification shared and tag in - nodes not present
2106 CALL fr_rlale(linale(k-4),linale(k-3),linale(k+1),itabm1,n)
2107 k=k+n+6
2108 !---------------------------------!
2109 ! ALE LINK DEFINED FROM GRNOD !
2110 !---------------------------------!
2111 ELSE !(N<=0)
2112 READ (iin,'(I10)', err=999,iostat=ierror)gr_id
2113 0999 IF(ierror/=0)THEN
2114 WRITE(istdo,*) ' ** ERROR IN ALE LINK: CANNOT READ GRNOD_ID VALUE'
2115 gr_id = 0
2116 CALL arret(2) !node group not read
2117 ENDIF
2118
2119 linale(k+1)=ngr2usr(gr_id,igrnod,ngrnod)
2120 IF(linale(k+1)==0)THEN
2121 WRITE(istdo,*) ' ** ERROR IN ALE LINK: INVALID GRNOD_ID'
2122 CALL arret(2) !node group not found
2123 ENDIF
2124
2125 IF(im==0.AND.igrnod(linale(k+1))%SORTED /= 1)THEN
2126 !option 0 needs to be defined with a /GRNOD/NODENS
2127 WRITE(istdo,*)' ** ERROR IN ALE LINK:'
2128 WRITE(istdo,*)' UNSORTABLE NODE GROUP REQUIRED WITH OPTION 0 (/GRNOD/NODENS)'
2129 WRITE(iout ,*)' ** ERROR IN ALE LINK:'
2130 WRITE(iout,*) ' UNSORTABLE NODE GROUP REQUIRED WITH OPTION 0 (/GRNOD/NODENS)'
2131 CALL arret(2)
2132 ENDIF
2133
2134 ic=i3+2*i2+4*i1
2135 IF(ic == 0)THEN
2136 ic=7
2137 i1=1
2138 i2=1
2139 i3=1
2140 ENDIF
2141 IF(ispmd == 0.AND.mcheck == 0)THEN
2142 WRITE(iout,2200)m1,m2,i1,i2,i3,im !N=1 with GRNOD option
2143 WRITE(iout,2202)gr_id
2144 ENDIF
2145 linale(k-4)=m1
2146 linale(k-3)=m2
2147 linale(k-2)=n
2148 linale(k-1)=ic
2149 linale(k)=im
2150 !shared verification & untag non present nodes
2151 CALL fr_rlale(linale(k-4),linale(k-3),linale(k+1),itabm1,n)
2152 k=k+1+6
2153 ENDIF !(N>0)
2154
2155 !---------------------------------!
2156 ! ALE LINK SET ON !
2157 !---------------------------------!
2158 ELSEIF(m1==-1)THEN ! ON
2159 ALLOCATE(alelin_on_off(iabs(n)))
2160 READ (iin,'(10I10)')(alelin_on_off(i),i=1,n)
2161 !L=1 !pos
2162 kk=0 ! number of ale link deactivated
2163 DO i1=1,n
2164 l=1 !pos
2165 m1=alelin_on_off(i1)
2166 DO i=1,nalelk_starter
2167 !---------------------------------!
2168 ! CHECK AND DEACTIVATE !
2169 !---------------------------------!
2170 IF(l>=slinale)THEN
2171 !warning : uID does not exist
2172 WRITE(istdo,*)' ** WARNING ALE LINK DOES NOT EXIST :',m1
2173 WRITE(iout,*) ' ** WARNING ALE LINK DOES NOT EXIST :',m1
2174 EXIT
2175 ELSE
2176 uid = linale(l+0)
2177 n = linale(l+3)
2178 IF(uid==-m1)THEN
2179 nalelk_removed = nalelk_removed+1 !counting deactivated links
2180 linale(l+0)=-linale(l+0) !setting negative uID to skip ALE LINK treatment (alelin.F)
2181 !printout :
2182 IF(ispmd == 0.AND.mcheck == 0)WRITE(iout,2211)m1 !N=1 with GRNOD option
2183 EXIT !next ALELIN_ON_OFF
2184 ELSEIF(uid==m1)THEN
2185 nalelk_removed = nalelk_removed+1 !counting deactivated links
2186 !warning uID already set OFF :
2187 WRITE(istdo,*)' ** WARNING ALE LINK ALREADY ACTIVATED : ',m1
2188 EXIT !next ALELIN_ON_OFF
2189 ENDIF
2190
2191 ENDIF
2192 l = l+6+iabs(n) !next uid position in linale()
2193 ENDDO
2194 ENDDO
2195 DEALLOCATE(alelin_on_off)
2196 !IF(NALELK_removed>0)NALELK_removed=1 !/ALE/LINK/ON or OFF card is no longer taken into account
2197
2198 !---------------------------------!
2199 ! ALE LINK SET OFF !
2200 !---------------------------------!
2201 ELSEIF(m1==-2)THEN ! OFF
2202 ALLOCATE(alelin_on_off(iabs(n)))
2203 READ (iin,'(10I10)')(alelin_on_off(i),i=1,n)
2204 !L=1 !pos
2205 kk=0 ! number of ale link deactivated
2206 DO i1=1,n
2207 l=1 !pos
2208 m1=alelin_on_off(i1)
2209 DO i=1,nalelk_starter
2210 !---------------------------------!
2211 ! CHECK AND DEACTIVATE !
2212 !---------------------------------!
2213 IF(l>=slinale)THEN
2214 !warning : uID does not exist
2215 WRITE(istdo,*) ' ** WARNING ALE LINK DOES NOT EXIST :',m1
2216 WRITE(iout,*) ' ** WARNING ALE LINK DOES NOT EXIST :',m1
2217 EXIT
2218 ELSE
2219 uid = linale(l+0)
2220 n = linale(l+3)
2221 IF(uid==m1)THEN
2222 nalelk_removed = nalelk_removed+1 !counting deactivated links
2223 linale(l+0)=-linale(l+0) !setting negative uID to skip ALE LINK treatment (alelin.F)
2224 !printout :
2225 IF(ispmd == 0.AND.mcheck == 0)WRITE(iout,2210)m1 !N=1 with GRNOD option
2226 EXIT !next ALELIN_ON_OFF
2227 ELSEIF(uid==-m1)THEN
2228 nalelk_removed = nalelk_removed+1 !counting deactivated links
2229 !warning uID already set OFF :
2230 WRITE(istdo,*)' ** WARNING ALE LINK ALREADY DEACTIVATED :',m1
2231 EXIT !next ALELIN_ON_OFF
2232 ENDIF
2233 ENDIF
2234 l = l+6+iabs(n) !next uID position in LINALE()
2235 ENDDO
2236 ENDDO
2237 DEALLOCATE(alelin_on_off)
2238 !IF(NALELK_removed>0)NALELK_removed=1 !/ALE/LINK/ON or OFF card is no longer taken into account
2239 ENDIF
2240
2241 ENDDO
2242 slinale = slinale + llinal
2243 nalelk = nalelk+nalelink-nalelk_removed
2244 ENDIF !(NALELINK /= 0.AND.MCHECK == 0)
2245
2246 IF(mcheck==0)nalelk=nalelk+nalelk_starter
2247
2248C-------------------------------------------
2249C ALE ON/ OFF
2250C-------------------------------------------
2251 IF(nbpartaleon > 0) ALLOCATE(aleon_partids_tmp(nbpartaleon))
2252 IF(nbpartaleoff > 0) ALLOCATE(aleoff_partids_tmp(nbpartaleoff))
2253 ion = 0
2254 ioff = 0
2255 DO i = 1, nbpartaleon + nbpartaleoff
2256 READ(iin, '(I10, I10)') partid, ison
2257 IF (ison == 1) THEN
2258 ion = ion + 1
2259 aleon_partids_tmp(ion) = partid
2260 ELSE IF(ison == 0) THEN
2261 ioff = ioff + 1
2262 aleoff_partids_tmp(ioff) = partid
2263 ENDIF
2264 ENDDO
2265C Checking for duplicates and removal
2266 DO i = 1, nbpartaleon
2267 partid = aleon_partids_tmp(i)
2268 IF (partid /= -1) THEN
2269 DO j = 1, nbpartaleon
2270 IF (j == i) cycle
2271 IF (aleon_partids_tmp(j) == partid) THEN
2272 nbpartaleon = nbpartaleon - 1
2273 aleon_partids_tmp(j) = -1
2274 ENDIF
2275 ENDDO
2276 ENDIF
2277 ENDDO
2278 DO i = 1, nbpartaleoff
2279 partid = aleoff_partids_tmp(i)
2280 IF (partid /= -1) THEN
2281 DO j = 1, nbpartaleoff
2282 IF (j == i) cycle
2283 IF (aleoff_partids_tmp(j) == partid) THEN
2284 nbpartaleoff = nbpartaleoff - 1
2285 aleoff_partids_tmp(j) = -1
2286 ENDIF
2287 ENDDO
2288C Check if part is also defined in ON
2289C By default, all parts are ON, if a part is set to OFF
2290C and ON at the same time, we choose to OFF it
2291 DO j = 1, nbpartaleon
2292 IF (aleon_partids_tmp(j) == partid) THEN
2293 !!! WARNING TO OUTPUT HERE
2294 aleon_partids_tmp(j) = -1
2295 nbpartaleon = nbpartaleon - 1
2296 CALL ancmsg(msgid = 272, anmode = aninfo, i1 = partid)
2297 ENDIF
2298 ENDDO
2299 ENDIF
2300 ENDDO
2301C Checking that part ids exist in IPART(4, *)
2302 nbaleoff_part = nbpartaleoff
2303 nbaleon_part = nbpartaleon
2304 DO i = 1, nbpartaleon
2305 partid = aleon_partids_tmp(i)
2306 ison = 0
2307 IF (partid /= -1) THEN
2308 DO j = 1, npart
2309 IF (ipart(4 + (j - 1) * lipart1) == partid) THEN
2310 ison = 1
2311 aleon_partids_tmp(i) = j
2312 EXIT
2313 ENDIF
2314 ENDDO
2315 ENDIF
2316 IF (ison == 0) THEN
2317 aleon_partids_tmp(i) = -1
2318 nbaleon_part = nbaleon_part - 1
2319 CALL ancmsg(msgid = 271, anmode = aninfo, i1 = partid)
2320 CALL arret(2)
2321 ENDIF
2322 ENDDO
2323 DO i = 1, nbpartaleoff
2324 partid = aleoff_partids_tmp(i)
2325 ison = 0
2326 IF (partid /= -1) THEN
2327 DO j = 1, npart
2328 IF (ipart(4 + (j - 1) * lipart1) == partid) THEN
2329 ison = 1
2330 aleoff_partids_tmp(i) = j
2331 EXIT
2332 ENDIF
2333 ENDDO
2334 ENDIF
2335 IF (ison == 0) THEN
2336 aleoff_partids_tmp(i) = -1
2337 nbaleoff_part = nbaleoff_part - 1
2338 CALL ancmsg(msgid = 271, anmode = aninfo, i1 = partid)
2339 CALL arret(2)
2340 ENDIF
2341 ENDDO
2342C Filling module values
2343 IF (nbaleoff_part > 0) ALLOCATE(aleoff_partids(nbaleoff_part))
2344 IF (nbaleon_part > 0) ALLOCATE(aleon_partids(nbaleon_part))
2345 DO i = 1, nbpartaleon
2346 partid = aleon_partids_tmp(i)
2347 ison = 0
2348 IF (partid /= -1) THEN
2349 ison = ison + 1
2350 aleon_partids(ison) = partid
2351 ENDIF
2352 ENDDO
2353 ison = 0
2354 DO i = 1, nbpartaleoff
2355 partid = aleoff_partids_tmp(i)
2356 IF (partid /= -1) THEN
2357 ison = ison + 1
2358 aleoff_partids(ison) = partid
2359 ENDIF
2360 ENDDO
2361C OUTPUT
2362 DO i = 1, nbaleoff_part
2363 WRITE(iout, 4800) aleoff_partids(i)
2364 ENDDO
2365C-------------------------------------------
2366C Fill IPARG
2367C-------------------------------------------
2368 k1 = 1 + lipart1 * (npart + nthpart) + 2 * 9 * (npart + nthpart)
2369 k2 = k1 + numels
2370 DO ng = 1, ngroup
2371 ity = iparg(5, ng)
2372 nel = iparg(2, ng)
2373 nft = iparg(3, ng)
2374C Initially, IPARG(76, *) = 0 --> ON
2375C When restarting a computation, IPARG(76, NG) = Whatever was put here during previous run
2376C Find part associated to the current group
2377C by checking the first element of the group
2378C NB : all elements in a same group belong necessarily to the same part
2379
2380 IF (ity == 1) THEN
2381C Solid groups
2382 partid = ipart(k1 + nft)
2383 ELSEIF (ity == 2) THEN
2384C Shells and coques
2385 partid = ipart(k2 + nft)
2386 ENDIF
2387C Check if this part is found in ALEOFF_PARTIDS
2388 DO i = 1, nbaleoff_part
2389 IF (aleoff_partids(i) == partid) THEN
2390 iparg(76, ng) = 1 ! --> OFF
2391 EXIT
2392 ENDIF
2393 ENDDO
2394C Check if this part is found in ALEON_PARTIDS
2395 DO i = 1, nbaleon_part
2396 IF (aleon_partids(i) == partid) THEN
2397 iparg(76, ng) = 0 ! --> ON
2398 EXIT
2399 ENDIF
2400 ENDDO
2401 ENDDO
2402 IF (ALLOCATED(aleon_partids_tmp)) DEALLOCATE(aleon_partids_tmp)
2403 IF (ALLOCATED(aleoff_partids_tmp)) DEALLOCATE(aleoff_partids_tmp)
2404 IF (ALLOCATED(aleon_partids)) DEALLOCATE(aleon_partids)
2405 IF (ALLOCATED(aleoff_partids)) DEALLOCATE(aleoff_partids)
2406C-------------------------------------------
2407C MODIFICATION DES CONDITIONS AUX LIMITES
2408C-------------------------------------------
2409 IF(nubcsn /= 0)THEN
2410 CALL lcbcsf(icode,iskew,nubcsn,itab,itabm1,npby ,iskwn,weight)
2411 IF(nspmd > 1) THEN
2412C necessary processing to recover limit cond nodes
2413 iwiout = 0
2414 IF (ispmd /= 0) CALL spmd_chkw(iwiout,iout)
2415 CALL spmd_glob_isum9(iwiout,1)
2416 CALL spmd_ibcast(iwiout,iwiout,1,1,0,2)
2417 IF (iwiout > 0) THEN
2418 CALL spmd_wiout(iout,iwiout)
2419 iwiout = 0
2420 END IF
2421 END IF
2422 END IF
2423C-------------------------------------------
2424C INTERFACE MODIFICATION
2425C-------------------------------------------
2426 IF(nintch /= 0)THEN
2427 DO i=1,nintch
2428 READ(iin,'(2I8,2F16.0)')noint,nsearch,tstart,tfin
2429 IF(tfin == 0.0)tfin=1.e30
2430 IF(tfin <= tstart) THEN
2431 CALL ancmsg(msgid=307,anmode=anstop)
2432 CALL arret(2)
2433 ENDIF
2434
2435 jpri = -1
2436 DO j=1,ninter
2437 IF(noint == ipari(15,j)) THEN
2438 IF(nsearch /= 0) ipari(13,j)=nsearch
2439 intbuf_tab(j)%VARIABLES(3) = tstart
2440 intbuf_tab(j)%VARIABLES(11) = tfin
2441 jpri=j
2442 ENDIF
2443 ENDDO
2444 IF(jpri/=-1) THEN
2445 IF(ispmd == 0.AND.mcheck == 0) WRITE(iout,2300)noint,ipari(13,jpri),tstart,tfin
2446 ELSE
2447 IF(ispmd == 0.AND.mcheck == 0) WRITE(iout,2301)noint
2448 ENDIF
2449 ENDDO
2450 ENDIF
2451C
2452 IF(mcheck /= 0)THEN
2453 DO i = 1,mx_outp
2454 nv_outp = nv_outp + outp_v(i)
2455 nss_outp = nss_outp + outp_ss(i)
2456 nst_outp = nst_outp + outp_st(i)
2457 ncs_outp = ncs_outp + outp_cs(i)
2458 nct_outp = nct_outp + outp_ct(i)
2459 nts_outp = nts_outp + outp_ts(i)
2460 nps_outp = nps_outp + outp_ps(i)
2461 npt_outp = npt_outp + outp_pt(i)
2462 nrs_outp = nrs_outp + outp_rs(i)
2463 nrt_outp = nrt_outp + outp_rt(i)
2464 ENDDO
2465 IF(outp_v(12) == 1)nv_outp=nv_outp+1
2466 ENDIF
2467C-------------------------------------------
2468C MODIFICATION OF RIGID BODY
2469C-------------------------------------------
2470 IF(nrbyof /= 0)
2471 1 CALL rbyonf(iparg,ipari ,ms ,in ,
2472 2 ixs ,ixq ,ixc ,ixt ,ixp ,
2473 3 ixr ,skew ,itab ,itabm1 ,iskwn ,
2474 4 npby ,0 ,nrbyof,wa ,lpby ,
2475 5 rby ,x ,v ,vr ,ixtg ,
2476 6 igrv ,ibgr ,weight,fr_rby2,partsav,
2477 7 ipart,elbuf_str,icfield,lcfield,tagslv_rby)
2478C
2479 IF(nrbyon /= 0)
2480 1 CALL rbyonf(iparg,ipari ,ms ,in ,
2481 2 ixs ,ixq ,ixc ,ixt ,ixp ,
2482 3 ixr ,skew ,itab ,itabm1 ,iskwn ,
2483 4 npby ,1 ,nrbyon,wa ,lpby ,
2484 5 rby ,x ,v ,vr ,ixtg ,
2485 6 igrv ,ibgr ,weight,fr_rby2,partsav,
2486 7 ipart,elbuf_str,icfield,lcfield,tagslv_rby)
2487C---------------------------------------------------------
2488C READING DATA FOR FLUX & ROTATION
2489C---------------------------------------------------------
2490 IF (nsflsw /= 0.AND.mcheck == 0) THEN
2491 CALL lecflsw (nsflsw,ntflsw,neflsw,nnflsw,crflsw,x,ixs,iparg,wa)
2492 ENDIF
2493C-------------------------------------------------
2494C MODIFICATION OF TIME FUNCTIONS
2495C-------------------------------------------------
2496 IF (nfct /= 0) THEN
2497 CALL lecfun (npc, pld, nfct, npts, table)
2498 END IF
2499C-------------------------------------------------
2500C RESET VELOCITIES
2501C-------------------------------------------------
2502 IF (niniv /= 0) THEN
2503 CALL lecinv (niniv,x,v,vr,itab,iframe,xframe,igrnod,fxbipm,fxbvit,fxbrpm)
2504 END IF
2505C---------------------------------------------------------
2506C CUTS
2507C---------------------------------------------------------
2508 IF(ncuts > 0.AND.mcheck == 0)CALL leccut(icut,xcut,itabm1)
2509C---------------------------------------------------------
2510C ANIMATION (DT ,DMAS)
2511C---------------------------------------------------------
2512! DO I=1,SANIN
2513! ANIN(I)=ZERO
2514! ENDDO
2515C---------------------------------------------------------
2516C FILTERED SAMPLED OUTPUT
2517C---------------------------------------------------------
2518 IF(nnoise > 0)CALL lecnoise(inoise,itabm1,names_and_titles)
2519 nnoiser=nnoise
2520C---------------------------------------------------------
2521C MADYMO COUPLING GENERAL DATA (CONVERSION FACTORS).
2522C---------------------------------------------------------
2523#ifdef DNC
2524 IF( imadcpl > 0.AND.invers>=40.AND.mcheck == 0)CALL leccpl()
2525#endif
2526C---------------------------------------------------------
2527C OPTIONS SPMD SPECIFIQUES
2528 IF(ipread > 0.AND.ispmd == 0.AND.mcheck == 0) WRITE(iout,4000)
2529 IF(iddw > 0.AND.ispmd == 0.AND.mcheck == 0) WRITE(iout,4100)
2530C---------------------------------------------------------
2531C RAYLEIGH DAMPING
2532C---------------------------------------------------------
2533 CALL lecdamp(ndampn, dampr, igrnod)
2534C---------------------------------------------------------
2535 istatcnd_sav= istatcnd
2536 IF(impl_s/=0) THEN
2537 IF(isprb==0.AND.nadmesh/=0)THEN
2538 IF (ispmd == 0) THEN
2539 CALL ancmsg(msgid=131,anmode=aninfo)
2540 ENDIF
2541 CALL arret(2)
2542 ENDIF
2543 istatcnd = 0
2544 END IF
2545C---------------------------------------------------------
2546 IF (impl_s == 1.OR.neig > 0) CALL lecimpl
2547 IF (nfxinp > 0) CALL lecfxinp(nfxinp)
2548 IF(ale%SUB%NODSUBDT /= 0)THEN
2549 IF(ispmd == 0) THEN
2550 CALL ancmsg(msgid=129,anmode=aninfo_blind)
2551 END IF
2552 CALL arret(2)
2553 END IF
2554 IF (neigoff > 0) CALL leceig(neigoff, neoff)
2555 IF (nfvmesh > 0) CALL lecfvbag(nfvmesh, monvol, volmon, x)
2556 IF (nfvmodi > 0) CALL lecfvbag1(nfvmodi, monvol, volmon)
2557C-------------------------------------------
2558C .sta files
2559C-------------------------------------------
2560 IF(mcheck /= 0)THEN
2561 nc_stat = 0
2562 DO i = 1,mx_stat
2563 nc_stat = nc_stat + stat_c(i)
2564 ENDDO
2565 ENDIF
2566 IF(nstatprt /= 0 .OR. nstatall /= 0)THEN
2567 k1=1+lipart1*(npart+nthpart)+2*9*(npart+nthpart)
2568 k2=k1+numels
2569 k3=k2+numelq
2570 k4=k3+numelc
2571 k5=k4+numelt
2572 k6=k5+numelp
2573 k7=k6+numelr
2574 k8=k7
2575 k9=k8+numeltg
2576 CALL lecstat(ipart,ipart_state,elbuf_str,ipm,iparg,ipart(k1),
2577 . ipart(k3),ipart(k8),mat_param)
2578 END IF
2579C-------------------------------------------
2580C .dynain files
2581C-------------------------------------------
2582 IF(dynain_data%NDYNAINPRT /= 0 .OR. dynain_data%NDYNAINALL /= 0)THEN
2583 k1=1+lipart1*(npart+nthpart)+2*9*(npart+nthpart)
2584 k2=k1+numels
2585 k3=k2+numelq
2586 k4=k3+numelc
2587 k5=k4+numelt
2588 k6=k5+numelp
2589 k7=k6+numelr
2590 k8=k7
2591 k9=k8+numeltg
2592 CALL read_dynain(ipart,dynain_data,ipart(k3),ipart(k8),ixc,ixtg)
2593 END IF
2594C-------------------------------------------
2595C .h3d files
2596C-------------------------------------------
2597 CALL prelech3d(numgeo ,npropgi ,npropmi ,nummat ,numply ,
2598 . igeo ,ipm ,h3d_data ,multi_fvm,mds_output_table,
2599 . mds_nmat ,max_depvar,mds_ndepsvar,mat_param,numsphg)
2600 CALL lech3d(output,geo,igeo,ipm,ipart,h3d_data,multi_fvm,ipari,iparg,tag_skins6,
2601 . mds_label,mds_output_table,mds_nmat,max_depvar,mds_ndepsvar,
2602 . elbuf_str,stack,ibcl,iloadp,lloadp,loads,mat_param,pblast,
2603 . igrpart,npc,pld,snpc,stf)
2604c
2605 IF(h3d_data%N_OUTP_H3D /= 0 .AND. ispmd == 0)THEN
2606 WRITE(iout,5000)h3d_data%TH3D,h3d_data%DTH3D
2607 WRITE(iout,*)' |'
2608 DO i=1,h3d_data%N_OUTP_H3D
2609 IF( h3d_data%OUTPUT_LIST(i)%ETYPE == 1 .OR.
2610 . h3d_data%OUTPUT_LIST(i)%ETYPE == 2 .AND. numelcg+numeltgg > 0 .OR.
2611 . h3d_data%OUTPUT_LIST(i)%ETYPE == 3 .AND. numelsg > 0 .OR.
2612 . h3d_data%OUTPUT_LIST(i)%ETYPE == 4 .AND. numeltrg+numelpg+numelrg > 0 .OR.
2613 . h3d_data%OUTPUT_LIST(i)%ETYPE == 5 .AND. numsphg > 0 .OR.
2614 . h3d_data%OUTPUT_LIST(i)%ETYPE == 6 .AND. numelqg > 0 .OR.
2615 . h3d_data%OUTPUT_LIST(i)%ETYPE == 7 .AND. numsking > 0) THEN
2616 char1=' '
2617 IF(h3d_data%OUTPUT_LIST(i)%ETYPE == 1)THEN
2618 char1='NODAL'
2619 ELSEIF(h3d_data%OUTPUT_LIST(i)%ETYPE == 2)THEN
2620 char1='SHELL'
2621 ELSEIF(h3d_data%OUTPUT_LIST(i)%ETYPE == 3)THEN
2622 char1='SOLID'
2623 ELSEIF(h3d_data%OUTPUT_LIST(i)%ETYPE == 4)THEN
2624 char1='ONED'
2625 ELSEIF(h3d_data%OUTPUT_LIST(i)%ETYPE == 5)THEN
2626 char1='SPH'
2627 ELSEIF(h3d_data%OUTPUT_LIST(i)%ETYPE == 6)THEN
2628 char1='QUAD'
2629 ELSEIF(h3d_data%OUTPUT_LIST(i)%ETYPE == 7)THEN
2630 char1='SKIN'
2631 ENDIF
2632c
2633 char2=' '
2634 IF(h3d_data%OUTPUT_LIST(i)%OUTP_TYPE == 1)THEN
2635 char2='SCALAR'
2636 ELSEIF(h3d_data%OUTPUT_LIST(i)%OUTP_TYPE == 2)THEN
2637 char2='VECTOR'
2638 ELSEIF(h3d_data%OUTPUT_LIST(i)%OUTP_TYPE == 3)THEN
2639 char2='TENSOR'
2640 ELSEIF(h3d_data%OUTPUT_LIST(i)%OUTP_TYPE == 4)THEN
2641 char2='TORSOR'
2642 ENDIF
2643c
2644 stext1 = h3d_data%OUTPUT_LIST(i)%S_STRING1
2645 WRITE(iout,*) ' |----'//char1//' '//char2//' : '
2646c WRITE(IOUT,*) ' |----'//H3D_DATA%OUTPUT_LIST(I)%KEYWORD
2647 WRITE(iout,*) ' | '//h3d_data%OUTPUT_LIST(i)%STRING1(1:stext1)
2648
2649 IF (h3d_data%OUTPUT_LIST(i)%IUVAR > 0)
2650 . WRITE(iout,*) ' | UVAR=', h3d_data%OUTPUT_LIST(i)%IUVAR
2651
2652 IF (h3d_data%OUTPUT_LIST(i)%PLY > 0)
2653 . WRITE(iout,*) ' | PLY=',h3d_data%OUTPUT_LIST(i)%PLY
2654
2655 IF (h3d_data%OUTPUT_LIST(i)%LAYER > 0)
2656 . WRITE(iout,*) ' | LAYER=',h3d_data%OUTPUT_LIST(i)%LAYER
2657
2658 IF (h3d_data%OUTPUT_LIST(i)%IPT > 0)
2659 . WRITE(iout,*) ' | IPT=',h3d_data%OUTPUT_LIST(i)%IPT
2660
2661 IF (h3d_data%OUTPUT_LIST(i)%IR > 0)
2662 . WRITE(iout,*) ' | IR=',h3d_data%OUTPUT_LIST(i)%IR
2663
2664 IF (h3d_data%OUTPUT_LIST(i)%IR > 0)
2665 . WRITE(iout,*) ' | IS=',h3d_data%OUTPUT_LIST(i)%IS
2666
2667 IF (h3d_data%OUTPUT_LIST(i)%IT > 0)
2668 . WRITE(iout,*) ' | IT=',h3d_data%OUTPUT_LIST(i)%IT
2669
2670 WRITE(iout,*)' |'
2671 ENDIF
2672 ENDDO
2673 WRITE(iout,*)' '
2674 WRITE(iout,*)' '
2675 ENDIF
2676C-------------------------------------------
2677C H3D/SENSOR
2678C-------------------------------------------
2679 IF (ispmd == 0.AND.mcheck == 0) THEN
2680 DO k=1,h3d_data%N_SENS_H3D
2681 iok = 0
2682 IF(h3d_data%LSENS_H3D(k) /= 0)THEN
2683 DO i=1,sensors%NSENSOR
2684 IF(h3d_data%LSENS_H3D(k) == sensor_tab(i)%SENS_ID)THEN
2685 h3d_data%LSENS_H3D(k)=i
2686 iok = 1
2687 EXIT
2688 ENDIF
2689 ENDDO
2690 ENDIF
2691 IF(iok == 0) THEN
2692 CALL ancmsg(msgid=283,anmode=aninfo,i1=h3d_data%LSENS_H3D(k))
2693 CALL arret(2)
2694 ENDIF
2695 ENDDO
2696 ENDIF !(ISPMD==0.AND.MCHECK==0)
2697C-------------------------------
2698C Output by interface :
2699C 1 - Skid lines
2700C-------------------------------
2701 IF(h3d_data%N_SCAL_SKID > 0.AND.nintskidold==0) THEN
2702 ninterskid = h3d_data%N_SCAL_SKID
2703 IF(nintstamp/=0) THEN
2704 ALLOCATE (pskids(ninterskid,numnodg))
2705 pskids(1:ninterskid,1:numnodg) = zero
2706 ELSE
2707 ALLOCATE (pskids(ninterskid,numnod))
2708 pskids(1:ninterskid,1:numnod) = zero
2709 ENDIF
2710 ELSEIF(h3d_data%N_SCAL_SKID == 0.AND.nintskidold==0)THEN
2711 ALLOCATE (pskids(0,0))
2712 ENDIF
2713C-------------------------------------------
2714C 2 - Frictionnal energy: allocs tabs
2715C--------------------------------------------
2716 IF(h3d_data%N_SCAL_CSE_FRIC >0.AND.output%DATA%S_EFRIC == 0)THEN
2717 ALLOCATE(output%DATA%EFRICG(numnod))
2718 output%DATA%EFRICG(1:numnod) = zero
2719 IF(nintstamp/=0) THEN
2720 ALLOCATE(output%DATA%EFRICG_STAMP(numnodg))
2721 output%DATA%EFRICG_STAMP(1:numnodg) = zero
2722 ELSE
2723 ALLOCATE(output%DATA%EFRICG_STAMP(0))
2724 ENDIF
2725 ELSEIF(output%DATA%S_EFRIC==0)THEN
2726 ALLOCATE(output%DATA%EFRICG(0))
2727 ALLOCATE(output%DATA%EFRICG_STAMP(0))
2728 ENDIF
2729
2730 IF(h3d_data%N_SCAL_CSE_FRICINT > 0.AND.output%DATA%S_EFRICINT==0) THEN
2731C
2732 ninefricg = h3d_data%N_SCAL_CSE_FRICINT
2733 output%DATA%NINEFRIC_STAMP = 0
2734 output%DATA%NINEFRIC = 0
2735 IF(nintstamp/=0) THEN
2736 DO n=1,ninter
2737 ni = h3d_data%N_CSE_FRIC_INTER (n)
2738 IF(ni/= 0.AND.ipari(7,n)==21) output%DATA%NINEFRIC_STAMP = output%DATA%NINEFRIC_STAMP + 1
2739 ENDDO
2740 ENDIF
2741 IF(output%DATA%NINEFRIC_STAMP==ninefricg) THEN
2742 CALL my_alloc(output%DATA%EFRIC_STAMP,ninefricg,numnodg)
2743 output%DATA%EFRIC_STAMP(1:ninefricg,1:numnodg) = zero
2744 CALL my_alloc(output%DATA%EFRIC,0,0)
2745 ELSEIF(output%DATA%NINEFRIC_STAMP==0) THEN
2746 output%DATA%NINEFRIC = ninefricg
2747 CALL my_alloc(output%DATA%EFRIC,output%DATA%NINEFRIC,numnod)
2748 output%DATA%EFRIC(1:ninefricg,1:numnod) = zero
2749 CALL my_alloc(output%DATA%EFRIC_STAMP,0,0)
2750 ELSE
2751 output%DATA%NINEFRIC = ninefricg-output%DATA%NINEFRIC_STAMP
2752 CALL my_alloc(output%DATA%EFRIC_STAMP,output%DATA%NINEFRIC_STAMP,numnodg)
2753 output%DATA%EFRIC_STAMP(1:output%DATA%NINEFRIC_STAMP,1:numnodg) = zero
2754 CALL my_alloc(output%DATA%EFRIC,output%DATA%NINEFRIC,numnod)
2755 output%DATA%EFRIC(1:output%DATA%NINEFRIC,1:numnod) = zero
2756 ns = 0
2757 nn= 0
2758 DO n=1,ninter
2759 ni = h3d_data%N_CSE_FRIC_INTER (n)
2760 IF(ni/= 0.AND.ipari(7,n)==21) THEN
2761 ns = ns+1
2762 h3d_data%N_CSE_FRIC_INTER (n) = output%DATA%NINEFRIC + ns
2763 ELSEIF (ni/=0) THEN
2764 nn = nn+1
2765 h3d_data%N_CSE_FRIC_INTER (n) = nn
2766 ENDIF
2767 ENDDO
2768 ENDIF
2769 ELSEIF(output%DATA%S_EFRICINT==0)THEN
2770 CALL my_alloc(output%DATA%EFRIC,0,0)
2771 CALL my_alloc(output%DATA%EFRIC_STAMP,0,0)
2772 ENDIF
2773C-----------------------------------
2774C-------------------------------------------
2775C Animations
2776C-------------------------------------------
2777
2778C------------------------------------------------
2779C Check function load of output frequencie and initial time step
2780C------------------------------------------------
2781 IF(output%nb_anim_frame >0.AND.output%DTANIM_FCT_ID==0) THEN
2782 IF(output%TANIM_STOP0 >zero.AND.output%TANIM_STOP0 /= ep20.AND.output%TANIM0 > 0) THEN
2783 output%DTANIM0 = (output%TANIM_STOP0 - output%TANIM0) / output%nb_anim_frame
2784 ELSEIF(output%TANIM_STOP0 >zero.AND.output%TANIM_STOP0 /= ep20) THEN
2785 output%DTANIM0 = (output%TANIM_STOP0-tt) / output%nb_anim_frame
2786 ELSEIF(output%TANIM0 >zero) THEN
2787 output%DTANIM0 = (tstop-output%TANIM0) / output%nb_anim_frame
2788 ELSE
2789 output%DTANIM0 = (tstop-tt) / output%nb_anim_frame
2790 ENDIF
2791 IF (output%DTANIM0 <= zero) THEN
2792 CALL ancmsg(msgid=293,anmode=aninfo,c1='ANIM',c2='ANIM')
2793 CALL arret(0)
2794 ENDIF
2795 ENDIF
2796
2797 ok = 0
2798 dtfct = 0
2799 IF (output%DTANIM_FCT_ID > 0) THEN
2800
2801 DO jj=1,nfunct
2802 IF(output%DTANIM_FCT_ID == npc(nfunct+2+jj-1)) THEN
2803 dtfct=jj
2804 ok = 1
2805 EXIT
2806 ENDIF
2807 ENDDO
2808 IF (ok == 0) THEN
2809 CALL ancmsg(msgid=315,
2810 . msgtype=msgerror,
2811 . anmode=aninfo_blind_1,
2812 . c1='ANIM',
2813 . i2=output%DTANIM_FCT_ID)
2814 ENDIF
2815 IF (ok >0) output%DTANIM_FCT_ID = dtfct
2816
2817 IF(output%DTANIM_FCT_ID > 0) THEN
2818 output%DTANIM0 = finter(output%DTANIM_FCT_ID,zero,npc,pld,dydx)
2819 !need to add the python also
2820 ENDIF
2821
2822 IF (output%DTANIM0 <= zero) THEN
2823 CALL ancmsg(msgid=293,anmode=aninfo,c1='ANIM',c2='ANIM')
2824 CALL arret(0)
2825 ENDIF
2826 ENDIF
2827
2828 IF(output%TANIM0 /= zero) output%TANIM = output%TANIM0
2829 IF(output%TANIM_STOP0/=zero .AND. output%TANIM_STOP0/=ep20) output%TANIM_STOP = output%TANIM_STOP0
2830 IF(output%DTANIM0 > zero) output%DTANIM= output%DTANIM0
2831 IF(output%DTANIM<=zero) output%TANIM = ep30
2832 IF (output%TANIM < tt-dt2.AND.output%DTANIM > zero)output%TANIM = output%TANIM
2833 . + int((tt-dt2-output%TANIM)/output%DTANIM)*output%DTANIM
2834 IF (output%TANIM < tt-dt2)output%TANIM = output%TANIM+output%DTANIM
2835C-----------------------
2836 IF(anim_v(14)+h3d_data%N_VECT_DROT > 0 .AND.
2837 . ((isecut == 0 .AND. iisrot == 0 .AND. impose_dr == 0 .AND. idrot == 0) .OR. iroddl == 0)) THEN
2838 IF(ispmd == 0) THEN
2839 WRITE(iout,*) ' ** WARNING ** : /ANIM/DROT OPTION USED',
2840 . ' WHILE ROTATIONAL DOF ARE NOT COMPUTED',
2841 . ' (IDROT = 0 IN /IOFLAG OPTION)'
2842 WRITE(istdo,*) ' ** WARNING ** : /ANIM/DROT OPTION USED',
2843 . ' WHILE ROTATIONAL DOF ARE NOT COMPUTED',
2844 . ' (IDROT = 0 IN /IOFLAG OPTION)'
2845 ENDIF
2846 IF(anim_v(14) == 1) THEN
2847 anim_v(14) = 0
2848 nv_ani = nv_ani - 1
2849 ENDIF
2850 ENDIF
2851C-----------------------
2852 IF(ispmd == 0.AND.mcheck == 0) THEN
2853C-----------------------
2854C Tracing Animations cards
2855C-----------------------
2856 WRITE(iout,1120)output%TANIM,output%DTANIM,output%TANIM_STOP,sensors%ANIM_ID,sensors%ANIM_DT,
2857 + anim_e(1),anim_e(2),anim_e(3),anim_e(25),
2858 + anim_e(4),anim_e(5),anim_e(6),
2859 + anim_e(7),anim_e(8),anim_e(9),
2860 + anim_e(10)+anim_e(4960)+anim_e(4961)+anim_e(4962)
2861 WRITE(iout,1129) anim_n(3),anim_n(4),anim_n(6)
2862 WRITE(iout,1130)
2863 + anim_v(1),anim_v(2),anim_v(3),anim_v(4),anim_v(5),
2864 + anim_v(6),anim_v(7),anim_v(9),anim_v(12),
2865 + anim_t(1),anim_t(2),anim_t(3),anim_t(4)
2866 WRITE(iout,1140)
2867 + anim_t(5),anim_t(6),anim_t(7),anim_t(8),
2868 + anim_m,anim_k
2869 ENDIF
2870C-------------------------------------------
2871c-----------------------------------------
2872c deallocate temporary SENSOR arrays
2873c-----------------------------------------
2874 IF (ALLOCATED(sensors%ANIM_TMP)) DEALLOCATE (sensors%ANIM_TMP)
2875 IF (ALLOCATED(sensors%STAT_TMP)) DEALLOCATE (sensors%STAT_TMP)
2876 IF (ALLOCATED(sensors%OUTP_TMP)) DEALLOCATE (sensors%OUTP_TMP)
2877 IF (ALLOCATED(sensors%STOP_TMP)) DEALLOCATE (sensors%STOP_TMP)
2878C---------------------------------------------------------
2879 IF(ierr == 0) RETURN
2880 CALL arret(0)
2881 1001 FORMAT(/
2882 & 1x,'ALE EULER SOLVERS'/
2883 & 1x,'-----------------'/)
2884 1002 FORMAT(
2885 & 1x,' +--STAGGERRED SCHEME'/
2886 & 1x,' | +--MOMENTUM : ',a16,' eta=',g14.7/
2887 & 1x,' | +--MASS : ',a16,' eta=',g14.7/
2888 & 1x,' | +--ENERGY : ',a16,' eta=',g14.7/
2889 & 1x,' | +--VOLUME FRACTION : ',a16)
2890 1003 FORMAT(
2891 & 1x,' |'/
2892 & 1x,' +--COLOCATED SCHEME (LAW151)'/
2893 & 1x,' | +--MOMENTUM : ',a16/
2894 & 1x,' | +--MASS : ',a16/
2895 & 1x,' | +--ENERGY : ',a16/
2896 & 1x,' | +--VOLUME FRACTION : ',a16 )
2897 1004 FORMAT(
2898 & 1x,' | +--LOW MACH OPTION : ENABLED' )
2899 1005 FORMAT(
2900 & 1x,' | +--MUSCL OPTION : ENABLED beta=',g14.7)
2901 1006 FORMAT(
2902 & 1x,' |'/
2903 & 1x,' +--TIME STEP'/
2904 & 1x,' | +--COURANT NUMBER : ',g14.7/
2905 & 1x,' | +--MINIMUM TIME STEP : ',g14.7)
2906 1007 FORMAT(
2907 & 1x,' |'/
2908 & 1x,' +--GRID SMOOTHING'/
2909 & 1x,' +--FORMULATION : ',a17)
2910 1008 FORMAT(
2911 & 1x,' +--PARAMETERS :')
2912
2913 1099 FORMAT(
2914 . ' MULTIDOMAINS COUPLING . . . . . . . . . . . . . .',g14.7//)
2915 1100 FORMAT(
2916 . ' FINAL TIME . . . . . . . . . . . . . . . . . . . ',g14.7//
2917 . ' TIME INTERVAL FOR TIME HISTORY PLOTS . . . . . . ',g14.7//
2918 . ' TIME STEP SCALE FACTOR . . . . . . . . . . . . . ',g14.7//
2919 . ' MINIMUM TIME STEP . . . . . . . . . . . . . . . ',g14.7//)
2920 1105 FORMAT(
2921 . ' BRICK TIME STEP SCALE FACTOR . . . . . . . . . . ',g14.7/
2922 . ' BRICK MINIMUM TIME STEP. . . . . . . . . . . . . ',g14.7/
2923 . ' MIN. TIME STEP FLAG (1:STOP RUN, 2:DELETE BRICK) ',i5//
2924 . ' QUAD TIME STEP SCALE FACTOR. . . . . . . . . . . ',g14.7/
2925 . ' QUAD MINIMUM TIME STEP . . . . . . . . . . . . . ',g14.7/
2926 . ' MIN. TIME STEP FLAG (1:STOP RUN, 2:DELETE QUAD). ',i5//
2927 . ' SHELL TIME STEP SCALE FACTOR . . . . . . . . . . ',g14.7/
2928 . ' SHELL MINIMUM TIME STEP. . . . . . . . . . . . . ',g14.7/
2929 . ' MIN. TIME FLAG (1:STOP, 2:DELETE, 3:SMALL STRAIN)',i5//
2930 . ' TRUSS TIME STEP SCALE FACTOR . . . . . . . . . . ',g14.7/
2931 . ' TRUSS MINIMUM TIME STEP. . . . . . . . . . . . . ',g14.7/
2932 . ' MIN. TIME FLAG (1:STOP, 2:DELETE). . . . . . . . ',i5//
2933 . ' BEAM TIME STEP SCALE FACTOR. . . . . . . . . . . ',g14.7/
2934 . ' BEAM MINIMUM TIME STEP . . . . . . . . . . . . . ',g14.7/
2935 . ' MIN. TIME FLAG (1:STOP, 2:DELETE). . . . . . . . ',i5//
2936 . ' SPRING TIME STEP SCALE FACTOR. . . . . . . . . . ',g14.7/
2937 . ' SPRING MINIMUM TIME STEP . . . . . . . . . . . . ',g14.7/
2938 . ' MIN. TIME FLAG (1:STOP, 2:DELETE). . . . . . . . ',i5//
2939 . ' AIRBAG TIME STEP SCALE FACTOR. . . . . . . . . . ',g14.7/
2940 . ' AIRBAG MINIMUM TIME STEP . . . . . . . . . . . . ',g14.7/
2941 . ' MIN. TIME FLAG (1:STOP). . . . . . . . . . . . . ',i5//
2942 . ' CONTACT TIME STEP SCALE FACTOR . . . . . . . . . ',g14.7/
2943 . ' CONTACT MINIMUM TIME STEP. . . . . . . . . . . . ',g14.7/
2944 . ' MIN. TIME FLAG(1:STOP, 2:REMOVE NODE FROM INTERF)',i5/ )
2945 1155 FORMAT('/DT/FVMBAG/0 OPTION')
2946 1156 FORMAT('/DT/FVMBAG/1 OPTION')
2947 1157 FORMAT('FVMBAG TIME STEP SYNTHESIS')
2948 1147 FORMAT(
2949 . ' FVMBAG ID. . . . . . . . . . . . . . . . . . . . ',i10)
2950 1148 FORMAT(
2951 . ' FVMBAG1 ID. . . . . . . . . . . . . . . . . . . . ',i10)
2952 1149 FORMAT(
2953 . ' FVMBAG2 ID. . . . . . . . . . . . . . . . . . . . ',i10)
2954 1151 FORMAT(
2955 . ' SMOOTH PARTICLES TIME STEP SCALE FACTOR. . . . . ',g14.7/
2956 . ' SMOOTH PARTICLES MINIMUM TIME STEP . . . . . . . ',g14.7/
2957 . ' MIN. TIME FLAG (1:STOP, 2:DELETE, 5:KILL). . . . ',i5 )
2958 1152 FORMAT(
2959 . ' FVMBAG TIME STEP SCALE FACTOR. . . . . . . . . . ',g14.7/
2960 . ' FVMBAG MINIMUM TIME STEP . . . . . . . . . . . . ',g14.7/
2961 . ' FVMBAG TIME STEP FLAG. . . . . . . . . . . . . . ',i5,/)
2962 1153 FORMAT(
2963 . ' . . TIME STEP SCALE FACTOR. . . . . . . . . . . . . . ',g14.7/
2964 . ' . . MINIMUM TIME STEP . . . . . . . . . . . . . . . . ',g14.7/
2965 . ' . . TIME STEP FLAG. . . . . . . . . . . . . . . . . . ',i5)
2966 1154 FORMAT(
2967 . ' . . CHARACTERISTIC LENGTH OPTION. . . . . . . . . . . ',i5/
2968 . ' . . TIME STEP SMOOTHING FACTOR. . . . . . . . . . . . ',g14.7/
2969 . ' . . PREVIOUS TIME STEP. . . . . . . . . . . . . . . . ',g14.7)
2970 1107 FORMAT(
2971 . ' CONTACT NODAL TIME STEP SCALE FACTOR . . . . . . ',g14.7/
2972 . ' CONTACT NODAL MINIMUM TIME STEP. . . . . . . . . ',g14.7/
2973 . ' CONTACT NODAL MIN. TIME FLAG . . . . . . . . . . ',i5/
2974 . ' 3:INCREASE MASS, ORIGINAL FORMULATION . . . . ',/
2975 . ' 8:INCREASE MASS, IMPROVED FORMULATION . . . . ',//)
2976 1209 FORMAT(
2977 . ' AMS CONTACT TIME STEP IS ON . . . . . . . . . . . . .',/
2978 . ' . . . . . .(ALL CONTACTS WILL BE CONCERNED).) . . . .',/
2979 . ' AMS CONTACT TIME STEP SCALE FACTOR. . . . . . . . . .',g14.7/
2980 . ' AMS CONTACT MINIMUM TIME STEP . . . . . . . . . . . .',g14.7//)
2981 1106 FORMAT(
2982 . ' NODAL TIME STEP SCALE FACTOR, . . . . . . . . . ',g14.7/
2983 . ' NODAL MINIMUM TIME STEP. . . . . . . . . . . . . ',g14.7/
2984 . ' MIN. TIME FLAG . . . . . . . . . . . . . . . . . ',i5/
2985 . ' 1:STOP. . . . . . . . . . . . . . . . . . . . ',/
2986 . ' 3:INCREASE MASS, ORIGINAL FORMULATION . . . . ',/
2987 . ' 8:INCREASE MASS, IMPROVED FORMULATION . . . . ',//)
2988 1116 FORMAT(
2989 . ' NODAL TIME STEP SCALE FACTOR, . . . . . . . . . ',g14.7/
2990 . ' NODAL MINIMUM TIME STEP. . . . . . . . . . . . . ',g14.7/)
2991 1206 FORMAT(
2992 . ' NODAL TIME STEP SCALE FACTOR, . . . . . . . . . ',g14.7/
2993 . ' NODAL MINIMUM TIME STEP. . . . . . . . . . . . . ',g14.7/
2994 . ' ADDED MASS RATIO (DM/M0) . . . . . . . . . . . . ',g14.7/
2995 . ' MIN. TIME FLAG . . . . . . . . . . . . . . . . . ',i5/
2996 . ' 1:STOP. . . . . . . . . . . . . . . . . . . . ',/
2997 . ' 3:INCREASE MASS, ORIGINAL FORMULATION . . . . ',/
2998 . ' 8:INCREASE MASS, IMPROVED FORMULATION . . . . ',//)
2999 1108 FORMAT(
3000 . ' ADVANCED MASS SCALING IS ON . . . . . . . . . . .',/
3001 . ' AMS NODAL TIME STEP SCALE FACTOR. . . . . . . . .',g14.7/
3002 . ' AMS NODAL MINIMUM TIME STEP . . . . . . . . . . .',g14.7/
3003 . ' AMS TOLERANCE ON CONVERGENCE. . . . . . . . . . .',g14.7/
3004 . ' AMS MAXIMUM NUMBER OF ITERATIONS FOR CONVERGENCE.',i10/
3005 . ' OUTPUT FREQUENCY OF INFORMATION / CONVERGENCE . .',i10/
3006 . ' PART GROUP ID (=0 ALL PARTS) . . . . . . . . . . ',i10//)
3007 1109 FORMAT(
3008 . ' ADVANCED MASS SCALING IS ON . . . . . . . . . . .',/
3009 . ' AMS TIME STEP SCALE FACTOR. . . . . . . . . . . .',g14.7/
3010 . ' AMS MINIMUM TIME STEP . . . . . . . . . . . . . .',g14.7/
3011 . ' AMS TOLERANCE ON CONVERGENCE. . . . . . . . . . .',g14.7/
3012 . ' AMS MAXIMUM NUMBER OF ITERATIONS FOR CONVERGENCE.',i10/
3013 . ' NUMBER OF PRECONDITIONING VECTORS FOR PCG . . . .',i10/
3014 . ' OUTPUT FREQUENCY OF INFORMATION / CONVERGENCE . .',i10/
3015 . ' PART GROUP ID (=0 ALL PARTS) . . . . . . . . . . ',i10//)
3016 2109 FORMAT(
3017 . ' ADVANCED MASS SCALING IS ON . . . . . . . . . . .',/
3018 . ' AMS TIME STEP SCALE FACTOR. . . . . . . . . . . .',g14.7/
3019 . ' AMS MINIMUM TIME STEP . . . . . . . . . . . . . .',g14.7/
3020 . ' AMS TOLERANCE ON CONVERGENCE. . . . . . . . . . .',g14.7/
3021 . ' AMS MAXIMUM NUMBER OF ITERATIONS FOR CONVERGENCE.',i10/
3022 . ' NUMBER OF PRECONDITIONING VECTORS FOR PCG . . . .',i10/
3023 . ' OUTPUT FREQUENCY OF INFORMATION / CONVERGENCE . .',i10/
3024 . ' AUTO ELEMENT SELECTION - TIME STEP CRITERIA . . .',g14.7/
3025 . ' PART GROUP ID (=0 ALL PARTS) . . . . . . . . . . ',i10//)
3026 1110 FORMAT(
3027 . ' PRINTOUT CYCLE FREQUENCY . . . . . . . . . . . . ',i5 //
3028 . ' RESTART CYCLE FREQUENCY. . . . . . . . . . . . . ',i10 //
3029 . ' MAXIMUM RESTART WRITE BEFORE OVERWRITE . . . . . ',i5 //
3030 . ' INVERS: INPUT DECK VERSION . . . . . . . . . . . ',i5/
3031 . ' ITTYP : FLAG FOR TYPE OF T-FILE. . . . . . . . . ' ,i5/)
3032 1120 FORMAT(/
3033 . ' TIME FOR FIRST ANIMATION-FILE PLOT . . . . . . . ',g14.7/
3034 . ' TIME INTERVAL FOR ANIMATION-FILE PLOTS . . . . . ',g14.7/
3035 . ' TIME TO STOP ANIMATION-FILE PLOTS. . . . . . . . ',g14.7/
3036 . ' SENSOR FOR ANIMATION-FILE PLOTS. . . . . . . . . ',i10/
3037 . ' TIME INTERVAL FOR SENSOR ANIMATION-FILE PLOTS. . ',g14.7//
3038 . ' VARIABLE SAVED ON ANIMATION FILES (1: YES) :',/
3039 . ' PLASTIC STRAIN. . . . . . ',i5/
3040 . ' DENSITY (SOLID ONLY). ',i5/
3041 . ' SPECIFIC ENERGY . . . . . ',i5/
3042 . ' HOURGLASS ENERGY. . . . . ',i5/
3043 . ' TEMPERATURE (SOLID ONLY). ',i5/
3044 . ' THICKNESS (SHELL ONLY). ',i5/
3045 . ' PRESSURE (SOLID ONLY). ',i5/
3046 . ' VON MISES . . . . . . . . ',i5/
3047 . ' TURBULENT ENERGY (FLUID). ',i5/
3048 . ' TURBULENT VISCOSITY(FLUID)',i5/
3049 . ' VORTICITY (FLUID) . . . . ',i5)
3050 1129 FORMAT(
3051 . ' NODAL VALUES SAVED ON ANIMATION FILES (1: YES) :',/
3052 . ' PRESSURE. . . . . . . . . ',i5/
3053 . ' DENSITY . . . . . . . . . ',i5/
3054 . ' TEMPERATURE . . . . . . . ',i5)
3055 1130 FORMAT(
3056 . ' VECTOR SAVED ON ANIMATION FILES (1: YES) :',/
3057 . ' VELOCITY VECTOR . . . . . ',i5/
3058 . ' DISPLACEMENT VECTOR . . . ',i5/
3059 . ' ACCELERATION VECTOR . . . ',i5/
3060 . ' CONTACT FORCES. . . . . . ',i5/
3061 . ' INTERNAL FORCES . . . . . ',i5/
3062 . ' EXTERNAL FORCES . . . . . ',i5/
3063 . ' SECTION RBY RWALL FORCES. ',i5/
3064 . ' ROTATIONAL VELOCITY VECTOR',i5/
3065 . ' CONTACT PRESSURE (VECTORS)',i5/
3066 . ' SHELL TENSOR SAVED ON ANIMATION FILES (1: YES) :',/
3067 . ' MEMBRANE STRESS . . . . . ',i5/
3068 . ' BENDING STRESS(MOMENT/T^2)',i5/
3069 . ' UPPER LAYER STRESS. . . . ',i5/
3070 . ' LOWER LAYER STRESS. . . . ',i5)
3071 1140 FORMAT(
3072 . ' MEMBRANE STRAIN . . . . . ',i5/
3073 . ' CURVATURE . . . . . . . . ',i5/
3074 . ' UPPER LAYER STRAIN. . . . ',i5/
3075 . ' LOWER LAYER STRAIN. . . . ',i5/
3076 . ' NODAL MASS SAVED ON ANIMATION FILES (1: YES) :',i5/
3077 . ' KEEP DELETED ELEMENT (1: YES) :',i5//)
3078 1150 FORMAT(
3079 . ' INITIAL TIME STEP. . . . . . . . . . . . . . . . ',g14.7/
3080 . ' MAXIMUM TIME STEP. . . . . . . . . . . . . . . . ',g14.7)
3081 1160 FORMAT(/
3082 . ' IMPLICIT : CONJUGATED GRADIENT ' /
3083 . ' GLOBAL CONVERGENCE PRECISION . . . . . . . . . . ',g14.7/
3084 . ' INCREMENTAL CONVERGENCE PRECISION. . . . . . . . ',g14.7/
3085 . ' MAXIMUM NUMBER OF ITERATIONS . . . . . . . . . . ',i5)
3086 1171 FORMAT(/
3087 . ' DYNAMIC RELAXATION ' /
3088 . ' NODE GROUP ID (=0 ALL NODES) . . . . . . . . . . ',i10/
3089 . ' BETA . . . . . . . . . . . . . . . . . . . . . . ',g14.7/
3090 . ' PERIOD . . . . . . . . . . . . . . . . . . . . . ',g14.7)
3091 1172 FORMAT(/
3092 . ' KINEMATIC RELAXATION ' /
3093 . ' NODE GROUP ID (=0 ALL NODES) . . . . . . . . . . ',i10)
3094 1180 FORMAT(/
3095 . ' PARALLEL ARITHMETIC OFF')
3096 1181 FORMAT(/
3097 . ' PARALLEL ARITHMETIC ON')
3098 1182 FORMAT(/
3099 . ' PARALLEL ARITHMETIC FLAG . . . . . . . . . . . . ',i5)
3100 1196 FORMAT(/
3101 . ' QUASI-COMPRESSIBLE FORMULATION ON (/INCMP)')
3102 1198 FORMAT(
3103 . ' MOMENTUM : MIXED INTEGRATION')
3104 1199 FORMAT(
3105 . ' MOMENTUM : VOLUME INTEGRATION')
3106 1200 FORMAT(
3107 . 28x,' ALPHA : DONEA COEFFICIENT. . . . . . . . . . . ',g14.7/
3108 . 28x,' GAMMA : GRID VELOCITY LIMITATION FACTOR. . . . ',g14.7/
3109 . 28x,' FscaleX : X-GRID VELOCITY SCALE FACTOR . . . . . ',g14.7/
3110 . 28x,' FscaleY : Y-GRID VELOCITY SCALE FACTOR . . . . . ',g14.7/
3111 . 28x,' FscaleZ : Z-GRID VELOCITY SCALE FACTOR . . . . . ',g14.7/
3112 . 28x,' VOLMIN : MINIMUM VOLUME FOR ELEMENT DELETION. . ',g14.7//)
3113 1220 FORMAT(
3114 . 28x,' UMAX : MAXIMUM ABSOLUTE GRID VELOCITY . . . . . ',g14.7/
3115 . 28x,' VMIN : MINIMUM VOLUME FOR ELEMENT DELETION. . . ',g14.7//)
3116 1250 FORMAT(
3117 . 28x,' DT0 : TYPICAL TIME STEP. . . . . . . . . . . . ',g14.7/
3118 . 28x,' DT0* : EFFECTIVE TIME STEP. . . . . . . . . . . ',g14.7/
3119 . 28x,' GAMMA : NON LINEARITY FACTOR . . . . . . . . . . ',g14.7/
3120 . 28x,' ETA : DAMPING COEFFICIENT . . . . . . . . . . ',g14.7/
3121 . 28x,' NU : SHEAR FACTOR . . . . . . . . . . . . . . ',g14.7/
3122 . 28x,' VOLMIN: MINIMUM VOLUME FOR ELEMENT DELETION. . . ',g14.7//)
3123 1254 FORMAT(
3124 . 28x,' ALPHA : STABILITY FACTOR . . . . . . . . . . . . ',g14.7/
3125 . 28x,' GAMMA : NON LINEARITY FACTOR . . . . . . . . . . ',g14.7/
3126 . 28x,' BETA : DAMPING COEFFICIENT. . . . . . . . . . . ',g14.7/
3127 . 28x,' LC : CHARACTERISTIC LENGTH. . . . . . . . . . ',g14.7//)
3128 1257 FORMAT(
3129 . 28x,' ENABLED DEFORMATION . . . . . . . . . . . . . . ',a3/
3130 . 28x,' ENABLED ROTATION . . . . . . . . . . . . . . . ',a3/
3131 . 28x,' SCALE FACTOR FOR DEFORMATION . . . . . . . . . . ',g14.7/
3132 . 28x,' SCALE FACTOR FOR ROTATION . . . . . . . . . . . ',g14.7//)
3133 1300 FORMAT(
3134 . ' NUMBER OF INTERFACES TO BE ELIMINATED. . . . . . ',i8//
3135 . ' NUMBER OF PARTS TO BE ELIMINATED . . . . . . . . ',i8//
3136 . ' NUMBER OF SOLID ELEMENT BLOCKS TO BE ELIMINATED ',i8//
3137 . ' NUMBER OF QUAD ELEMENT BLOCKS TO BE ELIMINATED ',i8//
3138 . ' NUMBER OF SHELL ELEMENT BLOCKS TO BE ELIMINATED ',i8//
3139 . ' NUMBER OF TRUSS ELEMENT BLOCKS TO BE ELIMINATED ',i8//
3140 . ' NUMBER OF BEAM ELEMENT BLOCKS TO BE ELIMINATED ',i8//
3141 . ' NUMBER OF SPRING ELEMENT BLOCKS TO BE ELIMINATED ',i8//
3142 . ' NUMBER OF SH_3N ELEMENT BLOCKS TO BE ELIMINATED ',i8//
3143 . ' NUMBER OF SPH PARTICLES BLOCKS TO BE ELIMINATED ',i8/)
3144 1400 FORMAT(
3145 . ' NUMBER OF RIGID LINKS. . . . . . . . . . . . . . ',i8/)
3146 1450 FORMAT(
3147 . ' NUMBER OF ALE LINKS. . . . . . . . . . . . . . . ',i8/)
3148 1500 FORMAT(
3149 . ' NUMBER OF NEW BOUNDARY CONDITIONS. . . . . . . . ',i8/)
3150 1550 FORMAT(
3151 . ' REMOVE INTER.7 SEGMENT AFTER SHELL FAILURE(1 YES)',i5/)
3152 1810 FORMAT(///' LIST OF ELIMINATED SOLID ELEMENTS ')
3153 1820 FORMAT(///' LIST OF ELIMINATED QUAD ELEMENTS ')
3154 1830 FORMAT(///' LIST OF ELIMINATED SHELL ELEMENTS ')
3155 1840 FORMAT(///' LIST OF ELIMINATED TRUSS ELEMENTS ')
3156 1850 FORMAT(///' LIST OF ELIMINATED BEAM ELEMENTS ')
3157 1860 FORMAT(///' LIST OF ELIMINATED SPRING ELEMENTS ')
3158 1870 FORMAT(///' LIST OF ELIMINATED SH_3N ELEMENTS ')
3159 1880 FORMAT(///' LIST OF ELIMINATED SPH PARTICLES ')
3160 1890 FORMAT(///' LIST OF ELIMINATED PARTS ')
3161 2000 FORMAT(/' INTERFACE NUMBER',i10,' IS ELIMINATED')
3162 2100 FORMAT(/' RIGID LINK:',i5,
3163 . ' TRANSLATION X,Y,Z',3(1x,i1),
3164 . ' ROTATION X,Y,Z',3(1x,i1),
3165 . ' SKEW ',i10,/
3166 . ' ---------- ',i10,' NODES :')
3167 2200 FORMAT(/' ALE LINK: ',/,
3168 . ' --------',/,
3169 . ' MAIN NODES : ',2i8,/,
3170 . ' (X,Y,Z) : (',i1,',',i1,',',i1,')',/,
3171 . ' TYPE : ',i2)
3172 2201 FORMAT( ' NODES :')
3173 2202 FORMAT( ' NODE GROUP : ',i2)
3174 2210 FORMAT(' DEACTIVATING ALE LINK ID:',i5)
3175 2211 FORMAT(' ACTIVATING ALE LINK ID:',i5)
3176 2300 FORMAT(/' INTERFACE CHANGES'/
3177 . ' INTERFACE NB . . . . . . . . . . . . . . . . . . ',i10/
3178 . ' SEARCH OF CLOSEST NODES EACH NSEARCH TIME STEPS. ',i5/
3179 . ' START TIME . . . . . . . . . . . . . . . . . . . ',g14.7/
3180 . ' STOP TIME . . . . . . . . . . . . . . . . . . . ',g14.7)
3181
3182 2301 FORMAT(/' ERROR IN INTERFACE CHANGES'/
3183 . ' INTERFACE NB IS NOT EXISTING . . . . . . . . . . ',i10)
3184
3185 4000 FORMAT(/
3186 . ' SPMD PARALLEL RESTART READING (PREAD) ACTIVATED')
3187 4100 FORMAT(/
3188 . ' SPMD ELEMENT WEIGHT ESTIMATION (DDW) ACTIVATED')
3189C
3190 4500 FORMAT(/
3191 & 1x,'EXTERNAL LIBRARY FOR USERS CODE INTERFACE ')
3192 4600 FORMAT(
3193 & 1x,'LIBRARY NAME . . . . . . . . . . . . . . . . . . . . ',a/
3194 & 1x,'RADIOSS USERS CODE INTERFACE VERSION . . . . . . . .',i10//)
3195 4700 FORMAT(
3196 . ' BRICK TIME STEP SCALE FACTOR . . . . . . . . . . ',g14.7/
3197 . ' BRICK MINIMUM TIME STEP. . . . . . . . . . . . . ',g14.7/
3198 . ' MIN. TIME STEP FLAG (1:STOP RUN, 2:DELETE 3:CST )',i5//
3199 . ' BRICK_CST MINIMUM ASPECT RATIO (Tet collapse). . ',g14.7/
3200 . ' BRICK_CST MINIMUM VOLUME CHANGE . . . . . . . . ',g14.7//
3201 . ' QUAD TIME STEP SCALE FACTOR. . . . . . . . . . . ',g14.7/
3202 . ' QUAD MINIMUM TIME STEP . . . . . . . . . . . . . ',g14.7/
3203 . ' MIN. TIME STEP FLAG (1:STOP RUN, 2:DELETE QUAD). ',i5//
3204 . ' SHELL TIME STEP SCALE FACTOR . . . . . . . . . . ',g14.7/
3205 . ' SHELL MINIMUM TIME STEP. . . . . . . . . . . . . ',g14.7/
3206 . ' MIN. TIME FLAG (1:STOP, 2:DELETE, 3:SMALL STRAIN)',i5//
3207 . ' TRUSS TIME STEP SCALE FACTOR . . . . . . . . . . ',g14.7/
3208 . ' TRUSS MINIMUM TIME STEP. . . . . . . . . . . . . ',g14.7/
3209 . ' MIN. TIME FLAG (1:STOP, 2:DELETE). . . . . . . . ',i5//
3210 . ' BEAM TIME STEP SCALE FACTOR. . . . . . . . . . . ',g14.7/
3211 . ' BEAM MINIMUM TIME STEP . . . . . . . . . . . . . ',g14.7/
3212 . ' MIN. TIME FLAG (1:STOP, 2:DELETE). . . . . . . . ',i5//
3213 . ' SPRING TIME STEP SCALE FACTOR. . . . . . . . . . ',g14.7/
3214 . ' SPRING MINIMUM TIME STEP . . . . . . . . . . . . ',g14.7/
3215 . ' MIN. TIME FLAG (1:STOP, 2:DELETE). . . . . . . . ',i5//
3216 . ' AIRBAG TIME STEP SCALE FACTOR. . . . . . . . . . ',g14.7/
3217 . ' AIRBAG MINIMUM TIME STEP . . . . . . . . . . . . ',g14.7/
3218 . ' MIN. TIME FLAG (1:STOP). . . . . . . . . . . . . ',i5//
3219 . ' INTERFACE TYPE 7 TIME STEP SCALE FACTOR. . . . . ',g14.7/
3220 . ' INTERFACE TYPE 7 MINIMUM TIME STEP . . . . . . . ',g14.7/
3221 . ' MIN. TIME FLAG(1:STOP, 2:REMOVE NODE FROM INTERF)',i5/ )
3222 4720 FORMAT(
3223 . ' ACCURATE TIME STEP FOR SHELLS IS USED. . . . . . ',/)
3224 4730 FORMAT(
3225 . ' ACCURATE TIME STEP FOR 4-NODE & 8-NODE SOLID IS USED',/)
3226 4740 FORMAT(
3227 . ' ACCURATE TIME STEP FOR 10-NODE TETRAHEDRA IS USED ',/)
3228 4800 FORMAT(
3229 . ' PART DEACTIVATED FOR ALE / EULER COMPUTATION',2x,i10)
3230 5000 FORMAT(/' H3D FILES : '/
3231 . ' TIME FOR FIRST H3D-FILE PLOT . . . . . . . . . . ',g14.7/
3232 . ' TIME INTERVAL FOR H3D-FILE PLOTS . . . . . . . . ',g14.7/
3233 . ' VARIABLE SAVED ON H3D FILES : ')
3234 5001 FORMAT(/
3235 . ' KINEMATIC RELAXATION ' /
3236 . ' NODE GROUP ID (=0 ALL NODES) . . . . . . . . . . ',i10/
3237 . ' START TIME. . . . . . . . . . . . . . . . . . . ',g14.7/
3238 . ' STOP TIME. . . . . . . . . . . . . . . . . . . ',g14.7)
3239 5010 FORMAT(/
3240 . ' ADAPTIVE DYNAMIC RELAXATION ' ,/
3241 . ' NODE GROUP ID (=0 ALL NODES) . . . . . . . . . . ',i10/)
3242 5011 FORMAT(/
3243 . ' ADAPTIVE DYNAMIC RELAXATION ' /
3244 . ' NODE GROUP ID (=0 ALL NODES) . . . . . . . . . . ',i10/
3245 . ' START TIME. . . . . . . . . . . . . . . . . . . ',g14.7/
3246 . ' STOP TIME. . . . . . . . . . . . . . . . . . . ',g14.7)
3247 5020 FORMAT(
3248 . ' BRICK_DEL MINIMUM COLLAPSE RATIO. . . . . . . . . ',g14.7/
3249 . ' BRICK_DEL MINIMUM VOLUME CHANGE . . . . . . . . . ',g14.7/
3250 . ' BRICK_DEL MAXIMUM ASPECT RATIO . . . . . . . . . ',g14.7/
3251 . ' BRICK_DEL MAXIMUM VOLUME CHANGE . . . . . . . . . ',g14.7//)
3252
3253 END
3254!||====================================================================
3255!|| prout_buck ../engine/source/input/lectur.F
3256!||--- called by ------------------------------------------------------
3257!|| lecimpl ../engine/source/input/lectur.F
3258!||--- uses -----------------------------------------------------
3259!|| imp_kbcs ../engine/share/modules/impbufdef_mod.F
3260!||====================================================================
3261 SUBROUTINE prout_buck(IP,NBUCK,IBUCK)
3262C-----------------------------------------------
3263C M o d u l e s
3264C-----------------------------------------------
3265 USE imp_kbcs
3266C-----------------------------------------------
3267C I m p l i c i t T y p e s
3268C-----------------------------------------------
3269#include "implicit_f.inc"
3270C-----------------------------------------------
3271C D u m m y A r g u m e n t s
3272C-----------------------------------------------
3273 INTEGER IP,NBUCK,IBUCK
3274C-----------------------------------------------
3275C L o c a l V a r i a b l e s
3276C-----------------------------------------------
3277 CHARACTER*25 MSG_BSOL(2)
3278 DATA
3279 . msg_bsol
3280 . / 'BCS',
3281 . '-' /
3282C-----------------------------------------------
3283 IF (ibuck > 0) THEN
3284 WRITE(ip,1000)
3285 ELSE
3286 WRITE(ip,2000)
3287 END IF
3288 WRITE(ip,3000)nbuck,shift_b,emin_b,emax_b,msgl_b,maxset_b, msg_bsol(1)
3289
3290 RETURN
3291 1000 FORMAT(
3292 . ' EULER BUCKLING ANALYSIS (RESTART):'/)
3293 2000 FORMAT(
3294 . ' EULER BUCKLING ANALYSIS :'/)
3295 3000 FORMAT(
3296 . ' NUMBER OF MODES TO BE COMPUTED: . . . . . ',2x,i5/
3297 . ' SHIFT IN BUCKLING MODES PENCIL: . . . . . ',2x,g14.7/
3298 . ' MINIMUM EIGENVALUE: . . . . . . . . . . . ',2x,g14.7/
3299 . ' MAXIMUM EIGENVALUE: . . . . . . . . . . . ',2x,g14.7/
3300 . ' OUTPUT MESSAGE LEVEL: . . . . . . . . . . ',2x,i5/
3301 . ' NUMBER OF VECTORS IN BLOCK OR SET: . . . ',2x,i5/
3302 . ' LINEAR SOLVER: . . . . . . . . . . . . . ',2x,a/)
3303 END
3304!||====================================================================
3305!|| lecimpl ../engine/source/input/lectur.F
3306!||--- called by ------------------------------------------------------
3307!|| lectur ../engine/source/input/lectur.F
3308!||--- calls -----------------------------------------------------
3309!|| ancmsg ../engine/source/output/message/message.F
3310!|| arret ../engine/source/system/arret.F
3311!|| floatmin ../common_source/tools/math/precision.c
3312!|| prout_buck ../engine/source/input/lectur.F
3313!|| spmd_max_s ../engine/source/mpi/implicit/imp_spmd.F
3314!||--- uses -----------------------------------------------------
3315!|| imp_kbcs ../engine/share/modules/impbufdef_mod.F
3316!|| imp_pcg_proj ../engine/share/modules/impbufdef_mod.F
3317!|| imp_spbrm ../engine/share/modules/impbufdef_mod.F
3318!|| message_mod ../engine/share/message_module/message_mod.F
3319!||====================================================================
3320 SUBROUTINE lecimpl
3321C-----------------------------------------------
3322C M o d u l e s
3323C-----------------------------------------------
3324 USE imp_kbcs
3325 USE imp_pcg_proj
3326 USE imp_spbrm
3327 USE message_mod
3328C-----------------------------------------------
3329C I m p l i c i t T y p e s
3330C-----------------------------------------------
3331#include "implicit_f.inc"
3332C-----------------------------------------------
3333C C o m m o n B l o c k s
3334C-----------------------------------------------
3335#include "com01_c.inc"
3336#include "com04_c.inc"
3337#include "com06_c.inc"
3338#include "units_c.inc"
3339#include "task_c.inc"
3340#include "impl1_c.inc"
3341#include "impl2_c.inc"
3342#include "buckcom.inc"
3343C-----------------------------------------------
3344C L o c a l V a r i a b l e s
3345C-----------------------------------------------
3346 INTEGER J,NTY,IPRINT,IBID,ISOL
3347 my_real CS1(2),S
3348 REAL FLMIN
3349 CHARACTER*25 MSG_TYPE(9),MSG_ISOL(9),MSG_INSOL(4),MSG_PRE(5),MSG_BSOL(3)
3350 DATA
3351 . MSG_TYPE
3352 . / 'STATIC LINEAR',
3353 . 'STATIC NON-LINEAR',
3354 . 'DYNAMIC NON-LINEAR',
3355 . 'MODAL ANALYSIS',
3356 . 'CHECK',
3357 . 'QUASI-STATIC LINEAR',
3358 . 'QUASI-STATIC NON-LINEAR',
3359 . 'BUCKLING',
3360 . 'RADIOSS-AcuSolve DC-FSI'/,
3361 . MSG_ISOL
3362 . / 'PREC. CONJUGATE GRADIENT',
3363 . 'DIRECT(MUMPS) ',
3364c .'Prec.Lanzos',
3365 . 'DIRECT',
3366 . 'MULTI-LEVEL CONDENSATION',
3367 . 'MIXE SOLVERS: 1 + 3',
3368 . 'MIXE SOLVERS: 1 + 4',
3369 . 'AUTO SELECT SOLVER ',
3370 . 'PCG(GPU) ',
3371 . 'PCGP(PROJECTION) '/,
3372 . msg_insol
3373 . / 'MODIFIED NEWTON',
3374 . 'QUASI-NEWTON',
3375 . 'SLOAN ELASTOPLAS',
3376 . 'FULL NEWTON'/,
3377 . msg_pre
3378 . / 'NONE',
3379 . 'DIAGONAL JACOBIEN',
3380 . 'IMCOMP. CHOLESKY(0)',
3381 . 'STABLILIZED IC(0)',
3382 . 'FACTORED APPRO. INVERSES'/,
3383 . msg_bsol
3384 . / 'SUPERLU',
3385 . 'MUMPS',
3386 . 'MULTI-LEVEL CONDENSATION'/
3387C-----------------------------------------------
3388 iprint=0
3389C IMPMV>0 isolv/=1 --> IMPMV=0
3390 IF (ispmd == 0) iprint=1
3391 IF (imp_chk > 0) THEN
3392 nty = 5
3393 IF (iprec /= 5) iprec=5
3394 IF (iline /= 1) iline=1
3395 IF (isolv /= 1) isolv=1
3396 IF (d_tol /= zero) d_tol = zero
3397 ENDIF
3398 IF(isolv >= 3 ) THEN
3399 CALL ancmsg(msgid=296,anmode=aninfo,i1=isolv)
3400 isolv = 2
3401 END IF
3402
3403#ifndef MUMPS5
3404 IF(isolv == 2 .AND. nspmd > 1) THEN
3405 WRITE(6,*) "Fatal error: MUMPS is required, but not available"
3406 CALL arret(5)
3407 ENDIF
3408#endif
3409 IF (neig > 0) THEN
3410 iline=1
3411 isprb=0
3412 iqstat=0
3413 IF (nspmd == 1) THEN
3414 iprec=1
3415 ELSE
3416 iprec=5
3417 ENDIF
3418 ENDIF
3419 IF (nbuck > 0) THEN
3420 iline=1
3421C----fix mono w/ BCS solver
3422 IF (nspmd == 1) isolv =2
3423 ENDIF
3424C--------solvers----
3425C machine precision minimum -simple
3426 CALL floatmin(cs1(1),cs1(2),flmin)
3427 p_mach = two*sqrt(flmin)
3428 IF (nspmd > 1)CALL spmd_max_s(p_mach)
3429C AUTO SELECT SOLVER
3430 IF (isolv == 7 ) THEN
3431
3432 IF (iline == 1) THEN
3433
3434 WRITE(iout,*) ' ** WARNING ** : SOLVER AUTO SELECT IS NOT ','COMPATIBLE WITH LINEAR RUN '
3435 WRITE(iout,*) ' ** RESETTING TO ** : DEFAULT ONE '
3436
3437 WRITE(istdo,*) ' ** WARNING ** : SOLVER AUTO SELECT IS NOT ','COMPATIBLE WITH LINEAR RUN '
3438 WRITE(istdo,*) ' ** RESETTING TO ** : DEFAULT ONE '
3439
3440 isolv = 0
3441
3442 END IF
3443
3444 ENDIF
3445C----- default solver-> MUMPS
3446 IF (isolv == 0) isolv = 2
3447C-------ISOLV=2 --> use MUMPS anyway
3448 imumpsv = 0
3449 IF (isolv > 1 .AND. isolv <9 ) imumpsv = 1
3450 IF (isolv == 2 .OR. isolv ==3 ) THEN
3451 imumpsv = 1
3452 isolv = 3
3453 END IF
3454 IF ((isolv == 3.OR.isolv == 4).AND.intp_c < 0) THEN
3455 IF(iprint==1) THEN
3456 WRITE(iout,*) ' ** WARNING ** : DIRECT SOLVER IS NOT ','COMPATIBLE WITH IMPL/INTER/KNONL OPTION '
3457 WRITE(iout,*) ' ** RESETTING TO ** : MIXE ONE '
3458 ENDIF
3459 isolv = isolv + 2
3460 ENDIF
3461C-------ISOLV=9 -> PCG w/ Projection
3462 IF (isolv == 9) THEN
3463 IF (m_vs ==0) m_vs=20
3464 IF (ipro_s0 ==0) ipro_s0=4
3465 END IF
3466 IF (m_vs > 0) THEN
3467C------case /IMPL/PROJV/n w/o /SOLV/9
3468 IF (isolv /= 9) THEN
3469 WRITE(iout,*) ' ** WARNING ** : /IMPL/PROJV IS ','ONLY COMPATIBLE WITH PCGP SOLVER '
3470 WRITE(iout,*) ' ** CHANGE TO ** : ISOLV=9 '
3471 isolv = 9
3472 END IF
3473 END IF
3474
3475 IF (nbuck > 0.AND.nspmd == 1.AND.isolv /= 3) THEN
3476 WRITE(iout,*) ' ** WARNING ** : /IMPL/BUCKL IS ','ONLY COMPATIBLE WITH BCS SOLVER '
3477 WRITE(iout,*) ' ** CHANGE TO ** : ISOLV=3 '
3478 isolv = 3
3479 ENDIF
3480
3481 IF (iprec == 0.OR.iprec > 6)iprec=5
3482
3483 IF (isolv == 2) THEN
3484 IF (itol == 0.OR.itol > 1) itol=2
3485 IF (l_tol == zero) l_tol=p_mach
3486 ELSEIF (isolv == 1.OR.isolv>=7) THEN
3487 IF (itol == 0) itol=3
3488 IF (itol > 4) THEN
3489 IF(iprint==1) THEN
3490 WRITE(iout,*) ' ** WARNING ** : WRONG LINEAR STOP ','CRITERION NUMBER, RESET TO DEFAULT ONE ** '
3491 ENDIF
3492 itol=3
3493 ENDIF
3494 IF (l_tol == zero) THEN
3495 IF (itol == 3) THEN
3496 l_tol=p_mach
3497 IF (n_pat > 1.OR.isolv == 9) l_tol=p_mach*em01
3498 ELSE
3499 l_tol=em5
3500 ENDIF
3501 ENDIF
3502C--------direct-----
3503 ELSEIF (isolv == 3.OR.isolv == 4) THEN
3504 iprec=1
3505 IF (nspmd>1) imumpsv = 1
3506C--------MIX-----
3507 ELSEIF (isolv == 5.OR.isolv == 6) THEN
3508 iprec=1
3509 imumpsv = 1
3510 IF (itol == 0) itol=1
3511 IF (itol > 4) THEN
3512 IF(iprint==1) THEN
3513 WRITE(iout,*) ' ** WARNING ** : WRONG LINEAR STOP ','CRITERION NUMBER, RESET TO DEFAULT ONE ** '
3514 ENDIF
3515 itol=3
3516 ENDIF
3517 IF (l_tol == zero) THEN
3518 IF (itol == 3) THEN
3519 l_tol=p_mach*em01
3520 ELSE
3521 l_tol=em5
3522 ENDIF
3523 ENDIF
3524 ELSE
3525 IF (ispmd == 0)THEN
3526 WRITE(iout,*) ' ** warning ** : solver non available '
3527 WRITE(IOUT,*) ' ** resetting to ** : default one '
3528 ENDIF
3529 ISOLV = 1
3530 ENDIF
3531
3532.AND. IF (IPREC > 2IPREC /= 5) THEN
3533 IF (ISPMD == 0)THEN
3534 WRITE(IOUT,*) ' ** warning ** : ',
3535 . 'this precondition method is no more supported '
3536 WRITE(IOUT,*) ' ** resetting to ** : default one '
3537 ENDIF
3538 IPREC = 5
3539 ENDIF
3540
3541 IF (N_PAT > 1) THEN
3542 IF (N_PAT > 4) THEN
3543 IF(IPRINT==1) THEN
3544 WRITE(IOUT,*) ' ** warning ** : unavailable precondition',' matrix pattern, reset to 4 ** '
3545 ENDIF
3546 N_PAT = 4
3547 ENDIF
3548 IF (IPREC /= 5) THEN
3549 IF(IPRINT==1) THEN
3550 WRITE(IOUT,*) ' ** warning ** : input precondition matrix',' pattern only available with iprec=5 : ignored **'
3551 ENDIF
3552 N_PAT = 1
3553 ENDIF
3554 ENDIF
3555
3556 P_MACH = TWO*FLMIN
3557.AND..AND. IF (NVOLU>0 IMPMV > 0 ISOLV/=1) THEN
3558 IF(IPRINT==1) THEN
3559 WRITE(IOUT,*) ' ** warning ** : direct solver is not ','compatible with monitored volume type3 '
3560 WRITE(IOUT,*) ' ** stiffness will be ignored ** '
3561 ENDIF
3562 IMPMV = 0
3563 END IF
3564
3565 IF (ILINE == 1) THEN
3566 NTY=1
3567 IF (IQSTAT > 0) NTY=6
3568 INSOLV = 0
3569 IF (ILINTF > 0) IMP_INT7 = 2
3570 IF (INTP_C < 0) ITTOFF = 1
3571 IF (IKT > 0) IKT=0
3572 NDTFIX = 0
3573 IF (NBUCK == 0) IKPRES = 0
3574 IDTC = 0
3575 IF (ISCAU > 0) THEN
3576 IF (ISMDISP > 0) THEN
3577 IF(IPRINT==1) THEN
3578 WRITE(IOUT,*) ' ** warning ** : small displacement is ignored',' incompatible with cauchy stress output option;'
3579 ENDIF
3580 ISMDISP = 0
3581 ENDIF
3582 ELSE
3583 ISMDISP = 1
3584 END IF
3585 ELSE
3586C--------nonlinear parametres-----
3587.AND. IF (ISPRB == 1IDYNA > 0) THEN
3588 IF(IPRINT==1) THEN
3589 WRITE(IOUT,*) ' ** warning ** : dynamic IMPLICIT is not ', 'compatible with impl/sprback option '
3590 WRITE(IOUT,*) ' ** resetting ** : static one '
3591 ENDIF
3592 IDYNA=0
3593 ENDIF
3594.AND. IF (ISPRB == 1ISIGINI > 0) THEN
3595 IF(IPRINT==1) THEN
3596 WRITE(IOUT,*) ' ** warning ** : pre-stresses option is not ', 'compatible with impl/sprback option '
3597 WRITE(IOUT,*) ' ** deactivating it '
3598 ENDIF
3599 ISIGINI=0
3600 ENDIF
3601
3602 IF (IQSTAT > 1) IQSTAT=1
3603
3604 IF (IDYNA > 0) THEN
3605 NTY=3
3606 IF (IQSTAT > 0) THEN
3607 IF(IPRINT==1) THEN
3608 WRITE(IOUT,*) ' ** warning ** : quasi-satic IMPLICIT is not ','compatible with dynamic option '
3609 WRITE(IOUT,*) ' ** resetting ** : dynamic one '
3610 ENDIF
3611 IQSTAT=0
3612 END IF
3613 IF (IDYNA == 2) THEN
3614 IF (NEWM_A < HALF) NEWM_A = HALF
3615 IF (NEWM_B < HALF*NEWM_A) NEWM_B = HALF*NEWM_A
3616 ELSE
3617 IF (HHT_A == ZERO) THEN
3618 HHT_A=-ZEP05
3619 ELSE
3620 HHT_A=MIN(ZERO,HHT_A)
3621 HHT_A=MAX(-THIRD,HHT_A)
3622 ENDIF
3623 ENDIF
3624 NEXP = 1
3625 ELSE
3626 NTY=2
3627 IF (IQSTAT > 0) NTY=7
3628 END IF !(IDYNA > 0)
3629C--------------Non linear solvers
3630 IF (INSOLV > 4)THEN
3631 IF(IPRINT==1) THEN
3632 WRITE(IOUT,*) ' ** warning ** : wrong nonlinear solver '
3633 WRITE(IOUT,*) ' ** resetting to** : default one '
3634 ENDIF
3635 INSOLV = 0
3636 ENDIF
3637 IF (INSOLV==4)THEN
3638 IKT=4
3639 N_LIM=1
3640.OR. IF (ISOLV==5 ISOLV==6) ISOLV=3
3641 ENDIF
3642
3643 IF (INSOLV == 0)INSOLV = 1
3644C--------stop criteria-------
3645 IF (NITOL < 12) THEN
3646 IF (NITOL > 3) THEN
3647 IF(IPRINT==1) THEN
3648 WRITE(IOUT,*) ' ** warning ** : wrong nonlinear stop ',
3649 . 'criterion number, reset to default one ** '
3650 ENDIF
3651 NITOL=0
3652 ENDIF
3653 IF (NITOL == 0) NITOL=2
3654 IF (NITOL == 1) THEN
3655 IF (N_TOL == ZERO) N_TOL=EM3
3656 ELSEIF (NITOL == 2) THEN
3657 IF (N_TOL == ZERO) N_TOL=FIVE*EM3
3658 ELSEIF (NITOL == 3) THEN
3659 IF (N_TOL == ZERO) N_TOL=EM3
3660 ENDIF
3661C IF (N_TOL == ZERO) N_TOL=EM3
3662 ELSEIF (NITOL == 12) THEN
3663 IF (N_TOLE == ZERO) N_TOLE=EM3
3664 IF (N_TOLF == ZERO) N_TOLF=EM02
3665 ELSEIF (NITOL == 23) THEN
3666 IF (N_TOLU == ZERO) N_TOLU=EM02
3667 IF (N_TOLF == ZERO) N_TOLF=EM02
3668 ELSEIF (NITOL == 13) THEN
3669 IF (N_TOLU == ZERO) N_TOLU=EM02
3670 IF (N_TOLE == ZERO) N_TOLE=EM3
3671 ELSEIF (NITOL == 123) THEN
3672 IF (N_TOLU == ZERO) N_TOLU=EM02
3673 IF (N_TOLE == ZERO) N_TOLE=EM3
3674 IF (N_TOLF == ZERO) N_TOLF=EM02
3675 ELSE
3676 CALL ANCMSG(MSGID=132,ANMODE=ANINFO,
3677 . I1=NITOL)
3678 CALL ARRET(2)
3679 END IF !(NITOL < 12) THEN
3680C
3681 IF (IDTC == 1) THEN
3682 IF (NL_DTP == 0) NL_DTP=4
3683 IF (NL_DTN == 0) NL_DTN=15
3684 IF (SCAL_DTN == ZERO) SCAL_DTN=TWO_THIRD
3685 IF (SCAL_DTP == ZERO) SCAL_DTP=ONEP01
3686 ELSEIF (IDTC == 2) THEN
3687 IF (NL_DTP == 0) NL_DTP=6
3688 IF (NL_DTN == 0) NL_DTN=20
3689 IF (IDYNA == 0) THEN
3690 IF (SCAL_DTN == ZERO) SCAL_DTN=TWO_THIRD
3691 IF (SCAL_DTP == ZERO) SCAL_DTP=ONEP1
3692 ELSE
3693 IF (SCAL_DTN == ZERO) SCAL_DTN=TWO_THIRD
3694 IF (SCAL_DTP == ZERO) SCAL_DTP=ONEP1
3695 ENDIF
3696C-------ALEN=ALEN0 IF ALEN0=0 otherwise ALEN is computed automatically
3697C-------RIKS ARC LENGTH METHOD ------
3698 ELSEIF (IDTC == 3) THEN
3699 IF (NL_DTP == 0) NL_DTP=12
3700 IF (NL_DTN == 0) NL_DTN=25
3701 IF (SCAL_DTN == ZERO) SCAL_DTN=TWO_THIRD
3702 IF (SCAL_DTP == ZERO) SCAL_DTP=ONEP2
3703 IF (IAL_M == 0) IAL_M = 2
3704 IF (NDTFIX > 0) THEN
3705 IF(IPRINT==1) THEN
3706 WRITE(ISTDO,*)
3707 . ' ** warning :riks method is not ','compatible with fixed time point '
3708 WRITE(ISTDO,*) ' ** fixed time point : deactivated '
3709 WRITE(IOUT,*)' ** warning :riks method is not ','compatible with fixed time point '
3710 WRITE(IOUT,*) ' ** fixed time point : deactivated '
3711 NDTFIX = 0
3712 IDTFIX = 0
3713 ENDIF
3714 ENDIF
3715 ELSE
3716.NOT..AND. IF((ISMDISP > 0ISOLV < 4)) THEN
3717 IF(IPRINT==1) THEN
3718 WRITE(ISTDO,*) ' ** warning: no timestep control method defined **'
3719 WRITE(ISTDO,*)' ** potential infinite loop IF no convergence is achieved **'
3720 WRITE(IOUT,*)' ** warning: no timestep control method defined **'
3721 WRITE(IOUT,*) ' ** potential infinite loop IF no convergence is achieved **'
3722 ENDIF
3723.NOT..AND. END IF !((ISMDISP > 0ISOLV < 4)) THEN
3724 IF (SCAL_DTN == ZERO) SCAL_DTN=HALF
3725 END IF !(IDTC == 1)
3726C
3727 IF (DT_MAX == ZERO) DT_MAX=EP10
3728 IF (DT_MIN == ZERO) DT_MIN=EM10
3729C
3730 IF (DT_IMP == ZERO) THEN
3731 IF(IPRINT==1) THEN
3732 WRITE(ISTDO,*)' ** warning: no initial timestep defined **'
3733 WRITE(IOUT,*)' ** warning: no initial timestep defined **'
3734 ENDIF
3735 ENDIF
3736C
3737 IF (RF_MAX == ZERO) RF_MAX=EP30
3738C
3739 IF (ISMDISP == 1 ) THEN
3740 IKG =0
3741 IKPRES = 0
3742 IF (IKPROJ == 0 ) IKPROJ =-1
3743C
3744.AND..AND..AND. IF (IDYNA > 0 IDTC > 0ISOLV < 4 SCAL_DTP /= ONE) THEN
3745 IF(IPRINT==1) THEN
3746 WRITE(ISTDO,*)' ** warning: const. time-step will be used with small disp. **'
3747 WRITE(IOUT,*)' ** warning: const. time-step will be used with small disp. **'
3748 ENDIF
3749 SCAL_DTP = ONE
3750 ENDIF
3751 END IF
3752C
3753 IF (N_LIM == 0) THEN
3754.OR. IF (ISOLV < 3IDTC == 3) THEN
3755 N_LIM=3
3756.OR. ELSEIF (ISOLV == 5ISOLV == 6) THEN
3757 N_LIM=3
3758 ELSE
3759 N_LIM=6
3760 END IF
3761 END IF
3762C
3763 END IF !IF (ILINE == 1) THEN
3764
3765 IF (IMP_CHK > 0) THEN
3766 NTY = 5
3767 IF (IQSTAT > 0) IQSTAT=1
3768 ENDIF
3769
3770 IF (NBUCK > 0) NTY = 8
3771.AND. IF (ISOLV == 3IMUMPSD == 0) IMUMPSD=1
3772C
3773 IF (ISOLV == 5) IMUMPSD=1
3774C
3775C-----------attention initialization should not be inside below--
3776 IF(IPRINT==1) THEN
3777.AND. IF (NSPMD == 1ISOLV == 3) THEN
3778 MSG_ISOL(ISOLV)='direct(bcs)'
3779 ELSEIF (ISOLV == 3) THEN
3780 MSG_ISOL(ISOLV)='direct(mumps)'
3781 ENDIF
3782 IF (NEIG == 0) THEN
3783.OR. IF(ISOLV == 3ISOLV == 4) THEN
3784 ISOL = ISOLV
3785 IF (IMUMPSV==1) ISOL = 2
3786 WRITE(IOUT,5010)MSG_TYPE(NTY),MSG_ISOL(ISOL),LPRINT
3787 ELSE
3788 ISOL=MIN(9,ISOLV)
3789 WRITE(IOUT,5000)MSG_TYPE(NTY),MSG_ISOL(ISOL),MSG_PRE(IPREC),ITOL,L_LIM,L_TOL,LPRINT
3790 ENDIF
3791 IBID =0
3792 IF (INTP_C < 0) IBID =1
3793 WRITE(IOUT,5020)IKG,IKPRES,IAUTSPC,ISPRB,IBID
3794 ELSE
3795 WRITE(IOUT,5050)MSG_TYPE(4)
3796 ENDIF
3797 IF (M_VS > 0) WRITE(IOUT,8400)M_VS,IPRO_S0
3798 IF (IKPROJ /= 0) WRITE(IOUT,7600)IKPROJ
3799C
3800 IF (ILINE /= 1) THEN
3801 IBID = INSOLV
3802 IF (INSOLV==4) INSOLV=1
3803 IF (IBID == 5) IBID = 3
3804 IF (NITOL < 12) THEN
3805 WRITE(IOUT,5100)MSG_INSOL(IBID),DT_IMP,NITOL,N_TOL
3806 ELSEIF (NITOL < 123) THEN
3807 IF (NITOL == 12) THEN
3808 WRITE(IOUT,5112)MSG_INSOL(IBID),DT_IMP,NITOL,N_TOLE,N_TOLF
3809 ELSEIF (NITOL == 23) THEN
3810 WRITE(IOUT,5123)MSG_INSOL(IBID),DT_IMP,NITOL,N_TOLF,N_TOLU
3811 ELSEIF (NITOL == 13) THEN
3812 WRITE(IOUT,5113)MSG_INSOL(IBID),DT_IMP,NITOL,N_TOLE,N_TOLU
3813 ENDIF
3814 ELSEIF (NITOL == 123) THEN
3815 WRITE(IOUT,5132)MSG_INSOL(IBID),DT_IMP,NITOL,N_TOLE,N_TOLF,
3816 . N_TOLU
3817 ENDIF
3818C
3819.OR. IF(ISOLV == 5ISOLV == 6) THEN
3820 WRITE(IOUT,5150)IPUPD,N_LIM,NPRINT,ISIGINI,IRREF,IDTC,DT_MIN,DT_MAX
3821 ELSE
3822 WRITE(IOUT,5180)N_LIM, NPRINT,ISIGINI,IRREF,IDTC,DT_MIN,DT_MAX
3823 ENDIF
3824C
3825 IF(IRIG_M == 1) THEN
3826 WRITE(IOUT,8600)
3827 ELSEIF(IRIG_M > 1) THEN
3828 WRITE(IOUT,8700) E_REF(1:3)
3829 ENDIF
3830C
3831 WRITE(IOUT,8010) ISMDISP
3832 WRITE(IOUT,8200) RF_MIN,RF_MAX
3833 IF (NCY_MAX > 0) WRITE(IOUT,7900) NCY_MAX
3834 IF (IDTC == 1) THEN
3835 WRITE(IOUT,7700)
3836 WRITE(IOUT,5200)NL_DTP,SCAL_DTP,NL_DTN,SCAL_DTN
3837 ELSEIF (IDTC == 2) THEN
3838 WRITE(IOUT,7700)
3839 WRITE(IOUT,5300)NL_DTP,NL_DTN,SCAL_DTN,SCAL_DTP,ALEN0
3840 ELSEIF (IDTC == 3) THEN
3841 S = SCAL_RIKS
3842 WRITE(IOUT,7800)
3843 WRITE(IOUT,7650)NL_DTP,NL_DTN,SCAL_DTN,SCAL_DTP,ALEN0,IAL_M,S
3844 SCAL_RIKS = S*S
3845 ENDIF
3846 END IF !(ILINE /= 1)
3847
3848
3849C
3850 IF (IDYNA == 2) THEN
3851 WRITE(IOUT,6300)NEWM_A,NEWM_B
3852 ELSEIF (IDYNA > 0) THEN
3853 WRITE(IOUT,5400)HHT_A
3854.AND. ELSEIF (IQSTAT > 0SCAL_DTQ /= ONE) THEN
3855 WRITE(IOUT,6400)SCAL_DTQ
3856 ENDIF
3857 IF (IDY_DAMP > 0) WRITE(IOUT,6200)DAMPA_IMP,DAMPB_IMP
3858 IF (N_PAT > 1) WRITE(IOUT,6100)N_PAT
3859 IF (ITTOFF > 0) WRITE(IOUT,6500)
3860 IF (ILINTF > 0) WRITE(IOUT,6600)
3861 IF (ISCAU > 0) WRITE(IOUT,6700)
3862 IF (IMP_LR > 0) WRITE(IOUT,6800)
3863C
3864 ENDIF ! IF(IPRINT==1) THEN
3865C
3866 IF (ILINE /= 1) THEN
3867 IF (ILINE_S == 0) ILINE_S=3
3868 IF (ILINE_S == 100) THEN
3869 ILINE_S=0
3870 IF(IPRINT==1)WRITE(IOUT,7000)
3871 ENDIF
3872 ELSE
3873 ILINE_S=0
3874 ENDIF
3875
3876.AND. IF (ILINE /= 1ILINE_S > 0) THEN
3877 IF (NLS_LIM == 0) NLS_LIM=20
3878 IF (LS_TOL == ZERO) LS_TOL=EM03
3879 IF(IPRINT==1)WRITE(IOUT,6900)ILINE_S,NLS_LIM,LS_TOL
3880 IF (ILINE_S == 3) THEN
3881.AND. IF (NITOL /= 2NITOL /= 4) THEN
3882 LS_TOL=FIVE*LS_TOL
3883 ELSE
3884 LS_TOL=HALF*LS_TOL
3885 ENDIF
3886 ELSEIF (ILINE_S == 2) THEN
3887 LS_TOL=HALF*LS_TOL
3888 ENDIF
3889 ENDIF
3890
3891.AND. IF (ILINE /= 1IREFI > 0) THEN
3892 IF(IPRINT==1)WRITE(IOUT,7100)IREFI
3893 END IF
3894
3895.AND. IF (ILINE /= 1NDIVER == 0) THEN
3896.AND. IF (IREFI > 2ILINE_S /= 1) THEN
3897 NDIVER=2
3898.AND. ELSEIF(IKT > 0ILINE_S /= 1) THEN
3899 NDIVER=0
3900 ELSE
3901 NDIVER=1
3902 END IF
3903 ELSEIF(NDIVER == -1) THEN
3904 NDIVER=0
3905 END IF
3906
3907.AND. IF (ILINE /= 1NDIVER > 0) THEN
3908 IF(IPRINT==1)WRITE(IOUT,7200)NDIVER
3909 NDIVER = NDIVER + 1
3910 END IF
3911
3912 IF (ILINE /= 1) THEN
3913 IF (TOL_DIV == ZERO ) THEN
3914 TOL_DIV = EP04
3915 IF (ILINE_S == 1) TOL_DIV=EP03
3916 IF (ISMDISP == 1) TOL_DIV=EP10
3917 ELSE
3918 IF(IPRINT==1)WRITE(IOUT,8500)TOL_DIV
3919 END IF
3920 END IF
3921
3922 IF (IKT == 1) THEN
3923 IF(IPRINT==1)WRITE(IOUT,7300)
3924 ELSEIF (IKT == 2) THEN
3925 IF(IPRINT==1)WRITE(IOUT,7400)
3926 ELSEIF (IKT == 3) THEN
3927 IF(IPRINT==1)WRITE(IOUT,8000)
3928 ELSEIF (IKT == 4) THEN
3929 IF(IPRINT==1)WRITE(IOUT,8100)
3930 END IF
3931
3932 IF (NDTFIX > 0) THEN
3933 IF(IPRINT==1) THEN
3934 WRITE(IOUT,7500) NDTFIX
3935 WRITE(IOUT,7510)(DTIMPF(J),J=1,NDTFIX)
3936 END IF
3937 END IF
3938
3939 IF (NBUCK > 0) THEN
3940 IF (IMUMPSV> 0) THEN
3941.AND. IF (BISOLV /= 1BISOLV /= 2) THEN
3942 WRITE(ISTDO,*) ' ** error ** unavailable solver for buckling analysis'
3943 WRITE(IOUT,*) ' ** error ** unavailable solver for buckling analysis'
3944 CALL ARRET(2)
3945 ENDIF
3946 IF (ISPMD == 0) THEN
3947 WRITE(IOUT,6000) NBUCK, SHFTBUCK, BNITER, BINCV, BMAXNCV, MSG_BSOL(2)
3948 ENDIF
3949 ELSE
3950 CALL PROUT_BUCK(IOUT,NBUCK,IBUCKL)
3951 END IF
3952 ENDIF
3953C
3954 RETURN
3955 5000 FORMAT(/
3956 . ' IMPLICIT options used :'//
3957 . ' IMPLICIT TYPE : . . . . . . . . . . . . . ',2X,A/
3958 . ' linear solver : . . . . . . . . . . . . . ',2X,A/
3959 . ' precondition method : . . . . . . . . . . ',2X,A/
3960 . ' stop criterion for linear solver . . . . .',2X,I5/
3961 . ' iteration num. limit for linear solver . .',2X,I5/
3962 . ' tolerance for linear solver . . . . . . . ',2X,G14.7/
3963 . ' printout frequency for linear solver . . .',2X,I5/)
3964 5010 FORMAT(/
3965 . ' IMPLICIT options used :'//
3966 . ' IMPLICIT TYPE : . . . . . . . . . . . . . ',2X,A/
3967 . ' linear solver : . . . . . . . . . . . . . ',2X,A/
3968 . ' printout frequency for linear solver . . .',2X,I5/)
3969 5020 FORMAT(/
3970 . ' geometrical stiffness flag . . .. . . . .',2X,I5/
3971 . ' load(pressure) stiffness flag . . . . . ',2X,I5/
3972 . ' autospc flag(0:off,1:on,2:all) . . .. . .',2X,I5/
3973 . ' spring-back option : . . . . . . . . . . .',2X,I5/
3974 . ' special pcg solver for contact . . . .. . ',2X,I5/)
3975 5050 FORMAT(/
3976 . ' IMPLICIT options used :'//
3977 . ' IMPLICIT TYPE : . . . . . . . . . . . . . ',2X,A//)
3978 5100 FORMAT(
3979 . ' non-linear solver : . . . . . . . . . . . ',2X,A/
3980 . ' initial time step . . . . . . . . . . . . ',2X,G14.7/
3981 . ' stop criterion for non-linear solver . . .',2X,I5/
3982 . ' tolerance for non-linear solver . . . . . ',2X,G14.7)
3983 5112 FORMAT(
3984 . ' non-linear solver : . . . . . . . . . . . ',2X,A/
3985 . ' initial time step . . . . . . . . . . . . ',2X,G14.7/
3986 . ' stop criterion for non-linear solver . . .',2X,I5/
3987 . ' energy tolerance for non-linear solver . .',2X,G14.7/
3988 . ' force tolerance for non-linear solver . .',2X,G14.7)
3989 5113 FORMAT(
3990 . ' non-linear solver : . . . . . . . . . . . ',2X,A/
3991 . ' initial time step . . . . . . . . . . . . ',2X,G14.7/
3992 . ' stop criterion for non-linear solver . . .',2X,I5/
3993 . ' energy tolerance for non-linear solver . .',2X,G14.7/
3994 . ' disp. tolerance for non-linear solver . .',2X,G14.7)
3995 5123 FORMAT(
3996 . ' non-linear solver : . . . . . . . . . . . ',2X,A/
3997 . ' initial time step . . . . . . . . . . . . ',2X,G14.7/
3998 . ' stop criterion for non-linear solver . . .',2X,I5/
3999 . ' force tolerance for non-linear solver . .',2X,G14.7/
4000 . ' disp. tolerance for non-linear solver . .',2X,G14.7)
4001 5132 FORMAT(
4002 . ' non-linear solver : . . . . . . . . . . . ',2X,A/
4003 . ' initial time step . . . . . . . . . . . . ',2X,G14.7/
4004 . ' stop criterion for non-linear solver . . .',2X,I5/
4005 . ' energy tolerance for non-linear solver . .',2X,G14.7/
4006 . ' force tolerance for non-linear solver . .',2X,G14.7/
4007 . ' disp. tolerance for non-linear solver . .',2X,G14.7)
4008 5150 FORMAT(
4009 . ' reforming flag in mix solver(0:auto,>0:iter)',2X,I5/
4010 . ' iteration num.limit for reforming matrix .',2X,I5/
4011 . ' printout frequency for non-linear solver .',2X,I5/
4012 . ' pre-stresses control flag . . . . . . . . ',2X,I5/
4013 . ' reference residual flag . . . . . . .. . .',2X,I5/
4014 . ' time step control method . . . . . . . . .',2X,I5/
4015 . ' minimum time step . . . . . . . . . . . . ',2X,G14.7/
4016 . ' maximum time step . . . . . . . . . . . . ',2X,G14.7//)
4017 5180 FORMAT(
4018 . ' iteration num.limit for reforming matrix .',2X,I5/
4019 . ' printout frequency for non-linear solver .',2X,I5/
4020 . ' pre-stresses control flag . . . . . . . . ',2X,I5/
4021 . ' reference residual flag . . . . . . .. . .',2X,I5/
4022 . ' time step control method . . . . . . . . .',2X,I5/
4023 . ' minimum time step . . . . . . . . . . . . ',2X,G14.7/
4024 . ' maximum time step . . . . . . . . . . . . ',2X,G14.7//)
4025 5200 FORMAT(
4026 . ' converge iteration num. for increasing dt. ',2X,I5/
4027 . ' increasing time step scale factor. . . . .',2X,G14.7/
4028 . ' converge iteration num. for decreasing dt .',2X,I5/
4029 . ' decreasing time step scale factor. . . . .',2X,G14.7/)
4030 5300 FORMAT(
4031 . ' desired converge iteration num. . . . . . .',2X,I5/
4032 . ' maximum converge iteration num. . . . . . .',2X,I5/
4033 . ' decreasing time step scale factor. . . . .',2X,G14.7/
4034 . ' maximum increasing time step scale factor ',2X,G14.7/
4035 . ' constant arc-length. . . . . . .. . . . . ',2X,G14.7//)
4036 5400 FORMAT(
4037 . ' time integration with hht-alpha constant ',2X,G14.7//)
4038 6000 FORMAT(
4039 . ' euler buckling analysis :'/
4040 . ' number of modes to be computed :. . . . . ',2X,I5/
4041 . ' shift in buckling modes pencil :. . . . . ',2X,G14.7/
4042 . ' maximum number of arnoldi iterations :. . ',2X,I5/
4043 . ' initial factor for subspace dimension : . ',2X,I5/
4044 . ' maximum factor for subspace dimension : . ',2X,I5/
4045 . ' linear solver : . . . . . . . . . . . . . ',2X,A/)
4046 6100 FORMAT(
4047 . ' precondition matrix sparse pattern(a^n) . ',2X,I5/)
4048 6200 FORMAT(
4049 . ' IMPLICIT rayleigh damping coefficients : ',2X,2G14.7/)
4050 6300 FORMAT(
4051 . ' time integration with newmark constants ',2X,2G14.7//)
4052 6400 FORMAT(
4053 . ' quasi-static inertia scale factor : ',2X,G14.7//)
4054 6500 FORMAT(
4055 . ' crossing contact node detection deactivated .',/)
4056 6600 FORMAT(
4057 . ' linear analyse taking into account contact . .',/)
4058 6700 FORMAT(
4059 . ' cauchy stress selected for linear analyse . .',/)
4060 6800 FORMAT(
4061 . ' taking into account large rigid rotation . .',/)
4062 6900 FORMAT(
4063 . ' line-search method :. . . . . . . . . . . .',2X,I5/
4064 . ' maximum line-search iteration num. . . . . ',2X,I5/
4065 . ' tolerance for line-search iteration. . . .',2X,G14.7/)
4066 7000 FORMAT(
4067 . ' line-search deactivated . . . . . . . . . . .',/)
4068 7100 FORMAT(
4069 . ' reference residual options for contact . . ',2X,I5/)
4070 7200 FORMAT(
4071 . ' divergence criterion numbers : . . . . . .',2X,I5/)
4072 7300 FORMAT(
4073 . ' simplified tangent stiffness activated . .',/)
4074 7400 FORMAT(
4075 . ' average continuum tangent stiffness activated',/)
4076 7500 FORMAT(
4077 . ' fixed time point number: . . . . . . . .',2X,I5)
4078 7510 FORMAT( /,3X,6G20.13,//)
4079 7600 FORMAT(
4080 . ' full projection for qeph(-1:off,1:on) . . .. . .',2X,I5/)
4081 7650 FORMAT(
4082 . ' desired converge iteration num. . . . . . .',2X,I5/
4083 . ' maximum converge iteration num. . . . . . .',2X,I5/
4084 . ' decreasing time step scale factor. . . . .',2X,G14.7/
4085 . ' maximum increasing time step scale factor ',2X,G14.7/
4086 . ' constant arc-length. . . . . . .. . . . . ',2X,G14.7/
4087 . ' constraint TYPE . . . .. . . . . . . . . . ',2X,I5/
4088 . ' loading contribution scale factor . . . . ',2X,G14.7//)
4089 7700 FORMAT('automatic time step control:'/)
4090 7800 FORMAT('automatic time step with riks method:'/)
4091 7900 FORMAT(
4092 . ' maximum increment(ncycle) num . . .. . .',2X,I5/)
4093 8000 FORMAT(
4094 . ' continuum tangent stiffness activated. .',/)
4095 8100 FORMAT(
4096 . ' consistent tangent stiffness activated. .',/)
4097 8200 FORMAT(
4098 . ' minimum reference force residual . . . .',2X,G14.7/
4099 . ' maximum reference force residual . . . . ',2X,G14.7//)
4100 8010 FORMAT(
4101 . ' small displacement option(0:off,1:on) . . .. . .',2X,I5/)
4102 8400 FORMAT(
4103 . ' number of projection vectors of pcg . . . ',2X,I5/
4104 . ' projection vector initialization method . ',2X,I5/)
4105 8500 FORMAT(
4106 . ' diverging tol. of relative force residual :',2X,G14.7/)
4107 8600 FORMAT(
4108 . ' free rigid motion /mrigm used . . . . . . .'/)
4109 8700 FORMAT(
4110 . ' free rigid motion /mrigm w/ ref_node_id:',2X,3I10/)
4111 END
4112
anim_build_index_all
subroutine anim_build_index_all(ispmd, mcheck, sensors, igeo, geo)
Definition anim_build_index_all.F:36
chkipari
subroutine chkipari(ipari)
Definition chkstfn3.F:211
damping
subroutine damping(nodft, nodlt, v, vr, a, ar, damp, ms, in, igrnod, dim, itask, weight, tagslv_rby, wfext)
Definition damping.F:882
prout_buck
subroutine prout_buck(ip, nbuck, ibuck)
Definition lectur.F:3262
lecimpl
subroutine lecimpl
Definition lectur.F:3321
lectur
subroutine lectur(icode, iskew, iskwn, ixtg, ixs, ixq, ixc, ixt, ixp, ixr, itab, itabm1, npc, iparg, igrv, ibgr, ipari, npby, lpby, nnlink, llink, linale, neflsw, nnflsw, icut, inoise, x, v, vr, ms, in, skew, pld, rby, wa, crflsw, xcut, dampr, igrnod, kxsp, weight, fr_rby2, fr_rl, partsav, ipart, pm, monvol, volmon, ipart_state, geo, table, iframe, xframe, elbuf_str, igeo, intbuf_tab, ipm, h3d_data, multi_fvm, igrpart, tag_skins6, icfield, lcfield, tagslv_rby, mds_label, mds_output_table, mds_nmat, max_depvar, mds_ndepsvar, stack, ibcl, iloadp, lloadp, sensors, dynain_data, dt, loads, output, names_and_titles, mat_param, glob_therm, pblast)
Definition lectur.F:105
lecfun
subroutine lecfun(npc, pld, nfct, nptst, table)
Definition lecfun.F:35
alpha
#define alpha
Definition eval.h:35
fr_rlink1
subroutine fr_rlink1(nod, itabm1, fr_rl, nsn)
Definition fr_rlink1.F:36
fr_rlale
subroutine fr_rlale(m1, m2, nod, itabm1, itag)
Definition fr_rlink1.F:131
spmd_max_s
subroutine spmd_max_s(s)
Definition imp_spmd.F:1195
lcbcsf
subroutine lcbcsf(icode, iskew, numbcsn, itab, itabm1, npby, iskwn, weight)
Definition lcbcsf.F:36
leccut
subroutine leccut(icut, xcut, itabm1)
Definition leccut.F:31
lecdamp
subroutine lecdamp(nd, dampr, igrnod)
Definition lecdamp.F:36
leceig
subroutine leceig(neigoff, neoff)
Definition leceig.F:31
lecflsw
subroutine lecflsw(nsflsw, ntflsw, neflsw, nnflsw, crflsw, x, ixs, iparg, itmp)
Definition lecflsw.F:37
lecfvbag1
subroutine lecfvbag1(nfvmodi, monvol, volmon)
Definition lecfvbag1.F:29
lecfvbag
subroutine lecfvbag(nfvmesh, monvol, volmon, x)
Definition lecfvbag.F:34
lecfxinp
subroutine lecfxinp(nfxinp)
Definition lecfxinp.F:33
lech3d
subroutine lech3d(output, geo, igeo, ipm, ipart, h3d_data, multi_fvm, ipari, iparg, tag_skins6, mds_label, mds_output_table, mds_nmat, max_depvar, mds_ndepsvar, elbuf_str, stack, ibcl, iloadp, lloadp, loads, mat_param, pblast, igrpart, npc, tf, snpc, stf)
Definition lech3d.F:97
lecinv
subroutine lecinv(niniv, x, v, vr, itab, iframe, xframe, igrnod, fxbipm, fxbvit, fxbrpm)
Definition lecinv.F:39
lecnoise
subroutine lecnoise(inoise, itabm1, names_and_titles)
Definition lecnoise.F:45
lecstat
subroutine lecstat(ipart, ipart_state, elbuf_tab, ipm, iparg, iparts, ipartc, ipartg, mat_param)
Definition lecstat.F:37
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
for
for(i8=*sizetab-1;i8 >=0;i8--)
Definition mumps_common.c:91
ale_mod
Definition ale_mod.F:150
ale_mod::ale
type(ale_) ale
Definition ale_mod.F:253
alemuscl_mod
Definition alemuscl_mod.F:55
alemuscl_mod::alemuscl_param
type(alemuscl_param_) alemuscl_param
Definition alemuscl_mod.F:76
check_mod
Definition check_mod.F:34
check_mod::restart_file
integer restart_file
Definition check_mod.F:52
check_mod::check_message
character(len=2048), dimension(check_message_size) check_message
Definition check_mod.F:54
dt_mod
Definition dt_mod.F:80
fvbag_mod
Definition fvbag_mod.F:44
fvbag_mod::fvdata_old
type(fvbag_data), dimension(:), allocatable fvdata_old
Definition fvbag_mod.F:193
fvbag_mod::num_opt_dt_fvmbag_1
integer num_opt_dt_fvmbag_1
Definition fvbag_mod.F:199
fvbag_mod::num_opt_dt_fvmbag_0
integer num_opt_dt_fvmbag_0
Definition fvbag_mod.F:198
fxb_mod
Definition fxb_mod.F:32
fxb_mod::fxbipm
integer, dimension(:,:), allocatable fxbipm
Definition fxb_mod.F:39
groupdef_mod
Definition groupdef_mod.F:668
h3d_mod
Definition h3d_mod.F:311
imp_kbcs
Definition impbufdef_mod.F:471
imp_kbcs::maxset_b
integer maxset_b
Definition impbufdef_mod.F:475
imp_kbcs::msgl_b
integer msgl_b
Definition impbufdef_mod.F:475
imp_pcg_proj
Definition impbufdef_mod.F:646
imp_pcg_proj::m_vs
integer m_vs
Definition impbufdef_mod.F:648
imp_spbrm
Definition impbufdef_mod.F:696
intstamp_glob_mod
Definition intstamp_glob_mod.F:30
message_mod
Definition message_mod.F:1257
names_and_titles_mod
Definition names_and_titles_mod.F:1002
stack_mod
Definition stack_mod.F:75
state_mod
Definition state_mod.F:54
table_mod
Definition table_mod.F:113
floatmin
void floatmin(int *a, int *b, float *flm)
Definition precision.c:71
rbyonf
subroutine rbyonf(iparg, ipari, ms, in, ixs, ixq, ixc, ixt, ixp, ixr, skew, itab, itabm1, iskwn, npby, onof, nrbynf, itag, lpby, rby, x, v, vr, ixtg, igrv, ibgr, weight, fr_rby2, partsav, ipart, elbuf_tab, icfield, lcfield, tagslv_rby)
Definition rbyonf.F:42
read_dynain
subroutine read_dynain(ipart, dynain_data, ipartc, iparttg, ixc, ixtg)
Definition read_dynain.F:39
rlink0
subroutine rlink0(v, vr, ms, in, nlik, idir, idrot, nod, nsn, ic, icr, iskw, isk, skew, iskwn, fr_rl, weight)
Definition rlink0.F:35
spmd_chkw
subroutine spmd_chkw(iwiout, iout)
Definition spmd_chkw.F:38
spmd_ibcast
subroutine spmd_ibcast(tabi, tabr, n1, n2, from, add)
Definition spmd_ibcast.F:57
spmd_glob_isum9
subroutine spmd_glob_isum9(v, len)
Definition spmd_th.F:520
spmd_wiout
subroutine spmd_wiout(iout, iwiout)
Definition spmd_wiout.F:40
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:895
arret
subroutine arret(nn)
Definition arret.F:86
static
subroutine static(v, vr, a, ar, ms, in, igrnod, weight_md, wfext)
Definition static.F:33
table_mod::ttable
Definition table_mod.F:126