Skip to content

Instantly share code, notes, and snippets.

@jthielen

jthielen/metpy_test_with_pint_master_20190726.txt

Created July 27, 2019 02:28
Show Gist options
  • Save jthielen/f24c9a953ff4edc3a795179638a77f5e to your computer and use it in GitHub Desktop.
Save jthielen/f24c9a953ff4edc3a795179638a77f5e to your computer and use it in GitHub Desktop.
Download ZIP
Results of a run of pytest with the current master branches of MetPy and pint (as of 2019-07-26)
Raw
metpy_test_with_pint_master_20190726.txt
============================= test session starts ==============================
platform linux -- Python 3.7.3, pytest-4.6.2, py-1.8.0, pluggy-0.12.0
Matplotlib: 3.0.2
Freetype: 2.6.1
Dependencies: Matplotlib (3.0.2), NumPy (1.16.4), Pint (0.10.dev0), SciPy (1.3.0)
rootdir: /home/jon/dev/MetPy, inifile: setup.cfg
plugins: flake8-1.0.4, cov-2.7.1, mpl-0.10
collected 878 items
metpy/calc/tests/test_basic.py ............................FF.F.....F... [ 4%]
...............F.F. [ 6%]
metpy/calc/tests/test_calc_tools.py ...............................F.F.. [ 10%]
...F............F................................................. [ 18%]
metpy/calc/tests/test_cross_sections.py ............ [ 19%]
metpy/calc/tests/test_indices.py ............... [ 21%]
metpy/calc/tests/test_kinematics.py .................................... [ 25%]
.....FF........F [ 27%]
metpy/calc/tests/test_thermo.py ........................................ [ 32%]
......................FF..............FF................F.....FFFF...... [ 40%]
[ 40%]
metpy/calc/tests/test_turbulence.py ........................... [ 43%]
metpy/interpolate/tests/test_geometry.py ............. [ 44%]
metpy/interpolate/tests/test_grid.py ................................... [ 48%]
...... [ 49%]
metpy/interpolate/tests/test_interpolate_tools.py ...... [ 50%]
metpy/interpolate/tests/test_one_dimension.py ................. [ 52%]
metpy/interpolate/tests/test_points.py ............... [ 53%]
metpy/interpolate/tests/test_slices.py ....... [ 54%]
metpy/io/tests/test_gini.py .......... [ 55%]
metpy/io/tests/test_nexrad.py .......................................... [ 60%]
........................................................................ [ 68%]
................ [ 70%]
metpy/plots/tests/test_cartopy_utils.py .... [ 70%]
metpy/plots/tests/test_ctables.py .......... [ 72%]
metpy/plots/tests/test_declarative.py .............. [ 73%]
metpy/plots/tests/test_mapping.py ................. [ 75%]
metpy/plots/tests/test_mpl.py . [ 75%]
metpy/plots/tests/test_skewt.py ........................x [ 78%]
metpy/plots/tests/test_station_plot.py .................. [ 80%]
metpy/plots/tests/test_util.py ............. [ 82%]
metpy/plots/tests/test_wx_symbols.py ... [ 82%]
metpy/tests/test_cbook.py . [ 82%]
metpy/tests/test_deprecation.py . [ 82%]
metpy/tests/test_gridding.py . [ 82%]
metpy/tests/test_units.py ..FF...................... [ 85%]
metpy/tests/test_xarray.py ......F...................................... [ 90%]
........................................................................ [ 99%]
........ [100%]
=================================== FAILURES ===================================
_________________________ test_height_to_geopotential __________________________
def test_height_to_geopotential():
"""Test conversion from height to geopotential."""
height = units.Quantity([0, 1000, 2000, 3000], units.m)
geopot = height_to_geopotential(height)
assert_array_almost_equal(geopot, units.Quantity([0., 9817, 19632,
> 29443], units('m**2 / second**2')), 0)
metpy/calc/tests/test_basic.py:293:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
actual = array([ 0. , 9818.14474002, 19633.20840782, 29445.19245349])
desired = array([ 0., 9817., 19632., 29443.]), decimal = 0
def assert_array_almost_equal(actual, desired, decimal=7):
"""Check that arrays are almost equal, including units.
Wrapper around :func:`numpy.testing.assert_array_almost_equal`
"""
actual, desired = check_and_drop_units(actual, desired)
> numpy.testing.assert_array_almost_equal(actual, desired, decimal)
E AssertionError:
E Arrays are not almost equal to 0 decimals
E
E Mismatch: 25%
E Max absolute difference: 2.19245349
E Max relative difference: nan
E x: array([ 0., 9818., 19633., 29445.])
E y: array([ 0., 9817., 19632., 29443.])
metpy/testing.py:160: AssertionError
______________________ test_height_to_geopotential_32bit _______________________
def test_height_to_geopotential_32bit():
"""Test conversion to geopotential with 32-bit values."""
heights = np.linspace(20597, 20598, 11, dtype=np.float32) * units.m
truth = np.array([201590.422, 201591.391, 201592.375, 201593.344,
201594.312, 201595.297, 201596.266, 201597.250,
201598.219, 201599.203, 201600.172], dtype=np.float32) * units('J/kg')
> assert_almost_equal(height_to_geopotential(heights), truth, 2)
metpy/calc/tests/test_basic.py:303:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
actual = array([201604.31, 201605.28, 201606.27, 201607.23, 201608.22, 201609.19,
201610.17, 201611.14, 201612.12, 201613.1 , 201614.06],
dtype=float32)
desired = array([201590.42, 201591.39, 201592.38, 201593.34, 201594.31, 201595.3 ,
201596.27, 201597.25, 201598.22, 201599.2 , 201600.17],
dtype=float32)
decimal = 2
def assert_almost_equal(actual, desired, decimal=7):
"""Check that values are almost equal, including units.
Wrapper around :func:`numpy.testing.assert_almost_equal`
"""
actual, desired = check_and_drop_units(actual, desired)
> numpy.testing.assert_almost_equal(actual, desired, decimal)
E AssertionError:
E Arrays are not almost equal to 2 decimals
E
E Mismatch: 100%
E Max absolute difference: 13.90625
E Max relative difference: 6.898136e-05
E x: array([201604.31, 201605.28, 201606.27, 201607.23, 201608.22, 201609.19,
E 201610.17, 201611.14, 201612.12, 201613.1 , 201614.06],
E dtype=float32)
E y: array([201590.42, 201591.39, 201592.38, 201593.34, 201594.31, 201595.3 ,
E 201596.27, 201597.25, 201598.22, 201599.2 , 201600.17],
E dtype=float32)
metpy/testing.py:151: AssertionError
______________________ test_geopotential_to_height_32bit _______________________
def test_geopotential_to_height_32bit():
"""Test conversion from geopotential to height with 32-bit values."""
geopot = np.arange(201590, 201600, dtype=np.float32) * units('J/kg')
truth = np.array([20596.957, 20597.059, 20597.162, 20597.266,
20597.365, 20597.469, 20597.570, 20597.674,
20597.777, 20597.881], dtype=np.float32) * units.m
> assert_almost_equal(geopotential_to_height(geopot), truth, 2)
metpy/calc/tests/test_basic.py:321:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
actual = array([20595.531, 20595.635, 20595.736, 20595.84 , 20595.941, 20596.045,
20596.146, 20596.25 , 20596.354, 20596.455], dtype=float32)
desired = array([20596.957, 20597.059, 20597.162, 20597.266, 20597.365, 20597.469,
20597.57 , 20597.674, 20597.777, 20597.88 ], dtype=float32)
decimal = 2
def assert_almost_equal(actual, desired, decimal=7):
"""Check that values are almost equal, including units.
Wrapper around :func:`numpy.testing.assert_almost_equal`
"""
actual, desired = check_and_drop_units(actual, desired)
> numpy.testing.assert_almost_equal(actual, desired, decimal)
E AssertionError:
E Arrays are not almost equal to 2 decimals
E
E Mismatch: 100%
E Max absolute difference: 1.4257812
E Max relative difference: 6.922291e-05
E x: array([20595.53, 20595.63, 20595.74, 20595.84, 20595.94, 20596.04,
E 20596.15, 20596.25, 20596.35, 20596.46], dtype=float32)
E y: array([20596.96, 20597.06, 20597.16, 20597.27, 20597.37, 20597.47,
E 20597.57, 20597.67, 20597.78, 20597.88], dtype=float32)
metpy/testing.py:151: AssertionError
_________________________ test_add_pressure_to_height __________________________
def test_add_pressure_to_height():
"""Test the height at pressure above height calculation."""
height = add_pressure_to_height(110.8286757 * units.m, 100 * units.hPa)
> assert_almost_equal(height, 988.0028867 * units.meter, 5)
metpy/calc/tests/test_basic.py:398:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
actual = 988.0029400097312, desired = 988.0028867, decimal = 5
def assert_almost_equal(actual, desired, decimal=7):
"""Check that values are almost equal, including units.
Wrapper around :func:`numpy.testing.assert_almost_equal`
"""
actual, desired = check_and_drop_units(actual, desired)
> numpy.testing.assert_almost_equal(actual, desired, decimal)
E AssertionError:
E Arrays are not almost equal to 5 decimals
E ACTUAL: 988.0029400097312
E DESIRED: 988.0028867
metpy/testing.py:151: AssertionError
___________________ test_altimeter_to_station_pressure_inhg ____________________
def test_altimeter_to_station_pressure_inhg():
"""Test the altimeter to station pressure function with inches of mercury."""
altim = 29.8 * units.inHg
elev = 500 * units.m
> res = altimeter_to_station_pressure(altim, elev)
metpy/calc/tests/test_basic.py:637:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
metpy/xarray.py:457: in wrapper
return func(*args, **kwargs)
metpy/calc/basic.py:919: in altimeter_to_station_pressure
return ((altimeter_value ** n - ((p0 ** n * gamma * height) / t0)) ** (1 / n) + (
../pint/pint/quantity.py:830: in __sub__
return self._add_sub(other, operator.sub)
../pint/pint/quantity.py:76: in wrapped
result = f(self, *args, **kwargs)
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
self = <Quantity(1.9076292329237703, 'inch_Hg ** 0.190267')>
other = <Quantity(0.04210930305337787, 'hectopascal ** 0.190267')>
op = <built-in function sub>
@check_implemented
def _add_sub(self, other, op):
"""Perform addition or subtraction operation and return the result.
:param other: object to be added to / subtracted from self
:type other: Quantity or any type accepted by :func:`_to_magnitude`
:param op: operator function (e.g. operator.add, operator.isub)
:type op: function
"""
if not self._check(other):
# other not from same Registry or not a Quantity
if _eq(other, 0, True):
# If the other value is 0 (but not Quantity 0)
# do the operation without checking units.
# We do the calculation instead of just returning the same
# value to enforce any shape checking and type casting due to
# the operation.
units = self._units
magnitude = op(self._magnitude,
_to_magnitude(other, self.force_ndarray))
elif self.dimensionless:
units = UnitsContainer()
magnitude = op(self.to(units)._magnitude,
_to_magnitude(other, self.force_ndarray))
else:
raise DimensionalityError(self._units, 'dimensionless')
return self.__class__(magnitude, units)
if not self.dimensionality == other.dimensionality:
raise DimensionalityError(self._units, other._units,
self.dimensionality,
> other.dimensionality)
E pint.errors.DimensionalityError: Cannot convert from 'inch_Hg ** 0.190267' ([mass] ** 0.190267 / [length] ** 0.190267 / [time] ** 0.380533) to 'hectopascal ** 0.190267' ([mass] ** 0.190267 / [length] ** 0.190267 / [time] ** 0.380533)
../pint/pint/quantity.py:736: DimensionalityError
__________________ test_altimiter_to_sea_level_pressure_inhg ___________________
def test_altimiter_to_sea_level_pressure_inhg():
"""Test the altimeter to sea level pressure function with inches of mercury."""
altim = 29.8 * units.inHg
elev = 500 * units.m
temp = 30 * units.degC
> res = altimeter_to_sea_level_pressure(altim, elev, temp)
metpy/calc/tests/test_basic.py:656:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
metpy/xarray.py:457: in wrapper
return func(*args, **kwargs)
metpy/calc/basic.py:985: in altimeter_to_sea_level_pressure
psfc = altimeter_to_station_pressure(altimeter_value, height)
metpy/xarray.py:457: in wrapper
return func(*args, **kwargs)
metpy/calc/basic.py:919: in altimeter_to_station_pressure
return ((altimeter_value ** n - ((p0 ** n * gamma * height) / t0)) ** (1 / n) + (
../pint/pint/quantity.py:830: in __sub__
return self._add_sub(other, operator.sub)
../pint/pint/quantity.py:76: in wrapped
result = f(self, *args, **kwargs)
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
self = <Quantity(1.9076292329237703, 'inch_Hg ** 0.190267')>
other = <Quantity(0.04210930305337787, 'hectopascal ** 0.190267')>
op = <built-in function sub>
@check_implemented
def _add_sub(self, other, op):
"""Perform addition or subtraction operation and return the result.
:param other: object to be added to / subtracted from self
:type other: Quantity or any type accepted by :func:`_to_magnitude`
:param op: operator function (e.g. operator.add, operator.isub)
:type op: function
"""
if not self._check(other):
# other not from same Registry or not a Quantity
if _eq(other, 0, True):
# If the other value is 0 (but not Quantity 0)
# do the operation without checking units.
# We do the calculation instead of just returning the same
# value to enforce any shape checking and type casting due to
# the operation.
units = self._units
magnitude = op(self._magnitude,
_to_magnitude(other, self.force_ndarray))
elif self.dimensionless:
units = UnitsContainer()
magnitude = op(self.to(units)._magnitude,
_to_magnitude(other, self.force_ndarray))
else:
raise DimensionalityError(self._units, 'dimensionless')
return self.__class__(magnitude, units)
if not self.dimensionality == other.dimensionality:
raise DimensionalityError(self._units, other._units,
self.dimensionality,
> other.dimensionality)
E pint.errors.DimensionalityError: Cannot convert from 'inch_Hg ** 0.190267' ([mass] ** 0.190267 / [length] ** 0.190267 / [time] ** 0.380533) to 'hectopascal ** 0.190267' ([mass] ** 0.190267 / [length] ** 0.190267 / [time] ** 0.380533)
../pint/pint/quantity.py:736: DimensionalityError
_____ test_get_bound_pressure_height[pressure4-bound4-None-True-expected4] _____
pressure = <Quantity([1000. 900. 800. 700. 600. 500. 400. 300. 200. 100.], 'hectopascal')>
bound = <Quantity(0.9880028, 'kilometer')>, hgts = None, interp = True
expected = (<Quantity(900, 'hectopascal')>, <Quantity(0.9880028, 'kilometer')>)
@pytest.mark.parametrize('pressure, bound, hgts, interp, expected', [
(get_bounds_data()[0], 900 * units.hPa, None, True,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 900 * units.hPa, None, False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 870 * units.hPa, None, True,
(870 * units.hPa, 1.2665298 * units.kilometer)),
(get_bounds_data()[0], 870 * units.hPa, None, False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 0.9880028 * units.kilometer, None, True,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 0.9880028 * units.kilometer, None, False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 1.2665298 * units.kilometer, None, True,
(870 * units.hPa, 1.2665298 * units.kilometer)),
(get_bounds_data()[0], 1.2665298 * units.kilometer, None, False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 900 * units.hPa, get_bounds_data()[1], True,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 900 * units.hPa, get_bounds_data()[1], False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 870 * units.hPa, get_bounds_data()[1], True,
(870 * units.hPa, 1.2643214 * units.kilometer)),
(get_bounds_data()[0], 870 * units.hPa, get_bounds_data()[1], False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 0.9880028 * units.kilometer, get_bounds_data()[1], True,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 0.9880028 * units.kilometer, get_bounds_data()[1], False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 1.2665298 * units.kilometer, get_bounds_data()[1], True,
(870.9869087 * units.hPa, 1.2665298 * units.kilometer)),
(get_bounds_data()[0], 1.2665298 * units.kilometer, get_bounds_data()[1], False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 0.98800289 * units.kilometer, get_bounds_data()[1], True,
(900 * units.hPa, 0.9880028 * units.kilometer))
])
def test_get_bound_pressure_height(pressure, bound, hgts, interp, expected):
"""Test getting bounds in layers with various parameter combinations."""
bounds = _get_bound_pressure_height(pressure, bound, heights=hgts, interpolate=interp)
> assert_array_almost_equal(bounds[0], expected[0], 5)
metpy/calc/tests/test_calc_tools.py:261:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
actual = array(900.00001612), desired = 900, decimal = 5
def assert_array_almost_equal(actual, desired, decimal=7):
"""Check that arrays are almost equal, including units.
Wrapper around :func:`numpy.testing.assert_array_almost_equal`
"""
actual, desired = check_and_drop_units(actual, desired)
> numpy.testing.assert_array_almost_equal(actual, desired, decimal)
E AssertionError:
E Arrays are not almost equal to 5 decimals
E
E Mismatch: 100%
E Max absolute difference: 1.61150167e-05
E Max relative difference: 1.79055741e-08
E x: array(900.00002)
E y: array(900)
metpy/testing.py:160: AssertionError
_____ test_get_bound_pressure_height[pressure6-bound6-None-True-expected6] _____
pressure = <Quantity([1000. 900. 800. 700. 600. 500. 400. 300. 200. 100.], 'hectopascal')>
bound = <Quantity(1.2665298, 'kilometer')>, hgts = None, interp = True
expected = (<Quantity(870, 'hectopascal')>, <Quantity(1.2665298, 'kilometer')>)
@pytest.mark.parametrize('pressure, bound, hgts, interp, expected', [
(get_bounds_data()[0], 900 * units.hPa, None, True,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 900 * units.hPa, None, False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 870 * units.hPa, None, True,
(870 * units.hPa, 1.2665298 * units.kilometer)),
(get_bounds_data()[0], 870 * units.hPa, None, False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 0.9880028 * units.kilometer, None, True,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 0.9880028 * units.kilometer, None, False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 1.2665298 * units.kilometer, None, True,
(870 * units.hPa, 1.2665298 * units.kilometer)),
(get_bounds_data()[0], 1.2665298 * units.kilometer, None, False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 900 * units.hPa, get_bounds_data()[1], True,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 900 * units.hPa, get_bounds_data()[1], False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 870 * units.hPa, get_bounds_data()[1], True,
(870 * units.hPa, 1.2643214 * units.kilometer)),
(get_bounds_data()[0], 870 * units.hPa, get_bounds_data()[1], False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 0.9880028 * units.kilometer, get_bounds_data()[1], True,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 0.9880028 * units.kilometer, get_bounds_data()[1], False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 1.2665298 * units.kilometer, get_bounds_data()[1], True,
(870.9869087 * units.hPa, 1.2665298 * units.kilometer)),
(get_bounds_data()[0], 1.2665298 * units.kilometer, get_bounds_data()[1], False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 0.98800289 * units.kilometer, get_bounds_data()[1], True,
(900 * units.hPa, 0.9880028 * units.kilometer))
])
def test_get_bound_pressure_height(pressure, bound, hgts, interp, expected):
"""Test getting bounds in layers with various parameter combinations."""
bounds = _get_bound_pressure_height(pressure, bound, heights=hgts, interpolate=interp)
> assert_array_almost_equal(bounds[0], expected[0], 5)
metpy/calc/tests/test_calc_tools.py:261:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
actual = array(870.00001703), desired = 870, decimal = 5
def assert_array_almost_equal(actual, desired, decimal=7):
"""Check that arrays are almost equal, including units.
Wrapper around :func:`numpy.testing.assert_array_almost_equal`
"""
actual, desired = check_and_drop_units(actual, desired)
> numpy.testing.assert_array_almost_equal(actual, desired, decimal)
E AssertionError:
E Arrays are not almost equal to 5 decimals
E
E Mismatch: 100%
E Max absolute difference: 1.70334024e-05
E Max relative difference: 1.95786234e-08
E x: array(870.00002)
E y: array(870)
metpy/testing.py:160: AssertionError
__ test_get_bound_pressure_height[pressure12-bound12-hgts12-True-expected12] ___
pressure = <Quantity([1000. 900. 800. 700. 600. 500. 400. 300. 200. 100.], 'hectopascal')>
bound = <Quantity(0.9880028, 'kilometer')>
hgts = <Quantity([ 0.11082868 0.98800295 1.94800727 3.01066437 4.20430412 5.57162493
7.18180872 9.15932612 11.76894158 15.78930577], 'kilometer')>
interp = True
expected = (<Quantity(900, 'hectopascal')>, <Quantity(0.9880028, 'kilometer')>)
@pytest.mark.parametrize('pressure, bound, hgts, interp, expected', [
(get_bounds_data()[0], 900 * units.hPa, None, True,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 900 * units.hPa, None, False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 870 * units.hPa, None, True,
(870 * units.hPa, 1.2665298 * units.kilometer)),
(get_bounds_data()[0], 870 * units.hPa, None, False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 0.9880028 * units.kilometer, None, True,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 0.9880028 * units.kilometer, None, False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 1.2665298 * units.kilometer, None, True,
(870 * units.hPa, 1.2665298 * units.kilometer)),
(get_bounds_data()[0], 1.2665298 * units.kilometer, None, False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 900 * units.hPa, get_bounds_data()[1], True,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 900 * units.hPa, get_bounds_data()[1], False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 870 * units.hPa, get_bounds_data()[1], True,
(870 * units.hPa, 1.2643214 * units.kilometer)),
(get_bounds_data()[0], 870 * units.hPa, get_bounds_data()[1], False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 0.9880028 * units.kilometer, get_bounds_data()[1], True,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 0.9880028 * units.kilometer, get_bounds_data()[1], False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 1.2665298 * units.kilometer, get_bounds_data()[1], True,
(870.9869087 * units.hPa, 1.2665298 * units.kilometer)),
(get_bounds_data()[0], 1.2665298 * units.kilometer, get_bounds_data()[1], False,
(900 * units.hPa, 0.9880028 * units.kilometer)),
(get_bounds_data()[0], 0.98800289 * units.kilometer, get_bounds_data()[1], True,
(900 * units.hPa, 0.9880028 * units.kilometer))
])
def test_get_bound_pressure_height(pressure, bound, hgts, interp, expected):
"""Test getting bounds in layers with various parameter combinations."""
bounds = _get_bound_pressure_height(pressure, bound, heights=hgts, interpolate=interp)
> assert_array_almost_equal(bounds[0], expected[0], 5)
metpy/calc/tests/test_calc_tools.py:261:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
actual = array([900.00001682]), desired = 900, decimal = 5
def assert_array_almost_equal(actual, desired, decimal=7):
"""Check that arrays are almost equal, including units.
Wrapper around :func:`numpy.testing.assert_array_almost_equal`
"""
actual, desired = check_and_drop_units(actual, desired)
> numpy.testing.assert_array_almost_equal(actual, desired, decimal)
E AssertionError:
E Arrays are not almost equal to 5 decimals
E
E Mismatch: 100%
E Max absolute difference: 1.68248043e-05
E Max relative difference: 1.8694227e-08
E x: array([900.00002])
E y: array(900)
metpy/testing.py:160: AssertionError
____ test_get_layer[pressure2-variable2-None-bottom2-depth2-True-expected2] ____
pressure = <Quantity([1000 900 800 700 600 500 400 300 200 100], 'hectopascal')>
variable = <Quantity([ 25. 16.66666667 8.33333333 0. -8.33333333
-16.66666667 -25. -33.33333333 -41.66666667 -50. ], 'degree_Celsius')>
heights = None, bottom = <Quantity(2, 'kilometer')>
depth = <Quantity(3, 'kilometer')>, interp = True
expected = (<Quantity([794.85264282 700. 600. 540.01696548], 'hectopascal')>, <Quantity([ 7.93049516 0. -8.33333333 -13.14758845], 'degree_Celsius')>)
@pytest.mark.parametrize('pressure, variable, heights, bottom, depth, interp, expected', [
(layer_test_data()[0], layer_test_data()[1], None, None, 150 * units.hPa, True,
(np.array([1000, 900, 850]) * units.hPa,
np.array([25.0, 16.666666, 12.62262]) * units.degC)),
(layer_test_data()[0], layer_test_data()[1], None, None, 150 * units.hPa, False,
(np.array([1000, 900]) * units.hPa, np.array([25.0, 16.666666]) * units.degC)),
(layer_test_data()[0], layer_test_data()[1], None, 2 * units.km, 3 * units.km, True,
(np.array([794.85264282, 700., 600., 540.01696548]) * units.hPa,
np.array([7.93049516, 0., -8.33333333, -13.14758845]) * units.degC))
])
def test_get_layer(pressure, variable, heights, bottom, depth, interp, expected):
"""Test get_layer functionality."""
p_layer, y_layer = get_layer(pressure, variable, heights=heights, bottom=bottom,
depth=depth, interpolate=interp)
> assert_array_almost_equal(p_layer, expected[0], 5)
metpy/calc/tests/test_calc_tools.py:349:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
actual = array([794.85265485, 700. , 600. , 540.01698666])
desired = array([794.85264282, 700. , 600. , 540.01696548])
decimal = 5
def assert_array_almost_equal(actual, desired, decimal=7):
"""Check that arrays are almost equal, including units.
Wrapper around :func:`numpy.testing.assert_array_almost_equal`
"""
actual, desired = check_and_drop_units(actual, desired)
> numpy.testing.assert_array_almost_equal(actual, desired, decimal)
E AssertionError:
E Arrays are not almost equal to 5 decimals
E
E Mismatch: 25%
E Max absolute difference: 2.11808643e-05
E Max relative difference: 3.92225906e-08
E x: array([794.85265, 700. , 600. , 540.01699])
E y: array([794.85264, 700. , 600. , 540.01697])
metpy/testing.py:160: AssertionError
____________________ test_q_vector_without_static_stability ____________________
q_vector_data = (<Quantity([[ 2.50000000e+01 8.68240888e+00 -8.68240888e+00 -2.50000000e+01]
[ 8.68240888e+00 -6.12323400e-15 -6.123...047447 185324.79047447 185324.79047447]
[185324.79047447 185324.79047447 185324.79047447 185324.79047447]], 'meter')>)
def test_q_vector_without_static_stability(q_vector_data):
"""Test the Q-vector function without using static stability."""
u, v, temp, p, dx, dy = q_vector_data
# Treating as 700 hPa data
q1, q2 = q_vector(u, v, temp[1], p[1], dx, dy)
q1_truth = (np.array([[-2.7454089e-14, -3.0194267e-13, -3.0194267e-13, -2.7454089e-14],
[-1.8952185e-13, -2.2269905e-14, -2.2269905e-14, -1.8952185e-13],
[-1.9918390e-13, -2.3370829e-14, -2.3370829e-14, -1.9918390e-13],
[-5.6160772e-14, -3.5145951e-13, -3.5145951e-13, -5.6160772e-14]])
* units('m^2 kg^-1 s^-1'))
q2_truth = (np.array([[-4.4976059e-14, -4.3582378e-13, 4.3582378e-13, 4.4976059e-14],
[-3.0124244e-13, -3.5724617e-14, 3.5724617e-14, 3.0124244e-13],
[3.1216232e-13, 3.6662900e-14, -3.6662900e-14, -3.1216232e-13],
[8.6038280e-14, 4.6968342e-13, -4.6968342e-13, -8.6038280e-14]])
* units('m^2 kg^-1 s^-1'))
> assert_almost_equal(q1, q1_truth, 20)
metpy/calc/tests/test_kinematics.py:1062:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
actual = array([[-2.74540907e-14, -3.01942689e-13, -3.01942689e-13,
-2.74540907e-14],
[-1.89521858e-13, -2.22699...07e-14,
-1.99183913e-13],
[-5.61607757e-14, -3.51459530e-13, -3.51459530e-13,
-5.61607757e-14]])
desired = array([[-2.7454089e-14, -3.0194267e-13, -3.0194267e-13, -2.7454089e-14],
[-1.8952185e-13, -2.2269905e-14, -2.22...3370829e-14, -2.3370829e-14, -1.9918390e-13],
[-5.6160772e-14, -3.5145951e-13, -3.5145951e-13, -5.6160772e-14]])
decimal = 20
def assert_almost_equal(actual, desired, decimal=7):
"""Check that values are almost equal, including units.
Wrapper around :func:`numpy.testing.assert_almost_equal`
"""
actual, desired = check_and_drop_units(actual, desired)
> numpy.testing.assert_almost_equal(actual, desired, decimal)
E AssertionError:
E Arrays are not almost equal to 20 decimals
E
E Mismatch: 25%
E Max absolute difference: 1.97713253e-20
E Max relative difference: 7.30728665e-08
E x: array([[-2.7454090723180448e-14, -3.0194268852999482e-13,
E -3.0194268852998139e-13, -2.7454090723230513e-14],
E [-1.8952185792609708e-13, -2.2269906036398466e-14,...
E y: array([[-2.7454089e-14, -3.0194267e-13, -3.0194267e-13, -2.7454089e-14],
E [-1.8952185e-13, -2.2269905e-14, -2.2269905e-14, -1.8952185e-13],
E [-1.9918390e-13, -2.3370829e-14, -2.3370829e-14, -1.9918390e-13],
E [-5.6160772e-14, -3.5145951e-13, -3.5145951e-13, -5.6160772e-14]])
metpy/testing.py:151: AssertionError
_____________________ test_q_vector_with_static_stability ______________________
q_vector_data = (<Quantity([[ 2.50000000e+01 8.68240888e+00 -8.68240888e+00 -2.50000000e+01]
[ 8.68240888e+00 -6.12323400e-15 -6.123...047447 185324.79047447 185324.79047447]
[185324.79047447 185324.79047447 185324.79047447 185324.79047447]], 'meter')>)
def test_q_vector_with_static_stability(q_vector_data):
"""Test the Q-vector function using static stability."""
u, v, temp, p, dx, dy = q_vector_data
sigma = static_stability(p[:, np.newaxis, np.newaxis], temp)
# Treating as 700 hPa data
q1, q2 = q_vector(u, v, temp[1], p[1], dx, dy, sigma[1])
q1_truth = (np.array([[-1.4158140e-08, -1.6197987e-07, -1.6875014e-07, -1.6010616e-08],
[-9.3971386e-08, -1.1252476e-08, -1.1252476e-08, -9.7617234e-08],
[-1.0785670e-07, -1.2403513e-08, -1.2403513e-08, -1.0364793e-07],
[-2.9186946e-08, -1.7577703e-07, -1.6937879e-07, -2.6112047e-08]])
* units('kg m^-2 s^-3'))
q2_truth = (np.array([[-2.3194263e-08, -2.3380160e-07, 2.4357380e-07, 2.6229040e-08],
[-1.4936626e-07, -1.8050836e-08, 1.8050836e-08, 1.5516129e-07],
[1.6903373e-07, 1.9457964e-08, -1.9457964e-08, -1.6243771e-07],
[4.4714390e-08, 2.3490489e-07, -2.2635441e-07, -4.0003646e-08]])
* units('kg m^-2 s^-3'))
> assert_almost_equal(q1, q1_truth, 14)
metpy/calc/tests/test_kinematics.py:1086:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
actual = array([[-1.41581383e-08, -1.61979842e-07, -1.68750114e-07,
-1.60106130e-08],
[-9.39713720e-08, -1.12524...11e-08,
-1.03647910e-07],
[-2.91869419e-08, -1.75777003e-07, -1.69378764e-07,
-2.61120434e-08]])
desired = array([[-1.4158140e-08, -1.6197987e-07, -1.6875014e-07, -1.6010616e-08],
[-9.3971386e-08, -1.1252476e-08, -1.12...2403513e-08, -1.2403513e-08, -1.0364793e-07],
[-2.9186946e-08, -1.7577703e-07, -1.6937879e-07, -2.6112047e-08]])
decimal = 14
def assert_almost_equal(actual, desired, decimal=7):
"""Check that values are almost equal, including units.
Wrapper around :func:`numpy.testing.assert_almost_equal`
"""
actual, desired = check_and_drop_units(actual, desired)
> numpy.testing.assert_almost_equal(actual, desired, decimal)
E AssertionError:
E Arrays are not almost equal to 14 decimals
E
E Mismatch: 37.5%
E Max absolute difference: 2.83183781e-14
E Max relative difference: 1.92526432e-07
E x: array([[-1.41581383039852e-08, -1.61979841681622e-07,
E -1.68750114288469e-07, -1.60106129998105e-08],
E [-9.39713719856350e-08, -1.12524745467649e-08,...
E y: array([[-1.4158140e-08, -1.6197987e-07, -1.6875014e-07, -1.6010616e-08],
E [-9.3971386e-08, -1.1252476e-08, -1.1252476e-08, -9.7617234e-08],
E [-1.0785670e-07, -1.2403513e-08, -1.2403513e-08, -1.0364793e-07],
E [-2.9186946e-08, -1.7577703e-07, -1.6937879e-07, -2.6112047e-08]])
metpy/testing.py:151: AssertionError
_______________________________ test_q_vector_4d _______________________________
data_4d = D_4D_Test_Data(height=<Quantity([[[[5661.239 5572.5986 5500.839 5511.0386]
[5682.6787 5607.1587 5550.7188 5547.27...222107.84922049 -222107.84922049]]]], 'meter')>, latitude=<Quantity([[[[46.]
[44.]
[42.]
[40.]]]], 'degree')>)
def test_q_vector_4d(data_4d):
"""Test q_vector on a 4D (time, pressure, y, x) grid."""
f = coriolis_parameter(data_4d.latitude)
u_g, v_g = geostrophic_wind(data_4d.height, f, data_4d.dx, data_4d.dy)
q1, q2 = q_vector(u_g, v_g, data_4d.temperature, data_4d.pressure, data_4d.dx, data_4d.dy)
q1_truth = np.array([[[[-8.98245364e-13, 2.03803219e-13, 2.88874668e-12, 2.18043424e-12],
[4.37446820e-13, 1.21145200e-13, 1.51859353e-12, 3.82803347e-12],
[-1.20538030e-12, 2.27477298e-12, 3.47570178e-12, 3.03123012e-12],
[-1.51597275e-12, 8.02915408e-12, 7.71292472e-12, -2.22078527e-12]],
[[5.72960497e-13, 1.04264321e-12, -1.75695523e-13, 1.20745997e-12],
[2.94807953e-13, 5.80261767e-13, 6.23668595e-13, 7.31474131e-13],
[-4.04218965e-14, 3.24794013e-13, 1.39539675e-12, 2.82242029e-12],
[3.27509076e-13, 5.61307677e-13, 1.13454829e-12, 4.63551274e-12]],
[[2.23877015e-13, 5.77177907e-13, 1.62133659e-12, 5.43858376e-13],
[2.65333917e-13, 2.41006445e-13, 3.72510595e-13, 7.35822030e-13],
[6.56644633e-14, 1.99773842e-13, 5.20573457e-13, 1.69706608e-12],
[4.15915138e-14, 1.19910880e-13, 1.03632944e-12, 1.99561829e-12]]],
[[[-2.68870846e-13, 1.35977736e-12, 4.17548337e-12, 1.50465522e-12],
[4.62457018e-14, 1.25888111e-13, 2.15928418e-12, 4.70656495e-12],
[-1.25393137e-12, 9.54737370e-13, 1.48443002e-12, 2.12375621e-12],
[-2.93284658e-12, 6.06555344e-12, 4.21151397e-12, -2.12250513e-12]],
[[4.23461674e-13, 1.39393686e-13, 5.89509120e-13, 2.55041326e-12],
[5.73125714e-13, 5.60965341e-13, 7.65040451e-13, 9.49571939e-13],
[2.17153819e-14, 3.97023968e-13, 1.09194718e-12, 1.90731542e-12],
[1.45101233e-13, 1.79588608e-13, 1.03018848e-12, 3.62186462e-12]],
[[5.32674437e-13, 5.13465061e-13, 1.15582657e-12, 1.04827520e-12],
[2.77261345e-13, 2.33645555e-13, 4.59592371e-13, 5.34293340e-13],
[1.47376125e-13, 1.95746242e-13, 3.45854003e-13, 7.47741411e-13],
[-2.14078421e-14, 1.75226662e-13, 4.85424103e-13, 1.10808035e-12]]],
[[[6.41348753e-13, 1.88256910e-12, 5.21213092e-12, 2.07707653e-12],
[1.30753737e-12, 4.77125469e-13, 2.15204760e-12, 3.07374453e-12],
[-2.30546806e-13, 2.49929428e-13, 8.82215204e-14, 2.45990265e-12],
[-7.25812141e-12, 8.47072439e-13, -2.06762495e-12,
-4.40132129e-12]],
[[6.03705941e-13, -6.71320661e-13, 9.10543636e-13, 5.82480651e-13],
[9.54081741e-13, 6.11781160e-13, 6.95995265e-13, 8.67169047e-13],
[7.86580678e-14, 5.27405484e-13, 7.45800341e-13, 1.33965768e-12],
[2.22480631e-13, -1.98920384e-13, 8.56608245e-13, 1.59793218e-12]],
[[4.47195537e-13, 2.18235390e-13, 3.30926531e-13, -4.06675908e-14],
[1.70517246e-13, 2.18234962e-13, 3.78622612e-13, 5.03962144e-13],
[2.59462161e-13, 2.65626826e-13, 2.04642555e-13, 6.02812047e-13],
[1.69339642e-13, 2.91716502e-13, -1.20043003e-14,
4.43770388e-13]]]]) * units('m^2 kg^-1 s^-1')
q2_truth = np.array([[[[3.33980776e-12, -1.32969763e-13, 1.01454470e-12, 6.02652581e-12],
[2.52898242e-13, -1.71069245e-13, -8.24708561e-13, 1.66384429e-13],
[-3.50646511e-12, -1.68929195e-12, 7.76215111e-13, 1.54486058e-12],
[-1.75492099e-12, -3.86524071e-12, -1.89368596e-12,
-5.14689517e-12]],
[[-2.09848775e-13, -6.25683634e-13, -1.40009292e-13, 1.08972893e-12],
[-2.58259284e-13, -2.67211578e-13, -6.41928957e-14, 5.90625597e-13],
[-2.73346325e-13, -2.28248227e-13, -4.76577835e-13,
-8.48559875e-13],
[1.21003124e-12, -5.10541546e-13, 6.35947149e-14, 2.44893915e-12]],
[[-6.72309334e-14, -3.56791270e-13, -4.13553842e-14, 3.81212108e-13],
[-3.55860413e-13, -1.22880634e-13, -3.19443665e-14,
-4.71232601e-14],
[-2.82289531e-13, -1.20965929e-13, 1.14160715e-13, -6.85113982e-14],
[5.17465531e-14, -4.61129211e-13, 5.33793701e-13, 1.28285338e-12]]],
[[[1.71894904e-12, -1.35675428e-12, 1.48328005e-13, 3.22454170e-12],
[-2.12666583e-13, -1.17551681e-13, -6.93968059e-13, 1.76659826e-12],
[-2.67906914e-12, -3.78250861e-13, -9.88730956e-13, 2.88200442e-12],
[1.48225123e-12, 2.15004833e-13, -4.84554577e-12, 2.77242999e-12]],
[[-3.09626209e-13, -2.52138997e-13, 4.58311589e-14, 2.03206766e-12],
[-3.95662347e-13, -2.99828956e-13, 1.08715446e-14, 1.06096030e-12],
[-2.46502471e-13, -2.43524217e-13, -3.81250581e-13,
-1.70270366e-13],
[8.12479206e-13, -1.38629628e-13, -8.05591138e-13,
-7.80286006e-13]],
[[-2.19626566e-13, -1.52852503e-13, 4.07706963e-13, 1.52323163e-12],
[-2.56235985e-13, -1.20817691e-13, 6.51260820e-15, 3.49591511e-13],
[-2.44063890e-13, -1.21871642e-13, -9.09798480e-14,
-1.59903476e-13],
[-2.47929201e-13, -1.77269110e-13, -1.12991330e-13,
-6.06795348e-13]]],
[[[-6.48288201e-13, -1.96951432e-12, -5.53508048e-13, 1.94507133e-12],
[-2.00769011e-12, -3.72469047e-13, -4.59116219e-13, 1.11322705e-13],
[-3.83507643e-12, 1.18054543e-13, -4.24001455e-13, -5.88688871e-13],
[-1.84528711e-12, 1.54974343e-12, -7.36123184e-13, 1.06256777e-13]],
[[-4.58487019e-13, -1.89124158e-13, 2.58416604e-13, 8.14652306e-13],
[-6.09664269e-13, -3.51509413e-13, 2.39576397e-13, 5.80539044e-13],
[-1.68850738e-13, -3.49553817e-13, -2.26470205e-13, 7.79989044e-13],
[2.23081718e-13, 1.20195366e-13, -1.01508013e-12, -2.15527487e-13]],
[[-1.68054338e-13, -5.06878852e-14, 2.77697698e-13, 8.37521961e-13],
[-1.39462599e-13, -1.36628363e-13, 3.13920124e-14, 4.55413406e-13],
[-1.06658890e-13, -2.19817426e-13, -8.35968065e-14, 1.88190788e-13],
[-2.27182863e-13, -2.74607819e-13, -1.10587309e-13,
-3.88915866e-13]]]]) * units('m^2 kg^-1 s^-1')
> assert_array_almost_equal(q1, q1_truth, 20)
metpy/calc/tests/test_kinematics.py:1702:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
actual = array([[[[-8.98245420e-13, 2.03803231e-13, 2.88874686e-12,
2.18043437e-12],
[ 4.37446848e-13, 1... 6.02812085e-13],
[ 1.69339653e-13, 2.91716520e-13, -1.20043010e-14,
4.43770416e-13]]]])
desired = array([[[[-8.98245364e-13, 2.03803219e-13, 2.88874668e-12,
2.18043424e-12],
[ 4.37446820e-13, 1... 6.02812047e-13],
[ 1.69339642e-13, 2.91716502e-13, -1.20043003e-14,
4.43770388e-13]]]])
decimal = 20
def assert_array_almost_equal(actual, desired, decimal=7):
"""Check that arrays are almost equal, including units.
Wrapper around :func:`numpy.testing.assert_array_almost_equal`
"""
actual, desired = check_and_drop_units(actual, desired)
> numpy.testing.assert_array_almost_equal(actual, desired, decimal)
E AssertionError:
E Arrays are not almost equal to 20 decimals
E
E Mismatch: 77.8%
E Max absolute difference: 5.00771366e-19
E Max relative difference: 6.59987194e-08
E x: array([[[[-8.9824541957213273e-13, 2.0380323138284710e-13,
E 2.8887468637244123e-12, 2.1804343746946653e-12],
E [ 4.3744684750934834e-13, 1.2114520784515458e-13,...
E y: array([[[[-8.98245364e-13, 2.03803219e-13, 2.88874668e-12,
E 2.18043424e-12],
E [ 4.37446820e-13, 1.21145200e-13, 1.51859353e-12,...
metpy/testing.py:160: AssertionError
________________________________ test_cape_cin _________________________________
def test_cape_cin():
"""Test the basic CAPE and CIN calculation."""
p = np.array([959., 779.2, 751.3, 724.3, 700., 269.]) * units.mbar
temperature = np.array([22.2, 14.6, 12., 9.4, 7., -38.]) * units.celsius
dewpoint = np.array([19., -11.2, -10.8, -10.4, -10., -53.2]) * units.celsius
parcel_prof = parcel_profile(p, temperature[0], dewpoint[0])
cape, cin = cape_cin(p, temperature, dewpoint, parcel_prof)
> assert_almost_equal(cape, 75.7340825 * units('joule / kilogram'), 6)
metpy/calc/tests/test_thermo.py:698:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
actual = 75.73380786652106, desired = 75.7340825, decimal = 6
def assert_almost_equal(actual, desired, decimal=7):
"""Check that values are almost equal, including units.
Wrapper around :func:`numpy.testing.assert_almost_equal`
"""
actual, desired = check_and_drop_units(actual, desired)
> numpy.testing.assert_almost_equal(actual, desired, decimal)
E AssertionError:
E Arrays are not almost equal to 6 decimals
E ACTUAL: 75.73380786652106
E DESIRED: 75.7340825
metpy/testing.py:151: AssertionError
_____________________________ test_cape_cin_no_el ______________________________
def test_cape_cin_no_el():
"""Test that CAPE works with no EL."""
p = np.array([959., 779.2, 751.3, 724.3]) * units.mbar
temperature = np.array([22.2, 14.6, 12., 9.4]) * units.celsius
dewpoint = np.array([19., -11.2, -10.8, -10.4]) * units.celsius
parcel_prof = parcel_profile(p, temperature[0], dewpoint[0]).to('degC')
cape, cin = cape_cin(p, temperature, dewpoint, parcel_prof)
assert_almost_equal(cape, 0.08610409 * units('joule / kilogram'), 6)
> assert_almost_equal(cin, -89.8179205 * units('joule / kilogram'), 6)
metpy/calc/tests/test_thermo.py:710:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
actual = -89.81795976062058, desired = -89.8179205, decimal = 6
def assert_almost_equal(actual, desired, decimal=7):
"""Check that values are almost equal, including units.
Wrapper around :func:`numpy.testing.assert_almost_equal`
"""
actual, desired = check_and_drop_units(actual, desired)
> numpy.testing.assert_almost_equal(actual, desired, decimal)
E AssertionError:
E Arrays are not almost equal to 6 decimals
E ACTUAL: -89.81795976062058
E DESIRED: -89.8179205
metpy/testing.py:151: AssertionError
_________________________ test_surface_based_cape_cin __________________________
def test_surface_based_cape_cin():
"""Test the surface-based CAPE and CIN calculation."""
p = np.array([959., 779.2, 751.3, 724.3, 700., 269.]) * units.mbar
temperature = np.array([22.2, 14.6, 12., 9.4, 7., -38.]) * units.celsius
dewpoint = np.array([19., -11.2, -10.8, -10.4, -10., -53.2]) * units.celsius
cape, cin = surface_based_cape_cin(p, temperature, dewpoint)
> assert_almost_equal(cape, 75.7340825 * units('joule / kilogram'), 6)
metpy/calc/tests/test_thermo.py:952:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
actual = 75.73380786652106, desired = 75.7340825, decimal = 6
def assert_almost_equal(actual, desired, decimal=7):
"""Check that values are almost equal, including units.
Wrapper around :func:`numpy.testing.assert_almost_equal`
"""
actual, desired = check_and_drop_units(actual, desired)
> numpy.testing.assert_almost_equal(actual, desired, decimal)
E AssertionError:
E Arrays are not almost equal to 6 decimals
E ACTUAL: 75.73380786652106
E DESIRED: 75.7340825
metpy/testing.py:151: AssertionError
_____________________ test_most_unstable_cape_cin_surface ______________________
def test_most_unstable_cape_cin_surface():
"""Test the most unstable CAPE/CIN calculation when surface is most unstable."""
pressure = np.array([959., 779.2, 751.3, 724.3, 700., 269.]) * units.mbar
temperature = np.array([22.2, 14.6, 12., 9.4, 7., -38.]) * units.celsius
dewpoint = np.array([19., -11.2, -10.8, -10.4, -10., -53.2]) * units.celsius
mucape, mucin = most_unstable_cape_cin(pressure, temperature, dewpoint)
> assert_almost_equal(mucape, 75.7340825 * units('joule / kilogram'), 6)
metpy/calc/tests/test_thermo.py:962:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
actual = 75.73380786652106, desired = 75.7340825, decimal = 6
def assert_almost_equal(actual, desired, decimal=7):
"""Check that values are almost equal, including units.
Wrapper around :func:`numpy.testing.assert_almost_equal`
"""
actual, desired = check_and_drop_units(actual, desired)
> numpy.testing.assert_almost_equal(actual, desired, decimal)
E AssertionError:
E Arrays are not almost equal to 6 decimals
E ACTUAL: 75.73380786652106
E DESIRED: 75.7340825
metpy/testing.py:151: AssertionError
__________________________ test_wet_bulb_temperature ___________________________
def test_wet_bulb_temperature():
"""Test wet bulb calculation with scalars."""
val = wet_bulb_temperature(1000 * units.hPa, 25 * units.degC, 15 * units.degC)
truth = 18.34345936 * units.degC # 18.59 from NWS calculator
> assert_almost_equal(val, truth)
metpy/calc/tests/test_thermo.py:1123:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
actual = 18.343459589135325, desired = 18.34345936, decimal = 7
def assert_almost_equal(actual, desired, decimal=7):
"""Check that values are almost equal, including units.
Wrapper around :func:`numpy.testing.assert_almost_equal`
"""
actual, desired = check_and_drop_units(actual, desired)
> numpy.testing.assert_almost_equal(actual, desired, decimal)
E AssertionError:
E Arrays are not almost equal to 7 decimals
E ACTUAL: 18.343459589135325
E DESIRED: 18.34345936
metpy/testing.py:151: AssertionError
____________________________ test_lfc_not_below_lcl ____________________________
def test_lfc_not_below_lcl():
"""Test sounding where LFC appears to be (but isn't) below LCL."""
levels = np.array([1002.5, 1001.7, 1001., 1000.3, 999.7, 999., 998.2, 977.9,
966.2, 952.3, 940.6, 930.5, 919.8, 909.1, 898.9, 888.4,
878.3, 868.1, 858., 848., 837.2, 827., 816.7, 805.4]) * units.hPa
temperatures = np.array([17.9, 17.9, 17.8, 17.7, 17.7, 17.6, 17.5, 16.,
15.2, 14.5, 13.8, 13., 12.5, 11.9, 11.4, 11.,
10.3, 9.7, 9.2, 8.7, 8., 7.4, 6.8, 6.1]) * units.degC
dewpoints = np.array([13.6, 13.6, 13.5, 13.5, 13.5, 13.5, 13.4, 12.5,
12.1, 11.8, 11.4, 11.3, 11., 9.3, 10., 8.7, 8.9,
8.6, 8.1, 7.6, 7., 6.5, 6., 5.4]) * units.degC
lfc_pressure, lfc_temp = lfc(levels, temperatures, dewpoints)
# Before patch, LFC pressure would show 1000.5912165339967 hPa
> assert_almost_equal(lfc_pressure, 811.8263397 * units.mbar, 6)
metpy/calc/tests/test_thermo.py:1201:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
actual = 811.8262615743613, desired = 811.8263397, decimal = 6
def assert_almost_equal(actual, desired, decimal=7):
"""Check that values are almost equal, including units.
Wrapper around :func:`numpy.testing.assert_almost_equal`
"""
actual, desired = check_and_drop_units(actual, desired)
> numpy.testing.assert_almost_equal(actual, desired, decimal)
E AssertionError:
E Arrays are not almost equal to 6 decimals
E ACTUAL: 811.8262615743613
E DESIRED: 811.8263397
metpy/testing.py:151: AssertionError
______________________________ test_multiple_lfcs ______________________________
def test_multiple_lfcs():
"""Test sounding with multiple LFCs.
If which='top', return lowest-pressure LFC.
If which='all', return all LFCs
"""
levels = np.array([966., 937.2, 925., 904.6, 872.6, 853., 850., 836., 821., 811.6, 782.3,
754.2, 726.9, 700., 648.9, 624.6, 601.1, 595., 587., 576., 555.7,
534.2, 524., 500., 473.3, 400., 384.5, 358., 343., 308.3, 300., 276.,
273., 268.5, 250., 244.2, 233., 200.]) * units.mbar
temperatures = np.array([18.2, 16.8, 16.2, 15.1, 13.3, 12.2, 12.4, 14., 14.4,
13.7, 11.4, 9.1, 6.8, 4.4, -1.4, -4.4, -7.3, -8.1,
-7.9, -7.7, -8.7, -9.8, -10.3, -13.5, -17.1, -28.1, -30.7,
-35.3, -37.1, -43.5, -45.1, -49.9, -50.4, -51.1, -54.1, -55.,
-56.7, -57.5]) * units.degC
dewpoints = np.array([16.9, 15.9, 15.5, 14.2, 12.1, 10.8, 8.6, 0., -3.6, -4.4,
-6.9, -9.5, -12., -14.6, -15.8, -16.4, -16.9, -17.1, -27.9, -42.7,
-44.1, -45.6, -46.3, -45.5, -47.1, -52.1, -50.4, -47.3, -57.1,
-57.9, -58.1, -60.9, -61.4, -62.1, -65.1, -65.6,
-66.7, -70.5]) * units.degC
lfc_pressure_top, lfc_temp_top = lfc(levels, temperatures, dewpoints)
lfc_pressure_bottom, lfc_temp_bottom = lfc(levels, temperatures, dewpoints,
which='bottom')
lfc_pressure_all, _ = lfc(levels, temperatures, dewpoints, which='all')
> assert_almost_equal(lfc_pressure_top, 705.4346277 * units.mbar, 6)
metpy/calc/tests/test_thermo.py:1230:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
actual = 705.4345895435765, desired = 705.4346277, decimal = 6
def assert_almost_equal(actual, desired, decimal=7):
"""Check that values are almost equal, including units.
Wrapper around :func:`numpy.testing.assert_almost_equal`
"""
actual, desired = check_and_drop_units(actual, desired)
> numpy.testing.assert_almost_equal(actual, desired, decimal)
E AssertionError:
E Arrays are not almost equal to 6 decimals
E ACTUAL: 705.4345895435765
E DESIRED: 705.4346277
metpy/testing.py:151: AssertionError
______________________________ test_multiple_els _______________________________
def test_multiple_els():
"""Test sounding with multiple ELs.
If which='top', return lowest-pressure EL.
If which='all', return all ELs
"""
levels = np.array([966., 937.2, 925., 904.6, 872.6, 853., 850., 836., 821., 811.6, 782.3,
754.2, 726.9, 700., 648.9, 624.6, 601.1, 595., 587., 576., 555.7,
534.2, 524., 500., 473.3, 400., 384.5, 358., 343., 308.3, 300., 276.,
273., 268.5, 250., 244.2, 233., 200.]) * units.mbar
temperatures = np.array([18.2, 16.8, 16.2, 15.1, 13.3, 12.2, 12.4, 14., 14.4,
13.7, 11.4, 9.1, 6.8, 4.4, -1.4, -4.4, -7.3, -8.1,
-7.9, -7.7, -8.7, -9.8, -10.3, -13.5, -17.1, -28.1, -30.7,
-35.3, -37.1, -43.5, -45.1, -49.9, -50.4, -51.1, -54.1, -55.,
-56.7, -57.5]) * units.degC
dewpoints = np.array([16.9, 15.9, 15.5, 14.2, 12.1, 10.8, 8.6, 0., -3.6, -4.4,
-6.9, -9.5, -12., -14.6, -15.8, -16.4, -16.9, -17.1, -27.9, -42.7,
-44.1, -45.6, -46.3, -45.5, -47.1, -52.1, -50.4, -47.3, -57.1,
-57.9, -58.1, -60.9, -61.4, -62.1, -65.1, -65.6,
-66.7, -70.5]) * units.degC
el_pressure_top, el_temp_top = el(levels, temperatures, dewpoints)
el_pressure_bottom, el_temp_bottom = el(levels, temperatures, dewpoints, which='bottom')
el_pressure_all, _ = el(levels, temperatures, dewpoints, which='all')
> assert_almost_equal(el_pressure_top, 228.0575059 * units.mbar, 6)
metpy/calc/tests/test_thermo.py:1261:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
actual = 228.05753786963396, desired = 228.0575059, decimal = 6
def assert_almost_equal(actual, desired, decimal=7):
"""Check that values are almost equal, including units.
Wrapper around :func:`numpy.testing.assert_almost_equal`
"""
actual, desired = check_and_drop_units(actual, desired)
> numpy.testing.assert_almost_equal(actual, desired, decimal)
E AssertionError:
E Arrays are not almost equal to 6 decimals
E ACTUAL: 228.05753786963396
E DESIRED: 228.0575059
metpy/testing.py:151: AssertionError
_________________________ test_cape_cin_custom_profile _________________________
def test_cape_cin_custom_profile():
"""Test the CAPE and CIN calculation with a custom profile passed to LFC and EL."""
p = np.array([959., 779.2, 751.3, 724.3, 700., 269.]) * units.mbar
temperature = np.array([22.2, 14.6, 12., 9.4, 7., -38.]) * units.celsius
dewpoint = np.array([19., -11.2, -10.8, -10.4, -10., -53.2]) * units.celsius
parcel_prof = parcel_profile(p, temperature[0], dewpoint[0]) + 5 * units.delta_degC
cape, cin = cape_cin(p, temperature, dewpoint, parcel_prof)
> assert_almost_equal(cape, 1443.505086499895 * units('joule / kilogram'), 6)
metpy/calc/tests/test_thermo.py:1275:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
actual = 1443.5040016400164, desired = 1443.505086499895, decimal = 6
def assert_almost_equal(actual, desired, decimal=7):
"""Check that values are almost equal, including units.
Wrapper around :func:`numpy.testing.assert_almost_equal`
"""
actual, desired = check_and_drop_units(actual, desired)
> numpy.testing.assert_almost_equal(actual, desired, decimal)
E AssertionError:
E Arrays are not almost equal to 6 decimals
E ACTUAL: 1443.5040016400164
E DESIRED: 1443.505086499895
metpy/testing.py:151: AssertionError
_________________________________ test_axhline _________________________________
self = <Quantity(930, 'millibar')>, key = 0
def __getitem__(self, key):
try:
> value = self._magnitude[key]
E TypeError: 'int' object is not subscriptable
../pint/pint/quantity.py:1465: TypeError
During handling of the above exception, another exception occurred:
@pytest.mark.skipif(pint.__version__ == '0.9', reason=('Currently broken upstream (see '
'pint#751'))
@pytest.mark.mpl_image_compare(tolerance=0, remove_text=True)
def test_axhline():
r"""Ensure that passing a quantity to axhline does not error."""
fig, ax = plt.subplots()
> ax.axhline(930 * units('mbar'))
metpy/tests/test_units.py:45:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
../../miniconda3/envs/dev/lib/python3.7/site-packages/matplotlib/axes/_axes.py:805: in axhline
yy = self.convert_yunits(y)
../../miniconda3/envs/dev/lib/python3.7/site-packages/matplotlib/artist.py:195: in convert_yunits
return ax.yaxis.convert_units(y)
../../miniconda3/envs/dev/lib/python3.7/site-packages/matplotlib/axis.py:1521: in convert_units
if munits.ConversionInterface.is_numlike(x):
../../miniconda3/envs/dev/lib/python3.7/site-packages/matplotlib/units.py:126: in is_numlike
for thisx in x:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
self = <Quantity(930, 'millibar')>, key = 0
def __getitem__(self, key):
try:
value = self._magnitude[key]
return self.__class__(value, self._units)
except TypeError:
raise TypeError("Neither Quantity object nor its magnitude ({})"
> "supports indexing".format(self._magnitude))
E TypeError: Neither Quantity object nor its magnitude (930)supports indexing
../pint/pint/quantity.py:1469: TypeError
_________________________________ test_axvline _________________________________
self = <Quantity(0, 'degree_Celsius')>, key = 0
def __getitem__(self, key):
try:
> value = self._magnitude[key]
E TypeError: 'int' object is not subscriptable
../pint/pint/quantity.py:1465: TypeError
During handling of the above exception, another exception occurred:
@pytest.mark.skipif(pint.__version__ == '0.9', reason=('Currently broken upstream (see '
'pint#751'))
@pytest.mark.mpl_image_compare(tolerance=0, remove_text=True)
def test_axvline():
r"""Ensure that passing a quantity to axvline does not error."""
fig, ax = plt.subplots()
> ax.axvline(0 * units('degC'))
metpy/tests/test_units.py:57:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
../../miniconda3/envs/dev/lib/python3.7/site-packages/matplotlib/axes/_axes.py:874: in axvline
xx = self.convert_xunits(x)
../../miniconda3/envs/dev/lib/python3.7/site-packages/matplotlib/artist.py:186: in convert_xunits
return ax.xaxis.convert_units(x)
../../miniconda3/envs/dev/lib/python3.7/site-packages/matplotlib/axis.py:1521: in convert_units
if munits.ConversionInterface.is_numlike(x):
../../miniconda3/envs/dev/lib/python3.7/site-packages/matplotlib/units.py:126: in is_numlike
for thisx in x:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
self = <Quantity(0, 'degree_Celsius')>, key = 0
def __getitem__(self, key):
try:
value = self._magnitude[key]
return self.__class__(value, self._units)
except TypeError:
raise TypeError("Neither Quantity object nor its magnitude ({})"
> "supports indexing".format(self._magnitude))
E TypeError: Neither Quantity object nor its magnitude (0)supports indexing
../pint/pint/quantity.py:1469: TypeError
______________________________ test_convert_units ______________________________
test_var = <xarray.DataArray 'Temperature' (time: 1, isobaric: 29, y: 118, x: 292)>
array([[[[ 18.440033, ..., 30.752533],
...e: Initialized analysis product
GRIB_level_type: 100
GRIB_VectorComponentFlag: gridRelative
def test_convert_units(test_var):
"""Test in-place conversion of units."""
test_var.metpy.convert_units('degC')
# Check that variable metadata is updated
> assert test_var.attrs['units'] == 'degC'
E AssertionError: assert 'degree_Celsius' == 'degC'
E - degree_Celsius
E + degC
metpy/tests/test_xarray.py:100: AssertionError
=============================== warnings summary ===============================
metpy/calc/tests/test_basic.py::test_windchill_invalid
/home/jon/miniconda3/envs/dev/lib/python3.7/site-packages/numpy/testing/_private/utils.py:1011: RuntimeWarning: invalid value encountered in less
return z < 1.5 * 10.0**(-decimal)
metpy/calc/tests/test_basic.py::test_heat_index_invalid
/home/jon/dev/pint/pint/quantity.py:1166: RuntimeWarning: invalid value encountered in sqrt
magnitude = new_self._magnitude ** _to_magnitude(other, self.force_ndarray)
metpy/plots/tests/test_declarative.py::test_latlon
/home/jon/miniconda3/envs/dev/lib/python3.7/site-packages/cartopy/mpl/geoaxes.py:1063: DeprecationWarning: Using or importing the ABCs from 'collections' instead of from 'collections.abc' is deprecated, and in 3.8 it will stop working
if not isinstance(regrid_shape, collections.Sequence):
-- Docs: https://docs.pytest.org/en/latest/warnings.html
======== 25 failed, 852 passed, 1 xfailed, 3 warnings in 58.55 seconds =========
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment