[MITgcm-devel] upcoming changes in seaice_growth.F

Sun May 29 02:14:06 EDT 2011

Gael, mea culpa.  I did not follow recent pkg/seaice changes very closely and I have not yet tried the EVOLUTION branch.

My understanding from latest devel list exchange is that areaMin will only be used for regularizing thickness in heat budget computations, i.e., same as in original Jinlun/Hibler code.

Do we need to hang on to the ALLOW_PRECLUDE_INFINITESIMAL_AREA bit of code?

Thanks for merging/updating the sea ice code.

Dimitris Menemenlis

On May 27, 2011, at 6:34 AM, Gael Forget wrote:

Hi Dimitris,

we are not going to remove areaMin because it is needed in the LEGACY branch.

rev.87 is when I added the merged branch (which I refer to as EVOLUTION) including
the set of pathological cases that I thought required some treatment. The one you are referring
to (area=max(area,areaMin) if ice is present) I included based on Martin's comment
C The default of A22 = 0.15 is a common threshold for defining
C the ice edge. This ice concentration usually does not occur
C due to thermodynamics but due to advection.

After I completed this phase of the merging in October, I sent a detailed email to the devel list, entitled

[MITgcm-devel] wrapping up seaice_growth.F  october 2010 revision

where I was referring to this pathological cases part specifically, and suggesting comparison runs.
I don't recall getting much feedback. So, to be on the safe side, SEAICE_GROWTH_LEGACY has
remained the default since then. I am still hoping to get more feedback...

Have you ever tried the EVOLUTION branch? Did you notice any big difference compared
with LEGACY branch results? (I didn't at the time)

Cheers,
Gael

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-devel/attachments/20110528/149abac4/attachment.htm>