Code Duplication    Length = 10-12 lines in 2 locations

metpy/calc/tests/test_indices.py 2 locations

@@ 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)