OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
intti2f.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!|| intti2f ../engine/source/interfaces/interf/intti2f.F
25!||--- called by ------------------------------------------------------
26!|| intti1 ../engine/source/interfaces/interf/intti1.F
27!||--- calls -----------------------------------------------------
28!|| i2curvf ../engine/source/interfaces/interf/i2curvf.F
29!|| i2curvfp ../engine/source/interfaces/interf/i2curvfp.F
30!|| i2fomo3 ../engine/source/interfaces/interf/i2for3.F
31!|| i2fomo3p ../engine/source/interfaces/interf/i2for3p.F
32!|| i2for10 ../engine/source/interfaces/interf/i2for10.F
33!|| i2for10p ../engine/source/interfaces/interf/i2for10p.F
34!|| i2for25 ../engine/source/interfaces/interf/i2for25.F
35!|| i2for25p ../engine/source/interfaces/interf/i2for25p.F
36!|| i2for26 ../engine/source/interfaces/interf/i2for26.F
37!|| i2for26p ../engine/source/interfaces/interf/i2for26p.F
38!|| i2for27 ../engine/source/interfaces/interf/i2for27.f
39!|| i2for27p ../engine/source/interfaces/interf/i2for27p.F
40!|| i2for28 ../engine/source/interfaces/interf/i2for28.F
41!|| i2for28p ../engine/source/interfaces/interf/i2for28p.F
42!|| i2for3 ../engine/source/interfaces/interf/i2for3.F
43!|| i2for3n ../engine/source/interfaces/interf/i2for3.F
44!|| i2for3o ../engine/source/interfaces/interf/i2for3.F
45!|| i2for3p ../engine/source/interfaces/interf/i2for3p.F
46!|| i2for3pn ../engine/source/interfaces/interf/i2for3p.F
47!|| i2for3po ../engine/source/interfaces/interf/i2for3p.F
48!|| i2mom10 ../engine/source/interfaces/interf/i2for10.F
49!|| i2mom10p ../engine/source/interfaces/interf/i2for10p.F
50!|| i2mom3 ../engine/source/interfaces/interf/i2for3.F
51!|| i2mom3n ../engine/source/interfaces/interf/i2for3.F
52!|| i2mom3p ../engine/source/interfaces/interf/i2for3p.F
53!|| i2mom3pn ../engine/source/interfaces/interf/i2for3p.F
54!|| i2mzerop ../engine/source/interfaces/interf/i2for3p.F
55!|| i2therm ../engine/source/interfaces/interf/i2therm.F
56!||--- uses -----------------------------------------------------
57!|| h3d_mod ../engine/share/modules/h3d_mod.F
58!|| intbufdef_mod ../common_source/modules/interfaces/intbufdef_mod.F90
59!||====================================================================
60 SUBROUTINE intti2f(
61 1 IPARI ,X ,V ,A ,
62 2 VR ,AR ,MS ,IN ,WEIGHT ,STIFN ,
63 3 STIFR ,FSKYI2 ,IADI2 ,I2MSCH ,DMAST ,ADM ,
64 4 I0 ,NIR ,I2SIZE ,ADI ,IGEO ,BUFGEO ,
65 5 FSAV ,FNCONT ,NODNX_SMS,DMINT2 ,SAV_FOR_PENA ,
66 6 MS_PENA,DT2T ,NELTST ,ITYPTST,INTBUF_TAB,TEMP,
67 7 FTHE ,FTHESKYI2,CONDN,CONDNSKYI2,ITAB,
68 8 SAV_INER_POFF ,H3D_DATA,T2FAC_SMS,FNCONTP ,
69 A FTCONTP,IDT_THERM ,THEACCFACT)
70C-----------------------------------------------
71C M o d u l e s
72C-----------------------------------------------
73 USE intbufdef_mod
74 USE h3d_mod
75C-----------------------------------------------
76C I m p l i c i t T y p e s
77C-----------------------------------------------
78#include "implicit_f.inc"
79C-----------------------------------------------
80C D u m m y A r g u m e n t s
81C-----------------------------------------------
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
87 my_real
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
98C-----------------------------------------------
99C C o m m o n B l o c k s
100C-----------------------------------------------
101#include "com01_c.inc"
102#include "parit_c.inc"
103#include "param_c.inc"
104C-----------------------------------------------
105C L o c a l V a r i a b l e s
106C-----------------------------------------------
107 integer
108 . nsn,nmn,nrts,nrtm,ilev,i0base,idel2,penflag,irot,noint,intth,ivisc
109 my_real
110 . kthe
111C=======================================================================
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)
120C
121 i0base = i0
122C
123 IF (iparit == 0) THEN
124 IF((ilev == 1.OR.ilev == 3).AND.iroddl /= 0)THEN
125C FORMULATION SANS MOMENT SUR MAIN
126 CALL i2fomo3(
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
135C FORMULATION AVEC MOMENT SUR MAIN
136 CALL i2for3(
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 )
141C
142 IF(iroddl /= 0)THEN
143 CALL i2mom3(
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
150C FORMULATION decouplant les ddl rot
151 CALL i2for3o(
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
158C FORMULATION AVEC RUPTURE
159 CALL i2for10(
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)
163C
164 IF(iroddl /= 0) THEN
165 CALL i2mom10(
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
171C
172 ELSEIF (ilev == 25) THEN
173c formulation en penalite avec rotations
174 penflag = ipari(43)
175 irot = ipari(48)
176C
177 CALL i2for25(
178 . x ,v ,vr ,a ,ar ,
179 . ms ,stifn ,stifr ,weight ,intbuf_tab%IRECTM,
180 . intbuf_tab%NSV,intbuf_tab%IRTLM,intbuf_tab%CSTS,intbuf_tab%SKEW,intbuf_tab%DSM,
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
187c formulation en penalite + rot
188
189 penflag = ipari(43)
190 irot = ipari(48)
191C
192 CALL i2for26(
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
201c formulation en penalite avec rotations
202 penflag = ipari(43)
203 irot = ipari(48)
204c
205 CALL i2for27(
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 )
216C
217 ELSEIF (ilev == 28) THEN
218c formulation en penalite avec rotations
219 penflag = ipari(43)
220 irot = ipari(48)
221c
222 CALL i2for28(
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 )
234C
235 ELSEIF (ilev == 30 .AND. iroddl /= 0) THEN
236c formulation avec la courbure du main
237 CALL i2curvf(
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,intbuf_tab%CSTS,
241 . intbuf_tab%NMAS,intbuf_tab%SMAS,intbuf_tab%SINER,fsav ,nsn ,
242 . nmn ,idel2 ,fncont ,h3d_data ,fncontp ,
243 . ftcontp)
244C
245 ELSEIF(ilev == 0 .OR. iroddl == 0)THEN
246C FORMULATION AVEC MOMENT SUR MAIN
247 CALL i2for3n(
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
254 CALL i2mom3n(
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
261C
262 ENDIF
263C----------- Fin Parith OFF
264 ELSE
265C
266C parith/on spmd
267C attention aux routines duales du parith/on spmd
268C
269 IF((ilev == 1.OR.ilev == 3).AND.iroddl /= 0)THEN
270C FORMULATION SANS MOMENT SUR MAIN
271 CALL i2fomo3p(
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
280C FORMULATION AVEC MOMENT SUR MAIN
281 i0base = i0
282 CALL i2for3p(
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 )
288C
289 IF(iroddl /= 0)THEN
290 i0 = i0base
291 CALL i2mom3p(
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
299C FORMULATION decouplant les ddl rot
300 i0base = i0
301 CALL i2for3po(
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
315C FORMULATION AVEC RUPTURE
316 i0base = i0
317 CALL i2for10p(
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)
321C
322 IF(iroddl /= 0)THEN
323 i0 = i0base
324 CALL i2mom10p(
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
331C
332 ELSEIF (ilev == 25) THEN
333 penflag = ipari(43)
334 irot = ipari(48)
335c formulation en penalite
336C
337 CALL i2for25p(
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 )
346C
347 ELSEIF (ilev == 26) THEN
348c formulation en penalite avec rotations
349 CALL i2for26p(
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 )
357C
358 ELSEIF (ilev == 27) THEN
359c formulation en cinematique - penalite
360 penflag = ipari(43)
361 irot = ipari(48)
362
363 CALL i2for27p(
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 )
375c
376 ELSEIF (ilev == 28) THEN
377c formulation en cinematique - penalite
378 penflag = ipari(43)
379 irot = ipari(48)
380 CALL i2for28p(
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 )
392c
393 ELSEIF (ilev == 30 .AND. iroddl /= 0) THEN
394c formulation avec la courbure du main
395 CALL i2curvfp(
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)
402C
403 ELSEIF(ilev == 0.OR.iroddl == 0)THEN
404C FORMULATION AVEC MOMENT SUR MAIN
405 i0base = i0
406 CALL i2for3pn(
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 )
413C
414 IF(iroddl /= 0)THEN
415 i0 = i0base
416 CALL i2mom3pn(
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
424C
425C---
426 ENDIF
427 ENDIF
428
429C---------Thermal modelling for interface Type 2--------
430C
431 intth = ipari(47)
432 kthe = intbuf_tab%VARIABLES(15)
433 IF( intth > 0 ) THEN
434 i0 = i0base
435 CALL i2therm (
436 1 x ,nsn ,intbuf_tab%NSV ,intbuf_tab%IRTLM ,ms ,
437 2 weight ,intbuf_tab%IRECTM ,intbuf_tab%CSTS ,iadi2 ,kthe ,
438 3 temp ,intbuf_tab%AREAS2 ,fthe ,ftheskyi2,condn ,
439 4 condnskyi2,i0 ,itab ,idt_therm ,theaccfact)
440 ENDIF
441C
442 RETURN
443 END
i2curvf
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)
Definition i2curvf.F:42
i2curvfp
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)
Definition i2curvfp.F:42
i2mom10
subroutine i2mom10(nsn, nmn, ar, irect, crst, msr, nsv, irtl, in, ms, a, x, weight, stifr, stifn, irupt, ilev, csts_bis)
Definition i2for10.F:147
i2for10
subroutine i2for10(a, ms, stifn, weight, irect, nsv, msr, irtl, irupt, crst, fsm, nsn, csts_bis)
Definition i2for10.F:31
i2mom10p
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)
Definition i2for10p.F:221
i2for10p
subroutine i2for10p(a, ms, stifn, weight, iadi2, fskyi2, nsv, irupt, crst, fsm, nir, nsn, i0, i2size, csts_bis)
Definition i2for10p.F:31
i2for25
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)
Definition i2for25.F:44
i2for25p
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)
Definition i2for25p.F:45
i2for26
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)
Definition i2for26.F:43
i2for26p
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)
Definition i2for26p.F:42
i2for27
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)
Definition i2for27.F:46
i2for27p
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)
Definition i2for27p.F:46
i2for28
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)
Definition i2for28.F:45
i2for28p
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)
Definition i2for28p.F:46
i2mom3
subroutine i2mom3(nsn, nmn, ar, irect, crst, msr, nsv, irtl, in, ms, a, x, weight, stifr, stifn, miner, csts_bis)
Definition i2for3.F:1064
i2for3
subroutine i2for3(nsn, nmn, a, irect, crst, msr, nsv, irtl, ms, weight, stifn, mmass, x, v, fsav, fncont, h3d_data, csts_bis, fncontp, ftcontp)
Definition i2for3.F:311
i2mom3n
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)
Definition i2for3.F:712
i2for3o
subroutine i2for3o(nsn, nmn, a, irect, crst, msr, nsv, irtl, ms, weight, stifn, mmass, x, v, fsav, fncont, csts_bis, h3d_data, fncontp, ftcontp)
Definition i2for3.F:516
i2fomo3
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)
Definition i2for3.F:1291
i2for3n
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)
Definition i2for3.F:38
i2mom3pn
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)
Definition i2for3p.F:726
i2mzerop
subroutine i2mzerop(fskyi2, i0, nir, i2size, iadi2, nsn, nsv, weight)
Definition i2for3p.F:1577
i2for3p
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)
Definition i2for3p.F:295
i2for3pn
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)
Definition i2for3p.F:40
i2fomo3p
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)
Definition i2for3p.F:1236
i2for3po
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)
Definition i2for3p.F:538
i2mom3p
subroutine i2mom3p(nsn, nmn, ar, irect, crst, msr, nsv, irtl, in, ms, a, x, weight, stifr, fskyi2, iadi2, i0, nir, i2size, stifn, csts_bis)
Definition i2for3p.F:1046
i2therm
subroutine i2therm(x, nsn, nsv, irtl, ms, weight, irect, crst, iadi2, kthe, temp, areas, fthe, ftheskyi, condn, condnskyi, i0, itab, idt_therm, theaccfact)
Definition i2therm.F:32
intti2f
subroutine intti2f(ipari, x, v, a, vr, ar, ms, in, weight, stifn, stifr, fskyi2, iadi2, i2msch, dmast, adm, i0, nir, i2size, adi, igeo, bufgeo, fsav, fncont, nodnx_sms, dmint2, sav_for_pena, ms_pena, dt2t, neltst, ityptst, intbuf_tab, temp, fthe, ftheskyi2, condn, condnskyi2, itab, sav_iner_poff, h3d_data, t2fac_sms, fncontp, ftcontp, idt_therm, theaccfact)
Definition intti2f.F:70
h3d_mod
Definition h3d_mod.F:308