Issue with breaching and resolving flats

[TheMalekCH]

Hello there

For my bachelor thesis I am comparing preprocessing tools. I also want to include RichDEM but have an issue, especially with the FlatResolve algorithm, as some others previously had.

Because by resolving flats new depressions are produced and while I was able to use epsilon filling instead for the fill algorithm there is no alternative with the breaching algorithm. I also tried only breaching, but then discovered that all breached channels are the same height and therefore again a flat.

Is there any workaround within RichDEM to either resolve flats without producing new depressions or alternatively breach without producing new flats?

Thanks in advance and kind regards