<div><font>Hi Bertrand,</font></div><div><font><br></font></div><div><font>I tried to close the momentum/energy budget using MITgcm one year ago and</font></div><div><font>suffered a lot as you. So I write two notebooks and collect all the information</font></div><div><font>I can find. You may find them useful:</font></div><div><font><a href="https://github.com/miniufo/notebooks/blob/master/Budget%20in%20MITgcm%20Part%20III.ipynb">https://github.com/miniufo/notebooks/blob/master/Budget%20in%20MITgcm%20Part%20III.ipynb</a></font></div><div><a href="https://github.com/miniufo/notebooks/blob/master/Budget%20in%20MITgcm%20Part%20IV.ipynb">https://github.com/miniufo/notebooks/blob/master/Budget%20in%20MITgcm%20Part%20IV.ipynb</a></div><div><br></div><div><font>However, I did not remember any clue on close w-momentum budget in</font></div><div><font>nonhydrostatic case. Maybe you are the first one trying to do it.</font></div><div><font><br></font></div><div><font>On energy budget, you are right that you have to take the dot product of</font></div><div><font>velocity vector and momentum equations. I remember that there is a major</font></div><div><font>change of the diagnostic package that some diagnostics have been re-grouped</font></div><div>to simplified the budget (you can find this in the reference part of the notebook).</div><div>So you have to make sure the version of the model that contains the change or</div><div>not.</div><div><sign signid="0"><div><font><br></font><font><br></font><font><br></font><font><br></font><div style="color:#909090;font-family:Arial Narrow;font-size:12px">------------------</div></div><div style="font-size:14px;font-family:Verdana;color:#000;"><div> <font size="2">Best regards <br> <br>Yu-Kun Qian (钱钰坤) <br>Center for Monsoon and Environment Research</font> </div> <font size="2">Department of Atmospheric Sciences<br>School of Environmental Science and Engineering<font size="2"> <br> </font>Sun Yat-sen University <br>No. 135 Xingang West Road, Haizhu District <br>Guangzhou, 510275, P.R. China <br>Tel; 020-84115227 <br>Email: <a href="mailto:qianyk@mail2.sysu.edu.cn">qianyk@mail3.sysu.edu.cn</a> </font>
<div> <font size="2"> </font> </div></div></sign></div><div> </div><div><tincludetail><div> </div><div> </div><div style="font:Verdana normal 14px;color:#000;"><div style="FONT-SIZE: 12px;FONT-FAMILY: Arial Narrow;padding:2px 0 2px 0;">------------------ 原始邮件 ------------------</div><div style="FONT-SIZE: 12px;background:#efefef;padding:8px;"><div id="menu_sender"><b>发件人:</b> "Bertrand Delorme"<bdelorme@stanford.edu>;</div><div><b>发送时间:</b> 2020年11月24日(星期二) 晚上9:31</div><div><b>收件人:</b> "mitgcm-support@mitgcm.org"<mitgcm-support@mitgcm.org>; <wbr></div><div></div><div><b>主题:</b> Re: [MITgcm-support] Non-Hydrostatic Energy Budget</div></div><div> </div><div style="position:relative;"><div id="tmpcontent_res"></div>
<style></style>
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt">Hi everyone,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">Follow-up to my previous post below. It seems that the diagnostics Wm_Advec, Wm_Diss and AB_gW are available with the code in NH mode. However, to have a complete budget, I would need the total
vertical pressure gradient (something like a </span><span lang="EN-US" style="font-size:11.0pt">Wm_dPhiZ diagnostic). It would also be nice to have the tendency of the vertical velocity (something like TOTWTEND). I can’t find any of these two diagnostics in
the code of the model. Is there a reason why those haven’t been included with the other diagnostics of the vertical momentum equation? Is there any ways to access them outside of an offline calculation?<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">Thank you,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">Bertrand<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><o:p> </o:p></span></p>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span style="color:black">De : </span></b><span style="color:black">Bertrand DELORME <bdelorme@stanford.edu><br>
<b>Date : </b>vendredi 13 novembre 2020 à 16:24<br>
<b>À : </b>"mitgcm-support@mitgcm.org" <mitgcm-support@mitgcm.org><br>
<b>Objet : </b>Non-Hydrostatic Momentum Budget</span><span style="color:black;mso-fareast-language:FR"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
</div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">Hi everyone,</span><o:p></o:p></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">I am trying to compute the fully non-hydrostatic energy budget of my internal wave simulation using the Diagnostics package and I have a few questions:</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"> </span><o:p></o:p></p>
<p class="MsoListParagraph" style="margin-left:54.0pt;text-indent:-36.0pt;mso-list:l0 level1 lfo2">
<!--[if !supportLists]--><span style="mso-list:Ignore">i)<span style=" ; ">
</span></span><span lang="EN-US" style="font-size:11.0pt">I read that the diagnostics Um_Advec and Vm_Advec already contain the Coriolis terms. Does that mean that in NH mode, Um_Advec contains both the -f*v and +fh*w, where f is the vertical component
of the Coriolis parameter and fh the horizontal component of the Coriolis parameter? When calculating the energy budget, one needs to take the dot product of the momentum equations with the velocity vector, effectively canceling out the Coriolis terms in the
energy equation. However, using only the horizontal momentum equations diagnostics (i.e., Um_Advec and Vm_Advec), we don’t get the -fh*u term coming from the vertical momentum equation that permits to cancel out the +fh*w term from the zonal momentum equation.
Is that really what is happening when we perform U . Um_Advec + V . Vm_Advec (i.e., we get a u*fh*w term that is not balanced)? If yes, is there a way to get a Wm_Advec term from the diagnostics package or do I have to calculate it offline when performing
my energy budget?</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"> </span><o:p></o:p></p>
<p class="MsoListParagraph" style="margin-left:54.0pt;text-indent:-36.0pt;mso-list:l0 level1 lfo2">
<!--[if !supportLists]--><span style="mso-list:Ignore">ii)<span style=" ; ">
</span></span><span lang="EN-US" style="font-size:11.0pt"> In a previous post from Jean-Michel, I read that the surface pressure horizontal gradient needs to be added to the Um_dPHdx and Vm_dPHdy diagnostics. However, when I run my simulation, I don’t
get the Um_dPHdx and Vm_dPHdy diagnostics but I get the Um_dPhiX and Vm_dPhiY diagnostics instead, which seem to already incorporate the surface pressure gradient. Am I right on that? Also, do these new diagnostics incorporate the NH pressure field?</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">Thank you,</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">Bertrand Delorme</span><o:p></o:p></p>
</div>
</div></div></tincludetail></div><!--<![endif]-->