Added PatchZonalStats for computing zonal statistics for patch-scale variables
from yearly RHESSys output