Add point_in_cell method to query call

Author: fluidnumerics-joe

tests/v4/test_spatialhash.py

@@ -1,6 +1,7 @@
 import numpy as np
 
 from parcels._datasets.structured.generic import datasets
+from parcels._index_search import curvilinear_point_in_cell
 from parcels.xgrid import XGrid
 
 
@@ -15,7 +16,7 @@ def test_invalid_positions():
     ds = datasets["2d_left_rotated"]
     grid = XGrid.from_dataset(ds)
 
-    j, i = grid.get_spatial_hash().query([np.nan, np.inf], [np.nan, np.inf])
+    j, i = grid.get_spatial_hash().query([np.nan, np.inf], [np.nan, np.inf], curvilinear_point_in_cell)
     assert np.all(j == -1)
     assert np.all(i == -1)
 
@@ -27,7 +28,7 @@ def test_mixed_positions():
     lon = grid.lon.mean()
     y = [lat, np.nan]
     x = [lon, np.nan]
-    j, i = grid.get_spatial_hash().query(y, x)
+    j, i = grid.get_spatial_hash().query(y, x, curvilinear_point_in_cell)
     assert j[0] == 29  # Actual value for 2d_left_rotated center
     assert i[0] == 14  # Actual value for 2d_left_rotated center
     assert j[1] == -1