<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=gb2312">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Helvetica;
panose-1:2 11 5 4 2 2 2 2 2 4;}
@font-face
{font-family:SimSun;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:DengXian;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:Lato;
panose-1:2 15 5 2 2 2 4 3 2 3;}
@font-face
{font-family:DengXian;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:SimSun;
panose-1:2 1 6 0 3 1 1 1 1 1;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
text-align:justify;
text-justify:inter-ideograph;
font-size:10.5pt;
font-family:DengXian;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
span.apple-converted-space
{mso-style-name:apple-converted-space;}
.MsoChpDefault
{mso-style-type:export-only;}
/* Page Definitions */
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 90.0pt 72.0pt 90.0pt;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body lang="ZH-CN" link="blue" vlink="#954F72" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US">Hi Martin,</span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">You are right (rhoNil is \rho_c, rhoConst is \rho_0), thank you very much!</span><span lang="EN-US" style="font-size:12.0pt"><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:12.0pt;font-family:SimSun">Cheers<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:12.0pt;font-family:SimSun"><o:p> </o:p></span></p>
<div style="mso-element:para-border-div;border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal" style="border:none;padding:0cm"><b>发件人<span lang="EN-US">: </span>
</b><span lang="EN-US"><a href="mailto:Martin.Losch@awi.de">Martin Losch</a><br>
</span><b>发送时间<span lang="EN-US">: </span></b><span lang="EN-US">2022</span>年<span lang="EN-US">10</span>月<span lang="EN-US">10</span>日<span lang="EN-US"> 16:51<br>
</span><b>收件人<span lang="EN-US">: </span></b><span lang="EN-US"><a href="mailto:mitgcm-support@mitgcm.org">MITgcm Support</a><br>
</span><b>主题<span lang="EN-US">: </span></b><span lang="EN-US">Re: [MITgcm-support] How to calculate the phiHyd (hydrostatic pressure potential anomaly) with model output ?</span></p>
</div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:12.0pt;font-family:SimSun"><o:p> </o:p></span></p>
<p class="MsoNormal" align="left" style="text-align:left"><span lang="EN-US">Hi,</span><span lang="EN-US" style="font-size:12.0pt"><o:p></o:p></span></p>
<div>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US">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).</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US">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.</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US">Hope this helps,</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US">Martin</span></p>
<div>
<p class="MsoNormal"><span lang="EN-US"><br>
<br>
</span></p>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<p class="MsoNormal"><span lang="EN-US">On 7. Oct 2022, at 08:39, Mike Fan <<a href="mailto:ifanliming@outlook.com">ifanliming@outlook.com</a>> wrote:</span></p>
</div>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<div>
<p class="MsoNormal"><span lang="EN-US">Hello all,</span><span lang="EN-US"><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> </span><span lang="EN-US"><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">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"><a href="https://mitgcm.readthedocs.io/en/latest/examples/baroclinic_gyre/baroclinic_gyre.html#output-files">4.2.
Baroclinic Ocean Gyre<span class="apple-converted-space"> </span><span lang="EN-US"><span lang="EN-US">― MITgcm checkpoint68k-6-g80404458a documentation</span></span></a>.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">In my model, setting a linaer EOS,<span class="apple-converted-space"> </span></span><span lang="EN-US" style="font-family:"Lato",sans-serif;color:black;background:#F3F6F6">implicitFreeSurface=’true’</span><span lang="EN-US">,
and rhoConst=999.8. Here is my question, there is a rho_0, how the model calculates it?<span class="apple-converted-space"> </span></span><span lang="EN-US"><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> </span><span lang="EN-US"><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">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).</span><span lang="EN-US"><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> </span><span lang="EN-US"><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> </span><span lang="EN-US"><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Leming</span><span lang="EN-US"><o:p></o:p></span></p>
<p class="MsoNormal" align="left" style="text-align:left"><span lang="EN-US"><</span>屏幕截图<span lang="EN-US"> 2022-09-29 095309.jpg><</span>屏幕截图<span lang="EN-US"> 2022-10-07 090528.jpg></span><span lang="EN-US" style="font-size:9.0pt;font-family:"Helvetica",sans-serif">_______________________________________________<br>
MITgcm-support mailing list<br>
</span><span lang="EN-US"><a href="mailto:MITgcm-support@mitgcm.org"><span style="font-size:9.0pt;font-family:"Helvetica",sans-serif">MITgcm-support@mitgcm.org</span></a></span><span lang="EN-US" style="font-size:9.0pt;font-family:"Helvetica",sans-serif"><br>
</span><span lang="EN-US"><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support"><span style="font-size:9.0pt;font-family:"Helvetica",sans-serif">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</span></a></span><span lang="EN-US" style="font-size:12.0pt;font-family:SimSun"><o:p></o:p></span></p>
</div>
</blockquote>
</div>
</div>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:12.0pt;font-family:SimSun"><o:p> </o:p></span></p>
</div>
</body>
</html>