Continued development of computeReachSTATS routine
added sub-routine for computing reach profile statistics
from a given subcatchment and reference reach shapefile