Function flow_direction_d8 

pub fn flow_direction_d8(dem: &Dem) -> Vec<D8Code> 

Compute D8 flow direction for every grid cell.

Returns VecD8Codewith lengthrows * cols`. Flat or edge cells that have no steeper neighbour receive code 0.