[MITgcm-support] pressure field from PHIHYD

[Matthew Mazloff]

Hi Martha

See Montgomery potential section on page 35 of http://scor-int.org/Publications/Thermophysical_Manual_09Jan09.pdf <http://scor-int.org/Publications/Thermophysical_Manual_09Jan09.pdf>

Its not completely straightforward, however, so if you want to compare to an observational estimate you may want to generate P, where P = -g*rC + PHIHYD, and then copy the method used by the estimate you are trying to reproduce….

Matt







> On May 4, 2016, at 8:10 AM, Martha Buckley <marthabuckley at gmail.com> wrote:
> 
> 
> I want to plot pressure on a given isopynal, and to be comparable to other papers (particularly observational estimates), I want to plot the full pressure field (with the depth dependence).   How do I do this using PHIHYD? I am using r* coordinates.  
> 
> Lets say for simplicity I have a profile of PHIHYD. Also rhoconst (constant density), RC (cell centered depths), and g (gravity).  
> 
> Panom=PHIHYD*rhoconst;
> P=Panom+(g*rhoconst*-1*mygrid.RC)
> 
> Is this correct?  Or does a term involving SSH needed to be added back on too?  I thought PHIHYD was the full dynamic part of hydrostatic pressure, but the manual indicates there is a surface term as well. 
> thanks,
> Martha
> 
> 
> thanks,
> Martha
> 
> -- 
> Martha W. Buckley
> marthabuckley at gmail.com <mailto:marthabuckley at gmail.com>
> http://sites.google.com/site/marthabuckley/home <http://sites.google.com/site/marthabuckley/home>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20160504/4ed7aa2b/attachment.htm>