[MITgcm-devel] recently added SHELFICEMassStepping in pkg/shelfice

[Jean-Michel Campin]

Hi Dan and Martin,

regarding the newly added option "SHELFICEMassStepping" 
and new S/R SHELFICE_STEP_ICEMASS:

I realized that the units of forcing field array 
 "shelfIceMassDynTendency" are not documented (in SHELFICE.h, partly my fault), 
and would like to change it from a thickness tendency contribution
to a mass tendency per unit area [kg/m^2/s], and taking
this opportunity to remove the hard coded "910." number
(which represents ice-shelf density, but is inconsitent with 
pkg/shelfice parameter "rhoShelfIce").
Dan, are you OK with this modification ?

An other minor thing: I am thinking of changing few field arrays
in pkg/shelfice (including shelficeMass, which can now be stepped forward)
from _RS to _RL ; this should have very little consequences since
almost nobody uses RS=real*4. I think that by upgrading few arrays
there from _RS to _RL I can get 1 cg2d matching digits (currently: zero digits)
when I run testreport with -r4 option, see e.g.:
http://mitgcm.org/testing/results/2015_01/tr_danton_20150103_2/summary.txt

Cheers,
Jean-Michel