OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i11mainf.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!|| i11mainf ../engine/source/interfaces/int11/i11mainf.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!|| frictionparts_model_isot ../engine/source/interfaces/int07/frictionparts_model.F
31!|| i11ass3 ../engine/source/interfaces/int11/i11ass3.F
32!|| i11cdcor3 ../engine/source/interfaces/int11/i11cdcor3.F
33!|| i11cor3 ../engine/source/interfaces/int11/i11cor3.F
34!|| i11dst3 ../engine/source/interfaces/int11/i11dst3.F
35!|| i11for3 ../engine/source/interfaces/int11/i11for3.F
36!|| i11therm ../engine/source/interfaces/int11/i11therm.F
37!|| startime ../engine/source/system/timer_mod.F90
38!|| stoptime ../engine/source/system/timer_mod.F90
39!|| sum_6_float_sens ../engine/source/system/parit.F
40!||--- uses -----------------------------------------------------
41!|| h3d_mod ../engine/share/modules/h3d_mod.F
42!|| intbuf_fric_mod ../common_source/modules/interfaces/intbuf_fric_mod.F90
43!|| intbufdef_mod ../common_source/modules/interfaces/intbufdef_mod.F90
44!|| message_mod ../engine/share/message_module/message_mod.F
45!|| my_alloc_mod ../common_source/tools/memory/my_alloc.F90
46!|| output_mod ../common_source/modules/output/output_mod.F90
47!|| timer_mod ../engine/source/system/timer_mod.F90
48!||====================================================================
49 SUBROUTINE i11mainf(OUTPUT, TIMERS,
50 1 IPARI ,X ,A ,
51 2 ICODT ,FSAV ,V ,MS ,DT2T ,
52 3 NELTST ,ITYPTST ,ITAB ,STIFN ,FSKYI ,
53 4 ISKY ,FCONT ,NIN ,LINDMAX ,JTASK ,
54 5 NB_JLT ,NB_JLT_NEW,NB_STOK_N,NISKYFI,NEWFRONT,
55 6 NSTRF ,SECFCUM ,VISCN ,NUM_IMP, NS_IMP ,
56 7 NE_IMP ,MSKYI_SMS ,ISKYI_SMS,NODNX_SMS,ICONTACT,
57 8 INTBUF_TAB,PM ,TEMP , FTHE ,FTHESKYI ,
58 9 NPC ,TF ,CONDN ,CONDNSKYI,FBSAV6,
59 A ISENSINT,DIMFB ,FSAVSUB ,H3D_DATA ,INTBUF_FRIC_TAB,
60 B ITASK ,TAGNCONT ,KLOADPINTER,LOADPINTER,LOADP_HYD_INTER,
61 C DGAPLOADINT,S_LOADPINTER,NODADT_THERM)
62C-----------------------------------------------
63C M o d u l e s
64C-----------------------------------------------
65 USE timer_mod
66 USE output_mod
67 USE intbufdef_mod
68 USE h3d_mod
69 USE intbuf_fric_mod
70 USE message_mod
71 USE my_alloc_mod
72C-----------------------------------------------
73C I m p l i c i t T y p e s
74C-----------------------------------------------
75#include "implicit_f.inc"
76C-----------------------------------------------
77C G l o b a l P a r a m e t e r s
78C-----------------------------------------------
79#include "mvsiz_p.inc"
80C-----------------------------------------------
81C C o m m o n B l o c k s
82C-----------------------------------------------
83#include "com04_c.inc"
84#include "com08_c.inc"
85#include "param_c.inc"
86#include "warn_c.inc"
87#include "task_c.inc"
88#include "parit_c.inc"
89#include "impl1_c.inc"
90#include "timeri_c.inc"
91#include "macro.inc"
92C-----------------------------------------------
93C D u m m y A r g u m e n t s
94C-----------------------------------------------
95 TYPE(output_), intent(inout) :: output
96 TYPE(TIMER_), INTENT(INOUT) :: TIMERS
97 INTEGER ,INTENT(IN) :: NODADT_THERM
98 INTEGER NELTST,ITYPTST,NIN,NEWFRONT,NSTRF(*)
99 INTEGER IPARI(*), ICODT(*),
100 . ITAB(*), ISKY(*),ICONTACT(*),TAGNCONT(NLOADP_HYD_INTER,NUMNOD)
101 INTEGER , INTENT(IN) :: S_LOADPINTER
102 INTEGER , INTENT(IN) :: KLOADPINTER(NINTER+1),LOADPINTER(S_LOADPINTER),
103 . LOADP_HYD_INTER(NLOADP_HYD)
104C Interface Statistics
105 INTEGER NB_JLT,NB_JLT_NEW,NB_STOK_N,JTASK,
106 . NISKYFI, LINDMAX,ITASK
107 INTEGER NUM_IMP,NS_IMP(*),NE_IMP(*),
108 . iskyi_sms(*), nodnx_sms(*),npc(*), isensint(*),dimfb
109C REAL
110 my_real , INTENT(IN) :: dgaploadint(s_loadpinter)
111 my_real dt2t,
112 . x(*), a(3,*), fsav(*), v(3,*),
113 . ms(*), stifn(*), fskyi(lskyi,nfskyi), fcont(3,*),
114 . secfcum(7,numnod,nsect), viscn(*), mskyi_sms(*),
115 . temp(*),fthe(*),ftheskyi(*),tf(*),condn(*),condnskyi(*),
116 . pm(npropm,*),fsavsub(*)
117
118 DOUBLE PRECISION FBSAV6(12,6,DIMFB)
119
120 TYPE(intbuf_struct_) INTBUF_TAB
121 TYPE(H3D_DATABASE) :: H3D_DATA
122 TYPE(intbuf_fric_struct_), TARGET, DIMENSION(NINTERFRIC) :: INTBUF_FRIC_TAB
123C-----------------------------------------------
124C L o c a l V a r i a b l e s
125C-----------------------------------------------
126 INTEGER
127 . I, J, H, IBC, NOINT, ISECIN, I_STOK,
128 . JLT , NFT, JLT_NEW, IGAP, IVIS2,
129 . NB_LOC, I_STOK_LOC,DEBUT, IGSTI,IFORM,INTTH,IKTHE,
130 . iformth,sfsavparit,nisub,intfric,nsetprts,npartfric,mfrot,
131 . ierror,iorthfric,ifric,jj,ninloadp
132 INTEGER N1(MVSIZ), N2(MVSIZ), M1(MVSIZ), M2(MVSIZ),
133 . cs_loc(mvsiz), cm_loc(mvsiz),
134 . nsms(mvsiz)
135 INTEGER
136 . IELECI(MVSIZ),IELESI(MVSIZ), IPARTFRICSI(MVSIZ),
137 . IPARTFRICMI(MVSIZ)
138 INTEGER,DIMENSION(:), ALLOCATABLE :: INDEX2
139C REAL
140 my_real
141 . STARTT, FRIC, GAP, STOPT,
142 . VISC,VISCF,STIGLO, GAPMIN, KMIN, KMAX,DTMINI,
143 . TINT,XTHE,KTHE,FRAD,DRAD,DGAPLOAD
144C-----------------------------------------------
145C REAL
146 my_real
147 . nx(mvsiz),ny(mvsiz),nz(mvsiz),
148 . hs1(mvsiz), hs2(mvsiz), hm1(mvsiz), hm2(mvsiz),
149 . stif(mvsiz),gapv(mvsiz),
150 . xxs1(mvsiz), xxs2(mvsiz), xys1(mvsiz), xys2(mvsiz),
151 . xzs1(mvsiz), xzs2(mvsiz), xxm1(mvsiz), xxm2(mvsiz),
152 . xym1(mvsiz), xym2(mvsiz), xzm1(mvsiz), xzm2(mvsiz),
153 . vxs1(mvsiz), vxs2(mvsiz), vys1(mvsiz), vys2(mvsiz),
154 . vzs1(mvsiz), vzs2(mvsiz), vxm1(mvsiz), vxm2(mvsiz),
155 . vym1(mvsiz), vym2(mvsiz), vzm1(mvsiz), vzm2(mvsiz),
156 . ms1(mvsiz), ms2(mvsiz), mm1(mvsiz), mm2(mvsiz),
157 . fni(mvsiz), tempi1(mvsiz),tempi2(mvsiz),tempm1(mvsiz),
158 . tempm2(mvsiz),phis1(mvsiz),phis2(mvsiz),phim1(mvsiz),
159 . phim2(mvsiz),areac(mvsiz), condints1(mvsiz),
160 . condints2(mvsiz),condintm1(mvsiz),condintm2(mvsiz),
161 . penrad(mvsiz),fx1(mvsiz), fx2(mvsiz),
162 . fx3(mvsiz), fx4(mvsiz),fy1(mvsiz), fy2(mvsiz),
163 . fy3(mvsiz), fy4(mvsiz),fz1(mvsiz), fz2(mvsiz),
164 . fz3(mvsiz), fz4(mvsiz),k1(mvsiz) , k2(mvsiz) ,
165 . k3(mvsiz) , k4(mvsiz) ,c1(mvsiz) , c2(mvsiz) ,
166 . c3(mvsiz) , c4(mvsiz)
167 my_real, DIMENSION(:,:,:), ALLOCATABLE :: fsavparit
168 my_real
169 . xfiltr_fric,fric_coefs(mvsiz,10),viscffric(mvsiz),fricc(mvsiz)
170 INTEGER, DIMENSION(:) ,POINTER :: TABCOUPLEPARTS_FRIC
171 INTEGER, DIMENSION(:) ,POINTER :: TABPARTS_FRIC
172 INTEGER, DIMENSION(:) ,POINTER :: ADPARTS_FRIC
173 my_real, DIMENSION(:) ,POINTER :: TABCOEF_FRIC
174
175 INTEGER,TARGET, DIMENSION(1):: TABCOUPLEPARTS_FRIC_BID
176 INTEGER,TARGET, DIMENSION(1):: TABPARTS_FRIC_BID
177 INTEGER,TARGET, DIMENSION(1):: ADPARTS_FRIC_BID
178 my_real,TARGET, DIMENSION(1):: TABCOEF_FRIC_BID
179 INTEGER :: NRTS, NTY
180C-----------------------------------------------
181 CALL my_alloc(index2,lindmax)
182C-----------------------------------------------
183C calculation of interface buffer addresses
184C addresses of buffers J10-JFI and K10-KFI
185C are systematically replaced by arrays
186C JD(i) et KD(i), en gardant les memes numeros d'indexes.
187C old direct addresses Jn, Kn are modified
188C in the routine without usual additional comments
189C----------------------------------------------------
190C
191 nrts =ipari(3)
192 nty =ipari(7)
193 ibc =ipari(11)
194 ivis2 =ipari(14)
195 noint =ipari(15)
196 igap =ipari(21)
197 isecin=ipari(28)
198 iform =ipari(30)
199 igsti =ipari(34)
200 nisub =ipari(36)
201C
202 stiglo=-intbuf_tab%STFAC(1)
203 startt=intbuf_tab%VARIABLES(3)
204 stopt =intbuf_tab%VARIABLES(11)
205 IF(startt>tt) RETURN
206 IF(tt>stopt) RETURN
207C
208 fric =intbuf_tab%VARIABLES(1)
209 gap =intbuf_tab%VARIABLES(2)
210 gapmin=intbuf_tab%VARIABLES(13)
211 visc =intbuf_tab%VARIABLES(14)
212 viscf =intbuf_tab%VARIABLES(15)
213 kmin =intbuf_tab%VARIABLES(17)
214 kmax =intbuf_tab%VARIABLES(18)
215 dtmini=intbuf_tab%VARIABLES(41)
216 dgapload=intbuf_tab%VARIABLES(46)
217 num_imp = 0
218 IF (impl_s==1) THEN
219 visc =zero
220 viscf =zero
221 ENDIF
222C heat interface
223 intth = ipari(47)
224 ikthe = ipari(43)
225 iformth =ipari(44)
226 kthe = intbuf_tab%VARIABLES(20)
227 xthe = intbuf_tab%VARIABLES(22)
228 tint = intbuf_tab%VARIABLES(21)
229 frad = intbuf_tab%VARIABLES(23)
230 drad = intbuf_tab%VARIABLES(24)
231 ifric = 0
232C--- Corresponding Friction model
233 intfric=ipari(72)
234 mfrot = 0
235 iorthfric = 0
236 npartfric = 0
237 xfiltr_fric = 0
238 nsetprts = 0
239 IF(intfric /= 0) THEN
240 tabcoupleparts_fric => intbuf_fric_tab(intfric)%TABCOUPLEPARTS_FRIC
241 tabcoef_fric => intbuf_fric_tab(intfric)%TABCOEF_FRIC
242 tabparts_fric => intbuf_fric_tab(intfric)%TABPARTS_FRIC
243 adparts_fric => intbuf_fric_tab(intfric)%ADPARTS_FRIC
244 xfiltr_fric = intbuf_fric_tab(intfric)%XFILTR_FRIC
245 nsetprts = intbuf_fric_tab(intfric)%NSETPRTS
246 npartfric = intbuf_fric_tab(intfric)%S_TABPARTS_FRIC
247 iorthfric = intbuf_fric_tab(intfric)%IORTHFRIC
248 ELSE
249 tabcoupleparts_fric => tabcoupleparts_fric_bid
250 tabparts_fric => tabparts_fric_bid
251 tabcoef_fric => tabcoef_fric_bid
252 adparts_fric => adparts_fric_bid
253 ENDIF
254C
255 ninloadp = ipari(95) ! load pressure related to inter
256C
257 i_stok = intbuf_tab%I_STOK(1)
258C this part is done in parallel after the calculation of element forces
259C static decoupage
260 nb_loc = i_stok / nthread
261 IF (jtask==nthread) THEN
262 i_stok_loc = i_stok-nb_loc*(nthread-1)
263 ELSE
264 i_stok_loc = nb_loc
265 ENDIF
266 debut = (jtask-1)*nb_loc
267 i_stok = 0
268C recalculation of istok
269 DO i = debut+1, debut+i_stok_loc
270 IF(intbuf_tab%CAND_N(i)<0) THEN
271 i_stok = i_stok + 1
272 index2(i_stok) = i
273C inbuf == cand_S
274 intbuf_tab%CAND_N(i) = -intbuf_tab%CAND_N(i)
275 ENDIF
276 ENDDO
277 IF (debug(3)>=1) THEN
278 nb_jlt = nb_jlt + i_stok_loc
279 nb_stok_n = nb_stok_n + i_stok
280 ENDIF
281C
282 sfsavparit = 0
283 DO i=1,nisub+1
284 IF(isensint(i)/=0) THEN
285 sfsavparit = sfsavparit + 1
286 ENDIF
287 ENDDO
288 IF (sfsavparit /= 0) THEN
289 ALLOCATE(fsavparit(nisub+1,11,i_stok),stat=ierror)
290 IF(ierror/=0) THEN
291 CALL ancmsg(msgid=19,anmode=aninfo,
292 . c1='(/INTER/TYPE11)')
293 CALL arret(2)
294 ENDIF
295 DO j=1,i_stok
296 DO i=1,11
297 DO h=1,nisub+1
298 fsavparit(h,i,j) = zero
299 ENDDO
300 ENDDO
301 ENDDO
302 ELSE
303 ALLOCATE(fsavparit(0,0,0),stat=ierror)
304 IF(ierror/=0) THEN
305 CALL ancmsg(msgid=19,anmode=aninfo,
306 . c1='(/INTER/TYPE11)')
307 CALL arret(2)
308 ENDIF
309 ENDIF
310C
311 DO nft = 0 , i_stok - 1 , nvsiz
312 jlt = min( nvsiz, i_stok - nft )
313C preparation CANDIDATES retenus
314 CALL i11cdcor3(
315 1 jlt,index2(nft+1),intbuf_tab%CAND_E,intbuf_tab%CAND_N,cm_loc,
316 2 cs_loc)
317
318 CALL i11cor3(
319 1 jlt ,intbuf_tab%IRECTS,intbuf_tab%IRECTM,x ,v ,
320 2 cs_loc ,cm_loc ,intbuf_tab%STFS ,intbuf_tab%STFM,gapmin ,
321 3 intbuf_tab%GAP_S,intbuf_tab%GAP_M,igap ,gapv ,ms ,
322 4 stif ,xxs1 ,xxs2 ,xys1 ,xys2 ,
323 5 xzs1 ,xzs2 ,xxm1 ,xxm2 ,xym1 ,
324 6 xym2 ,xzm1 ,xzm2 ,vxs1 ,vxs2 ,
325 7 vys1 ,vys2 ,vzs1 ,vzs2 ,vxm1 ,
326 8 vxm2 ,vym1 ,vym2 ,vzm1 ,vzm2 ,
327 9 ms1 ,ms2 ,mm1 ,mm2 ,n1 ,
328 a n2 ,m1 ,m2 ,nrts ,nin ,
329 b igsti ,kmin ,kmax ,nodnx_sms ,nsms ,
330 c intbuf_tab%GAP_SL, intbuf_tab%GAP_ML,intth,temp ,tempi1 ,
331 d tempi2 ,tempm1 ,tempm2,intbuf_tab%AREAS,intbuf_tab%AREAM,
332 e areac ,ieleci ,ielesi,intbuf_tab%IELEC,intbuf_tab%IELES,
333 f iformth , itab ,intfric ,intbuf_tab%IPARTFRICS,ipartfricsi,
334 g intbuf_tab%IPARTFRICM,ipartfricmi)
335
336 CALL i11dst3(
337 1 jlt ,cs_loc ,cm_loc ,hs1 ,hs2 ,
338 2 hm1 , hm2 ,nx ,ny ,nz ,
339 3 stif ,n1 ,n2 ,m1 ,m2 ,
340 4 jlt_new ,xxs1 ,xxs2 ,xys1 ,xys2 ,
341 5 xzs1 ,xzs2 ,xxm1 ,xxm2 ,xym1 ,
342 6 xym2 ,xzm1 ,xzm2 ,vxs1 ,vxs2 ,
343 7 vys1 ,vys2 ,vzs1 ,vzs2 ,vxm1 ,
344 8 vxm2 ,vym1 ,vym2 ,vzm1 ,vzm2 ,
345 9 ms1 ,ms2 ,mm1 ,mm2 ,gapv ,
346 a nsms ,index2(nft+1),drad , intfric ,ipartfricsi,
347 b ipartfricmi,dgapload)
348
349 jlt = jlt_new
350 IF(jlt_new/=0) THEN
351 IF (imonm > 0 .AND. jtask == 1) CALL startime(timers,20)
352 ipari(29) = 1
353 IF (debug(3)>=1)
354 . nb_jlt_new = nb_jlt_new + jlt_new
355
356C-------------------------------------------------------------------------------
357C Friction model : computation of friction coefficients based on Material of connected Parts
358C-------------------------------------------------------------------------------
359 IF(itask==1) CALL startime(timers,macro_timer_fric)
360 jj = 0
361 CALL frictionparts_model_isot(
362 1 intfric ,jlt ,ipartfricsi ,ipartfricmi ,adparts_fric,
363 2 nsetprts ,tabcoupleparts_fric,npartfric ,tabparts_fric,tabcoef_fric,
364 3 fric ,viscf ,intbuf_tab%FRIC_P,fric_coefs , fricc ,
365 4 viscffric ,nty ,mfrot ,iorthfric , ifric ,
366 5 jj , tint ,tempi1 ,npc ,tf ,
367 6 temp , hs1 ,hs2 ,hm1 ,hm2 ,
368 7 n1 , n2 ,m1 ,m2 ,iform )
369 IF(itask==1) CALL stoptime(timers,macro_timer_fric)
370
371 CALL i11for3(output,
372 1 jlt ,fsav ,gap ,fric ,ms ,
373 2 visc ,viscf ,noint ,itab ,cs_loc ,
374 3 cm_loc ,stif ,dt2t ,hs1 ,hs2 ,
375 4 hm1 ,hm2 ,n1 ,n2 , m1 ,
376 5 m2 ,ivis2 ,neltst ,ityptst ,nx ,
377 6 ny ,nz ,gapv,intbuf_tab%PENIS,intbuf_tab%PENIM ,
378 7 ipari(22) ,newfront,nrts ,ms1 ,ms2 ,
379 8 mm1 ,mm2 ,vxs1 ,vys1 ,vzs1 ,
380 9 vxs2 ,vys2 ,vzs2 ,vxm1 ,vym1 ,
381 a vzm1 ,vxm2 ,vym2 ,vzm2 ,nin ,
382 b dtmini,iform ,intbuf_tab%FTSAVX,intbuf_tab%FTSAVY,intbuf_tab%FTSAVZ,
383 c index2(nft+1),intbuf_tab%IFPEN ,intbuf_tab%STFS,fni ,
384 e fx1 ,fy1 ,fz1 ,fx2 ,fy2 ,
385 f fz2 ,fx3 ,fy3 ,fz3 ,fx4 ,
386 g fy4 ,fz4 ,k1 ,k2 ,k3 ,
387 h k4 ,c1 ,c2 ,c3 ,c4 ,
388 i intth ,drad ,penrad ,isensint ,fsavparit ,
389 j nisub ,nft ,intbuf_tab%ADDSUBS ,intbuf_tab%ADDSUBM,
390 k intbuf_tab%LISUBS,intbuf_tab%LISUBM,intbuf_tab%LISUB,fsavsub,fricc ,
391 l viscffric ,tagncont ,kloadpinter,loadpinter,loadp_hyd_inter ,
392 m intbuf_tab%TYPSUB,intbuf_tab%INFLG_SUBS,intbuf_tab%INFLG_SUBM ,
393 . ninloadp ,dgaploadint,
394 n s_loadpinter )
395
396C
397 IF( intth > 0 ) THEN
398 CALL i11therm(
399 1 jlt ,pm ,intth ,penrad , kthe ,
400 2 tempi1 ,tempi2 ,tempm1 ,tempm2 ,phis1 ,
401 3 phis2 ,tint ,areac ,ieleci ,ielesi ,
402 4 frad ,gapv ,fni ,ikthe ,xthe ,
403 5 npc ,drad ,tf ,hs1 ,hs2 ,
404 6 hm1 ,hm2 ,condints1 ,condints2,phim1 ,
405 7 phim2 ,condintm1,condintm2 ,iformth )
406 ENDIF
407C
408 CALL i11ass3 (
409 1 jlt ,a ,nin ,noint ,cs_loc ,
410 2 stifn ,stif ,fskyi ,isky ,fcont ,
411 3 hs1 ,hs2 ,hm1 ,hm2 ,n1 ,
412 4 n2 ,m1 ,m2 ,niskyfi ,isecin ,
413 5 nstrf ,secfcum ,viscn ,nrts ,iskyi_sms,
414 6 nsms ,icontact ,mskyi_sms ,fx1 ,fy1 ,
415 7 fz1 ,fx2 ,fy2 ,fz2 ,fx3 ,
416 8 fy3 ,fz3 ,fx4 ,fy4 ,fz4 ,
417 9 k1 ,k2 ,k3 ,k4 ,c1 ,
418 a c2 ,c3 ,c4 ,intth ,phis1 ,
419 b phis2 ,phim1 ,phim2 ,fthe ,ftheskyi ,
420 c condints1 ,condints2 ,condintm1 ,condintm2 ,condn ,
421 d condnskyi ,jtask ,h3d_data ,nodadt_therm)
422C
423 IF(impl_s==1) THEN
424 DO i = 1 ,jlt_new
425 ns_imp(i+num_imp)=cs_loc(i)
426 ne_imp(i+num_imp)=cm_loc(i)
427 ENDDO
428 num_imp=num_imp+jlt_new
429 ENDIF
430 IF (imonm > 0 .AND. jtask == 1) CALL stoptime(timers,20)
431 ENDIF
432 ENDDO
433 IF (sfsavparit /= 0)THEN
434 CALL sum_6_float_sens(fsavparit, nisub+1, 11, i_stok,1,i_stok,
435 . fbsav6, 12, 6, dimfb, isensint )
436 ENDIF
437 IF(ALLOCATED(fsavparit)) DEALLOCATE (fsavparit)
438 IF(ALLOCATED(index2)) DEALLOCATE (index2)
439
440C
441 RETURN
442 END
my_real
#define my_real
Definition cppsort.cpp:32
frictionparts_model_isot
subroutine frictionparts_model_isot(intfric, jlt, ipartfricsi, ipartfricmi, adparts_fric, nset, tabcoupleparts_fric, npartfric, tabparts_fric, tabcoef_fric, fric, viscf, frot_p, fric_coefs, fricc, viscffric, nty, mfrot, iorthfric, ifric, jlt_tied, tint, tempi, npc, tf, temp, h1, h2, h3, h4, ix1, ix2, ix3, ix4, iform)
Definition frictionparts_model.F:350
i11ass3
subroutine i11ass3(jlt, a, nin, noint, cs_loc, stifn, stif, fskyi, isky, fcont, hs1, hs2, hm1, hm2, n1, n2, m1, m2, niskyfi, isecin, nstrf, secfcum, viscn, nrts, iskyi_sms, nsms, icontact, mskyi_sms, fx1, fy1, fz1, fx2, fy2, fz2, fx3, fy3, fz3, fx4, fy4, fz4, k1, k2, k3, k4, c1, c2, c3, c4, intth, phis1, phis2, phim1, phim2, fthe, ftheskyi, condints1, condints2, condintm1, condintm2, condn, condnskyi, jtask, h3d_data, nodadt_therm)
Definition i11ass3.F:51
i11cdcor3
subroutine i11cdcor3(jlt, index, cand_m, cand_s, cand_m_n, cand_s_n)
Definition i11cdcor3.F:32
i11cor3
subroutine i11cor3(jlt, irects, irectm, x, v, cand_s, cand_m, stfs, stfm, gap, gap_s, gap_m, igap, gapv, ms, stif, xxs1, xxs2, xys1, xys2, xzs1, xzs2, xxm1, xxm2, xym1, xym2, xzm1, xzm2, vxs1, vxs2, vys1, vys2, vzs1, vzs2, vxm1, vxm2, vym1, vym2, vzm1, vzm2, ms1, ms2, mm1, mm2, n1, n2, m1, m2, nrts, nin, igsti, kmin, kmax, nodnx_sms, nsms, gap_s_l, gap_m_l, intth, temp, tempi1, tempi2, tempm1, tempm2, areas, aream, areac, ieleci, ielesi, ielec, ieles, iform, itab, intfric, ipartfrics, ipartfricsi, ipartfricm, ipartfricmi)
Definition i11cor3.F:50
i11for3
subroutine i11for3(output, jlt, fsav, gap, fric, ms, visc, viscf, noint, itab, cs_loc, cm_loc, stif, dt2t, hs1, hs2, hm1, hm2, n1, n2, m1, m2, ivis2, neltst, ityptst, nx, ny, nz, gapv, penis, penim, inacti, newfront, nrts, ms1, ms2, mm1, mm2, vxs1, vys1, vzs1, vxs2, vys2, vzs2, vxm1, vym1, vzm1, vxm2, vym2, vzm2, nin, dtmini, iform, cand_fx, cand_fy, cand_fz, index, ifpen, stfs, fni, fx1, fy1, fz1, fx2, fy2, fz2, fx3, fy3, fz3, fx4, fy4, fz4, k1, k2, k3, k4, c1, c2, c3, c4, intth, drad, penrad, isensint, fsavparit, nisub, nft, addsubs, addsubm, lisubs, lisubm, lisub, fsavsub, fricc, viscffric, tagncont, kloadpinter, loadpinter, loadp_hyd_inter, typsub, inflg_subs, inflg_subm, ninloadp, dgaploadint, s_loadpinter)
Definition i11for3.F:56
i11mainf
subroutine i11mainf(output, timers, ipari, x, a, icodt, fsav, v, ms, dt2t, neltst, ityptst, itab, stifn, fskyi, isky, fcont, nin, lindmax, jtask, nb_jlt, nb_jlt_new, nb_stok_n, niskyfi, newfront, nstrf, secfcum, viscn, num_imp, ns_imp, ne_imp, mskyi_sms, iskyi_sms, nodnx_sms, icontact, intbuf_tab, pm, temp, fthe, ftheskyi, npc, tf, condn, condnskyi, fbsav6, isensint, dimfb, fsavsub, h3d_data, intbuf_fric_tab, itask, tagncont, kloadpinter, loadpinter, loadp_hyd_inter, dgaploadint, s_loadpinter, nodadt_therm)
Definition i11mainf.F:62
i11therm
subroutine i11therm(jlt, pm, intth, penrad, kthe, tempi1, tempi2, tempm1, tempm2, phis1, phis2, tint, areac, ieleci, ielesi, frad, gapv, fni, ifunctk, xthe, npc, drad, tf, hs1, hs2, hm1, hm2, condints1, condints2, phim1, phim2, condintm1, condintm2, iform)
Definition i11therm.F:39
min
#define min(a, b)
Definition macros.h:20
h3d_mod
Definition h3d_mod.F:311
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
i11dst3
subroutine i11dst3(jlt, gap, cand_s, cand_m, irects, irectm, nx, ny, nz, n1, n2, m1, m2, jlt_new, x, igap, gap_s, gap_m, gapv2, gap_s_l, gap_m_l, drad, dgapload)
Definition i11dst3.F:33
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
startime
subroutine startime(event, itask)
Definition timer.F:93
stoptime
subroutine stoptime(event, itask)
Definition timer.F:135