Esope source files

   1 : C ROVAP0    SOURCE    TTMF3     12/07/05    21:16:12     7425           
   2 :       FUNCTION ROVAP0(P,T)
   3 : C--------------------------------------------------------------------
   4 : C Densité vapeur fonction de P et de T (issu de RALOC)
   5 : C--------------------------------------------------------------------
   6 : C
   7 : C  E/  P      : Pression partielle de vapeur (en Pa)
   8 : C  E/  T      : Température (en K)
   9 : C   /S ROVAP0 : Densité de la vapeur (en kg/m3)
  10 : C
  11 : C--------------------------------------------------------------------
  12 :       IMPLICIT INTEGER(I-N)
  13 :       IMPLICIT REAL*8 (A-H,O-Z)                 
  14 : C
  15 :       B0 = -4512.70756D0  
  16 :       B1 = 4245585.9D0    
  17 :       B2 = -1584624342.D0 
  18 :       R  = 461.523D0      
  19 : *
  20 :       DNUM1 = B0*T*T + (B1*T) + B2
  21 :       VAL1   = LOG((R*R/4.D0) + (P*DNUM1/T/T/T/T))
  22 :       VAL2   = EXP(VAL1/2.D0)
  23 :       ROVAP0 = 1.D0/((VAL2 + R/2.D0) * T / P)
  24 : C
  25 :       END
  26 :  
  27 :  
  28 :