feat: add backwards compatibility for variable query param

Author: hrodmn

tests/test_dependencies.py

@@ -27,7 +27,7 @@ def test_granule_search_temporal_interval_unchanged():
 def test_interpolated_xarray_params_single_datetime():
     """Test InterpolatedXarrayParams with single datetime interpolation."""
     xarray_params = InterpolatedXarrayParams(
-        variable="temperature", sel=["time={datetime}", "lev=1000"]
+        variables=["temperature"], sel=["time={datetime}", "lev=1000"]
     )
 
     single_datetime = datetime(2025, 9, 23, 0, 0, 0, tzinfo=timezone.utc)
@@ -39,13 +39,13 @@ def test_interpolated_xarray_params_single_datetime():
     result = dependencies.interpolated_xarray_ds_params(xarray_params, granule_search)
 
     assert result.sel == [f"time={single_datetime.isoformat()}", "lev=1000"]
-    assert result.variable == "temperature"
+    assert result.variables == ["temperature"]
 
 
 def test_interpolated_xarray_params_datetime_range():
     """Test InterpolatedXarrayParams with datetime range (uses start datetime)."""
     xarray_params = InterpolatedXarrayParams(
-        variable="temperature", sel=["time={datetime}"]
+        variables=["temperature"], sel=["time={datetime}"]
     )
 
     start_datetime = datetime(2025, 9, 23, 0, 0, 0, tzinfo=timezone.utc)
@@ -63,7 +63,7 @@ def test_interpolated_xarray_params_datetime_range():
 def test_interpolated_xarray_params_no_datetime_template():
     """Test InterpolatedXarrayParams when sel doesn't contain datetime template."""
     xarray_params = InterpolatedXarrayParams(
-        variable="temperature",
+        variables=["temperature"],
         sel=["time=2025-01-01T00:00:00Z", "lev=1000"],
     )
 
@@ -80,7 +80,7 @@ def test_interpolated_xarray_params_no_datetime_template():
 
 def test_interpolated_xarray_params_no_sel():
     """Test InterpolatedXarrayParams when sel is None or empty."""
-    xarray_params = InterpolatedXarrayParams(variable="temperature", sel=None)
+    xarray_params = InterpolatedXarrayParams(variables=["temperature"], sel=None)
 
     single_datetime = datetime(2025, 9, 23, 0, 0, 0, tzinfo=timezone.utc)
     granule_search = GranuleSearch(
@@ -91,13 +91,13 @@ def test_interpolated_xarray_params_no_sel():
     result = dependencies.interpolated_xarray_ds_params(xarray_params, granule_search)
 
     assert result.sel is None
-    assert result.variable == "temperature"
+    assert result.variables == ["temperature"]
 
 
 def test_interpolated_xarray_params_multiple_templates():
     """Test InterpolatedXarrayParams with multiple datetime templates."""
     xarray_params = InterpolatedXarrayParams(
-        variable="temperature",
+        variables=["temperature"],
         sel=["time={datetime}", "start_time={datetime}", "lev=1000"],
     )
 

titiler/cmr/dependencies.py

@@ -115,7 +115,9 @@ def __post_init__(self):
 class XarrayDsParams(DefaultDependency):
     """Xarray Dataset Options."""
 
-    variables: Annotated[list[str], Query(description="Xarray Variable name.")]
+    variables: Annotated[
+        list[str], Query(alias="variable", description="Xarray Variable name.")
+    ]
 
     sel: Annotated[
         list[SelDimStr] | None,

titiler/cmr/reader.py

@@ -134,7 +134,7 @@ def get_variables(
         xarray.DataArray: 2D or 3D DataArray.
 
     """
-    da = xr.concat([ds[variable] for variable in variables], dim="band")
+    da = xr.concat([ds[variable] for variable in variables], dim="band").squeeze()
 
     for selector in _parse_dsl(sel):
         dimension = selector["dimension"]