C TOCON1 SOURCE CHAT 05/01/13 03:42:39 5004 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