Reduce DeprecationWarnings in Testing

docs/release-notes.rst

@@ -5,6 +5,14 @@
 Release Notes
 =============
 
+v2024.03.0 (Unreleased)
+---------------------------
+Upcoming release
+
+Internal Changes
+^^^^^^^^^^^^^^^^
+* Reduce DeprecationWarnings in testing by `Cora Schneck`_ in (:pr:`582`)
+
 v2024.03.0 (March 29, 2024)
 ---------------------------
 This release includes a bug fix for ``delta_pressure``.

geocat/comp/interpolation.py

@@ -459,18 +459,17 @@ def interp_hybrid_to_pressure(data: xr.DataArray,
     #     }},
     # )
 
-    # If an unchunked Xarray input is given, chunk it just with its dims
-    if data.chunks is None:
-        data_chunk = dict([
-            (k, v) for (k, v) in zip(list(data.dims), list(data.shape))
-        ])
-        data = data.chunk(data_chunk)
+    # set chunked Xarray input to chunk with dims
+    data_chunk = dict([
+        (k, v) for (k, v) in zip(list(data.dims), list(data.shape))
+    ])
+    data = data.chunk(data_chunk)
 
     # Chunk pressure equal to data's chunks
-    pressure = pressure.chunk(data.chunks)
+    pressure = pressure.chunk(data_chunk)
 
     # Output data structure elements
-    out_chunks = list(data.chunks)
+    out_chunks = list(zip((data_chunk.values())))
     out_chunks[interp_axis] = (new_levels.size,)
     out_chunks = tuple(out_chunks)
     # ''' end of boilerplate

test/test_dask.py

@@ -45,6 +45,7 @@ def test_dewtemp_dask(self):
         assert isinstance((out - 273.15).data, dask.array.Array)
         assert np.allclose(out - 273.15, dt_2, atol=0.1)
 
+    @pytest.mark.filterwarnings("ignore:invalid value encountered in sqrt")
     def test_heat_index_dask(self):
 
         ncl_gt_1 = [

test/test_interpolation.py

@@ -172,6 +172,8 @@ def test_interp_hybrid_to_pressure_extrap_temp(self, temp_in, press_in,
         temp_extrap_expected = ds_out.Tpx.rename(lev_p='plev')
         xr.testing.assert_allclose(temp_extrap_expected, result)
 
+    @pytest.mark.filterwarnings(
+        "ignore: Interpolation point out of data bounds encountered")
     def test_interp_hybrid_to_pressure_extrap_geopotential(
             self, geopotential_in, press_in, _hyam, _hybm, t_bot, phis,
             ds_out) -> None:

test/test_meteorology.py

@@ -96,6 +96,9 @@ def test_alt_coef(self) -> None:
                            self.ncl_gt_2,
                            atol=0.005)
 
+    @pytest.mark.filterwarnings(
+        "ignore:invalid value encountered in sqrt"
+    )  # ignores RuntimeWarning sqrt is negative, np.sqrt produces nan
     def test_xarray_alt_coef(self) -> None:
         assert np.allclose(heat_index(xr.DataArray(self.t2),
                                       xr.DataArray(self.rh2), True),
@@ -110,6 +113,9 @@ def test_list_input(self) -> None:
                            self.ncl_gt_1,
                            atol=0.005)
 
+    @pytest.mark.filterwarnings(
+        "ignore:invalid value encountered in sqrt"
+    )  # ignores RuntimeWarning sqrt is negative, np.sqrt produces nan
     def test_xarray_input(self) -> None:
         t = xr.DataArray(self.t1)
         rh = xr.DataArray(self.rh1)