MOM6 Vertical Layer Plotting

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

MOM6 Vertical Layer Plotting #222

Open claireyung opened 1 year ago

claireyung commented 1 year ago

I'd like to make a general MOM6 notebook that can plot layers (in depth space) of the grid the model ran on e.g HYCOM (like a general version of https://github.com/COSIMA/mom6-diagnostics/blob/main/PanAntarctic/HYCOM1_Transects.ipynb with correct interpolation) or for rho or sigma layers. However, this requires e output (h output is also useful but doesn't work with with ice-shelf cavities) in the native coordinate output, i.e. not remapped to z or rho, which I couldn't find in any currently public COSIMA MOM6 output.

navidcy commented 1 year ago

Say no more!

navidcy commented 1 year ago

Tag some people here who you think they might know where you can get the output w the diags you need?

adele-morrison commented 1 year ago

We want a config with interesting vertical coords (e.g. hybrid or rho not z*). Possibilities I can think of:

OM4. What vertical coord is this? Claire said it doesn't have the right diags needed.
panant-01-hycom1-v13, but we'd need to run it a bit longer to output the necessary diags.
@pedrocol's Amery config. What vertical coord does this use? SIGMA_SHELF_ZSTAR? Does it have requried diags?
One of the anu-tub cases until we get output from a more realistic case? @angus-g @AndyHoggANU?

pedrocol commented 1 year ago

Amery config uses z*, which implies vanishing layers at the top of the cavity, and the topography, as usual. The diag_table is set to save the "e" field both in the model grid and the interpolated grid.

navidcy commented 1 year ago

@claireyung this is still pending... is it a good idea to pursue it during hackathon v3.0?