[MITgcm-support] How to output b (buoyancy) using MITgcm?

[Martin Losch]

Brainf,

buoyancy is b=gravity*rho/rhoConst, correct?

 From what I see in the internal_wave example:
EOStype = 'linear',
talpha = 2.e-4,
sbeta = 0.,
so that rho = rho0*(1-talpha*(theta-theta_ref)); that is your density.

If you want to use the diagnostics package, use the variable RHOanom  
(density anomaly = rho-rhoConst, rhoConst is set in set_defaults.F or  
your "data" file), and you compute buoyancy like this b = gravity* 
(RHOanom+rhoConst)/rhoConst.

Martin
On 4 Dec 2006, at 02:23, Brainf wrote:

> Hi mitgcm,
>     I want to simulate internal waves using MITgcm. I am very  
> interested in what
> Sonya Legg have done like "Internal wave breaking at concave and  
> convex continental" research.
>     I found he plotted b(buoyancy) contour to see internal waves,  
> and said the contour interval is △b=9.81*10^(-6)m/s^2 ,but how to  
> output this?  The traditional output  files are  U V W T S and ETA;  
> buoyancy is not included. I checked available diagnostics, can't  
> find variable buoyancy also. Do I need to edit the mdsio package  
> output source code?
>     I start from  verification/internal_wave experiment,and have  
> change the resolution to 640x100.
>
>
> Best regards.
>
> Brainf
> brainf at 163.com
> 2006-12-04
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support