Open changliao1025 opened 2 years ago
Depression filling must be domain scale.
After compset_calculate_flow_accumulation
, there are possibly multiple watersheds.
Starting from this point, all the remaining algorithms should consider multiple watersheds.
It is also possible that some cells does not belong to any watershed, therefore it is desired to use the domain container so the program won't miss any cell.
Once all the cells are processed, we can then group them into watersheds, and output individual watershed as well, just like the pyflowline component.
Stream segment topology information maybe added later. And the program needs to sync both domain-based and watershed-based.
compset_find_index_by_cellid
can only be used at domain.
Some algorithms are better handled at watershed scale. However, the model is expected to run with multiple watersheds. As a result, some redesign may be needed.