qbilan.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	qbilan (partsav, off, eint, rho, rk, vol, vy1, vy2, vy3, vy4, vz1, vz2, vz3, vz4, vnew, ipartq, ehou, r22, r23, r32, r33, gresav, grth, igrth, iexpan, eintth, itask, nel, jtur, jcvt, igre, sensors, g_wpla, wpla)

Function/Subroutine Documentation

qbilan()

subroutine qbilan	(		partsav,
			off,
			eint,
			rho,
			rk,
			vol,
			vy1,
			vy2,
			vy3,
			vy4,
			vz1,
			vz2,
			vz3,
			vz4,
			vnew,
		integer, dimension(*)	ipartq,
			ehou,
			r22,
			r23,
			r32,
			r33,
			gresav,
		integer, dimension(*)	grth,
		integer, dimension(*)	igrth,
		integer	iexpan,
			eintth,
		integer	itask,
		integer, intent(in)	nel,
		integer, intent(in)	jtur,
		integer, intent(in)	jcvt,
		integer, intent(in)	igre,
		type (sensors_), intent(inout)	sensors,
		integer, intent(in)	g_wpla,
		intent(in)	wpla )

Definition at line 35 of file qbilan.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----6------------------------------------------
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-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      INTEGER, INTENT(IN) :: JTUR
      INTEGER, INTENT(IN) :: JCVT
      INTEGER, INTENT(IN) :: IGRE
      my_real
     .   partsav(npsav,*), off(*), eint(*), rho(*), rk(*), vol(*),
     .   vy1(*),vy2(*),vy3(*),vy4(*),vz1(*),vz2(*),vz3(*),vz4(*),
     .   vnew(*), ehou(*), r22(*), r23(*), r32(*), r33(*),gresav(*),
     .   eintth(*)
      INTEGER IEXPAN,ITASK,
     .   IPARTQ(*),GRTH(*),IGRTH(*)
      type (sensors_),INTENT(INOUT) :: SENSORS
      INTEGER, INTENT(IN) :: G_WPLA
      my_real,DIMENSION(NEL*G_WPLA), INTENT(IN) :: wpla
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      my_real
     .  vya(mvsiz), vza(mvsiz), va2(mvsiz), xmas(mvsiz), vy, vz,
     .  off_sav(mvsiz),rbidon(1),ei_sav(mvsiz),ek_sav(mvsiz),
     .  xm_sav(mvsiz),ym_sav(mvsiz),zm_sav(mvsiz),ei(mvsiz),ek(mvsiz)
      INTEGER I,M, FLAG
C-----------------------------------------------
      rbidon = zero
      off_sav= zero
      ei_sav = zero
      ek_sav = zero
      xm_sav = zero
      ym_sav = zero
      zm_sav = zero
      
      DO i=1,nel
        vya(i)=vy1(i)+vy2(i)+vy3(i)+vy4(i)
        vza(i)=vz1(i)+vz2(i)+vz3(i)+vz4(i)
      ENDDO
 
      IF(jcvt/=0)THEN
        DO i=1,nel
          vy=r22(i)*vya(i)+r23(i)*vza(i)
          vz=r32(i)*vya(i)+r33(i)*vza(i)
          vya(i)=vy
          vza(i)=vz
        END DO
      END IF
 
      DO i=1,nel
        va2(i)=
     1         vy1(i)*vy1(i)+vy2(i)*vy2(i)+vy3(i)*vy3(i)+vy4(i)*vy4(i)
     2        +vz1(i)*vz1(i)+vz2(i)*vz2(i)+vz3(i)*vz3(i)+vz4(i)*vz4(i)
      ENDDO
 
      DO i=1,nel
        xmas(i)=rho(i)*vnew(i)*off(i)  
      ENDDO
 
      flag = 0
      IF (igre /= 0) THEN
          DO i=1,nel
            off_sav(i) = one
            ei_sav(i)  = eint(i)*vol(i)
            ek_sav(i)  = xmas(i)*va2(i)*one_over_8
            xm_sav(i)  = zero
            ym_sav(i)  = xmas(i)*vya(i)*fourth
            zm_sav(i)  = xmas(i)*vza(i)*fourth
          ENDDO
          CALL grelem_sav(1   ,nel   ,gresav,igrth ,grth  ,
     2                    off_sav,ei_sav,ek_sav,xm_sav,ym_sav,
     3                    zm_sav,xmas  ,rbidon,rbidon,rbidon,
     4                    rbidon,rbidon,rbidon,rbidon,rbidon,
     5                    rbidon,rbidon,rbidon,rbidon,rbidon,
     6                    rbidon,flag)
      ENDIF
 
      DO i=1,nel
        m=ipartq(i)
        partsav(1,m)=partsav(1,m) + eint(i)*vol(i)*off(i)
        partsav(2,m)=partsav(2,m) + xmas(i)*va2(i)*one_over_8
        partsav(3,m)=zero
        partsav(4,m)=partsav(4,m) + xmas(i)*vya(i)*fourth
        partsav(5,m)=partsav(5,m) + xmas(i)*vza(i)*fourth
        partsav(6,m)=partsav(6,m) + xmas(i)
        IF (g_wpla > 0)  partsav(29,m)=partsav(29,m) + wpla(i)
      ENDDO
 
      IF(jtur/=0)THEN
       DO i=1,nel
         m=ipartq(i)
         partsav(7,m)=partsav(7,m) + rk(i)*vol(i)*off(i)
       ENDDO
      ENDIF
 
      IF(iexpan/=0)THEN
       DO i=1,nel
        m=ipartq(i)
        partsav(27,m)=partsav(27,m) + eintth(i)*vol(i)*off(i)
       ENDDO
      ENDIF
 
      DO i=1,nel
        m=ipartq(i)
        ei(i) = eint(i)*vol(i)*off(i)
        ek(i) = xmas(i)*va2(i)*one_over_8
      ENDDO
      CALL sensor_energy_bilan(1,nel,ei,ek,off,ipartq,itask,sensors)  
 
      RETURN