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

[Dimitris Menemenlis]

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