Heat and salt budget notebook

COSIMA / cosima-recipes

A cookbook of recipes (i.e., examples) for analysing ocean and sea ice model output

https://cosima-recipes.readthedocs.io

Apache License 2.0

46 stars 66 forks source link

Heat and salt budget notebook #192

Open adele-morrison opened 1 year ago

adele-morrison commented 1 year ago

Would it be useful to have a notebook showing how to close heat and salt budgets? We don’t have 3D interior heat / salt budget terms for most simulations, but we could at least show how to close interior heat / salt content changes with surface fluxes. This information can potentially be figured out from the Surface Watermass Transformation notebook, but probably that’s not obvious to all?

rmholmes commented 1 year ago

Yes! I can do this. I've got lots of stuff lying around on various aspects of this which I've been meaning to document more precisely. I can cover:

Closing globally-integrated budgets (i.e. what are the external sources, also see https://forum.access-hive.org.au/t/non-closure-of-heat-budget-in-access-cm2-esmf-cmip-output/244/2)
Closing vertically-integrated budgets
Closing full 3D budgets (including which terms are on in the different resolutions)
How the budget diagnostics relate to changes in the properties themselves (e.g. snapshots versus averages, accumulation of rounding errors etc).

Much of my code is in Matlab, but this is a good opportunity to python-ify it.

ChrisC28 commented 1 year ago

Is there any appetite to look at closing tracer budgets in MOM6 with the new vertical coordinates? While I'm aware of the GFDL example by Graeme McGilchrist, i think a ACCESS-OM3 version could be really useful.

adele-morrison commented 1 year ago

Yes that would be awesome. I'd imagine we'd need to rerun a MOM6 simulation and output extra diagnostics?

ChrisC28 commented 1 year ago

I can try to get a year or two of the EAC config run by then with the appropriate diagnostics added, but Setonix is still down for maintenance. Could we ask @angus-g or @AndyHoggANU very nicely to run another year of the global 0.1 degree MOM6 with the additional diagnostics?

rmholmes commented 1 year ago

I've made a start on this from the MOM5/ACCESS-OM2 perspective in this notebook: https://github.com/rmholmes/cosima-recipes/blob/heat_budget/DocumentedExamples/Closing_Heat_Budgets_in_ACCESS-OM2.ipynb I will continue!