[MITgcm-devel] more seaice

[Jinlun Zhang]

Dimitris, looks like it's better to remove it and avoid confusion.
Jinlun


Dimitris Menemenlis wrote:

> > 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.