Esope source files

   1 : C DRVDP1    SOURCE    CHAT      05/01/12    22:55:33     5004           
   2 :       SUBROUTINE DRVDP1(MPOVA1,MPOVA2,MPOVA3,ISAT)
   3 : C--------------------------------------------------------------------
   4 : C Densité de la vapeur d'eau 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 l'enthalpie vapeur spécifique à saturation : les data sont des CHPOS
  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 la densité de la vapeur d'eau
  16 : C                 par rapport à la pression partielle de vapeur
  17 : C                 (en kg/m3/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 = DRVDP0(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 = DRVDP0(X1,X2)
  42 :             MPOVA3.VPOCHA(I,1) = X3
  43 :  20      CONTINUE
  44 :       ENDIF
  45 : C
  46 :       RETURN
  47 :       END
  48 :  
  49 :  
  50 :