multi_bilan.F File Reference

#include "implicit_f.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	multi_bilan (partsav, vol, multi_fvm, nel, nft, iparts, gresav, igrth, grth)

Function/Subroutine Documentation

multi_bilan()

subroutine multi_bilan	(	dimension(npsav, *), intent(inout)	partsav,
		dimension(*), intent(in)	vol,
		type(multi_fvm_struct), intent(in)	multi_fvm,
		integer, intent(in)	nel,
		integer, intent(in)	nft,
		integer, dimension(*), intent(in)	iparts,
		dimension(*), intent(inout)	gresav,
		integer, dimension(*), intent(in)	igrth,
		integer, dimension(*), intent(in)	grth )

Definition at line 31 of file multi_bilan.F.

C-----------------------------------------------
C     M o d u l e s
C-----------------------------------------------
      USE multi_fvm_mod
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      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      my_real, INTENT(INOUT) :: partsav(npsav, *), gresav(*)
      my_real, INTENT(IN) :: vol(*)
      TYPE(MULTI_FVM_STRUCT), INTENT(IN) :: MULTI_FVM
      INTEGER, INTENT(IN) :: NEL, NFT, IPARTS(*), IGRTH(*), GRTH(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: I, M
      my_real :: ei_sav(nel), ek_sav(nel), 
     .     xm_sav(nel), ym_sav(nel), zm_sav(nel), rbidon(1), va2(nel),
     .     mas(nel), off(nel), vx(nel), vy(nel), vz(nel)
      
      rbidon(1) = zero
      off(1:nel) = one
      DO i = 1, nel
         mas(i) = vol(i) * multi_fvm%RHO(i + nft)
         vx(i) = multi_fvm%VEL(1, i + nft)
         vy(i) = multi_fvm%VEL(2, i + nft)
         vz(i) = multi_fvm%VEL(3, i + nft)
         va2(i) = vx(i)**2 + vy(i)**2 + vz(i)**2
      ENDDO
      m = iparts(1)
      DO i=1,nel
         partsav(1, m) = partsav(1, m) + multi_fvm%EINT(i + nft) * vol(i)
         partsav(2, m) = partsav(2, m) + half * mas(i) * va2(i)
         partsav(6, m) = partsav(6, m) + mas(i)
         partsav(3, m) = partsav(3, m) + mas(i) * vx(i)
         partsav(4, m) = partsav(4, m) + mas(i) * vy(i)
         partsav(5, m) = partsav(5, m) + mas(i) * vz(i)
      ENDDO