Trouble with pit filling and slope applying

[rbunel]

import richdem as rd
import numpy as np
dem              = nd.array([[250,260,265,270],[240,200,200,210],[230,220,225,205]])
dem              = rd.rdarray(dem, no_data=-9999)
dem_filled       = rd.FillDepressions(dem, epsilon=False, in_place=False,topology='D4')
dem_flat_resolve = rd.ResolveFlats(dem_filled, in_place=False)

print(dem)
print(dem_filled)
print(dem_flat_resolve)

>>> print(dem)
[[250 260 265 270]
 [240 200 200 210]
 [230 220 225 205]]
>>> print(dem_filled)
[[250 260 265 270]
 [240 210 210 210]
 [230 220 225 205]]
>>> print(dem_flat_resolve)
[[250 260 265 270]
 [240 206 208 210]
 [230 220 225 205]]
>>> 

I realize that when I try to correct flat areas with the ResolveFlats function, my flat area turns back into depression. The slope gradient for the flat area is in the opposite direction to what might be expected.

Could you please enlighten me on that point?

Thanks in advance

[r-barnes]

@rbunel: This looks like a bug. I'll try to get around to fixing it soon.