SumUpstreamArea function

Calculate upstream area sums