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

[Martin Losch]

Hi,
I would put
>> HICE(I,J)    = MIN(HICE(I,J),MAX_HICE)
and let MAX_HICE default to zero, the other is OK with me.

Martin

On 16 Jul 2007, at 09:50, 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
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel