[MITgcm-support] C grid and pressure

[Jody Klymak]

Hi,


On Jan 7, 2012, at  8:37 AM, Jean-Michel Campin wrote:

>> From those output files, the pressure P can be obtained as:
> P(i,j,k)/rhoConst = -g*rC(k) + PH(i,j,k) + PNH(i,j,k) 
> 
> And PHL would be the equivalent of PH (hydrostatic pressure potential anomaly)
> but at the bottom, so that the hydrostatic pressure at the bottom Pb
> would be: Pb(i,j)/rhoConst = g*Depth(i,j) + PHL(i,j)

I wonder if this could be documented better somewhere? I too spent a lot of time tracking this down, and still wasn't sure if it was correct. Please forgive me, if it is documented, and I didn't know where to look.

Perhaps in 3.5.1.1:

PH.0000nIter and PNH.0000nIter - ocean: hydrostatic and non-hydrostatic pressure anomallies P(i,j,k) = -g*rC(k)*rhoConst + rhoConst*PH(i,j,k)+rhoConst*PNH(i,j,k).


Thanks,  Jody



--
Jody Klymak    
http://web.uvic.ca/~jklymak/