safer way to find close incices

geoscixyz · Nov 1, 2023 · 5e48e6e · 5e48e6e
1 parent 6419180
commit 5e48e6e
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/geoscilabs/dcip/DCWidgetResLayer2_5D.py b/geoscilabs/dcip/DCWidgetResLayer2_5D.py
@@ -435,8 +435,8 @@ def PLOT(
         VNprim = 0.0
 
     else:
-        MInd = np.where(xSurface == M)
-        NInd = np.where(xSurface == N)
+        MInd = np.abs(xSurface - M).argmin(axis=0)
+        NInd = np.abs(xSurface - N).argmin(axis=0)
 
         VM = phiTotalSurface[MInd[0]]
         VN = phiTotalSurface[NInd[0]]