def test_mixed_layer():
"""Test the mixed layer calculation."""
pressure = np.array([959., 779.2, 751.3, 724.3, 700., 269.]) * units.hPa
temperature = np.array([22.2, 14.6, 12., 9.4, 7., -38.]) * units.degC
mixed_layer_temperature = mixed_layer(pressure,
temperature,
depth=250 * units.hPa)[0]
assert_almost_equal(mixed_layer_temperature, 16.4024930 * units.degC, 6)
def test_mixed_layer():
"""Test the mixed layer calculation."""
pressure = np.array([959., 779.2, 751.3, 724.3, 700., 269.]) * units.hPa
temperature = np.array([22.2, 14.6, 12., 9.4, 7., -38.]) * units.degC
mixed_layer_temperature = mixed_layer(pressure, temperature, depth=250 * units.hPa)[0]
assert_almost_equal(mixed_layer_temperature, 16.4024930 * units.degC, 6)