C DHVDP3 SOURCE CHAT 05/01/12 22:49:08 5004 SUBROUTINE DHVDP3(MLREE1,MLREE2,MLREE3,ISAT) C-------------------------------------------------------------------- C Enthalpie vapeur spécifique en fonction de P et de T: C dérivée partielle par rapport à la pression partielle de vapeur C C Dans le cas ou seule la pression est fournie, on considère que la C pression P transmise est Psat et on calcule Tsat afin d'évaluer C la fonction derivée C-------------------------------------------------------------------- C C E/ MLREE1 : Pression partielle de vapeur (en Pa) C E/ MLREE2 : Température (en K) C E/ ISAT : 0 si T est donné, 1 sinon (T=TSAT(P)) C /S MLREE3 : Dérivée partielle de l'enthalpie spécifique de la C vapeur par rapport à la pression partielle de C vapeur (en J/kg/Pa) C----------------------------------------------------------------------- C C Langage : ESOPE + FORTRAN77 C 10/2002 C C----------------------------------------------------------------------- IMPLICIT INTEGER(I-N) IMPLICIT REAL*8 (A-H,O-Z) C -INC SMLREEL C JG1 = MLREE1.PROG(/1) IF (ISAT.EQ.0) THEN DO 10 I=1,JG1 X1 = MLREE1.PROG(I) X2 = MLREE2.PROG(I) X3 = DHVDP0(X1,X2) MLREE3.PROG(I) = X3 10 CONTINUE ELSE DO 20 I=1,JG1 X1 = MLREE1.PROG(I) X2 = TSATP0(X1) X3 = DHVDP0(X1,X2) MLREE3.PROG(I) = X3 20 CONTINUE ENDIF C RETURN END