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

[Martin Losch]

Hi Chris,
good point, Dimitris, could you do this while you are at changing  
seaice_growth.F?
But beware, this flag is used twice. The first instance refers to  
Hibler's original "ridging scheme" c=min(c,1); and we definitely  
NEED_THIS. Only the second instance refers to capping the effective  
thickness.

M.
> C NOW TRUNCATE AREA
> #ifdef DO_WE_NEED_THIS
>           AREA(I,J,1,bi,bj)=MIN(ONE,AREA(I,J,1,bi,bj))
>          ENDDO
>         ENDDO
> #ifdef ALLOW_AUTODIFF_TAMC
> CADJ STORE area(:,:,:,bi,bj) = comlev1_bibj,
> CADJ &                         key = iicekey, byte = isbyte
> CADJ STORE hsnow(:,:,bi,bj) = comlev1_bibj,
> CADJ &                         key = iicekey, byte = isbyte
> #endif /* ALLOW_AUTODIFF_TAMC */
>         DO J=1,sNy
>          DO I=1,sNx
>           AREA(I,J,1,bi,bj) = MAX(ZERO,AREA(I,J,1,bi,bj))
>           HSNOW(I,J,bi,bj)  = MAX(ZERO,HSNOW(I,J,bi,bj))
> #endif
>           AREA(I,J,1,bi,bj) = AREA(I,J,1,bi,bj)*HEFFM(I,J,bi,bj)
>           HEFF(I,J,1,bi,bj) = HEFF(I,J,1,bi,bj)*HEFFM(I,J,bi,bj)
> #ifdef DO_WE_NEED_THIS
> c          HEFF(I,J,1,bi,bj)=MIN(MAX_HEFF,HEFF(I,J,1,bi,bj))
> #endif

On 16 Jul 2007, at 13:38, Chris Hill wrote:

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