C ROVAP0    SOURCE    TTMF3     12/07/05    21:16:12     7425           
      FUNCTION ROVAP0(P,T)
C--------------------------------------------------------------------
C Densité vapeur fonction de P et de T (issu de RALOC)
C--------------------------------------------------------------------
C
C  E/  P      : Pression partielle de vapeur (en Pa)
C  E/  T      : Température (en K)
C   /S ROVAP0 : Densité de la vapeur (en kg/m3)
C
C--------------------------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)                 
C
      B0 = -4512.70756D0  
      B1 = 4245585.9D0    
      B2 = -1584624342.D0 
      R  = 461.523D0      
*
      DNUM1 = B0*T*T + (B1*T) + B2
      VAL1   = LOG((R*R/4.D0) + (P*DNUM1/T/T/T/T))
      VAL2   = EXP(VAL1/2.D0)
      ROVAP0 = 1.D0/((VAL2 + R/2.D0) * T / P)
C
      END