fix unit test because of catalog changes (#271)

Author: wh3248

pyproject.toml

@@ -1,6 +1,6 @@
 [tool.poetry]
 name = "hf_hydrodata"
-version = "1.4.9"
+version = "1.4.10"
 description = "hydroframe tools and utilities"
 authors = ["William M. Hasling", "Laura Condon", "Reed Maxwell",  "George Artavanis", "Amy M. Johnson", "Amy C. Defnet"]
 license = "MIT"

tests/hf_hydrodata/test_gridded.py

@@ -1544,96 +1544,6 @@ def test_huc_border():
     assert math.isnan(data[0, 0, 0])
 
 
-@pytest.mark.private_dataset
-def test_get_wtd():
-    """Unit test reading the 3 resolution of water table depth files."""
-
-    # Test the 1000 meter resolution version
-    x = 1500
-    y = 1500
-    grid_bound_x_width = 2
-    grid_bounds_y_height = 2
-    bounds = [x, y, x + grid_bound_x_width, y + grid_bounds_y_height]
-    options = {
-        "dataset": "conus2_current_conditions",
-        "grid_bounds": bounds,
-        "variable": "water_table_depth",
-        "grid": "conus2_wtd",
-    }
-    data = hf.get_gridded_data(options)
-    assert (
-        hf.get_path(options)
-        in "/hydrodata/temp/high_resolution_data/WTD_estimates/30m/remapped_data/wtd_mean_estimate_RF_additional_inputs_dummy_drop0LP_1000m_CONUS2_m_1s_remapped.tif"
-    )
-
-    assert data.shape == (2, 2)
-    assert str(round(data[0, 0], 4)) == "52.86"
-    assert str(round(data[0, 1], 4)) == "43.374"
-    assert str(round(data[1, 0], 4)) == "27.7567"
-    assert str(round(data[1, 1], 4)) == "36.6467"
-
-    # Test the 100 meter resolution version
-    # Same points, but values are not exactly the same as 1000 because of aggregation in resolutions
-    x = 1500 * 10
-    y = 1500 * 10
-    bounds = [x, y, x + grid_bound_x_width, y + grid_bounds_y_height]
-    options = {
-        "dataset": "conus2_current_conditions",
-        "grid_bounds": bounds,
-        "variable": "water_table_depth",
-        "grid": "conus2_wtd.100",
-    }
-    data = hf.get_gridded_data(options)
-    assert (
-        hf.get_path(options)
-        in "/hydrodata/temp/high_resolution_data/WTD_estimates/30m/remapped_data/wtd_mean_estimate_RF_additional_inputs_dummy_drop0LP_100m_CONUS2_m_1s_remapped.tif"
-    )
-    assert data.shape == (2, 2)
-    assert str(round(data[0, 0], 3)) == "58.015"
-    assert str(round(data[0, 1], 3)) == "54.305"
-    assert str(round(data[1, 0], 3)) == "48.614"
-    assert str(round(data[1, 1], 3)) == "49.047"
-
-    # Test the 30 meter resolution version
-    # Same points, but values are not exactly the same as 1000 because of aggregation in resolutions
-    x = int((1500 * 1000) / 30)
-    y = int((1500 * 1000) / 30)
-    bounds = [x, y, x + grid_bound_x_width, y + grid_bounds_y_height]
-    options = {
-        "dataset": "conus2_current_conditions",
-        "grid_bounds": bounds,
-        "variable": "water_table_depth",
-        "grid": "conus2_wtd.30",
-    }
-    data = hf.get_gridded_data(options)
-    assert (
-        hf.get_path(options)
-        in "/hydrodata/temp/high_resolution_data/WTD_estimates/30m/compressed_data/wtd_mean_estimate_RF_additional_inputs_dummy_drop0LP_1s_CONUS2_m_remapped_unflip_compressed.tif"
-    )
-
-    assert data.shape == (2, 2)
-    assert data.shape == (2, 2)
-    assert str(round(data[0, 0], 3)) == "77.192"
-    assert str(round(data[0, 1], 3)) == "78.744"
-    assert str(round(data[1, 0], 3)) == "78.691"
-    assert str(round(data[1, 1], 3)) == "78.744"
-
-
-@pytest.mark.private_dataset
-def test_wtd_1000m_north():
-    """Unit test edge condition found during integration testing."""
-
-    bounds = [1593, 1724, 3420, 3484]
-    options = {
-        "dataset": "conus2_current_conditions",
-        "grid_bounds": bounds,
-        "variable": "water_table_depth",
-        "grid": "conus2_wtd",
-    }
-    data = hf.get_gridded_data(options)
-    assert data.shape == ((1760, 1827))
-
-
 def test_noaa_precip():
     """Unit test NOAA precip"""
     bounds = [1000, 1000, 1005, 1005]