OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i22mainf.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!|| i22mainf ../engine/source/interfaces/int22/i22mainf.F
25!||--- called by ------------------------------------------------------
26!|| intfop2 ../engine/source/interfaces/interf/intfop2.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../engine/source/output/message/message.F
29!|| arret ../engine/source/system/arret.F
30!|| i22cor3 ../engine/source/interfaces/int22/i22cor3.F
31!|| i22for3 ../engine/source/interfaces/int22/i22for3.F
32!|| i22wetsurf ../engine/source/interfaces/int22/i22wetsurf.F
33!|| my_barrier ../engine/source/system/machine.F
34!|| startime ../engine/source/system/timer_mod.F90
35!|| stoptime ../engine/source/system/timer_mod.F90
36!|| sum_6_float_sens ../engine/source/system/parit.f
37!||--- uses -----------------------------------------------------
38!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
39!|| groupdef_mod ../common_source/modules/groupdef_mod.F
40!|| h3d_mod ../engine/share/modules/h3d_mod.F
41!|| i22tri_mod ../common_source/modules/interfaces/cut-cell-search_mod.F
42!|| intbufdef_mod ../common_source/modules/interfaces/intbufdef_mod.f90
43!|| message_mod ../engine/share/message_module/message_mod.F
44!|| output_mod ../common_source/modules/output/output_mod.f90
45!|| timer_mod ../engine/source/system/timer_mod.F90
46!||====================================================================
47 SUBROUTINE i22mainf(OUTPUT,TIMERS,
48 1 IPARI ,INTBUF_TAB ,X ,A ,
49 2 ICODT ,FSAV ,V ,MS ,DT2T ,
50 3 NELTST ,ITYPTST ,ITAB ,STIFN ,FSKYI ,
51 4 ISKY ,FCONT ,NIN ,LINDMAX ,KINET ,
52 5 JTASK ,NB_JLT ,NB_JLT_NEW ,NB_STOK_N ,
53 6 NISKYFI ,NEWFRONT ,NSTRF ,SECFCUM ,ICONTACT ,
54 7 VISCN ,NUM_IMP,
55 9 NS_IMP ,NE_IMP ,IND_IMP ,FSAVSUB ,NRTMDIM ,
56 A FSAVBAG ,IGRBRIC ,
57 B EMINX ,IXS ,IXS16 ,IXS20 ,FNCONT ,
58 C FTCONT ,IAD_ELEM ,FR_ELEM ,RCONTACT ,ACONTACT ,
59 D PCONTACT ,TEMP ,FTHE ,FTHESKYI,
60 E PM ,IPARG ,IAD17 ,MSKYI_SMS ,ISKYI_SMS,
61 F NODNX_SMS ,MS0 ,ELBUF_TAB ,NV46 ,FBSAV6 ,
62 G ISENSINT ,DIMFB ,H3D_DATA )
63C=======================================================================
64C-----------------------------------------------
65C D e s c r i p t i o n
66C-----------------------------------------------
67C Interface Type22 (/INTER/TYPE22) is an FSI coupling method based on cut cell method.
68C This experimental cut cell method is not completed, abandoned, and is not an official option.
69C
70C-----------------------------------------------
71C M o d u l e s
72C-----------------------------------------------
73 USE output_mod
74 USE timer_mod
75 USE elbufdef_mod
76 USE i22tri_mod
77 USE intbufdef_mod
78 USE h3d_mod
79 USE message_mod
80 USE groupdef_mod
81C-----------------------------------------------
82C I m p l i c i t T y p e s
83C-----------------------------------------------
84#include "implicit_f.inc"
85C-----------------------------------------------
86C G l o b a l P a r a m e t e r s
87C-----------------------------------------------
88#include "mvsiz_p.inc"
89C-----------------------------------------------
90C C o m m o n B l o c k s
91C-----------------------------------------------
92#include "com01_c.inc"
93#include "com04_c.inc"
94#include "com08_c.inc"
95#include "param_c.inc"
96#include "warn_c.inc"
97#include "task_c.inc"
98#include "parit_c.inc"
99#include "timeri_c.inc"
100C-----------------------------------------------
101C D u m m y A r g u m e n t s
102C-----------------------------------------------
103 TYPE(output_), INTENT(INOUT) :: OUTPUT
104 TYPE (ELBUF_STRUCT_), TARGET, DIMENSION(NGROUP) :: ELBUF_TAB
105 TYPE(TIMER_), INTENT(INOUT) :: TIMERS
106 INTEGER NELTST,ITYPTST,NIN,NEWFRONT,NSTRF(*),
107 . NRTMDIM, IAD17,NV46, ISENSINT(*), DIMFB
108 INTEGER IPARI(NPARI,NINTER), ICODT(*),ICONTACT(*),
109 . ITAB(*), ISKY(*), KINET(*),
110 . IPARG(NPARG,*)
111 INTEGER NB_JLT,NB_JLT_NEW,NB_STOK_N,JTASK,
112 . NISKYFI, LINDMAX
113 INTEGER NUM_IMP,NS_IMP(*),NE_IMP(*),IND_IMP(*)
114 INTEGER IXS(*) ,IXS16(*) ,IXS20(*)
115 INTEGER IAD_ELEM(2,*),FR_ELEM(*),
116 . ISKYI_SMS(*), NODNX_SMS(*)
117 my_real
118 . eminx(*)
119 my_real dt2t,
120 . x(*), a(3,*), fsav(*), v(3,*),fsavbag(*),
121 . ms(*),stifn(*),fskyi(lskyi,4),fcont(3,*),ms0(*),
122 . secfcum(7,numnod,nsect),viscn(*), fsavsub(*),
123 . fncont(3,*), ftcont(3,*), rcontact(*), acontact(*),
124 . pcontact(*),
125 . temp(*),fthe(*),ftheskyi(lskyi),pm(npropm,*),
126 . mskyi_sms(*)
127
128 DOUBLE PRECISION FBSAV6(12,6,DIMFB)
129
130 TYPE(intbuf_struct_) INTBUF_TAB
131 TYPE(H3D_DATABASE) :: H3D_DATA
132C-----------------------------------------------
133 TYPE (GROUP_) , DIMENSION(NGRBRIC) :: IGRBRIC
134C-----------------------------------------------
135C L o c a l V a r i a b l e s
136C-----------------------------------------------
137 INTEGER I, I_STOK, JLT_NEW, JLT , NFT, IVIS2,
138 . IBC, NOINT, NSEG, ISECIN, IBAG, IADM,
139 . IGAP, INACTI, IFQ, MFROT, IGSTI, NISUB,
140 . NB_LOC, I_STOK_LOC,DEBUT,
141 . ILAGM, LENR, LENT, MAXCC,INTTH,I22GRSH3N,SFSAVPARIT,
142 . IERROR,ISU1
143 INTEGER IX1(MVSIZ), IX2(MVSIZ), IX3(MVSIZ), IX4(MVSIZ),
144 . NSVG(MVSIZ), CB_LOC(MVSIZ),CE_LOC(MVSIZ),
145 . cand_b_n(mvsiz),cand_e_n(mvsiz),kini(mvsiz),
146 . index2(lindmax),
147 . isdsiz(nspmd+1),ircsiz(nspmd+1),itag(numnod),
148 . ieleci(mvsiz), nsms(mvsiz), iad, j, h
149 my_real
150 . startt, fric, gap, stopt,
151 . visc,viscf,stiglo,gapmin,
152 . kmin, kmax, gapmax,rstif,fheat,tint,rhoh
153C-----------------------------------------------
154 my_real
155 . nx1(mvsiz), nx2(mvsiz), nx3(mvsiz), nx4(mvsiz),
156 . ny1(mvsiz), ny2(mvsiz), ny3(mvsiz), ny4(mvsiz),
157 . nz1(mvsiz), nz2(mvsiz), nz3(mvsiz), nz4(mvsiz),
158 . lb1(mvsiz), lb2(mvsiz), lb3(mvsiz), lb4(mvsiz),
159 . lc1(mvsiz), lc2(mvsiz), lc3(mvsiz), lc4(mvsiz),
160 . p1(mvsiz), p2(mvsiz), p3(mvsiz), p4(mvsiz),
161 . x1(mvsiz), x2(mvsiz), x3(mvsiz), x4(mvsiz),
162 . y1(mvsiz), y2(mvsiz), y3(mvsiz), y4(mvsiz),
163 . z1(mvsiz), z2(mvsiz), z3(mvsiz), z4(mvsiz),
164 . xi(mvsiz), yi(mvsiz), zi(mvsiz), stif(mvsiz),
165 . n1(mvsiz), n2(mvsiz), n3(mvsiz), pene(mvsiz),
166 . h1(mvsiz), h2(mvsiz), h3(mvsiz), h4(mvsiz),
167 . gapv(mvsiz),vxi(mvsiz),vyi(mvsiz),vzi(mvsiz),msi(mvsiz),
168 . tempi(mvsiz),phi(mvsiz),areasi(mvsiz)
169 my_real
170 . , DIMENSION(:,:),ALLOCATABLE :: surf
171 my_real
172 . , DIMENSION(:), ALLOCATABLE :: pres
173 SAVE surf,pres
174 my_real
175 . anglt, padm
176 INTEGER NRTMFT, NRTMLT, NMNFT, NMNLT, NRADM
177 my_real
178 . NNX1(MVSIZ), NNX2(MVSIZ), NNX3(MVSIZ), NNX4(MVSIZ),
179 . NNY1(MVSIZ), NNY2(MVSIZ), NNY3(MVSIZ), NNY4(MVSIZ),
180 . NNZ1(MVSIZ), NNZ2(MVSIZ), NNZ3(MVSIZ), NNZ4(MVSIZ),
181 . N_SCUT(1:3,NBCUT_MAX,MVSIZ), !normal for cut surface (1:nbcut<8)
182 . CoG(3,NBCUT_MAX,MVSIZ),
183 . Seff(NBCUT_MAX,MVSIZ), Delta(4,NBCUT_MAX,MVSIZ)
184 INTEGER ICURV
185 my_real, DIMENSION(:,:,:), ALLOCATABLE,TARGET :: FSAVPARIT
186 my_real, POINTER :: pfsavparit
187 INTEGER :: NRTM, NSN, NTY
188
189 iadm = 0
190
191 nrtm = ipari(4,nin)
192 nsn = ipari(5,nin)
193 nty = ipari(7,nin)
194 noint = ipari(15,nin)
195 inacti = ipari(22,nin) ! unused
196 ibag = ipari(32,nin) ! Pressure Correction
197 nisub = ipari(36,nin)
198 isu1 = ipari(45,nin)
199 intth = ipari(47,nin)
200
201 !Geometrical tracking
202 i22grsh3n = ipari(48,nin)
203C
204 stiglo = -intbuf_tab%STFAC(1)
205 startt = intbuf_tab%VARIABLES(3)
206 stopt = intbuf_tab%VARIABLES(11)
207 fric = intbuf_tab%VARIABLES(1)
208 gap = intbuf_tab%VARIABLES(2)
209 gapmin = intbuf_tab%VARIABLES(13)
210 visc = intbuf_tab%VARIABLES(14)
211 viscf = intbuf_tab%VARIABLES(15)
212 gapmax = intbuf_tab%VARIABLES(16)
213 kmin = intbuf_tab%VARIABLES(17)
214 kmax = intbuf_tab%VARIABLES(18)
215 rstif = intbuf_tab%VARIABLES(20)
216 fheat = intbuf_tab%VARIABLES(21)
217 tint = intbuf_tab%VARIABLES(22)
218 i_stok = intbuf_tab%I_STOK(1)
219
220 debut = 0
221
222 !-------------------------------!
223 !UNPLUG CONDITIONS !
224 !-------------------------------!
225 IF(startt>tt) RETURN
226 IF(tt>stopt) RETURN
227
228 !-------------------------------!
229 !STATIC CUTTING !
230 !-------------------------------!
231 nb_loc = i_stok / nthread
232 IF (jtask==nthread) THEN
233 i_stok_loc = i_stok-nb_loc*(nthread-1)
234 ELSE
235 i_stok_loc = nb_loc
236 ENDIF
237 debut = (jtask-1)*nb_loc
238 i_stok = 0
239
240 !-------------------------------!
241 !recalculate du i_stok
242 !-------------------------------!
243 DO i = debut+1, debut+i_stok_loc
244 IF(intbuf_tab%CAND_N(i)/=0) THEN
245 i_stok = i_stok + 1
246 index2(i_stok) = i
247 ENDIF
248 ENDDO
249
250 !-------------------------------!
251 !/DEBUG/CAND/NN !
252 !-------------------------------!
253 IF (debug(3)>=1) THEN
254 nb_jlt = nb_jlt + i_stok_loc
255 nb_stok_n = nb_stok_n + i_stok
256 ENDIF
257
258c------------------------------------------------
259 CALL my_barrier()
260 IF(jtask==1)THEN
261 ALLOCATE(surf(3,nrtmdim))
262 ALLOCATE(pres(nrtmdim))
263 DO i = 1, nrtm
264 pres(i) = zero
265 surf(1,i) = zero
266 surf(2,i) = zero
267 surf(3,i) = zero
268 ENDDO
269
270!PARITH/ON TO BE DONE LATER
271
272! IF(IPARIT/=0)THEN
273! !C sort depending on user node ids
274! !C calculate parith/on pressure
275! !C for a given face, node order must be the same one
276! CALL I18TRI(
277! 1 INDEX2 ,I_STOK,INTBUF_TAB%CAND_N,ITAB ,NSN ,
278! 2 INTBUF_TAB%NSV,NIN )
279! END IF
280 ENDIF
281 CALL my_barrier()
282c------------------------------------------------
283C
284 sfsavparit = 0
285 NULLIFY(pfsavparit)
286 DO i=1,nisub+1
287 IF(isensint(i)/=0) THEN
288 sfsavparit = sfsavparit + 1
289 ENDIF
290 ENDDO
291 IF (sfsavparit /= 0) THEN
292 ALLOCATE(fsavparit(nisub+1,11,i_stok),stat=ierror)
293 IF(ierror/=0) THEN
294 CALL ancmsg(msgid=19,anmode=aninfo,
295 . c1='(/INTER/TYPE22)')
296 CALL arret(2)
297 ENDIF
298 DO j=1,i_stok
299 DO i=1,11
300 DO h=1,nisub+1
301 fsavparit(h,i,j) = zero
302 ENDDO
303 ENDDO
304 ENDDO
305 ELSE
306 ALLOCATE(fsavparit(0,0,0),stat=ierror)
307 IF(ierror/=0) THEN
308 CALL ancmsg(msgid=19,anmode=aninfo,
309 . c1='(/INTER/TYPE22)')
310 CALL arret(2)
311 ENDIF
312 ENDIF
313C
314 DO nft = 0 , i_stok - 1 , nvsiz
315
316 jlt = min( nvsiz, i_stok - nft )
317
318! if(ibug22_trivox==1)then
319! print *, "================================================"
320! print *, " I22MAINF : "
321! print *, " NFT = ", NFT
322! print *, " JLT = ", JLT
323! print *, "================================================"
324! endif
325
326 ! cand_b & cand_e replaced with cand_b_n & cand_e_n
327 CALL i22cor3(
328 1 jlt ,x , intbuf_tab%IRECTM ,intbuf_tab%NSV ,cand_e_n ,
329 2 cand_b_n ,intbuf_tab%STFM , intbuf_tab%STFNS ,x1 ,x2 ,
330 3 x3 ,x4 , y1 ,y2 ,y3 ,
331 4 y4 ,z1 , z2 ,z3 ,z4 ,
332 5 xi ,yi , zi ,stif ,ix1 ,
333 6 ix2 ,ix3 , ix4 ,nsvg ,igap ,
334 7 gap ,intbuf_tab%GAP_S , intbuf_tab%GAP_M ,gapv ,ms ,
335 8 vxi ,vyi ,
336 a vzi ,msi , nsn ,v ,kinet ,
337 b kini ,nty , nin ,igsti ,kmin ,
338 c kmax ,gapmax , gapmin ,iadm ,index2(nft+1) ,
339 d intth ,temp , intbuf_tab%CAND_E(1) ,intbuf_tab%CAND_N(1) ,
340 e tempi ,phi , intbuf_tab%AREAS ,intbuf_tab%IELEC ,areasi ,
341 f ieleci ,nodnx_sms , nsms ,intbuf_tab%GAP_SL ,intbuf_tab%GAP_ML,
342 g igrbric(isu1)%ENTITY,jtask)
343
344 CALL i22wetsurf(
345 1 jlt ,cand_b_n ,cand_e_n ,cb_loc ,ce_loc ,
346 2 x1 ,x2 ,x3 ,x4 ,y1 ,
347 3 y2 ,y3 ,y4 ,z1 ,z2 ,
348 4 z3 ,z4 ,xi ,yi ,zi ,
349 5 nx1 ,nx2 ,nx3 ,nx4 ,ny1 ,
350 6 ny2 ,ny3 ,ny4 ,nz1 ,nz2 ,
351 7 nz3 ,nz4 ,lb1 ,lb2 ,lb3 ,
352 8 lb4 ,lc1 ,lc2 ,lc3 ,lc4 ,
353 9 p1 ,p2 ,p3 ,p4 ,ix1 ,
354 a ix2 ,ix3 ,ix4 ,nsvg ,stif ,
355 b jlt_new ,gapv ,inacti ,intbuf_tab%CAND_P ,n_scut ,
356 c index2(nft+1) ,vxi ,vyi ,
357 d vzi ,msi ,kini ,surf ,ibag ,
358 e itab ,intbuf_tab%IRECTM ,intbuf_tab%I_STOK(1) ,ixs ,nft ,
359 f cog ,seff ,delta ,x)
360
361 !jlt = jlt_new
362
363
364 !tmp
365 jlt_new =1
366 !-------!
367
368 IF (imonm > 0 .AND. jtask == 1) CALL startime(timers,20)
369 IF(jlt_new/=0) THEN
370 ipari(29,nin) = 1
371 IF (debug(3)>=1)nb_jlt_new = nb_jlt_new + jlt_new
372 IF (sfsavparit /= 0) pfsavparit => fsavparit(1,1,nft+1)
373 CALL i22for3(output,
374 1 jlt ,a ,v ,ibc ,icodt ,
375 2 fsav ,gap ,fric ,ms ,visc ,
376 3 viscf ,noint ,intbuf_tab%STFNS ,itab ,cb_loc ,
377 4 stiglo ,stifn ,stif ,fskyi ,isky ,
378 5 nx1 ,nx2 ,nx3 ,nx4 ,ny1 ,
379 6 ny2 ,ny3 ,ny4 ,nz1 ,nz2 ,
380 7 nz3 ,nz4 ,lb1 ,lb2 ,lb3 ,
381 8 lb4 ,lc1 ,lc2 ,lc3 ,lc4 ,
382 9 p1 ,p2 ,p3 ,p4 ,fcont ,
383 b ix1 ,ix2 ,ix3 ,ix4 ,nsvg ,
384 c ivis2 ,neltst ,ityptst ,dt2t ,intth ,
385 d gapv ,inacti ,intbuf_tab%CAND_P ,index2(nft+1) ,niskyfi ,
386 e kinet ,newfront ,isecin ,nstrf ,secfcum ,
387 f x ,intbuf_tab%IRECTM ,ce_loc ,mfrot ,ifq ,
388 g intbuf_tab%FRIC_P ,intbuf_tab%FTSAVX ,intbuf_tab%FTSAVY ,intbuf_tab%FTSAVZ ,
389 + intbuf_tab%XFILTR ,
390 h intbuf_tab%IFPEN ,ibag ,icontact ,
391 j viscn ,vxi ,vyi ,vzi ,msi ,
392 k kini ,nin ,nisub ,intbuf_tab%LISUB ,intbuf_tab%ADDSUBS ,
393 l intbuf_tab%ADDSUBM ,intbuf_tab%LISUBS ,intbuf_tab%LISUBM ,fsavsub ,
394 . intbuf_tab%CAND_N ,
395 m ipari(33,nin) ,ipari(39,nin) ,pres ,fncont ,ms0 ,
396 n n_scut ,surf ,cog ,cand_e_n ,seff ,
397 o elbuf_tab ,x1 ,x2 ,x3 ,x4 ,
398 3 y1 ,y2 ,y3 ,y4 ,z1 ,
399 4 z2 ,z3 ,z4 ,ixs ,nv46 ,
400 5 delta ,isensint ,pfsavparit ,iparg ,h3d_data )
401 ENDIF
402 IF (imonm > 0 .AND. jtask == 1) CALL stoptime(timers,20)
403 ENDDO
404c
405c
406 IF (sfsavparit /= 0)THEN
407 CALL sum_6_float_sens(fsavparit, nisub+1, 11, i_stok,1,i_stok,
408 . fbsav6, 12, 6, dimfb, isensint )
409 ENDIF
410 IF (ALLOCATED(fsavparit)) DEALLOCATE (fsavparit)
411C
412 IF(inacti==7.AND.ibag/=0)THEN
413 CALL my_barrier()
414 IF(jtask==1) THEN
415! CALL ASS18(NRTM ,INTBUF_TAB%IRECTM,A ,X ,SURF ,
416! 2 PRES ,FCONT ,IBAG,FSAV,FSAVBAG,
417! 3 FR_I18,ISKY ,FSKYI )
418 END IF
419 CALL my_barrier()
420 IF(jtask == 1) DEALLOCATE(surf,pres)
421 ENDIF
422
423
424 CALL my_barrier()
425
426 IF(jtask == 1) THEN
427 IF(ALLOCATED(surf))DEALLOCATE(surf)
428 IF(ALLOCATED(pres))DEALLOCATE(pres)
429 ENDIF
430
431C
432 RETURN
433 END
434
i22cor3
subroutine i22cor3(jlt, x, irect, nsv, cand_e, cand_b, stf, stfn, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, xi, yi, zi, stif, ix1, ix2, ix3, ix4, nsvg, igap, gap, gap_s, gap_m, gapv, ms, vxi, vyi, vzi, msi, nsn, v, kinet, kini, ity, nin, igsti, kmin, kmax, gapmax, gapmin, iadm, index, intth, temp, cand__e, cand__b, tempi, phi, areas, ielec, areasi, ieleci, nodnx_sms, nsms, gap_s_l, gap_m_l, bufbric, jtask)
Definition i22cor3.F:49
i22for3
subroutine i22for3(output, jlt, a, v, ibc, icodt, fsav, gap, fric, ms, visc, viscf, noint, stfn, itab, cb_loc, stiglo, stifn, stif, fskyi, isky, nx1, nx2, nx3, nx4, ny1, ny2, ny3, ny4, nz1, nz2, nz3, nz4, lb1, lb2, lb3, lb4, lc1, lc2, lc3, lc4, p1, p2, p3, p4, fcont, ix1, ix2, ix3, ix4, nsvg, ivis2, neltst, ityptst, dt2t, intth, gapv, inacti, cand_p, index, niskyfi, kinet, newfront, isecin, nstrf, secfcum, x, irect, ce_loc, mfrot, ifq, frot_p, cand_fx, cand_fy, cand_fz, alpha0, ifpen, ibag, icontact, viscn, vxi, vyi, vzi, msi, kini, nin, nisub, lisub, addsubs, addsubm, lisubs, lisubm, fsavsub, cand_n, ilagm, icurv, pres, fncont, ms0, n_scut, n_surf, cog, cand_e, swet, elbuf_tab, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, ixs, nv46, delta, isensint, fsavparit, iparg, h3d_data)
Definition i22for3.F:65
i22mainf
subroutine i22mainf(output, timers, ipari, intbuf_tab, x, a, icodt, fsav, v, ms, dt2t, neltst, ityptst, itab, stifn, fskyi, isky, fcont, nin, lindmax, kinet, jtask, nb_jlt, nb_jlt_new, nb_stok_n, niskyfi, newfront, nstrf, secfcum, icontact, viscn, num_imp, ns_imp, ne_imp, ind_imp, fsavsub, nrtmdim, fsavbag, igrbric, eminx, ixs, ixs16, ixs20, fncont, ftcont, iad_elem, fr_elem, rcontact, acontact, pcontact, temp, fthe, ftheskyi, pm, iparg, iad17, mskyi_sms, iskyi_sms, nodnx_sms, ms0, elbuf_tab, nv46, fbsav6, isensint, dimfb, h3d_data)
Definition i22mainf.F:63
i22wetsurf
subroutine i22wetsurf(jlt, cand_b, cand_e, cb_loc, ce_loc, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, xi, yi, zi, nx1, nx2, nx3, nx4, ny1, ny2, ny3, ny4, nz1, nz2, nz3, nz4, lb1, lb2, lb3, lb4, lc1, lc2, lc3, lc4, p1, p2, p3, p4, ix1, ix2, ix3, ix4, nsvg, stif, jlt_new, gapv, inacti, cand_p, n_scut, index, vxi, vyi, vzi, msi, kini, surf, ibag, itab, irect, i_stok, ixs, nft, cog, seff, delta, x)
Definition i22wetsurf.F:52
min
#define min(a, b)
Definition macros.h:20
groupdef_mod
Definition groupdef_mod.F:668
h3d_mod
Definition h3d_mod.F:311
i22tri_mod
Definition cut-cell-search_mod.F:212
message_mod
Definition message_mod.F:1257
sum_6_float_sens
subroutine sum_6_float_sens(f, a, b, c, jft, jlt, f6, d, e, g, isensint)
Definition parit.F:540
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:895
arret
subroutine arret(nn)
Definition arret.F:86
my_barrier
subroutine my_barrier
Definition machine.F:31
startime
subroutine startime(event, itask)
Definition timer.F:93
stoptime
subroutine stoptime(event, itask)
Definition timer.F:135