Sources Esope | Cast3M

zvap0
C ZVAP0     SOURCE    SERRE     16/11/25    21:15:06     9165           
      FUNCTION ZVAP0(RO,T)
C
C---------------------------------------------------------------------
C Correction gaz reel pour la vapeur (issue de RALOC)
C---------------------------------------------------------------------
C
C  E/  RO    : Densité de vapeur (en kg/m3)
C  E/  T     : Température (en K)
C   /S ZVAP0 : Correction gaz reel pour la 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.523D0
      ZVAP1 = ((B0*T + B1) * T + B2) * RO / (R * T**3.D0) + 1.D0
      ZVAP0 = MAX(ZVAP1,0.5D0)
C
      END