[MITgcm-devel] freezing in the interior

[Dimitris Menemenlis]

Martin, that's a good question.
I had not been thinking that far ahead.
Maybe for time being I exclude when allow_autodiff_tamc
is defined until we gain some experience with forward code?

Dimitris

On Wed, 2011-04-20 at 08:48 -0700, Martin Losch wrote:
> Hi Dimitris,
> 
> is this supposed to work with the adjoint? if so you need to update
> model_ad_diff.list and probably add some store directive to
> do_ocean_physics.F. If you do not want it to work with the adjoint
> then you should exclude the call when allow_autodiff_tamc is defined.
> 
> Martin
> On Apr 20, 2011, at 3:59 PM, Dimitris Menemenlis wrote:
> 
> > Gael, Steph, JM, I have checked in a first draft of
> freeze_interior.F
> > subroutine.  It is enabled with flag allowInteriorFreezing in
> PARM01.
> > Right now all the negative heat comes to the surface level, where it
> > would contribute to ice shelf or sea ice thickness or to cooling the
> > surface level.
> > 
> > I like Gael's suggestion of distributing the negative heat
> throughout
> > the water column as the ice rises.  It's a bit like the inverse of
> > a salt plume.  The question is how much heat such a process would
> > give up on the way up.  Would it use up all the available heat in
> > each level down to the freezing point or would it mostly rise to the
> > surface at the sub-grid scale?
> > 
> > Before coding anything more complicated, I want to experiment with
> > this basic first version in our setups, see what happens and
> > correct any obvious bugs.
> 
>