segmentation error with fill-spill-merge on a small DEM

[jkingslake]

Hello!

I am getting an error that I do not understand running fill-spill-merge on a very small example DEM.

The DEM is here https://raw.githubusercontent.com/r-barnes/Barnes2019-DepressionHierarchy/master/test_cases/kerry_test3.dem. I saved it as downloaded.dem.

Richdem compiled without any errors (i followed the instructions in the README).

I ran the following

$ apps/rd_fill_spill_merge.exe downloaded.dem simple_test_2 -9 --swl 1 --save_dh my_dh

and got

m

Input DEM = downloaded.dem
m Output prefix = simple_test_2
m Surface water level = 1
m Surface water file =
m Ocean level = -9
p Reading topography...
m Data width = 10
m Data height = 10
m Data cells = 100
p Performing bucket fill...
p Getting depression hierarchy...
m Loading DH from my_dh
p Performing FillSpillMerge...
zsh: segmentation fault apps/rd_fill_spill_merge.exe downloaded.dem simple_test_2 -9 --swl 1 --save_d

I get the same behavior if I run it with fewer inputs:

apps/rd_fill_spill_merge.exe downloaded.dem simple_test_2 -9

or

apps/rd_fill_spill_merge.exe downloaded.dem simple_test_2 -9 --swl 1

The wtd raster is created but it is empty and the -hydrologic-surface-height raster is missing.

Any help would be very much appreciated!

cc @KCallaghan