sbilan43.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	sbilan43 (nel, iparts, partsav, eint, rho, areat, vxa, vya, vza, va2, fill)

Function/Subroutine Documentation

sbilan43()

subroutine sbilan43	(	integer	nel,
		integer, dimension(*)	iparts,
			partsav,
			eint,
			rho,
			areat,
			vxa,
			vya,
			vza,
			va2,
			fill )

Definition at line 28 of file sbilan43.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NEL
      INTEGER IPARTS(*)
C     REAL
      my_real
     .  partsav(npsav,*), eint(*), rho(*), areat(*),
     .  vxa(*), vya(*), vza(*), va2(*), fill(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,M
C     REAL
      my_real
     .   ei_sav(nel),ek_sav(nel),xmas(nel),
     .   xm_sav(nel),ym_sav(nel),zm_sav(nel)
C-----------------------------------------------
      ei_sav = zero
      ek_sav = zero
      xm_sav = zero
      ym_sav = zero
      zm_sav = zero
      DO i=1,nel
        xmas(i)=fill(i)*rho(i)*areat(i)
      ENDDO
C
      DO i=1,nel
        m = iparts(i)
        partsav(1,m)=partsav(1,m) + fill(i)*eint(i)
        partsav(2,m)=partsav(2,m) + xmas(i)*va2(i)*one_over_16
        partsav(6,m)=partsav(6,m) + xmas(i)
        partsav(3,m)=partsav(3,m) + xmas(i)*vxa(i)*one_over_8
        partsav(4,m)=partsav(4,m) + xmas(i)*vya(i)*one_over_8
        partsav(5,m)=partsav(5,m) + xmas(i)*vza(i)*one_over_8
      ENDDO
C-----------
      RETURN