[MITgcm-support] T & rho!

Van Thinh Nguyen vtnguyen at moisie.math.uwaterloo.ca
Wed Sep 6 12:29:13 EDT 2006


Hi Martin & others,

Thanks Martin for your explain.
It's true, below 18.5 then N2=0 due to rho=const. This density profile 
looks like 2 layers, I would like to study solitary waves generating by 
tidal flows over a depression bump (tidal current u=1.5*sin(omega*t), 
omega=M2). The model has OBCs at west & east, here are some specific 
options what I used to setup to the model:

-Domain: delX=10000*20., delY=200., delZ=120*1. (m) (2D)
  Gausian bump at centre (after X axis) and h(x,z)=-80.-40./(cosh(x/L));

-OBCs at west & East, using obcs_cal.F to define the tidal current and
  OBCs for U,V W and Eta; using obcs_sponge.F to reduce reflecting waves
  from OBCs

-using OBCs, KPP packages

-applying viscC2Leith=4. and viscC2LeithD=4. due to viscAz=1.E-5 (or
  1.E-6)

I could reduce a huge noise, however a bit of noises and reflecting waves
from OBCs are still remained.

If someone may have any experience about such model, please help me an 
advice.

Thanks a lot

Vth

------------------------------------------------------

On Wed, 6 Sep 2006, Martin Losch wrote:

> Hi VanThinh,
>
> if hydrogThetaFile is not set, then tRef is used to initialize the 
> temperature.
> If sBeta=0 (no salinity), you then treat pot temperature as if it were 
> density.
> You can do as you described, but your density profile leads to N2 = 0 below 
> 18.5 meters, doesn't it? I suggest Alistair's approach, see
> http://forge.csail.mit.edu/pipermail/mitgcm-support/2004-June/002269.html
>
> Martin
>
> On Sep 5, 2006, at 9:23 PM, Van Thinh Nguyen wrote:
>
>> Hi Martin,
>> 
>> I looked in your discuss on T & rho in MITgcm-support archives (June 2004), 
>> but I am really not sure what I am using in my model:
>> 
>> I have a vertical density profile as:
>> 
>> rho(z)=0.001*(18.5-2.75*(1+tanh((z+10)/2)))
>> 
>> From exp1 using linear EOS, we have:
>>
>>    d(rho)/dz = - rho_0*tAlpha*d(T)/dz
>> 
>> Integrating this relation, I get T(z)= - rho(z)/(rho_0*tAlpha)
>> (here rho_0=1000kg/m3, tAlpha=2.E-4 and rho(z) given as above)
>> 
>> This T(z) can be used as tRef in data?
>> 
>> Thanks so much for any help,
>> 
>> Van Thinh
>> ------------------------------------------------------
>> 
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support



More information about the MITgcm-support mailing list