[MITgcm-support] linear instability

[mlosch at awi-bremerhaven.de]

Hi Zoe,
as far as I know, this (Dirichlet) type of boundary condition for 
temperature is not availble for the top nor the bottom because it 
is not normally used in ocean modeling (it is for open lateral 
boundaries). However, you could hack restoring boundary 
conditions in model/src/external_forcing.F for the bottom (for the 
top they are already coded), so that the temperature is nearly 
your desired values. The alternative is to find all the places in the 
code where the open boundary conditions are applied (call 
obcs_apply_ts and call obcs_apply_t) and apply the Dirichlet 
BCs at the surface and the bottom there.

Martin

Martin Losch
Alfred Wegener Institute 
Postfach 120161, 27515 Bremerhaven, Germany; 
Tel./Fax: ++49(0471)4831-1872/1797



----- Original Message -----
From: zlr at noc.soton.ac.uk
Date: Thursday, March 9, 2006 7:24 pm
Subject: [MITgcm-support] linear instability

> Hello!
> 
> I have been using MIT as a comparison for a non-hydro. finite 
> element adaptive 
> mesh ocean model.  I am trying to set the top and bottom bcs 
on 
> temp to 
> preserve the initial values eg 0C at the top and 50C at the 
bottom. 
> Any advice 
> would be greatly appreciated.
> 
> Thanks,
> 
> Zoe Roberts (PhD student)
> National Oceanography Centre, Southampton
> UK
> 
> --------------------------------------------------------------------
> ---
> National Oceanography Centre, Southampton :: 
> http://www.noc.soton.ac.uk
> 
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>