[MITgcm-devel] freezing in the interior

[Menemenlis, Dimitris (3248)]

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?

Dimitris