OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
wrcomm.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!|| wrcomr ../starter/source/restart/ddsplit/wrcomm.F
25!||--- called by ------------------------------------------------------
26!|| ddsplit ../starter/source/restart/ddsplit/ddsplit.F
27!||--- uses -----------------------------------------------------
28!||====================================================================
29 SUBROUTINE wrcomr(LMXVREA,LVARREA,TABVREA,MS_ND,P,INTER_PARAMETERS,
30 . UNITAB ,GLOB_THERM, OUTPUT)
31C-----------------------------------------------
32C M o d u l e s
33C-----------------------------------------------
34 USE anim_mod
35 USE outputs_mod
36 USE ale_mod
37 USE parameters_mod
38 USE unitab_mod
39 USE multimat_param_mod , ONLY : m51_iflg6
40 use glob_therm_mod
41 USE output_mod , ONLY : output_
42C-----------------------------------------------
43C I m p l i c i t T y p e s
44C-----------------------------------------------
45#include "implicit_f.inc"
46C-----------------------------------------------
47C D u m m y A r g u m e n t s
48C-----------------------------------------------
49 INTEGER LMXVREA,LVARREA, I
50 my_real tabvrea(*),ms_nd
51 TYPE(unit_type_), INTENT(IN) :: UNITAB
52 TYPE(glob_therm_) ,intent(in) :: glob_therm
53 TYPE(output_),INTENT(IN) :: OUTPUT
54C-----------------------------------------------
55C C o m m o n B l o c k s
56C-----------------------------------------------
57#include "com06_c.inc"
58#include "com08_c.inc"
59#include "lagmult.inc"
60#include "remesh_c.inc"
61#include "scr18_c.inc"
62#include "sms_c.inc"
63#include "sphcom.inc"
64#include "inter22.inc"
65C-----------------------------------------------
66#include "tablen_c.inc"
67#include "ddspmd_c.inc"
68C-----------------------------------------------
69C L o c a l V a r i a b l e s
70C-----------------------------------------------
71 INTEGER IVAR,J,RVARINT,P
72 TYPE (PARAMETERS_) ,INTENT(IN):: INTER_PARAMETERS
73C--------------------------------------
74C PREPARATION DES VARIABLES REELLES POUR ECRITURE SUR RFILE.
75C--------------------------------------
76C-----
77 DO ivar=1,lmxvrea
78 tabvrea(ivar)=zero
79 ENDDO
80C-----
81 tabvrea(1)=tt
82 tabvrea(2)=dt1
83 tabvrea(3)=dt2
84 tabvrea(4)=dt12
85 tabvrea(5)=dt2old
86 tabvrea(6)=tstop
87 tabvrea(7)=0 !no longer used
88 tabvrea(8)=this
89C-----
90 tabvrea(9)=0 !no longer used
91 tabvrea(10)=dthis
92 tabvrea(11)=dtfac
93 tabvrea(12)=dtmin
94 tabvrea(13)=volmin
95 tabvrea(14)=output%TH%WFEXT
96 tabvrea(15)=reint
97 tabvrea(16)=0 !no longer used
98 tabvrea(17)=0 !no longer used
99 tabvrea(18)=dtanim
100 tabvrea(19)=tanim
101 tabvrea(20)=ureint
102 tabvrea(21)=econtv
103 tabvrea(22)=ehour
104 tabvrea(23)=dtoutp
105 tabvrea(24)=toutp
106 tabvrea(25)=t1s
107 tabvrea(26)=dt2s
108 tabvrea(27)=dt12s
109 tabvrea(28)=zero !old USREINT
110 tabvrea(29)=r2rfx1
111 tabvrea(30)=r2rfx2
112 tabvrea(31)=dampa
113 tabvrea(32)=dampb
114 tabvrea(33)=dampw
115C-----
116 tabvrea(34)=ale%GRID%ALPHA
117 tabvrea(35)=ale%GRID%GAMMA
118 tabvrea(36)=ale%GRID%VGX
119 tabvrea(37)=ale%GRID%VGY
120 tabvrea(38)=ale%GRID%VGZ
121C-----
122 tabvrea(39)=spasort
123 tabvrea(40)=spatrue
124 DO i=1,9
125 tabvrea(40+i) = dthis1(i)
126 tabvrea(49+i) = this1(i)
127 ENDDO
128 tabvrea(59) = zero
129 tabvrea(60) = zero
130 tabvrea(61)=lagm_tol
131 tabvrea(62)=lag_alph
132 tabvrea(63)=lag_alphs
133 lvarrea=63
134C
135C---- debut check point restart
136C
137 DO i = 1, 42
138 tabvrea(63+i)=zero
139 ENDDO
140C
141C---- fin check point restart
142C
143 lvarrea= 105
144 tabvrea(106)=dwmad
145 lvarrea=106
146 tabvrea(107)=unitab%FAC_MASS
147 tabvrea(108)=unitab%FAC_LENGTH
148 tabvrea(109)=unitab%FAC_TIME
149 tabvrea(110)=t1sh
150 tabvrea(111) = ale%UPWIND%UPWSM
151 rvarint = 111
152 rvarint = rvarint + 1
153 tabvrea(rvarint) = tpsref
154C THERMAL ANALYSIS
155 rvarint = rvarint + 1
156 tabvrea(rvarint) = glob_therm%HEAT_STORED
157 rvarint = rvarint + 1
158 tabvrea(rvarint) = glob_therm%HEAT_MECA
159 rvarint = rvarint + 1
160 tabvrea(rvarint) = glob_therm%HEAT_CONV
161 rvarint = rvarint + 1
162 tabvrea(rvarint) = glob_therm%HEAT_RADIA
163 rvarint = rvarint + 1
164 tabvrea(rvarint) = glob_therm%HEAT_FFLUX
165C TABLEAU SOL1TNL
166 DO i=0, maxlaw
167 rvarint = rvarint + 1
168 tabvrea(rvarint) = sol1tnl(i,1)
169 rvarint = rvarint + 1
170 tabvrea(rvarint) = sol1tnl(i,2)
171 ENDDO
172C TABLEAU SOL8TNL
173 DO i=0, maxlaw
174 rvarint = rvarint + 1
175 tabvrea(rvarint) = sol8tnl(i,1)
176 rvarint = rvarint + 1
177 tabvrea(rvarint) = sol8tnl(i,2)
178 ENDDO
179C TABLEAU TET4TNL
180 DO i=0, maxlaw
181 rvarint = rvarint + 1
182 tabvrea(rvarint) = tet4tnl(i,1)
183 ENDDO
184C TABLEAU TET10TNL
185 DO i=0, maxlaw
186 rvarint = rvarint + 1
187 tabvrea(rvarint) = tet10tnl(i,1)
188 ENDDO
189C TABLEAU SHTNL
190 DO i=0, maxlaw
191 DO j=0,3
192 rvarint = rvarint + 1
193 tabvrea(rvarint) = shtnl(i,j,1)
194 rvarint = rvarint + 1
195 tabvrea(rvarint) = shtnl(i,j,2)
196 rvarint = rvarint + 1
197 tabvrea(rvarint) = shtnl(i,j,3)
198 ENDDO
199 ENDDO
200C TABLEAU TRITNL
201 DO i=0, maxlaw
202 DO j=0,3
203 rvarint = rvarint + 1
204 tabvrea(rvarint) = tritnl(i,j,1)
205 rvarint = rvarint + 1
206 tabvrea(rvarint) = tritnl(i,j,2)
207 rvarint = rvarint + 1
208 tabvrea(rvarint) = tritnl(i,j,3)
209 ENDDO
210 ENDDO
211C TABLEAU SOLTELT
212 DO i=1,8
213 rvarint = rvarint + 1
214 tabvrea(rvarint) = soltelt(i)
215 ENDDO
216 rvarint = rvarint + 1
217 tabvrea(rvarint) = tet4telt(1)
218 rvarint = rvarint + 1
219 tabvrea(rvarint) = tet10telt(1)
220C TABLEAU SHTELT
221 DO i=1,10
222 rvarint = rvarint + 1
223 tabvrea(rvarint) = shtelt(i)
224 ENDDO
225C TABLEAU TRITELT
226 DO i=1,10
227 rvarint = rvarint + 1
228 tabvrea(rvarint) = tritelt(i)
229 ENDDO
230C
231C .sta files
232 rvarint = rvarint + 1
233 tabvrea(rvarint) = dtstat
234 rvarint = rvarint + 1
235 tabvrea(rvarint) = tstat
236C .dynain files
237 rvarint = rvarint + 1
238 tabvrea(rvarint) = zero !DTDYNAIN
239 rvarint = rvarint + 1
240 tabvrea(rvarint) = zero !TDYNAIN
241C
242 rvarint = rvarint + 1
243 tabvrea(rvarint) = tadmesh
244 rvarint = rvarint + 1
245 tabvrea(rvarint) = dtadmesh
246C
247 rvarint = rvarint + 1
248 tabvrea(rvarint) = dtfacs
249 rvarint = rvarint + 1
250 tabvrea(rvarint) = dtmins
251 rvarint = rvarint + 1
252 tabvrea(rvarint) = tol_sms
253C abf files
254 DO i=1,10
255 rvarint = rvarint + 1
256 tabvrea(rvarint) = tabfis(i)
257 rvarint = rvarint + 1
258 tabvrea(rvarint) = dtabf(i)
259 rvarint = rvarint + 1
260 tabvrea(rvarint) = dtabfwr(i)
261 rvarint = rvarint + 1
262 tabvrea(rvarint) = tabfwr(i)
263 ENDDO
264C
265C
266 rvarint = rvarint + 1
267 tabvrea(rvarint) = ratio22
268C
269 rvarint = rvarint + 1
270 tabvrea(rvarint) = m51_iflg6
271C
272 rvarint = rvarint + 1
273 tabvrea(rvarint) = ms_nd ! MASS Middle nodes Itet=2
274C
275 rvarint = rvarint + 1
276 tabvrea(rvarint) = zero ! DMSI2 Itet=2
277 rvarint = rvarint + 1
278 tabvrea(rvarint) = zero ! MSI20 Itet=2
279C
280C---- debut check point restart
281C
282C DTFAC1(51) and DTMIN1(51)
283 rvarint = rvarint + 1
284 tabvrea(rvarint) = zero
285 rvarint = rvarint + 1
286 tabvrea(rvarint) = zero
287C DTFAC1(52) and DTMIN1(52)
288 rvarint = rvarint + 1
289 tabvrea(rvarint) = zero
290 rvarint = rvarint + 1
291 tabvrea(rvarint) = zero
292C
293 !DTFAC1(102) - DTMIN1(102) - /DT/ALE
294 rvarint = rvarint + 1
295 tabvrea(rvarint) = zero
296 rvarint = rvarint + 1
297 tabvrea(rvarint) = zero
298C
299 rvarint = rvarint + 1
300 tabvrea(rvarint) = zero ! DTFACS_INT
301 rvarint = rvarint + 1
302 tabvrea(rvarint) = zero ! DTMINS_INT
303C
304C BETATE_0, EINT_0,ENCIN_0,PIMAX ,PCMAX,PINT ,PCIN ,FILtr+8,NFIRST
305 DO i = 1, 16
306 rvarint = rvarint + 1
307 tabvrea(rvarint) = zero
308 END DO
309C
310 rvarint = rvarint + 2
311C ENCONTD , ECONT_CUMU
312 rvarint = rvarint + 1
313 tabvrea(rvarint) = econtd
314 rvarint = rvarint + 1
315 tabvrea(rvarint) = econt_cumu
316C INITIAL MASS FROM STARTER
317 rvarint = rvarint + 1
318 tabvrea(rvarint) = mass0_start
319 rvarint = rvarint + 1
320 tabvrea(rvarint) = zero
321C ANIM FILES STOP TIME
322 rvarint = rvarint + 1
323 tabvrea(rvarint) = tanim_stop
324C INITIAL COORDS OF SPH RESERVE
325 rvarint = rvarint + 1
326 tabvrea(rvarint) = xi_res
327 rvarint = rvarint + 1
328 tabvrea(rvarint) = yi_res
329 rvarint = rvarint + 1
330 tabvrea(rvarint) = zi_res
331C
332 rvarint = rvarint + 1
333 tabvrea(rvarint) = inter_parameters%DT_STIFINT
334C KINETIC ENERGY
335 rvarint = rvarint + 1
336 tabvrea(rvarint) = zero !ENCIN
337 rvarint = rvarint + 1
338 tabvrea(rvarint) = zero !ENROT
339C
340C---- fin check point restart
341C
342 lvarrea=rvarint
343C----------------------------------------------------------------
344C warning to radioss developers : Check LVARREA < LMXVREA
345C----------------------------------------------------------------
346 IF(lvarrea>lmxvrea)stop 123
347C
348 RETURN
349 END
my_real
#define my_real
Definition cppsort.cpp:32
ale_mod
Definition ale_mod.F:149
ale_mod::ale
type(ale_) ale
Definition ale_mod.F:249
anim_mod
Definition anim_mod.F:96
outputs_mod
Definition outputs_mod.F:59
parameters_mod
Definition parameters_mod.F:49
unitab_mod
Definition unitab_mod.F:513
wrcomr
subroutine wrcomr(lmxvrea, lvarrea, tabvrea, ms_nd, p, inter_parameters, unitab, glob_therm, output)
Definition wrcomm.F:31
unitab_mod::unit_type_
Definition unitab_mod.F:520