[MITgcm-devel] experiment to check conservation in seaice

[Martin Losch]

Hi Jean-Michel,

I finally set up a simple experiment based on the verification/global_ocean.cs32x15 with seaice and seaice thermodynamics (no winton). I failed miserably in trying to check the salt and heat conservation even for the ocean. The only thing that i managed to get right is volume conservation in the ocean, which is almost trivial of course. The seaice model has a drift of about 4mm/y as expected from earlier experiment.

I put the configuration on faulks:~mlosch/cs32conserve.tgz (just the code and input directories).

In the input directory there's a matlab script that tries to make use of the diagnostics output in netcdf (check_conserve.m). Please teach me how to close the heat and salt budgets.

Martin