Index: wflow-py/Scripts/wflow_flood_lib.py
===================================================================
diff -u -rf6d346c8a0883f44e74b99d2dce672a3c881b30e -r34850193a8dc11a1d40383cc39212dfad7b371b2
--- wflow-py/Scripts/wflow_flood_lib.py	(.../wflow_flood_lib.py)	(revision f6d346c8a0883f44e74b99d2dce672a3c881b30e)
+++ wflow-py/Scripts/wflow_flood_lib.py	(.../wflow_flood_lib.py)	(revision 34850193a8dc11a1d40383cc39212dfad7b371b2)
@@ -572,3 +572,8 @@
 
     else:
         return x, y, data, fill_val
+
+def define_max_strahler(stream_file, logging=logging):
+    xax, yax, stream_data, fill_value = gdal_readmap(stream_file, 'GTiff', logging=logging)
+    return stream_data.max()
+