<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="text-align: left; text-indent: 0px; background-color: rgb(255, 255, 255); margin: 0px; font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: black;" class="elementToProof">
Hi Jody,</div>
<div style="text-align: left; text-indent: 0px; background-color: rgb(255, 255, 255); margin: 0px; font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: black;" class="elementToProof">
<br>
</div>
<div style="text-align: left; text-indent: 0px; background-color: rgb(255, 255, 255); margin: 0px; font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: black;" class="elementToProof">
I have been reading your reference materials which is useful. Finally I was able to output the energy dignostics in latest version based on your code, but several questions have arisen:</div>
<div style="text-align: left; text-indent: 0px; background-color: rgb(255, 255, 255); margin: 0px; font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: black;" class="elementToProof">
<br>
</div>
<div style="text-align: justify; text-indent: 0px; background-color: rgb(255, 255, 255); margin: 0px; font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: black;" class="elementToProof">
1, For <b>baroclinic</b> <b>energy</b> equation <img style="width: 274px; height: 32px; max-width: 753px; min-width: auto; min-height: auto; margin: 0px;" size="3516" height="32" width="274" id="x_image_0" data-outlook-trace="F:2|T:2" src="cid:592a7939-f6bc-48b9-b311-cb17b6089067"> how
did you calculate the first three terms from dignostics outputs, while it seems the equation used in dignostics code is the
<b>total energy</b> equation? I can only find your github pages cause the link in your paper seems no longer valid. I'm confused that Fbc should be <u'p'>, but the dignostics code use the original velosity
<b>uVel</b> and <b>vVel</b> directly to get uPbc and so on.</div>
<div style="text-align: justify; text-indent: 0px; background-color: rgb(255, 255, 255); margin: 0px; font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: black;" class="elementToProof">
<br>
</div>
<div style="text-align: justify; text-indent: 0px; background-color: rgb(255, 255, 255); margin: 0px; font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: black;" class="elementToProof">
2, For Conversion term, the W is derived from <img style="width: 218px; height: 35px; max-width: 753px; min-width: auto; min-height: auto; margin: 0px;" size="6154" height="35" width="218" id="x_image_2" data-outlook-trace="F:2|T:2" src="cid:791c1654-63b9-4bb1-9f91-ec3c854ef76b">in
Kang's 2011 paper. So we don't need the W velosity from the model output? I'm not sure if they are different or my code was wrong.</div>
<div style="text-align: justify; text-indent: 0px; background-color: rgb(255, 255, 255); margin: 0px; font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: black;" class="elementToProof">
<br>
</div>
<div style="text-align: justify; text-indent: 0px; background-color: rgb(255, 255, 255); margin: 0px; font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: black;" class="elementToProof">
Sorry if it’s a lot to explain, but I could really use your help with this. </div>
<div style="text-align: left; text-indent: 0px; background-color: rgb(255, 255, 255); margin: 0px; font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: black;" class="elementToProof">
<br>
</div>
<div style="text-align: left; text-indent: 0px; background-color: rgb(255, 255, 255); margin: 0px; font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: black;" class="elementToProof">
Thanks and cheers!</div>
<div style="text-align: left; text-indent: 0px; background-color: rgb(255, 255, 255); margin: 0px; font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: black;" class="elementToProof">
Sean</div>
<div style="text-align: left; text-indent: 0px; background-color: rgb(255, 255, 255); margin: 0px; font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: black;" class="elementToProof">
University of East Anglia</div>
<div style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<br>
</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>发件人:</b> MITgcm-support <mitgcm-support-bounces@mitgcm.org> 代表 Jody Klymak <jklymak@uvic.ca><br>
<b>发送时间:</b> 2025年7月19日 1:18<br>
<b>收件人:</b> mitgcm-support@mitgcm.org <mitgcm-support@mitgcm.org><br>
<b>主题:</b> Re: [MITgcm-support] Question about internal tide energy flux balance</font>
<div> </div>
</div>
<style>
<!--
@font-face
{font-family:"Cambria Math"}
@font-face
{font-family:Calibri}
@font-face
{font-family:Aptos}
@font-face
{font-family:wf_segoe-ui_normal}
p.x_MsoNormal, li.x_MsoNormal, div.x_MsoNormal
{margin:0cm;
font-size:12.0pt;
font-family:"Aptos",sans-serif}
a:link, span.x_MsoHyperlink
{color:blue;
text-decoration:underline}
span.x_EmailStyle19
{font-family:"Aptos",sans-serif;
color:windowtext}
.x_MsoChpDefault
{font-size:10.0pt}
@page WordSection1
{margin:72.0pt 72.0pt 72.0pt 72.0pt}
div.x_WordSection1
{}
-->
</style>
<div lang="EN-CA" link="blue" vlink="purple" style="word-wrap:break-word">
<p style="border:4px; border-style:solid; border-color:#FF0000; padding:0.5em; text-align:center">
<b>Warning:</b> This email is from outside the UEA system. Do not click on links or attachments unless you expect them from the sender and know the content is safe.
</p>
<div></div>
<div>
<div class="x_WordSection1">
<p class="x_MsoNormal"><span style="font-size:11.0pt">Hi Sean</span></p>
<p class="x_MsoNormal"><span style="font-size:11.0pt"> </span></p>
<p class="x_MsoNormal"><span style="font-size:11.0pt">This is probably outdated, but the writeup and the diagnostic may help you with your problem:</span></p>
<p class="x_MsoNormal"><span style="font-size:11.0pt"> </span></p>
<p class="x_MsoNormal"><span style="font-size:11.0pt"><a href="https://ocean-physics.seos.uvic.ca/~jklymak/MITgcm/" originalsrc="https://ocean-physics.seos.uvic.ca/~jklymak/MITgcm/">https://ocean-physics.seos.uvic.ca/~jklymak/MITgcm/</a></span></p>
<p class="x_MsoNormal"><span style="font-size:11.0pt"> </span></p>
<p class="x_MsoNormal"><span style="font-size:11.0pt">in particular <a href="https://ocean-physics.seos.uvic.ca/~jklymak/MITgcm/energy_diagnostics_fill.F" originalsrc="https://ocean-physics.seos.uvic.ca/~jklymak/MITgcm/energy_diagnostics_fill.F">
https://ocean-physics.seos.uvic.ca/~jklymak/MITgcm/energy_diagnostics_fill.F</a></span></p>
<p class="x_MsoNormal"><span style="font-size:11.0pt">and </span></p>
<p class="x_MsoNormal"><span style="font-size:11.0pt"><a href="https://ocean-physics.seos.uvic.ca/~jklymak/MITgcm/EnergyDiagnostic.pdf" originalsrc="https://ocean-physics.seos.uvic.ca/~jklymak/MITgcm/EnergyDiagnostic.pdf">https://ocean-physics.seos.uvic.ca/~jklymak/MITgcm/EnergyDiagnostic.pdf</a></span></p>
<p class="x_MsoNormal"><span style="font-size:11.0pt"> </span></p>
<p class="x_MsoNormal"><span style="font-size:11.0pt">This was used extensively in</span>
<span style="font-size:11.0pt"><a href="https://doi.org/10.1175/2010JPO4396.1" originalsrc="https://doi.org/10.1175/2010JPO4396.1" target="_blank">10.1175/2010JPO4396.1</a> and 10.1175/JPO-D-16-0061.1</span></p>
<p class="x_MsoNormal"><span style="font-size:11.0pt"> </span></p>
<p class="x_MsoNormal"><span style="font-size:11.0pt">Hopefully that helps, if not, I can probably provide more guidance.
</span></p>
<p class="x_MsoNormal"><span style="font-size:11.0pt"> </span></p>
<p class="x_MsoNormal"><span style="font-size:11.0pt">Cheers Jody</span></p>
<p class="x_MsoNormal"><span style="font-size:11.0pt"> </span></p>
<div style="border:none; border-top:solid #B5C4DF 1.0pt; padding:3.0pt 0cm 0cm 0cm">
<p class="x_MsoNormal"><b><span style="font-family:"Calibri",sans-serif; color:black">From:
</span></b><span style="font-family:"Calibri",sans-serif; color:black">MITgcm-support <mitgcm-support-bounces@mitgcm.org> on behalf of "Shiqiang Hu (ENV - Postgraduate Researcher)" <Shiqiang.Hu@uea.ac.uk><br>
<b>Reply-To: </b>"mitgcm-support@mitgcm.org" <mitgcm-support@mitgcm.org><br>
<b>Date: </b>Friday, July 18, 2025 at 9:58 AM<br>
<b>To: </b>"mitgcm-support@mitgcm.org" <mitgcm-support@mitgcm.org><br>
<b>Subject: </b>[MITgcm-support] Question about internal tide energy flux balance</span></p>
</div>
<div>
<p class="x_MsoNormal"> </p>
</div>
<table class="x_MsoNormalTable" border="0" cellspacing="0" cellpadding="0" align="left" width="100%" style="width:100.0%; background:revert!important; border:revert!important; bottom:revert!important; color:revert!important; direction:revert!important; display:revert!important; font-size:revert!important; height:revert!important; letter-spacing:revert!important; line-height:revert!important; margin:revert!important; opacity:revert!important; order:revert!important; outline:revert!important; overflow:revert!important; padding:revert!important; position:revert!important; tab-size:revert!important; table-layout:revert!important; text-align:revert!important; text-indent:revert!important; text-orientation:revert!important; text-overflow:revert!important; text-transform:revert!important; top:revert!important; vertical-align:revert!important; visibility:revert!important; white-space:revert!important; width:revert!important; word-break:revert!important; word-spacing:revert!important; writing-mode:revert!important; zoom:revert!important; display:table!important; border-collapse:seperate!important; border-spacing:0px 0px!important; float:none">
<tbody>
<tr style="background:revert!important; border:revert!important; bottom:revert!important; color:revert!important; direction:revert!important; display:revert!important; font-size:revert!important; height:revert!important; letter-spacing:revert!important; line-height:revert!important; margin:revert!important; opacity:revert!important; order:revert!important; outline:revert!important; overflow:revert!important; padding:revert!important; position:revert!important; tab-size:revert!important; table-layout:revert!important; text-align:revert!important; text-indent:revert!important; text-orientation:revert!important; text-overflow:revert!important; text-transform:revert!important; top:revert!important; vertical-align:revert!important; visibility:revert!important; white-space:revert!important; width:revert!important; word-break:revert!important; word-spacing:revert!important; writing-mode:revert!important; zoom:revert!important">
<td width="0" style="width:.3pt; background:#A6A6A6; padding:5.25pt 1.5pt 5.25pt 1.5pt; background:revert!important; border:revert!important; bottom:revert!important; color:revert!important; direction:revert!important; display:revert!important; font-size:revert!important; height:revert!important; letter-spacing:revert!important; line-height:revert!important; margin:revert!important; opacity:revert!important; order:revert!important; outline:revert!important; overflow:revert!important; padding:revert!important; position:revert!important; tab-size:revert!important; table-layout:revert!important; text-align:revert!important; text-indent:revert!important; text-orientation:revert!important; text-overflow:revert!important; text-transform:revert!important; top:revert!important; vertical-align:revert!important; visibility:revert!important; white-space:revert!important; width:revert!important; word-break:revert!important; word-spacing:revert!important; writing-mode:revert!important; zoom:revert!important">
</td>
<td width="100%" style="width:100.0%; background:#EAEAEA; padding:5.25pt 3.75pt 5.25pt 11.25pt; background:revert!important; border:revert!important; bottom:revert!important; color:revert!important; direction:revert!important; display:revert!important; font-size:revert!important; height:revert!important; letter-spacing:revert!important; line-height:revert!important; margin:revert!important; opacity:revert!important; order:revert!important; outline:revert!important; overflow:revert!important; padding:revert!important; position:revert!important; tab-size:revert!important; table-layout:revert!important; text-align:revert!important; text-indent:revert!important; text-orientation:revert!important; text-overflow:revert!important; text-transform:revert!important; top:revert!important; vertical-align:revert!important; visibility:revert!important; white-space:revert!important; width:revert!important; word-break:revert!important; word-spacing:revert!important; writing-mode:revert!important; zoom:revert!important">
<div>
<p class="x_MsoNormal" style=""><span style="font-size:9.0pt; font-family:"wf_segoe-ui_normal",serif; color:#212121">You don't often get email from shiqiang.hu@uea.ac.uk.
<a href="https://aka.ms/LearnAboutSenderIdentification">Learn why this is important</a>
</span></p>
</div>
</td>
<td width="75" style="width:56.25pt; background:#EAEAEA; padding:5.25pt 3.75pt 5.25pt 3.75pt; background:revert!important; border:revert!important; bottom:revert!important; color:revert!important; direction:revert!important; display:revert!important; font-size:revert!important; height:revert!important; letter-spacing:revert!important; line-height:revert!important; margin:revert!important; opacity:revert!important; order:revert!important; outline:revert!important; overflow:revert!important; padding:revert!important; position:revert!important; tab-size:revert!important; table-layout:revert!important; text-align:revert!important; text-indent:revert!important; text-orientation:revert!important; text-overflow:revert!important; text-transform:revert!important; top:revert!important; vertical-align:revert!important; visibility:revert!important; white-space:revert!important; width:revert!important; word-break:revert!important; word-spacing:revert!important; writing-mode:revert!important; zoom:revert!important">
</td>
</tr>
</tbody>
</table>
<div>
<div>
<p class="x_MsoNormal"><span style="color:black">Hi MITgcm commuity,</span></p>
</div>
<div>
<p class="x_MsoNormal"><span style="color:black"> </span></p>
</div>
<div>
<p class="x_MsoNormal"><span style="color:black">Hope you all well.</span></p>
</div>
<div>
<p class="x_MsoNormal"><span style="color:black"> </span></p>
</div>
<div>
<p class="x_MsoNormal"><span style="color:black">I built an idealized 3D internal tide model based on
<a href="https://github.com/MITgcm/MITgcm/blob/master/verification/internal_wave" originalsrc="https://github.com/MITgcm/MITgcm/blob/master/verification/internal_wave" title="https://github.com/MITgcm/MITgcm/blob/master/verification/internal_wave">
internal_wave experiment</a> to calculate the energy flux F=<u'p'>(<a href="https://doi.org/10.1175/JTECH1784.1" originalsrc="https://doi.org/10.1175/JTECH1784.1" title="https://doi.org/10.1175/JTECH1784.1">Nash et al., 2005</a>). Periodic 1-mode internal tides
propagate to the right and interact with a simple topography. The model works but I got stuck at the energy flux balance part. Sometimes when I chose a small area within the topography to calculate the net energy flux across its boundary, the result is positive,
which means there are rather dissipation but genaration of enerny. So I wanna ask if someone has experience with this. </span></p>
</div>
<div>
<p class="x_MsoNormal"><span style="color:black">Basically, I use <i>RHOAnoma</i>+<i>rhoConst</i> to get original density
<i>rho</i>, cause the paper gets density anomaly from <i><span style="background:white">rho(z,t)-rho_mean(z),
</span></i><span style="background:white">where<i> rho_mean(z) </i>is the time mean vertical density profile (I guess I could also just use the
<i>Rhoref(z) </i>as the reference profile?). The question is the C-grids, should I calculate these values on T-points or U-V points so the energy could be more balenced? Or maybe some other factors I didn't realize? </span></span></p>
</div>
<div>
<p class="x_MsoNormal"><span style="color:black; background:white">By the way, the non-hydrostatic pressure is quite small, so I ignored it.</span><span style="color:black"></span></p>
</div>
<div>
<p class="x_MsoNormal"><span style="color:black"> </span></p>
</div>
<div>
<p class="x_MsoNormal"><span style="color:black; background:white">Thanks for any suggestions.</span><span style="color:black"></span></p>
</div>
<div>
<p class="x_MsoNormal"><span style="color:black; background:white">Kind regards,</span><span style="color:black"></span></p>
</div>
<div>
<p class="x_MsoNormal"><span style="color:black"> </span></p>
</div>
<div>
<p class="x_MsoNormal"><span style="color:black; background:white">Sean</span><span style="color:black"></span></p>
</div>
<div>
<p class="x_MsoNormal"><span style="color:black">University of East Anglia</span></p>
</div>
</div>
</div>
</div>
</div>
</body>
</html>