def test_water_at_freezing(self):
    '\n        Reproduce verification results from IAPWS-IF97 for water at 0C\n        http://www.iapws.org/relguide/supsat.pdf\n        '
    water = SaturatedWater()
    steam = SaturatedSteam()
    Tk = 273.16
    ref_vapor_pressure = 611.657
    ref_dp_dT = 44.436693
    ref_saturated_water_rho = 999.789
    ref_saturated_steam_rho = 0.00485426
    ref_alpha = (- 11.529101)
    ref_saturated_water_enthalpy = 0.611786
    ref_saturated_steam_enthalpy = 2500500.0
    ref_phi = (- 0.04)
    ref_saturated_water_entropy = 0
    ref_saturated_steam_entropy = 9154.0
    self.assertAlmostEqual(ref_vapor_pressure, water.vaporPressure(Tk=Tk), 3)
    self.assertAlmostEqual(ref_vapor_pressure, steam.vaporPressure(Tk=Tk), 3)
    self.assertAlmostEqual(ref_dp_dT, water.vaporPressurePrime(Tk=Tk), 3)
    self.assertAlmostEqual(ref_dp_dT, steam.vaporPressurePrime(Tk=Tk), 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), 3)
    self.assertAlmostEqual(ref_alpha, steam.auxiliaryQuantitySpecificEnthalpy(Tk=Tk), 3)
    self.assertAlmostEqual(ref_saturated_water_enthalpy, water.enthalpy(Tk=Tk), 2)
    self.assertAlmostEqual((ref_saturated_steam_enthalpy / steam.enthalpy(Tk=Tk)), 1, 2)
    self.assertAlmostEqual(ref_phi, water.auxiliaryQuantitySpecificEntropy(Tk=Tk), 2)
    self.assertAlmostEqual(ref_phi, steam.auxiliaryQuantitySpecificEntropy(Tk=Tk), 2)
    self.assertAlmostEqual(ref_saturated_water_entropy, water.entropy(Tk=Tk), 3)
    self.assertAlmostEqual((ref_saturated_steam_entropy / steam.entropy(Tk=Tk)), 1, 3)
Reproduce verification results from IAPWS-IF97 for water at 0C
http://www.iapws.org/relguide/supsat.pdf
def test_water_at_freezing(self):
    '\n        Reproduce verification results from IAPWS-IF97 for water at 0C\n        http://www.iapws.org/relguide/supsat.pdf\n        '
    water = SaturatedWater()
    steam = SaturatedSteam()
    Tk = 273.16
    ref_vapor_pressure = 611.657
    ref_dp_dT = 44.436693
    ref_saturated_water_rho = 999.789
    ref_saturated_steam_rho = 0.00485426
    ref_alpha = (- 11.529101)
    ref_saturated_water_enthalpy = 0.611786
    ref_saturated_steam_enthalpy = 2500500.0
    ref_phi = (- 0.04)
    ref_saturated_water_entropy = 0
    ref_saturated_steam_entropy = 9154.0
    self.assertAlmostEqual(ref_vapor_pressure, water.vaporPressure(Tk=Tk), 3)
    self.assertAlmostEqual(ref_vapor_pressure, steam.vaporPressure(Tk=Tk), 3)
    self.assertAlmostEqual(ref_dp_dT, water.vaporPressurePrime(Tk=Tk), 3)
    self.assertAlmostEqual(ref_dp_dT, steam.vaporPressurePrime(Tk=Tk), 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), 3)
    self.assertAlmostEqual(ref_alpha, steam.auxiliaryQuantitySpecificEnthalpy(Tk=Tk), 3)
    self.assertAlmostEqual(ref_saturated_water_enthalpy, water.enthalpy(Tk=Tk), 2)
    self.assertAlmostEqual((ref_saturated_steam_enthalpy / steam.enthalpy(Tk=Tk)), 1, 2)
    self.assertAlmostEqual(ref_phi, water.auxiliaryQuantitySpecificEntropy(Tk=Tk), 2)
    self.assertAlmostEqual(ref_phi, steam.auxiliaryQuantitySpecificEntropy(Tk=Tk), 2)
    self.assertAlmostEqual(ref_saturated_water_entropy, water.entropy(Tk=Tk), 3)
    self.assertAlmostEqual((ref_saturated_steam_entropy / steam.entropy(Tk=Tk)), 1, 3)