[MITgcm-devel] freezing in the interior

[Dimitris Menemenlis]

Jean-Michel,
I was thinking of adding to main code as this could also
be useful for pkg/icefront, and in some circumstances
for pkg/seaice.

I would follow the
allowFreezing/freeze_surface.F
useOldFreezing/freeze.F
template and add a flag and subroutine in main code say:
allowFreezingInterior/freeze_interior.F

Is this OK?

On Tue, 2011-04-19 at 05:48 -0700, Jean-Michel Campin wrote:
> Hi Dimitris,
> 
> Would this be part of the pkg/shelfice ? I know there is
> a depth dependent freezing temperature used in this pkg.
> 
> Cheers,
> Jean-Michel
> 
> On Mon, Apr 18, 2011 at 05:32:24PM -0700, Menemenlis, Dimitris (3248) wrote:
> > MITgcm ice czars, in configurations with ice shelf cavities we often end up with temperatures below the in situ freezing point at mid depth.
> > 
> > I am thinking of adding a subroutine similar to freeze.F but that will move the excess negative heat from the ocean interior to the surface level (or to the level right below the ice shelf) in order to conserve heat in the model.
> > 
> > Has anyone else given some thought to, or already written some code to deal with, this mid-depth freezing problem?