Would be great to have a function to combine routing products from two or more folders into one product. This would be useful in instances where you have two sections of catchments that you ultimately want in one model, but with two distinct outlets.
Something like:
basinmaker.basinmaker.delineate.Combine_Routing_Products(routing_product_folders, gis_platform='qgis')
routing_product_folders = ['folder1','folder2']
gis_platform='qgis'
Output is a single routing product determined from combining one or more routing products. Key challenge would be to confirm that the individual routing products have non-overlapping extents. Note that this can currently be done manually with GIS software, and subsequent operations seem to work ok.