OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
wrcomm.F File Reference
#include "implicit_f.inc"
#include "com06_c.inc"
#include "com08_c.inc"
#include "lagmult.inc"
#include "remesh_c.inc"
#include "scr18_c.inc"
#include "sms_c.inc"
#include "sphcom.inc"
#include "inter22.inc"
#include "tablen_c.inc"
#include "ddspmd_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine wrcomr (lmxvrea, lvarrea, tabvrea, ms_nd, p, inter_parameters, unitab, glob_therm, output)

Function/Subroutine Documentation

◆ wrcomr()

subroutine wrcomr ( integer lmxvrea,
integer lvarrea,
tabvrea,
ms_nd,
integer p,
type (parameters_), intent(in) inter_parameters,
type(unit_type_), intent(in) unitab,
type(glob_therm_), intent(in) glob_therm,
type(output_), intent(in) output )

Definition at line 30 of file wrcomm.F.

32C-----------------------------------------------
33C M o d u l e s
34C-----------------------------------------------
35 USE output_mod
36 USE sph_mod, ONLY : xi_res, yi_res, zi_res
37 USE ale_mod
38 USE parameters_mod
39 USE unitab_mod
40 USE multimat_param_mod , ONLY : m51_iflg6
41 use glob_therm_mod
42 USE output_mod , ONLY : output_
43C-----------------------------------------------
44C I m p l i c i t T y p e s
45C-----------------------------------------------
46#include "implicit_f.inc"
47C-----------------------------------------------
48C D u m m y A r g u m e n t s
49C-----------------------------------------------
50 INTEGER LMXVREA,LVARREA, I
51 my_real tabvrea(*),ms_nd
52 TYPE(UNIT_TYPE_), INTENT(IN) :: UNITAB
53 TYPE(glob_therm_) ,intent(in) :: glob_therm
54 TYPE(OUTPUT_),INTENT(IN) :: OUTPUT
55C-----------------------------------------------
56C C o m m o n B l o c k s
57C-----------------------------------------------
58#include "com06_c.inc"
59#include "com08_c.inc"
60#include "lagmult.inc"
61#include "remesh_c.inc"
62#include "scr18_c.inc"
63#include "sms_c.inc"
64#include "sphcom.inc"
65#include "inter22.inc"
66C-----------------------------------------------
67#include "tablen_c.inc"
68#include "ddspmd_c.inc"
69C-----------------------------------------------
70C L o c a l V a r i a b l e s
71C-----------------------------------------------
72 INTEGER IVAR,J,RVARINT,P
73 TYPE (PARAMETERS_) ,INTENT(IN):: INTER_PARAMETERS
74C--------------------------------------
75C PREPARATION OF REAL VARIABLES FOR WRITING TO RFILE.
76C--------------------------------------
77C-----
78 DO IVAR=1,LMXVREA
79 TABVREA(IVAR)=ZERO
80 ENDDO
81C-----
82 TABVREA(1)=TT
83 TABVREA(2)=DT1
84 TABVREA(3)=DT2
85 TABVREA(4)=DT12
86 TABVREA(5)=DT2OLD
87 TABVREA(6)=TSTOP
88 TABVREA(7)=0 !no longer used
89 TABVREA(8)=THIS
90C-----
91 TABVREA(9)=0 !no longer used
92 TABVREA(10)=DTHIS
93 TABVREA(11)=DTFAC
94 TABVREA(12)=DTMIN
95 TABVREA(13)=VOLMIN
96 TABVREA(14)=OUTPUT%TH%WFEXT
97 TABVREA(15)=REINT
98 TABVREA(16)=0 !no longer used
99 TABVREA(17)=0 !no longer used
100 TABVREA(18)=OUTPUT%DTANIM
101 TABVREA(19)=OUTPUT%TANIM
102 TABVREA(20)=UREINT
103 TABVREA(21)=ECONTV
104 TABVREA(22)=EHOUR
105 TABVREA(23)=DTOUTP
106 TABVREA(24)=TOUTP
107 TABVREA(25)=T1S
108 TABVREA(26)=DT2S
109 TABVREA(27)=DT12S
110 TABVREA(28)=ZERO !old USREINT
111 TABVREA(29)=R2RFX1
112 TABVREA(30)=R2RFX2
113 TABVREA(31)=DAMPA
114 TABVREA(32)=DAMPB
115 TABVREA(33)=DAMPW
116C-----
117 TABVREA(34)=ALE%GRID%ALPHA
118 TABVREA(35)=ALE%GRID%GAMMA
119 TABVREA(36)=ALE%GRID%VGX
120 TABVREA(37)=ALE%GRID%VGY
121 TABVREA(38)=ALE%GRID%VGZ
122C-----
123 TABVREA(39)=SPASORT
124 TABVREA(40)=SPATRUE
125 DO I=1,9
126 TABVREA(40+I) = DTHIS1(I)
127 TABVREA(49+I) = THIS1(I)
128 ENDDO
129 TABVREA(59) = ZERO
130 TABVREA(60) = ZERO
131 TABVREA(61)=LAGM_TOL
132 TABVREA(62)=LAG_ALPH
133 TABVREA(63)=LAG_ALPHS
134 LVARREA=63
135C
136C---- Debut Check Point Restart
137C
138 DO I = 1, 42
139 TABVREA(63+I)=ZERO
140 ENDDO
141C
142C---- end check point restart
143C
144 LVARREA= 105
145 TABVREA(106)=DWMAD
146 LVARREA=106
147 TABVREA(107)=UNITAB%FAC_MASS
148 TABVREA(108)=UNITAB%FAC_LENGTH
149 TABVREA(109)=UNITAB%FAC_TIME
150 TABVREA(110)=T1SH
151 TABVREA(111) = ALE%UPWIND%UPWSM
152 RVARINT = 111
153 RVARINT = RVARINT + 1
154 TABVREA(RVARINT) = TPSREF
155C THERMAL ANALYSIS
156 RVARINT = RVARINT + 1
157 TABVREA(RVARINT) = GLOB_THERM%HEAT_STORED
158 RVARINT = RVARINT + 1
159 TABVREA(RVARINT) = GLOB_THERM%HEAT_MECA
160 RVARINT = RVARINT + 1
161 TABVREA(RVARINT) = GLOB_THERM%HEAT_CONV
162 RVARINT = RVARINT + 1
163 TABVREA(RVARINT) = GLOB_THERM%HEAT_RADIA
164 RVARINT = RVARINT + 1
165 TABVREA(RVARINT) = GLOB_THERM%HEAT_FFLUX
166C TABLEAU SOL1TNL
167 DO I=0, MAXLAW
168 RVARINT = RVARINT + 1
169 TABVREA(RVARINT) = SOL1TNL(I,1)
170 RVARINT = RVARINT + 1
171 TABVREA(RVARINT) = SOL1TNL(I,2)
172 ENDDO
173C TABLEAU SOL8TNL
174 DO I=0, MAXLAW
175 RVARINT = RVARINT + 1
176 TABVREA(RVARINT) = SOL8TNL(I,1)
177 RVARINT = RVARINT + 1
178 TABVREA(RVARINT) = SOL8TNL(I,2)
179 ENDDO
180C TABLEAU TET4TNL
181 DO I=0, MAXLAW
182 RVARINT = RVARINT + 1
183 TABVREA(RVARINT) = TET4TNL(I,1)
184 ENDDO
185C TABLEAU TET10TNL
186 DO I=0, MAXLAW
187 RVARINT = RVARINT + 1
188 TABVREA(RVARINT) = TET10TNL(I,1)
189 ENDDO
190C TABLEAU SHTNL
191 DO I=0, MAXLAW
192 DO J=0,3
193 RVARINT = RVARINT + 1
194 TABVREA(RVARINT) = SHTNL(I,J,1)
195 RVARINT = RVARINT + 1
196 TABVREA(RVARINT) = SHTNL(I,J,2)
197 RVARINT = RVARINT + 1
198 TABVREA(RVARINT) = SHTNL(I,J,3)
199 ENDDO
200 ENDDO
201C TABLEAU TRITNL
202 DO I=0, MAXLAW
203 DO J=0,3
204 RVARINT = RVARINT + 1
205 TABVREA(RVARINT) = TRITNL(I,J,1)
206 RVARINT = RVARINT + 1
207 TABVREA(RVARINT) = TRITNL(I,J,2)
208 RVARINT = RVARINT + 1
209 TABVREA(RVARINT) = TRITNL(I,J,3)
210 ENDDO
211 ENDDO
212C TABLEAU SOLTELT
213 DO I=1,8
214 RVARINT = RVARINT + 1
215 TABVREA(RVARINT) = SOLTELT(I)
216 ENDDO
217 RVARINT = RVARINT + 1
218 TABVREA(RVARINT) = TET4TELT(1)
219 RVARINT = RVARINT + 1
220 TABVREA(RVARINT) = TET10TELT(1)
221C TABLEAU SHTELT
222 DO I=1,10
223 RVARINT = RVARINT + 1
224 TABVREA(RVARINT) = SHTELT(I)
225 ENDDO
226C TABLEAU TRITELT
227 DO I=1,10
228 RVARINT = RVARINT + 1
229 TABVREA(RVARINT) = TRITELT(I)
230 ENDDO
231C
232C .sta files
233 RVARINT = RVARINT + 1
234 TABVREA(RVARINT) = DTSTAT
235 RVARINT = RVARINT + 1
236 TABVREA(RVARINT) = TSTAT
237C .dynain files
238 RVARINT = RVARINT + 1
239 TABVREA(RVARINT) = ZERO !DTDYNAIN
240 RVARINT = RVARINT + 1
241 TABVREA(RVARINT) = ZERO !TDYNAIN
242C
243 RVARINT = RVARINT + 1
244 TABVREA(RVARINT) = TADMESH
245 RVARINT = RVARINT + 1
246 TABVREA(RVARINT) = DTADMESH
247C
248 RVARINT = RVARINT + 1
249 TABVREA(RVARINT) = DTFACS
250 RVARINT = RVARINT + 1
251 TABVREA(RVARINT) = DTMINS
252 RVARINT = RVARINT + 1
253 TABVREA(RVARINT) = TOL_SMS
254C abf files
255 DO I=1,10
256 RVARINT = RVARINT + 1
257 TABVREA(RVARINT) = TABFIS(I)
258 RVARINT = RVARINT + 1
259 TABVREA(RVARINT) = DTABF(I)
260 RVARINT = RVARINT + 1
261 TABVREA(RVARINT) = DTABFWR(I)
262 RVARINT = RVARINT + 1
263 TABVREA(RVARINT) = TABFWR(I)
264 ENDDO
265C
266C
267 RVARINT = RVARINT + 1
268 TABVREA(RVARINT) = RATIO22
269C
270 RVARINT = RVARINT + 1
271 TABVREA(RVARINT) = M51_IFLG6
272C
273 RVARINT = RVARINT + 1
274 TABVREA(RVARINT) = MS_ND ! MASS Middle nodes Itet=2
275C
276 RVARINT = RVARINT + 1
277 TABVREA(RVARINT) = ZERO ! DMSI2 Itet=2
278 RVARINT = RVARINT + 1
279 TABVREA(RVARINT) = ZERO ! MSI20 Itet=2
280C
281C---- Debut Check Point Restart
282C
283C DTFAC1(51) and DTMIN1(51)
284 RVARINT = RVARINT + 1
285 TABVREA(RVARINT) = ZERO
286 RVARINT = RVARINT + 1
287 TABVREA(RVARINT) = ZERO
288C DTFAC1(52) and DTMIN1(52)
289 RVARINT = RVARINT + 1
290 TABVREA(RVARINT) = ZERO
291 RVARINT = RVARINT + 1
292 TABVREA(RVARINT) = ZERO
293C
294 !DTFAC1(102) - DTMIN1(102) - /DT/ALE
295 RVARINT = RVARINT + 1
296 TABVREA(RVARINT) = ZERO
297 RVARINT = RVARINT + 1
298 TABVREA(RVARINT) = ZERO
299C
300 RVARINT = RVARINT + 1
301 TABVREA(RVARINT) = ZERO ! DTFACS_INT
302 RVARINT = RVARINT + 1
303 TABVREA(RVARINT) = ZERO ! DTMINS_INT
304C
305C BETATE_0, EINT_0,ENCIN_0,PIMAX ,PCMAX,PINT ,PCIN ,FILtr+8,NFIRST
306 DO I = 1, 16
307 RVARINT = RVARINT + 1
308 TABVREA(RVARINT) = ZERO
309 END DO
310C
311 RVARINT = RVARINT + 2
312C ENCONTD , ECONT_CUMU
313 RVARINT = RVARINT + 1
314 TABVREA(RVARINT) = ECONTD
315 RVARINT = RVARINT + 1
316 TABVREA(RVARINT) = ECONT_CUMU
317C INITIAL MASS FROM STARTER
318 RVARINT = RVARINT + 1
319 TABVREA(RVARINT) = MASS0_START
320 RVARINT = RVARINT + 1
321 TABVREA(RVARINT) = ZERO
322C ANIM FILES STOP TIME
323 RVARINT = RVARINT + 1
324 TABVREA(RVARINT) = OUTPUT%TANIM_STOP
325C INITIAL COORDS OF SPH RESERVE
326 RVARINT = RVARINT + 1
327 TABVREA(RVARINT) = XI_RES
328 RVARINT = RVARINT + 1
329 TABVREA(RVARINT) = YI_RES
330 RVARINT = RVARINT + 1
331 TABVREA(RVARINT) = ZI_RES
332C
333 RVARINT = RVARINT + 1
334 TABVREA(RVARINT) = INTER_PARAMETERS%DT_STIFINT
335C KINETIC ENERGY
336 RVARINT = RVARINT + 1
337 TABVREA(RVARINT) = ZERO !ENCIN
338 RVARINT = RVARINT + 1
339 TABVREA(RVARINT) = ZERO !ENROT
340C
341C---- end check point restart
342C
343 LVARREA=RVARINT
344C----------------------------------------------------------------
345C warning to radioss developers : Check LVARREA < LMXVREA
346C----------------------------------------------------------------
347 IF(LVARREA>LMXVREA)STOP 123
348C
349 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
ale_mod
Definition ale_mod.F:150
parameters_mod
Definition parameters_mod.F:49
sph_mod
Definition sph_mod.F:31
unitab_mod
Definition unitab_mod.F:517