|
@@ 38-49 (lines=12) @@
|
| 35 |
|
data.variables['v_wind'][:], |
| 36 |
|
heights=data.variables['height'][:], |
| 37 |
|
depth=6000 * units('meter')) |
| 38 |
|
assert_almost_equal(u, 6.0208700094534775 * units('m/s'), 7) |
| 39 |
|
assert_almost_equal(v, 7.966031839967931 * units('m/s'), 7) |
| 40 |
|
|
| 41 |
|
|
| 42 |
|
def test_mean_pressure_weighted_elevated(): |
| 43 |
|
"""Test pressure-weighted mean wind function with a base above the surface.""" |
| 44 |
|
with UseSampleData(): |
| 45 |
|
data = get_upper_air_data(datetime(2016, 5, 22, 0), 'DDC', source='wyoming') |
| 46 |
|
u, v = mean_pressure_weighted(data.variables['pressure'][:], |
| 47 |
|
data.variables['u_wind'][:], |
| 48 |
|
data.variables['v_wind'][:], |
| 49 |
|
heights=data.variables['height'][:], |
| 50 |
|
depth=3000 * units('meter'), |
| 51 |
|
bottom=data.variables['height'][0] + 3000 * units('meter')) |
| 52 |
|
assert_almost_equal(u, 8.270829843626476 * units('m/s'), 7) |
|
@@ 64-73 (lines=10) @@
|
| 61 |
|
data.variables['u_wind'][:], data.variables['v_wind'][:], |
| 62 |
|
data.variables['height'][:])) |
| 63 |
|
truth = [1.4537892577864744, 2.0169333025630616, 10.587950761120482, 13.915130377372801, |
| 64 |
|
6.0208700094534775, 7.9660318399679308] * units('m/s') |
| 65 |
|
assert_almost_equal(motion.flatten(), truth, 8) |
| 66 |
|
|
| 67 |
|
|
| 68 |
|
def test_bulk_shear(): |
| 69 |
|
"""Test bulk shear with observed sounding.""" |
| 70 |
|
with UseSampleData(): |
| 71 |
|
data = get_upper_air_data(datetime(2016, 5, 22, 0), 'DDC', source='wyoming') |
| 72 |
|
u, v = bulk_shear(data.variables['pressure'][:], data.variables['u_wind'][:], |
| 73 |
|
data.variables['v_wind'][:], heights=data.variables['height'][:], |
| 74 |
|
depth=6000 * units('meter')) |
| 75 |
|
truth = [29.899581266946115, -14.389225800205509] * units('knots') |
| 76 |
|
assert_almost_equal(u.to('knots'), truth[0], 8) |
Unidata /
MetPy
