import topotoolbox
import matplotlib.pyplot as plt
dem = topotoolbox.load_dem('bigtujunga')
fd = topotoolbox.FlowObject(dem)
_= fd.ezgetnal(dem).plot()
plt.show()