OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
rdresb.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!|| rdresb ../engine/source/output/restart/rdresb.F
25!||--- called by ------------------------------------------------------
26!|| radioss2 ../engine/source/engine/radioss2.F
27!||--- calls -----------------------------------------------------
28!|| alelag_rrest ../engine/source/output/restart/rdresb.F
29!|| ancmsg ../engine/source/output/message/message.F
30!|| arret ../engine/source/system/arret.F
31!|| cfieldrest ../engine/source/output/restart/rdresb.F
32!|| checksum_restart_read ../common_source/modules/output/checksum_mod.F90
33!|| convrrest ../engine/source/output/restart/rdresb.F
34!|| crkxfem_rrest ../engine/source/output/restart/rdresb.F
35!|| crkxfem_rrestanim ../engine/source/output/restart/rdresb.F
36!|| cur_fil_c ../common_source/tools/input_output/write_routines.c
37!|| decompress_i_nnz ../common_source/comm/decompress_nnz.F
38!|| decompress_r_nnz ../common_source/comm/decompress_nnz.F
39!|| drape_rrest ../engine/source/output/restart/rdresb.F
40!|| eigrrest ../engine/source/output/restart/rdresb.F
41!|| elbuf_ini ../engine/source/elements/elbuf/elbuf_ini.F
42!|| fillxdp ../engine/source/output/restart/fillxdp.F
43!|| fvrrest ../engine/source/output/restart/rdresb.F
44!|| fxbresb ../engine/source/output/restart/rdresb.F
45!|| fxfluxrrest ../engine/source/constraints/thermic/fxfluxrrest.F
46!|| fxtemprrest ../engine/source/output/restart/rdresb.F
47!|| get_nrbe3pen_l ../engine/source/output/restart/restart_rbe3pen.F90
48!|| get_nrwallpen_l ../engine/source/output/restart/restart_rwallpen.F90
49!|| group_ini ../engine/source/groups/group_ini.F
50!|| imprrest ../engine/source/output/restart/rdresb.F
51!|| init_global_node_id ../common_source/modules/nodal_arrays.F90
52!|| init_global_shell_id ../common_source/modules/connectivity.F90
53!|| intbuf_ini ../common_source/interf/intbuf_ini.F
54!|| intfric_rresti ../common_source/modules/interfaces/intbuf_fric_mod.F90
55!|| intfric_rrestr ../common_source/modules/interfaces/intbuf_fric_mod.F90
56!|| intstamp_rresti ../engine/share/modules/intstamp_mod.F
57!|| intstamp_rrestr ../engine/share/modules/intstamp_mod.F
58!|| knotlocelrest ../engine/source/output/restart/rdresb.F
59!|| knotlocpcrest ../engine/source/output/restart/rdresb.F
60!|| knotrest ../engine/source/output/restart/rdresb.F
61!|| line_ini ../engine/source/groups/line_ini.F
62!|| loadprest ../engine/source/output/restart/rdresb.F
63!|| nfrrest ../engine/source/output/restart/rdresb.F
64!|| nitscherrest ../engine/source/output/restart/rdresb.F
65!|| pblast_init_tables ../common_source/modules/loads/pblast_mod.F90
66!|| pblast_load ../common_source/modules/loads/pblast_mod.F90
67!|| ply_info_rest ../engine/source/output/restart/rdresb.F
68!|| plyxfem_ravuply ../engine/source/output/restart/rdresb.F
69!|| plyxfem_rrest ../engine/source/output/restart/rdresb.F
70!|| plyxfem_rrestanim ../engine/source/output/restart/rdresb.F
71!|| python_deserialize ../common_source/modules/python_mod.F90
72!|| r_bufbric_22 ../engine/source/interfaces/int22/r_bufbric_22.F
73!|| radiarrest ../engine/source/output/restart/rdresb.F
74!|| rdcomr ../engine/source/output/restart/rdcomm.F
75!|| read_ale_grid ../engine/source/output/restart/read_ale_grid.F90
76!|| read_bcs_nrf ../engine/source/output/restart/read_bcs_nrf.F90
77!|| read_bcs_wall ../engine/source/output/restart/read_bcs_wall.F90
78!|| read_cluster ../engine/source/output/cluster/read_cluster.F
79!|| read_db ../common_source/tools/input_output/read_db.F
80!|| read_dpdb ../common_source/tools/input_output/read_db.F
81!|| read_elgroup_param ../engine/source/output/restart/read_elgroup_param.F
82!|| read_failwave ../engine/source/output/restart/read_failwave.F
83!|| read_i_c ../common_source/tools/input_output/write_routines.c
84!|| read_inivel ../engine/source/output/restart/read_inivel.F90
85!|| read_joint ../engine/source/output/restart/read_joint.F
86!|| read_matparam ../engine/source/output/restart/read_matparam.F
87!|| read_nloc_struct ../engine/source/output/restart/read_nloc_struct.F
88!|| read_pcyl ../engine/source/output/restart/read_pcyl.F
89!|| read_rrbe3pen ../engine/source/output/restart/restart_rbe3pen.F90
90!|| read_rrwallpen ../engine/source/output/restart/restart_rwallpen.F90
91!|| read_sensors ../engine/source/output/restart/read_sensors.F
92!|| read_th_restart ../engine/source/output/th/read_th_restart.F
93!|| read_units ../common_source/comm/write_units.F
94!|| rmatrrest ../engine/source/output/restart/rdresb.F
95!|| rthbufrest ../engine/source/output/restart/rdresb.F
96!|| spmd_inisph ../engine/source/mpi/sph/spmd_sptool.F
97!|| spmd_initfi ../engine/source/mpi/interfaces/spmd_i7tool.f
98!|| stack_rrest ../engine/source/output/restart/rdresb.F
99!|| subset_ini ../engine/source/groups/subset_ini.F
100!|| surf_ini ../engine/source/groups/surf_ini.F
101!|| table_rresti ../engine/source/tools/curve/table_tools.F
102!|| table_rrestr ../engine/source/tools/curve/table_tools.F
103!|| thcrrest ../engine/source/output/restart/rdresb.F
104!|| thsurf_read_restart ../common_source/modules/interfaces/th_surf_mod.F
105!|| userwi_read ../engine/source/user_interface/userwindow_interface_routines.f
106!|| wigerest ../engine/source/output/restart/rdresb.F
107!||--- uses -----------------------------------------------------
108!|| ale_connectivity_mod ../common_source/modules/ale/ale_connectivity_mod.F
109!|| ale_mod ../common_source/modules/ale/ale_mod.F
110!|| alefvm_mod ../common_source/modules/ale/alefvm_mod.F
111!|| bcs_mod ../common_source/modules/boundary_conditions/bcs_mod.F90
112!|| checksum_output_option_mod ../common_source/modules/output/checksum_mod.F90
113!|| cluster_mod ../engine/share/modules/cluster_mod.F
114!|| connectivity_mod ../common_source/modules/connectivity.F90
115!|| damp_mod ../engine/share/modules/damp_mod.F
116!|| drape_mod ../engine/share/modules/drape_mod.f
117!|| ebcs_mod ../common_source/modules/boundary_conditions/ebcs_mod.F90
118!|| ecnd_mod ../engine/share/modules/ecdn_mod.F
119!|| element_mod ../common_source/modules/elements/element_mod.F90
120!|| glob_therm_mod ../common_source/modules/mat_elem/glob_therm_mod.f90
121!|| h3d_mod ../engine/share/modules/h3d_mod.F
122!|| interfaces_mod ../common_source/modules/interfaces/interfaces_mod.f90
123!|| intstamp_glob_mod ../engine/share/modules/intstamp_glob_mod.F
124!|| loads_mod ../common_source/modules/loads/loads_mod.F90
125!|| mat_elem_mod ../common_source/modules/mat_elem/mat_elem_mod.F90
126!|| message_mod ../engine/share/message_module/message_mod.F
127!|| monvol_struct_mod ../engine/share/modules/monvol_struct_mod.F
128!|| multi_fvm_mod ../common_source/modules/ale/multi_fvm_mod.F90
129!|| nlocal_reg_mod ../common_source/modules/nlocal_reg_mod.f
130!|| nodal_arrays_mod ../common_source/modules/nodal_arrays.F90
131!|| outmax_mod ../common_source/modules/outmax_mod.F
132!|| output_mod ../common_source/modules/output/output_mod.F90
133!|| pblast_mod ../common_source/modules/loads/pblast_mod.F90
134!|| pinchtype_mod ../common_source/modules/pinchtype_mod.F
135!|| pload_cyl_mod ../common_source/modules/loads/pload_cyl_mod.F
136!|| python_funct_mod ../common_source/modules/python_mod.F90
137!|| rbe3_mod ../common_source/modules/constraints/rbe3_mod.F90
138!|| read_ale_grid_mod ../engine/source/output/restart/read_ale_grid.F90
139!|| read_bcs_nrf_mod ../engine/source/output/restart/read_bcs_nrf.F90
140!|| read_bcs_wall_mod ../engine/source/output/restart/read_bcs_wall.F90
141!|| read_inivel_mod ../engine/source/output/restart/read_inivel.F90
142!|| restart_rbe3pen_mod ../engine/source/output/restart/restart_rbe3pen.F90
143!|| restart_rwallpen_mod ../engine/source/output/restart/restart_rwallpen.f90
144!|| restmod ../engine/share/modules/restart_mod.F
145!|| rwall_mod ../common_source/modules/constraints/rwall_mod.F90
146!|| seatbelt_mod ../common_source/modules/seatbelt_mod.F
147!|| sensor_mod ../common_source/modules/sensor_mod.F90
148!|| skew_mod ../common_source/modules/skew_mod.F90
149!|| sph_mod ../engine/share/modules/sph_mod.F
150!|| stack_mod ../engine/share/modules/stack_mod.F
151!|| state_mod ../common_source/modules/state_mod.F
152!|| table_glob_mod ../engine/share/modules/table_glob_mod.f
153!|| time_mod ../engine/share/modules/time_mod.F
154!|| unitab_mod ../common_source/modules/unitab_mod.F
155!|| user_windows_mod ../common_source/modules/user_windows_mod.F
156!||====================================================================
157 SUBROUTINE rdresb( AF ,IAF ,LVARREA ,IRUNN ,
158 . MULTI_FVM ,H3D_DATA ,PINCH_DATA ,ALE_CONNECTIVITY ,T_MONVOL ,
159 . SENSORS ,EBCS_TAB ,DYNAIN_DATA ,USER_WINDOWS ,OUTPUT ,
160 . INTERFACES ,LOADS ,MAT_ELEM ,PYTHON ,IFLOW ,
161 . SKEWS ,RFLOW ,LIFLOW ,LRFLOW ,IMPL_S0 ,
162 . FORNEQS ,UNITAB ,
163 . STACK ,DRAPE_SH4N ,DRAPE_SH3N ,DRAPEG ,NDRAPE ,
164 . GLOB_THERM ,PBLAST ,ELEMENT ,NODES ,RBE3 ,
165 . RWALL )
166C-----------------------------------------------
167C M o d u l e s
168C-----------------------------------------------
169 USE connectivity_mod
170 USE nodal_arrays_mod
171 USE python_funct_mod
172 USE restmod
173 USE intstamp_glob_mod
174 USE table_glob_mod
175 USE message_mod
176 USE sph_mod
177 USE cluster_mod
178 USE ecnd_mod
179 USE multi_fvm_mod
180 USE h3d_mod
181 USE nlocal_reg_mod
182 USE mat_elem_mod
183 USE pinchtype_mod
184 USE ale_connectivity_mod
185 USE monvol_struct_mod
186 USE sensor_mod
187 USE ebcs_mod
188 USE seatbelt_mod
189 USE pblast_mod
190 USE outmax_mod
191 USE pload_cyl_mod
192 USE state_mod
193 USE user_windows_mod
194 USE ale_mod
195 USE output_mod
196 USE interfaces_mod
197 USE sensor_mod
198 USE loads_mod
199 USE damp_mod
200 USE bcs_mod , ONLY : bcs
201 USE skew_mod
202 USE unitab_mod
203 USE alefvm_mod
204 USE stack_mod
205 USE drape_mod
206 USE time_mod
207 use read_inivel_mod, only : read_inivel
208 use glob_therm_mod
209 use read_ale_grid_mod, only : read_ale_grid
210 use read_bcs_wall_mod, only : read_bcs_wall
211 use read_bcs_nrf_mod, only : read_bcs_nrf
212 use rbe3_mod
213 use restart_rbe3pen_mod, only : get_nrbe3pen_l,read_rrbe3pen
214 use checksum_output_option_mod
215 use element_mod , only : nixc
216 use rwall_mod
217 use restart_rwallpen_mod, only : get_nrwallpen_l,read_rrwallpen
218C-----------------------------------------------
219C I m p l i c i t T y p e s
220C-----------------------------------------------
221#include "implicit_f.inc"
222C-----------------------------------------------
223C G l o b a l P a r a m e t e r s
224C-----------------------------------------------
225#include "mvsiz_p.inc"
226C-----------------------------------------------
227C C o m m o n B l o c k s
228C-----------------------------------------------
229#include "com01_c.inc"
230#include "com04_c.inc"
231#include "com_xfem1.inc"
232#include "couple_c.inc"
233#include "intstamp_c.inc"
234#include "param_c.inc"
235#include "parit_c.inc"
236#include "remesh_c.inc"
237#include "scr06_c.inc"
238#include "scr03_c.inc"
239#include "scr05_c.inc"
240#include "scr07_c.inc"
241#include "scr14_c.inc"
242#include "scr16_c.inc"
243#include "scr19_c.inc"
244#include "scrcut_c.inc"
245#include "scrfs_c.inc"
246#include "scrnoi_c.inc"
247#include "scr_fac_c.inc"
248#include "sms_c.inc"
249#include "sphcom.inc"
250#include "spmd_c.inc"
251#include "tabsiz_c.inc"
252#include "task_c.inc"
253C-----------------------------------------------
254C D u m m y A r g u m e n t s
255C-----------------------------------------------
256 INTEGER IAF(*),LVARREA,LEN_G,LEN_M, IRUNN
257 my_real AF(*)
258 TYPE(MULTI_FVM_STRUCT) :: MULTI_FVM
259 TYPE(H3D_DATABASE) :: H3D_DATA
260 TYPE(PINCH) :: PINCH_DATA
261 TYPE(t_ale_connectivity), INTENT(INOUT) :: ALE_CONNECTIVITY
262 TYPE(MONVOL_STRUCT_),DIMENSION(NVOLU), INTENT(INOUT) :: T_MONVOL
263 TYPE(t_ebcs_tab), INTENT(INOUT) :: EBCS_TAB
264 TYPE(dynain_database) ,INTENT(INOUT) :: DYNAIN_DATA
265 TYPE(user_windows_) ,INTENT(INOUT) :: USER_WINDOWS
266 TYPE(output_) ,INTENT(INOUT) :: OUTPUT !< output structure
267 TYPE(interfaces_) ,INTENT(INOUT) :: INTERFACES
268 TYPE(sensors_) ,INTENT(INOUT) :: SENSORS
269 TYPE(loads_) ,INTENT(INOUT) :: LOADS
270 TYPE(mat_elem_) ,INTENT(INOUT) :: MAT_ELEM
271 TYPE(python_) ,INTENT(INOUT) :: PYTHON
272 TYPE(skew_) ,INTENT(INOUT) :: SKEWS
273 TYPE(unit_type_) ,INTENT(INOUT) :: UNITAB !< units conversion
274 TYPE (STACK_PLY) ,INTENT(INOUT) :: STACK
275 TYPE(drape_) ,INTENT(INOUT) :: DRAPE_SH3N(NUMELTG_DRAPE),DRAPE_SH4N(NUMELC_DRAPE)
276 TYPE(drapeg_) ,INTENT(INOUT) :: DRAPEG
277 TYPE(glob_therm_) ,INTENT(INOUT) :: GLOB_THERM
278 TYPE(pblast_) ,INTENT(INOUT) :: PBLAST
279 TYPE(connectivity_) ,INTENT(INOUT) :: ELEMENT
280 TYPE(nodal_arrays_) ,INTENT(INOUT) :: NODES
281 TYPE(rbe3_) ,INTENT(INOUT) :: RBE3
282 TYPE(rwall_) ,INTENT(INOUT) :: RWALL
283 INTEGER, INTENT(IN) :: LIFLOW
284 INTEGER, INTENT(IN) :: LRFLOW
285 my_real, INTENT(INOUT) :: forneqs(3,numnod)
286 INTEGER, INTENT(INOUT) :: IFLOW(LIFLOW)
287 my_real, INTENT(INOUT) :: rflow(lrflow)
288 INTEGER, INTENT(IN) :: IMPL_S0
289 INTEGER, INTENT(IN) :: NDRAPE
290C-----------------------------------------------
291C L o c a l V a r i a b l e s
292C-----------------------------------------------
293 INTEGER IDUMM(MVSIZ)
294 INTEGER I,J,LMXVREA, LEN,MY_ILEN,MY_RLEN,ISTAT,II,NS_DIFF
295 parameter(lmxvrea=1000+30*maxlaw+30)
296 my_real tabvrea(lmxvrea)
297 INTEGER, DIMENSION(:), ALLOCATABLE::IBID
298 my_real, DIMENSION(:), ALLOCATABLE::RBID
299 DOUBLE PRECISION, DIMENSION(:),ALLOCATABLE :: XDPDUM
300 INTEGER, DIMENSION(NVOLU) :: NTRI
301 INTEGER, DIMENSION(:), ALLOCATABLE :: INTEGER_DATA
302 INTEGER :: INT1, INT2, INT3, INT0, ITMP
303 INTEGER :: MY_SIZE
304 INTEGER :: PYTHON_LEN
305 INTEGER, DIMENSION(:), ALLOCATABLE :: BUFFER
306 INTEGER :: SKEW_LEN
307 LOGICAL :: IS_EBCS_PARALLEL
308 INTEGER :: NPT,NRBE3PEN_L
309C-----------------------------------------------
310C S o u r c e L i n e s
311C-----------------------------------------------
312 DO i=1,lmxvrea
313 tabvrea(i)=zero
314 ENDDO
315 IF(lvarrea>lmxvrea)THEN
316 CALL ancmsg(msgid=33,anmode=aninfo,i1=lvarrea,i2=lmxvrea)
317 CALL arret(2)
318 ENDIF
319C--------------------------------------
320C READING INTEGERS
321C--------------------------------------
322C Restart file in binary format ieee
323 len_g = npropgi*numgeo
324 len_m = npropmi*nummat
325C
326 CALL cur_fil_c(0)
327
328 CALL checksum_restart_read(output%CHECKSUM)
329
330 CALL read_i_c(nodes%NODGLOB,snodglob)
331
332 CALL read_i_c(nodes%MAIN_PROC,sweight)
333
334 CALL read_i_c(nodes%WEIGHT,sweight)
335
336 CALL read_i_c(nodes%ICODE,sicode)
337
338 CALL read_i_c(nodes%ISKEW,siskew)
339
340 CALL read_i_c(iskwn,siskwn)
341
342 CALL read_i_c(iframe,siframe)
343
344 CALL read_i_c(neth,sneth)
345
346 CALL read_i_c(ibcslag,sibcslag)
347
348 CALL read_i_c(ipart,sipart)
349
350 CALL read_i_c(nom_opt,snom_opt)
351
352 CALL read_i_c(python_len,1)
353
354 ! read python data in form of integer buffer
355 ALLOCATE(buffer(python_len))
356 buffer(1) = python_len
357 CALL read_i_c(buffer(2:python_len), python_len-1)
358 CALL python_deserialize(python, buffer) ! fills python structure
359 DEALLOCATE(buffer)
360
361 CALL read_i_c(npc,snpc)
362
363 CALL read_i_c(ixtg,sixtg)
364
365 CALL read_i_c(ixtg1,sixtg1)
366
367 CALL read_i_c(ixs,sixs)
368
369 CALL read_i_c(ixq,sixq)
370
371 CALL read_i_c(element%shell%ixc,sixc)
372 DO i = 1, sixc / nixc ! loop over the number of shell
373 element%shell%nodes(1:4,i) = element%shell%ixc(2:5, i)
374 element%shell%pid(i) = element%shell%ixc(6, i)
375 element%shell%matid(i) = element%shell%ixc(1, i)
376 element%shell%user_id(i) = element%shell%ixc(7, i)
377 ENDDO
378 CALL init_global_shell_id(element%SHELL)
379
380 CALL read_i_c(ixt,sixt)
381
382 CALL read_i_c(ixp,sixp)
383
384 CALL read_i_c(ixr,sixr)
385
386 CALL read_i_c(nodes%ITAB,sitab)
387 ! Create the map from global to local ids
388 CALL init_global_node_id(nodes,numnod)
389
390 CALL read_i_c(nodes%ITABM1,sitabm1)
391
392 CALL read_i_c(gjbufi,sgjbufi)
393
394 CALL read_i_c(ale_connectivity%NALE,ale%GLOBAL%SNALE)
395
396 IF (iale + ieuler + ialelag + glob_therm%ITHERM /= 0) THEN
397 CALL read_i_c(itmp, 1)
398 ALLOCATE(ale_connectivity%ee_connect%iad_connect(itmp))
399 CALL read_i_c(ale_connectivity%ee_connect%iad_connect, itmp)
400
401 CALL read_i_c(itmp, 1)
402 ALLOCATE(ale_connectivity%ee_connect%connected(itmp))
403 CALL read_i_c(ale_connectivity%ee_connect%connected, itmp)
404
405 CALL read_i_c(itmp, 1)
406 ALLOCATE(ale_connectivity%ee_connect%type(itmp))
407 CALL read_i_c(ale_connectivity%ee_connect%type, itmp)
408
409 CALL read_i_c(itmp, 1)
410 ALLOCATE(ale_connectivity%ee_connect%iface2(itmp))
411 CALL read_i_c(ale_connectivity%ee_connect%iface2, itmp)
412 ENDIF
413
414 CALL read_i_c(id_global_vois,ale%GLOBAL%SIELVS)
415
416 CALL read_i_c(face_vois,ale%GLOBAL%SIELVS)
417
418 IF(ale_connectivity%has_ne_connect) THEN
419 CALL read_i_c(ne_nercvois, nspmd + 1)
420 ALLOCATE(ne_lercvois(ne_nercvois(nspmd+1)))
421 CALL read_i_c(ne_lercvois, ne_nercvois(nspmd+1))
422 CALL read_i_c(ne_nesdvois, nspmd + 1)
423 ALLOCATE(ne_lesdvois(ne_nesdvois(nspmd+1)))
424 CALL read_i_c(ne_lesdvois, ne_nesdvois(nspmd+1))
425 ELSE
426 ALLOCATE(ne_lercvois(0))
427 ALLOCATE(ne_lesdvois(0))
428 ENDIF
429 IF (iale > 0 .AND. ale%GRID%NWALE_RST >= 0) THEN
430 IF(ale_connectivity%has_nn_connect)THEN
431 ALLOCATE(ale_connectivity%NN_CONNECT%IAD_CONNECT(numnod + 1))
432 CALL read_i_c(ale_connectivity%NN_CONNECT%IAD_CONNECT, numnod + 1)
433 ALLOCATE(ale_connectivity%NN_CONNECT%CONNECTED(ale_connectivity%NN_CONNECT%IAD_CONNECT(numnod + 1)))
434 CALL read_i_c(ale_connectivity%NN_CONNECT%CONNECTED, ale_connectivity%NN_CONNECT%IAD_CONNECT(numnod + 1))
435 ENDIF
436 ENDIF
437 IF(ale_connectivity%has_ne_connect)THEN
438 ALLOCATE(ale_connectivity%NE_CONNECT%IAD_CONNECT(numnod + 1))
439 CALL read_i_c(ale_connectivity%NE_CONNECT%IAD_CONNECT, numnod + 1)
440 ALLOCATE(ale_connectivity%NE_CONNECT%CONNECTED(ale_connectivity%NE_CONNECT%IAD_CONNECT(numnod + 1)))
441 CALL read_i_c(ale_connectivity%NE_CONNECT%CONNECTED, ale_connectivity%NE_CONNECT%IAD_CONNECT(numnod + 1))
442 ENDIF
443
444 CALL read_i_c(ns_diff, 1)
445 multi_fvm%NS_DIFF = .false.
446 IF (ns_diff == 1) THEN
447 multi_fvm%NS_DIFF = .true.
448 ale_connectivity%HAS_IDGLOB = .true.
449 IF (.NOT. ALLOCATED(ale_connectivity%IDGLOB%ID)) THEN
450 ALLOCATE(ale_connectivity%IDGLOB%ID(numels + nsvois + numelq + nqvois + numeltg + ntgvois))
451 ENDIF
452 IF (.NOT. ALLOCATED(ale_connectivity%IDGLOB%UID)) THEN
453 ALLOCATE(ale_connectivity%IDGLOB%UID(numels + nsvois + numelq + nqvois + numeltg + ntgvois))
454 ENDIF
455 IF (.NOT. ALLOCATED(multi_fvm%ELEM_DATA%CENTROID)) THEN
456 ALLOCATE(multi_fvm%ELEM_DATA%CENTROID(3, numels + nsvois + numelq + nqvois + numeltg + ntgvois))
457 ENDIF
458 IF (.NOT. ALLOCATED(multi_fvm%FACE_DATA%CENTROID)) THEN
459 ALLOCATE(multi_fvm%FACE_DATA%CENTROID(3, 6, numels + nsvois + numelq + nqvois + numeltg + ntgvois))
460 ENDIF
461 IF (.NOT. ALLOCATED(multi_fvm%VOL)) THEN
462 ALLOCATE(multi_fvm%VOL(numels + nsvois + numelq + nqvois + numeltg + ntgvois))
463 ENDIF
464 CALL read_i_c(ale_connectivity%IDGLOB%ID, numels + nsvois + numelq + nqvois + numeltg + ntgvois)
465 CALL read_i_c(ale_connectivity%IDGLOB%UID, numels + nsvois + numelq + nqvois + numeltg + ntgvois)
466 ENDIF
467
468 CALL read_i_c(ifill,sifill)
469
470 CALL read_i_c(ims,sims)
471
472 CALL subset_ini
473
474 CALL group_ini
475
476 CALL surf_ini
477
478 CALL line_ini
479
480 CALL read_i_c(kxx,skxx)
481
482 CALL read_i_c(ixx,sixx)
483
484 CALL read_i_c(kxsp,skxsp)
485
486 CALL read_i_c(ixsp,sixsp)
487
488 CALL read_i_c(nod2sp,snod2sp)
489
490 CALL read_i_c(ispsym,sispsym)
491
492 CALL read_i_c(ispcond,sispcond)
493
494 CALL read_i_c(isphio,sisphio)
495
496 CALL read_i_c(lprtsph,slprtsph)
497
498 CALL read_i_c(lonfsph,slonfsph)
499
500 CALL read_i_c(ibufssg_io,sibufssg_io)
501
502 IF(nsphsol/=0)THEN
503
504 CALL read_i_c(sph2sol,numsph)
505
506 CALL read_i_c(sol2sph,2*numels8)
507
508 CALL read_i_c(irst,3*nsphsol)
509
510 CALL read_i_c(sol2sph_typ,numels8)
511 END IF
512
513 CALL read_i_c(ibcl,sibcl)
514
515 IF (nintloadp>0) THEN
516 CALL read_i_c(kloadpinter, ninter+1)
517 CALL read_i_c(loadpinter,ninter*nloadp_hyd)
518 ENDIF
519
520 CALL read_i_c(ibfv,sibfv)
521
522 CALL read_i_c(iecran,siecran)
523
524 CALL read_i_c(ilas,silas)
525
526 ! -----------------------------
527 ! EBCS option
528 CALL read_i_c(int0, 1) ! total number of EBCS
529 CALL read_i_c(int1, 1) ! number of multifluid ebcs
530 CALL read_i_c(int2, 1) ! number of parallel ebcs
531 CALL read_i_c(int3, 1) ! number of other ebcs
532
533 ebcs_tab%nebcs = int0 !(INT0=INT1+INT2+INT3)
534 ebcs_tab%nebcs_fvm = int1
535 ebcs_tab%nebcs_parallel = int2
536 ebcs_tab%nebcs_loc = int3
537
538 IF (int0 > 0) THEN
539 ALLOCATE(integer_data(int0))
540 CALL ebcs_tab%read_type_data(int0, integer_data)
541 ALLOCATE(ebcs_tab%my_typ(int0))
542 ALLOCATE(ebcs_tab%need_to_compute(int0))
543 ! save the type of /EBCS ( /EBCS/NRF --> typ = 10) ( /EBCS/PROPELLANT --> typ = 11)
544 ebcs_tab%my_typ(1:int0) = integer_data(1:int0)
545 ! -------------
546 ! check if the current proc needs to compute an /EBCS
547 ! multifluid /EBCS, /EBCS/NRF, and /EBCS/PROPELLANT are parallelized, can be computed by several proc
548 ! the others are sequential, only the main proc (ispmd=0) computes the other /EBCS
549 DO ii=1,int0
550 ebcs_tab%need_to_compute(ii) = .false.
551 IF(ispmd /= 0) THEN
552 is_ebcs_parallel = .false.
553 IF(ebcs_tab%my_typ(ii) == 10 .or. ebcs_tab%my_typ(ii) == 11)is_ebcs_parallel = .true.
554 IF(is_ebcs_parallel) ebcs_tab%need_to_compute(ii) = .true. ! /EBCS/NRF and /EBCS/PROPELLANT are //
555 ELSE
556 ebcs_tab%need_to_compute(ii) = .true. ! other /EBCS sequential
557 ENDIF
558 ENDDO
559 ! -------------
560 ! create the data of /EBCS options
561 CALL ebcs_tab%create_from_types(int1 + int2 + int3 , integer_data(1:int1 + int2 + int3 ))
562 ENDIF
563 ! read the data of /EBCS options
564 IF (int1 + int2 + int3 > 0) THEN
565 DO ii = 1, int1 + int2 + int3
566 CALL ebcs_tab%tab(ii)%poly%read_common_data()
567 CALL ebcs_tab%tab(ii)%poly%read_data()
568 ENDDO
569 ENDIF
570 ! array allocation : adress for parith/on option
571 ALLOCATE( ebcs_parithon(ebcs_tab%nebcs) )
572 ! -----------------------------
573
574 CALL read_i_c(laccelm,slaccelm)
575
576 CALL read_i_c(lgauge,3*nbgauge)
577
578 CALL read_i_c(nnlink,snnlink)
579
580 CALL read_i_c(lnlink,slnlink)
581
582 CALL userwi_read(user_windows,ispmd,nspmd,numnod)
583
584 CALL read_i_c(iparg,siparg)
585
586 CALL read_i_c(dd_iad,sdd_iad)
587
588 CALL read_i_c(igrv,sigrv)
589
590 CALL read_i_c(lgrav,slgrav)
591
592 CALL read_i_c(iexlnk,siexlnk)
593
594 CALL read_i_c(ibvel,sibvel)
595
596 CALL read_i_c(lbvel,slbvel)
597
598 CALL read_i_c(irbe2,sirbe2)
599
600 CALL read_i_c(lrbe2,slrbe2)
601
602 CALL read_i_c(rbe3%IRBE3,rbe3%nrbe3*irbe3_variables)
603
604 CALL get_nrbe3pen_l(rbe3%nrbe3,irbe3_variables,rbe3%IRBE3,nrbe3pen_l)
605
606 CALL read_i_c(rbe3%LRBE3,rbe3%lrbe3_sz)
607
608 CALL read_i_c(iactiv,siactiv)
609
610 CALL read_i_c(ibmpc,sibmpc)
611
612 CALL read_i_c(nodes%KINET,numnod)
613
614 CALL read_i_c(fasolfr,sfasolfr)
615
616 CALL read_i_c(segquadfr,ssegquadfr)
617
618 CALL read_i_c(ipari,sipari)
619
620c CALL READ_I_C(INBUF,SINBUF)
621
622 CALL r_bufbric_22 ! I22LEN |-> BRICK_LIST(...); CALL READ_I_C(); CALL READ_DB(); END SUBROUTINE
623
624 !read new structure INTBUF_TAB
625 CALL intbuf_ini(interfaces%INTBUF_TAB)
626
627 CALL read_i_c(rwall%NPRW,rwall%NNPRW*rwall%NRWALL)
628
629 CALL read_i_c(rwall%LPRW,rwall%SZ_LPRW)
630
631 CALL read_i_c(iconx,siconx)
632
633 CALL read_i_c(madprt,smadprt)
634
635 CALL read_i_c(madsh4,smadsh4)
636
637 CALL read_i_c(madsh3 ,smadsh3)
638
639 CALL read_i_c(madsol,smadsol)
640
641 CALL read_i_c(madnod,smadnod)
642
643 CALL read_i_c(madclnod,madclnods)
644
645 CALL read_i_c(madfail,smadfail)
646
647 CALL read_i_c(madidx,smadidx)
648
649 CALL read_i_c(npby,snpby)
650
651 CALL read_i_c(npbyl,snpbyl)
652
653 CALL read_i_c(lpby,slpby)
654
655 CALL read_i_c(lpbyl,slpbyl)
656
657 CALL read_i_c(lrivet,slrivet)
658
659 CALL read_i_c(nstrf,snstrf)
660
661 CALL read_i_c(nom_sect,snom_sect)
662
663 CALL read_joint()
664
665 CALL read_i_c(nodpor,snodpor)
666
667 CALL read_i_c(monvol,smonvol)
668
669 CALL read_i_c(ntri, nvolu)
670 DO ii = 1, nvolu
671 t_monvol(ii)%NB_FILL_TRI = ntri(ii)
672 IF (ntri(ii) > 0) THEN
673 ALLOCATE(t_monvol(ii)%FILL_TRI(3 * ntri(ii)))
674 CALL read_i_c(t_monvol(ii)%FILL_TRI(1:3 * ntri(ii)), 3 * ntri(ii))
675 ENDIF
676 ENDDO
677
678 CALL read_i_c(lagbuf,slagbuf)
679
680 CALL read_i_c(icnds10,3*ns10e)
681
682 ! Read Time History
683 CALL read_th_restart(output%TH)
684
685 CALL thsurf_read_restart(output%TH%TH_SURF)
686
687 CALL read_i_c(fr_mv,sfr_mv)
688
689 CALL read_i_c(nodes%BOUNDARY_ADD,siad_elem)
690
691 CALL read_i_c(iad_rby,siad_rby)
692
693 CALL read_i_c(iad_rby2,siad_rby2)
694
695 CALL read_i_c(iad_i2m,siad_i2m)
696
697 CALL read_i_c(iadcj,siadcj)
698
699 CALL read_i_c(iad_rbm,siad_rbm)
700
701 CALL read_i_c(iad_rbm2,siad_rbm2)
702
703 CALL read_i_c(iad_rbe2,siad_rbe2)
704
705 CALL read_i_c(rbe3%mpi%IAD_RBE3,nspmd+1)
706
707 CALL read_i_c(iad_sec,siad_sec)
708
709 CALL read_i_c(iad_cut,siad_cut)
710C
711 CALL read_i_c(iad_rbym,siad_rbym)
712
713 CALL read_i_c(iad_rbym2,siad_rbym2)
714C
715 CALL read_i_c(nodes%BOUNDARY,sfr_elem)
716
717 CALL read_i_c(fr_rby,sfr_rby)
718
719 CALL read_i_c(rwall%FR_WALL,sfr_wall)
720
721 CALL read_i_c(fr_rby2,sfr_rby2)
722
723 CALL read_i_c(fr_i2m,sfr_i2m)
724
725 CALL read_i_c(fr_ll,sfr_ll)
726
727 CALL read_i_c(fr_cj,sfr_cj)
728
729 CALL read_i_c(fr_rbm,sfr_rbm)
730
731 CALL read_i_c(fr_rbm2,sfr_rbm2)
732
733 CALL read_i_c(fr_rbe2,sfr_rbe2)
734
735 CALL read_i_c(rbe3%mpi%FR_RBE3,rbe3%mpi%fr_rbe3_sz)
736
737 CALL read_i_c(rbe3%mpi%FR_RBE3MP,rbe3%mpi%fr_rbe3_sz)
738
739 CALL read_i_c(fr_sec,sfr_sec)
740
741 CALL read_i_c(fr_cut,sfr_cut)
742
743 CALL read_i_c(rg_cut,srg_cut)
744
745 CALL read_i_c(fr_mad,sfr_mad)
746
747 CALL read_i_c(fr_i18,sfr_i18)
748
749 CALL read_i_c(dd_r2r,sdd_r2r)
750C
751 CALL read_i_c(dd_r2r_elem,sdd_r2r_elem)
752C
753 CALL read_i_c(fr_rbym,sfr_rbym)
754
755 CALL read_i_c(fr_rbym2,sfr_rbym2)
756C
757 CALL read_i_c(addcsrect,numnor+1)
758
759 CALL read_i_c(fr_nor,nbddnort)
760
761 CALL read_i_c(iad_frnor,(nspmd+1)*ninter25)
762
763 CALL read_i_c(procnor,nbccnor)
764
765 CALL read_i_c(interfaces%SPMD_ARRAYS%FR_EDG,2*nbddedgt)
766
767 CALL read_i_c(interfaces%SPMD_ARRAYS%IAD_FREDG,(nspmd+1)*ninter25)
768
769C
770 CALL read_i_c(iskwp,siskwp)
771
772 CALL read_i_c(nskwp,snskwp)
773
774 IF(snskwp>0) CALL read_i_c(iskwp_l,nskwp(ispmd+1))
775
776 CALL read_i_c(isensp,sisensp)
777
778 CALL read_i_c(nsensp,snsensp)
779
780 CALL read_i_c(iaccp,siaccp)
781
782 CALL read_i_c(naccp,snaccp)
783
784 CALL read_i_c(igaup,sigaup)
785
786 CALL read_i_c(ngaup,sngaup)
787
788 CALL read_i_c(fr_lagf,sfr_lagf)
789
790 CALL read_i_c(newfront,snewfront)
791
792 CALL read_i_c(nbrcvois,snbrcvois)
793
794 CALL read_i_c(lnrcvois,slnrcvois)
795
796 CALL read_i_c(nbsdvois,snbsdvois)
797
798 CALL read_i_c(lnsdvois,slnsdvois)
799
800 CALL read_i_c(nercvois,snercvois)
801
802 CALL read_i_c(lercvois,slercvois)
803
804 CALL read_i_c(nesdvois,snesdvois)
805
806 CALL read_i_c(lesdvois,slesdvois)
807
808 CALL read_i_c(npsegcom,snpsegcom)
809
810 CALL read_i_c(lsegcom,slsegcom)
811
812 CALL read_i_c(nporgeo,snporgeo)
813
814 CALL read_i_c(lnodpor,slnodpor)
815
816 CALL read_i_c(llagf,sllagf)
817C
818 nodes%WEIGHT_MD = nodes%WEIGHT
819
820C
821 IF(icrack3d > 0)THEN
822 CALL read_i_c(iad_edge,siad_edge)
823 CALL read_i_c(fr_edge,sfr_edge)
824 CALL read_i_c(fr_nbedge,sfr_nbedge)
825 ENDIF
826C
827 CALL read_i_c(iad_cndm,siad_cndm)
828 CALL read_i_c(fr_cndm,sfr_cndm)
829C
830C Part Parith/We read in SPM in all cases
831C
832C direct reading of ADDCNE
833 IF(iparit==1) THEN
834 CALL read_i_c(element%PON%ADSKY,numnod+1)
835C direct reading of PROCNE
836 CALL read_i_c(element%PON%PROCNE,lcne0)
837C lecture directe structure interface type2
838 IF (i2nsnt > 0) THEN
839 CALL read_i_c(addcni2,numnod+1)
840 ENDIF
841 CALL read_i_c(procni2,lcni2)
842C---Itet2 of S10
843 IF(ns10e>0) THEN
844 CALL read_i_c(addcncnd,saddcncnd)
845 ENDIF
846 CALL read_i_c(procncnd,sprocncnd)
847C Direct reading of IADS+IADWAL+IADRBK+IADI2+IADMV2+IADMV3+IADLL+IADRBM+IADI18+
848C IADRBMK
849 CALL read_i_c(element%PON%IADS ,SIZE(element%PON%IADS)) !< 1 ; 8xNUMELS solid indexes to FSKY
850 CALL read_i_c(element%PON%IADS10 ,SIZE(element%PON%IADS10 )) !< 6* NUMELS10
851 CALL read_i_c(element%PON%IADS20 ,SIZE(element%PON%IADS20 )) ! 12*NUMELS20
852 CALL read_i_c(element%PON%IADS16 ,SIZE(element%PON%IADS16 )) ! 8*NUMELS16
853 CALL read_i_c(element%PON%IADQ ,SIZE(element%PON%IADQ )) !<i87b ; quad i87b
854 CALL read_i_c(element%PON%IADC ,SIZE(element%PON%IADC )) !<i87C shell (4 nodes) indexes to FSKY
855 CALL read_i_c(element%PON%IAD_TRUSS ,SIZE(element%PON%IAD_TRUSS)) !< I87D 2xNUMELT
856 CALL read_i_c(element%PON%IAD_BEAM ,SIZE(element%PON%IAD_BEAM )) !< I87E 2xNUMELP
857 CALL read_i_c(element%PON%IAD_SPRING ,SIZE(element%PON%IAD_SPRING )) !<F 3xNUMELR
858 CALL read_i_c(element%PON%IAD_TG ,SIZE(element%PON%IAD_TG )) !<G 3xNUMELTG
859 CALL read_i_c(element%PON%IAD_TG6 ,SIZE(element%PON%IAD_TG6 )) !<H 3xNUMELTG6
860 CALL read_i_c(element%PON%IAD_MV ,SIZE(element%PON%IAD_MV )) !I 4xNSKYMV0
861 CALL read_i_c(element%PON%IAD_CONLD ,SIZE(element%PON%IAD_CONLD)) !<J 4xNCONLD
862 CALL read_i_c(element%PON%IAD_CONV ,SIZE(element%PON%IAD_CONV)) !<K 4x glob_therm%NCONV
863 CALL read_i_c(element%PON%IAD_RADIA ,SIZE(element%PON%IAD_RADIA)) !<L 4x glob_therm%Numrada
864 CALL read_i_c(element%PON%IAD_LOADP ,SIZE(element%PON%IAD_LOADP)) !<M SLLOADP
865 CALL read_i_c(element%PON%IAD_FXFLUX ,SIZE(element%PON%IAD_FXFLUX)) !<N 4x glob_therm%nfxflux
866
867 CALL read_i_c(iadwal,siadwal)
868
869 CALL read_i_c(iadrbk,siadrbk)
870
871 CALL read_i_c(iadi2,siadi2)
872
873 CALL read_i_c(iadcnd,siadcnd)
874
875 CALL read_i_c(iadmv2,siadmv2)
876
877 CALL read_i_c(iadmv3,siadmv3)
878
879 CALL read_i_c(iadll,siadll)
880
881 CALL read_i_c(iadrbm,siadrbm)
882
883 CALL read_i_c(iadi18,siadi18)
884C
885 CALL read_i_c(iadrbmk,siadrbmk)
886C
887C ply xfem
888C
889 IF(iplyxfem > 0) THEN
890 CALL read_i_c(adsky_pxfem,nplyxfe+1 )
891C direct reading of PROCNE
892 CALL read_i_c(procne_pxfem,lcnepxfem)
893 CALL read_i_c(iadc_pxfem,4*eplyxfe)
894 ENDIF
895C
896C xfem for layered shell (cracks)
897C
898 IF (icrack3d > 0) THEN
899 CALL read_i_c(adsky_crkxfem,ncrkxfe+1)
900 CALL read_i_c(cne_crkxfem,lcnecrkxfem)
901 CALL read_i_c(procne_crkxfem,lcnecrkxfem)
902 CALL read_i_c(iadc_crkxfem,4*ecrkxfec+3*ecrkxfetg)
903 CALL read_i_c(crknodiad,lcnecrkxfem)
904 ENDIF
905C
906 ELSEIF(ipari0==1)THEN
907C skip reading
908 len = numnod+1+lcne0+(numnod+1)*min(1,i2nsnt)+lcni2
909 + + nisky0+nskyrw0+nskyrbk0+niskyi2+nskymv0+nskymvc0
910 + + nskyll0+nskyrbm0+nskyi18+nskyrbe30+nskyrbmk0
911C
912 IF(ns10e>0) len = len +saddcncnd+sprocncnd+siadcnd
913 IF(iplyxfem > 0)len = len + lcnepxfem + 4*eplyxfe+ nplyxfe+1
914 DO i = 1,len,mvsiz
915 CALL read_i_c(idumm,min(mvsiz,len-i+1))
916 ENDDO
917 ENDIF
918c ENDIF
919
920 ! ---------------
921 ! ebcs option : adress of FSKY array
922 IF(ebcs_tab%nebcs_parallel > 0) THEN
923 DO i=1,ebcs_tab%nebcs
924 is_ebcs_parallel = .false.
925 IF(ebcs_tab%my_typ(i)==10 .OR. ebcs_tab%my_typ(i)==11)is_ebcs_parallel = .true.
926 IF(is_ebcs_parallel) THEN
927 my_size = ebcs_tab%tab(i)%poly%nb_elem
928 ALLOCATE(ebcs_parithon(i)%ELEM_ADRESS(4,my_size))
929 CALL read_i_c(ebcs_parithon(i)%ELEM_ADRESS,4*my_size)
930 ENDIF
931 ENDDO
932 ENDIF
933 ! ---------------
934
935 CALL decompress_i_nnz(igeo,sigeo)
936
937 CALL decompress_i_nnz(ipm, sipm)
938
939 CALL read_i_c(icontact ,sicontact)
940c
941 CALL read_i_c(ipart_state ,npart)
942C--------------------------------------
943 IF(nadmesh /= 0)THEN
944 CALL read_i_c(sh4tree,ksh4tree*numelc)
945 CALL read_i_c(sh3tree,ksh3tree*numeltg)
946 CALL read_i_c(ipadmesh,kipadmesh*npart)
947 len=abs(lsh4trim)
948 CALL read_i_c(sh4trim,len)
949 len=abs(lsh3trim)
950 CALL read_i_c(sh3trim,len)
951 END IF
952C--------------------------------------
953C Allowance + Reading Structure Interface Entiere Party
954C--------------------------------------
955 CALL spmd_initfi(ipari,1,h3d_data,interfaces%PARAMETERS,
956 . glob_therm%IDT_THERM,glob_therm%INTHEAT)
957C
958C--------------------------------------
959C Reading + Intit Frontiere SPMD for SPH (integer part)
960C--------------------------------------
961 CALL spmd_inisph(1)
962C--------------------------------------
963 IF(nintstamp /= 0)THEN
964 CALL intstamp_rresti(intstamp)
965 END IF
966C--------------------------------------
967C Interface friction buffer (reak )
968C--------------------------------------
969 IF(ninterfric /= 0)THEN
970 CALL intfric_rresti(interfaces%INTBUF_FRIC_TAB,ninterfric)
971 END IF
972C--------------------------------------
973C Table structure
974C--------------------------------------
975 IF(ntable /= 0)THEN
976 CALL table_rresti(table)
977 END IF
978C--------------------------------------
979 CALL read_i_c(ithvar,sithvar)
980C--------------------------------------
981 IF (nloadc>0)THEN
982 CALL read_i_c(icfield,sizfield*nloadc)
983 CALL read_i_c(lcfield,slcfield)
984 ENDIF
985 IF (nloadp>0)THEN
986 CALL read_i_c(iloadp,sizloadp*nloadp)
987 CALL read_i_c(lloadp,slloadp)
988 ENDIF
989 IF (pblast%NLOADP_B>0)THEN
990 CALL pblast_load(pblast)
991 CALL pblast_init_tables(pblast%PBLAST_DATA) !PBLAST_DATA defined in PBLAST_MOD
992 ENDIF
993c
994 IF (loads%NLOAD_CYL > 0) THEN
995 CALL read_pcyl(loads)
996 ENDIF
997C----------------------------------
998 IF(icrack3d > 0)THEN
999 CALL read_i_c(ibordnode,sibordnode)
1000 CALL read_i_c(iedgesh,siedgesh)
1001 CALL read_i_c(nodedge,2*numedges)
1002 CALL read_i_c(iedge,numedges)
1003 CALL read_i_c(nodglobxfe,snodglobxfe)
1004 END IF
1005C----------------------------------
1006 CALL read_i_c(kxig3d,skxig3d)
1007 CALL read_i_c(ixig3d,sixig3d)
1008C----------------------------------
1009 CALL read_i_c(tag_skins6,numels)
1010C----------------------------------
1011C /BCS/CYCLIC
1012C----------------------------------
1013 IF(nbcscyc > 0)THEN
1014 CALL read_i_c(ibcscyc,sibcscyc)
1015 CALL read_i_c(lbcscyc,slbcscyc)
1016 END IF
1017C----------------------------------
1018C /BCS/WALL
1019C----------------------------------
1020 CALL read_i_c(itmp,1)
1021 bcs%NUM_WALL = itmp
1022 IF(bcs%NUM_WALL > 0)THEN
1023 CALL read_bcs_wall()
1024 ENDIF
1025C----------------------------------
1026C /BCS/NRF
1027C----------------------------------
1028 CALL read_i_c(itmp,1)
1029 bcs%NUM_NRF = itmp
1030 IF(bcs%NUM_NRF > 0)THEN
1031 CALL read_bcs_nrf()
1032 ENDIF
1033C----------------------------------
1034C reading real variables
1035C--------------------------------------
1036C warning to radioss developers :
1037C IF(LVARREA>LMXVREA)THEN
1038C WRITE(*,*) ' ** RADIOSS INTERNAL ERROR'
1039C CALL ARRET(7)
1040C ENDIF
1041C-----
1042 CALL read_db(tabvrea,lvarrea)
1043C
1044C-----
1045 CALL rdcomr(lmxvrea,tabvrea,dynain_data,interfaces%PARAMETERS,
1046 . output,glob_therm)
1047C--------------------------------------
1048C READING REALS
1049C--------------------------------------
1050 CALL read_db(nodes%X,sx)
1051 CALL read_db(nodes%D,sd)
1052 CALL read_db(nodes%V,sv)
1053 CALL read_db(nodes%VR,svr)
1054 CALL read_db(nodes%DR,sdr)
1055
1056 CALL read_db(thke,sthke)
1057
1058 CALL read_db(dampr,sdampr)
1059
1060 CALL read_db(damp,sdamp)
1061
1062 CALL read_db(nodes%MS,sms)
1063
1064 IF (n2d >0) THEN
1065 CALL read_db(ms_2d,numnod)
1066 ENDIF
1067
1068 CALL read_db(nodes%IN,sin)
1069
1070 CALL read_db(tf,stf)
1071
1072 CALL read_db(pm,spm)
1073
1074 CALL read_matparam(mat_elem)
1075
1076 CALL read_elgroup_param(mat_elem)
1077
1078 skew_len = skews%N_SKEW_VAR*skews%TOTAL_SKEW_NUMBER
1079 CALL read_db(skews%SKEW,skew_len)
1080
1081 CALL read_db(xframe,sxframe)
1082
1083 CALL decompress_r_nnz(geo,sgeo)
1084
1085 CALL read_db(eani,seani)
1086
1087 CALL read_db(bufmat,sbufmat)
1088
1089 CALL read_db(bufgeo,sbufgeo)
1090
1091 CALL read_db(bufsf,sbufsf)
1092
1093 CALL read_db(rbmpc,srbmpc)
1094
1095 CALL read_db(gjbufr,sgjbufr)
1096
1097 CALL read_db(w,sw)
1098
1099 CALL read_db(veul,sveul)
1100
1101 IF (multi_fvm%IS_USED) THEN
1102 CALL read_db(multi_fvm%PRES_SHIFT, 1)
1103 IF (n2d == 0) THEN
1104 CALL read_db(multi_fvm%VEL(1, :), numels)
1105 CALL read_db(multi_fvm%VEL(2, :), numels)
1106 CALL read_db(multi_fvm%VEL(3, :), numels)
1107 ELSE
1108 CALL read_db(multi_fvm%VEL(1, :), numelq + numeltg)
1109 CALL read_db(multi_fvm%VEL(2, :), numelq + numeltg)
1110 CALL read_db(multi_fvm%VEL(3, :), numelq + numeltg)
1111 ENDIF
1112 ENDIF
1113
1114 CALL read_db(fill,sfill)
1115
1116 CALL read_db(dfill,sdfill)
1117
1118 CALL read_db(alph,salph)
1119
1120 CALL read_db(wb,swb)
1121
1122 CALL read_db(dsave,sdsave)
1123
1124 CALL read_db(asave,sasave)
1125
1126 CALL read_db(spbuf,sspbuf)
1127
1128 CALL read_db(vsphio,svsphio)
1129
1130 CALL read_db(sphveln,ssphveln)
1131
1132 CALL read_db(msnf,smsnf)
1133
1134 CALL read_db(forc,sforc)
1135
1136 CALL read_db(vel,svel)
1137
1138 CALL read_db(fsav,sfsav)
1139
1140 CALL read_i_c(tab_ump,7*taille)
1141
1142 CALL read_i_c(poin_ump,nummat)
1143
1144 CALL read_db(tab_mat,stab_mat)
1145
1146 CALL read_db(fzero,sfzero)
1147
1148 CALL read_db(xlas,sxlas)
1149
1150 CALL read_db(accelm,saccelm)
1151
1152 CALL read_db(gauge,llgauge*nbgauge)
1153
1154 CALL read_sensors(sensors, python)
1155
1156 CALL read_db(fbvel,sfbvel)
1157
1158 CALL read_db(rbe3%FRBE3,rbe3%frbe3_sz)
1159
1160 CALL read_rrbe3pen(rbe3%PEN,nrbe3pen_l)
1161
1162 CALL read_db(factiv,lractiv*nactiv)
1163
1164 CALL read_db(grav,sgrav)
1165
1166 CALL read_db(fr_wave,sfr_wave)
1167
1168 CALL read_failwave(failwave)
1169
1170 CALL read_nloc_struct(nloc_dmg)
1171
1172 CALL read_db(parts0,sparts0)
1173
1174 CALL read_db(elbuf,selbuf)
1175
1176 CALL elbuf_ini(mat_elem,ngroup, iparg) ! allocate, read and initialize element buffer
1177
1178 CALL read_cluster(cluster)
1179
1180 CALL read_db(rwall%RWBUF,rwall%nrwlp*rwall%nrwall)
1181
1182 CALL read_db(rwall%RWSAV,rwall%SZ_RWSAV)
1183
1184 CALL get_nrwallpen_l(rwall%nrwall,rwall%nnprw,rwall%nprw,rwall%nrwall_pen,rwall%pen%lnspen)
1185
1186 CALL read_rrwallpen(rwall%nrwall_pen,rwall%pen,rwall%pen%lnspen)
1187
1188 CALL read_db(rby,srby)
1189
1190 CALL read_db(rbyl,srbyl)
1191
1192 CALL read_db(rivet,srivet)
1193
1194 CALL read_db(secbuf,ssecbuf)
1195
1196 CALL read_db(volmon,svolmon)
1197
1198 CALL read_db(lambda,slambda)
1199
1200 CALL read_db(rconx,srconx)
1201
1202 CALL read_db(rcontact ,srcontact)
1203 CALL read_db(acontact ,srcontact)
1204 CALL read_db(pcontact ,srcontact)
1205
1206C lecture struture interface partie reelle si irunn > 1
1207 IF(irunn>1.OR.mcheck/=0) CALL spmd_initfi(ipari,2,h3d_data,interfaces%PARAMETERS,
1208 . glob_therm%IDT_THERM,glob_therm%INTHEAT)
1209C--------------------------------------
1210C Reading + Intit Frontier SPMD for SPH (real part)
1211C--------------------------------------
1212 CALL spmd_inisph(2)
1213C--------------------------------------
1214C STOCKAGE FACULTATIFS
1215C--------------------------------------
1216 IF (iresmd==0.AND.npsav>=25) THEN
1217 CALL read_db(partsav,npsav*npart)
1218 ENDIF
1219
1220 IF(nnoise>0.AND. (irunn>1 .OR. mcheck/=0))THEN
1221c IF(NNOISE>0.AND.IRUNN>1)THEN
1222 CALL read_i_c (iaf(if01),sinoise)
1223 CALL read_db(af(mf01),sfnoise)
1224 ENDIF
1225
1226 IF (iresp == 1) THEN
1227 IF (irxdp == 1) THEN
1228 CALL read_dpdb(nodes%XDP,3*numnod)
1229 CALL read_dpdb(nodes%DDP,3*numnod)
1230 ELSE
1231 CALL fillxdp(nodes%X,nodes%XDP,nodes%D,nodes%DDP)
1232C next restart, XDP will be on restart head !
1233 irxdp=1
1234 ENDIF
1235C Cas double precision
1236 ELSE
1237 IF (irxdp ==1)THEN
1238C Reading Dummy if XDP is in Restart file. XDP written during next restart.
1239 ALLOCATE(xdpdum(3*numnod))
1240C Read XDP
1241 CALL read_db(xdpdum,3*numnod)
1242C Read DDP
1243 CALL read_db(xdpdum,3*numnod)
1244 DEALLOCATE(xdpdum)
1245 irxdp=0
1246 ENDIF
1247 ENDIF
1248 npsav = 29
1249C ----------------------------------------------------------------------
1250C used for checkpoint restarts.
1251 IF(mcheck/=0)THEN
1252 CALL read_i_c(my_ilen,1)
1253 CALL read_i_c(my_rlen,1)
1254 IF(mod(irform,5)>=2.AND.mod(irform,5)<=4) THEN
1255 CALL read_i_c(ilink,4*nrlink)
1256 CALL read_i_c(fr_rl,(nspmd+2)*nrlink)
1257 CALL read_i_c(llink,lllink)
1258 CALL read_i_c(iadrl,min(1,iparit)*lllink)
1259 !CALL READ_I_C(LINALE,LLINAL)
1260 CALL read_i_c(neflsw,nsflsw)
1261 CALL read_i_c(nnflsw,8*ntflsw)
1262 CALL read_i_c(icut,44*ncuts)
1263 CALL read_db(crflsw,9*nsflsw)
1264 CALL read_db(xcut,7*ncuts)
1265 ENDIF
1266 ELSE
1267 CALL read_i_c(my_ilen,1)
1268 CALL read_i_c(my_rlen,1)
1269 ALLOCATE(ibid(my_ilen),stat=istat)
1270 IF(istat/=0)THEN
1271 CALL ancmsg(msgid=20,anmode=aninfo)
1272 CALL arret(2)
1273 END IF
1274 ALLOCATE(rbid(my_rlen),stat=istat)
1275 IF(istat/=0)THEN
1276 CALL ancmsg(msgid=20,anmode=aninfo)
1277 CALL arret(2)
1278 END IF
1279 CALL read_i_c(ibid,my_ilen)
1280 CALL read_db (rbid,my_rlen)
1281 ENDIF
1282C ----------------------------------------------------------------------
1283 IF(nadmesh /= 0 .OR. irest_mselt /= 0)THEN
1284 CALL read_db(msc,numelc)
1285 CALL read_db(inc,numelc)
1286 CALL read_db(mstg,numeltg)
1287 CALL read_db(intg,numeltg)
1288 CALL read_db(ptg,3*numeltg)
1289 END IF
1290
1291 IF(irest_mselt /= 0)THEN
1292 CALL read_db(mssa,numels)
1293 CALL read_db(mstr,numelt)
1294 CALL read_db(msp,numelp)
1295 CALL read_db(msrt,numelr)
1296 END IF
1297
1298 IF(nadmesh /= 0)THEN
1299 CALL read_db(padmesh,kpadmesh*npart)
1300 IF(glob_therm%ITHERM_FE > 0)THEN
1301 CALL read_db(mcpc,numelc)
1302 CALL read_db(mcptg,numeltg)
1303 END IF
1304 END IF
1305C pinching
1306 IF(npinch > 0) THEN
1307 CALL read_db(pinch_data%VPINCH,3*npinch)
1308 CALL read_db(pinch_data%DPINCH,3*npinch)
1309 CALL read_db(pinch_data%XPINCH,3*npinch)
1310 CALL read_db(pinch_data%MSPINCH,npinch)
1311 ENDIF
1312C ----------------------------------------------------------------------
1313 IF(istatcnd /= 0)THEN
1314 CALL read_db(mscnd,numnod)
1315 CALL read_db(incnd,numnod)
1316 ENDIF
1317C ----------------------------------------------------------------------
1318 IF(nintstamp /= 0)THEN
1319 CALL intstamp_rrestr(intstamp)
1320 END IF
1321 IF(nintskidold > 0) THEN
1322 IF(nintstamp/=0) THEN
1323 CALL read_db(pskids, nintskidold*numnodg)
1324 ELSE
1325 CALL read_db(pskids, nintskidold*numnod)
1326 ENDIF
1327 ENDIF
1328 IF(ifcontmax >0)THEN
1329 IF(nintstamp/=0.AND.nspmd > 1 ) THEN
1330 IF(ispmd == 0) CALL read_db(fcont_max, 3*numnodg)
1331 ELSE
1332 CALL read_db(fcont_max, 3*numnod)
1333 ENDIF
1334 ENDIF
1335 IF(ifcontpmax >0)THEN
1336 CALL read_db(fncont_max, 3*spcont_max)
1337 CALL read_db(ftcont_max, 3*spcont_max)
1338 ENDIF
1339 IF(ifcont2max >0)THEN
1340 CALL read_db(fcont2_max, 3*numnod)
1341 ENDIF
1342 IF(ifcontp2max >0)THEN
1343 CALL read_db(fncont2_max, 3*numnod)
1344 CALL read_db(ftcont2_max, 3*numnod)
1345 CALL read_db(npcont2_max, 3*numnod)
1346 ENDIF
1347 IF(ifcont2min >0)THEN
1348 CALL read_db(fcont2_min, 3*numnod)
1349 ENDIF
1350 IF(ifcontp2min >0)THEN
1351 CALL read_db(fncont2_min, 3*numnod)
1352 CALL read_db(ftcont2_min, 3*numnod)
1353 CALL read_db(npcont2_min, 3*numnod)
1354 ENDIF
1355C--------Frictional energy output-------
1356 IF(output%DATA%S_EFRICINT >0)THEN
1357 IF(output%DATA%NINEFRIC > 0) CALL read_db(output%DATA%EFRIC, output%DATA%NINEFRIC*numnod)
1358 IF(output%DATA%NINEFRIC_STAMP > 0) CALL read_db(output%DATA%EFRIC_STAMP, output%DATA%NINEFRIC_STAMP*numnodg)
1359 ENDIF
1360 IF(output%DATA%S_EFRIC >0)THEN
1361 CALL read_db(output%DATA%EFRICG, numnod)
1362 IF(nintstamp/=0) CALL read_db(output%DATA%EFRICG_STAMP, numnodg)
1363 ENDIF
1364C--------------------------------------
1365C Interface friction buffer (Reals )
1366C--------------------------------------
1367 IF(ninterfric /= 0)THEN
1368 CALL intfric_rrestr(interfaces%INTBUF_FRIC_TAB,ninterfric)
1369 END IF
1370C-------------
1371C always
1372 CALL read_db(nodes%MS0,numnod)
1373 CALL read_db(nodes%IN0,sin)
1374 IF(idtmins_old==1) THEN
1375 CALL read_db(admsms,numnod)
1376 ELSEIF(idtmins_old==2) THEN
1377 CALL read_db(dmelc ,numelc)
1378 CALL read_db(dmeltg,numeltg)
1379 CALL read_db(dmels ,numels)
1380 CALL read_db(dmeltr,numelt)
1381 CALL read_db(dmelp ,numelp)
1382 CALL read_db(dmelrt,numelr)
1383 END IF
1384 IF(idtmins_old/=0.OR.idtmins_int_old/=0)THEN
1385 CALL read_db(res_sms,3*numnod)
1386 END IF
1387 IF(idtmins_old==2.OR.idtmins_int_old/=0) THEN
1388Cfor TYPE20
1389 CALL read_db(diag_sms ,numnod)
1390Cfor TYPE2
1391 CALL read_db(dmint2 ,4*i2nsn25)
1392 END IF
1393 IF (isms_selec /= 0) THEN
1394 CALL read_i_c(nativ0_sms,numnod)
1395 ENDIF
1396C--------------------------------------
1397C Structure Tables Real Party
1398C--------------------------------------
1399 IF(ntable /= 0)THEN
1400 CALL table_rrestr(table)
1401 END IF
1402C------------------------
1403C Read MDS Parameters
1404C------------------------
1405#ifdef DNC
1406 CALL eng_read_mds()
1407#endif
1408C--------------------------------------
1409C ALE LINK
1410C--------------------------------------
1411! print *,"rdresb.F: READ_I_C, SLINALE=", SLINALE
1412 IF (slinale>0)CALL read_i_c(linale,slinale)
1413C--------------------------------------
1414C ALE LINK
1415C--------------------------------------
1416 CALL read_ale_grid()
1417C--------------------------------------
1418C FXBODY (moved from radioss2.F)
1419C--------------------------------------
1420 IF (nfxbody>0) THEN
1421 CALL fxbresb()
1422 ENDIF
1423C--------------------------------------
1424C EIGEN MODES (moved from radioss2.F)
1425C--------------------------------------
1426 IF (neig>0) THEN
1427 CALL eigrrest()
1428 ENDIF
1429C--------------------------------------
1430C SEATBELTS - SLIPRING / RETRACTOR-
1431C--------------------------------------
1432 DO i = 1, nslipring
1433C
1434 CALL read_i_c(slipring(i)%ID, 1)
1435 CALL read_i_c(slipring(i)%IDG, 1)
1436 CALL read_i_c(slipring(i)%NFRAM, 1)
1437 CALL read_i_c(slipring(i)%IFUNC, 4)
1438 CALL read_i_c(slipring(i)%SENSID, 1)
1439 CALL read_i_c(slipring(i)%FL_FLAG, 1)
1440 CALL read_i_c(slipring(i)%RBODY, 1)
1441C
1442 CALL read_db(slipring(i)%DC, 1)
1443 CALL read_db(slipring(i)%A, 1)
1444 CALL read_db(slipring(i)%FRIC, 1)
1445 CALL read_db(slipring(i)%FAC_D, 3)
1446 CALL read_db(slipring(i)%FRICS, 1)
1447 CALL read_db(slipring(i)%FAC_S, 3)
1448C
1449 ALLOCATE(slipring(i)%FRAM(slipring(i)%NFRAM))
1450C
1451 DO j=1,slipring(i)%NFRAM
1452C
1453 CALL read_i_c(slipring(i)%FRAM(j)%UPDATE, 1)
1454 CALL read_i_c(slipring(i)%FRAM(j)%ANCHOR_NODE, 1)
1455 CALL read_i_c(slipring(i)%FRAM(j)%NODE, 3)
1456 CALL read_i_c(slipring(i)%FRAM(j)%NODE_NEXT, 3)
1457 CALL read_i_c(slipring(i)%FRAM(j)%NODE2_PREV, 1)
1458 CALL read_i_c(slipring(i)%FRAM(j)%N_REMOTE_PROC, 1)
1459 CALL read_i_c(slipring(i)%FRAM(j)%ORIENTATION_NODE, 1)
1460 CALL read_i_c(slipring(i)%FRAM(j)%STRAND_DIRECTION, 2)
1461 CALL read_i_c(slipring(i)%FRAM(j)%LOCKED, 1)
1462C
1463 CALL read_db(slipring(i)%FRAM(j)%VECTOR, 6)
1464 CALL read_db(slipring(i)%FRAM(j)%ORIENTATION_ANGLE, 1)
1465 CALL read_db(slipring(i)%FRAM(j)%MATERIAL_FLOW, 1)
1466 CALL read_db(slipring(i)%FRAM(j)%MATERIAL_FLOW_OLD, 1)
1467 CALL read_db(slipring(i)%FRAM(j)%DFS, 1)
1468 CALL read_db(slipring(i)%FRAM(j)%RESIDUAL_LENGTH, 2)
1469 CALL read_db(slipring(i)%FRAM(j)%CURRENT_LENGTH, 2)
1470 CALL read_db(slipring(i)%FRAM(j)%RINGSLIP, 1)
1471 CALL read_db(slipring(i)%FRAM(j)%BETA, 1)
1472 CALL read_db(slipring(i)%FRAM(j)%SLIP_FORCE, 3)
1473 CALL read_db(slipring(i)%FRAM(j)%PREV_REF_LENGTH, 1)
1474 CALL read_db(slipring(i)%FRAM(j)%INTVAR_STR1, 8)
1475 CALL read_db(slipring(i)%FRAM(j)%INTVAR_STR2, 8)
1476C
1477 ENDDO
1478C
1479 ENDDO
1480C
1481 DO i = 1, nretractor
1482 CALL read_i_c(retractor(i)%ID, 1)
1483 CALL read_i_c(retractor(i)%IDG, 1)
1484 CALL read_i_c(retractor(i)%UPDATE, 1)
1485 CALL read_i_c(retractor(i)%ANCHOR_NODE, 1)
1486 CALL read_i_c(retractor(i)%NODE, 2)
1487 CALL read_i_c(retractor(i)%NODE_NEXT, 2)
1488 CALL read_i_c(retractor(i)%STRAND_DIRECTION, 1)
1489 CALL read_i_c(retractor(i)%IFUNC, 3)
1490 CALL read_i_c(retractor(i)%ISENS, 2)
1491 CALL read_i_c(retractor(i)%TENS_TYP, 1)
1492 CALL read_i_c(retractor(i)%LOCKED, 1)
1493 CALL read_i_c(retractor(i)%LOCKED_FREEZE, 1)
1494 CALL read_i_c(retractor(i)%PRETENS_ACTIV, 1)
1495 CALL read_i_c(retractor(i)%INACTI_NNOD, 1)
1496 CALL read_i_c(retractor(i)%INACTI_NNOD_MAX, 1)
1497 ALLOCATE(retractor(i)%INACTI_NODE(retractor(i)%INACTI_NNOD_MAX))
1498 CALL read_i_c(retractor(i)%INACTI_NODE, retractor(i)%INACTI_NNOD)
1499 CALL read_i_c(retractor(i)%N_REMOTE_PROC, 1)
1500 CALL read_i_c(retractor(i)%S_TABLE, 2)
1501 CALL read_db(retractor(i)%VECTOR, 3)
1502 CALL read_db(retractor(i)%ELEMENT_SIZE, 1)
1503 CALL read_db(retractor(i)%FORCE, 1)
1504 CALL read_db(retractor(i)%MATERIAL_FLOW, 1)
1505 CALL read_db(retractor(i)%RESIDUAL_LENGTH, 1)
1506 CALL read_db(retractor(i)%FAC, 4)
1507 CALL read_db(retractor(i)%PULLOUT, 1)
1508 CALL read_db(retractor(i)%UNLOCK_FORCE, 1)
1509 CALL read_db(retractor(i)%LOCK_PULL, 1)
1510 CALL read_db(retractor(i)%LOCK_PULL_SAV, 1)
1511 CALL read_db(retractor(i)%LOCK_OFFSET, 1)
1512 CALL read_db(retractor(i)%LOCK_YIELD_FORCE, 1)
1513 CALL read_db(retractor(i)%RINGSLIP, 1)
1514 CALL read_db(retractor(i)%PRETENS_TIME, 1)
1515 CALL read_db(retractor(i)%PRETENS_PULL, 1)
1516 CALL read_db(retractor(i)%PRETENS_PULLMAX, 1)
1517 CALL read_db(retractor(i)%RET_FORCE, 1)
1518 DO j=1,2
1519 IF (retractor(i)%S_TABLE(j) > 0) THEN
1520 npt = retractor(i)%S_TABLE(j)
1521 retractor(i)%TABLE(j)%NDIM = 1
1522 ALLOCATE (retractor(i)%TABLE(j)%X(1))
1523 ALLOCATE (retractor(i)%TABLE(j)%X(1)%VALUES(npt))
1524 ALLOCATE (retractor(i)%TABLE(j)%Y)
1525 ALLOCATE (retractor(i)%TABLE(j)%Y%VALUES(npt))
1526 CALL read_db(retractor(i)%TABLE(j)%X(1)%VALUES(1:npt),npt)
1527 CALL read_db(retractor(i)%TABLE(j)%Y%VALUES(1:npt),npt)
1528 ENDIF
1529 ENDDO
1530 ENDDO
1531C
1532 IF (n_anchor_remote > 0) THEN
1533 CALL read_i_c(anchor_remote%ADD_PROC, nspmd+1)
1534 CALL read_i_c(anchor_remote%NODE, n_anchor_remote)
1535 ENDIF
1536C
1537 IF (n_anchor_remote_send > 0) THEN
1538 CALL read_i_c(anchor_remote_send%ADD_PROC, nspmd+1)
1539 CALL read_i_c(anchor_remote_send%NODE, n_anchor_remote_send)
1540 ENDIF
1541C
1542 IF ((nslipring_g + nretractor_g >0).AND.(ispmd == 0)) THEN
1543 ALLOCATE(th_slipring(nslipring_g,6))
1544 ALLOCATE(th_retractor(nretractor_g,3))
1545 th_slipring = 0
1546 th_retractor = 0
1547 CALL read_i_c(nseatbelt_th_proc, 1)
1548 IF ((nspmd >1).AND.(nseatbelt_th_proc > 0)) THEN
1549 ALLOCATE(seatbelt_th_exch(nseatbelt_th_proc))
1550 DO i=1,nseatbelt_th_proc
1551 CALL read_i_c(seatbelt_th_exch(i)%ID_PROC, 1)
1552 CALL read_i_c(seatbelt_th_exch(i)%ADD_PROC, 1)
1553 CALL read_i_c(seatbelt_th_exch(i)%NSLIPRING, 1)
1554 CALL read_i_c(seatbelt_th_exch(i)%NRETRACTOR, 1)
1555 ENDDO
1556 ENDIF
1557 ENDIF
1558C--------------------------------------
1559C /H3D/?/TMAX
1560C--------------------------------------
1561 CALL read_db(tm_dis,lmax_dis*numnod)
1562 CALL read_db(tm_vel,lmax_vel*numnod)
1563 CALL read_db(tm_nsig1,lmax_nsig*numnod)
1564 CALL read_db(tm_nsig3,lmax_nsig*numnod)
1565 CALL read_db(tm_nstra1,lmax_nstra*numnod)
1566 CALL read_db(tm_nstra3,lmax_nstra*numnod)
1567C--------------------------------------
1568 IF (nconld > 0) THEN
1569 CALL read_db(dpl0cld,6*nconld)
1570 CALL read_db(vel0cld,6*nconld)
1571 ENDIF
1572C--------------------------------------
1573C /DAMP/VREL
1574C--------------------------------------
1575 ndamp_vrel_rbyg = 0
1576 IF (ndamp_vrel > 0) THEN
1577 CALL read_i_c(id_damp_vrel,ndamp_vrel)
1578 CALL read_i_c(fr_damp_vrel,ndamp_vrel*(nspmd+2))
1579 ENDIF
1580 IF (ndamp > 0) THEN
1581 DO i=1,ndamp
1582 IF (nint(dampr(nrdamp*(i-1)+25)) > 0) ndamp_vrel_rbyg = ndamp_vrel_rbyg + 1
1583 ENDDO
1584 ENDIF
1585
1586 ! Flow
1587 IF (nflow>0) THEN
1588 CALL nfrrest(iflow, rflow)
1589 ENDIF
1590
1591 ! FVMBAG Structures
1592 CALL fvrrest(monvol)
1593
1594 ! Implicit
1595 IF (impl_s0>0) CALL imprrest(impl_s0)
1596
1597 ! heat transfer
1598 IF(glob_therm%ITHERM_FE > 0 ) THEN
1599 CALL thcrrest(nodes%MCP,nodes%TEMP)
1600 ENDIF
1601
1602 ! Nitsche formulation
1603 IF(nitsche > 0 ) THEN
1604 CALL nitscherrest(forneqs)
1605 ENDIF
1606
1607 ! Units
1608 CALL read_units(unitab)
1609
1610 ! Rigid MAterials
1611 IF(irigid_mat > 0) CALL rmatrrest(rbym,irbym, lnrbym,weight_rm)
1612
1613 ! thermal bc
1614 IF (glob_therm%NUMCONV > 0) CALL convrrest(ibcv, fconv ,glob_therm)
1615 IF (glob_therm%NUMRADIA > 0) CALL radiarrest(ibcr, fradia,glob_therm)
1616 IF (glob_therm%NFXFLUX > 0) CALL fxfluxrrest(ibfflux, fbfflux,glob_therm)
1617 IF (glob_therm%NFXTEMP > 0) CALL fxtemprrest(ibftemp, fbftemp,glob_therm)
1618
1619 ! PlyXFEM
1620 IF (iplyxfem > 0) THEN
1621 CALL plyxfem_rrest(ms_ply,zi_ply,inod_pxfem,iel_pxfem,
1622 . icode_ply,iskew_ply,msz2)
1623 ENDIF
1624 CALL plyxfem_rrestanim()
1625 IF(iplyxfem > 0) CALL plyxfem_ravuply()
1626
1627 ! Loads
1628 IF(nloadc > 0) CALL cfieldrest(cfield)
1629 IF(nloadp > 0) CALL loadprest(loadp)
1630
1631 ! Interface GAP Load
1632 IF (nintloadp>0) THEN
1633 CALL read_db(dgaploadint, ninter*nloadp_hyd )
1634 ENDIF
1635
1636! /INIVEL w/ T_start or senor_id
1637 IF (loads%NINIVELT>0) THEN
1638 ALLOCATE(loads%INIVELT(loads%NINIVELT))
1639 CALL read_inivel(loads%NINIVELT,loads%INIVELT)
1640 ENDIF
1641 ! CrackXFEM
1642 IF(icrack3d > 0) CALL crkxfem_rrest(inod_crkxfem,iel_crkxfem,nodlevxf)
1643 CALL crkxfem_rrestanim()
1644
1645 ! Finite Volume Method / ALE
1646 IF(alefvm_param%IEnabled > 0) CALL read_db(alefvm_buffer%FCELL ,6*numels)
1647
1648 ! ALELAG
1649 IF(ialelag > 0) CALL alelag_rrest()
1650
1651 IF(srthbuf > 0) CALL rthbufrest(rthbuf,srthbuf)
1652
1653 IF(sknot > 0) CALL knotrest(knot)
1654 IF(sknotlocpc > 0) CALL knotlocpcrest(knotlocpc)
1655 IF(sknotlocel > 0) CALL knotlocelrest(knotlocel)
1656 IF(numelig3d > 0) CALL wigerest(wige)
1657
1658 IF(ipart_stack >0) THEN
1659 CALL stack_rrest(stack%IGEO,stack%GEO, stack%PM )
1660 ENDIF
1661
1662 IF (ndrape > 0) THEN
1663 CALL drape_rrest(drape_sh4n,drape_sh3n, drapeg)
1664 ENDIF
1665
1666 IF (numply > 0) CALL ply_info_rest(ply_info)
1667
1668 ! read the starter elapsed time
1669 CALL read_dpdb(global_comp_time%STARTER_TIME,1)
1670 ! read the previous engine elapsed time
1671 IF(global_comp_time%RUN_NBR>1) THEN
1672 CALL read_dpdb(global_comp_time%ENGINE_TIME,global_comp_time%RUN_NBR-1)
1673 ENDIF
1674C
1675 RETURN
1676 END
1677!||====================================================================
1678!|| fxbresb ../engine/source/output/restart/rdresb.F
1679!||--- called by ------------------------------------------------------
1680!|| rdresb ../engine/source/output/restart/rdresb.F
1681!||--- calls -----------------------------------------------------
1682!|| my_getpid ../engine/source/system/pid.c
1683!|| read_db ../common_source/tools/input_output/read_db.F
1684!|| read_i_c ../common_source/tools/input_output/write_routines.c
1685!|| redsqi ../engine/source/output/tools/redsqi.F
1686!|| redsqr ../engine/source/output/tools/redsqr.F
1687!|| tmpenvf ../engine/source/system/tmpenv_c.c
1688!||--- uses -----------------------------------------------------
1689!|| fxb_mod ../engine/share/modules/fxb_mod.F
1690!|| scratchfile_mod ../engine/share/modules/scratchfiles_mod.F
1691!||====================================================================
1692 SUBROUTINE fxbresb()
1693C-----------------------------------------------
1694C M o d u l e s
1695C-----------------------------------------------
1696 USE scratchfile_mod
1697 USE fxb_mod
1698C-----------------------------------------------
1699C I m p l i c i t T y p e s
1700C-----------------------------------------------
1701#include "implicit_f.inc"
1702C-----------------------------------------------
1703C C o m m o n B l o c k s
1704C-----------------------------------------------
1705#include "com04_c.inc"
1706#include "scr05_c.inc"
1707#include "units_c.inc"
1708#include "fxbcom.inc"
1709#include "chara_c.inc"
1710#include "task_c.inc"
1711C-----------------------------------------------
1712C L o c a l V a r i a b l e s
1713C-----------------------------------------------
1714 INTEGER LEN_IPM,LEN_MOD, NRECM, NRECS, IRCM, IRCM0, IRCM1,
1715 . IRCS, IRCS0, IRCS1, LREC, I, J, RCLEN,TLEN
1716 my_real
1717 . VV(6)
1718 INTEGER MY_PID
1719 CHARACTER(LEN=10) :: CPID
1720 CHARACTER(LEN=6) :: CISPMD
1721 CHARACTER(LEN=4096) :: TMPDIR
1722C--------------------------------------
1723C READING INTEGERS
1724C--------------------------------------
1725 len_ipm=nbipm*nfxbody
1726 IF (irform/5<=1) THEN
1727 CALL redsqi (fxbipm,len_ipm,irform)
1728 IF (lennod>0) CALL redsqi (fxbnod,lennod,irform)
1729 IF (lenelm>0) CALL redsqi (fxbelm,lenelm,irform)
1730 IF (lengrvi>0) CALL redsqi (fxbgrvi,lengrvi,irform)
1731 ELSE
1732 CALL read_i_c(fxbipm,len_ipm)
1733 IF (lennod>0) CALL read_i_c(fxbnod,lennod)
1734 IF (lenelm>0) CALL read_i_c(fxbelm,lenelm)
1735 IF (lengrvi>0) CALL read_i_c(fxbgrvi,lengrvi)
1736 ENDIF
1737C--------------------------------------
1738C READING REALS
1739C--------------------------------------
1740 len_mod=lenmod*6
1741 IF (irform/5<=1) THEN
1742 IF (len_mod>0) CALL redsqr (fxbmod,len_mod,irform)
1743 IF (lenglm>0) CALL redsqr (fxbglm,lenglm,irform)
1744 IF (lencp>0) CALL redsqr (fxbcpm,lencp ,irform)
1745 IF (lencp>0) CALL redsqr (fxbcps,lencp ,irform)
1746 IF (lenlm>0) CALL redsqr (fxblm, lenlm ,irform)
1747 IF (lenfls>0) CALL redsqr (fxbfls,lenfls,irform)
1748 IF (lendls>0) CALL redsqr (fxbdls,lendls,irform)
1749 CALL redsqr (fxbdep,lenvar,irform)
1750 CALL redsqr (fxbvit,lenvar,irform)
1751 CALL redsqr (fxbacc,lenvar,irform)
1752 CALL redsqr (fxbrpm,lenrpm,irform)
1753 IF (lensig>0) CALL redsqr (fxbsig,lensig,irform)
1754 IF (lengrvr>0) CALL redsqr (fxbgrvr,lengrvr,irform)
1755 ELSE
1756 IF (len_mod>0) CALL read_db(fxbmod,len_mod)
1757 IF (lenglm>0) CALL read_db(fxbglm,lenglm)
1758 IF (lencp>0) CALL read_db(fxbcpm,lencp )
1759 IF (lencp>0) CALL read_db(fxbcps,lencp )
1760 IF (lenlm>0) CALL read_db(fxblm, lenlm )
1761 IF (lenfls>0) CALL read_db(fxbfls,lenfls)
1762 IF (lendls>0) CALL read_db(fxbdls,lendls)
1763 CALL read_db(fxbdep,lenvar)
1764 CALL read_db(fxbvit,lenvar)
1765 CALL read_db(fxbacc,lenvar)
1766 CALL read_db(fxbrpm,lenrpm)
1767 IF (lensig>0) CALL read_db(fxbsig,lensig)
1768 IF (lengrvr>0) CALL read_db(fxbgrvr,lengrvr)
1769 ENDIF
1770C reading mode and constraint files
1771 INQUIRE(iolength=rclen) vv
1772
1773 tlen=4096
1774 CALL tmpenvf(tmpdir,tlen)
1775 CALL my_getpid(my_pid)
1776 WRITE(cpid,'(I10.10)') my_pid
1777 WRITE(cispmd,'(I6.6)') ispmd
1778 ifxm_fn=tmpdir(1:tlen)//'/'//'25_'//rootn(1:lenrootn)//'_'//cpid//'_'//cispmd//'.tmp'
1779 ifxs_fn=tmpdir(1:tlen)//'/'//'26_'//rootn(1:lenrootn)//'_'//cpid//'_'//cispmd//'.tmp'
1780 len_ifxm_fn=len_trim(ifxm_fn)
1781 len_ifxs_fn=len_trim(ifxs_fn)
1782
1783 OPEN(unit=ifxm,file=trim(ifxm_fn),access='DIRECT',recl=rclen)
1784 OPEN(unit=ifxs,file=trim(ifxs_fn),access='DIRECT',recl=rclen)
1785
1786
1787 use_ifxm = 1
1788 use_ifxs = 1
1789 nrecm=0
1790 nrecs=0
1791 DO i=1,nfxbody
1792 ircm0=fxbipm(30,i)
1793 ircs0=fxbipm(31,i)
1794 ircm1=fxbipm(32,i)
1795 ircs1=fxbipm(33,i)
1796 nrecm=nrecm+ircm1-ircm0
1797 nrecs=nrecs+ircs1-ircs0
1798 ENDDO
1799 ircm=0
1800 ircs=0
1801 lrec=6
1802 DO i=1,nrecm
1803 ircm=ircm+1
1804 CALL read_db(vv,lrec)
1805 WRITE(ifxm,rec=ircm) (vv(j),j=1,lrec)
1806 ENDDO
1807 DO i=1,nrecs
1808 ircs=ircs+1
1809 CALL read_db(vv,lrec)
1810 WRITE(ifxs,rec=ircs) (vv(j),j=1,lrec)
1811 ENDDO
1812C
1813 RETURN
1814 END
1815!||====================================================================
1816!|| eigrrest ../engine/source/output/restart/rdresb.F
1817!||--- called by ------------------------------------------------------
1818!|| rdresb ../engine/source/output/restart/rdresb.F
1819!||--- calls -----------------------------------------------------
1820!|| my_getpid ../engine/source/system/pid.c
1821!|| read_db ../common_source/tools/input_output/read_db.F
1822!|| read_i_c ../common_source/tools/input_output/write_routines.c
1823!|| redsqi ../engine/source/output/tools/redsqi.F
1824!|| redsqr ../engine/source/output/tools/redsqr.F
1825!|| tmpenvf ../engine/source/system/tmpenv_c.c
1826!||--- uses -----------------------------------------------------
1827!|| eig_mod ../engine/share/modules/eig_mod.F
1828!|| scratchfile_mod ../engine/share/modules/scratchfiles_mod.f
1829!||====================================================================
1830 SUBROUTINE eigrrest()
1831C-----------------------------------------------
1832C M o d u l e s
1833C-----------------------------------------------
1834 USE eig_mod
1835 USE scratchfile_mod
1836C-----------------------------------------------
1837C I m p l i c i t T y p e s
1838C-----------------------------------------------
1839#include "implicit_f.inc"
1840C-----------------------------------------------
1841C C o m m o n B l o c k s
1842C-----------------------------------------------
1843#include "units_c.inc"
1844#include "com04_c.inc"
1845#include "scr05_c.inc"
1846#include "eigcom.inc"
1847#include "chara_c.inc"
1848#include "task_c.inc"
1849C-----------------------------------------------
1850C L o c a l V a r i a b l e s
1851C-----------------------------------------------
1852 INTEGER NRECM, IRCM, NBN, NBM, LREC, I, J, RCLEN,TLEN
1853 my_real
1854 . VV(6)
1855 CHARACTER(LEN=10) :: CPID
1856 CHARACTER(LEN=6) :: CISPMD
1857 CHARACTER(LEN=4096) :: TMPDIR
1858 INTEGER MY_PID
1859C--------------------------------------
1860C READING INTEGERS
1861C--------------------------------------
1862 IF (irform/5<=1) THEN
1863 CALL redsqi (eigipm,neig*neipm,irform)
1864 CALL redsqi (eigibuf,leibuf,irform)
1865 ELSE
1866 CALL read_i_c(eigipm,neig*neipm)
1867 CALL read_i_c(eigibuf,leibuf)
1868 ENDIF
1869C--------------------------------------
1870C READING REALS
1871C--------------------------------------
1872 IF (irform/5<=1) THEN
1873 CALL redsqr (eigrpm,neig*nerpm,irform)
1874 ELSE
1875 CALL read_db(eigrpm,neig*nerpm)
1876 ENDIF
1877C reading the additional modes file
1878 INQUIRE(iolength=rclen) vv
1879
1880 CALL my_getpid(my_pid)
1881 WRITE(cpid,'(I10.10)') my_pid
1882 WRITE(cispmd,'(I6.6)') ispmd
1883 tlen=4096
1884 CALL tmpenvf(tmpdir,tlen)
1885 ieigm_fn=tmpdir(1:tlen)//'/'//'27_'//rootn(1:lenrootn)//'_'//cpid//'_'//cispmd//'.tmp'
1886 len_ieigm_fn=len_trim(ieigm_fn)
1887 OPEN(unit=ieigm,file=trim(ieigm_fn),access='DIRECT',recl=rclen)
1888 use_ieigm = 1
1889
1890 nrecm=0
1891 DO i=1,neig
1892 nbn=eigipm(10,i)
1893 nbm=eigipm(14,i)
1894 nrecm=nrecm+nbn*nbm
1895 ENDDO
1896C
1897 ircm=0
1898 lrec=6
1899 DO i=1,nrecm
1900 ircm=ircm+1
1901 CALL read_db(vv,lrec)
1902 WRITE(ieigm,rec=ircm) (vv(j),j=1,lrec)
1903 ENDDO
1904C
1905 RETURN
1906 END
1907
1908!||====================================================================
1909!|| dsrrest ../engine/source/output/restart/rdresb.f
1910!||--- calls -----------------------------------------------------
1911!|| read_db ../common_source/tools/input_output/read_db.F
1912!|| read_i_c ../common_source/tools/input_output/write_routines.c
1913!|| redsqi ../engine/source/output/tools/redsqi.F
1914!|| redsqr ../engine/source/output/tools/redsqr.F
1915!||--- uses -----------------------------------------------------
1916!|| dsgraph_mod ../engine/share/modules/dsgraph_mod.F
1917!||====================================================================
1918 SUBROUTINE dsrrest(GRAPHE)
1919C-----------------------------------------------
1920C M o d u l e s
1921C-----------------------------------------------
1922 USE dsgraph_mod
1923C-----------------------------------------------
1924C I m p l i c i t T y p e s
1925C-----------------------------------------------
1926#include "implicit_f.inc"
1927C-----------------------------------------------
1928C C o m m o n B l o c k s
1929C-----------------------------------------------
1930#include "com01_c.inc"
1931#include "scr05_c.inc"
1932C-----------------------------------------------
1933C D u m m y A r g u m e n t s
1934C-----------------------------------------------
1935 TYPE(prgraph) :: GRAPHE(*)
1936C-----------------------------------------------
1937C L o c a l V a r i a b l e s
1938C-----------------------------------------------
1939 INTEGER I, J, TN(7), N, NDDL, DSNDDL, DSLEN, K, L, NSDEC
1940 my_real
1941 . CUTFREQ
1942C
1943 DO i=1,nsproc
1944C--------------------------------------
1945C reading integer parameters
1946C--------------------------------------
1947 n=7
1948 IF (irform/5<=1) THEN
1949 CALL redsqi(tn, n, irform)
1950 ELSE
1951 CALL read_i_c(tn, n)
1952 ENDIF
1953 dsnddl=tn(1)
1954 dslen=tn(2)
1955 nsdec=tn(4)
1956 graphe(i)%NDDL=dsnddl
1957 graphe(i)%NSUPEL=dslen
1958 graphe(i)%NSLEVEL=tn(3)
1959 graphe(i)%NSDEC=nsdec
1960 graphe(i)%NSVMAX=tn(5)
1961 graphe(i)%IPRI=tn(6)
1962 graphe(i)%NDDL_GLOB=tn(7)
1963C--------------------------------------
1964C READING REALS
1965C--------------------------------------
1966 n=1
1967 IF (irform/5<=1) THEN
1968 CALL redsqr(cutfreq, n, irform)
1969 ELSE
1970 CALL read_db(cutfreq, n)
1971 ENDIF
1972 graphe(i)%CUTFREQ=cutfreq
1973C--------------------------------------
1974C READING ARRAY DOF
1975C--------------------------------------
1976 ALLOCATE(graphe(i)%LSDDL(2,dsnddl),
1977 . graphe(i)%LSDDL_INI(2,dsnddl),
1978 . graphe(i)%LSDDL_GLOB(dsnddl),
1979 . graphe(i)%LSDDL_GLOB_INI(dsnddl))
1980 DO j=1,2
1981 IF (irform/5<=1) THEN
1982 CALL redsqi(graphe(i)%LSDDL(j,1:dsnddl), dsnddl, irform)
1983 ELSE
1984 CALL read_i_c(graphe(i)%LSDDL(j,1:dsnddl), dsnddl)
1985 ENDIF
1986 ENDDO
1987 IF (irform/5<=1) THEN
1988 CALL redsqi(graphe(i)%LSDDL_GLOB, dsnddl, irform)
1989 ELSE
1990 CALL read_i_c(graphe(i)%LSDDL_GLOB, dsnddl)
1991 ENDIF
1992 graphe(i)%NDDL_INI=graphe(i)%NDDL
1993 DO j=1,2
1994 DO k=1,dsnddl
1995 graphe(i)%LSDDL_INI(j,k)=graphe(i)%LSDDL(j,k)
1996 ENDDO
1997 ENDDO
1998 DO j=1,dsnddl
1999 graphe(i)%LSDDL_GLOB_INI(j)=graphe(i)%LSDDL_GLOB(j)
2000 ENDDO
2001C--------------------------------------
2002C READING SUPERELEMENTS
2003C--------------------------------------
2004 ALLOCATE(graphe(i)%DGRAPH(dslen))
2005 DO j=1,dslen
2006 n=3
2007 IF (irform/5<=1) THEN
2008 CALL redsqi(tn, n, irform)
2009 ELSE
2010 CALL read_i_c(tn, n)
2011 ENDIF
2012 graphe(i)%DGRAPH(j)%NDDL_I=tn(1)
2013 graphe(i)%DGRAPH(j)%NDDL_F=tn(2)
2014 graphe(i)%DGRAPH(j)%NSDMAX=tn(3)
2015 nddl=tn(1)+tn(2)
2016 ALLOCATE(graphe(i)%DGRAPH(j)%CHILD(nsdec))
2017 ALLOCATE(graphe(i)%DGRAPH(j)%DDLS(nddl),
2018 . graphe(i)%DGRAPH(j)%DDLS_INI(nddl))
2019 ALLOCATE(graphe(i)%DGRAPH(j)%IFAC(tn(3)+1,tn(2)),
2020 . graphe(i)%DGRAPH(j)%IFACM(tn(2)))
2021 ALLOCATE(graphe(i)%DGRAPH(j)%IFAC_INI(tn(3)+1,tn(2)),
2022 . graphe(i)%DGRAPH(j)%IFACM_INI(tn(2)))
2023 n=tn(3)+1
2024 IF (irform/5<=1) THEN
2025 CALL redsqi(graphe(i)%DGRAPH(j)%CHILD, nsdec, irform)
2026 CALL redsqi(graphe(i)%DGRAPH(j)%DDLS, nddl, irform)
2027 DO k=1,tn(2)
2028 CALL redsqi(graphe(i)%DGRAPH(j)
2029 . %IFAC(1:n,k), n, irform)
2030 ENDDO
2031 CALL redsqi(graphe(i)%DGRAPH(j)%IFACM, tn(2), irform)
2032 ELSE
2033 CALL read_i_c(graphe(i)%DGRAPH(j)%CHILD, nsdec)
2034 CALL read_i_c(graphe(i)%DGRAPH(j)%DDLS, nddl)
2035 DO k=1,tn(2)
2036 CALL read_i_c(graphe(i)%DGRAPH(j)
2037 . %IFAC(1:n,k), n)
2038 ENDDO
2039 CALL read_i_c(graphe(i)%DGRAPH(j)%IFACM, tn(2))
2040 ENDIF
2041 graphe(i)%DGRAPH(j)%NDDL_I_INI=graphe(i)%DGRAPH(j)%NDDL_I
2042 graphe(i)%DGRAPH(j)%NDDL_F_INI=graphe(i)%DGRAPH(j)%NDDL_F
2043 DO k=1,nddl
2044 graphe(i)%DGRAPH(j)%DDLS_INI(k)=
2045 . graphe(i)%DGRAPH(j)%DDLS(k)
2046 ENDDO
2047 DO k=1,tn(2)
2048 DO l=1,tn(3)+1
2049 graphe(i)%DGRAPH(j)%IFAC_INI(l,k)=
2050 . graphe(i)%DGRAPH(j)%IFAC(l,k)
2051 ENDDO
2052 graphe(i)%DGRAPH(j)%IFACM_INI(k)=
2053 . graphe(i)%DGRAPH(j)%IFACM(k)
2054 ENDDO
2055 ENDDO
2056C
2057 ENDDO
2058C-----------
2059 RETURN
2060 END
2061
2062!||====================================================================
2063!|| nfrrest ../engine/source/output/restart/rdresb.F
2064!||--- called by ------------------------------------------------------
2065!|| rdresb ../engine/source/output/restart/rdresb.F
2066!||--- calls -----------------------------------------------------
2067!|| read_db ../common_source/tools/input_output/read_db.F
2068!|| read_i_c ../common_source/tools/input_output/write_routines.c
2069!|| redsqi ../engine/source/output/tools/redsqi.F
2070!|| redsqr ../engine/source/output/tools/redsqr.F
2071!||====================================================================
2072 SUBROUTINE nfrrest(IFLOW, RFLOW)
2073C-----------------------------------------------
2074C I m p l i c i t T y p e s
2075C-----------------------------------------------
2076#include "implicit_f.inc"
2077C-----------------------------------------------
2078C C o m m o n B l o c k s
2079C-----------------------------------------------
2080#include "scr05_c.inc"
2081#include "flowcom.inc"
2082C-----------------------------------------------
2083C D u m m y A r g u m e n t s
2084C-----------------------------------------------
2085 INTEGER IFLOW(*)
2086 my_real
2087 . RFLOW(*)
2088C-----------------------------------------------
2089C L o c a l V a r i a b l e s
2090C-----------------------------------------------
2091C
2092C--------------------------------------
2093C READING INTEGERS
2094C--------------------------------------
2095 IF (irform/5<=1) THEN
2096 CALL redsqi (iflow,liflow,irform)
2097 ELSE
2098 CALL read_i_c(iflow,liflow)
2099 ENDIF
2100C--------------------------------------
2101C READING REALS
2102C--------------------------------------
2103 IF (irform/5<=1) THEN
2104 CALL redsqr (rflow,lrflow,irform)
2105 ELSE
2106 CALL read_db(rflow,lrflow)
2107 ENDIF
2108C
2109 RETURN
2110 END
2111
2112!||====================================================================
2113!|| imprrest ../engine/source/output/restart/rdresb.F
2114!||--- called by ------------------------------------------------------
2115!|| rdresb ../engine/source/output/restart/rdresb.F
2116!||--- calls -----------------------------------------------------
2117!|| cp_real ../engine/source/implicit/produt_v.F
2118!|| read_db ../common_source/tools/input_output/read_db.F
2119!||--- uses -----------------------------------------------------
2120!|| imp_dyna ../engine/share/modules/impbufdef_mod.F
2121!|| imp_rest ../engine/share/modules/impbufdef_mod.F
2122!||====================================================================
2123 SUBROUTINE imprrest(NIMPR)
2124C-----------------------------------------------
2125C M o d u l e s
2126C-----------------------------------------------
2127 USE imp_rest
2128 USE imp_dyna
2129C-----------------------------------------------
2130C I m p l i c i t T y p e s
2131C-----------------------------------------------
2132#include "implicit_f.inc"
2133C-----------------------------------------------
2134C C o m m o n B l o c k s
2135C-----------------------------------------------
2136#include "com01_c.inc"
2137#include "com04_c.inc"
2138#include "impl1_c.inc"
2139C-----------------------------------------------
2140C D u m m y A r g u m e n t s
2141C-----------------------------------------------
2142 INTEGER NIMPR
2143C-----------------------------------------------
2144C L o c a l V a r i a b l e s
2145C-----------------------------------------------
2146 INTEGER ND,IDYNA0
2147 my_real
2148 . A(3*NUMNOD),AR(3*NUMNOD)
2149C--------------------------------------
2150 IF(ALLOCATED(imp_rr)) DEALLOCATE(imp_rr)
2151 ALLOCATE(imp_rr(nimpr))
2152C--------------------------------------
2153C READING REALS
2154C--------------------------------------
2155 nd = numnod*3
2156 idyna0 = nimpr - 2
2157
2158 CALL read_db(imp_rr,nimpr)
2159
2160 IF (idyna0>0) THEN
2161 CALL read_db(a,nd)
2162 IF (iroddl/=0) CALL read_db(ar,nd)
2163 ENDIF
2164
2165C
2166 IF ((idyna+idyna0)>0) THEN
2167 ALLOCATE(dy_a(3,numnod))
2168 IF (iroddl/=0) ALLOCATE(dy_ar(3,numnod))
2169 CALL cp_real(nd,a,dy_a)
2170 IF (iroddl/=0) CALL cp_real(nd,ar,dy_ar)
2171 ENDIF
2172C
2173 RETURN
2174 END
2175
2176!||====================================================================
2177!|| fvrrest ../engine/source/output/restart/rdresb.F
2178!||--- called by ------------------------------------------------------
2179!|| rdresb ../engine/source/output/restart/rdresb.F
2180!||--- calls -----------------------------------------------------
2181!|| ancmsg ../engine/source/output/message/message.F
2182!|| arret ../engine/source/system/arret.F
2183!|| read_db ../common_source/tools/input_output/read_db.F
2184!|| read_i_c ../common_source/tools/input_output/write_routines.c
2185!|| spmd_comm_split ../engine/source/mpi/generic/spmd_comm_split.F
2186!|| spmd_fvb_comm_pattern ../engine/source/mpi/airbags/spmd_fvb_comm_pattern.F
2187!||--- uses -----------------------------------------------------
2188!|| fvbag_mod ../engine/share/modules/fvbag_mod.F
2189!|| fvmbag_meshcontrol_mod ../common_source/modules/airbag/fvmbag_meshcontrol_mod.F
2190!|| message_mod ../engine/share/message_module/message_mod.F
2191!||====================================================================
2192 SUBROUTINE fvrrest(MONVOL)
2193C-----------------------------------------------
2194C M o d u l e s
2195C-----------------------------------------------
2196 USE fvbag_mod
2197 USE fvmbag_meshcontrol_mod
2198 USE message_mod
2199C-----------------------------------------------
2200C I m p l i c i t T y p e s
2201C-----------------------------------------------
2202#include "implicit_f.inc"
2203C-----------------------------------------------
2204C C o m m o n B l o c k s
2205C-----------------------------------------------
2206#include "com01_c.inc"
2207#include "com04_c.inc"
2208#include "task_c.inc"
2209#include "param_c.inc"
2210#include "tabsiz_c.inc"
2211#include "units_c.inc"
2212C-----------------------------------------------
2213C D u m m y A r g u m e n t s
2214C-----------------------------------------------
2215 INTEGER, INTENT(IN) :: MONVOL(SMONVOL)
2216C-----------------------------------------------
2217C L o c a l V a r i a b l e s
2218C-----------------------------------------------
2219 INTEGER LEN, I, II, KK, TABI(11), NNS, NNTR, LENP, NPOLY, LENH, NPOLH,
2220 . LENP_ANIM, NPOLH_ANIM, LENH_ANIM, NNS_ANIM, NN_L, NNA_L,
2221 . NNSA_L, NSA, NELSA, NNI_L, IFV, ID_DT_OPTION, K1, ITYP,
2222 . ID_BAG, ID_BAG_INPUT, N
2223 INTEGER COLOR,KEY, IH3D_FLAG
2224 LOGICAL lFOUND
2225C-----------------------------------------------
2226C S o u r c e L i n e s
2227C-----------------------------------------------
2228
2229C----------------------------------------------------------
2230C CONSTANT INIT. : VISUALIZATION OF FVBAG RESULTS
2231C----------------------------------------------------------
2232 !---AIRBAGS_TOTAL_FVM_IN_H3D
2233 airbags_total_fvm_in_h3d = 0 !number of polyhedra to display for visualization subroutines (Sum FVDATA%NPOLH_ANIM)
2234
2235C--------------------------------------
2236C READING SIZES
2237C--------------------------------------
2238 len=1
2239 CALL read_i_c(nfvbag,len)
2240 ALLOCATE(fvdata_old(nfvbag))
2241 DO i=1,nfvbag
2242 NULLIFY(fvdata_old(i)%BRIC,fvdata_old(i)%TBRIC,
2243 * fvdata_old(i)%XB,fvdata_old(i)%SFAC,
2244 * fvdata_old(i)%IFVPOLY_ANIM,fvdata_old(i)%IFVTADR_ANIM,
2245 * fvdata_old(i)%IFVPOLH_ANIM,fvdata_old(i)%IFVPADR_ANIM,
2246 * fvdata_old(i)%REDIR_ANIM,fvdata_old(i)%IFVTRI_ANIM,
2247 * fvdata_old(i)%NOD_ANIM,fvdata_old(i)%IFVPOLY,
2248 * fvdata_old(i)%IFVTADR,fvdata_old(i)%IFVPOLH,
2249 * fvdata_old(i)%IFVPADR,fvdata_old(i)%IDPOLH,
2250 * fvdata_old(i)%IBPOLH,fvdata_old(i)%IFVNOD,
2251 * fvdata_old(i)%IFVTRI,fvdata_old(i)%RFVNOD,
2252 * fvdata_old(i)%MPOLH,fvdata_old(i)%EPOLH,
2253 * fvdata_old(i)%PPOLH,fvdata_old(i)%RPOLH,
2254 * fvdata_old(i)%GPOLH,fvdata_old(i)%TPOLH,
2255 * fvdata_old(i)%CPAPOLH,
2256 * fvdata_old(i)%CPBPOLH,fvdata_old(i)%CPCPOLH,
2257 * fvdata_old(i)%CPDPOLH,fvdata_old(i)%CPEPOLH,
2258 * fvdata_old(i)%CPFPOLH,fvdata_old(i)%DTPOLH,
2259 * fvdata_old(i)%RMWPOLH,fvdata_old(i)%VPOLH_INI,
2260 * fvdata_old(i)%QPOLH,
2261 * fvdata_old(i)%SSPPOLH,
2262 * fvdata_old(i)%CENTROID_POLH )
2263
2264 fvdata_old(i)%L_TYPE = 0
2265 fvdata_old(i)%ID_DT_OPTION = -1
2266 fvdata_old(i)%LAMBDA = zero
2267 fvdata_old(i)%DTOLD = zero
2268 fvdata_old(i)%CFL_COEF = zero
2269 fvdata_old(i)%DTMIN = zero
2270 fvdata_old(i)%PDISP_OLD = zero
2271 fvdata_old(i)%PDISP = zero
2272 ENDDO
2273
2274 ALLOCATE(fvspmd(nfvbag))
2275 DO i=1,nfvbag
2276 NULLIFY(fvspmd(i)%IBUF_L)
2277 NULLIFY(fvspmd(i)%IBUFA_L)
2278 NULLIFY(fvspmd(i)%IBUFSA_L)
2279 NULLIFY(fvspmd(i)%IXSA)
2280 NULLIFY(fvspmd(i)%ELEMSA)
2281 NULLIFY(fvspmd(i)%ITAB)
2282 ENDDO
2283C
2284
2285 ALLOCATE(kmesh(nfvbag))
2286 DO i=1,nfvbag
2287
2288 CALL read_i_c(fvdata_old(i)%L_TYPE,1)
2289 CALL read_i_c(fvdata_old(i)%ID_DT_OPTION,1)
2290 CALL read_db(fvdata_old(i)%LAMBDA,1)
2291 CALL read_db(fvdata_old(i)%DTOLD,1)
2292 CALL read_db(fvdata_old(i)%CFL_COEF,1)
2293 CALL read_db(fvdata_old(i)%DTMIN,1)
2294
2295 len=8
2296 CALL read_i_c(tabi,len)
2297 nn_l=tabi(1)
2298 nni_l=tabi(8)
2299 nna_l=tabi(2)
2300 nnsa_l=tabi(3)
2301 nsa=tabi(4)
2302 nelsa=tabi(6)
2303 fvspmd(i)%NN_L=nn_l
2304 fvspmd(i)%NNI_L=nni_l
2305 fvspmd(i)%NNA_L=nna_l
2306 fvspmd(i)%NNSA_L=nnsa_l
2307 fvspmd(i)%NSA=nsa
2308 fvspmd(i)%NNSA=tabi(5)
2309 fvspmd(i)%NELSA=nelsa
2310 fvspmd(i)%PMAIN=tabi(7)
2311 fvspmd(i)%NNA_L_GLOB = nna_l
2312
2313 color = min(1,nn_l+nni_l+nna_l+nnsa_l)
2314 IF(ispmd == fvspmd(i)%PMAIN - 1) THEN
2315 ! PMAIN will be zero in the subcommunicator
2316 key = 0
2317 color = 1
2318 ELSE
2319 key = 1
2320 ENDIF
2321
2322
2323 ! COLOR == 0 => ISPMD does not belong to the fvm
2324 ! COLOR == 1 => ISPMD belongs to the fvm
2325 fvspmd(i)%MPI_COMM = 0
2326 fvspmd(i)%RANK = 0
2327 fvspmd(i)%NSPMD = 1
2328 CALL spmd_comm_split(color,key,fvspmd(i)%MPI_COMM,
2329 . fvspmd(i)%RANK, fvspmd(i)%NSPMD)
2330
2331
2332 ALLOCATE(fvspmd(i)%ITAB(4,nspmd-1))
2333
2334 ! Compute the size of the comm. involved in SPMD_FVB_GATHER
2335 CALL spmd_fvb_comm_pattern(i)
2336
2337 ALLOCATE(fvspmd(i)%IBUF_L(2,nn_l+nni_l))
2338 len=2*(nn_l+nni_l)
2339 IF (len>0) THEN
2340 CALL read_i_c(fvspmd(i)%IBUF_L,len)
2341 ENDIF
2342 ALLOCATE(fvspmd(i)%IBUFA_L(2,nna_l))
2343 len=2*nna_l
2344 IF (len>0) THEN
2345 CALL read_i_c(fvspmd(i)%IBUFA_L,len)
2346 ENDIF
2347 ALLOCATE(fvspmd(i)%IBUFSA_L(2,nnsa_l))
2348 len=2*nnsa_l
2349 IF (len>0) THEN
2350 CALL read_i_c(fvspmd(i)%IBUFSA_L,len)
2351 ENDIF
2352 len = 1
2353 CALL read_i_c(kmesh(i), len)
2354 IF (ispmd/=fvspmd(i)%PMAIN-1) cycle
2355C
2356 ALLOCATE(fvspmd(i)%IXSA(8,nsa))
2357 len=8*nsa
2358 IF (len>0) THEN
2359 CALL read_i_c(fvspmd(i)%IXSA,len)
2360 ENDIF
2361 ALLOCATE(fvspmd(i)%ELEMSA(3,nelsa))
2362 len=3*nelsa
2363 IF (len>0) THEN
2364 CALL read_i_c(fvspmd(i)%ELEMSA,len)
2365 ENDIF
2366 len=11
2367 CALL read_i_c(tabi,len)
2368C
2369 nns=tabi(1)
2370 nntr=tabi(2)
2371 lenp=tabi(3)
2372 npoly=tabi(4)
2373 lenh=tabi(5)
2374 npolh=tabi(6)
2375 lenp_anim=tabi(7)
2376 lenh_anim=tabi(8)
2377 npolh_anim=tabi(9)
2378 nns_anim=tabi(10)
2379C
2380 fvdata_old(i)%NNS=nns
2381 fvdata_old(i)%NNTR=nntr
2382 fvdata_old(i)%LENP=lenp
2383 fvdata_old(i)%NPOLY=npoly
2384 fvdata_old(i)%LENH=lenh
2385 fvdata_old(i)%NPOLH=npolh
2386 fvdata_old(i)%NPOLH_ANIM=npolh_anim
2387 fvdata_old(i)%NNS_ANIM=nns_anim
2388 fvdata_old(i)%ID=tabi(11)
2389C
2390 ALLOCATE(fvdata_old(i)%IFVNOD(3,nns),
2391 . fvdata_old(i)%RFVNOD(2,nns),
2392 . fvdata_old(i)%IFVTRI(6,nntr),
2393 . fvdata_old(i)%IFVPOLY(lenp),
2394 . fvdata_old(i)%IFVTADR(npoly+1),
2395 . fvdata_old(i)%IFVPOLH(lenh),
2396 . fvdata_old(i)%IFVPADR(npolh+1),
2397 . fvdata_old(i)%IDPOLH(npolh),
2398 . fvdata_old(i)%IBPOLH(npolh),
2399 . fvdata_old(i)%MPOLH(npolh),
2400 . fvdata_old(i)%QPOLH(3,npolh),
2401 . fvdata_old(i)%EPOLH(npolh),
2402 . fvdata_old(i)%PPOLH(npolh),
2403 . fvdata_old(i)%RPOLH(npolh),
2404 . fvdata_old(i)%GPOLH(npolh),
2405 . fvdata_old(i)%TPOLH(npolh),
2406 . fvdata_old(i)%CPAPOLH(npolh),
2407 . fvdata_old(i)%CPBPOLH(npolh),
2408 . fvdata_old(i)%CPCPOLH(npolh),
2409 . fvdata_old(i)%CPDPOLH(npolh),
2410 . fvdata_old(i)%CPEPOLH(npolh),
2411 . fvdata_old(i)%CPFPOLH(npolh),
2412 . fvdata_old(i)%RMWPOLH(npolh),
2413 . fvdata_old(i)%VPOLH_INI(npolh),
2414 . fvdata_old(i)%DTPOLH(npolh),
2415 . fvdata_old(i)%SSPPOLH(npolh),
2416 . fvdata_old(i)%CENTROID_POLH(3,npolh) )
2417
2418 IF (npolh_anim>0)
2419 . ALLOCATE(fvdata_old(i)%IFVPOLY_ANIM(lenp_anim),
2420 . fvdata_old(i)%IFVTADR_ANIM(npoly+1),
2421 . fvdata_old(i)%IFVPOLH_ANIM(lenh_anim),
2422 . fvdata_old(i)%IFVPADR_ANIM(npolh_anim+1),
2423 . fvdata_old(i)%IFVTRI_ANIM(6,nntr),
2424 . fvdata_old(i)%REDIR_ANIM(nns_anim),
2425 . fvdata_old(i)%NOD_ANIM(3,nns_anim) )
2426
2427
2428C--------------------------------------
2429C READING INTEGERS
2430C--------------------------------------
2431 len=3*nns
2432 CALL read_i_c(fvdata_old(i)%IFVNOD,len)
2433 len=6*nntr
2434 CALL read_i_c(fvdata_old(i)%IFVTRI,len)
2435 CALL read_i_c(fvdata_old(i)%IFVPOLY,lenp)
2436 CALL read_i_c(fvdata_old(i)%IFVTADR,npoly+1)
2437 CALL read_i_c(fvdata_old(i)%IFVPOLH,lenh)
2438 CALL read_i_c(fvdata_old(i)%IFVPADR,npolh+1)
2439 CALL read_i_c(fvdata_old(i)%IDPOLH,npolh)
2440 CALL read_i_c(fvdata_old(i)%IBPOLH,npolh)
2441 IF (npolh_anim>0) THEN
2442 CALL read_i_c(fvdata_old(i)%IFVPOLY_ANIM,lenp_anim)
2443 CALL read_i_c(fvdata_old(i)%IFVTADR_ANIM,npoly+1)
2444 CALL read_i_c(fvdata_old(i)%IFVPOLH_ANIM,lenh_anim)
2445 CALL read_i_c(fvdata_old(i)%IFVPADR_ANIM,npolh_anim+1)
2446 len=6*nntr
2447 CALL read_i_c(fvdata_old(i)%IFVTRI_ANIM,len)
2448 CALL read_i_c(fvdata_old(i)%REDIR_ANIM,nns_anim)
2449 ENDIF
2450C--------------------------------------
2451C READING REALS
2452C--------------------------------------
2453 len = 1
2454 CALL read_db(fvdata_old(i)%PDISP_OLD,len)
2455 CALL read_db(fvdata_old(i)%PDISP,len)
2456 len=2*nns
2457 CALL read_db(fvdata_old(i)%RFVNOD,len)
2458 CALL read_db(fvdata_old(i)%MPOLH,npolh)
2459 len=3*npolh
2460 CALL read_db(fvdata_old(i)%QPOLH,len)
2461 CALL read_db(fvdata_old(i)%EPOLH,npolh)
2462 CALL read_db(fvdata_old(i)%PPOLH,npolh)
2463 CALL read_db(fvdata_old(i)%RPOLH,npolh)
2464 CALL read_db(fvdata_old(i)%GPOLH,npolh)
2465 CALL read_db(fvdata_old(i)%TPOLH,npolh)
2466 CALL read_db(fvdata_old(i)%CPAPOLH,npolh)
2467 CALL read_db(fvdata_old(i)%CPBPOLH,npolh)
2468 CALL read_db(fvdata_old(i)%CPCPOLH,npolh)
2469 CALL read_db(fvdata_old(i)%CPDPOLH,npolh)
2470 CALL read_db(fvdata_old(i)%CPEPOLH,npolh)
2471 CALL read_db(fvdata_old(i)%CPFPOLH,npolh)
2472 CALL read_db(fvdata_old(i)%RMWPOLH,npolh)
2473 CALL read_db(fvdata_old(i)%VPOLH_INI,npolh)
2474 CALL read_db(fvdata_old(i)%DLH,1)
2475 CALL read_db(fvdata_old(i)%DTPOLH,npolh)
2476 IF (npolh_anim>0) THEN
2477 len=3*nns_anim
2478 CALL read_db(fvdata_old(i)%NOD_ANIM,len)
2479 ENDIF
2480 !FVDATA_OLD(I)%SSPPOLH not read : H3D only. Generated if needed.
2481 !FVDATA_OLD(I)%CENTROID_POLH not read : H3D only. Generated if needed.
2482C
2483 k1 = 1
2484 ih3d_flag = 0
2485 !Looping over all existing bags to check which one is matching : collect IHED_FLAG
2486 DO n=1,nvolu
2487 ifv = monvol(k1-1 + 45)
2488 IF(ifv==i)THEN
2489 ih3d_flag = monvol(k1-1 + 75)
2490 IF(ih3d_flag==1)THEN
2491 airbags_total_fvm_in_h3d = airbags_total_fvm_in_h3d + npolh_anim
2492 ENDIF
2493 EXIT
2494 ENDIF
2495 k1=k1+nimv
2496 ENDDO !next n
2497
2498C
2499 ENDDO !next FVBAG
2500
2501
2502C--------------------------------------
2503C /MONVOL/FVMBAG1(typ=8)
2504C /DT/FVMBAG/1 if relevant for ITYP8 only
2505C--------------------------------------
2506C FVMBAG_INPUT_OPTIONS is a data structure from FREFORM (Engine Reader) which contains data from /DT/FVMBAG options
2507C these parameters are treated once all options are read and after initialization of parameters from restart file
2508 !looping over all /DT/FVMBAG options
2509 DO ii=num_opt_dt_fvmbag,1,-1
2510 id_dt_option = fvmbag_input_options(ii)%ID_DT_OPTION
2511 ! /DT/FVMBAG/2 otherwise next cycle
2512 IF(id_dt_option==1)THEN
2513 k1=1
2514 !Looping over all existing bags to check if one id is matching
2515 DO n=1,nvolu
2516 id_bag=monvol(k1-1+1)
2517 ityp=monvol(k1-1+2)
2518 IF (ityp == 8) THEN
2519 ifv = monvol(k1 -1 +45)
2520 !IF(FVDATA_OLD(IFV)%ID_DT_OPTION == 1)THEN
2521 fvdata_old(ifv)%L_TYPE = -1
2522 fvdata_old(ifv)%ID_DT_OPTION = 1
2523 fvdata_old(ifv)%LAMBDA = zero
2524 fvdata_old(ifv)%CFL_COEF = fvmbag_input_options(ii)%CFL_COEF
2525 fvdata_old(ifv)%DTMIN = fvmbag_input_options(ii)%DTMIN
2526 !ENDIF
2527 ENDIF
2528 k1=k1+nimv
2529 ENDDO !next N
2530 ENDIF !ID_DT_OPTION
2531 enddo!next II
2532
2533C--------------------------------------
2534C /MONVOL/FVMBAG1(typ=8) & /MONVOL/FVMBAG2(typ=11)
2535C UPDATING DATA WITH ENGINE FILE PARAMETERS(/DT/FVMBAG/2 )
2536C--------------------------------------
2537 ! FILL FVDATA BUFFER (NOW ALLOCATED) WITH PARAMETERS READ IN ENGINE FILE (/DT/FVMBAG OPTION)
2538 !looping over all /DT/FVMBAG options
2539 DO ii=1,num_opt_dt_fvmbag
2540 id_dt_option = fvmbag_input_options(ii)%ID_DT_OPTION
2541 ! /DT/FVMBAG/2 otherwise next cycle
2542 IF(id_dt_option == 2)THEN
2543 k1=1
2544 id_bag_input = fvmbag_input_options(ii)%ID_BAG
2545 lfound = .false.
2546 IF(id_bag_input > 0)THEN
2547 !Looping over all existing bags to check if one id is matching
2548 DO n=1,nvolu
2549 id_bag=monvol(k1-1+1)
2550 ityp=monvol(k1-1+2)
2551 IF(id_bag == id_bag_input)THEN
2552 ifv = monvol(k1 -1 +45)
2553 IF (ityp == 8 .OR. ityp == 11) THEN
2554 lfound=.true.
2555 fvdata_old(ifv)%L_TYPE = fvmbag_input_options(ii)%L_TYPE
2556 fvdata_old(ifv)%ID_DT_OPTION = fvmbag_input_options(ii)%ID_DT_OPTION
2557 fvdata_old(ifv)%LAMBDA = fvmbag_input_options(ii)%LAMBDA
2558 fvdata_old(ifv)%CFL_COEF = fvmbag_input_options(ii)%CFL_COEF
2559 fvdata_old(ifv)%DTMIN = fvmbag_input_options(ii)%DTMIN
2560 ENDIF
2561 ENDIF
2562 k1=k1+nimv
2563 IF(lfound)EXIT
2564 ENDDO !next N
2565 IF(.NOT.lfound)THEN
2566 CALL ancmsg(msgid=299,anmode=aninfo,i1=id_bag_input)
2567 CALL arret(2)
2568 ENDIF
2569 ENDIF !(ID_BAG_INPUT > 0)
2570 ENDIF !IF(ID_DT_OPTION == 2)
2571 enddo!next II
2572
2573
2574
2575
2576
2577C--------------------------------------
2578C /MONVOL/FVMBAG1(typ=6) & /MONVOL/FVMBAG1(typ=8) & /MONVOL/FVMBAG2(typ=11)
2579C DEFAULT
2580C--------------------------------------
2581 k1=1
2582 DO n=1,nvolu
2583 id_bag=monvol(k1-1+1)
2584 ityp=monvol(k1-1+2)
2585 IF (ityp == 6 .OR. ityp == 8 .OR. ityp == 11) THEN
2586 ifv = monvol(k1 -1 +45)
2587 IF(fvdata_old(ifv)%CFL_COEF == zero)THEN
2588 fvdata_old(ifv)%CFL_COEF = zep9
2589 ENDIF
2590 IF(fvdata_old(ifv)%DTMIN == zero)THEN
2591 fvdata_old(ifv)%DTMIN = em20
2592 ENDIF
2593 IF(fvdata_old(ifv)%LAMBDA == zero)THEN
2594 fvdata_old(ifv)%LAMBDA = em02
2595 ENDIF
2596 IF(fvdata_old(ifv)%L_TYPE == 0)THEN
2597 fvdata_old(ifv)%L_TYPE = 1
2598 ENDIF
2599 ENDIF
2600 k1=k1+nimv
2601 ENDDO !next N
2602
2603C--------------------------------------
2604C CHECK INCOMPATIBLE TIME STEP OPTION /DT/FVMBAG
2605C before deallocating FVMBAG_INPUT_OPTIONS
2606C--------------------------------------
2607 !Warn user if incompatible options are defined
2608 IF( num_opt_dt_fvmbag_0+num_opt_dt_fvmbag_1 > 0 .AND. num_opt_dt_fvmbag_2 > 0)THEN
2609 ! if /DT/FVMBAG/2 is used on a given airbag, then general option /DT/FVMBAG/1 is ignored for this airbag.
2610 DO kk=1,num_opt_dt_fvmbag
2611 !display WARNING DUPLICATED DEFINITION
2612 IF(fvmbag_input_options(kk)%ID_DT_OPTION == 2)THEN
2613 id_bag = fvmbag_input_options(kk)%ID_BAG
2614 CALL ancmsg(msgid=302,anmode=aninfo,i1=id_bag)
2615 ENDIF
2616 ENDDO
2617 ENDIF
2618
2619 !Warn user if incompatible options are defined
2620 IF( num_opt_dt_fvmbag_1 > 1 )THEN
2621 ! if /DT/FVMBAG/2 is used on a given airbag, then general option /DT/FVMBAG/1 is ignored for this airbag.
2622 CALL ancmsg(msgid=303,anmode=aninfo,i1=id_bag)
2623 CALL arret(2)
2624 ENDIF
2625
2626C--------------------------------------
2627 IF(ALLOCATED(fvmbag_input_options))DEALLOCATE(fvmbag_input_options)
2628C--------------------------------------
2629
2630
2631 RETURN
2632 END
2633
2634C shell heat transfer
2635!||====================================================================
2636!|| thcrrest ../engine/source/output/restart/rdresb.F
2637!||--- called by ------------------------------------------------------
2638!|| rdresb ../engine/source/output/restart/rdresb.F
2639!||--- calls -----------------------------------------------------
2640!|| read_db ../common_source/tools/input_output/read_db.F
2641!|| redsqr ../engine/source/output/tools/redsqr.F
2642!||====================================================================
2643 SUBROUTINE thcrrest(MCP, TEMP)
2644C-----------------------------------------------
2645C I m p l i c i t T y p e s
2646C-----------------------------------------------
2647#include "implicit_f.inc"
2648C-----------------------------------------------
2649C C o m m o n B l o c k s
2650C-----------------------------------------------
2651#include "com04_c.inc"
2652#include "scr05_c.inc"
2653C-----------------------------------------------
2654C D u m m y A r g u m e n t s
2655C-----------------------------------------------
2656 my_real
2657 . mcp(*),temp(*)
2658C-----------------------------------------------
2659C L o c a l V a r i a b l e s
2660C-----------------------------------------------
2661C
2662C--------------------------------------
2663C READING INTEGERS
2664C--------------------------------------
2665
2666C--------------------------------------
2667C READING REALS
2668C--------------------------------------
2669 IF (irform/5<=1) THEN
2670 CALL redsqr (mcp,numnod,irform)
2671 CALL redsqr (temp,numnod,irform)
2672 ELSE
2673 CALL read_db(mcp,numnod)
2674 CALL read_db(temp,numnod)
2675 ENDIF
2676C
2677 RETURN
2678 END
2679C
2680!||====================================================================
2681!|| nitscherrest ../engine/source/output/restart/rdresb.F
2682!||--- called by ------------------------------------------------------
2683!|| rdresb ../engine/source/output/restart/rdresb.F
2684!||--- calls -----------------------------------------------------
2685!|| read_db ../common_source/tools/input_output/read_db.F
2686!|| redsqr ../engine/source/output/tools/redsqr.F
2687!||====================================================================
2688 SUBROUTINE nitscherrest(FORNEQS)
2689C-----------------------------------------------
2690C I m p l i c i t T y p e s
2691C-----------------------------------------------
2692#include "implicit_f.inc"
2693C-----------------------------------------------
2694C C o m m o n B l o c k s
2695C-----------------------------------------------
2696#include "com04_c.inc"
2697#include "scr05_c.inc"
2698C-----------------------------------------------
2699C D u m m y A r g u m e n t s
2700C-----------------------------------------------
2701 my_real
2702 . forneqs(*)
2703C-----------------------------------------------
2704C L o c a l V a r i a b l e s
2705C-----------------------------------------------
2706C
2707C--------------------------------------
2708C READING INTEGERS
2709C--------------------------------------
2710
2711C--------------------------------------
2712C READING REALS
2713C--------------------------------------
2714 IF (irform/5<=1) THEN
2715 CALL redsqr (forneqs,3*numnod,irform)
2716 ELSE
2717 CALL read_db(forneqs,3*numnod)
2718 ENDIF
2719C
2720 RETURN
2721 END
2722C
2723C convection flow
2724!||====================================================================
2725!|| convrrest ../engine/source/output/restart/rdresb.F
2726!||--- called by ------------------------------------------------------
2727!|| rdresb ../engine/source/output/restart/rdresb.F
2728!||--- calls -----------------------------------------------------
2729!|| read_db ../common_source/tools/input_output/read_db.F
2730!|| read_i_c ../common_source/tools/input_output/write_routines.c
2731!|| redsqi ../engine/source/output/tools/redsqi.F
2732!|| redsqr ../engine/source/output/tools/redsqr.F
2733!||--- uses -----------------------------------------------------
2734!|| glob_therm_mod ../common_source/modules/mat_elem/glob_therm_mod.F90
2735!||====================================================================
2736 SUBROUTINE convrrest(IBCV, FCONV,GLOB_THERM)
2737
2738 use glob_therm_mod
2739C-----------------------------------------------
2740C I m p l i c i t T y p e s
2741C-----------------------------------------------
2742#include "implicit_f.inc"
2743C-----------------------------------------------
2744C C o m m o n B l o c k s
2745C-----------------------------------------------
2746#include "com04_c.inc"
2747#include "scr05_c.inc"
2748#include "param_c.inc"
2749C-----------------------------------------------
2750C D u m m y A r g u m e n t s
2751C-----------------------------------------------
2752 INTEGER IBCV(*)
2753 my_real :: FCONV(*)
2754 type (glob_therm_) :: GLOB_THERM
2755C-----------------------------------------------
2756C L o c a l V a r i a b l e s
2757C-----------------------------------------------
2758C
2759C--------------------------------------
2760C READING INTEGERS
2761C--------------------------------------
2762
2763C--------------------------------------
2764C READING REALS
2765C--------------------------------------
2766 IF (irform/5<=1) THEN
2767 CALL redsqr (fconv,glob_therm%LFACTHER*glob_therm%NUMCONV,irform)
2768 CALL redsqi (ibcv,glob_therm%NICONV*glob_therm%NUMCONV,irform)
2769 ELSE
2770 CALL read_db(fconv,glob_therm%LFACTHER*glob_therm%NUMCONV)
2771 CALL read_i_c(ibcv,glob_therm%NICONV*glob_therm%NUMCONV)
2772 ENDIF
2773C
2774 RETURN
2775 END
2776C radiative flux
2777!||====================================================================
2778!|| radiarrest ../engine/source/output/restart/rdresb.F
2779!||--- called by ------------------------------------------------------
2780!|| rdresb ../engine/source/output/restart/rdresb.F
2781!||--- calls -----------------------------------------------------
2782!|| read_db ../common_source/tools/input_output/read_db.F
2783!|| read_i_c ../common_source/tools/input_output/write_routines.c
2784!|| redsqi ../engine/source/output/tools/redsqi.f
2785!|| redsqr ../engine/source/output/tools/redsqr.F
2786!||--- uses -----------------------------------------------------
2787!|| glob_therm_mod ../common_source/modules/mat_elem/glob_therm_mod.F90
2788!||====================================================================
2789 SUBROUTINE radiarrest(IBCR, FRADIA, GLOB_THERM)
2790
2791 use glob_therm_mod
2792C-----------------------------------------------
2793C I m p l i c i t T y p e s
2794C-----------------------------------------------
2795#include "implicit_f.inc"
2796C-----------------------------------------------
2797C C o m m o n B l o c k s
2798C-----------------------------------------------
2799#include "com04_c.inc"
2800#include "scr05_c.inc"
2801#include "param_c.inc"
2802C-----------------------------------------------
2803C D u m m y A r g u m e n t s
2804C-----------------------------------------------
2805 INTEGER IBCR(*)
2806 my_real :: FRADIA(*)
2807 type (glob_therm_) ,intent(inout) :: GLOB_THERM
2808C-----------------------------------------------
2809C L o c a l V a r i a b l e s
2810C-----------------------------------------------
2811C
2812C--------------------------------------
2813C READING INTEGERS
2814C--------------------------------------
2815
2816C--------------------------------------
2817C READING REALS
2818C--------------------------------------
2819 IF (irform/5<=1) THEN
2820 CALL redsqr (fradia,glob_therm%LFACTHER*glob_therm%NUMRADIA,irform)
2821 CALL redsqi (ibcr,glob_therm%NIRADIA*glob_therm%NUMRADIA,irform)
2822 ELSE
2823 CALL read_db(fradia,glob_therm%LFACTHER*glob_therm%NUMRADIA)
2824 CALL read_i_c(ibcr,glob_therm%NIRADIA*glob_therm%NUMRADIA)
2825 ENDIF
2826C
2827 RETURN
2828 END
2829C imposed temperature
2830!||====================================================================
2831!|| fxtemprrest ../engine/source/output/restart/rdresb.F
2832!||--- called by ------------------------------------------------------
2833!|| rdresb ../engine/source/output/restart/rdresb.F
2834!||--- calls -----------------------------------------------------
2835!|| read_db ../common_source/tools/input_output/read_db.F
2836!|| read_i_c ../common_source/tools/input_output/write_routines.c
2837!|| redsqi ../engine/source/output/tools/redsqi.F
2838!|| redsqr ../engine/source/output/tools/redsqr.F
2839!||--- uses -----------------------------------------------------
2840!|| glob_therm_mod ../common_source/modules/mat_elem/glob_therm_mod.F90
2841!||====================================================================
2842 SUBROUTINE fxtemprrest(IBFTEMP, FBFTEMP,GLOB_THERM)
2843
2844 use glob_therm_mod
2845C-----------------------------------------------
2846C I m p l i c i t T y p e s
2847C-----------------------------------------------
2848#include "implicit_f.inc"
2849C-----------------------------------------------
2850C C o m m o n B l o c k s
2851C-----------------------------------------------
2852#include "com04_c.inc"
2853#include "scr05_c.inc"
2854#include "param_c.inc"
2855C-----------------------------------------------
2856C D u m m y A r g u m e n t s
2857C-----------------------------------------------
2858 INTEGER IBFTEMP(*)
2859 my_real :: FBFTEMP(*)
2860 type (glob_therm_) ,intent(inout) :: GLOB_THERM
2861C-----------------------------------------------
2862C L o c a l V a r i a b l e s
2863C-----------------------------------------------
2864C
2865C--------------------------------------
2866C READING INTEGERS
2867C--------------------------------------
2868
2869C--------------------------------------
2870C READING REALS
2871C--------------------------------------
2872 IF (irform/5<=1) THEN
2873 CALL redsqr (fbftemp,glob_therm%LFACTHER*glob_therm%NFXTEMP,irform)
2874 CALL redsqi (ibftemp,glob_therm%NIFT*glob_therm%NFXTEMP,irform)
2875 ELSE
2876 CALL read_db(fbftemp,glob_therm%LFACTHER*glob_therm%NFXTEMP)
2877 CALL read_i_c(ibftemp,glob_therm%NIFT*glob_therm%NFXTEMP)
2878 ENDIF
2879C
2880 RETURN
2881 END
2882C rigid material
2883!||====================================================================
2884!|| rmatrrest ../engine/source/output/restart/rdresb.F
2885!||--- called by ------------------------------------------------------
2886!|| rdresb ../engine/source/output/restart/rdresb.F
2887!||--- calls -----------------------------------------------------
2888!|| read_db ../common_source/tools/input_output/read_db.F
2889!|| read_i_c ../common_source/tools/input_output/write_routines.c
2890!|| redsqi ../engine/source/output/tools/redsqi.F
2891!|| redsqr ../engine/source/output/tools/redsqr.F
2892!||====================================================================
2893 SUBROUTINE rmatrrest(RBYM,IRBYM, LNRBYM,WEIGHT)
2894C-----------------------------------------------
2895C I m p l i c i t T y p e s
2896C-----------------------------------------------
2897#include "implicit_f.inc"
2898C-----------------------------------------------
2899C C o m m o n B l o c k s
2900C-----------------------------------------------
2901#include "com04_c.inc"
2902#include "scr05_c.inc"
2903C-----------------------------------------------
2904C D u m m y A r g u m e n t s
2905C-----------------------------------------------
2906 INTEGER IRBYM(*), LNRBYM(*),WEIGHT(*)
2907 my_real
2908 . RBYM(*)
2909C-----------------------------------------------
2910C L o c a l V a r i a b l e s
2911C-----------------------------------------------
2912C
2913C--------------------------------------
2914C READING INTEGERS
2915C--------------------------------------
2916
2917C--------------------------------------
2918C READING REALS
2919C--------------------------------------
2920 IF (irform/5<=1) THEN
2921 CALL redsqr (rbym,nrbym*nfrbym,irform)
2922 CALL redsqi(irbym,nrbym*nirbym ,irform)
2923 CALL redsqi(lnrbym,ngslnrbym,irform)
2924 CALL redsqi(weight,nrbym,irform)
2925 ELSE
2926 CALL read_db(rbym,nrbym*nfrbym)
2927 CALL read_i_c(irbym,nrbym*nirbym )
2928 CALL read_i_c(lnrbym,ngslnrbym )
2929 CALL read_i_c(weight,nrbym )
2930 ENDIF
2931C
2932 RETURN
2933 END
2934C shell xfem
2935!||====================================================================
2936!|| plyxfem_rrest ../engine/source/output/restart/rdresb.F
2937!||--- called by ------------------------------------------------------
2938!|| rdresb ../engine/source/output/restart/rdresb.F
2939!||--- calls -----------------------------------------------------
2940!|| read_db ../common_source/tools/input_output/read_db.F
2941!|| read_i_c ../common_source/tools/input_output/write_routines.c
2942!||====================================================================
2943 SUBROUTINE plyxfem_rrest(MS_LAYER,ZI_LAYER,INOD,IEL,ICODE,ISKEW,
2944 . MSZ2)
2945C-----------------------------------------------
2946C I m p l i c i t T y p e s
2947C-----------------------------------------------
2948#include "implicit_f.inc"
2949C-----------------------------------------------
2950C C o m m o n B l o c k s
2951C-----------------------------------------------
2952#include "com01_c.inc"
2953#include "com04_c.inc"
2954#include "param_c.inc"
2955C-----------------------------------------------
2956C D u m m y A r g u m e n t s
2957C-----------------------------------------------
2958 my_real
2959 . ms_layer(*),zi_layer(*),msz2(*)
2960 integer
2961 . icode(*),iskew(*),inod(*),iel(*)
2962C-----------------------------------------------
2963C L o c a l V a r i a b l e s
2964C-----------------------------------------------
2965C
2966C--------------------------------------
2967C READING INTEGERS
2968C--------------------------------------
2969
2970C--------------------------------------
2971C READING REALS
2972C--------------------------------------
2973 CALL read_db(ms_layer,nplyxfe*nplymax)
2974 CALL read_db(zi_layer,nplyxfe*nplymax)
2975 CALL read_db(msz2,nplyxfe)
2976C
2977 CALL read_i_c(inod,numnod)
2978 CALL read_i_c(iel,numelc)
2979C BCS ply
2980 IF(iplybcs > 0) THEN
2981 CALL read_i_c(icode,numnod)
2982 CALL read_i_c(iskew,numnod)
2983 ENDIF
2984 RETURN
2985 END
2986C shell xfem
2987!||====================================================================
2988!|| plyxfem_rrestanim ../engine/source/output/restart/rdresb.F
2989!||--- called by ------------------------------------------------------
2990!|| rdresb ../engine/source/output/restart/rdresb.F
2991!||--- calls -----------------------------------------------------
2992!|| read_i_c ../common_source/tools/input_output/write_routines.c
2993!||--- uses -----------------------------------------------------
2994!|| plyxfem_mod ../engine/share/modules/plyxfem_mod.F
2995!||====================================================================
2996 SUBROUTINE plyxfem_rrestanim()
2997C-----------------------------------------------
2998 USE plyxfem_mod
2999C-----------------------------------------------
3000C I m p l i c i t T y p e s
3001C-----------------------------------------------
3002#include "implicit_f.inc"
3003C-----------------------------------------------
3004C C o m m o n B l o c k s
3005C-----------------------------------------------
3006#include "com01_c.inc"
3007#include "param_c.inc"
3008C-----------------------------------------------
3009C L o c a l V a r i a b l e s
3010C-----------------------------------------------
3011 INTEGER K,ELSZ,NDSZ
3012C-----------------------------------------------
3013 IF(IPLYXFEM == 0) then
3014! Those array be used as dummy argument
3015! they MUST be allocated in all cases
3016 ALLOCATE(indx_ply(0))
3017 ALLOCATE(plyshell(0))
3018 ALLOCATE(plynod(0))
3019 ELSE
3020 CALL read_i_c(nplypart,1)
3021 ALLOCATE(indx_ply(nplypart))
3022 CALL read_i_c(indx_ply,nplypart)
3023
3024 ALLOCATE(plyshell(nplymax))
3025 ALLOCATE(plynod(nplymax))
3026
3027 DO k=1,nplymax
3028 CALL read_i_c(elsz,1)
3029 plyshell(k)%PLYNUMSHELL=elsz
3030 ALLOCATE(plyshell(k)%SHID(elsz))
3031 ALLOCATE(plyshell(k)%SHELLIPT(elsz))
3032 ALLOCATE(plyshell(k)%SHELLID(elsz))
3033
3034 CALL read_i_c(plyshell(k)%SHID, elsz)
3035 CALL read_i_c(plyshell(k)%SHELLIPT, elsz)
3036 CALL read_i_c(plyshell(k)%SHELLID, elsz)
3037
3038
3039 CALL read_i_c(ndsz,1)
3040 plynod(k)%PLYNUMNODS=ndsz
3041 ALLOCATE(plynod(k)%NODES(ndsz))
3042 ALLOCATE(plynod(k)%PLYNODID(ndsz))
3043
3044 CALL read_i_c(plynod(k)%NODES, ndsz)
3045 CALL read_i_c(plynod(k)%PLYNODID, ndsz)
3046 ENDDO
3047
3048 ALLOCATE(plysizg(nplymax))
3049 CALL read_i_c(plysizg,nplymax)
3050C
3051 ALLOCATE (idpid_ply(nplypart))
3052 CALL read_i_c(idpid_ply,nplypart)
3053 ENDIF
3054 RETURN
3055 END
3056C centrifugal force
3057!||====================================================================
3058!|| cfieldrest ../engine/source/output/restart/rdresb.F
3059!||--- called by ------------------------------------------------------
3060!|| rdresb ../engine/source/output/restart/rdresb.F
3061!||--- calls -----------------------------------------------------
3062!|| read_db ../common_source/tools/input_output/read_db.F
3063!||--- uses -----------------------------------------------------
3064!|| message_mod ../engine/share/message_module/message_mod.F
3065!||====================================================================
3066 SUBROUTINE cfieldrest(CFIELD)
3067C-----------------------------------------------
3068C M o d u l e s
3069C-----------------------------------------------
3070 USE message_mod
3071C-----------------------------------------------
3072C I m p l i c i t T y p e s
3073C-----------------------------------------------
3074#include "implicit_f.inc"
3075C-----------------------------------------------
3076C C o m m o n B l o c k s
3077C-----------------------------------------------
3078#include "com04_c.inc"
3079#include "param_c.inc"
3080C-----------------------------------------------
3081C D u m m y A r g u m e n t s
3082C-----------------------------------------------
3083 my_real
3084 . cfield(*)
3085C-----------------------------------------------
3086C L o c a l V a r i a b l e s
3087C-----------------------------------------------
3088C
3089C--------------------------------------
3090C READING INTEGERS
3091C--------------------------------------
3092
3093C--------------------------------------
3094C READING REALS
3095C--------------------------------------
3096 CALL read_db(cfield,lfacload*nloadc)
3097C
3098 RETURN
3099 END
3100!||====================================================================
3101!|| loadprest ../engine/source/output/restart/rdresb.F
3102!||--- called by ------------------------------------------------------
3103!|| rdresb ../engine/source/output/restart/rdresb.F
3104!||--- calls -----------------------------------------------------
3105!|| read_db ../common_source/tools/input_output/read_db.F
3106!||--- uses -----------------------------------------------------
3107!|| message_mod ../engine/share/message_module/message_mod.F
3108!||====================================================================
3109 SUBROUTINE loadprest(LOADP)
3110C-----------------------------------------------
3111C M o d u l e s
3112C-----------------------------------------------
3113 USE message_mod
3114C-----------------------------------------------
3115C I m p l i c i t T y p e s
3116C-----------------------------------------------
3117#include "implicit_f.inc"
3118C-----------------------------------------------
3119C C o m m o n B l o c k s
3120C-----------------------------------------------
3121#include "com04_c.inc"
3122#include "param_c.inc"
3123C-----------------------------------------------
3124C D u m m y A r g u m e n t s
3125C-----------------------------------------------
3126 my_real loadp(*)
3127C-----------------------------------------------
3128C L o c a l V a r i a b l e s
3129C-----------------------------------------------
3130C
3131C--------------------------------------
3132C READING INTEGERS
3133C--------------------------------------
3134
3135C--------------------------------------
3136C READING REALS
3137C--------------------------------------
3138 CALL read_db(loadp,lfacload*nloadp)
3139C
3140 RETURN
3141 END
3142
3143!||====================================================================
3144!|| crkxfem_rrest ../engine/source/output/restart/rdresb.F
3145!||--- called by ------------------------------------------------------
3146!|| rdresb ../engine/source/output/restart/rdresb.F
3147!||--- calls -----------------------------------------------------
3148!|| read_i_c ../common_source/tools/input_output/write_routines.c
3149!||====================================================================
3150 SUBROUTINE crkxfem_rrest(INOD_CRKXFEM,IEL_CRKXFEM,NODLEVXF)
3151C-----------------------------------------------
3152C I m p l i c i t T y p e s
3153C-----------------------------------------------
3154#include "implicit_f.inc"
3155C-----------------------------------------------
3156C C o m m o n B l o c k s
3157C-----------------------------------------------
3158#include "com04_c.inc"
3159#include "com_xfem1.inc"
3160C-----------------------------------------------
3161C D u m m y A r g u m e n t s
3162C-----------------------------------------------
3163 INTEGER INOD_CRKXFEM(*),IEL_CRKXFEM(*),NODLEVXF(*)
3164C-----------------------------------------------
3165C L o c a l V a r i a b l e s
3166C-----------------------------------------------
3167C=======================================================================
3168 CALL READ_I_C(INOD_CRKXFEM,NUMNOD)
3169 CALL READ_I_C(IEL_CRKXFEM,NUMELC+NUMELTG)
3170 CALL READ_I_C(NODLEVXF,NCRKXFE)
3171C---
3172 RETURN
3173 END
3174!||====================================================================
3175!|| crkxfem_rrestanim ../engine/source/output/restart/rdresb.F
3176!||--- called by ------------------------------------------------------
3177!|| rdresb ../engine/source/output/restart/rdresb.F
3178!||--- calls -----------------------------------------------------
3179!|| read_db ../common_source/tools/input_output/read_db.F
3180!|| read_i_c ../common_source/tools/input_output/write_routines.c
3181!||--- uses -----------------------------------------------------
3182!|| restmod ../engine/share/modules/restart_mod.F
3183!|| xfem2vars_mod ../engine/share/modules/xfem2vars_mod.F
3184!||====================================================================
3185 SUBROUTINE crkxfem_rrestanim()
3186C-----------------------------------------------
3187 USE xfem2vars_mod
3188 USE restmod
3189C-----------------------------------------------
3190C I m p l i c i t T y p e s
3191C-----------------------------------------------
3192#include "implicit_f.inc"
3193C-----------------------------------------------
3194C C o m m o n B l o c k s
3195C-----------------------------------------------
3196#include "com01_c.inc"
3197#include "com04_c.inc"
3198#include "com_xfem1.inc"
3199C-----------------------------------------------
3200C L o c a l V a r i a b l e s
3201C-----------------------------------------------
3202 INTEGER K,LEN,ELSZ,NENR,ELSZC,ELSZTG,NCRKXFE_G,
3203 . CRKNUMNODS
3204C=======================================================================
3205 IF(ICRACK3D == 0) then
3206 ALLOCATE(indx_crk(0))
3207 ALLOCATE(crkshell(0))
3208 ALLOCATE(crknod(0))
3209 ALLOCATE(crklvset(0))
3210 ALLOCATE(crkavx(0))
3211 ALLOCATE(crksky(0))
3212 ALLOCATE(xfem_phantom(0))
3213 ALLOCATE(crkedge(0))
3214 ALLOCATE(xedge4n(4,0))
3215 ALLOCATE(xedge3n(3,0))
3216 ELSE
3217 CALL read_i_c(ncrkpart,1)
3218 CALL read_i_c(ncrkxfe_g,1)
3219 ALLOCATE(indx_crk(ncrkpart))
3220 CALL read_i_c(indx_crk,ncrkpart)
3221C
3222 ALLOCATE(crkshell(nlevmax))
3223 ALLOCATE(crknod(nlevmax))
3224 ALLOCATE(crklvset(nlevmax))
3225 ALLOCATE(crkavx(nlevmax))
3226 ALLOCATE(crksky(nlevmax))
3227C
3228 DO k=1,nlevmax
3229 CALL read_i_c(elszc ,1)
3230 CALL read_i_c(elsztg ,1)
3231 CALL read_i_c(elsz ,1)
3232c
3233c-------- CRKSHELL
3234 crkshell(k)%CRKNUMSHELL = elsz
3235 crkshell(k)%CRKNUMSH4 = elszc
3236 crkshell(k)%CRKNUMSH3 = elsztg
3237C
3238 ALLOCATE(crkshell(k)%CRKSHELLID(elsz))
3239 ALLOCATE(crkshell(k)%ELTYPE(elsz))
3240 CALL read_i_c(crkshell(k)%CRKSHELLID ,elsz)
3241 CALL read_i_c(crkshell(k)%ELTYPE ,elsz)
3242c-------- CRKNOD
3243 len = 4*elsz
3244 ALLOCATE(crknod(k)%XFECRKNODID(len))
3245 ALLOCATE(crknod(k)%XFENODES(len))
3246 CALL read_i_c(crknod(k)%XFECRKNODID ,len)
3247 CALL read_i_c(crknod(k)%XFENODES ,len)
3248c-------- CRKSHELL
3249 len = 4*elsz
3250 ALLOCATE(crkshell(k)%XNODEL(4 ,elsz))
3251 CALL read_i_c(crkshell(k)%XNODEL, len)
3252c-------- CRKLVSET
3253 len = lcnecrkxfem
3254 ALLOCATE(crklvset(k)%ENR0(2,len))
3255 ALLOCATE(crklvset(k)%AREA(elsz))
3256 CALL read_i_c(crklvset(k)%ENR0(1,1:len) ,len)
3257 CALL read_i_c(crklvset(k)%ENR0(2,1:len) ,len)
3258 CALL read_db (crklvset(k)%AREA ,elsz)
3259c-------- CRKAVX
3260 len = lcnecrkxfem
3261 ALLOCATE(crkavx(k)%A(3,len))
3262 ALLOCATE(crkavx(k)%AR(3,len))
3263 ALLOCATE(crkavx(k)%V(3,len))
3264 ALLOCATE(crkavx(k)%VR(3,len))
3265 ALLOCATE(crkavx(k)%X(3,len))
3266 ALLOCATE(crkavx(k)%U(3,len))
3267 len = lcnecrkxfem*3
3268 CALL read_db(crkavx(k)%A ,len)
3269 CALL read_db(crkavx(k)%AR ,len)
3270 CALL read_db(crkavx(k)%V ,len)
3271 CALL read_db(crkavx(k)%VR ,len)
3272 CALL read_db(crkavx(k)%X ,len)
3273 CALL read_db(crkavx(k)%U ,len)
3274 len = 4*elsz
3275 ALLOCATE(crkavx(k)%XX(3,len))
3276 CALL read_db(crkavx(k)%XX, 3*len)
3277c-------- CRKSKY
3278 ALLOCATE(crksky(k)%FSKY(8,lcnecrkxfem))
3279 CALL read_db(crksky(k)%FSKY, 8*lcnecrkxfem)
3280c-------- CRKNOD
3281 len = 4*elsz
3282 ALLOCATE(crknod(k)%NOD2IAD(len))
3283 CALL read_i_c(crknod(k)%NOD2IAD, len)
3284 CALL read_i_c(crknumnods, 1)
3285 crknod(k)%CRKNUMNODS = crknumnods
3286c
3287 ENDDO ! K=1,NLEVMAX
3288c
3289c--- read XFEM_PHANTOM ----------------------------
3290c
3291
3292 nenr = int(ienrnod/nlevmax)
3293 ALLOCATE(xfem_phantom(nxlaymax))
3294 DO k=1,nxlaymax
3295 ALLOCATE(xfem_phantom(k)%ELCUT(elsz))
3296 ALLOCATE(xfem_phantom(k)%IFI(lcnecrkxfem))
3297 ALLOCATE(xfem_phantom(k)%TAGXP(5,ncrkxfe,ienrnod)) ! IENRNOD -> NENR
3298 ALLOCATE(xfem_phantom(k)%ITRI(2,elsz))
3299c
3300 CALL read_i_c(xfem_phantom(k)%ELCUT ,elsz)
3301 CALL read_i_c(xfem_phantom(k)%IFI ,lcnecrkxfem)
3302 CALL read_i_c(xfem_phantom(k)%TAGXP ,ncrkxfe*ienrnod*5)
3303 CALL read_i_c(xfem_phantom(k)%ITRI ,elsz*2)
3304 END DO
3305c
3306c--- read CRKEDGE ----------------------------
3307c
3308 len = numedges
3309 ALLOCATE(crkedge(nxlaymax))
3310 DO k=1,nxlaymax
3311 ALLOCATE(crkedge(k)%LAYCUT(elsz))
3312 ALLOCATE(crkedge(k)%IEDGEC(4,elszc))
3313 ALLOCATE(crkedge(k)%IEDGETG(3,elsztg))
3314 ALLOCATE(crkedge(k)%EDGEICRK(numedges))
3315 ALLOCATE(crkedge(k)%EDGEIFI(2,numedges))
3316 ALLOCATE(crkedge(k)%EDGEENR(2,numedges))
3317 ALLOCATE(crkedge(k)%EDGETIP(2,numedges))
3318 ALLOCATE(crkedge(k)%IBORDEDGE(numedges))
3319 ALLOCATE(crkedge(k)%ICUTEDGE(numedges))
3320 ALLOCATE(crkedge(k)%RATIO(numedges))
3321c
3322 CALL read_i_c(crkedge(k)%LAYCUT ,elsz)
3323 CALL read_i_c(crkedge(k)%IEDGEC ,elszc*4)
3324 CALL read_i_c(crkedge(k)%IEDGETG ,elsztg*3)
3325 CALL read_i_c(crkedge(k)%EDGEICRK ,numedges)
3326 CALL read_i_c(crkedge(k)%EDGEIFI ,numedges*2)
3327 CALL read_i_c(crkedge(k)%EDGEENR ,numedges*2)
3328 CALL read_i_c(crkedge(k)%EDGETIP ,numedges*2)
3329 CALL read_i_c(crkedge(k)%IBORDEDGE ,numedges)
3330 CALL read_i_c(crkedge(k)%ICUTEDGE ,numedges)
3331 CALL read_db (crkedge(k)%RATIO ,numedges)
3332 END DO
3333c-------------------------------
3334 ALLOCATE(xedge4n(4,elszc))
3335 ALLOCATE(xedge3n(3,elsztg))
3336 CALL read_i_c(xedge4n, 4*elszc)
3337 CALL read_i_c(xedge3n, 3*elsztg)
3338c-------------------------------
3339 len = 2*(numelc+numeltg)
3340 CALL read_i_c(elcutc,len)
3341 len = ncrkxfe
3342 CALL read_i_c(nodenr,len)
3343 len = ncrkxfe
3344 CALL read_i_c(kxfenod2elc,len)
3345 len = numnod*ienrnod
3346 CALL read_i_c(enrtag,len)
3347C---
3348 ENDIF
3349 RETURN
3350 END
3351
3352C shell xfem
3353!||====================================================================
3354!|| plyxfem_ravuply ../engine/source/output/restart/rdresb.F
3355!||--- called by ------------------------------------------------------
3356!|| rdresb ../engine/source/output/restart/rdresb.F
3357!||--- calls -----------------------------------------------------
3358!|| read_db ../common_source/tools/input_output/read_db.F
3359!||--- uses -----------------------------------------------------
3360!|| plyxfem_mod ../engine/share/modules/plyxfem_mod.F
3361!||====================================================================
3362 SUBROUTINE plyxfem_ravuply()
3363C-----------------------------------------------
3364 USE plyxfem_mod
3365C-----------------------------------------------
3366C I m p l i c i t T y p e s
3367C-----------------------------------------------
3368#include "implicit_f.inc"
3369C-----------------------------------------------
3370C C o m m o n B l o c k s
3371C-----------------------------------------------
3372#include "param_c.inc"
3373C-----------------------------------------------
3374C L o c a l V a r i a b l e s
3375C-----------------------------------------------
3376 INTEGER K,NOSZ
3377C-----------------------------------------------
3378 ALLOCATE(PLY(NPLYMAX))
3379 DO K=1,nplymax
3380 NULLIFY(ply(k)%V,ply(k)%U )
3381 ALLOCATE(ply(k)%V(3,nplyxfe), ply(k)%U(3,nplyxfe))
3382 nosz = 3*nplyxfe
3383 CALL read_db(ply(k)%U ,nosz)
3384 CALL read_db(ply(k)%V ,nosz)
3385 ENDDO
3386C
3387 RETURN
3388 END
3389C air + faom
3390!||====================================================================
3391!|| alelag_rrest ../engine/source/output/restart/rdresb.F
3392!||--- called by ------------------------------------------------------
3393!|| rdresb ../engine/source/output/restart/rdresb.F
3394!||--- calls -----------------------------------------------------
3395!|| read_db ../common_source/tools/input_output/read_db.F
3396!||--- uses -----------------------------------------------------
3397!|| aleflow_mod ../engine/share/modules/aleflow_mod.F
3398!||====================================================================
3399 SUBROUTINE alelag_rrest()
3400C-----------------------------------------------
3401 USE aleflow_mod
3402C-----------------------------------------------
3403C I m p l i c i t T y p e s
3404C-----------------------------------------------
3405#include "implicit_f.inc"
3406C-----------------------------------------------
3407C C o m m o n B l o c k s
3408C-----------------------------------------------
3409#include "com04_c.inc"
3410C-----------------------------------------------
3411C L o c a l V a r i a b l e s
3412C-----------------------------------------------
3413 INTEGER NOSZ
3414C-----------------------------------------------
3415 ALLOCATE (VFLOW(3*NUMNOD), DFLOW(3*NUMNOD),WFLOW(3*NUMNOD))
3416 NOSZ = 3*numnod
3417 CALL read_db(dflow, nosz)
3418 CALL read_db(vflow ,nosz)
3419 CALL read_db(wflow ,nosz)
3420C
3421 RETURN
3422 END
3423!||====================================================================
3424!|| rthbufrest ../engine/source/output/restart/rdresb.F
3425!||--- called by ------------------------------------------------------
3426!|| rdresb ../engine/source/output/restart/rdresb.F
3427!||--- calls -----------------------------------------------------
3428!|| read_db ../common_source/tools/input_output/read_db.F
3429!||--- uses -----------------------------------------------------
3430!|| message_mod ../engine/share/message_module/message_mod.F
3431!||====================================================================
3432 SUBROUTINE rthbufrest(RTHBUF,SRTHBUF)
3433C-----------------------------------------------
3434C M o d u l e s
3435C-----------------------------------------------
3436 USE message_mod
3437C-----------------------------------------------
3438C I m p l i c i t T y p e s
3439C-----------------------------------------------
3440#include "implicit_f.inc"
3441C-----------------------------------------------
3442C D u m m y A r g u m e n t s
3443C-----------------------------------------------
3444 my_real
3445 . rthbuf(*)
3446 integer
3447 . srthbuf
3448C-----------------------------------------------
3449C L o c a l V a r i a b l e s
3450C-----------------------------------------------
3451C
3452C--------------------------------------
3453C READING INTEGERS
3454C--------------------------------------
3455
3456C--------------------------------------
3457C READING REALS
3458C--------------------------------------
3459 CALL read_db(rthbuf,srthbuf)
3460C
3461 RETURN
3462 END
3463
3464!||====================================================================
3465!|| knotrest ../engine/source/output/restart/rdresb.F
3466!||--- called by ------------------------------------------------------
3467!|| rdresb ../engine/source/output/restart/rdresb.F
3468!||--- calls -----------------------------------------------------
3469!|| read_db ../common_source/tools/input_output/read_db.F
3470!||--- uses -----------------------------------------------------
3471!|| message_mod ../engine/share/message_module/message_mod.F
3472!||====================================================================
3473 SUBROUTINE knotrest(KNOT)
3474C-----------------------------------------------
3475C M o d u l e s
3476C-----------------------------------------------
3477 USE message_mod
3478C-----------------------------------------------
3479C I m p l i c i t T y p e s
3480C-----------------------------------------------
3481#include "implicit_f.inc"
3482C-----------------------------------------------
3483C C o m m o n B l o c k s
3484C-----------------------------------------------
3485#include "tabsiz_c.inc"
3486C-----------------------------------------------
3487C D u m m y A r g u m e n t s
3488C-----------------------------------------------
3489 my_real
3490 . knot(*)
3491C-----------------------------------------------
3492C L o c a l V a r i a b l e s
3493C-----------------------------------------------
3494C
3495C--------------------------------------
3496C READING INTEGERS
3497C--------------------------------------
3498C--------------------------------------
3499C READING REALS
3500C--------------------------------------
3501 CALL read_db(knot,sknot)
3502C
3503 RETURN
3504 END
3505!||====================================================================
3506!|| knotlocpcrest ../engine/source/output/restart/rdresb.F
3507!||--- called by ------------------------------------------------------
3508!|| rdresb ../engine/source/output/restart/rdresb.F
3509!||--- calls -----------------------------------------------------
3510!|| read_db ../common_source/tools/input_output/read_db.F
3511!||--- uses -----------------------------------------------------
3512!|| message_mod ../engine/share/message_module/message_mod.F
3513!||====================================================================
3514 SUBROUTINE knotlocpcrest(KNOTLOCPC)
3515C-----------------------------------------------
3516C M o d u l e s
3517C-----------------------------------------------
3518 USE message_mod
3519C-----------------------------------------------
3520C I m p l i c i t T y p e s
3521C-----------------------------------------------
3522#include "implicit_f.inc"
3523C-----------------------------------------------
3524C C o m m o n B l o c k s
3525C-----------------------------------------------
3526#include "tabsiz_c.inc"
3527C-----------------------------------------------
3528C D u m m y A r g u m e n t s
3529C-----------------------------------------------
3530 my_real
3531 . knotlocpc(*)
3532C-----------------------------------------------
3533C L o c a l V a r i a b l e s
3534C-----------------------------------------------
3535C
3536C--------------------------------------
3537C READING INTEGERS
3538C--------------------------------------
3539C--------------------------------------
3540C READING REALS
3541C--------------------------------------
3542 CALL read_db(knotlocpc,sknotlocpc)
3543C
3544 RETURN
3545 END
3546!||====================================================================
3547!|| knotlocelrest ../engine/source/output/restart/rdresb.F
3548!||--- called by ------------------------------------------------------
3549!|| rdresb ../engine/source/output/restart/rdresb.F
3550!||--- calls -----------------------------------------------------
3551!|| read_db ../common_source/tools/input_output/read_db.F
3552!||--- uses -----------------------------------------------------
3553!|| message_mod ../engine/share/message_module/message_mod.F
3554!||====================================================================
3555 SUBROUTINE knotlocelrest(KNOTLOCEL)
3556C-----------------------------------------------
3557C M o d u l e s
3558C-----------------------------------------------
3559 USE message_mod
3560C-----------------------------------------------
3561C I m p l i c i t T y p e s
3562C-----------------------------------------------
3563#include "implicit_f.inc"
3564C-----------------------------------------------
3565C C o m m o n B l o c k s
3566C-----------------------------------------------
3567#include "tabsiz_c.inc"
3568C-----------------------------------------------
3569C D u m m y A r g u m e n t s
3570C-----------------------------------------------
3571 my_real knotlocel(*)
3572C-----------------------------------------------
3573C L o c a l V a r i a b l e s
3574C-----------------------------------------------
3575C
3576C--------------------------------------
3577C READING INTEGERS
3578C--------------------------------------
3579C--------------------------------------
3580C READING REALS
3581C--------------------------------------
3582 CALL read_db(knotlocel,sknotlocel)
3583C
3584 RETURN
3585 END
3586!||====================================================================
3587!|| wigerest ../engine/source/output/restart/rdresb.F
3588!||--- called by ------------------------------------------------------
3589!|| rdresb ../engine/source/output/restart/rdresb.f
3590!||--- calls -----------------------------------------------------
3591!|| read_db ../common_source/tools/input_output/read_db.F
3592!||--- uses -----------------------------------------------------
3593!|| message_mod ../engine/share/message_module/message_mod.F
3594!||====================================================================
3595 SUBROUTINE wigerest(WIGE)
3596C-----------------------------------------------
3597C M o d u l e s
3598C-----------------------------------------------
3599 USE message_mod
3600C-----------------------------------------------
3601C I m p l i c i t T y p e s
3602C-----------------------------------------------
3603#include "implicit_f.inc"
3604C-----------------------------------------------
3605C C o m m o n B l o c k s
3606C-----------------------------------------------
3607#include "com04_c.inc"
3608C-----------------------------------------------
3609C D u m m y A r g u m e n t s
3610C-----------------------------------------------
3611 my_real wige(*)
3612C-----------------------------------------------
3613C L o c a l V a r i a b l e s
3614C-----------------------------------------------
3615C--------------------------------------
3616C READING INTEGERS
3617C--------------------------------------
3618C--------------------------------------
3619C READING REALS
3620C--------------------------------------
3621 CALL read_db(wige,numnod)
3622C
3623 RETURN
3624 END
3625C STACk PART
3626!||====================================================================
3627!|| stack_rrest ../engine/source/output/restart/rdresb.F
3628!||--- called by ------------------------------------------------------
3629!|| rdresb ../engine/source/output/restart/rdresb.F
3630!||--- calls -----------------------------------------------------
3631!|| read_db ../common_source/tools/input_output/read_db.F
3632!|| read_i_c ../common_source/tools/input_output/write_routines.c
3633!||====================================================================
3634 SUBROUTINE stack_rrest(IGEO,GEO,PM)
3635C-----------------------------------------------
3636C I m p l i c i t T y p e s
3637C-----------------------------------------------
3638#include "implicit_f.inc"
3639C-----------------------------------------------
3640C C o m m o n B l o c k s
3641C-----------------------------------------------
3642#include "param_c.inc"
3643C-----------------------------------------------
3644C D u m m y A r g u m e n t s
3645C-----------------------------------------------
3646 INTEGER IGEO(*)
3647 my_real GEO(*),PM(*)
3648C-----------------------------------------------
3649C L o c a l V a r i a b l e s
3650C-----------------------------------------------
3651 INTEGER LSIZ
3652
3653C-----------------------------------------------
3654 LSIZ =(4*npt_stack +2)*ns_stack
3655 CALL read_i_c(igeo, lsiz)
3656 lsiz =(6*npt_stack + 1)*ns_stack
3657 CALL read_db(geo, lsiz)
3658 lsiz =20*ns_stack
3659 CALL read_db(pm, lsiz)
3660 RETURN
3661 END
3662!||====================================================================
3663!|| drape_rrest ../engine/source/output/restart/rdresb.f
3664!||--- called by ------------------------------------------------------
3665!|| rdresb ../engine/source/output/restart/rdresb.f
3666!||--- calls -----------------------------------------------------
3667!|| read_db ../common_source/tools/input_output/read_db.F
3668!|| read_i_c ../common_source/tools/input_output/write_routines.c
3669!||--- uses -----------------------------------------------------
3670!|| drape_mod ../engine/share/modules/drape_mod.F
3671!||====================================================================
3672 SUBROUTINE drape_rrest(DRAPE_SH4N,DRAPE_SH3N,DRAPEG)
3673C
3674 USE drape_mod
3675C-----------------------------------------------
3676C I m p l i c i t T y p e s
3677C-----------------------------------------------
3678#include "implicit_f.inc"
3679C-----------------------------------------------
3680C C o m m o n B l o c k s
3681C-----------------------------------------------
3682#include "com04_c.inc"
3683C-----------------------------------------------
3684C D u m m y A r g u m e n t s
3685C-----------------------------------------------
3686 TYPE(drape_), DIMENSION(*) :: DRAPE_SH4N,DRAPE_SH3N
3687 TYPE(DRAPEG_) :: DRAPEG
3688C-----------------------------------------------
3689C L o c a l V a r i a b l e s
3690C-----------------------------------------------
3691 INTEGER LEN,OFF,II,J,IS,NSLICE,NPT,NDIM_INDX_PLY,NDIM,NDIM_SLICE,
3692 . NPT_DRAPE
3693 my_real, DIMENSION(:,:), ALLOCATABLE :: RTMP
3694 my_real, DIMENSION(:) , ALLOCATABLE :: THK
3695 INTEGER, DIMENSION(:,:), ALLOCATABLE :: ITMP
3696 INTEGER, DIMENSION(:) , ALLOCATABLE :: NPTDRP_ELT,ISLICE,INDX_PLY,
3697 . NPT_ELT
3698C--------------------------------------
3699C READING INTEGERS
3700C--------------------------------------
3701C 4 node Shell
3702 scdrape = 0
3703 stdrape = 0
3704 IF(numelc_drape > 0) THEN
3705 CALL read_i_c(ndim ,1)
3706 CALL read_i_c(ndim_slice,1)
3707 CALL read_i_c(ndim_indx_ply,1)
3708 drapeg%NDIM_SH4N = ndim
3709 drapeg%NDIM_SLICE_SH4N = ndim_slice
3710 drapeg%NDIM_INDX_PLY_SH4N = ndim_indx_ply
3711 !
3712 ALLOCATE(npt_elt(numelc_drape),itmp(ndim,2),rtmp(ndim,2),
3713 . nptdrp_elt(numelc_drape),islice(ndim_slice),indx_ply(ndim_indx_ply),
3714 . thk(numelc_drape) )
3715 thk = zero
3716 CALL read_i_c(drapeg%INDX_SH4N,numelc)
3717 CALL read_i_c(nptdrp_elt,numelc_drape)
3718 CALL read_i_c(npt_elt,numelc_drape)
3719 CALL read_i_c(islice,ndim_slice )
3720 CALL read_i_c(indx_ply,ndim_indx_ply)
3721 CALL read_i_c(itmp,2*ndim)
3722 CALL read_db(rtmp,2*ndim)
3723 CALL read_db(thk,numelc_drape)
3724 ndim_indx_ply = 0
3725 ndim_slice = 0
3726 ndim = 0
3727 DO ii= 1,numelc_drape
3728 drape_sh4n(ii)%NPLY_DRAPE = nptdrp_elt(ii)
3729 drape_sh4n(ii)%NPLY = npt_elt(ii)
3730 drape_sh4n(ii)%THICK = thk(ii)
3731 npt = drape_sh4n(ii)%NPLY
3732 ALLOCATE(drape_sh4n(ii)%INDX_PLY(npt))
3733 drape_sh4n(ii)%INDX_PLY = 0
3734 DO j = 1,npt
3735 drape_sh4n(ii)%INDX_PLY(j)=indx_ply(ndim_indx_ply + j)
3736 ENDDO
3737 ndim_indx_ply = ndim_indx_ply + npt
3738 npt_drape = drape_sh4n(ii)%NPLY_DRAPE
3739 ALLOCATE(drape_sh4n(ii)%DRAPE_PLY(npt_drape))
3740 DO j=1 ,npt_drape
3741 nslice = islice(ndim_slice + j)
3742 drape_sh4n(ii)%DRAPE_PLY(j)%NSLICE = nslice
3743 ALLOCATE(drape_sh4n(ii)%DRAPE_PLY(j)%IDRAPE(nslice,2),drape_sh4n(ii)%DRAPE_PLY(j)%RDRAPE(nslice,2))
3744 drape_sh4n(ii)%DRAPE_PLY(j)%IDRAPE = 0
3745 drape_sh4n(ii)%DRAPE_PLY(j)%RDRAPE = zero
3746 DO is= 1,nslice
3747 drape_sh4n(ii)%DRAPE_PLY(j)%IDRAPE(is,1) = itmp(ndim + is,1)
3748 drape_sh4n(ii)%DRAPE_PLY(j)%IDRAPE(is,2) = itmp(ndim + is,2)
3749 drape_sh4n(ii)%DRAPE_PLY(j)%RDRAPE(is,1) = rtmp(ndim + is,1)
3750 drape_sh4n(ii)%DRAPE_PLY(j)%RDRAPE(is,2) = rtmp(ndim + is,2)
3751 ENDDO
3752 ndim = ndim + nslice
3753 ENDDO
3754 ndim_slice = ndim_slice + npt_drape
3755 ENDDO ! NUMELC_DRAPE
3756 DEALLOCATE(npt_elt,itmp,rtmp,nptdrp_elt,islice,indx_ply, thk)
3757 scdrape = numelc
3758 ENDIF
3759C
3760C 3 Node Shell
3761 IF(numeltg_drape > 0) THEN
3762 CALL read_i_c(ndim ,1)
3763 CALL read_i_c(ndim_slice,1)
3764 Call read_i_c(ndim_indx_ply,1)
3765 drapeg%NDIM_SH3N = ndim
3766 drapeg%NDIM_SLICE_SH3N = ndim_slice
3767 drapeg%NDIM_INDX_PLY_SH3N = ndim_indx_ply
3768 ALLOCATE( npt_elt(numeltg_drape),itmp(ndim,2),rtmp(ndim,2),
3769 . nptdrp_elt(numeltg_drape),islice(ndim_slice),indx_ply(ndim_indx_ply),
3770 . thk(numeltg_drape))
3771 thk = zero
3772 CALL read_i_c(drapeg%INDX_SH3N,numeltg)
3773 CALL read_i_c(nptdrp_elt,numeltg_drape)
3774 CALL read_i_c(npt_elt,numeltg_drape)
3775 CALL read_i_c(islice,ndim_slice)
3776 CALL read_i_c(indx_ply,ndim_indx_ply)
3777 CALL read_i_c(itmp,2*ndim)
3778 CALL read_db(rtmp,2*ndim)
3779 CALL read_db(thk,numeltg_drape)
3780 ndim_indx_ply = 0
3781 ndim_slice = 0
3782 ndim = 0
3783 DO ii= 1,numeltg_drape
3784 drape_sh3n(ii)%NPLY_DRAPE = nptdrp_elt(ii)
3785 drape_sh3n(ii)%NPLY = npt_elt(ii)
3786 drape_sh3n(ii)%THICK = thk(ii)
3787 npt = drape_sh3n(ii)%NPLY
3788 ALLOCATE(drape_sh3n(ii)%INDX_PLY(npt))
3789 drape_sh3n(ii)%INDX_PLY = 0
3790 DO j = 1,npt
3791 drape_sh3n(ii)%INDX_PLY(j)=indx_ply(ndim_indx_ply + j)
3792 ENDDO
3793 ndim_indx_ply = ndim_indx_ply + npt
3794 npt_drape = drape_sh3n(ii)%NPLY_DRAPE
3795 ALLOCATE(drape_sh3n(ii)%DRAPE_PLY(npt_drape))
3796 DO j = 1,npt_drape
3797 nslice = islice(ndim_slice + j)
3798 drape_sh3n(ii)%DRAPE_PLY(j)%NSLICE = nslice
3799 ALLOCATE(drape_sh3n(ii)%DRAPE_PLY(j)%IDRAPE(nslice,2),drape_sh3n(ii)%DRAPE_PLY(j)%RDRAPE(nslice,2))
3800 drape_sh3n(ii)%DRAPE_PLY(j)%IDRAPE = 0
3801 drape_sh3n(ii)%DRAPE_PLY(j)%RDRAPE = zero
3802 DO is= 1,nslice
3803 drape_sh3n(ii)%DRAPE_PLY(j)%IDRAPE(is,1) = itmp(ndim + is,1)
3804 drape_sh3n(ii)%DRAPE_PLY(j)%IDRAPE(is,2) = itmp(ndim + is,2)
3805 drape_sh3n(ii)%DRAPE_PLY(j)%RDRAPE(is,1) = rtmp(ndim + is,1)
3806 drape_sh3n(ii)%DRAPE_PLY(j)%RDRAPE(is,2) = rtmp(ndim + is,2)
3807 ENDDO
3808 ndim = ndim + nslice
3809 ENDDO
3810 ndim_slice = ndim_slice + npt_drape
3811 ENDDO ! NUMELTG_DRAPE
3812 DEALLOCATE( npt_elt,itmp,rtmp,nptdrp_elt,islice,indx_ply,thk)
3813 stdrape = numeltg
3814 ENDIF
3815C
3816 RETURN
3817 END
3818!||====================================================================
3819!|| ply_info_rest ../engine/source/output/restart/rdresb.F
3820!||--- called by ------------------------------------------------------
3821!|| rdresb ../engine/source/output/restart/rdresb.F
3822!||--- calls -----------------------------------------------------
3823!|| read_i_c ../common_source/tools/input_output/write_routines.c
3824!||--- uses -----------------------------------------------------
3825!|| message_mod ../engine/share/message_module/message_mod.F
3826!||====================================================================
3827 SUBROUTINE ply_info_rest(PLY_INFO)
3828C-----------------------------------------------
3829C M o d u l e s
3830C-----------------------------------------------
3831 USE message_mod
3832C-----------------------------------------------
3833C I m p l i c i t T y p e s
3834C-----------------------------------------------
3835#include "implicit_f.inc"
3836C-----------------------------------------------
3837C C o m m o n B l o c k s
3838C-----------------------------------------------
3839#include "com04_c.inc"
3840C-----------------------------------------------
3841C D u m m y A r g u m e n t s
3842C-----------------------------------------------
3843 INTEGER PLY_INFO(*)
3844C-----------------------------------------------
3845C L o c a l V a r i a b l e s
3846C-----------------------------------------------
3847C
3848C--------------------------------------
3849C READING INTEGERS
3850C--------------------------------------
3851 CALL READ_I_C(PLY_INFO,3*NUMPLY)
3852C--------------------------------------
3853C READING REALS
3854C--------------------------------------
3855C
3856 RETURN
3857 END
my_real
#define my_real
Definition cppsort.cpp:32
decompress_i_nnz
subroutine decompress_i_nnz(array, len)
Definition decompress_nnz.F:32
decompress_r_nnz
subroutine decompress_r_nnz(array, len)
Definition decompress_nnz.F:71
table_rresti
subroutine table_rresti(table)
Definition table_tools.F:163
table_rrestr
subroutine table_rrestr(table)
Definition table_tools.F:227
fillxdp
subroutine fillxdp(x, xdp, d, ddp)
Definition fillxdp.F:29
fxfluxrrest
subroutine fxfluxrrest(ibfflux, fbfflux, glob_therm)
Definition fxfluxrrest.F:36
intbuf_ini
subroutine intbuf_ini(intbuf_tab)
Definition intbuf_ini.F:39
line_ini
subroutine line_ini
Definition line_ini.F:35
min
#define min(a, b)
Definition macros.h:20
ale_connectivity_mod
Definition ale_connectivity_mod.F:225
ale_mod
Definition ale_mod.F:150
ale_mod::ale
type(ale_) ale
Definition ale_mod.F:253
aleflow_mod
Definition aleflow_mod.F:30
alefvm_mod
Definition alefvm_mod.F:88
alefvm_mod::alefvm_buffer
type(alefvm_buffer_), target alefvm_buffer
Definition alefvm_mod.F:120
alefvm_mod::alefvm_param
type(alefvm_param_), target alefvm_param
Definition alefvm_mod.F:121
cluster_mod
Definition cluster_mod.F:40
damp_mod
Definition damp_mod.F:33
damp_mod::id_damp_vrel
integer, dimension(:), allocatable id_damp_vrel
Definition damp_mod.F:41
damp_mod::fr_damp_vrel
integer, dimension(:), allocatable fr_damp_vrel
Definition damp_mod.F:42
damp_mod::ndamp_vrel
integer ndamp_vrel
Definition damp_mod.F:44
damp_mod::ndamp_vrel_rbyg
integer ndamp_vrel_rbyg
Definition damp_mod.F:45
drape_mod
Definition drape_mod.F:54
drape_mod::numeltg_drape
integer numeltg_drape
Definition drape_mod.F:92
drape_mod::scdrape
integer scdrape
Definition drape_mod.F:92
drape_mod::stdrape
integer stdrape
Definition drape_mod.F:92
drape_mod::numelc_drape
integer numelc_drape
Definition drape_mod.F:92
dsgraph_mod
Definition dsgraph_mod.F:31
ecnd_mod
Definition ecdn_mod.F:34
ecnd_mod::iadcnd
integer, dimension(:), pointer iadcnd
Definition ecdn_mod.F:49
ecnd_mod::fr_cndm
integer, dimension(:), pointer fr_cndm
Definition ecdn_mod.F:47
ecnd_mod::procncnd
integer, dimension(:), pointer procncnd
Definition ecdn_mod.F:49
ecnd_mod::icnds10
integer, dimension(:), pointer icnds10
Definition ecdn_mod.F:42
ecnd_mod::addcncnd
integer, dimension(:), pointer addcncnd
Definition ecdn_mod.F:49
ecnd_mod::iad_cndm
integer, dimension(:), pointer iad_cndm
Definition ecdn_mod.F:47
eig_mod
Definition eig_mod.F:32
eig_mod::eigibuf
integer, dimension(:), allocatable eigibuf
Definition eig_mod.F:41
eig_mod::eigipm
integer, dimension(:,:), allocatable eigipm
Definition eig_mod.F:39
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_2
integer num_opt_dt_fvmbag_2
Definition fvbag_mod.F:200
fvbag_mod::fvspmd
type(fvbag_spmd), dimension(:), allocatable fvspmd
Definition fvbag_mod.F:129
fvbag_mod::airbags_total_fvm_in_h3d
integer airbags_total_fvm_in_h3d
Definition fvbag_mod.F:191
fvbag_mod::fvmbag_input_options
type(fvmbag_input_options_), dimension(:), allocatable fvmbag_input_options
Definition fvbag_mod.F:196
fvbag_mod::num_opt_dt_fvmbag_1
integer num_opt_dt_fvmbag_1
Definition fvbag_mod.F:199
fvbag_mod::num_opt_dt_fvmbag
integer num_opt_dt_fvmbag
Definition fvbag_mod.F:197
fvbag_mod::num_opt_dt_fvmbag_0
integer num_opt_dt_fvmbag_0
Definition fvbag_mod.F:198
fvbag_mod::nfvbag
integer nfvbag
Definition fvbag_mod.F:127
fvmbag_meshcontrol_mod
Definition fvmbag_meshcontrol_mod.F:46
fvmbag_meshcontrol_mod::kmesh
integer, dimension(:), allocatable kmesh
Definition fvmbag_meshcontrol_mod.F:48
fxb_mod
Definition fxb_mod.F:32
fxb_mod::fxbipm
integer, dimension(:,:), allocatable fxbipm
Definition fxb_mod.F:39
fxb_mod::fxbnod
integer, dimension(:), allocatable fxbnod
Definition fxb_mod.F:41
fxb_mod::fxbelm
integer, dimension(:), allocatable fxbelm
Definition fxb_mod.F:41
fxb_mod::fxbgrvi
integer, dimension(:), allocatable fxbgrvi
Definition fxb_mod.F:41
h3d_mod
Definition h3d_mod.F:311
imp_dyna
Definition impbufdef_mod.F:314
imp_rest
Definition impbufdef_mod.F:422
intstamp_glob_mod
Definition intstamp_glob_mod.F:30
intstamp_glob_mod::intstamp
type(intstamp_data), dimension(:), allocatable intstamp
Definition intstamp_glob_mod.F:40
message_mod
Definition message_mod.F:1257
monvol_struct_mod
Definition monvol_struct_mod.F:66
nlocal_reg_mod
Definition nlocal_reg_mod.F:119
outmax_mod
Definition outmax_mod.F:52
outmax_mod::ifcontp2max
integer ifcontp2max
Definition outmax_mod.F:69
outmax_mod::lmax_vel
integer lmax_vel
Definition outmax_mod.F:61
outmax_mod::ifcontmax
integer ifcontmax
Definition outmax_mod.F:69
outmax_mod::ifcont2max
integer ifcont2max
Definition outmax_mod.F:69
outmax_mod::ifcontp2min
integer ifcontp2min
Definition outmax_mod.F:69
outmax_mod::lmax_nstra
integer lmax_nstra
Definition outmax_mod.F:63
outmax_mod::ifcontpmax
integer ifcontpmax
Definition outmax_mod.F:69
outmax_mod::lmax_dis
integer lmax_dis
Definition outmax_mod.F:60
outmax_mod::lmax_nsig
integer lmax_nsig
Definition outmax_mod.F:62
outmax_mod::spcont_max
integer spcont_max
Definition outmax_mod.F:68
outmax_mod::ifcont2min
integer ifcont2min
Definition outmax_mod.F:69
pinchtype_mod
Definition pinchtype_mod.F:69
pinchtype_mod::npinch
integer npinch
Definition pinchtype_mod.F:73
pload_cyl_mod
Definition pload_cyl_mod.F:35
plyxfem_mod
Definition plyxfem_mod.F:31
plyxfem_mod::plynod
type(plynods), dimension(:), allocatable plynod
Definition plyxfem_mod.F:44
plyxfem_mod::indx_ply
integer, dimension(:), allocatable indx_ply
Definition plyxfem_mod.F:60
plyxfem_mod::nplypart
integer nplypart
Definition plyxfem_mod.F:59
plyxfem_mod::idpid_ply
integer, dimension(:), allocatable idpid_ply
Definition plyxfem_mod.F:61
plyxfem_mod::plysizg
integer, dimension(:), allocatable plysizg
Definition plyxfem_mod.F:124
plyxfem_mod::plyshell
type(plyshells), dimension(:), allocatable plyshell
Definition plyxfem_mod.F:56
restmod::poin_ump
integer, dimension(:), allocatable poin_ump
Definition restart_mod.F:151
restmod::iconx
integer, dimension(:), allocatable iconx
Definition restart_mod.F:60
restmod::igrv
integer, dimension(:), allocatable, target igrv
Definition restart_mod.F:60
restmod::fr_sec
integer, dimension(:), allocatable fr_sec
Definition restart_mod.F:83
restmod::iad_rby
integer, dimension(:), allocatable iad_rby
Definition restart_mod.F:83
restmod::id_global_vois
integer, dimension(:), allocatable id_global_vois
Definition restart_mod.F:86
restmod::fr_nbedge
integer, dimension(:), allocatable fr_nbedge
Definition restart_mod.F:83
restmod::iadrbmk
integer, dimension(:), allocatable iadrbmk
Definition restart_mod.F:90
restmod::fr_mad
integer, dimension(:), allocatable fr_mad
Definition restart_mod.F:83
restmod::ibcv
integer, dimension(:), allocatable ibcv
Definition restart_mod.F:106
restmod::lagbuf
integer, dimension(:), allocatable lagbuf
Definition restart_mod.F:60
restmod::ixx
integer, dimension(:), allocatable ixx
Definition restart_mod.F:60
restmod::lpby
integer, dimension(:), allocatable, target lpby
Definition restart_mod.F:60
restmod::fr_rl
integer, dimension(:), allocatable fr_rl
Definition restart_mod.F:83
restmod::ixs
integer, dimension(:), allocatable, target ixs
Definition restart_mod.F:60
restmod::iad_rbym2
integer, dimension(:), allocatable iad_rbym2
Definition restart_mod.F:90
restmod::iad_i2m
integer, dimension(:), allocatable iad_i2m
Definition restart_mod.F:83
restmod::iad_cut
integer, dimension(:), allocatable iad_cut
Definition restart_mod.F:83
restmod::lgrav
integer, dimension(:), allocatable lgrav
Definition restart_mod.F:83
restmod::npby
integer, dimension(:), allocatable, target npby
Definition restart_mod.F:60
restmod::kxig3d
integer, dimension(:), allocatable kxig3d
Definition restart_mod.F:60
restmod::iadmv3
integer, dimension(:), allocatable iadmv3
Definition restart_mod.F:83
restmod::face_vois
integer, dimension(:), allocatable face_vois
Definition restart_mod.F:86
restmod::iframe
integer, dimension(:), pointer iframe
Definition restart_mod.F:108
restmod::madfail
integer, dimension(:), pointer madfail
Definition restart_mod.F:108
restmod::lesdvois
integer, dimension(:), allocatable lesdvois
Definition restart_mod.F:83
restmod::lnrcvois
integer, dimension(:), allocatable lnrcvois
Definition restart_mod.F:83
restmod::nodenr
integer, dimension(:), allocatable nodenr
Definition restart_mod.F:83
restmod::nativ0_sms
integer, dimension(:), allocatable nativ0_sms
Definition restart_mod.F:90
restmod::newfront
integer, dimension(:), allocatable newfront
Definition restart_mod.F:83
restmod::iadc_crkxfem
integer, dimension(:), allocatable iadc_crkxfem
Definition restart_mod.F:57
restmod::nodpor
integer, dimension(:), allocatable nodpor
Definition restart_mod.F:60
restmod::ilink
integer, dimension(:), allocatable ilink
Definition restart_mod.F:83
restmod::llink
integer, dimension(:), allocatable llink
Definition restart_mod.F:83
restmod::madclnod
integer, dimension(:), allocatable madclnod
Definition restart_mod.F:57
restmod::ipadmesh
integer, dimension(:,:), allocatable ipadmesh
Definition restart_mod.F:181
restmod::lbvel
integer, dimension(:), allocatable lbvel
Definition restart_mod.F:60
restmod::lprtsph
integer, dimension(:), allocatable lprtsph
Definition restart_mod.F:60
restmod::nbsdvois
integer, dimension(:), allocatable nbsdvois
Definition restart_mod.F:83
restmod::lnodpor
integer, dimension(:), allocatable lnodpor
Definition restart_mod.F:83
restmod::ibcr
integer, dimension(:), allocatable ibcr
Definition restart_mod.F:106
restmod::ne_nercvois
integer, dimension(:), allocatable ne_nercvois
Definition restart_mod.F:90
restmod::iadmv2
integer, dimension(:), allocatable iadmv2
Definition restart_mod.F:83
restmod::neflsw
integer, dimension(:), allocatable neflsw
Definition restart_mod.F:83
restmod::ixig3d
integer, dimension(:), allocatable ixig3d
Definition restart_mod.F:60
restmod::linale
integer, dimension(:), allocatable linale
Definition restart_mod.F:83
restmod::cluster
type(cluster_), dimension(:), allocatable cluster
Definition restart_mod.F:151
restmod::iactiv
integer, dimension(:), allocatable iactiv
Definition restart_mod.F:60
restmod::crknodiad
integer, dimension(:), allocatable crknodiad
Definition restart_mod.F:57
restmod::ne_lercvois
integer, dimension(:), allocatable ne_lercvois
Definition restart_mod.F:90
restmod::ibcslag
integer, dimension(:), allocatable ibcslag
Definition restart_mod.F:60
restmod::ibufssg_io
integer, dimension(:), allocatable ibufssg_io
Definition restart_mod.F:57
restmod::sh4tree
integer, dimension(:,:), allocatable sh4tree
Definition restart_mod.F:175
restmod::fr_lagf
integer, dimension(:), allocatable fr_lagf
Definition restart_mod.F:83
restmod::ispsym
integer, dimension(:), allocatable ispsym
Definition restart_mod.F:60
restmod::sh4trim
integer, dimension(:), allocatable sh4trim
Definition restart_mod.F:185
restmod::addcsrect
integer, dimension(:), allocatable addcsrect
Definition restart_mod.F:83
restmod::ipm
integer, dimension(:), allocatable ipm
Definition restart_mod.F:83
restmod::ipart
integer, dimension(:), allocatable, target ipart
Definition restart_mod.F:60
restmod::iadwal
integer, dimension(:), allocatable iadwal
Definition restart_mod.F:83
restmod::fr_nor
integer, dimension(:), allocatable fr_nor
Definition restart_mod.F:83
restmod::isphio
integer, dimension(:), allocatable isphio
Definition restart_mod.F:60
restmod::fr_i18
integer, dimension(:), allocatable fr_i18
Definition restart_mod.F:83
restmod::ipari
integer, dimension(:), allocatable, target ipari
Definition restart_mod.F:60
restmod::fr_rbym2
integer, dimension(:), allocatable fr_rbym2
Definition restart_mod.F:90
restmod::igaup
integer, dimension(:), allocatable igaup
Definition restart_mod.F:83
restmod::iad_rbm
integer, dimension(:), allocatable iad_rbm
Definition restart_mod.F:83
restmod::iskew_ply
integer, dimension(:), allocatable iskew_ply
Definition restart_mod.F:90
restmod::nnflsw
integer, dimension(:), allocatable nnflsw
Definition restart_mod.F:83
restmod::nercvois
integer, dimension(:), allocatable nercvois
Definition restart_mod.F:83
restmod::failwave
type(failwave_str_) failwave
Definition restart_mod.F:155
restmod::ispcond
integer, dimension(:), allocatable ispcond
Definition restart_mod.F:60
restmod::ibordnode
integer, dimension(:), allocatable ibordnode
Definition restart_mod.F:83
restmod::sh3trim
integer, dimension(:), allocatable sh3trim
Definition restart_mod.F:185
restmod::iadrbk
integer, dimension(:), allocatable iadrbk
Definition restart_mod.F:83
restmod::iecran
integer, dimension(:), allocatable iecran
Definition restart_mod.F:57
restmod::iedgesh
integer, dimension(:), allocatable, target iedgesh
Definition restart_mod.F:83
restmod::fr_i2m
integer, dimension(:), allocatable fr_i2m
Definition restart_mod.F:83
restmod::ixt
integer, dimension(:), allocatable ixt
Definition restart_mod.F:60
restmod::lnlink
integer, dimension(:), allocatable lnlink
Definition restart_mod.F:60
restmod::ibftemp
integer, dimension(:), allocatable ibftemp
Definition restart_mod.F:106
restmod::ibfv
integer, dimension(:), allocatable ibfv
Definition restart_mod.F:60
restmod::iaccp
integer, dimension(:), allocatable iaccp
Definition restart_mod.F:83
restmod::iel_crkxfem
integer, dimension(:), allocatable, target iel_crkxfem
Definition restart_mod.F:83
restmod::iadi18
integer, dimension(:), allocatable iadi18
Definition restart_mod.F:83
restmod::lsegcom
integer, dimension(:), allocatable lsegcom
Definition restart_mod.F:83
restmod::iskwp_l
integer, dimension(:), allocatable iskwp_l
Definition restart_mod.F:90
restmod::inod_pxfem
integer, dimension(:), allocatable inod_pxfem
Definition restart_mod.F:83
restmod::dd_r2r_elem
integer, dimension(:), allocatable dd_r2r_elem
Definition restart_mod.F:57
restmod::kloadpinter
integer, dimension(:), allocatable kloadpinter
Definition restart_mod.F:60
restmod::ixr
integer, dimension(:), allocatable ixr
Definition restart_mod.F:60
restmod::sh3tree
integer, dimension(:,:), allocatable sh3tree
Definition restart_mod.F:175
restmod::lonfsph
integer, dimension(:), allocatable lonfsph
Definition restart_mod.F:60
restmod::madnod
integer, dimension(:), pointer madnod
Definition restart_mod.F:108
restmod::madidx
integer, dimension(:), allocatable madidx
Definition restart_mod.F:57
restmod::iexlnk
integer, dimension(:), allocatable iexlnk
Definition restart_mod.F:60
restmod::ixtg
integer, dimension(:), allocatable, target ixtg
Definition restart_mod.F:60
restmod::lpbyl
integer, dimension(:), pointer lpbyl
Definition restart_mod.F:108
restmod::nnlink
integer, dimension(:), allocatable nnlink
Definition restart_mod.F:60
restmod::ibcl
integer, dimension(:), allocatable, target ibcl
Definition restart_mod.F:60
restmod::adsky_crkxfem
integer, dimension(:), allocatable adsky_crkxfem
Definition restart_mod.F:57
restmod::madprt
integer, dimension(:), pointer madprt
Definition restart_mod.F:108
restmod::fr_rbym
integer, dimension(:), allocatable fr_rbym
Definition restart_mod.F:90
restmod::monvol
integer, dimension(:), allocatable monvol
Definition restart_mod.F:60
restmod::ifill
integer, dimension(:), allocatable ifill
Definition restart_mod.F:60
restmod::kxfenod2elc
integer, dimension(:), allocatable kxfenod2elc
Definition restart_mod.F:83
restmod::iskwp
integer, dimension(:), allocatable iskwp
Definition restart_mod.F:83
restmod::isensp
integer, dimension(:), allocatable isensp
Definition restart_mod.F:83
restmod::fr_rbe2
integer, dimension(:), allocatable fr_rbe2
Definition restart_mod.F:83
restmod::irbe2
integer, dimension(:), allocatable irbe2
Definition restart_mod.F:60
restmod::inod_crkxfem
integer, dimension(:), allocatable inod_crkxfem
Definition restart_mod.F:83
restmod::nporgeo
integer, dimension(:), allocatable nporgeo
Definition restart_mod.F:83
restmod::procne_crkxfem
integer, dimension(:), allocatable procne_crkxfem
Definition restart_mod.F:57
restmod::fr_rbm2
integer, dimension(:), allocatable fr_rbm2
Definition restart_mod.F:83
restmod::kxsp
integer, dimension(:), allocatable kxsp
Definition restart_mod.F:60
restmod::neth
integer, dimension(:), allocatable neth
Definition restart_mod.F:60
restmod::enrtag
integer, dimension(:), allocatable enrtag
Definition restart_mod.F:83
restmod::nodlevxf
integer, dimension(:), allocatable nodlevxf
Definition restart_mod.F:83
restmod::loadpinter
integer, dimension(:), allocatable loadpinter
Definition restart_mod.F:60
restmod::elcutc
integer, dimension(:), allocatable elcutc
Definition restart_mod.F:83
restmod::fr_ll
integer, dimension(:), allocatable fr_ll
Definition restart_mod.F:83
restmod::iad_sec
integer, dimension(:), allocatable iad_sec
Definition restart_mod.F:83
restmod::nsensp
integer, dimension(:), allocatable nsensp
Definition restart_mod.F:83
restmod::dd_iad
integer, dimension(:), allocatable dd_iad
Definition restart_mod.F:60
restmod::gjbufi
integer, dimension(:), allocatable gjbufi
Definition restart_mod.F:60
restmod::madsh3
integer, dimension(:), pointer madsh3
Definition restart_mod.F:108
restmod::icut
integer, dimension(:), allocatable icut
Definition restart_mod.F:83
restmod::fr_cj
integer, dimension(:), allocatable fr_cj
Definition restart_mod.F:83
restmod::iskwn
integer, dimension(:), allocatable, target iskwn
Definition restart_mod.F:60
restmod::nesdvois
integer, dimension(:), allocatable nesdvois
Definition restart_mod.F:83
restmod::cne_crkxfem
integer, dimension(:), allocatable cne_crkxfem
Definition restart_mod.F:57
restmod::iloadp
integer, dimension(:), allocatable, target iloadp
Definition restart_mod.F:60
restmod::iadcj
integer, dimension(:), allocatable iadcj
Definition restart_mod.F:83
restmod::ngaup
integer, dimension(:), allocatable ngaup
Definition restart_mod.F:83
restmod::lnrbym
integer, dimension(:), allocatable lnrbym
Definition restart_mod.F:83
restmod::nod2sp
integer, dimension(:), allocatable nod2sp
Definition restart_mod.F:60
restmod::adsky_pxfem
integer, dimension(:), allocatable adsky_pxfem
Definition restart_mod.F:57
restmod::nodglobxfe
integer, dimension(:), allocatable nodglobxfe
Definition restart_mod.F:83
restmod::weight_rm
integer, dimension(:), allocatable weight_rm
Definition restart_mod.F:57
restmod::ixp
integer, dimension(:), allocatable ixp
Definition restart_mod.F:60
restmod::laccelm
integer, dimension(:), allocatable laccelm
Definition restart_mod.F:60
restmod::nom_opt
integer, dimension(:), allocatable, target nom_opt
Definition restart_mod.F:60
restmod::iad_rbe2
integer, dimension(:), allocatable iad_rbe2
Definition restart_mod.F:83
restmod::bufgeo
double precision, dimension(:), allocatable bufgeo
Definition restart_mod.F:116
restmod::madsol
integer, dimension(:), pointer madsol
Definition restart_mod.F:108
restmod::fasolfr
integer, dimension(:), allocatable fasolfr
Definition restart_mod.F:83
restmod::iadi2
integer, dimension(:), allocatable iadi2
Definition restart_mod.F:83
restmod::npc
integer, dimension(:), allocatable, target npc
Definition restart_mod.F:60
restmod::igeo
integer, dimension(:), allocatable igeo
Definition restart_mod.F:83
restmod::ibmpc
integer, dimension(:), allocatable, target ibmpc
Definition restart_mod.F:60
restmod::ixtg1
integer, dimension(:), allocatable ixtg1
Definition restart_mod.F:60
restmod::fr_mv
integer, dimension(:), allocatable fr_mv
Definition restart_mod.F:83
restmod::ims
integer, dimension(:), allocatable ims
Definition restart_mod.F:60
restmod::fr_edge
integer, dimension(:), allocatable fr_edge
Definition restart_mod.F:83
restmod::lercvois
integer, dimension(:), allocatable lercvois
Definition restart_mod.F:83
restmod::iadrl
integer, dimension(:), allocatable iadrl
Definition restart_mod.F:83
restmod::addcni2
integer, dimension(:), allocatable addcni2
Definition restart_mod.F:83
restmod::lbcscyc
integer, dimension(:), allocatable lbcscyc
Definition restart_mod.F:57
restmod::fr_rby2
integer, dimension(:), allocatable fr_rby2
Definition restart_mod.F:83
restmod::iad_frnor
integer, dimension(:), allocatable iad_frnor
Definition restart_mod.F:83
restmod::iadc_pxfem
integer, dimension(:), allocatable iadc_pxfem
Definition restart_mod.F:57
restmod::ne_lesdvois
integer, dimension(:), allocatable ne_lesdvois
Definition restart_mod.F:90
restmod::ibvel
integer, dimension(:), allocatable ibvel
Definition restart_mod.F:60
restmod::iadrbm
integer, dimension(:), allocatable iadrbm
Definition restart_mod.F:83
restmod::lrivet
integer, dimension(:), allocatable lrivet
Definition restart_mod.F:60
restmod::npbyl
integer, dimension(:), pointer npbyl
Definition restart_mod.F:108
restmod::iad_edge
integer, dimension(:), allocatable iad_edge
Definition restart_mod.F:83
restmod::icfield
integer, dimension(:), allocatable, target icfield
Definition restart_mod.F:60
restmod::lgauge
integer, dimension(:), allocatable lgauge
Definition restart_mod.F:60
restmod::nstrf
integer, dimension(:), allocatable nstrf
Definition restart_mod.F:60
restmod::ibcscyc
integer, dimension(:), allocatable ibcscyc
Definition restart_mod.F:57
restmod::procnor
integer, dimension(:), allocatable procnor
Definition restart_mod.F:83
restmod::tag_skins6
integer, dimension(:), allocatable tag_skins6
Definition restart_mod.F:57
restmod::fr_rby
integer, dimension(:), allocatable fr_rby
Definition restart_mod.F:83
restmod::iad_rbm2
integer, dimension(:), allocatable iad_rbm2
Definition restart_mod.F:83
restmod::irbym
integer, dimension(:), allocatable irbym
Definition restart_mod.F:83
restmod::ixsp
integer, dimension(:,:), allocatable ixsp
Definition restart_mod.F:81
restmod::madsh4
integer, dimension(:), pointer madsh4
Definition restart_mod.F:108
restmod::iparg
integer, dimension(:), allocatable iparg
Definition restart_mod.F:60
restmod::ixq
integer, dimension(:), allocatable ixq
Definition restart_mod.F:60
restmod::iedge
integer, dimension(:), allocatable iedge
Definition restart_mod.F:83
restmod::ibfflux
integer, dimension(:), allocatable ibfflux
Definition restart_mod.F:106
restmod::nodedge
integer, dimension(:), allocatable nodedge
Definition restart_mod.F:83
restmod::ilas
integer, dimension(:), allocatable ilas
Definition restart_mod.F:60
restmod::iad_rby2
integer, dimension(:), allocatable iad_rby2
Definition restart_mod.F:83
restmod::icode_ply
integer, dimension(:), allocatable icode_ply
Definition restart_mod.F:90
restmod::lloadp
integer, dimension(:), allocatable lloadp
Definition restart_mod.F:83
restmod::fr_cut
integer, dimension(:), allocatable fr_cut
Definition restart_mod.F:83
restmod::segquadfr
integer, dimension(:), allocatable segquadfr
Definition restart_mod.F:83
restmod::procni2
integer, dimension(:), allocatable procni2
Definition restart_mod.F:83
restmod::lcfield
integer, dimension(:), allocatable lcfield
Definition restart_mod.F:83
restmod::kxx
integer, dimension(:), allocatable kxx
Definition restart_mod.F:60
restmod::nskwp
integer, dimension(:), allocatable nskwp
Definition restart_mod.F:83
restmod::ne_nesdvois
integer, dimension(:), allocatable ne_nesdvois
Definition restart_mod.F:90
restmod::iad_rbym
integer, dimension(:), allocatable iad_rbym
Definition restart_mod.F:90
restmod::nom_sect
integer, dimension(:), allocatable nom_sect
Definition restart_mod.F:60
restmod::npsegcom
integer, dimension(:), allocatable npsegcom
Definition restart_mod.F:83
restmod::nloc_dmg
type(nlocal_str_) nloc_dmg
Definition restart_mod.F:156
restmod::rg_cut
integer, dimension(:), allocatable rg_cut
Definition restart_mod.F:83
restmod::iel_pxfem
integer, dimension(:), allocatable iel_pxfem
Definition restart_mod.F:83
restmod::lrbe2
integer, dimension(:), allocatable lrbe2
Definition restart_mod.F:60
restmod::fr_rbm
integer, dimension(:), allocatable fr_rbm
Definition restart_mod.F:83
restmod::llagf
integer, dimension(:), allocatable llagf
Definition restart_mod.F:83
restmod::nbrcvois
integer, dimension(:), allocatable nbrcvois
Definition restart_mod.F:83
restmod::lnsdvois
integer, dimension(:), allocatable lnsdvois
Definition restart_mod.F:83
restmod::tab_ump
integer, dimension(:,:), allocatable tab_ump
Definition restart_mod.F:152
restmod::naccp
integer, dimension(:), allocatable naccp
Definition restart_mod.F:83
restmod::procne_pxfem
integer, dimension(:), allocatable procne_pxfem
Definition restart_mod.F:57
restmod::ithvar
integer, dimension(:), allocatable ithvar
Definition restart_mod.F:60
restmod::dd_r2r
integer, dimension(:), allocatable dd_r2r
Definition restart_mod.F:83
restmod::iadll
integer, dimension(:), allocatable iadll
Definition restart_mod.F:83
restmod::icontact
integer, dimension(:), allocatable icontact
Definition restart_mod.F:83
restmod::ipart_state
integer, dimension(:), allocatable ipart_state
Definition restart_mod.F:60
scratchfile_mod
Definition scratchfiles_mod.F:34
scratchfile_mod::use_ieigm
integer use_ieigm
Definition scratchfiles_mod.F:37
scratchfile_mod::len_ifxm_fn
integer len_ifxm_fn
Definition scratchfiles_mod.F:38
scratchfile_mod::use_ifxm
integer use_ifxm
Definition scratchfiles_mod.F:37
scratchfile_mod::len_ifxs_fn
integer len_ifxs_fn
Definition scratchfiles_mod.F:38
scratchfile_mod::ifxm_fn
character(len=10192) ifxm_fn
Definition scratchfiles_mod.F:35
scratchfile_mod::ieigm_fn
character(len=10192) ieigm_fn
Definition scratchfiles_mod.F:35
scratchfile_mod::len_ieigm_fn
integer len_ieigm_fn
Definition scratchfiles_mod.F:38
scratchfile_mod::ifxs_fn
character(len=10192) ifxs_fn
Definition scratchfiles_mod.F:35
scratchfile_mod::use_ifxs
integer use_ifxs
Definition scratchfiles_mod.F:37
seatbelt_mod
Definition seatbelt_mod.F:55
seatbelt_mod::nseatbelt_th_proc
integer nseatbelt_th_proc
Definition seatbelt_mod.F:185
seatbelt_mod::retractor
type(retractor_struct), dimension(:), allocatable retractor
Definition seatbelt_mod.F:155
seatbelt_mod::anchor_remote_send
type(seatbelt_remote_nodes_struct) anchor_remote_send
Definition seatbelt_mod.F:183
seatbelt_mod::anchor_remote
type(seatbelt_remote_nodes_struct) anchor_remote
Definition seatbelt_mod.F:183
seatbelt_mod::seatbelt_th_exch
type(seatbelt_th_exch_struct), dimension(:), allocatable seatbelt_th_exch
Definition seatbelt_mod.F:194
seatbelt_mod::slipring
type(slipring_struct), dimension(:), allocatable slipring
Definition seatbelt_mod.F:110
sph_mod
Definition sph_mod.F:31
sph_mod::sph2sol
integer, dimension(:), allocatable sph2sol
Definition sph_mod.F:34
sph_mod::irst
integer, dimension(:), allocatable irst
Definition sph_mod.F:34
sph_mod::sol2sph_typ
integer, dimension(:), allocatable sol2sph_typ
Definition sph_mod.F:38
sph_mod::sol2sph
integer, dimension(:), allocatable sol2sph
Definition sph_mod.F:34
stack_mod
Definition stack_mod.F:75
stack_mod::ply_info
integer, dimension(:,:), allocatable ply_info
Definition stack_mod.F:133
state_mod
Definition state_mod.F:54
table_glob_mod
Definition table_glob_mod.F:34
table_glob_mod::table
type(ttable), dimension(:), allocatable table
Definition table_glob_mod.F:44
time_mod
Definition time_mod.F:34
time_mod::global_comp_time
type(time_type) global_comp_time
Definition time_mod.F:56
unitab_mod
Definition unitab_mod.F:517
user_windows_mod
Definition user_windows_mod.F:48
xfem2vars_mod
Definition xfem2vars_mod.F:32
xfem2vars_mod::crksky
type(xfem_sky_), dimension(:), allocatable crksky
Definition xfem2vars_mod.F:44
xfem2vars_mod::xedge4n
integer, dimension(:,:), allocatable xedge4n
Definition xfem2vars_mod.F:39
xfem2vars_mod::indx_crk
integer, dimension(:), allocatable indx_crk
Definition xfem2vars_mod.F:38
xfem2vars_mod::crkedge
type(xfem_edge_), dimension(:), allocatable crkedge
Definition xfem2vars_mod.F:42
xfem2vars_mod::xedge3n
integer, dimension(:,:), allocatable xedge3n
Definition xfem2vars_mod.F:40
cp_real
subroutine cp_real(n, x, xc)
Definition produt_v.F:871
r_bufbric_22
subroutine r_bufbric_22()
Definition r_bufbric_22.F:39
rdcomr
subroutine rdcomr(lmxvrea, tabvrea, dynain_data, inter_parameters, output, glob_therm)
Definition rdcomm.F:1864
alelag_rrest
subroutine alelag_rrest()
Definition rdresb.F:3400
plyxfem_ravuply
subroutine plyxfem_ravuply()
Definition rdresb.F:3363
thcrrest
subroutine thcrrest(mcp, temp)
Definition rdresb.F:2644
imprrest
subroutine imprrest(nimpr)
Definition rdresb.F:2124
drape_rrest
subroutine drape_rrest(drape_sh4n, drape_sh3n, drapeg)
Definition rdresb.F:3673
rthbufrest
subroutine rthbufrest(rthbuf, srthbuf)
Definition rdresb.F:3433
knotlocpcrest
subroutine knotlocpcrest(knotlocpc)
Definition rdresb.F:3515
cfieldrest
subroutine cfieldrest(cfield)
Definition rdresb.F:3067
knotlocelrest
subroutine knotlocelrest(knotlocel)
Definition rdresb.F:3556
stack_rrest
subroutine stack_rrest(igeo, geo, pm)
Definition rdresb.F:3635
fvrrest
subroutine fvrrest(monvol)
Definition rdresb.F:2193
eigrrest
subroutine eigrrest()
Definition rdresb.F:1831
loadprest
subroutine loadprest(loadp)
Definition rdresb.F:3110
dsrrest
subroutine dsrrest(graphe)
Definition rdresb.F:1919
knotrest
subroutine knotrest(knot)
Definition rdresb.F:3474
nitscherrest
subroutine nitscherrest(forneqs)
Definition rdresb.F:2689
wigerest
subroutine wigerest(wige)
Definition rdresb.F:3596
fxtemprrest
subroutine fxtemprrest(ibftemp, fbftemp, glob_therm)
Definition rdresb.F:2843
fxbresb
subroutine fxbresb()
Definition rdresb.F:1693
ply_info_rest
subroutine ply_info_rest(ply_info)
Definition rdresb.F:3828
rmatrrest
subroutine rmatrrest(rbym, irbym, lnrbym, weight)
Definition rdresb.F:2894
rdresb
subroutine rdresb(af, iaf, lvarrea, irunn, multi_fvm, h3d_data, pinch_data, ale_connectivity, t_monvol, sensors, ebcs_tab, dynain_data, user_windows, output, interfaces, loads, mat_elem, python, iflow, skews, rflow, liflow, lrflow, impl_s0, forneqs, unitab, stack, drape_sh4n, drape_sh3n, drapeg, ndrape, glob_therm, pblast, element, nodes, rbe3, rwall)
Definition rdresb.F:166
convrrest
subroutine convrrest(ibcv, fconv, glob_therm)
Definition rdresb.F:2737
plyxfem_rrestanim
subroutine plyxfem_rrestanim()
Definition rdresb.F:2997
plyxfem_rrest
subroutine plyxfem_rrest(ms_layer, zi_layer, inod, iel, icode, iskew, msz2)
Definition rdresb.F:2945
nfrrest
subroutine nfrrest(iflow, rflow)
Definition rdresb.F:2073
crkxfem_rrest
subroutine crkxfem_rrest(inod_crkxfem, iel_crkxfem, nodlevxf)
Definition rdresb.F:3151
crkxfem_rrestanim
subroutine crkxfem_rrestanim()
Definition rdresb.F:3186
radiarrest
subroutine radiarrest(ibcr, fradia, glob_therm)
Definition rdresb.F:2790
read_cluster
subroutine read_cluster(cluster)
Definition read_cluster.F:34
read_dpdb
subroutine read_dpdb(a, n)
Definition read_db.F:253
read_db
subroutine read_db(a, n)
Definition read_db.F:90
read_elgroup_param
subroutine read_elgroup_param(mat_elem)
Definition read_elgroup_param.F:33
read_failwave
subroutine read_failwave(failwave)
Definition read_failwave.F:33
read_joint
subroutine read_joint()
Definition read_joint.F:34
read_matparam
subroutine read_matparam(mat_elem)
Definition read_matparam.F:46
read_nloc_struct
subroutine read_nloc_struct(nloc_dmg)
Definition read_nloc_struct.F:34
read_pcyl
subroutine read_pcyl(loads)
Definition read_pcyl.F:34
read_sensors
subroutine read_sensors(sensors, python)
Definition read_sensors.F:40
read_th_restart
subroutine read_th_restart(th)
Definition read_th_restart.F:34
redsqi
subroutine redsqi(ia, l, iform)
Definition redsqi.F:39
redsqr
subroutine redsqr(a, l, iform)
Definition redsqr.F:39
spmd_comm_split
subroutine spmd_comm_split(color, key, subcomm, rank, size_l)
Definition spmd_comm_split.F:33
spmd_fvb_comm_pattern
subroutine spmd_fvb_comm_pattern(ifv)
Definition spmd_fvb_comm_pattern.F:33
spmd_initfi
subroutine spmd_initfi(ipari, iflag, h3d_data, parameters, idt_therm, intheat)
Definition spmd_i7tool.F:3412
spmd_inisph
subroutine spmd_inisph(iflg)
Definition spmd_sptool.F:37
elbuf_ini
subroutine elbuf_ini(elbuf_tab, mat_param, mlaw_tag, prop_tag, fail_tag, igeo, ipm, iparg, ipart, ipartsp, ixs, ixq, ixc, ixtg, flag_xfem, ipartig3d, stack, igeo_stack, ixt, ixp, ixr, kxx, geo, eos_tag, istr_24, print_flag, defaults)
Definition elbuf_ini.F:45
group_ini
subroutine group_ini(igrnod, igrbric, igrquad, igrsh4n, igrsh3n, igrtruss, igrbeam, igrspring, igrpart)
Definition group_ini.F:33
subset_ini
subroutine subset_ini(subset)
Definition subset_ini.F:30
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
drape_mod::drape_
Definition drape_mod.F:71
drape_mod::drapeg_
Definition drape_mod.F:79
dsgraph_mod::prgraph
Definition dsgraph_mod.F:45
state_mod::dynain_database
Definition state_mod.F:64
unitab_mod::unit_type_
Definition unitab_mod.F:524
user_windows_mod::user_windows_
Definition user_windows_mod.F:58
surf_ini
subroutine surf_ini
Definition surf_ini.F:35
my_getpid
#define my_getpid
Definition tmpenv_c.c:47
tmpenvf
void tmpenvf(char *tmpdir, int *tmplen)
Definition tmpenv_c.c:149
userwi_read
subroutine userwi_read(user_windows, ispmd, nspmd, numnod)
Definition userwindow_interface_routines.F:35
read_i_c
void read_i_c(int *w, int *len)
Definition write_routines.c:805
cur_fil_c
void cur_fil_c(int *nf)
Definition write_routines.c:477
read_units
subroutine read_units(unitab)
Definition write_units.F:92