add domain aggregation
Refactored existing aggregation function to be used for time aggregation (returns a mesh) and for domain aggregation (returns numpy array). Also several small plotting fixes
Refactored existing aggregation function to be used for time aggregation (returns a mesh) and for domain aggregation (returns numpy array). Also several small plotting fixes