[Mitgcm-support] Overiding dynamics to set temperature.

[mitgcm-support at dev.mitgcm.org]

Timour,

 Temp. fixing can be imposed in the same place as the
resetting of sub-freezing temps. Not sure what the dynamical
response will be, but the easiest thing is to try it.
 Sub-freezing temp. fixing code is

soyinka{p38_head}% find . -name '*.F' -exec grep -Hi 'tfreezing' {} \;
./model/src/freeze.F:      _RL Tfreezing
./model/src/freeze.F:      Tfreezing=-1.9
./model/src/freeze.F:         IF (gTNm1(i,j,k,bi,bj) .LT. Tfreezing)
THEN
./model/src/freeze.F:          gTNm1(i,j,k,bi,bj)=Tfreezing


soyinka{p38_head}% find . -name '*.F' -exec grep -Hi 'LL FREEZE' {} \;
./model/src/dynamics.F:            CALL FREEZE( bi, bj, iMin, iMax,
jMin, jMax, k, myThid )

soyinka{p38_head}% find . -name '*.F' -exec grep -Hi 'allowFreezing' {}
\;
./model/src/dynamics.F:         IF (allowFreezing) THEN
./model/src/ini_parms.F:     & allowFreezing, ivdc_kappa,
./model/src/set_defaults.F:      allowFreezing       = .FALSE.

 This code can be modified to impose fixed temps anywhere.
 You need allowFreezing=.TRUE. in data (PARM01).

Chris
 
Timour Radko wrote:
> 
> Hi, Chris,
> 
> Talked to Alistair yesterday.  He made a general suggestion
> that it would be easier and safer to play around with
> temperature rather than with the velocities.  So I am now
> inclined to just relax the temperature at the base (!)
> of the ekman layer to what I need.  Would you
> approve this approach?
> 
> I also asked Alistair whether I can just set T (at each
> time step) at the first 3 vertical layers to some specific
> T(x,y).  On this he responded that this will either work
> or not work.  Is it easy to implement something like
> that in the model?
> 
> Thanks,
> Timour.
> 
> Chris Hill wrote:
> 
> > Timor,
> >
> >  Whats the plan? Did AListair and Jean-Michel have any
> > suggestions?
> >
> > Chris