31
32
33
39 USE multimat_param_mod , ONLY : m51_iflg6
40 use glob_therm_mod
41 USE output_mod , ONLY : output_
42
43
44
45#include "implicit_f.inc"
46
47
48
49 INTEGER LMXVREA,LVARREA, I
51 TYPE(UNIT_TYPE_), INTENT(IN) :: UNITAB
52 TYPE(glob_therm_) ,intent(in) :: glob_therm
53 TYPE(OUTPUT_),INTENT(IN) :: OUTPUT
54
55
56
57#include "com06_c.inc"
58#include "com08_c."
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"
65
66#include "tablen_c.inc"
67#include "ddspmd_c.inc"
68
69
70
71 INTEGER IVAR,J,RVARINT,P
72 TYPE (PARAMETERS_) ,INTENT(IN):: INTER_PARAMETERS
73
74
75
76
77 DO IVAR=1,LMXVREA
78 TABVREA(IVAR)=ZERO
79 ENDDO
80
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
89
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
115
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
121
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
134
135
136
137 DO I = 1, 42
138 TABVREA(63+I)=ZERO
139 ENDDO
140
141
142
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
154
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
165
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
172
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
179
180 DO I=0, MAXLAW
181 RVARINT = RVARINT + 1
182 TABVREA(RVARINT) = TET4TNL(I,1)
183 ENDDO
184
185 DO I=0, MAXLAW
186 RVARINT = RVARINT + 1
187 TABVREA(RVARINT) = TET10TNL(I,1)
188 ENDDO
189
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
200
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
211
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)
220
221 DO I=1,10
222 RVARINT = RVARINT + 1
223 TABVREA(RVARINT) = SHTELT(I)
224 ENDDO
225
226 DO I=1,10
227 RVARINT = RVARINT + 1
228 TABVREA(RVARINT) = TRITELT(I)
229 ENDDO
230
231
232 RVARINT = RVARINT + 1
233 TABVREA(RVARINT) = DTSTAT
234 RVARINT = RVARINT + 1
235 TABVREA(RVARINT) = TSTAT
236
237 RVARINT = RVARINT + 1
238 TABVREA(RVARINT) = ZERO !DTDYNAIN
239 RVARINT = RVARINT + 1
240 TABVREA(RVARINT) = ZERO !TDYNAIN
241
242 RVARINT = RVARINT + 1
243 TABVREA(RVARINT) = TADMESH
244 RVARINT = RVARINT + 1
245 TABVREA(RVARINT) = DTADMESH
246
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
253
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
264
265
266 RVARINT = RVARINT + 1
267 TABVREA(RVARINT) = RATIO22
268
269 RVARINT = RVARINT + 1
270 TABVREA(RVARINT) = M51_IFLG6
271
272 RVARINT = RVARINT + 1
273 TABVREA(RVARINT) = MS_ND ! MASS Middle nodes Itet=2
274
275 RVARINT = RVARINT + 1
276 TABVREA(RVARINT) = ZERO ! DMSI2 Itet=2
277 RVARINT = RVARINT + 1
278 TABVREA(RVARINT) = ZERO ! MSI20 Itet=2
279
280
281
282
283 RVARINT = RVARINT + 1
284 TABVREA(RVARINT) = ZERO
285 RVARINT = RVARINT + 1
286 TABVREA(RVARINT) = ZERO
287
288 RVARINT = RVARINT + 1
289 TABVREA(RVARINT) = ZERO
290 RVARINT = RVARINT + 1
291 TABVREA(RVARINT) = ZERO
292
293 !DTFAC1(102) - DTMIN1(102) - /DT/ALE
294 RVARINT = RVARINT + 1
295 TABVREA(RVARINT) = ZERO
296 RVARINT = RVARINT + 1
297 TABVREA(RVARINT) = ZERO
298
299 RVARINT = RVARINT + 1
300 TABVREA(RVARINT) = ZERO ! DTFACS_INT
301 RVARINT = RVARINT + 1
302 TABVREA(RVARINT) = ZERO ! DTMINS_INT
303
304
305 DO I = 1, 16
306 RVARINT = RVARINT + 1
307 TABVREA(RVARINT) = ZERO
308 END DO
309
310 RVARINT = RVARINT + 2
311
312 RVARINT = RVARINT + 1
313 TABVREA(RVARINT) = ECONTD
314 RVARINT = RVARINT + 1
315 TABVREA(RVARINT) = ECONT_CUMU
316
317 RVARINT = RVARINT + 1
318 TABVREA(RVARINT) = MASS0_START
319 RVARINT = RVARINT + 1
320 TABVREA(RVARINT) = ZERO
321
322 RVARINT = RVARINT + 1
323 TABVREA(RVARINT) = TANIM_STOP
324
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
331
332 RVARINT = RVARINT + 1
333 TABVREA(RVARINT) = INTER_PARAMETERS%DT_STIFINT
334
335 RVARINT = RVARINT + 1
336 TABVREA(RVARINT) = ZERO !ENCIN
337 RVARINT = RVARINT + 1
338 TABVREA(RVARINT) = ZERO !ENROT
339
340
341
342 LVARREA=RVARINT
343
344
345
346 IF(LVARREA>LMXVREA)STOP 123
347
348 RETURN