diff --git a/parcels/xgrid.py b/parcels/xgrid.py index 75140d0a9..6979d6f11 100644 --- a/parcels/xgrid.py +++ b/parcels/xgrid.py @@ -281,7 +281,7 @@ def search(self, z, y, x, ei=None): if "Z" in self.axes: zi, zeta = _search_1d_array(ds.depth.values, z) else: - zi, zeta = 0, 0.0 + zi, zeta = np.zeros(z.shape, dtype=int), np.zeros(z.shape, dtype=float) if ds.lon.ndim == 1: yi, eta = _search_1d_array(ds.lat.values, y)