def test_critical_angle_units():
"""Test critical angle with observed sounding and different storm motion units."""
data = get_upper_air_data(datetime(2016, 5, 22, 0), 'DDC')
# Set storm motion in m/s
ca_ms = critical_angle(data['pressure'], data['u_wind'],
data['v_wind'], data['height'],
stormu=10 * units('m/s'), stormv=10 * units('m/s'))
# Set same storm motion in kt and m/s
ca_kt_ms = critical_angle(data['pressure'], data['u_wind'],
data['v_wind'], data['height'],
stormu=10 * units('m/s'), stormv=19.4384449244 * units('kt'))
# Make sure the resulting critical angles are equal
assert_almost_equal(ca_ms, ca_kt_ms, 8)
def test_critical_angle_units():
"""Test critical angle with observed sounding and different storm motion units."""
data = get_upper_air_data(datetime(2016, 5, 22, 0), 'DDC')
# Set storm motion in m/s
ca_ms = critical_angle(data['pressure'], data['u_wind'],
data['v_wind'], data['height'],
u_storm=10 * units('m/s'), v_storm=10 * units('m/s'))
# Set same storm motion in kt and m/s
ca_kt_ms = critical_angle(data['pressure'], data['u_wind'],
data['v_wind'], data['height'],
u_storm=10 * units('m/s'), v_storm=19.4384449244 * units('kt'))
# Make sure the resulting critical angles are equal
assert_almost_equal(ca_ms, ca_kt_ms, 8)
def test_critical_angle():
"""Test critical angle with observed sounding."""
data = get_upper_air_data(datetime(2016, 5, 22, 0), 'DDC')
ca = critical_angle(data['pressure'], data['u_wind'],
data['v_wind'], data['height'],
stormu=0 * units('m/s'), stormv=0 * units('m/s'))
truth = [140.0626637513269] * units('degrees')
assert_almost_equal(ca, truth, 8)
def test_critical_angle():
"""Test critical angle with observed sounding."""
data = get_upper_air_data(datetime(2016, 5, 22, 0), 'DDC')
ca = critical_angle(data['pressure'], data['u_wind'],
data['v_wind'], data['height'],
u_storm=0 * units('m/s'), v_storm=0 * units('m/s'))
truth = [140.0626637513269] * units('degrees')
assert_almost_equal(ca, truth, 8)