<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body dir="auto">
Hi Josephine, would you by any chance have sea ice entering or exiting your domain?  Sea ice lateral fluxes would not be included in TOTTEND.
<div><br id="lineBreakAtBeginningOfSignature">
<div dir="ltr">Dimitris Menemenlis
<div>818-625-6498</div>
</div>
<div dir="ltr"><br>
<blockquote type="cite">On Sep 17, 2024, at 9:15 AM, Josephine Anselin - BAS <joslin53@bas.ac.uk> wrote:<br>
<br>
</blockquote>
</div>
<blockquote type="cite">
<div dir="ltr">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
 <style>@font-face { font-family: "Cambria Math"; }
@font-face { font-family: Calibri; }
p.MsoNormal, li.MsoNormal, div.MsoNormal { margin: 0cm; font-size: 11pt; font-family: Calibri, sans-serif; }
a:link, span.MsoHyperlink { color: rgb(5, 99, 193); text-decoration: underline; }
span.EmailStyle17 { font-family: Calibri, sans-serif; color: windowtext; }
.MsoChpDefault { font-family: Calibri, sans-serif; }
@page WordSection1 { size: 612pt 792pt; margin: 72pt; }
div.WordSection1 { page: WordSection1; }</style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
<div class="WordSection1">
<p class="MsoNormal">Hi everyone,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I am evaluating the heat budget in a 3D domain with nx=640 grid points in the x-direction; ny=1200 grid points in the y-direction and nz=150 grid points in the z direction. The domain includes open ocean, ice shelf and grounded ice cells.
 Temperature and salinity are restored at the open lateral boundaries. <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">The following parameter settings are used:<o:p></o:p></p>
<p class="MsoNormal">implicitFreeSurface=.TRUE.<o:p></o:p></p>
<p class="MsoNormal">nonlinFreeSurf=4.<o:p></o:p></p>
<p class="MsoNormal">useRealFreshWaterFlux=.FALSE.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">To begin with, I am trying to close the budget for the cells away from the ice shelf-ocean interface and away from the ocean surface. I was assuming that for these cells, the heat budget should be:<o:p></o:p></p>
<p class="MsoNormal">TOTTEND/86400 = Adv_tend + Dif_tend,<o:p></o:p></p>
<p class="MsoNormal">with Adv_tend and Dif_tend calculated as per page 3 in this document:
<a href="https://urldefense.us/v3/__https://mitgcm.org/download/daily_snapshot/MITgcm/doc/Heat_Salt_Budget_MITgcm.pdf__;!!PvBDto6Hs4WbVuu7!KF_zBQLnOlt1MolinKrV4tXovMauCTPoVUwB5SsQNc40LpaVGixzCmrmMV2UUGZHAgtYAwjB0Eiul4alXunzGNgbNnaKIMY$">
https://mitgcm.org/download/daily_snapshot/MITgcm/doc/Heat_Salt_Budget_MITgcm.pdf</a>
<o:p></o:p></p>
<p class="MsoNormal">(KPPg_TH and oceQsw are zero in the model configuration).<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">However, the heat budget does not close using this calculation.
<o:p></o:p></p>
<p class="MsoNormal">I would like to investigate the potential impact of the lateral heat flux due to temperature restoring to check if this might be causing the large residual.
<o:p></o:p></p>
<p class="MsoNormal">Is there a diagnostics analogous to TRELAX but for lateral heat flux instead of surface heat flux?<o:p></o:p></p>
<p class="MsoNormal">If not, does anyone have any suggestions on how this could be calculated using built-in diagnostics?<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thank you<o:p></o:p></p>
<p class="MsoNormal">Josephine</p>
</div>
</div>
</blockquote>
</div>
</body>
</html>