Esope source files

   1 : C DHVDP1    SOURCE    CHAT      05/01/12    22:49:05     5004           
   2 :       SUBROUTINE DHVDP1(MPOVA1,MPOVA2,MPOVA3,ISAT)
   3 : C--------------------------------------------------------------------
   4 : C Enthalpie vapeur spécifique en fonction de P et de T:
   5 : C dérivée partielle par rapport à la pression partielle de vapeur
   6 : C
   7 : C Dans le cas ou seule la pression est fournie, on considère que la
   8 : C pression P transmise est Psat et on calcule Tsat afin d'évaluer
   9 : C la fonction derivée
  10 : C--------------------------------------------------------------------
  11 : C
  12 : C  E/  MPOVA1   : Pression partielle de vapeur (en Pa)
  13 : C  E/  MPOVA2   : Température (en K)
  14 : C  E/  ISAT     : 0 si T est donné, 1 sinon (T=TSAT(P))
  15 : C   /S MPOVA3   : Dérivée partielle de l'enthalpie spécifique de la
  16 : C                 vapeur par rapport à la pression partielle de
  17 : C                 vapeur (en J/kg/Pa)
  18 : C-----------------------------------------------------------------------
  19 : C
  20 : C Langage : ESOPE + FORTRAN77
  21 : C   10/2002
  22 : C
  23 : C-----------------------------------------------------------------------
  24 :       IMPLICIT INTEGER(I-N)
  25 :       IMPLICIT REAL*8 (A-H,O-Z)                 
  26 : C
  27 : -INC SMCHPOI
  28 : C
  29 :       N1 = MPOVA1.VPOCHA(/1)
  30 :       IF (ISAT.EQ.0) THEN
  31 :          DO 10 I=1,N1
  32 :             X1 = MPOVA1.VPOCHA(I,1)
  33 :             X2 = MPOVA2.VPOCHA(I,1)
  34 :             X3 = DHVDP0(X1,X2)
  35 :             MPOVA3.VPOCHA(I,1) = X3
  36 :  10      CONTINUE
  37 :       ELSE
  38 :          DO 20 I=1,N1
  39 :             X1 = MPOVA1.VPOCHA(I,1)
  40 :             X2 = TSATP0(X1)
  41 :             X3 = DHVDP0(X1,X2)
  42 :             MPOVA3.VPOCHA(I,1) = X3
  43 :  20      CONTINUE
  44 :       ENDIF
  45 : C
  46 :       RETURN
  47 :       END
  48 :  
  49 :  
  50 :