[MITgcm-support] pressure from MITgcm

[Patrick Heimbach]

Sergey,

your calculation using PH=phiHydC (!)
p (PH) = rhoConst*PH(k) + g*rhoConst*D(k)

would be missing the part from the surface elevation
rho0*g*eta

This part is included in the field totPhiHyd,
but not in phiHydC (which is what you mentioned).
If that's the case, you would be missing all
the Eta variability which would explain the difference.
But maybe you meant totPhiHyd?

For details, see S/R diags_phi_hyd.F
in model/src/

-Patrick



On Fri, 2005-11-18 at 16:23, Sergey Vinogradov wrote:
> I had no clue about time-varying thickness of layers, and so that, and 
> the Pbar variability should have gone into my pressure-from-density 
> computation.
> Thanks, and have a good weekend!
> Regards
> Sergey
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
-- 
--------------------------------------------------------
Patrick Heimbach   Massachusetts Institute of Technology
FON: +1/617/253-5259                  EAPS, Room 54-1518
FAX: +1/617/253-4464             77 Massachusetts Avenue
mailto:heimbach at mit.edu               Cambridge MA 02139
http://www.mit.edu/~heimbach/                        USA