import topotoolbox
dem = topotoolbox.load_dem('bigtujunga')
fd = topotoolbox.FlowObject(dem)
print(fd.flowpathextract(12345)) # doctest: +SKIP