<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,<div class=""><br class=""></div><div class="">there are couple of constant “rho”s in the code, e.g. rhoConst, rhoNil, rhoConstFresh, which all have a default, which can be changed. I usually check, which value the model actually uses (because I always forget, how the defaults are determined), in the stanard output (or in an MPI run in STDOUT.0000).</div><div class=""><br class=""></div><div class="">For the EOS calculation, checkout how “find_rho_2d” (in find_rho.F) does it: it uses “rhoNil” to scale the density (this should be “\rho_c”, while rhoConst = “\rho_0”). I am not sure if the 1D field rhoRef is used in this context.</div><div class=""><br class=""></div><div class="">Hope this helps,</div><div class=""><br class=""></div><div class="">Martin<br class=""><div><br class=""><blockquote type="cite" class=""><div class="">On 7. Oct 2022, at 08:39, Mike Fan <<a href="mailto:ifanliming@outlook.com" class="">ifanliming@outlook.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><meta charset="UTF-8" class=""><div class="WordSection1" style="page: WordSection1; caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;"><div style="margin: 0cm; text-align: justify; font-size: 10.5pt; font-family: DengXian;" class=""><span lang="EN-US" style="font-size: 12pt;" class="">Hello all,<o:p class=""></o:p></span></div><div style="margin: 0cm; text-align: justify; font-size: 10.5pt; font-family: DengXian;" class=""><span lang="EN-US" style="font-size: 12pt;" class=""><o:p class=""> </o:p></span></div><div style="margin: 0cm; text-align: justify; font-size: 10.5pt; font-family: DengXian;" class=""><span lang="EN-US" style="font-size: 12pt;" class="">I want to use the model output (temp & salt & eta) to calculate the phiHyd, Through reading the guidebook, I learned that phiHyd is called hydrostatic pressure potential anomaly, using the equation like<span class="Apple-converted-space"> </span></span><span lang="EN-US" class=""><a href="https://mitgcm.readthedocs.io/en/latest/examples/baroclinic_gyre/baroclinic_gyre.html#output-files" style="color: blue; text-decoration: underline;" class="">4.2. Baroclinic Ocean Gyre<span class="Apple-converted-space"> </span><span lang="EN-US" class=""><span lang="EN-US" class="">― MITgcm checkpoint68k-6-g80404458a documentation</span></span></a>.</span></div><div style="margin: 0cm; text-align: justify; font-size: 10.5pt; font-family: DengXian;" class=""><span lang="EN-US" class=""><o:p class=""> </o:p></span></div><div style="margin: 0cm; text-align: justify; font-size: 10.5pt; font-family: DengXian;" class=""><span lang="EN-US" style="font-size: 12pt;" class="">In my model, setting a linaer EOS,<span class="Apple-converted-space"> </span></span><span lang="EN-US" style="font-family: Lato, sans-serif; background-color: rgb(243, 246, 246);" class="">implicitFreeSurface=’true’</span><span lang="EN-US" style="font-size: 12pt;" class="">, and rhoConst=999.8. Here is my question, there is a rho_0, how the model calculates it?<span class="Apple-converted-space"> </span><o:p class=""></o:p></span></div><div style="margin: 0cm; text-align: justify; font-size: 10.5pt; font-family: DengXian;" class=""><span lang="EN-US" style="font-size: 12pt;" class=""><o:p class=""> </o:p></span></div><div style="margin: 0cm; text-align: justify; font-size: 10.5pt; font-family: DengXian;" class=""><span lang="EN-US" style="font-size: 12pt;" class="">What I have done. I set rho_0=rhoref (got from RhoRef.data & RhoRef.meta), but the result is not the same as the phiHyd (got from phiHyd.00000******.t0**.nc).<o:p class=""></o:p></span></div><div style="margin: 0cm; text-align: justify; font-size: 10.5pt; font-family: DengXian;" class=""><span lang="EN-US" style="font-size: 12pt;" class=""><o:p class=""> </o:p></span></div><div style="margin: 0cm; text-align: justify; font-size: 10.5pt; font-family: DengXian;" class=""><span lang="EN-US" style="font-size: 12pt;" class=""><o:p class=""> </o:p></span></div><div style="margin: 0cm; text-align: justify; font-size: 10.5pt; font-family: DengXian;" class=""><span lang="EN-US" style="font-size: 12pt;" class="">Leming<o:p class=""></o:p></span></div></div><span id="cid:0B3B1199-8052-48F2-8171-72461775CC4C"><屏幕截图 2022-09-29 095309.jpg></span><span id="cid:7389E875-37D2-4D6B-A3A3-7DA7558EBF1D"><屏幕截图 2022-10-07 090528.jpg></span><span style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;" class="">_______________________________________________</span><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class=""><span style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;" class="">MITgcm-support mailing list</span><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class=""><a href="mailto:MITgcm-support@mitgcm.org" style="color: blue; text-decoration: underline; font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;" class="">MITgcm-support@mitgcm.org</a><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class=""><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" style="color: blue; text-decoration: underline; font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a></div></blockquote></div><br class=""></div></body></html>