[MITgcm-support] difference between PH and PHL

Jody Klymak jklymak at uvic.ca
Mon Aug 5 22:56:01 EDT 2013


On Aug 5, 2013, at  16:27 PM, Jean-Michel Campin <jmc at ocean.mit.edu> wrote:

> Hi Jody,
> 
> I will try to clarify:
> 1) for atmospheric setup (buoyancyRelation='ATMOSPHERIC', fluidIsAir=T),
> you are right, phiRef contains the geopotential height computed from
> tRef (sRef dependence is ignored, assuming dry air).
> 2) with oceanic set up in z type coordinate (z or z*, but not pressure)
> the reference state is homogenous rhoConst, and phiRef is
> just -gravity*rC or rF, so that Nicolas was right in this case.
>  P(i,j,k)/rhoConst = -g*rC(k) + PH(i,j,k) + PNH(i,j,k)
> or:
>  P(i,j,k)/rhoConst = PHrefC(k) + PH(i,j,k) + PNH(i,j,k)

OK, well, I guess the confusion is in find_rho.F, in which different equations of state appear to use different conventions?

calc_phi_hyd.F uses "rhoInSitu" which so far as I can tell is calculated in do_oceanic_phys.F L520 or so by a call to find_rho_2d.  For that function, for a linear equation of state the code is:

       refTemp=tRef(kRef)
       refSalt=sRef(kRef)

       dRho = rhoNil-rhoConst

       DO j=jMin,jMax
        DO i=iMin,iMax
         rhoLoc(i,j)=rhoNil*(
     &     sBeta*(sFld(i,j)-refSalt)
     &   -tAlpha*(tFld(i,j)-refTemp) )
     &        + dRho
        ENDDO
       ENDDO

So perhaps an inconsistency has crept in between the meaning for Linear Equation of state (which use, a lot), and the non-linear ones?  I note that "POLY3" at least has a different coding and indeed appears to add in rho_ref.  

Thanks,   Jody




> 3) with oceanic set up in pressure coordinate, it would be more logical
> to use tRef,sRef (similar to atmosphere in pressure coords) 
> but instead we continue to use rhoConst to define the reference state.
> see in model/src/set_ref_state.F, lines 255-272.
> 
> And the way the hydrostatic pressure anomaly is computed is consistent with
> the definition of PhiRef (i.e., uses tRef only if atmosphere in p,
> otherwise based on rhoConst).
> 
> I hope I did not make too much confusion.
> 
> Cheers,
> Jean-Michel
> 
> On Mon, Aug 05, 2013 at 01:42:42PM -0700, Jody Klymak wrote:
>> Hi All,
>> 
>> On Aug 5, 2013, at  13:24 PM, Nicolas Grisouard <nicolas.grisouard at stanford.edu> wrote:
>> 
>>> I'm going to quote Jean-Michel, who answered this question when I asked it a couple of years ago (Jan. 7 2012 to be specific):
>>> 
>>> "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) 
>> 
>> Except this isn't quite correct.  
>>> P(i,j,k)/rhoConst = -g*rC(k) + PH(i,j,k) + PNH(i,j,k) + PHref(k)
>> where PHref(k) is the hydrostatic pressure due to Tref and Sref.
>> see calc_phi_hyd.F and find_rho.F.  Unless I am grossly misunderstanding the code.
>> Thanks,   Jody
>> 
>> 
>> 
>>> 
>>> 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)
>>> 
>>> Cheers,
>>> Jean-Michel"
>>> 
>>> It's pretty useful, because you don't have to interpolate the hydrostatic pressure on the topography.
>>> Cheers,
>>> Nico.
>>> 
>>> On 05/08/13 06:12, Robert Nazarian wrote:
>>>> Hi,
>>>> 
>>>> I'm using the MITgcm internal wave code and analyzing the pressure perturbations. I noticed in the output directory that         there are PH and PHL files, which are both hydrostatic pressure files. Can someone elaborate on the difference between the two?
>>>> 
>>>> Thanks so much,
>>>> Rob
>>>> 
>>>> Robert Nazarian
>>>> Department of Atmospheric and Oceanic Sciences
>>>> Princeton University
>>>> rn2 at princeton.edu
>>>> 
>>>> 
>>>> _______________________________________________
>>>> MITgcm-support mailing list
>>>> MITgcm-support at mitgcm.org
>>>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>>> 
>>> 
>>> -- 
>>> Nicolas Grisouard - Postdoctoral Scholar
>>> Stanford School of Earth Sciences
>>> 473 via Ortega - room M11
>>> Stanford, CA 94305-4216
>>> Dept. fax: (+1) 650-498-5099 
>>> Office Skype ID: nicolas.grisouard
>>> http://nicolas.grisouard.free.fr
>>> http://www.mendeley.com/profiles/nicolas-grisouard/
>>> _______________________________________________
>>> MITgcm-support mailing list
>>> MITgcm-support at mitgcm.org
>>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>> 
>> --
>> Jody Klymak    
>> http://web.uvic.ca/~jklymak/
>> 
>> 
>> 
>> 
> 
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-support
> 
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support

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







More information about the MITgcm-support mailing list