def test_water_at_critcalPoint(self):
    '\n        Reproduce verification results from IAPWS-IF97 for water at 647.096K\n        http://www.iapws.org/relguide/supsat.pdf\n        '
    water = SaturatedWater()
    steam = SaturatedSteam()
    Tk = 647.096
    ref_vapor_pressure = 22064000.0
    ref_dp_dT = 268000.0
    ref_saturated_water_rho = 322
    ref_saturated_steam_rho = 322
    ref_alpha = 1548000.0
    ref_saturated_water_enthalpy = 2086600.0
    ref_saturated_steam_enthalpy = 2086600.0
    ref_phi = 3578.0
    ref_saturated_water_entropy = 4410.0
    ref_saturated_steam_entropy = 4410.0
    self.assertAlmostEqual((ref_vapor_pressure / water.vaporPressure(Tk=Tk)), 1, 3)
    self.assertAlmostEqual((ref_vapor_pressure / steam.vaporPressure(Tk=Tk)), 1, 3)
    self.assertAlmostEqual((ref_dp_dT / water.vaporPressurePrime(Tk=Tk)), 1, 3)
    self.assertAlmostEqual((ref_dp_dT / steam.vaporPressurePrime(Tk=Tk)), 1, 3)
    self.assertAlmostEqual(ref_saturated_water_rho, water.densityKgM3(Tk=Tk), 0)
    self.assertAlmostEqual(ref_saturated_steam_rho, steam.densityKgM3(Tk=Tk), 0)
    self.assertAlmostEqual((ref_alpha / water.auxiliaryQuantitySpecificEnthalpy(Tk=Tk)), 1, 3)
    self.assertAlmostEqual((ref_alpha / steam.auxiliaryQuantitySpecificEnthalpy(Tk=Tk)), 1, 3)
    self.assertAlmostEqual((ref_saturated_water_enthalpy / water.enthalpy(Tk=Tk)), 1, 2)
    self.assertAlmostEqual((ref_saturated_steam_enthalpy / steam.enthalpy(Tk=Tk)), 1, 2)
    self.assertAlmostEqual((ref_phi / water.auxiliaryQuantitySpecificEntropy(Tk=Tk)), 1, 3)
    self.assertAlmostEqual((ref_phi / steam.auxiliaryQuantitySpecificEntropy(Tk=Tk)), 1, 3)
    self.assertAlmostEqual((ref_saturated_water_entropy / water.entropy(Tk=Tk)), 1, 3)
    self.assertAlmostEqual((ref_saturated_steam_entropy / steam.entropy(Tk=Tk)), 1, 3)
Reproduce verification results from IAPWS-IF97 for water at 647.096K
http://www.iapws.org/relguide/supsat.pdf
def test_water_at_critcalPoint(self):
    '\n        Reproduce verification results from IAPWS-IF97 for water at 647.096K\n        http://www.iapws.org/relguide/supsat.pdf\n        '
    water = SaturatedWater()
    steam = SaturatedSteam()
    Tk = 647.096
    ref_vapor_pressure = 22064000.0
    ref_dp_dT = 268000.0
    ref_saturated_water_rho = 322
    ref_saturated_steam_rho = 322
    ref_alpha = 1548000.0
    ref_saturated_water_enthalpy = 2086600.0
    ref_saturated_steam_enthalpy = 2086600.0
    ref_phi = 3578.0
    ref_saturated_water_entropy = 4410.0
    ref_saturated_steam_entropy = 4410.0
    self.assertAlmostEqual((ref_vapor_pressure / water.vaporPressure(Tk=Tk)), 1, 3)
    self.assertAlmostEqual((ref_vapor_pressure / steam.vaporPressure(Tk=Tk)), 1, 3)
    self.assertAlmostEqual((ref_dp_dT / water.vaporPressurePrime(Tk=Tk)), 1, 3)
    self.assertAlmostEqual((ref_dp_dT / steam.vaporPressurePrime(Tk=Tk)), 1, 3)
    self.assertAlmostEqual(ref_saturated_water_rho, water.densityKgM3(Tk=Tk), 0)
    self.assertAlmostEqual(ref_saturated_steam_rho, steam.densityKgM3(Tk=Tk), 0)
    self.assertAlmostEqual((ref_alpha / water.auxiliaryQuantitySpecificEnthalpy(Tk=Tk)), 1, 3)
    self.assertAlmostEqual((ref_alpha / steam.auxiliaryQuantitySpecificEnthalpy(Tk=Tk)), 1, 3)
    self.assertAlmostEqual((ref_saturated_water_enthalpy / water.enthalpy(Tk=Tk)), 1, 2)
    self.assertAlmostEqual((ref_saturated_steam_enthalpy / steam.enthalpy(Tk=Tk)), 1, 2)
    self.assertAlmostEqual((ref_phi / water.auxiliaryQuantitySpecificEntropy(Tk=Tk)), 1, 3)
    self.assertAlmostEqual((ref_phi / steam.auxiliaryQuantitySpecificEntropy(Tk=Tk)), 1, 3)
    self.assertAlmostEqual((ref_saturated_water_entropy / water.entropy(Tk=Tk)), 1, 3)
    self.assertAlmostEqual((ref_saturated_steam_entropy / steam.entropy(Tk=Tk)), 1, 3)