[MITgcm-devel] sea ice diagnostics

[Martin Losch]

Hi Jean-Michel, and others,

we are working on mechanical energy diagnostics for the sea ice model (following Bouchat and Tremblay 2014).

While doing so I realized that maybe the diagnostics of seaice are not in the right place:
seaice_diagnostics_state is called from do_statevars_diags when “seqFlag=0”. this is probably OK for fields like HEFF, AREA, UICE, VICE, but not for derived quantities like SIzeta, SIeta, SIdelta, because they are all 0 in the first time step (after a pickup, say). With that all of SIsig1/2, etc are also not correct because they require zeta/eta from the previous iteration. The same is true for (some of) the energy diagnostics.

I am not sure what the correct philosophy is here. The model state variables (diagnostics_fill_state) are all computed at the beginning of the time step, but not always, so it’s a bit difficult for me to diagnose, when what type of variable should be passed to diagnostics_fill.

Could you give me some advice?

Martin