address usage warnings

src/ravenpy/config/rvs.py

@@ -29,9 +29,9 @@ class RVI(RV):
 
     run_name: Optional[str] = optfield(alias="RunName")
     calendar: Optional[o.Calendar] = optfield(alias="Calendar")
-    start_date: Optional[date] = optfield(alias="StartDate")
+    start_date: Optional[Union[str, date]] = optfield(alias="StartDate")
     assimilation_start_time: Optional[date] = optfield(alias="AssimilationStartTime")
-    end_date: Optional[date] = optfield(alias="EndDate")
+    end_date: Optional[Union[str, date]] = optfield(alias="EndDate")
     duration: Optional[float] = optfield(alias="Duration")
     time_step: Optional[Union[float, str]] = optfield(alias="TimeStep")
     interpolation: Optional[o.Interpolation] = optfield(alias="Interpolation")

src/ravenpy/extractors/routing_product.py

@@ -14,6 +14,7 @@
     from osgeo import ogr, osr
     from shapely import wkt
 
+    osr.UseExceptions()
 except (ImportError, ModuleNotFoundError) as e:
     msg = gis_import_error_message.format(Path(__file__).stem)
     raise ImportError(msg) from e

src/ravenpy/utilities/analysis.py

@@ -10,8 +10,10 @@
 
 try:
     import rasterio
-    from osgeo.gdal import Dataset, DEMProcessing
+    from osgeo.gdal import Dataset, DEMProcessing, UseExceptions
     from shapely.geometry import GeometryCollection, MultiPolygon, Polygon, shape
+
+    UseExceptions()
 except (ImportError, ModuleNotFoundError) as e:
     msg = gis_import_error_message.format(Path(__file__).stem)
     raise ImportError(msg) from e

src/ravenpy/utilities/forecasting.py

@@ -374,7 +374,9 @@ def ensemble_prediction(
 
     # Run the model for each year
     ensemble = []
-    forecast_ds = xr.open_dataset(forecast, use_cftime=True)
+
+    time_coder = xr.coders.CFDatetimeCoder(use_cftime=True)
+    forecast_ds = xr.open_dataset(forecast, decode_times=time_coder)
 
     for member in range(0, len(forecast_ds[ens_dim])):
         # Prepare model instance

tests/test_forecasting.py

@@ -1,4 +1,4 @@
-import datetime as dt
+from datetime import datetime as dt
 
 import numpy as np
 import pytest
@@ -75,7 +75,7 @@ def test_forecasting_GEPS(numeric_config, get_local_testdata):
     # Prepare a RAVEN model run using historical data, GR4JCN in this case.
     # This is a dummy run to get initial states. In a real forecast situation,
     # this run would end on the day before the forecast, but process is the same.
-    wup.start_date = dt.datetime(2000, 1, 1)
+    wup.start_date = dt(2000, 1, 1)
     wup.duration = 30
 
     e = Emulator(wup)

tests/test_geoserver.py

@@ -76,7 +76,7 @@ def test_hydrobasins_upstream_aggregate(self, tmp_path):
         aggregated = self.geoserver.hydrobasins_aggregate(gdf_upstream)
 
         assert len(aggregated) == 1
-        assert float(aggregated.SUB_AREA.values) == 4977.8
+        assert aggregated.SUB_AREA.values[0] == 4977.8
         np.testing.assert_equal(
             aggregated.geometry.bounds.values,
             np.array([[-83.8167, 8.7625, -82.7125, 9.5875]]),