<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Hi CJ,<div class=""><br class=""></div><div class="">the composition of the diagnostics PHIHYD (totPhiHyd in the code) depends on the model configuration. For ECCO, I am guessing that the non-linear free surface with r-star coordinates is turned on and then totPhiHyd does not include the eta-contribution explicitly, but only impliclty via the variable r-star coordinate (this is more accurate, because then the free surface changes the layer thicknesses). Please have a look at model/src/diags_phi_hyd.F, starting line 107, which overwrites the default version at the top of the subroutine (which corresponds to your understanding of totPhiHyd). It’s probably going to be very tricky to reconstruct totPhiHyd exactly in this context (but I don’t know).</div><div class=""><br class=""></div><div class="">BTW, the conversion between pressure and geopotential is done with rhoConst (which happens to be 1029 in the ECCO context).</div><div class=""><br class=""></div><div class="">Martin<br class=""><div><br class=""><blockquote type="cite" class=""><div class="">On 19. Oct 2022, at 11:21, 秋菊陳 <<a href="mailto:ccj0315@gmail.com" class="">ccj0315@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><meta http-equiv="Content-Type" content="text/html; charset=utf-8" class=""><div dir="ltr" class=""><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Greetings,<br class=""></div><div dir="ltr" class=""><div class=""><br class=""></div><div class="">I would like to get some clarification about the variable PHIHYD, both physical meaning and to relate the value with other variables (SSHDYN, RHOAnoma and EXFpress). </div><div class="">(e.g. <a href="https://data.nas.nasa.gov/ecco/data.php?dir=/eccodata/llc_90/ECCOv4/Release4/interp_daily" class="">https://data.nas.nasa.gov/ecco/data.php?dir=/eccodata/llc_90/ECCOv4/Release4/interp_daily</a>. Yes, I am a ECCO data user.)</div><div class=""><br class=""></div><div class="">To my understanding, PHIHYD is the hydrostatic-pressure-anomaly/1029. The other way round, we could consider "PHIHYD*1029" as the sum of hydrostatic-pressure-anomaly ABOVE.</div><div class="">From the schematic figure-1 below, the hydrostatic-pressure-anomaly of the first grid PHIHYD(1)*1029 could be treated as sum of 3 components:</div><div class="">(1) atmospheric pressure (EXFpress), </div><div class="">(2) sea water pressure from sea surface (eta, SSHDYN) to z=0, i.e. g*( RHOAnoma+1029 )*SSHDYN </div><div class="">(3) sea water pressure from z=0 to first data grid z(1), i.e. g*RHOAnoma*abs(z(1))</div><div class=""><br class=""></div><div class="">However, I cannot get the result [sum of 3 components, fig2a green line] close to PHIHYD(1) [fig2a black line], unless ignoring component (2) [fig2a red line]. Did I misunderstand something?</div><div class="">Any help would be appreciated.<br class=""></div><div class=""><br class=""></div><div class="">Best,</div><div class=""><br class=""></div><div class="">CJ </div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">Figure-1 Figure-2</div><div class=""><div class=""><span id="cid:ii_l9f5cs1t2"><schematic.png></span> <span id="cid:ii_l9fcr3ol3"><check_PHIHYD1_nasnasa.png></span></div><div class=""></div></div></div></div>
</div></div>
_______________________________________________<br class="">MITgcm-support mailing list<br class=""><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br class=""></div></blockquote></div><br class=""></div></body></html>