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