Add option to delineate subbasins using either streamorder of Pfafstetter thresholds

[markhegnauer]

In regions with no / or limited number of gauges, it's sometimes still useful to derive subbasins for analysis and recognizability by users / clients. In Pyflwdir there are already options to do this (https://deltares.github.io/pyflwdir/latest/basins.html), but these functions are not yet available through HydroMT.

I would suggest implementing this in a way that the user, through the inifile can set this up:

[setup_subbasin]
subbasin_method = gauges / pfafstetter / streamorder / minarea / ... (so either user specified location, or automatic generation)
subbasin_method_threshold = 5 (for the automatic methods)

[Starphire1]

I am also very interested in this catchment delineation option.