[MITgcm-devel] description of seaice parameters

[Dimitris Menemenlis]

Hi Martin,

> I'll try to add documentation for these parameters: 
> SEAICE_sensHeat/latentWater/latentIce to SEAICE_PARAMS.h and 
> seaice_readparms.F.

Thanks.

> pstar = SEAICE_strength=2.75e4 is believed to be very high by our ice
> modelers (Ralph Timmermann in particular). He said that he wouldn't
> question anything between 1.5e4 and 2e4 (or 2.5e4), based on his 
> experience. Also, this value in particular seems to be very different
> between Artic and Antartic applications. Also, a freezing point of 
> Tf=-1.86 seems to make more sense. In the end, having Tf=f(S,p) does
> make most sense.

My suggestion is that information about range believed to be reasonable 
be added to SEAICE_PARAMS.h but that we leave default values unchanged 
for time being, until we have a bit more state estimation experience 
with pkg/seaice.  But happy to be overriden if you feel strongly about 
changing these values.

> HO (or H0=h naught, as they say here) is effectively 0.5m in
> pkg/seaice, Ralph uses 1.0m (derived from personal experience).

Same comment as above.  But should we remove the "cancelling" factors of 
two from the code?

Jinlun, what is your opinion?

Dimitris

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