avoid extrapolation artefacts

GliderToolsCommunity · Apr 11, 2024 · 0a6a5b8 · 0a6a5b8
1 parent 4bcfee8
commit 0a6a5b8
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/glidertools/mapping.py b/glidertools/mapping.py
@@ -974,7 +974,9 @@ def grid_data(
     gridded = gridded.reindex(labels.astype(float))
 
     if interp_lim > 0:
-        gridded = gridded.interpolate(limit=interp_lim).bfill(limit=interp_lim)
+        gridded = gridded.interpolate(limit=interp_lim, limit_area="inside").bfill(
+            limit=interp_lim
+        )
 
     if not return_xarray:
         return gridded