Esope source files

   1 : C TOCON1    SOURCE    CHAT      05/01/13    03:42:39     5004           
   2 :       SUBROUTINE TOCON1(FLOT1,FLOT2,MPOVA1,MPOVA2,MPOVA3,
   3 :      S                  MPOVA4,MPOVA5,MPOVA6,MPOVA7)
   4 : C--------------------------------------------------------------------
   5 : C Répartition vapeur/eau liquide au sein de chaque compartiment et
   6 : C données thermodynamiques associées à la vapeur
   7 : C--------------------------------------------------------------------
   8 : C  E/   FLOT1  : Température max assurant la validité de PSATT
   9 : C  E/   FLOT2  : Relaxation pour le titre vapeur
  10 : C  E/   MPOVA1 : Température (en K)
  11 : C  E/   MPOVA2 : Densité H2O (en kg/m3)
  12 : C  E/   MPOVA3 : Titre vapeur initial (-)
  13 : C   /S  MPOVA4 : Titre vapeur final (-)
  14 : C   /S  MPOVA5 : Densité vapeur (en kg/m3)
  15 : C   /S  MPOVA6 : Pression partielle de vapeur (en Pa)
  16 : C   /S  MPOVA7 : Température de saturation (en K)
  17 : C--------------------------------------------------------------------
  18 :       IMPLICIT INTEGER(I-N)
  19 :       IMPLICIT REAL*8 (A-H,O-Z)                 
  20 : C
  21 : -INC SMCHPOI
  22 :       POINTEUR MPOVA7.MPOVAL
  23 : C
  24 :       N1 = MPOVA1.VPOCHA(/1)
  25 :       SEGINI, MPOVA4=MPOVA1
  26 :       SEGINI, MPOVA5=MPOVA1
  27 :       SEGINI, MPOVA6=MPOVA1
  28 :       SEGINI, MPOVA7=MPOVA1
  29 :       DO 10 I=1,N1
  30 :       TG  = MIN(FLOT1,MPOVA1.VPOCHA(I,1))
  31 :       PS  = PSATT0(TG)
  32 :       RS  = ROVAP0(PS,TG)
  33 :       XV0 = MPOVA3.VPOCHA(I,1)
  34 :       XV1 = MIN(RS/MPOVA2.VPOCHA(I,1),1.D0)
  35 :       MPOVA4.VPOCHA(I,1) = XV0 + FLOT2*(XV1-XV0)
  36 :       MPOVA5.VPOCHA(I,1) = MPOVA4.VPOCHA(I,1) * MPOVA2.VPOCHA(I,1)
  37 :       MPOVA6.VPOCHA(I,1) = PVAP0(MPOVA5.VPOCHA(I,1),MPOVA1.VPOCHA(I,1))
  38 :       MPOVA7.VPOCHA(I,1) = TSATP0(MPOVA6.VPOCHA(I,1))
  39 :  10   CONTINUE
  40 : C
  41 :       RETURN
  42 :       END
  43 :  
  44 :