Sources Esope | Cast3M

rovap0
C ROVAP0    SOURCE    SERRE     16/11/25    21:15:04     9165           
      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
C
      RAC1 = ((B0*T + B1) * T + B2) * P / (T**4.D0) + (R*R/4.D0)
      IF (RAC1 .GT. 0.D0) THEN
         RAC1 = RAC1 ** 0.5D0
      ELSE
         RAC1 = 0.D0
      ENDIF
      ZVAP1 = RAC1/R + 0.5D0
C
      ROVAP0 = P / (R * T * ZVAP1)
C
      END