[MITgcm-devel] more seaice

[Martin Losch]

Hi,
On Thursday, April 29, 2004, at 06:28 PM, 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?

I will check in masking theta, surfaceTendencyTice and EmPmR. Qnet/Qsw, 
wouldn't make any difference, so I'll leave it for now.

The timestep issue, I am not that clear about. I have made all changes 
in growth, but I think all deltaT should replaced by either deltaTmom, 
deltaTfreesurf or deltaTtracers. What do the code czars think?
(I won't have time to do this until next week, anyway)

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

What does npseudo do for lsr, because I noticed that if npseudo > 1, 
lsr will be called npseudo times, and the "error"-print out in stdout 
reduces with each call.

Martin

> 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
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://dev.mitgcm.org/mailman/listinfo/mitgcm-devel