Estimating surface flow paths on a digital elevation model using a triangular facet network