Esope source files

   1 : C DRVDT3    SOURCE    CHAT      05/01/12    22:56:05     5004           
   2 :       SUBROUTINE DRVDT3(MLREE1,MLREE2,MLREE3,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 température
   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/  MLREE1   : Pression partielle de vapeur (en Pa)
  13 : C  E/  MLREE2   : Température (en K)
  14 : C  E/  ISAT     : 0 si T est donné, 1 sinon (T=TSAT(P))
  15 : C   /S MLREE3   : Dérivée partielle de la densité de la vapeur d'eau
  16 : C                 par rapport à la température (en kg/m3/K)
  17 : C
  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 SMLREEL
  28 : C
  29 :       JG1 = MLREE1.PROG(/1)
  30 :       IF (ISAT.EQ.0) THEN
  31 :          DO 10 I=1,JG1
  32 :             X1 = MLREE1.PROG(I)
  33 :             X2 = MLREE2.PROG(I)
  34 :             X3 = DRVDT0(X1,X2)
  35 :             MLREE3.PROG(I) = X3
  36 :  10      CONTINUE
  37 :       ELSE
  38 :          DO 20 I=1,JG1
  39 :             X1 = MLREE1.PROG(I)
  40 :             X2 = TSATP0(X1)
  41 :             X3 = DRVDT0(X1,X2)
  42 :             MLREE3.PROG(I) = X3
  43 :  20      CONTINUE
  44 :       ENDIF
  45 : C
  46 :       RETURN
  47 :       END
  48 :  
  49 :  
  50 :