70
71
72
73 USE intbufdef_mod
75
76
77
78#include "implicit_f.inc"
79
80
81
82 INTEGER IPARI(*), IADI2(*),WEIGHT(*),IGEO(NPROPGI,*),
83 . NODNX_SMS(*),ITAB(*)
84 INTEGER I2MSCH, I0, NIR, I2SIZE,NELTST,ITYPTST
85 INTEGER ,intent(in) :: IDT_THERM
86 my_real ,
intent(in) :: theaccfact
88 . x(*), v(*), a(*),dmast,adm(*),
89 . ms(*),in(*),ar(*),vr(*),stifn(*),stifr(*),fskyi2(*),
90 . adi(*), bufgeo(*), fsav(*),fncont(3,*),
91 . dmint2(*),sav_for_pena(*),ms_pena(*),dt2t,temp(*),
92 . fthe(*),ftheskyi2(*),condn(*),condnskyi2(*),
93 . sav_iner_poff(*),t2fac_sms(*),
94 . fncontp(3,*),ftcontp(3,*)
95
96 TYPE(INTBUF_STRUCT_) INTBUF_TAB
97 TYPE (H3D_DATABASE) :: H3D_DATA
98
99
100
101#include "com01_c.inc"
102#include "parit_c.inc"
103#include "param_c.inc"
104
105
106
107 integer
108 . nsn,nmn,nrts,nrtm,ilev,i0base,idel2,penflag,irot,noint,intth,ivisc
110 . kthe
111
112 nrts =ipari(3)
113 nrtm =ipari(4)
114 nsn =ipari(5)
115 nmn =ipari(6)
116 noint =ipari(15)
117 idel2 = ipari(17)
118 ilev =ipari(20)
119 ivisc = ipari(72)
120
121 i0base = i0
122
123 IF (iparit == 0) THEN
124 IF((ilev == 1.OR.ilev == 3).AND.iroddl /= 0)THEN
125
127 1 nsn ,nmn ,a ,intbuf_tab%IRECTM,intbuf_tab%DPARA,
128 2 intbuf_tab%MSR,intbuf_tab%NSV,intbuf_tab%IRTLM,ms ,weight ,
129 3 ar ,in ,x ,stifn ,stifr ,
130 4 ilev ,dmast ,adm ,intbuf_tab%NMAS,idel2 ,
131 5 intbuf_tab%SMAS,intbuf_tab%SINER,v ,intbuf_tab%CSTS,fsav ,
132 6 fncont ,h3d_data ,fncontp ,ftcontp )
133
134 ELSEIF(ilev == 2)THEN
135
137 1 nsn ,nmn ,a ,intbuf_tab%IRECTM,intbuf_tab%CSTS,
138 2 intbuf_tab%MSR,intbuf_tab%NSV,intbuf_tab%IRTLM,ms ,weight ,
139 3 stifn ,intbuf_tab%NMAS,x ,v ,fsav ,
140 4 fncont ,h3d_data ,intbuf_tab%CSTS_BIS ,fncontp ,ftcontp )
141
142 IF(iroddl /= 0)THEN
144 1 nsn ,nmn ,ar ,intbuf_tab%IRECTM,intbuf_tab%CSTS,
145 2 intbuf_tab%MSR,intbuf_tab%NSV,intbuf_tab%IRTLM,in ,ms ,
146 3 a ,x ,weight ,stifr ,stifn ,
147 4 sav_iner_poff,intbuf_tab%CSTS_BIS)
148 ENDIF
149 ELSEIF (ilev == 4)THEN
150
152 1 nsn ,nmn ,a ,intbuf_tab%IRECTM,intbuf_tab%CSTS,
153 2 intbuf_tab%MSR,intbuf_tab%NSV,intbuf_tab%IRTLM,ms ,weight ,
154 3 stifn ,intbuf_tab%NMAS,x ,v ,fsav ,
155 5 fncont ,intbuf_tab%CSTS_BIS,h3d_data ,fncontp,ftcontp)
156 ELSEIF (ilev == 10 .OR. ilev == 11 .OR. ilev == 12 .OR.
157 . ilev == 20 .OR. ilev == 21 .OR. ilev == 22) THEN
158
160 1 a ,ms ,stifn ,weight ,intbuf_tab%IRECTM,
161 2 intbuf_tab%NSV,intbuf_tab%MSR,intbuf_tab%IRTLM,intbuf_tab%IRUPT,intbuf_tab%CSTS,
162 3 intbuf_tab%FSM,nsn,intbuf_tab%CSTS_BIS)
163
164 IF(iroddl /= 0) THEN
166 1 nsn ,nmn ,ar ,intbuf_tab%IRECTM,intbuf_tab%CSTS,
167 2 intbuf_tab%MSR,intbuf_tab%NSV,intbuf_tab%IRTLM,in ,ms ,
168 3 a ,x ,weight ,stifr ,stifn ,
169 4 intbuf_tab%IRUPT,ilev,intbuf_tab%CSTS_BIS)
170 ENDIF
171
172 ELSEIF (ilev == 25) THEN
173
174 penflag = ipari(43)
175 irot = ipari(48)
176
178 . x ,v ,vr ,a ,ar ,
179 . ms ,stifn ,stifr ,weight
180 . intbuf_tab%NSV,intbuf_tab%IRTLM,intbuf_tab%CSTS,intbuf_tab%SKEW
181 . intbuf_tab%FSM,intbuf_tab%FINI,fsav ,fncont ,nsn ,
182 . intbuf_tab%SPENALTY,intbuf_tab%STFR_PENALTY,intbuf_tab%VARIABLES(14),penflag,irot ,
183 . noint ,nodnx_sms ,dmint2 , sav_for_pena ,
184 . ms_pena ,dt2t ,neltst ,ityptst ,
185 . ivisc ,h3d_data ,fncontp ,ftcontp )
186 ELSEIF (ilev == 26) THEN
187
188
189 penflag = ipari(43)
190 irot = ipari(48)
191
193 . x ,v ,vr ,a ,ar ,
194 . ms ,stifn ,stifr ,weight ,intbuf_tab%IRECTM,
195 . intbuf_tab%NSV,intbuf_tab%IRTLM,intbuf_tab%DSM,intbuf_tab%FSM,intbuf_tab%FINI,
196 . fsav ,fncont ,nsn ,intbuf_tab%SPENALTY,intbuf_tab%STFR_PENALTY,
197 . intbuf_tab%VARIABLES(14),noint ,nodnx_sms ,dmint2 ,sav_for_pena,
198 . ms_pena ,in ,dt2t ,neltst ,ityptst ,
199 . sav_iner_poff,h3d_data ,fncontp ,ftcontp )
200 ELSEIF (ilev == 27) THEN
201
202 penflag = ipari(43)
203 irot = ipari(48)
204
206 . x ,v ,vr ,a ,ar ,
207 . ms ,in ,stifn ,stifr ,weight ,
208 . intbuf_tab%IRECTM,intbuf_tab%NSV,intbuf_tab%IRTLM,intbuf_tab%DSM,intbuf_tab%FSM,
209 . intbuf_tab%FINI,intbuf_tab%SPENALTY,intbuf_tab%STFR_PENALTY,fsav ,fncont ,
210 . intbuf_tab%IRUPT,intbuf_tab%VARIABLES(14),intbuf_tab%CSTS,intbuf_tab%MSR,adi ,
211 . intbuf_tab%SMAS,intbuf_tab%SINER,intbuf_tab%NMAS(1+nmn),sav_for_pena,ms_pena ,
212 . noint ,nodnx_sms ,dmint2 ,dt2t ,neltst ,
213 . ityptst ,nsn ,nmn ,idel2 , penflag ,irot ,intbuf_tab%SKEW,
214 . sav_iner_poff, h3d_data ,intbuf_tab%DPARA,
215 . intbuf_tab%MSEGTYP2,intbuf_tab%CSTS_BIS,t2fac_sms ,fncontp ,ftcontp )
216
217 ELSEIF (ilev == 28) THEN
218
219 penflag = ipari(43)
220 irot = ipari(48)
221
223 . x ,v ,vr ,a ,ar ,
224 . ms ,in ,stifn ,stifr ,weight ,
225 . intbuf_tab%IRECTM,intbuf_tab%NSV,intbuf_tab%IRTLM,intbuf_tab%DSM,intbuf_tab%FSM,
226 . intbuf_tab%FINI,intbuf_tab%SPENALTY,intbuf_tab%STFR_PENALTY,fsav ,fncont ,
227 . intbuf_tab%IRUPT,intbuf_tab%VARIABLES(14),intbuf_tab%CSTS,intbuf_tab%MSR,adi ,
228 . intbuf_tab%SMAS,intbuf_tab%SINER,intbuf_tab%NMAS(1+nmn),sav_for_pena,ms_pena ,
229 . noint ,nodnx_sms ,dmint2 ,dt2t ,neltst ,
230 . ityptst ,nsn ,nmn ,idel2 ,penflag ,
231 . irot,intbuf_tab%SKEW,
232 . intbuf_tab%DPARA,intbuf_tab%NMAS,dmast,adm,sav_iner_poff,
233 . h3d_data ,fncontp ,ftcontp )
234
235 ELSEIF (ilev == 30 .AND. iroddl /= 0) THEN
236
238 . ms ,in ,x ,v ,vr ,
239 . a ,ar ,stifn ,stifr ,weight ,
240 . intbuf_tab%NSV,intbuf_tab%MSR,intbuf_tab%IRTLM,intbuf_tab%IRECTM
241 . intbuf_tab%NMAS,intbuf_tab%SMAS,intbuf_tab%SINER,fsav ,nsn ,
242 . nmn ,idel2 ,fncont ,h3d_data ,fncontp ,
243 . ftcontp)
244
245 ELSEIF(ilev == 0 .OR. iroddl == 0)THEN
246
248 1 nsn ,nmn ,a ,intbuf_tab%IRECTM,intbuf_tab%CSTS,
249 2 intbuf_tab%MSR,intbuf_tab%NSV,intbuf_tab%IRTLM,ms ,weight ,
250 3 stifn ,intbuf_tab%NMAS,idel2 ,intbuf_tab%SMAS,x ,
251 4 v ,fsav ,fncont ,h3d_data ,intbuf_tab%CSTS_BIS ,
252 5 fncontp ,ftcontp )
253 IF(iroddl /= 0)THEN
255 1 nsn ,nmn ,ar ,intbuf_tab%IRECTM, intbuf_tab%CSTS,
256 2 intbuf_tab%MSR,intbuf_tab%NSV,intbuf_tab%IRTLM,in ,ms ,
257 3 a ,x ,weight ,stifr ,stifn ,
258 4 idel2 ,intbuf_tab%SMAS,intbuf_tab%SINER,intbuf_tab%NMAS(1+nmn),adi,
259 5 sav_iner_poff,h3d_data,intbuf_tab%CSTS_BIS)
260 ENDIF
261
262 ENDIF
263
264 ELSE
265
266
267
268
269 IF((ilev == 1.OR.ilev == 3).AND.iroddl /= 0)THEN
270
272 1 nsn ,nmn ,a ,intbuf_tab%IRECTM,intbuf_tab%DPARA,
273 2 intbuf_tab%MSR,intbuf_tab%NSV,intbuf_tab%IRTLM,ms ,weight ,
274 3 ar ,in ,x ,stifn ,stifr ,
275 4 fskyi2 ,iadi2 ,ilev ,dmast ,adm ,
276 5 intbuf_tab%NMAS,i0 ,nir ,i2size ,idel2 ,
277 6 intbuf_tab%SMAS,intbuf_tab%SINER,v ,intbuf_tab%CSTS,fsav ,
278 7 fncont ,h3d_data ,fncontp ,ftcontp )
279 ELSEIF(ilev == 2)THEN
280
281 i0base = i0
283 1 nsn ,nmn ,a ,intbuf_tab%CSTS,intbuf_tab%MSR,
284 2 intbuf_tab%NSV,ms ,weight ,stifn ,intbuf_tab%NMAS,
285 3 fskyi2 ,iadi2 ,i0 ,nir ,i2size ,
286 4 intbuf_tab%IRECTM,x ,v ,fsav ,fncont ,
287 5 intbuf_tab%IRTLM ,h3d_data,intbuf_tab%CSTS_BIS ,fncontp ,ftcontp )
288
289 IF(iroddl /= 0)THEN
290 i0 = i0base
292 1 nsn ,nmn ,ar ,intbuf_tab%IRECTM,intbuf_tab%CSTS,
293 2 intbuf_tab%MSR,intbuf_tab%NSV,intbuf_tab%IRTLM,in ,ms ,
294 3 a ,x ,weight ,stifr ,fskyi2 ,
295 4 iadi2 ,i0 ,nir ,i2size ,stifn ,
296 5 intbuf_tab%CSTS_BIS)
297 ENDIF
298 ELSEIF(ilev == 4)THEN
299
300 i0base = i0
302 1 nsn ,nmn ,a ,intbuf_tab%CSTS,intbuf_tab%MSR,
303 2 intbuf_tab%NSV,ms ,weight ,stifn ,intbuf_tab%NMAS,
304 3 fskyi2 ,iadi2 ,i0 ,nir ,i2size ,
305 4 intbuf_tab%IRECTM,x ,v ,fsav ,fncont ,
306 5 intbuf_tab%IRTLM ,h3d_data ,intbuf_tab%CSTS_BIS ,
307 6 fncontp ,ftcontp )
308 IF(iroddl /= 0) THEN
309 i0 = i0base
310 CALL i2mzerop(fskyi2,i0 ,nir ,i2size,iadi2,
311 2 nsn ,intbuf_tab%NSV,weight)
312 ENDIF
313 ELSEIF (ilev == 10 .OR. ilev == 11 .OR. ilev == 12 .OR.
314 . ilev == 20 .OR. ilev == 21 .OR. ilev == 22) THEN
315
316 i0base = i0
318 1 a ,ms ,stifn ,weight ,iadi2 ,
319 2 fskyi2 ,intbuf_tab%NSV,intbuf_tab%IRUPT,intbuf_tab%CSTS,intbuf_tab%FSM,
320 3 nir ,nsn ,i0 ,i2size,intbuf_tab%CSTS_BIS)
321
322 IF(iroddl /= 0)THEN
323 i0 = i0base
325 1 x ,a ,ar ,ms ,in ,
326 2 stifn ,stifr ,weight ,intbuf_tab%IRECTM,intbuf_tab%NSV,
327 3 intbuf_tab%MSR,intbuf_tab%IRTLM,intbuf_tab%IRUPT,intbuf_tab%CSTS,iadi2 ,
328 4 fskyi2 ,nsn ,nmn ,nir ,i0 ,
329 5 i2size ,idel2 ,intbuf_tab%CSTS_BIS)
330 ENDIF
331
332 ELSEIF (ilev == 25) THEN
333 penflag = ipari(43)
334 irot = ipari(48)
335
336
338 . x ,v ,vr ,a ,ar ,
339 . ms ,stifn ,stifr ,weight ,intbuf_tab%IRECTM,
340 . intbuf_tab%NSV,intbuf_tab%IRTLM,intbuf_tab%CSTS,intbuf_tab%SKEW,intbuf_tab%DSM,
341 . intbuf_tab%FSM,intbuf_tab%FINI,fsav ,fncont ,nsn ,
342 . i0 ,i2size ,iadi2 ,fskyi2 ,intbuf_tab%SPENALTY,
343 . intbuf_tab%STFR_PENALTY,intbuf_tab%VARIABLES(14),penflag,irot ,noint ,
344 . nodnx_sms ,dmint2 ,dt2t ,neltst ,ityptst ,
345 . ivisc ,h3d_data ,fncontp ,ftcontp )
346
347 ELSEIF (ilev == 26) THEN
348
350 . x ,v ,vr ,a ,ar ,
351 . ms ,stifn ,stifr ,weight ,intbuf_tab%IRECTM,
352 . intbuf_tab%NSV,intbuf_tab%IRTLM,intbuf_tab%DSM,intbuf_tab%FSM,intbuf_tab%FINI,
353 . fsav ,fncont ,nsn ,i0 ,i2size ,
354 . iadi2 ,fskyi2 ,intbuf_tab%SPENALTY,intbuf_tab%STFR_PENALTY,intbuf_tab%VARIABLES(14),
355 . noint ,nodnx_sms ,dmint2 ,in ,dt2t ,neltst ,ityptst , h3d_data,
356 . fncontp ,ftcontp )
357
358 ELSEIF (ilev == 27) THEN
359
360 penflag = ipari(43)
361 irot = ipari(48)
362
364 . x ,v ,vr ,a ,ar ,
365 . ms ,in ,stifn ,stifr ,weight ,
366 . intbuf_tab%IRECTM,intbuf_tab%NSV,intbuf_tab%IRTLM,intbuf_tab%CSTS,intbuf_tab%SKEW
367 . intbuf_tab%DSM,intbuf_tab%FSM,intbuf_tab%FINI,fsav ,fncont ,
368 . nsn ,i0 ,i2size ,iadi2 ,fskyi2 ,
369 . intbuf_tab%SPENALTY,intbuf_tab%STFR_PENALTY,intbuf_tab%VARIABLES(14),penflag,irot ,
370 . noint ,nodnx_sms ,dmint2 ,dt2t ,neltst ,
371 . ityptst ,intbuf_tab%NMAS,intbuf_tab%SMAS,intbuf_tab%SINER,intbuf_tab%NMAS(1+nmn),
372 . intbuf_tab%MSR,intbuf_tab%IRUPT,adi ,idel2,nir ,
373 . nmn, h3d_data ,intbuf_tab%DPARA,
374 . intbuf_tab%MSEGTYP2 ,intbuf_tab%CSTS_BIS,t2fac_sms,fncontp ,ftcontp )
375
376 ELSEIF (ilev == 28) THEN
377
378 penflag = ipari(43)
379 irot = ipari(48)
381 . x ,v ,vr ,a ,ar ,
382 . ms ,in ,stifn ,stifr ,weight ,
383 . intbuf_tab%IRECTM,intbuf_tab%NSV,intbuf_tab%IRTLM,intbuf_tab%CSTS,intbuf_tab%SKEW,
384 . intbuf_tab%DSM,intbuf_tab%FSM,intbuf_tab%FINI,fsav ,fncont ,
385 . nsn ,i0 ,i2size ,iadi2 ,fskyi2 ,
386 . intbuf_tab%SPENALTY,intbuf_tab%STFR_PENALTY,intbuf_tab%VARIABLES(14),penflag,irot ,
387 . noint ,nodnx_sms ,dmint2 ,dt2t ,neltst ,
388 . ityptst ,intbuf_tab%NMAS,intbuf_tab%SMAS,intbuf_tab%SINER,intbuf_tab%NMAS(1+nmn),
389 . intbuf_tab%MSR,intbuf_tab%IRUPT,adi ,idel2,nir ,
390 . nmn ,intbuf_tab%DPARA,
391 . intbuf_tab%NMAS,dmast,adm,h3d_data,fncontp ,ftcontp )
392
393 ELSEIF (ilev == 30 .AND. iroddl /= 0) THEN
394
396 . ms ,in ,x ,v ,vr ,
397 . a ,ar ,stifn ,stifr ,weight ,
398 . intbuf_tab%NSV,intbuf_tab%MSR,intbuf_tab%IRTLM,intbuf_tab%IRECTM,intbuf_tab%CSTS,
399 . nsn ,nmn ,idel2 ,i0 ,i2size ,
400 . iadi2 ,intbuf_tab%NMAS,intbuf_tab%SMAS,intbuf_tab%SINER,fskyi2,
401 . fsav ,fncont ,h3d_data ,fncontp ,ftcontp)
402
403 ELSEIF(ilev == 0.OR.iroddl == 0)THEN
404
405 i0base = i0
407 1 nsn ,nmn ,a ,intbuf_tab%CSTS,intbuf_tab%NSV,
408 2 ms ,weight ,stifn ,intbuf_tab%NMAS,fskyi2 ,
409 3 iadi2 ,i0 ,nir ,i2size ,idel2 ,
410 4 intbuf_tab%SMAS,intbuf_tab%IRECTM,x ,v ,fsav ,
411 5 fncont ,intbuf_tab%IRTLM ,h3d_data ,intbuf_tab%CSTS_BIS ,
412 6 fncontp ,ftcontp )
413
414 IF(iroddl /= 0)THEN
415 i0 = i0base
417 1 nsn ,nmn ,ar ,intbuf_tab%IRECTM,intbuf_tab%CSTS,
418 2 intbuf_tab%MSR,intbuf_tab%NSV,intbuf_tab%IRTLM,in ,ms ,
419 3 a ,x ,weight ,stifr ,fskyi2 ,
420 4 stifn ,iadi2 ,i0 ,nir ,i2size ,
421 5 idel2 ,intbuf_tab%SMAS,intbuf_tab%SINER,intbuf_tab%NMAS(1+nmn),adi ,
422 6 h3d_data ,intbuf_tab%CSTS_BIS )
423 ENDIF
424
425
426 ENDIF
427 ENDIF
428
429
430
431 intth = ipari(47)
432 kthe = intbuf_tab%VARIABLES(15)
433 IF( intth > 0 ) THEN
434 i0 = i0base
436 1 x ,nsn ,intbuf_tab%NSV
437 2 weight ,intbuf_tab%IRECTM ,intbuf_tab%CSTS ,iadi2
438 3 temp ,intbuf_tab%AREAS2 ,fthe ,ftheskyi2,condn ,
439 4 condnskyi2,i0 ,itab ,idt_therm ,theaccfact)
440 ENDIF
441
442 RETURN
subroutine i2curvf(ms, in, x, v, vr, a, ar, stifn, stifr, weight, nsv, msr, irtl, irect, crst, mmass, smass, siner, fsav, nsn, nmn, idel2, fncont, h3d_data, fncontp, ftcontp)
subroutine i2curvfp(ms, in, x, v, vr, a, ar, stifn, stifr, weight, nsv, msr, irtl, irect, crst, nsn, nmn, idel2, i0, i2size, iadi2, mmass, smass, siner, fskyi2, fsav, fncont, h3d_data, fncontp, ftcontp)
subroutine i2mom10(nsn, nmn, ar, irect, crst, msr, nsv, irtl, in, ms, a, x, weight, stifr, stifn, irupt, ilev, csts_bis)
subroutine i2for10(a, ms, stifn, weight, irect, nsv, msr, irtl, irupt, crst, fsm, nsn, csts_bis)
subroutine i2mom10p(x, a, ar, ms, in, stifn, stifr, weight, irect, nsv, msr, irtl, irupt, crst, iadi2, fskyi2, nsn, nmn, nir, i0, i2size, idel2, csts_bis)
subroutine i2for10p(a, ms, stifn, weight, iadi2, fskyi2, nsv, irupt, crst, fsm, nir, nsn, i0, i2size, csts_bis)
subroutine i2for25(x, v, vr, a, ar, ms, stifn, stifr, weight, irect, nsv, irtl, crst, skew, xini, dx, fini, fsav, fncont, nsn, stfn, stfr, visc, penflag, irot, noint, nodnx_sms, dmint2, sav_for_pena, ms_pena, dt2t, neltst, ityptst, ivisc, h3d_data, fncontp, ftcontp)
subroutine i2for25p(x, v, vr, a, ar, ms, stifn, stifr, weight, irect, nsv, irtl, crst, skew, xini, dx, fini, fsav, fncont, nsn, i0, i2size, iadi2, fskyi2, stfn, stfr, visc, penflag, irot, noint, nodnx_sms, dmint2, dt2t, neltst, ityptst, ivisc, h3d_data, fncontp, ftcontp)
subroutine i2for26(x, v, vr, a, ar, ms, stifn, stifr, weight, irect, nsv, irtl, dr, dl, fini, fsav, fncont, nsn, stfn, stfr, visc, noint, nodnx_sms, dmint2, sav_for_pena, ms_pena, in, dt2t, neltst, ityptst, miner, h3d_data, fncontp, ftcontp)
subroutine i2for26p(x, v, vr, a, ar, ms, stifn, stifr, weight, irect, nsv, irtl, dr, dl, fini, fsav, fncont, nsn, i0, i2size, iadi2, fskyi2, stfn, stfr, visc, noint, nodnx_sms, dmint2, in, dt2t, neltst, ityptst, h3d_data, fncontp, ftcontp)
subroutine i2for27(x, v, vr, a, ar, ms, in, stifn, stifr, weight, irectm, nsv, irtlm, dr, dl, fini, stfn_pen, stfr_pen, fsav, fncont, ipena, visc, csts, msr, adi, smas, siner, mmas, sav_for_pena, ms_pena, noint, nodnx_sms, dmint2, dt2t, neltst, ityptst, nsn, nmn, idel2, penflag, irot, skew, sav_iner_poff, h3d_data, dpara, msegtyp2, csts_bis, t2fac_sms, fncontp, ftcontp)
subroutine i2for27p(x, v, vr, a, ar, ms, in, stifn, stifr, weight, irectm, nsv, irtlm, csts, skew, dr, dl, fini, fsav, fncont, nsn, i0, i2size, iadi2, fskyi2, stfn_pen, stfr_pen, visc, penflag, irot, noint, nodnx_sms, dmint2, dt2t, neltst, ityptst, mmas, smas, siner, miner, msr, ipena, adi, idel2, nir, nmn, h3d_data, dpara, msegtyp2, csts_bis, t2fac_sms, fncontp, ftcontp)
subroutine i2for28(x, v, vr, a, ar, ms, in, stifn, stifr, weight, irectm, nsv, irtlm, dr, dl, fini, stfn_pen, stfr_pen, fsav, fncont, ipena, visc, csts, msr, adi, smas, siner, mmas, sav_for_pena, ms_pena, noint, nodnx_sms, dmint2, dt2t, neltst, ityptst, nsn, nmn, idel2, penflag, irot, skew, dpara, nmas, dmast, adm, sav_iner_poff, h3d_data, fncontp, ftcontp)
subroutine i2for28p(x, v, vr, a, ar, ms, in, stifn, stifr, weight, irectm, nsv, irtlm, csts, skew, dr, dl, fini, fsav, fncont, nsn, i0, i2size, iadi2, fskyi2, stfn_pen, stfr_pen, visc, penflag, irot, noint, nodnx_sms, dmint2, dt2t, neltst, ityptst, mmas, smas, siner, miner, msr, ipena, adi, idel2, nir, nmn, dpara, nmas, dmast, adm, h3d_data, fncontp, ftcontp)
subroutine i2mom3(nsn, nmn, ar, irect, crst, msr, nsv, irtl, in, ms, a, x, weight, stifr, stifn, miner, csts_bis)
subroutine i2for3(nsn, nmn, a, irect, crst, msr, nsv, irtl, ms, weight, stifn, mmass, x, v, fsav, fncont, h3d_data, csts_bis, fncontp, ftcontp)
subroutine i2mom3n(nsn, nmn, ar, irect, crst, msr, nsv, irtl, in, ms, a, x, weight, stifr, stifn, idel2, smass, siner, nmas, adi, miner, h3d_data, csts_bis)
subroutine i2for3o(nsn, nmn, a, irect, crst, msr, nsv, irtl, ms, weight, stifn, mmass, x, v, fsav, fncont, csts_bis, h3d_data, fncontp, ftcontp)
subroutine i2fomo3(nsn, nmn, a, irect, dpara, msr, nsv, irtl, ms, weight, ar, in, x, stifn, stifr, ilev, dmast, adm, mmass, idel2, smass, siner, v, crst, fsav, fncont, h3d_data, fncontp, ftcontp)
subroutine i2for3n(nsn, nmn, a, irect, crst, msr, nsv, irtl, ms, weight, stifn, mmass, idel2, smass, x, v, fsav, fncont, h3d_data, csts_bis, fncontp, ftcontp)
subroutine i2mom3pn(nsn, nmn, ar, irect, crst, msr, nsv, irtl, in, ms, a, x, weight, stifr, fskyi2, stifn, iadi2, i0, nir, i2size, idel2, smass, siner, miner, adi, h3d_data, csts_bis)
subroutine i2mzerop(fskyi2, i0, nir, i2size, iadi2, nsn, nsv, weight)
subroutine i2for3p(nsn, nmn, a, crst, msr, nsv, ms, weight, stifn, mmass, fskyi2, iadi2, i0, nir, i2size, irect, x, v, fsav, fncont, irtl, h3d_data, csts_bis, fncontp, ftcontp)
subroutine i2for3pn(nsn, nmn, a, crst, nsv, ms, weight, stifn, mmass, fskyi2, iadi2, i0, nir, i2size, idel2, smass, irect, x, v, fsav, fncont, irtl, h3d_data, csts_bis, fncontp, ftcontp)
subroutine i2fomo3p(nsn, nmn, a, irect, dpara, msr, nsv, irtl, ms, weight, ar, in, x, stifn, stifr, fskyi2, iadi2, ilev, dmast, adm, mmass, i0, nir, i2size, idel2, smass, siner, v, crst, fsav, fncont, h3d_data, fncontp, ftcontp)
subroutine i2for3po(nsn, nmn, a, crst, msr, nsv, ms, weight, stifn, mmass, fskyi2, iadi2, i0, nir, i2size, irect, x, v, fsav, fncont, irtl, h3d_data, csts_bis, fncontp, ftcontp)
subroutine i2mom3p(nsn, nmn, ar, irect, crst, msr, nsv, irtl, in, ms, a, x, weight, stifr, fskyi2, iadi2, i0, nir, i2size, stifn, csts_bis)
subroutine i2therm(x, nsn, nsv, irtl, ms, weight, irect, crst, iadi2, kthe, temp, areas, fthe, ftheskyi, condn, condnskyi, i0, itab, idt_therm, theaccfact)