skip msftmz in decimation

Author: nocollier

src/ref_sample_data/data_request/cmip6.py

@@ -106,14 +106,17 @@ def decimate_dataset(self, dataset: xr.Dataset) -> xr.Dataset | None:
         has_latlon = "lat" in dataset.dims and "lon" in dataset.dims
         has_ij = "i" in dataset.dims and "j" in dataset.dims
 
+        # The AMOC variable `msftmz` has these strange dims and we do not want to decimate
+        skip_decimate = {"time", "basin", "lev", "lat"}.issubset(dataset.dims)
+
         if has_latlon:
             assert len(dataset.lat.dims) == 1 and len(dataset.lon.dims) == 1
 
             result = decimate_rectilinear(dataset)
         elif has_ij:
             # 2d curvilinear grid (generally ocean variables)
             result = decimate_curvilinear(dataset)
-        else:
+        elif not skip_decimate:
             raise ValueError("Cannot decimate this grid: too many dimensions")
 
         if "time" in dataset.dims and self.time_span is not None: