adding checks to index_search to validate bcoord between 0 and 1

Author: erikvansebille

parcels/_index_search.py

@@ -67,6 +67,8 @@ def _search_time_index(field: Field, time: datetime):
             if field.data.time.data[ti] != field.data.time.data[ti + 1]
             else 0
         )
+    if tau < 0 or tau > 1:  # TODO only for debugging; test can go?
+        raise ValueError(f"Time {time} is out of bounds for field time data {field.data.time.data}.")
     return tau, ti
 
 

parcels/xgrid.py

@@ -399,6 +399,8 @@ def _search_1d_array(
         return 0, 0.0
     i = np.argmin(arr <= x) - 1
     bcoord = (x - arr[i]) / (arr[i + 1] - arr[i])
+    if bcoord < 0 or bcoord > 1:  # TODO only for debugging; test can go?
+        raise ValueError(f"Position {x} is out of bounds for array {arr}.")
     return i, bcoord