[MITgcm-support] phiHydLow

[Dimitris Menemenlis]

Does anyone have experience in using phiHydLow ?

Dimitris

----------  Forwarded Message  ----------

Subject: Re: [Fwd: Re: bottom pressure]
Date: Friday 13 August 2004 06:11 pm
From: Ichiro Fukumori <fukumori at jpl.nasa.gov>
To: Dimitris Menemenlis <menemenlis at jpl.nasa.gov>
Cc: Benny Cheng <bcheng at jpl.nasa.gov>, christopher Hill <cnh at mit.edu>

Dimitris,

phiHydLow seems to have a bug for it to be ocean bottom pressure.
Benny was able to reproduce phiHydLow offline from TS and ETA but
we found that that calculation had a bug; density is non-zero below
ocean bottom and pressure calculation sums up the contribution
from these non-zero densities below topography.  The problem is in
how values below ocean bottom are masked.  A similar bug may
be present in MITGCM itself.

-- Ichiro

>From: Dimitris Menemenlis <menemenlis at jpl.nasa.gov>
>Organization: JPL
>To: Benny Cheng <bcheng at jpl.nasa.gov>
>Subject: Re: bottom pressure
>Date: Fri, 6 Aug 2004 13:39:47 -0700
>User-Agent: KMail/1.6.2
>Cc: Ichiro Fukumori <fukumori at jpl.nasa.gov>, dimitri at pacific.jpl.nasa.gov
>References: <5.1.0.14.2.20040806101117.02854ec0 at pacific.jpl.nasa.gov>
>         <1091823500.12635.81.camel at lochness.jpl.nasa.gov>
>In-Reply-To: <1091823500.12635.81.camel at lochness.jpl.nasa.gov>
>MIME-Version: 1.0
>Content-Disposition: inline
>Content-Type: text/plain; charset="utf-8"
>Message-Id: <200408061339.47797.menemenlis at jpl.nasa.gov>
>X-UIDL: 1950fbdc7128ef9d4dfd71bc77ad4489
>Content-Transfer-Encoding: 7bit
>
>On Friday 06 August 2004 01:18 pm, Benny Cheng wrote:
> > how do you enable bottom pressure output on the MITgcm ?
>
>Benny, you need to save the variable
>
>C     phiHydLow  :: Phi-Hydrostatic at r-lower boundary
>C                  (bottom in z-coordinates, top in p-coordinates)
>
>or
>
>C     phiHydLowtave:: Hydrostatic (ocean) pressure / (atmos) geo- Potential
>C                     at the fixed boundary: (ocean) bottom pressure
>C                     (atmos) geo- Potential
>
>I had disabled the dump of phiHydLowtave in pkg/timeave using CPP flag
>HRCUBE.  So you can either reenable it there for time-averaged values, or
>you can dump instantaneous values of phiHydLow directly.
>
>D.
>
>--
>Dimitris Menemenlis <menemenlis at jpl.nasa.gov>
>Jet Propulsion Lab, California Institute of Technology
>MS 300-323, 4800 Oak Grove Dr, Pasadena CA 91109-8099
>tel: 818-354-1656;  fax: 818-393-6720

===============================================================
Ichiro Fukumori               | Jet Propulsion Laboratory
e-mail: fukumori at jpl.nasa.gov | Mail Stop 300-323
phone: +1 (818) 354-6965      | 4800 Oak Grove Drive
fax: +1 (818) 393-6720        | Pasadena, California 91109, U.S.A.
===============================================================



-------------------------------------------------------

-- 
Dimitris Menemenlis <menemenlis at jpl.nasa.gov>
Jet Propulsion Lab, California Institute of Technology
MS 300-323, 4800 Oak Grove Dr, Pasadena CA 91109-8099
tel: 818-354-1656;  fax: 818-393-6720