Merge pull request #2159 from OceanParcels/minor_bug_xgrid_search

Fixing a small vectorization bug in xgrid search

Author: erikvansebille

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)