subcatch function

Identify a hydrological catchment based on a single direction flow algorithm