def test_mean_pressure_weighted_temperature():
"""Test pressure-weighted mean temperature function with vertical interpolation."""
data = get_upper_air_data(datetime(2016, 5, 22, 0), 'DDC')
t, = mean_pressure_weighted(data['pressure'],
data['temperature'],
height=data['height'],
depth=6000 * units('meter'))
assert_almost_equal(t, 281.535035296836 * units('kelvin'), 7)
def test_mean_pressure_weighted():
"""Test pressure-weighted mean wind function with vertical interpolation."""
data = get_upper_air_data(datetime(2016, 5, 22, 0), 'DDC')
u, v = mean_pressure_weighted(data['pressure'],
data['u_wind'],
data['v_wind'],
height=data['height'],
depth=6000 * units('meter'))
assert_almost_equal(u, 6.0208700094534775 * units('m/s'), 7)
assert_almost_equal(v, 7.966031839967931 * units('m/s'), 7)
def test_mean_pressure_weighted():
"""Test pressure-weighted mean wind function with vertical interpolation."""
data = get_upper_air_data(datetime(2016, 5, 22, 0), 'DDC')
u, v = mean_pressure_weighted(data['pressure'],
data['u_wind'],
data['v_wind'],
heights=data['height'],
depth=6000 * units('meter'))
assert_almost_equal(u, 6.0208700094534775 * units('m/s'), 7)
assert_almost_equal(v, 7.966031839967931 * units('m/s'), 7)
def test_mean_pressure_weighted_elevated():
"""Test pressure-weighted mean wind function with a base above the surface."""
data = get_upper_air_data(datetime(2016, 5, 22, 0), 'DDC')
u, v = mean_pressure_weighted(data['pressure'],
data['u_wind'],
data['v_wind'],
heights=data['height'],
depth=3000 * units('meter'),
bottom=data['height'][0] + 3000 * units('meter'))
assert_almost_equal(u, 8.270829843626476 * units('m/s'), 7)
assert_almost_equal(v, 1.7392601775853547 * units('m/s'), 7)
def test_mean_pressure_weighted_elevated():
"""Test pressure-weighted mean wind function with a base above the surface."""
data = get_upper_air_data(datetime(2016, 5, 22, 0), 'DDC')
u, v = mean_pressure_weighted(data['pressure'],
data['u_wind'],
data['v_wind'],
height=data['height'],
depth=3000 * units('meter'),
bottom=data['height'][0] + 3000 * units('meter'))
assert_almost_equal(u, 8.270829843626476 * units('m/s'), 7)
assert_almost_equal(v, 1.7392601775853547 * units('m/s'), 7)
