prelecflow.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
C            CFLUID                          /resol/prelecfluid.F
!||====================================================================
!||    prelecflow     ../engine/source/elements/solid/solide/prelecflow.f
!||--- called by ------------------------------------------------------
!||    forint         ../engine/source/elements/forint.F
!||--- calls      -----------------------------------------------------
!||    initbuf        ../engine/share/resol/initbuf.F
!||    my_barrier     ../engine/source/system/machine.F
!||    spmd_e1vois    ../engine/source/mpi/fluid/spmd_cfd.F
!||    startimeg      ../engine/source/system/timer.F
!||    stoptimeg      ../engine/source/system/timer.F
!||    vfluid         ../engine/source/elements/solid/solide/vfluid.F
!||--- uses       -----------------------------------------------------
!||    elbufdef_mod   ../common_source/modules/mat_elem/elbufdef_mod.F90
!||    initbuf_mod    ../engine/share/resol/initbuf.F
!||====================================================================
      SUBROUTINE prelecflow(ELBUF_TAB,NGROUNC,IGROUNC,IPARG  ,NBSDVOIS,
     1                      NERCVOIS,NESDVOIS,LERCVOIS,LESDVOIS,PHI1 ,
     2                      PHI2  ,POR)
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE initbuf_mod
      USE elbufdef_mod            
C----6---------------------------------------------------------------7---------8
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "param_c.inc"
#include      "vect01_c.inc"
#include      "task_c.inc"
C-----------------------------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IPARG(NPARG,*),NBSDVOIS(*),IGROUNC(*),
     .   NERCVOIS(*),NESDVOIS(*),LERCVOIS(*),LESDVOIS(*),NGROUNC
C     REAL
      my_real
     .   phi1(*),phi2(*),por(*)
      TYPE(elbuf_struct_), TARGET, DIMENSION(NGROUP) :: ELBUF_TAB
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER INDXOF(MVSIZ)
      INTEGER I,II,J,N, NG, NVC, MLW, JFT, JLT,ISOLNOD,ITHK,IPLA,
     .        lencom,kad,nel,nf1,ig
      my_real,
     .   DIMENSION(:), POINTER :: uvar
C======================================================================|
C   
      lencom = 0    
      IF(nspmd > 1)lencom=nercvois(nspmd+1)+nesdvois(nspmd+1) 
C Boucle parallele dynamique SMP
C
!$OMP DO 
 
       DO ig = 1, ngrounc
            ng = igrounc(ig)
            IF(iparg(8,ng)==1)GOTO 150
            ity   =iparg(5,ng)
            IF (iddw>0) CALL startimeg(ng)
            mlw     = iparg(1,ng)
C MLW = 777 ---->cell open foam       
           IF (mlw /= 77 ) GOTO 150
C           
           CALL initbuf(iparg    ,ng      ,
     2        mlw     ,nel     ,nft     ,kad     ,ity     ,
     3        npt     ,jale    ,ismstr  ,jeul    ,jtur    ,
     4        jthe    ,jlag    ,jmult   ,jhbe    ,jivf    ,
     5        nvaux   ,jpor    ,jcvt    ,jclose  ,ipla    ,
     6        irep    ,iint    ,igtyp   ,israt   ,isrot   ,
     7        icsen   ,isorth  ,isorthg ,ifailure,jsms)
C   
            isolnod = iparg(28,ng)
            
C---------------------------------
            lft   = 1 
            llt   = min(nvsiz,nel)
            jft=lft
            jlt=llt
            nf1 = nft+1
            iad = kad
            jsph=0
C----6---------------------------------------------------------------7---------8
        IF(ity==1 .AND. jlag==1)THEN
C-------------------------
C          STANDARD 8 NODE SOLID ELEMENT (JHBE = 1,2,101,102)
C-------------------------
           uvar => elbuf_tab(ng)%BUFLY(1)%MAT(1,1,1)%VAR
           CALL vfluid(
     1   uvar,    nel,     mlw,     phi1,
     2   phi2,    por,     nft)
        ENDIF
C-------------------------
        IF (iddw>0) CALL stoptimeg(ng)
 
 150     CONTINUE 
      ENDDO           
C     
      CALL my_barrier
C
C-----------------------------
C       ECHANGES SPMD  SUR ELTS VOISINS
C-----------------------------       
      IF (nspmd > 1) THEN       
cc        LENCOMN=NBRCVOIS(MAXPROC+1)+NBSDVOIS(MAXPROC+1)
C
C Partie non parallele smt
C
!$OMP SINGLE
       CALL spmd_e1vois(phi1     ,nercvois,nesdvois,lercvois,
     1                     lesdvois,lencom  )
C       
       CALL spmd_e1vois(phi2     ,nercvois,nesdvois,lercvois,
     1                     lesdvois,lencom  ) 
C       
       CALL spmd_e1vois(por     ,nercvois,nesdvois,lercvois,
     1                     lesdvois,lencom  )
C Fin Partie non parallele smt
!$OMP END SINGLE       
      ENDIF
      END