fvinjt8_1.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	fvinjt8_1 (njet, ibagjet, rbagjet, igeo, geo, pm, ivolu, rvolu, dmout, dhout)

Function/Subroutine Documentation

fvinjt8_1()

subroutine fvinjt8_1	(	integer	njet,
		integer, dimension(nibjet,*)	ibagjet,
			rbagjet,
		integer, dimension(npropgi,*)	igeo,
			geo,
			pm,
		integer, dimension(*)	ivolu,
			rvolu,
			dmout,
			dhout )

Definition at line 29 of file fvinjt8_1.F.

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      "com08_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NJET, IBAGJET(NIBJET,*), IGEO(NPROPGI,*), IVOLU(*)
C     REAL
      my_real
     .   dmout, dhout, rbagjet(nrbjet,*),
     .   geo(npropg,*), pm(npropm,*), rvolu(*)
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER IJ, IK, KK, I_INJ, I_TYPINJ, NGASES, I_GAS
      my_real
     .   r_igc1, rmwi, rmwg, rnm, rnmi, rnmg, mw, gmtot, gmi, dgmout, dgeout
C----------------------------------------------------------------------------
      r_igc1=pm(27,ivolu(66))
      rmwi  = rvolu(10)
C--------------------------
C     NOMBRE TOTAL DE MOLES
C--------------------------
      gmi = rvolu(11)
      rnm = gmi*rmwi
      DO ij=1,njet
        i_inj = iabs(ibagjet(13,ij))
        i_typinj = igeo(22,i_inj)
        ngases   = igeo(23,i_inj)
        DO ik=1,ngases
          IF (i_typinj==1) THEN
            i_gas = igeo(100+(ik-1)*3+1,i_inj)
          ELSE IF (i_typinj==2) THEN
            i_gas = igeo(100+(ik-1)*2+1,i_inj)
          END IF
          mw   = pm(20,i_gas)
          rmwg = r_igc1/mw
          kk=20+(ik-1)*4
          gmtot= rbagjet(kk+1,ij)
          rnmg = gmtot*rmwg
          rnm  = rnm+rnmg
        ENDDO
      ENDDO
      rvolu(14)=rnm
C-------------------------------
C     MASSE et ENTHALPY par GAZ
C-------------------------------
      rnmi=gmi*rmwi
      dgmout=dmout*rnmi/rnm
      dgeout=dhout*rnmi/rnm
C      DGEIN =DGMIN*RVOLU(63)
      rvolu(22)=dgeout
      rvolu(24)=dgmout
C      RVOLU(64)=DGMIN
C      RVOLU(65)=DGEIN
      rvolu(11)=rvolu(11)+dgmout*dt1
C
      DO ij=1,njet
         i_inj = iabs(ibagjet(13,ij))
         i_typinj = igeo(22,i_inj)
         ngases   = igeo(23,i_inj)
         rbagjet( 9,ij)=zero
         rbagjet(10,ij)=zero
         DO ik=1,ngases
          IF (i_typinj==1) THEN
            i_gas = igeo(100+(ik-1)*3+1,i_inj)
          ELSE IF (i_typinj==2) THEN
            i_gas = igeo(100+(ik-1)*2+1,i_inj)
          END IF
          mw   = pm(20,i_gas)
          rmwg = r_igc1/mw
          kk=20+(ik-1)*4
          gmtot= rbagjet(kk+1,ij)
          rnmg  =gmtot*rmwg
          dgmout=dmout*rnmg/rnm
          dgeout=dhout*rnmg/rnm
          rbagjet(kk+1,ij)=gmtot+dgmout*dt1
          rbagjet(kk+3,ij)=dgmout
          rbagjet(kk+4,ij)=dgeout
          rbagjet( 8,ij)=rbagjet( 8,ij)+dgmout*dt1
          rbagjet( 9,ij)=rbagjet( 9,ij)+dgmout
          rbagjet(10,ij)=rbagjet(10,ij)+dgeout
         ENDDO
      ENDDO
      RETURN