What is your issue?

For regridding (rectilinear) geospatial data, we are working on an xarray plugin "xarray-regrid".

As accessor we use regrid.

Currently we have a few methods implemented:

• linear, nearest neighbor, cubic (simple wrappers around xr.interp)
• conservative (custom method, based on a comment by Stephan Hoyer here)
  • It is compatible with dask.distributed and outperforms xESMF, on par with CDO.
• A "most common value" zonal statistics method.

On our repository we also have some notebooks benchmarking the results against xESMF and CDO.

I read I should let you know about using the accessor, hope you are interested!