Function compute_flow_accumulation 

pub fn compute_flow_accumulation(
    flow_dir: &[u8],
    width: usize,
    height: usize,
) -> Vec<u32>

Compute flow accumulation from flow direction map