OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
aconve.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!|| aconve ../engine/source/ale/aconve.f
25!||--- called by ------------------------------------------------------
26!|| alethe ../engine/source/ale/alethe.F
27!||--- calls -----------------------------------------------------
28!|| a22conv3 ../engine/source/ale/alefvm/cut_cells/a22conv3.F
29!|| a4conv3 ../engine/source/ale/ale3d/a4conv3.F
30!|| aconv2 ../engine/source/ale/ale2d/aconv2.F
31!|| aconv3 ../engine/source/ale/ale3d/aconv3.F
32!|| bconv2 ../engine/source/ale/ale2d/bconv2.f
33!|| initbuf ../engine/share/resol/initbuf.F
34!|| my_barrier ../engine/source/system/machine.F
35!|| spmd_e1vois ../engine/source/mpi/fluid/spmd_cfd.F
36!|| varcondec ../engine/source/interfaces/interf/varcondec.F
37!||--- uses -----------------------------------------------------
38!|| ale_connectivity_mod ../common_source/modules/ale/ale_connectivity_mod.F
39!|| ale_mod ../common_source/modules/ale/ale_mod.f
40!|| alefvm_mod ../common_source/modules/ale/alefvm_mod.F
41!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
42!|| i22bufbric_mod ../common_source/modules/interfaces/cut-cell-search_mod.F
43!|| i22tri_mod ../common_source/modules/interfaces/cut-cell-search_mod.f
44!|| initbuf_mod ../engine/share/resol/initbuf.F
45!|| multimat_param_mod ../common_source/modules/multimat_param_mod.F90
46!|| segvar_mod ../engine/share/modules/segvar_mod.F
47!||====================================================================
48 SUBROUTINE aconve(
49 1 IPARG ,ELBUF_TAB ,FLUX ,FLU1 ,PHI ,
50 2 ALE_CONNECT ,NVAR ,ITASK ,NERCVOIS ,NESDVOIS,
51 3 LERCVOIS ,LESDVOIS ,LENCOM ,SEGVAR ,BHOLE ,
52 4 ITRIMAT ,QMV ,IFLG ,IXS ,IXQ ,
53 5 PM ,X )
54C-----------------------------------------------
55C M o d u l e s
56C-----------------------------------------------
57 USE initbuf_mod
58 USE elbufdef_mod
59 USE i22tri_mod
60 USE i22bufbric_mod
61 USE segvar_mod
62 USE ale_connectivity_mod
63 USE alefvm_mod , only:alefvm_param
64 USE ale_mod
65 USE multimat_param_mod , ONLY : m51_nvphas, m51_n0phas
66C-----------------------------------------------
67C I m p l i c i t T y p e s
68C-----------------------------------------------
69#include "implicit_f.inc"
70C-----------------------------------------------
71C C o m m o n B l o c k s
72C-----------------------------------------------
73#include "com01_c.inc"
74#include "com04_c.inc"
75#include "vect01_c.inc"
76#include "param_c.inc"
77#include "task_c.inc"
78#include "spmd_c.inc"
79#include "inter22.inc"
80#include "warn_c.inc"
81C-----------------------------------------------
82C D u m m y A r g u m e n t s
83C-----------------------------------------------
84 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
85
86 INTEGER IPARG(NPARG,NGROUP),NVAR,ITRIMAT,
87 . NERCVOIS(*),NESDVOIS(*),LERCVOIS(*), LESDVOIS(*),
88 . BHOLE(*),IXS(NIXS,NUMELS),IXQ(NIXQ,NUMELQ), LENCOM,
89 . iflg, itask
90
91 my_real flux(*), flu1(*) , phi(*) ,
92 . qmv(*) , pm(npropm,nummat), x(3, numnod)
93
94 TYPE(t_segvar) :: SEGVAR
95 TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECT
96C-----------------------------------------------
97C L o c a l V a r i a b l e s
98C-----------------------------------------------
99 INTEGER NMN, NM, NG, JMUL, IADR, I, J, K, NF1,ISILENT,NFX, IOFF, IMAT
100 INTEGER JCODV(ALE%GLOBAL%LCONV),CODTOT,NGSEG,ISEG,ISOLNOD
101 INTEGER ADD0, ADD
102
103 TYPE(L_BUFEL_) ,POINTER :: LBUF
104 TYPE(g_bufel_) ,POINTER :: GBUF
105 TYPE(buf_mat_) ,POINTER :: MBUF
106
107 my_real, DIMENSION(:), POINTER :: var, prho , pvol , peint, piad22
108 INTEGER :: ICELLv,IB,IBv,NIN,NUM, MCELL, IDX, NDIM
109C-----------------------------------------------
110C S o u r c e L i n e s
111C-----------------------------------------------
112C--------------------
113 CALL my_barrier
114C--------------------
115
116 NULLIFY (var)
117
118C=======================================================================
119C PHI IS DEFINED DEPENDING ON NVAR VALUE
120C=======================================================================
121
122 nmn=max(1,nmult)
123
124 DO nm=1,nmn
125 DO ng=itask+1,ngroup,nthread
126C ALE ON / OFF
127 IF (iparg(76, ng) == 1) cycle ! --> OFF
128 CALL initbuf(iparg ,ng ,
129 2 mtn ,llt ,nft ,iadr ,ity ,
130 3 npt ,jale ,ismstr ,jeul ,jtur ,
131 4 jthe ,jlag ,jmul ,jhbe ,jivf ,
132 5 nvaux ,jpor ,jcvt ,jclose ,jplasol ,
133 6 irep ,iint ,igtyp ,israt ,isrot ,
134 7 icsen ,isorth ,isorthg ,ifailure,jsms )
135 isilent = iparg(64,ng)
136 !-----------------------------!
137 ! UNPLUG CONDITIONS !
138 !-----------------------------!
139 IF (jale+jeul == 0) cycle
140 IF (iparg(8,ng) == 1) cycle
141 IF (max(1,jmul) < nm) cycle
142 IF (itrimat /= 0 .AND. mtn /= 51) cycle
143 !-----------------------------!
144 ! POINTERS FOR CURRENT GROUP !
145 !-----------------------------!
146 gbuf => elbuf_tab(ng)%GBUF
147 lbuf => elbuf_tab(ng)%BUFLY(nm)%LBUF(1,1,1)
148 mbuf => elbuf_tab(ng)%BUFLY(nm)%MAT(1,1,1)
149
150 CALL varcondec(jcodv,iparg(34,ng),codtot)
151 IF (jcodv(nvar) /= 0) THEN
152 isolnod = iparg(28,ng)
153 IF (jmul /= 0) mtn =iparg(24+nm,ng)
154 lft=1
155 !----------------------------!
156 ! N V A R = 1 !
157 ! in : RHO !
158 ! out : RHO * VOL !
159 !----------------------------!
160 IF (nvar == 1) THEN
161 IF(itrimat==0)THEN
162 prho => lbuf%RHO(1:llt)
163 pvol => lbuf%VOL(1:llt)
164 ELSE
165 !USE PHASIS DATA
166 add0 = m51_n0phas + (itrimat-1)*m51_nvphas
167 add = add0 + 9 ! ADD+9 => RHO
168 k = llt*(add-1) ! UVAR(I,ADD) = UVAR(K+I)
169 prho => mbuf%VAR(k+1:k+llt)
170 add = add0 + 11 ! ADD+11 => VOL
171 k = llt*(add-1) ! UVAR(I,ADD) = UVAR(K+I)
172 pvol => mbuf%VAR(k+1:k+llt)
173 !dbKEY="R/V TRIMAT "
174 END IF
175#include "vectorize.inc"
176 DO i=lft,llt
177 j=i+nft
178 phi(j)=prho(i)
179 ENDDO
180 DO i=lft,llt
181 prho(i) = prho(i)*pvol(i)
182 ENDDO
183 !----------------------------!
184 ! N V A R = 2 !
185 ! in : E (J/m3)!
186 ! out : EINT (J) !
187 !----------------------------!
188 ELSEIF (nvar == 2) THEN
189 IF(itrimat == 0)THEN
190 peint=> lbuf%EINT(1:llt)
191 pvol => lbuf%VOL(1:llt)
192 ELSE
193 !USE PHASIS DATA
194 add0 = m51_n0phas + (itrimat-1)*m51_nvphas
195 add = add0 + 8 ! ADD+8 => EINT
196 k = llt*(add-1) ! UVAR(I,ADD) = UVAR(K+I)
197 peint=> mbuf%VAR(k+1:k+llt)
198 add = add0 + 11 ! ADD+11 => VOL
199 k = llt*(add-1) ! UVAR(I,ADD) = UVAR(K+I)
200 pvol => mbuf%VAR(k+1:k+llt)
201 END IF
202#include "vectorize.inc"
203 DO i=lft,llt
204 j=i+nft
205 phi(j)=peint(i)
206 ENDDO
207 DO i=lft,llt
208 peint(i) = peint(i)*pvol(i)
209 ENDDO
210 !----------------------------!
211 ! N V A R = 3 !
212 !----------------------------!
213 ELSEIF (nvar == 3) THEN ! JTUR > 0
214#include "vectorize.inc"
215 DO i=lft,llt
216 j=i+nft
217 phi(j)=lbuf%RK(i)
218 ENDDO
219 DO i=lft,llt
220 lbuf%RK(i) = lbuf%RK(i)*lbuf%VOL(i)
221 ENDDO
222 !----------------------------!
223 ! N V A R = 4 !
224 !----------------------------!
225 ELSEIF (nvar == 4) THEN ! JTUR > 0
226#include "vectorize.inc"
227 DO i=lft,llt
228 j=i+nft
229 phi(j)=lbuf%RE(i)
230 ENDDO
231 DO i=lft,llt
232 lbuf%RE(i) = lbuf%RE(i)*lbuf%VOL(i)
233 ENDDO
234 !----------------------------!
235 ! N V A R = 5 !
236 !----------------------------!
237 !UVAR(I,1) : rho_liq*V_liq/V
238 !UVAR(I,2) : density of gas
239 !UVAR(I,3) : density of liquid
240 !UVAR(I,4) : volumetric fraction of liquid
241 !UVAR(I,5) : volumetric fraction of gas
242 ELSEIF (nvar == 5) THEN ! Law 37 or law 41
243 var => elbuf_tab(ng)%BUFLY(1)%MAT(1,1,1)%VAR(1:llt) ! <---- UVAR(1:LLT,1)
244 IF (mtn == 41) THEN
245 prho => gbuf%RHO
246#include "vectorize.inc"
247 DO i=lft,llt
248 j=i+nft
249 phi(j) = var(i) * prho(i)
250 ENDDO
251 DO i=lft,llt
252 var(i) = var(i) * prho(i) * lbuf%VOL(i)
253 ENDDO
254 ELSE
255#include "vectorize.inc"
256 DO i=lft,llt
257 j=i+nft
258 phi(j) = var(i)
259 ENDDO
260 DO i=lft,llt
261 var(i) = var(i)*lbuf%VOL(i)
262 ENDDO
263 ENDIF
264 !-------------------------------!
265 ! N V A R = 6 !
266 ! in : RHO*v_x !
267 ! out : RHO*v_x *VOL !
268 !-------------------------------!
269 ELSEIF (nvar == 6) THEN
270 IF(itrimat == 0)THEN
271 idx = 1
272 ndim = 3
273 var => gbuf%MOM(1:llt*ndim)
274#include "vectorize.inc"
275 DO i=lft,llt
276 j = i+nft
277 k = llt*(idx-1) + i
278 phi(j) = var(k)
279 var(k) = var(k) *lbuf%VOL(i) ![rhoU]*VOL, U needed for fluxes calculation
280 ENDDO
281 ENDIF
282
283 !-----------------------------!
284 ! N V A R = 7 !
285 ! in : RHO*v_y !
286 ! out : RHO*v_y *VOL !
287 !-----------------------------!
288 ELSEIF (nvar == 7) THEN
289 IF(itrimat == 0)THEN
290 idx = 2
291 ndim = 3
292 var => gbuf%MOM(1:llt*ndim)
293#include "vectorize.inc"
294 DO i=lft,llt
295 j = i+nft
296 k = llt*(idx-1) + i
297 phi(j) = var(k)
298 var(k) = var(k) *lbuf%VOL(i) ![rhoU]*VOL, U needed for fluxes calculation
299 ENDDO
300 ENDIF
301
302 !-----------------------------!
303 ! N V A R = 8 !
304 ! in : RHO*v_z !
305 ! out : RHO*v_z *VOL !
306 !-----------------------------!
307 ELSEIF (nvar == 8) THEN
308 IF(itrimat == 0)THEN
309 idx = 3
310 ndim = 3
311 var => gbuf%MOM(1:llt*ndim)
312#include "vectorize.inc"
313 DO i=lft,llt
314 j = i+nft
315 k = llt*(idx-1) + i
316 phi(j) = var(k)
317 var(k) = var(k) *lbuf%VOL(i) ![rhoU]*VOL, U needed for fluxes calculation
318 ENDDO
319 ENDIF
320
321 ELSEIF (nvar == 9 .AND. isilent == 1) THEN
322 ELSEIF (nvar == 10 .AND. isilent == 1) THEN
323 ELSE
324 !--------------------------!
325 ! DEFAULT 3D CASE !
326 !--------------------------!
327 IF (n2d == 0) THEN
328#include "vectorize.inc"
329 DO i=lft,llt
330 j=i+nft
331 imat=ixs(1,j)
332 phi(j)=pm(180+nvar,imat)*lbuf%RHO(i)
333 END DO
334 !--------------------------!
335 ! DEFAULT 2D CASE !
336 !--------------------------!
337 ELSE
338#include "vectorize.inc"
339 DO i=lft,llt
340 j=i+nft
341 imat=ixq(1,j)
342 phi(j)=pm(180+nvar,imat)*lbuf%RHO(i)
343 END DO
344 END IF
345 END IF
346 ELSE
347 DO i=lft,llt
348 j=i+nft
349 phi(j)=zero
350 ENDDO
351 ENDIF
352
353 !--------------------------!
354 ! INTERFACE 22 !
355 !--------------------------!
356 !Fill Secnd cells linked to current main cell with main cell value.
357 IF(int22 > 0)THEN
358 nin = 1
359 piad22 => elbuf_tab(ng)%GBUF%TAG22(lft:llt)
360 DO i=lft,llt
361 j = i+nft
362 ib = nint(piad22(i))
363 IF(ib==0)cycle
364 num = brick_list(nin,ib)%SecndList%Num
365 mcell = brick_list(nin,ib)%mainID
366 brick_list(nin,ib)%POLY(mcell)%PHI = phi(j)
367 DO k=1,num
368 ibv = brick_list(nin,ib)%SecndList%IBV(k)
369 icellv = brick_list(nin,ib)%SecndList%ICELLv(k)
370 brick_list(nin,ibv)%POLY(icellv)%PHI = phi(j)
371 ENDDO
372 enddo!next I
373 ENDIF
374
375 ENDDO !next NG
376
377
378 !--------------------------!
379 ! EBCS **** !
380 !--------------------------!
381 ioff = 0
382 IF(nsegflu > 0)THEN
383 ioff = numels+numelq+numeltg
384 IF(nspmd > 1) THEN
385 ioff = ioff + nsvois
386 ENDIF
387 ngseg=nsegflu/nvsiz
388 IF(nsegflu-ngseg*nvsiz > 0)ngseg=ngseg+1
389 DO i=itask+1,ngseg,nthread
390 iseg=(i-1)*nvsiz
391 SELECT CASE(nvar)
392 !RHO
393 CASE(1)
394 IF(itrimat==0)THEN
395 DO j=iseg+1,min(iseg+nvsiz,nsegflu)
396 phi(ioff+j)=segvar%RHO(j)
397 ENDDO
398 ELSE
399 DO j=iseg+1,min(iseg+nvsiz,nsegflu)
400 phi(ioff+j)=segvar%PHASE_RHO(itrimat,j)
401 ENDDO
402 ENDIF
403 !EINT
404 CASE(2)
405 IF(itrimat==0)THEN
406 DO j=iseg+1,min(iseg+nvsiz,nsegflu)
407 phi(ioff+j)=segvar%EINT(j)
408 ENDDO
409 ELSE
410 DO j=iseg+1,min(iseg+nvsiz,nsegflu)
411 phi(ioff+j)=segvar%PHASE_EINT(itrimat,j)
412 ENDDO
413 ENDIF
414 !RK
415 CASE(3)
416 DO j=iseg+1,min(iseg+nvsiz,nsegflu)
417 phi(ioff+j)=segvar%RK(j)
418 ENDDO
419 !RE
420 CASE(4)
421 DO j=iseg+1,min(iseg+nvsiz,nsegflu)
422 phi(ioff+j)=segvar%RE(j)
423 ENDDO
424 !UVAR
425 CASE(5)
426 DO j=iseg+1,min(iseg+nvsiz,nsegflu)
427 phi(ioff+j)=segvar%UVAR(j)
428 ENDDO
429 END SELECT
430 ENDDO
431 ENDIF
432C----------------------
433 CALL my_barrier
434C----------------------
435 !-----------------------------------------------!
436 ! SPMD EXCHANGES !
437 !-----------------------------------------------!
438 IF (nspmd > 1) THEN
439!$OMP SINGLE
440 CALL spmd_e1vois(phi,nercvois,nesdvois,lercvois, lesdvois,lencom )
441!$OMP END SINGLE
442 END IF
443
444
445C=======================================================================
446C CONVECTION DEPENDING ON NVAR VALUE
447C=======================================================================
448
449 ! CONVECTION OFF WITH DEBUG OPTION /DEBUG/CONV0/[cycle]
450 IF(debug(10) /= 0)THEN
451 IF(ncycle >= debug(10))THEN
452 !write (*,*) "ADVECTION SET OFF"
453 cycle
454 ENDIF
455 ENDIF
456
457 !----------------------------!
458 ! CUT CELL (INTER22) !
459 !----------------------------!
460 !TREAT FIRST POLYHEDRA FROM CUT CELL BUFFER
461 !ELEMENTARY BUFFER IS ALSO UPDATED AND USED LATER
462 !TO TREAT USUAL BRICK ELEMS WHILE DEALING WITH ADJACENT VALUES
463 IF(int22 > 0)THEN
464 nf1=nft+1+(nm-1)*numels
465 nfx=nft+(nm-1)*numels
466 CALL a22conv3(
467 1 phi ,
468 2 iflg ,
469 3 itrimat , nvar , itask ,
470 4 elbuf_tab, ixs , iparg)
471 ENDIF
472
473 DO ng=itask+1,ngroup,nthread
474C ALE ON / OFF
475 IF (iparg(76, ng) == 1) cycle ! --> OFF
476 CALL varcondec(jcodv,iparg(34,ng),codtot)
477 IF (jcodv(nvar) == 0) cycle
478 CALL initbuf( iparg ,ng ,
479 2 mtn ,llt ,nft ,iadr ,ity ,
480 3 npt ,jale ,ismstr ,jeul ,jtur ,
481 4 jthe ,jlag ,jmul ,jhbe ,jivf ,
482 5 nvaux ,jpor ,jcvt ,jclose ,jplasol ,
483 6 irep ,iint ,igtyp ,israt ,isrot ,
484 7 icsen ,isorth ,isorthg ,ifailure,jsms )
485 isilent = iparg(64,ng)
486 IF (isilent == 1) cycle !silent boundaries law 11 & 51
487 IF (iparg(8,ng) == 1) cycle
488 IF (max(1,jmul) < nm) cycle
489 IF (itrimat /= 0 .AND. mtn /= 51) cycle
490
491 isolnod = iparg(28,ng)
492
493 !-----------------------------!
494 ! pointers from current group !
495 !-----------------------------!
496 gbuf => elbuf_tab(ng)%GBUF
497 lbuf => elbuf_tab(ng)%BUFLY(nm)%LBUF(1,1,1)
498 mbuf => elbuf_tab(ng)%BUFLY(nm)%MAT(1,1,1)
499 !-----------------------------!
500 IF (jmul /= 0) THEN
501 mtn =iparg(24+nm,ng)
502 ENDIF
503 lft=1
504
505 !----------------------------!
506 ! N V A R = 1 !
507 !----------------------------!
508 IF (nvar == 1) THEN
509 IF(itrimat == 0)THEN
510 prho => lbuf%RHO(1:llt)
511 ELSE
512 !USE PHASIS DATA
513 add0 = m51_n0phas + (itrimat-1)*m51_nvphas
514 add = add0 + 9 ! ADD+9 => RHO
515 k = llt*(add-1) ! UVAR(I,ADD) = UVAR(K+I)
516 prho => mbuf%VAR(k+1:k+llt)
517 END IF
518 var => prho
519 !----------------------------!
520 ! N V A R = 2 !
521 !----------------------------!
522 ELSEIF (nvar == 2) THEN
523 IF(itrimat == 0)THEN
524 peint=> lbuf%EINT(1:llt)
525 ELSE
526 !USE PHASIS DATA
527 add0 = m51_n0phas + (itrimat-1)*m51_nvphas
528 add = add0 + 8 ! ADD+9 => RHO
529 k = llt*(add-1) ! UVAR(I,ADD) = UVAR(K+I)
530 peint => mbuf%VAR(k+1:k+llt)
531 END IF
532 var => peint
533 !----------------------------!
534 ! N V A R = 3 !
535 !----------------------------!
536 ELSEIF (nvar == 3) THEN
537 var => elbuf_tab(ng)%BUFLY(nm)%LBUF(1,1,1)%RK(1:llt)
538 !----------------------------!
539 ! N V A R = 4 !
540 !----------------------------!
541 ELSEIF (nvar == 4) THEN
542 var => elbuf_tab(ng)%BUFLY(nm)%LBUF(1,1,1)%RE(1:llt)
543 !----------------------------!
544 ! N V A R = 5 !
545 !----------------------------!
546 ELSEIF (nvar == 5) THEN
547 var => elbuf_tab(ng)%BUFLY(1)%MAT(1,1,1)%VAR(1:llt)
548 !----------------------------!
549 ! N V A R = 6 !
550 !----------------------------!
551 ELSEIF (nvar == 6) THEN
552 IF(alefvm_param%IEnabled==0)THEN
553 IF (mtn == 51 .AND. itrimat /= 0) THEN
554 var => elbuf_tab(ng)%BUFLY(1)%MAT(1,1,1)%VAR(5*llt+1:6*llt)
555 ELSE
556 var => elbuf_tab(ng)%BUFLY(1)%MAT(1,1,1)%VAR(llt+1:2*llt)
557 ENDIF
558 ELSE
559 var => elbuf_tab(ng)%GBUF%MOM( 1 : llt )
560 ENDIF
561 !----------------------------!
562 ! N V A R = 7 !
563 !----------------------------!
564 ELSEIF (nvar == 7) THEN
565 IF(alefvm_param%IEnabled == 0)THEN
566 var => elbuf_tab(ng)%BUFLY(1)%MAT(1,1,1)%VAR(2*llt+1:3*llt)
567 ELSE
568 var => elbuf_tab(ng)%GBUF%MOM( llt*1+1 : llt*1+llt )
569 ENDIF
570 !----------------------------!
571 ! N V A R = 8 !
572 !----------------------------!
573 ELSEIF (nvar == 8) THEN
574 IF(alefvm_param%IEnabled == 0)THEN
575 var => elbuf_tab(ng)%BUFLY(1)%MAT(1,1,1)%VAR(3*llt+1:4*llt)
576 ELSE
577 var => elbuf_tab(ng)%GBUF%MOM( llt*2+1 : llt*2+llt )
578 ENDIF
579 !----------------------------!
580 ! N V A R = 9 !
581 !----------------------------!
582 ELSEIF (nvar == 9) THEN
583 var => elbuf_tab(ng)%BUFLY(1)%MAT(1,1,1)%VAR(4*llt+1:5*llt)
584 ENDIF
585
586 !----------------------------!
587 ! CALL ACONVE3 !
588 !----------------------------!
589 piad22 => elbuf_tab(ng)%GBUF%TAG22(lft:llt)
590 IF (n2d == 0) THEN
591 nf1=nft+1+(nm-1)*numels
592 nfx=nft+(nm-1)*numels
593 pvol => lbuf%VOL(1:llt)
594 IF (isolnod /= 4) THEN
595 CALL aconv3(
596 1 var , phi ,flux(6*nfx+1), flu1(nf1) ,ixs ,
597 2 ale_connect , ioff ,qmv(12*nfx+1), iflg ,
598 3 piad22 , nvar ,itask)
599 ELSE
600 CALL a4conv3(
601 1 var ,phi,flux(6*nfx+1),flu1(nf1),
602 2 ale_connect ,ioff )
603 ENDIF
604 !----------------------------!
605 ! CALL ACONVE2 !
606 !----------------------------!
607 ELSE
608 nf1=nft+1+(nm-1)*numelq
609 nfx=nft+(nm-1)*numelq
610 IF (nmult == 0) THEN
611 CALL aconv2(var ,phi ,flux(4*nfx+1),flu1(nf1),
612 . ale_connect ,qmv(8*nfx+1),iflg ,ixq ,
613 . x ,ioff )
614 ELSE
615 CALL bconv2(var, phi, flux(4*nfx+1), flu1(nf1), ale_connect ,bhole ,nm)
616 ENDIF
617 ENDIF !(N2D == 0)
618 !----------------------------!
619 ENDDO ! next NG
620
621
622C----------------------
623 CALL my_barrier
624C----------------------
625 END DO !next NM
626
627C-----------------------------------------------
628 RETURN
629 END
a22conv3
subroutine a22conv3(phi, iflg, itrimat, nvar, itask, elbuf_tab, ixs, iparg)
Definition a22conv3.F:40
a4conv3
subroutine a4conv3(vtot, phi, flux, flu1, ale_connect, ioff)
Definition a4conv3.F:32
aconv2
subroutine aconv2(vtot, phi, flux, flu1, ale_connect, qmv, iflg, ixq, x, ioff)
Definition aconv2.F:31
aconv3
subroutine aconv3(vtot, phi, flux, flu1, ixs, ale_connect, ioff, qmv, iflg, tag22, nvar, itask)
Definition aconv3.F:37
aconve
subroutine aconve(iparg, elbuf_tab, flux, flu1, phi, ale_connect, nvar, itask, nercvois, nesdvois, lercvois, lesdvois, lencom, segvar, bhole, itrimat, qmv, iflg, ixs, ixq, pm, x)
Definition aconve.F:54
bconv2
subroutine bconv2(vtot, phi, flux, flu1, ale_connect, bhole, nm)
Definition bconv2.F:32
my_real
#define my_real
Definition cppsort.cpp:32
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
ale_connectivity_mod
Definition ale_connectivity_mod.F:223
ale_mod
Definition ale_mod.F:149
alefvm_mod
Definition alefvm_mod.F:88
alefvm_mod::alefvm_param
type(alefvm_param_), target alefvm_param
Definition alefvm_mod.F:121
i22bufbric_mod
Definition cut-cell-search_mod.F:95
i22bufbric_mod::brick_list
type(brick_entity), dimension(:,:), allocatable, target brick_list
Definition cut-cell-search_mod.F:105
i22tri_mod
Definition cut-cell-search_mod.F:212
initbuf_mod
Definition initbuf.F:160
initbuf_mod::initbuf
subroutine initbuf(iparg, ng, mtn, llt, nft, iad, ity, npt, jale, ismstr, jeul, jtur, jthe, jlag, jmult, jhbe, jivf, mid, jpor, jcvt, jclose, jpla, irep, iint, igtyp, israt, isrot, icsen, isorth, isorthg, ifailure, jsms)
Definition initbuf.F:261
segvar_mod
Definition segvar_mod.F:61
spmd_e1vois
subroutine spmd_e1vois(phi, nercvois, nesdvois, lercvois, lesdvois, lencom)
Definition spmd_cfd.F:375
my_barrier
subroutine my_barrier
Definition machine.F:31
segvar_mod::t_segvar
Definition segvar_mod.F:69
varcondec
subroutine varcondec(icodv, varconv, codtot)
Definition varcondec.F:33