[MITgcm-support] Buoyancy anomaly

[Martin Losch]

Hi Jody,

I guess I was suggesting that the appropriate diagnostics to look at  
stability is DRHODR, rather than diagnosing the density anomaly and  
infering drho/dr from that. As far as i know it's (the code version  
of) DRHODR that is relevant to vertical acceleration/mixing. The  
vertical gradient of the anomaly is tricky because it also depends on  
the reference stratification tRef/sRef.

I think you are right and I have confused one with the other in my  
previous comment.
and I was under the impresstion that Lee was plotting the density  
anomaly (not the gradient).

Martin

On 22 Mar 2007, at 23:03, Jody Klymak wrote:

> Hi Martin,
>
> On Mar 22, 2007, at 1:45 AM, Martin Losch wrote:
>
>>
>> 1. I am not sure, if you are looking at the appropriate  
>> diagnostics. the density anomaly for a linear equation of state is
>> rho-rhoConst = rhoNil*tAlpha*(theta(i,j,k,bi,bj)-tRef(k))
>> Your stratification is the vertical gradient of this.  I believe  
>> the diagnostic for  this is:
>>       diagName  = 'DRHODR  '
>>       diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)'
>>       diagUnits = 'kg/m^4
>
> I'm sorry, but I'm not sure what you are saying here.  Are you  
> saying the diagnostic that Lee was plotting was the gradient of the  
> density anomaly \rho'?
>
> \rho' = \rho-\rho0-rho_ref(z)?
>
> i.e. d\rho/dz -d\rho_ref/dz?
>
> That would jibe with your equation above.  Otherwise, I'm a little  
> confused.
>
> Thanks a lot for any clarification,
>
> Cheers, Jody
>
>
> --
> Jody Klymak
> http://web.uvic.ca/~jklymak/
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support