Add extra test

codecov.yaml

@@ -5,5 +5,5 @@ coverage:
   status:
     project:
       default:
-        target: 80%
+        target: 75%
         threshold: 0.5%

eo_tides/model.py

@@ -885,10 +885,9 @@ def tag_tides(
     **model_tides_kwargs,
 ):
     """
-    Model tide heights for every timestep in a multi-
-    dimensional dataset, and add them as a new `tide_height`
-    variable giving the height of the tide at the exact
-    moment of each acquisition.
+    Model tide heights for every timestep in a multi-dimensional
+    dataset, and add them as a new `tide_height` (and optionally,
+    `ebb_flow`) variable that "tags" each observation with tide data.
 
     The function models tides at the centroid of the dataset
     by default, but a custom tidal modelling location can

tests/test_model.py

@@ -506,6 +506,22 @@ def test_tag_tides(satellite_ds, measured_tides_ds, ebb_flow, swap_dims, tidepos
         assert abs(val_stats["Bias"]) < 0.20
 
 
+def test_tag_tides_multiple(satellite_ds):
+    # Model multiple models at once
+    tagged_tides_ds = tag_tides(satellite_ds, model=["FES2014", "HAMTIDE11"], ebb_flow=True)
+
+    assert "tide_model" in tagged_tides_ds.dims
+    assert tagged_tides_ds.tide_height.dims == ("time", "tide_model")
+    assert tagged_tides_ds.ebb_flow.dims == ("time", "tide_model")
+
+    # Test that multiple tide models are correlated
+    val_stats = eval_metrics(
+        x=tagged_tides_ds.sel(tide_model="FES2014").tide_height,
+        y=tagged_tides_ds.sel(tide_model="HAMTIDE11").tide_height,
+    )
+    assert val_stats["Correlation"] > 0.99
+
+
 # Run tests for default and custom resolutions
 @pytest.mark.parametrize("resolution", [None, "custom"])
 def test_pixel_tides(satellite_ds, measured_tides_ds, resolution):