Sources Esope | Cast3M

tocon1
C TOCON1    SOURCE    CB215821  20/11/25    13:41:06     10792          
      SUBROUTINE TOCON1(FLOT1,FLOT2,MPOVA1,MPOVA2,MPOVA3,
     S                  MPOVA4,MPOVA5,MPOVA6,MPOVA7)
C--------------------------------------------------------------------
C Répartition vapeur/eau liquide au sein de chaque compartiment et
C données thermodynamiques associées à la vapeur
C--------------------------------------------------------------------
C  E/   FLOT1  : Température max assurant la validité de PSATT
C  E/   FLOT2  : Relaxation pour le titre vapeur
C  E/   MPOVA1 : Température (en K)
C  E/   MPOVA2 : Densité H2O (en kg/m3)
C  E/   MPOVA3 : Titre vapeur initial (-)
C   /S  MPOVA4 : Titre vapeur final (-)
C   /S  MPOVA5 : Densité vapeur (en kg/m3)
C   /S  MPOVA6 : Pression partielle de vapeur (en Pa)
C   /S  MPOVA7 : Température de saturation (en K)
C--------------------------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
C
-INC SMCHPOI
      POINTEUR MPOVA7.MPOVAL
C
      N1 = MPOVA1.VPOCHA(/1)
      SEGINI, MPOVA4=MPOVA1
      SEGINI, MPOVA5=MPOVA1
      SEGINI, MPOVA6=MPOVA1
      SEGINI, MPOVA7=MPOVA1
      DO 10 I=1,N1
      TG  = MIN(FLOT1,MPOVA1.VPOCHA(I,1))
      PS  = PSATT0(TG)
      RS  = ROVAP0(PS,TG)
      XV0 = MPOVA3.VPOCHA(I,1)
      XV1 = MIN(RS/MPOVA2.VPOCHA(I,1),1.D0)
      MPOVA4.VPOCHA(I,1) = XV0 + FLOT2*(XV1-XV0)
      MPOVA5.VPOCHA(I,1) = MPOVA4.VPOCHA(I,1) * MPOVA2.VPOCHA(I,1)
      MPOVA6.VPOCHA(I,1) = PVAP0(MPOVA5.VPOCHA(I,1),MPOVA1.VPOCHA(I,1))
      MPOVA7.VPOCHA(I,1) = TSATP0(MPOVA6.VPOCHA(I,1))
 10   CONTINUE
C
      RETURN
      END