Esope source files

   1 : C HVS0      SOURCE    TTMF3     12/07/05    21:15:55     7425           
   2 :       FUNCTION HVS0(P,T)
   3 : C
   4 : C--------------------------------------------------------------------
   5 : * Enthalpie vapeur spécifique fonction de P et de T (issue de RALOC).
   6 : C--------------------------------------------------------------------
   7 : C
   8 : C  E/  P    : Pression partielle de vapeur (en Pa)
   9 : C  E/  T    : Température (en K)
  10 : C   /S HVS0 : Enthalpie spécifique de la vapeur (en J/kg)
  11 : C
  12 : C--------------------------------------------------------------------
  13 :       IMPLICIT INTEGER(I-N)
  14 :       IMPLICIT REAL*8 (A-H,O-Z)                 
  15 : C
  16 :       A0 = 2009.9595D3     
  17 :       A1 = 1.748383D3      
  18 :       A2 = 0.277768D0/2.D0 
  19 :       A3 = 0.275D-3/3.D0   
  20 :       B0 = -4512.70756D0
  21 :       B1 = 4245585.9D0
  22 :       B2 = -1584624342.D0 
  23 :       R  = 461.523D0
  24 : C
  25 :       DNUM1  = (2.D0*B0*T + (3.D0*B1)) * T + 4.D0*B2
  26 :       DENOM1 = (B0*T +  B1) * T + B2 
  27 :       VAL1   = LOG(R*R/4.D0 + P*DENOM1/T/T/T/T)
  28 :       VAL2   = EXP(VAL1/2.D0)
  29 :       VPMRT1 = (VAL2 - R/2.D0) * T  
  30 :       HVS0   = ((A3*T+A2)*T+A1)*T+A0 + DNUM1/DENOM1*VPMRT1
  31 : C
  32 :       END
  33 :  
  34 :  
  35 :