It would be nice to have global sums in netcdf output: e.g. for each hemisphere, save total area, volume etc as a 1d variable straight out of the model.