prelecflow.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"
#include "mvsiz_p.inc"
#include "com01_c.inc"
#include "param_c.inc"
#include "vect01_c.inc"
#include "task_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	prelecflow (elbuf_tab, ngrounc, igrounc, iparg, nbsdvois, nercvois, nesdvois, lercvois, lesdvois, phi1, phi2, por)

Function/Subroutine Documentation

prelecflow()

subroutine prelecflow	(	type(elbuf_struct_), dimension(ngroup), target	elbuf_tab,
		integer	ngrounc,
		integer, dimension(*)	igrounc,
		integer, dimension(nparg,*)	iparg,
		integer, dimension(*)	nbsdvois,
		integer, dimension(*)	nercvois,
		integer, dimension(*)	nesdvois,
		integer, dimension(*)	lercvois,
		integer, dimension(*)	lesdvois,
			phi1,
			phi2,
			por )

Definition at line 39 of file prelecflow.F.

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