[MITgcm-support] T & rho!

[Martin Losch]

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