Update geospatial_utilities.py

ImperialCollegeLondon · Feb 25, 2024 · 682bbe4 · 682bbe4
1 parent 2bfedaf
commit 682bbe4
Showing 1 changed file with 6 additions and 1 deletion.
diff --git a/swmmanywhere/geospatial_utilities.py b/swmmanywhere/geospatial_utilities.py
@@ -422,7 +422,12 @@ def delineate_catchment(grid: pysheds.sgrid.sGrid,
                                 y=y_snap, 
                                 fdir=flow_dir, 
                                 dirmap=dirmap, 
-                                xytype='coordinate')
+                                xytype='coordinate',
+                                algorithm = 'recursive'
+                                )
+        # n.b. recursive algorithm is not recommended, but crashes with a seg 
+        # fault occasionally otherwise.
+
         grid_.clip_to(catch)
 
         # Polygonize the catchment