[MITgcm-devel] more seaice

[Dimitris Menemenlis]

> independently from all of the previous emails about seaice I would like 
> to make a suggestion for the seaice-package: I would like to introduce 
> masks (maskC) into the routine GROWTH for:
> ---> surfaceTendencyTice (line 139), because of neatness
> ---> theta (after line 139), because of the substraction of two "large" 
> numbers to prevent non-zero values over land
> ---> EmPmR (line 287), because the nonzero emp over land will lead to 
> nonzero sealevel over land under certain circumstances (e.g. if 
> convertFW2Salt = 35, or useRealFreshWaterFlux, or other up-to-now 
> unphysical situtations). I am not sure whether this is a real problem, 
> but it affects the values of cg2d_b, which are not masked. Over time, 
> however, Eta can become very large over land (60m after 10yrs, in my 
> case). I don't know if it makes sense to have a nonzero EmPmR over 
> land? Certainly not if it's computed from a SEA ice model (o:
> (---> Qnet is another candidate, but that has not lead to any visible 
> problems.)

OK, for all of the above.  Do you want me to have a go at
them, or are you already working on these?

Forgot to mention, NPSEUDO is only used for ADI, and as I said
ADI does not work properly.  Jinlun, should we just get rid of
it from the main branch of the code to avoid confusion?
It will still be there in the CVS archive.

One comparison/report (paper?) that would be very useful would
be to compare the thice and seaice (with and without dynamics)
packages.  I did some very preliminary tests some time ago:
http://ecco.jpl.nasa.gov/sea_ice/cube_sphere/32x32ice/
Figures are I think after a 5-year spin-up.
Hibler refers to pkg/seaice and Winton to pkg/thice
But it would be good to do something a bit more systematic
and complete.

D.

-- 
Dimitris Menemenlis <menemenlis at jpl.nasa.gov>
Jet Propulsion Lab, California Insitute of Technology
MS 300-323, 4800 Oak Grove Dr, Pasadena CA 91109-8099
tel: 818-354-1656;  fax: 818-393-6720