[MITgcm-support] T & rho!

[Van Thinh Nguyen]

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
------------------------------------------------------