dhdt currently masks out elevations < 0 m

[whyjz]

https://github.com/whyjz/CARST/blob/master/carst/libdhdt.py#L383

Not an issue for most areas, but if users are stacking files with valid pixel values < 0, these values will be lost in the stacked file.