aggregation_diagnostics function

Compute multiple aggregation diagnostics for a set of variables