[MITgcm-support] Stratification simulations

[Jordan Dawe]

David Hebert wrote:
> Hi all,
>
> I have run two identical simulations, one with a constant temperature 
> field, and one with temperature gradient specified via tRef array. The 
> simulations are run in non-hydrostatic mode. When I compute the 
> vertical velocity fields, and consequently vertical KE = 0.5 w^2, the 
> results are identical.
In your CPP_OPTIONS.h file, do you have this defined?

C o Include/exclude phi_hyd calculation code
#define INCLUDE_PHIHYD_CALCULATION_CODE

If INCLUDE_PHIHYD_CALCULATION_CODE is undefined, the model is compiled 
without hydrostatic pressure code, I think.  Tripped me up for a while once.

Jordan