Sources Esope | Cast3M

pvap0
C PVAP0     SOURCE    SERRE     16/11/25    21:15:04     9165           
      FUNCTION PVAP0(RO,T)
C
C---------------------------------------------------------------------
C Pression partielle de vapeur fonction de RO et de T (issue de RALOC)
C---------------------------------------------------------------------
C
C  E/  RO    : Densité de vapeur (en kg/m3)
C  E/  T     : Température (en K)
C   /S PVAP0 : Pression partielle de vapeur (p = ro R Z T)
C
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.513D0
C
      ZVAP1 = ((B0*T + B1) * T + B2) * RO / (R * T**3.D0) + 1.D0
      ZVAP0 = MAX(ZVAP1,0.5D0)
      PVAP0 = ZVAP0 * RO * R * T
C
      END