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

Chris Hill cnh at mit.edu
Mon Jul 16 07:38:57 EDT 2007


One request (since JMC is in the Norwegian mountains this week).
Could someone please change "DO_WE_NEED_THIS" to a phrase that explains 
the role of the option e.g. "SEAICE_CAP_HEFF" or something!

Thanks,

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




More information about the MITgcm-devel mailing list