[MITgcm-devel] seaice+thSIce changes

[Jean-Michel Campin]

Martin,

I started to make changes in seaice advection (pkg/seaice)
of thSIce variables, but realized that more EXCH calls
are needed to implement correctly the advection of Qice
(instead of Qice*hEff as it is now).
I left this piece of code within #ifdef THSICE_NEW_ADVECT
(undef in pkg/seaice/seaice_advdiff.F).

Also, I am going to try to enforce a minimum ice volume (or hEff)
together with a minimum ice Area after the advection.
(lines > C jmc: new version: > C jmc: comment out the old version:
in pkg/seaice/seaice_advdiff.F ); it might prevent
some random blow-up (as the one you reported earlier ?).

I also moved around the diffusion (mainly because it's
a pain to change 3 S/R, seaice_calc_rhs, seaice_advdiff &
seaice_advection instead of only 2, every time I had to add an
argument to SEAICE_ADVECTION), disable it for hSnow & Qice
(since we still need to use the 1rst Order upwind scheme ...),
but should be OK (I did short tests, both with & without thSIce,
and nothing has changed).

Jean-Michel