Using xgcm and axis properties to automatically determine indices offsets in Interpolators

[erikvansebille]

Now that we have xdcm logic in the grid class, we should use that to automatically determine which indices to use as corner-nodes in Interpolators.

See e.g. the code below
https://github.com/OceanParcels/Parcels/blob/cebda92bf6168b359a369b39b16f3e5cfa9cffdb/parcels/application_kernels/interpolation.py#L232-L258

Originally posted by @erikvansebille in #2152 (comment)