aflux0.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "com08_c.inc"
#include "vect01_c.inc"
#include "param_c.inc"
#include "task_c.inc"
#include "inter22.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	aflux0 (iparg, elbuf_str, pm, ixs, ixq, v, w, x, flux, flu2, veul, fill, dfill, vol, ale_connect, alph, itask, nv46, ipm)

Function/Subroutine Documentation

aflux0()

subroutine aflux0	(	integer, dimension(nparg,ngroup)	iparg,
		type (elbuf_struct_), dimension(ngroup), target	elbuf_str,
			pm,
		integer, dimension(nixs,numels)	ixs,
		integer, dimension(7,numelq)	ixq,
			v,
			w,
			x,
			flux,
			flu2,
			veul,
			fill,
			dfill,
			vol,
		type(t_ale_connectivity), intent(in)	ale_connect,
			alph,
		integer	itask,
		integer	nv46,
		integer, dimension(npropmi,nummat)	ipm )

Definition at line 46 of file aflux0.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE initbuf_mod
      USE elbufdef_mod            
      USE i22tri_mod        
      USE ale_connectivity_mod
      USE alefvm_mod , only:alefvm_param
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "com04_c.inc"
#include      "com08_c.inc"
#include      "vect01_c.inc"
#include      "param_c.inc"
#include      "task_c.inc"
#include      "inter22.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      my_real pm, v, w, x
      INTEGER IPARG(NPARG,NGROUP), IXS(NIXS,NUMELS), IXQ(7,NUMELQ), IPM(NPROPMI,NUMMAT)
      INTEGER ITASK,NV46
      my_real flux(*), flu2(*), veul(*), fill(numnod,*),dfill(numnod,*), vol(*), alph(2,*)
      TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_STR    
      TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECT
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NG, I, NF1,ISOLNOD,NSG,NEL
      my_real udt
      TYPE(G_BUFEL_)  ,POINTER :: GBUF     
      TYPE(L_BUFEL_)  ,POINTER :: LBUF1,LBUF2    
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
      IF(itask == 0) THEN
       IF (nmult /= 0) THEN
         DO ng=1,ngroup
C     ALE ON / OFF
            IF (iparg(76, ng)  ==  1) cycle ! --> OFF
           CALL initbuf(iparg    ,ng      ,
     2        mtn     ,llt     ,nft     ,iad     ,ity     ,
     3        npt     ,jale    ,ismstr  ,jeul    ,jtur    ,
     4        jthe    ,jlag    ,jmult   ,jhbe    ,jivf    ,
     5        nvaux   ,jpor    ,jcvt    ,jclose  ,jplasol ,
     6        irep    ,iint    ,igtyp   ,israt   ,isrot   ,
     7        icsen   ,isorth  ,isorthg ,ifailure,jsms    )
          IF (jmult == 0) cycle
          gbuf  => elbuf_str(ng)%GBUF
          lbuf1 => elbuf_str(ng)%BUFLY(1)%LBUF(1,1,1)
          lbuf2 => elbuf_str(ng)%BUFLY(2)%LBUF(1,1,1)
          lft=1
          IF (dt1 > zero)THEN
            udt=one/dt1
                ELSE
                 udt=zero
                ENDIF
                DO i=lft,llt
                  alph(1,i+nft)=lbuf1%FRAC(i)
                  alph(2,i+nft)=lbuf2%FRAC(i)
                  vol(i+nft)=gbuf%VOL(i)*udt
          ENDDO
        ENDDO
c
       ENDIF
      ENDIF
C--------------------
      CALL my_barrier
C--------------------
      DO ng=itask+1,ngroup,nthread
C     ALE ON / OFF 
         IF (iparg(76, ng)  ==  1) cycle ! --> OFF
           CALL initbuf(iparg    ,ng      ,
     2        mtn     ,llt     ,nft     ,iad     ,ity     ,
     3        npt     ,jale    ,ismstr  ,jeul    ,jtur    ,
     4        jthe    ,jlag    ,jmult   ,jhbe    ,jivf    ,
     5        nvaux   ,jpor    ,jcvt    ,jclose  ,jplasol ,
     6        irep    ,iint    ,igtyp   ,israt   ,isrot   ,
     7        icsen   ,isorth  ,isorthg ,ifailure,jsms    )
        IF (jale+jeul == 0)   cycle
        IF (iparg(8,ng) == 1) cycle
        gbuf => elbuf_str(ng)%GBUF
        nsg = iparg(10,ng)
        isolnod = iparg(28,ng)
        lft=1
        nf1=nft+1
        nel = llt
c
        IF(n2d == 0.AND.jeul == 0.AND. isolnod == 4)THEN
         CALL a4flux3(pm,ixs,v,w,x,flux(6*nft+1),flu2(nf1),ale_connect)
        ELSEIF(n2d == 0.AND.jeul == 0)THEN
         IF(alefvm_param%IEnabled==0)THEN        
           CALL aflux3(pm,ixs,v,w,x,flux(6*nft+1),flu2(nf1),ale_connect,nsg,gbuf%TAG22)
         ELSE
           CALL alefvm_aflux3(pm       ,ixs         ,w  ,flux(6*nft+1),
     2                        flu2(nf1),ale_connect,
     3                        ipm      ,nv46       ,x   ,
     4                        nel      )        
         ENDIF
        ELSEIF(n2d /= 0.AND.jeul == 0)THEN
         CALL aflux2(pm,ixq,v,w,x,
     .               flux(4*nft+1),flu2(nf1),
     .               fill,dfill,vol,ale_connect,alph)
        ELSEIF(n2d == 0.AND.jeul /= 0)THEN
         IF(alefvm_param%IEnabled==0)THEN
           CALL eflux3(pm,ixs,v,flux(6*nft+1),flu2(nf1),veul,ale_connect,gbuf%TAG22)
         ELSE
           CALL alefvm_eflux3(pm   ,ixs          ,flux(6*nft+1), flu2(nf1),
     2                        veul ,ale_connect  ,ipm          ,
     3                        nv46 ,x            ,nel      )
         ENDIF
        ELSEIF(n2d /= 0.AND.jeul /= 0)THEN
         CALL eflux2(pm,ixq,v,x,flux(4*nft+1),flu2(nf1),fill,dfill,vol,ale_connect,alph)
        ENDIF
      ENDDO
      
C-----------
      CALL my_barrier
C-----------      
      
      IF(int22 > 0)THEN
        IF(i22_aleul == 2)THEN
          CALL eflux3_int22_fvm(
     .      pm         , ixs          , flux(1)   , flu2(1)   ,
     .      iparg      , elbuf_str    , itask     ,
     .      nv46       , ipm          , x)         
        ELSE
          CALL aflux3_int22_fvm(
     .      pm         , ixs          , flux(1)   , flu2(1)   ,
     .      iparg      , elbuf_str    , itask     ,
     .      nv46       , ipm          , x            , w)         
       ENDIF
      ENDIF
      
C-----------      
      CALL my_barrier      
C-----------
 
      RETURN