Request to add function regarding HRU creation through pySWATPlus

[debpal]

pySWATPlus is a excellent module to interact without the need of QSWAT+. For my study area, watershed delineation can't be handled by QSWAT+ due to 900 million cells. I use the 'Use existing watershed' option in the SWAT+ and prepare the input files by pyflwdir module (https://github.com/Deltares/pyflwdir) since it can handle large raster and give fast result. For 900 million cells, it take only 20 minutes to prepare all files. The SWAT+ developer Chris George helped me to avoid Dinf flow direction calculation when using the 'Use existing watershed' option?

Can you please add a function of HRU creation so that I can create a automatic simulation workflow from start to end without QSWAT+ intervention? It would be helpful for me to do multi-objective optimization of drainage filling in the DEM. We would be happy to collaborate also.