The assignment of drainage direction over flat surfaces in raster digital elevation models