Fix regridding

Author: bouweandela

pyproject.toml

@@ -23,7 +23,6 @@ matplotlib = ">=3.10.3,<4"
 scipy = ">=1.16.0,<2"
 xarray = ">=2025.7.1,<2026"
 typer = ">=0.16.0,<0.17"
-xcdat = ">=0.9.1,<0.10"
 xesmf = ">=0.8.10,<0.9"
 
 [tool.pixi.pypi-dependencies]

src/ref_sample_data/resample.py

@@ -1,6 +1,5 @@
 import numpy as np
 import xarray as xr
-import xcdat
 import xesmf
 
 
@@ -44,7 +43,40 @@ def decimate_rectilinear(dataset: xr.Dataset) -> xr.Dataset:
     """
     # Decimate the dataset, but update the bounds
     # 10x10 degree grid
-    output_grid = xcdat.create_uniform_grid(-90, 90, 10, 0, 359, 10)
+    grid_size = 10  # degrees
+    nlat = int(180 / grid_size) + 1
+    nlon = int(360 / grid_size)
+    output_grid = xr.Dataset.from_dict(
+        {
+            "attrs": {},
+            "coords": {
+                "lat": {
+                    "attrs": {
+                        "axis": "Y",
+                        "coordinate": "latitude",
+                        "units": "degrees_north",
+                    },
+                    "data": np.linspace(-90, 90, nlat, endpoint=True),
+                    "dims": ("lat",),
+                },
+                "lon": {
+                    "attrs": {
+                        "axis": "X",
+                        "coordinate": "longitude",
+                        "units": "degrees_east",
+                    },
+                    "data": np.linspace(0, 360, nlon, endpoint=False),
+                    "dims": ("lon",),
+                },
+            },
+            "data_vars": {},
+            "dims": {
+                "lat": nlat,
+                "lon": nlon,
+            },
+        }
+    )
+
     regrid = xesmf.Regridder(dataset, output_grid, "bilinear", periodic=True)
     result = regrid(dataset.copy())
     result = result.bounds.add_bounds("Y").bounds.add_bounds("X")