[MITgcm-devel] freezing in the interior

[Gael Forget]

Hi Dimitris,

I also had thought of doing such a hack, so I am glad you are doing it. 

One thought : given that you will take some negative heat from a typically thick layer at depth,
putting it all at once in the surface layer may not prove very numerically stable; another option 
would be to spread that negative heat over the above layers (looping from the bottom layer up?).

Cheers,
Gael





On Apr 19, 2011, at 9:18 AM, Dimitris Menemenlis wrote:

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