Consequent domain and time masking in meshseries aggregate functions

Added masking in both dimension in aggregate functions (except aggregate_over_time which returns the whole Mesh, so only time_mask there).

Made use of data instead of values for efficient reading of the data.