sr8bilan.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"
#include "vect01_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	sr8bilan (partsav, eint, rho, rk, vol, vnew, vxa, vya, vza, va2, iparts, gresav, grth, igrth, xx, yy, zz, xx2, yy2, zz2, xy, yz, zx, itask, iparg, sensors)

Function/Subroutine Documentation

sr8bilan()

subroutine sr8bilan	(		partsav,
			eint,
			rho,
			rk,
			vol,
			vnew,
			vxa,
			vya,
			vza,
			va2,
		integer, dimension(*)	iparts,
			gresav,
		integer, dimension(*)	grth,
		integer, dimension(*)	igrth,
			xx,
			yy,
			zz,
			xx2,
			yy2,
			zz2,
			xy,
			yz,
			zx,
		integer	itask,
		integer, dimension(*)	iparg,
		type (sensors_), intent(inout)	sensors )

Definition at line 33 of file sr8bilan.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE sensor_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.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      "param_c.inc"
#include      "vect01_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
C     REAL
      my_real
     .   partsav(npsav,*), eint(*), rho(*), rk(*), vol(*),
     .   vnew(mvsiz), vxa(mvsiz), vya(mvsiz), vza(mvsiz), va2(mvsiz),
     .   gresav(*), xx(mvsiz), yy(mvsiz), zz(mvsiz),
     .   xx2(mvsiz), yy2(mvsiz), zz2(mvsiz),
     .   xy(mvsiz), yz(mvsiz), zx(mvsiz)
      INTEGER 
     .   IPARTS(*),GRTH(*),IGRTH(*),IPARG(*),ITASK
      type (sensors_),INTENT(INOUT) :: SENSORS
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, M, FLAG
C     REAL
      my_real
     .   xmas(mvsiz),
     .   off(mvsiz), ei(mvsiz) , ek(mvsiz),
     .   xm(mvsiz) , ym(mvsiz) , zm(mvsiz),
     .   xcg(mvsiz), ycg(mvsiz), zcg(mvsiz),
     .   xxm(mvsiz), yym(mvsiz), zzm(mvsiz),
     .   ixx(mvsiz), iyy(mvsiz), izz(mvsiz),
     .   ixy(mvsiz), iyz(mvsiz), izx(mvsiz),
     .   rei(mvsiz), rek(mvsiz), off_l(mvsiz)
C-----------------------------------------------
      flag = iparg(80) ! Flag pour sorties additionnelles
C-----------------------------------------------
C     CALCUL BILAN SORTIES DEF 
C-----------------------------------------------
      DO i=lft,llt
        vxa(i) = vxa(i)*one_over_8
        vya(i) = vya(i)*one_over_8
        vza(i) = vza(i)*one_over_8
        va2(i) = va2(i)*one_over_8
      ENDDO
C
      DO i=lft,llt
        xmas(i)= rho(i)*vnew(i)
        ei(i)  = eint(i)*vol(i)
        ek(i)  = xmas(i)*va2(i)*half
        xm(i)  = xmas(i)*vxa(i)
        ym(i)  = xmas(i)*vya(i)
        zm(i)  = xmas(i)*vza(i)
      ENDDO
C
      DO i=lft,llt
        m=iparts(i)
        partsav(1,m)=partsav(1,m) + ei(i)
        partsav(2,m)=partsav(2,m) + ek(i)
        partsav(3,m)=partsav(3,m) + xm(i)
        partsav(4,m)=partsav(4,m) + ym(i)
        partsav(5,m)=partsav(5,m) + zm(i)
        partsav(6,m)=partsav(6,m) + xmas(i)
      ENDDO
C
C-----------------------------------------------
C     CALCUL BILAN SORTIES ADDITIONNELLES 
C-----------------------------------------------
      IF(flag==1) THEN
        DO i=lft,llt
          xx(i) = xx(i)*one_over_8
          yy(i) = yy(i)*one_over_8
          zz(i) = zz(i)*one_over_8
          xy(i) = xy(i)*one_over_8
          yz(i) = yz(i)*one_over_8
          zx(i) = zx(i)*one_over_8
          xx2(i)= xx2(i)*one_over_8
          yy2(i)= yy2(i)*one_over_8
          zz2(i)= zz2(i)*one_over_8
        ENDDO
C
        DO i=lft,llt
          xcg(i)= xmas(i)*xx(i)
          ycg(i)= xmas(i)*yy(i)
          zcg(i)= xmas(i)*zz(i)
          ixy(i)= -xmas(i)*xy(i)
          iyz(i)= -xmas(i)*yz(i)
          izx(i)= -xmas(i)*zx(i)
          ixx(i)= xmas(i)*(yy2(i) + zz2(i))
          iyy(i)= xmas(i)*(zz2(i) + xx2(i))
          izz(i)= xmas(i)*(xx2(i) + yy2(i))
          xxm(i)= vza(i)*ycg(i)-vya(i)*zcg(i)
          yym(i)= vxa(i)*zcg(i)-vza(i)*xcg(i)
          zzm(i)= vya(i)*xcg(i)-vxa(i)*ycg(i)
          rei(i)= zero
          rek(i)= zero
        ENDDO
C
        DO i=lft,llt
          m=iparts(i)
          partsav(9,m) =partsav(9,m)  + xcg(i)
          partsav(10,m)=partsav(10,m) + ycg(i)
          partsav(11,m)=partsav(11,m) + zcg(i)
          partsav(12,m)=partsav(12,m) + xxm(i)
          partsav(13,m)=partsav(13,m) + yym(i)
          partsav(14,m)=partsav(14,m) + zzm(i)
          partsav(15,m)=partsav(15,m) + ixx(i)
          partsav(16,m)=partsav(16,m) + iyy(i)
          partsav(17,m)=partsav(17,m) + izz(i)
          partsav(18,m)=partsav(18,m) + ixy(i)
          partsav(19,m)=partsav(19,m) + iyz(i)
          partsav(20,m)=partsav(20,m) + izx(i)
          partsav(21,m)=partsav(21,m) + rei(i)
          partsav(22,m)=partsav(22,m) + rek(i) 
        ENDDO
      ENDIF
C
C-----------------------------------------------
      IF (igre /= 0) THEN
          DO i=lft,llt
            off(i)     = one
          ENDDO
          CALL grelem_sav(lft   ,llt   ,gresav,igrth ,grth  ,
     2                    off   ,ei    ,ek    ,xm    ,ym    ,
     3                    zm    ,xmas  ,xcg   ,ycg   ,zcg   ,
     4                    xxm   ,yym   ,zzm   ,ixx   ,iyy   ,
     5                    izz   ,ixy   ,iyz   ,izx   ,rei   ,
     6                    rek   ,flag  )
      ENDIF
C-----------------------------------------------
C
      IF(jtur/=0)THEN
       DO 50 i=lft,llt
       m=iparts(i)
       partsav(7,m)=partsav(7,m) + rk(i)*vol(i)
   50  CONTINUE
      ENDIF
 
      DO i=lft,llt
            off_l(i) = one
      ENDDO
      CALL sensor_energy_bilan(lft,llt,ei,ek,off_l,iparts,itask,sensors) 
C
      RETURN