[MITgcm-devel] heff_max in pkg/seaice/seaice_growth.F

[Jinlun Zhang]

It is not really a good idea to put a cap on ice thickness. Something 
else is wrong if a cap has been use to prevent blowup.
Jinlun

Dimitris Menemenlis wrote:

> I plan to make following two changes to seaice_growth:
>
> 1. Comment the following line.
>
>> C     Capping the actual ice thickness effectively enforces a
>> C     minimum of heat flux through the ice and helps getting rid of 
>> C     very thick ice.
>>           HICE(I,J)    = MIN(HICE(I,J),9.0 _d +00)
>
> An ran a test and she finds that this limit actually causes thicker 
> rather than thinner ice, i.e., the limit on HICE, allows more cooling 
> of the ocean water to take place hence growing thicker ice.
>
> 2. Uncomment following line.
>
>> #ifdef DO_WE_NEED_THIS
>> c          HEFF(I,J,1,bi,bj)=MIN(MAX_HEFF,HEFF(I,J,1,bi,bj))
>> #endif
>
> Without a cap on HEFF_MAX, the CS510 configuration crashes.
>
> OK?
>
> Dimitris