Esope source files

   1 : C DHLDP1    SOURCE    CHAT      05/01/12    22:48:51     5004           
   2 :       SUBROUTINE DHLDP1(MPOVA1,MPOVA2,MPOVA3,ISAT)
   3 : C--------------------------------------------------------------------
   4 : C Enthalpie liquide spécifique en fonction de P et de T:
   5 : C dérivée partielle par rapport à la pression
   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 liquide spécifique à saturation: les data sont des CHPOS
  10 : C--------------------------------------------------------------------
  11 : C
  12 : C  E/  MPOVA1   : Pression (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
  16 : C                 liquide par rapport à la pression (en J/kg/Pa)
  17 : C-----------------------------------------------------------------------
  18 : C
  19 : C Langage : ESOPE + FORTRAN77
  20 : C   10/2002
  21 : C
  22 : C-----------------------------------------------------------------------
  23 :       IMPLICIT INTEGER(I-N)
  24 :       IMPLICIT REAL*8 (A-H,O-Z)                 
  25 : C
  26 : -INC SMCHPOI
  27 : C
  28 :       N1 = MPOVA1.VPOCHA(/1)
  29 :       IF (ISAT.EQ.0) THEN
  30 :          DO 10 I=1,N1
  31 :             X1 = MPOVA1.VPOCHA(I,1)
  32 :             X2 = MPOVA2.VPOCHA(I,1)
  33 :             X3 = DHLDP0(X1,X2)
  34 :             MPOVA3.VPOCHA(I,1) = X3
  35 :  10      CONTINUE
  36 :       ELSE
  37 :          DO 20 I=1,N1
  38 :             X1 = MPOVA1.VPOCHA(I,1)
  39 :             X2 = TSATP0(X1)
  40 :             X3 = DHLDP0(X1,X2)
  41 :             MPOVA3.VPOCHA(I,1) = X3
  42 :  20      CONTINUE
  43 :       ENDIF
  44 : C
  45 :       RETURN
  46 :       END
  47 :  
  48 :  
  49 :