<html><head><meta http-equiv="content-type" content="text/html; charset=utf-8"></head><body style="overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">Hi Paola,<div><br></div><div>I never tried this, but I think it is important to know how the Coriolis terms are evaluted in the simulation. For this you need to check the relevant parameters in your STDOUT.0000 (do you use the mom-fluxform or the vector invariant form? parameters like: selectCoriScheme, useJamartWetPoints, useEnergyConservingCoriolis, selectVortScheme, SadournyCoriolis), and then it will be a bit tedious to find out, what the model actually does for the specific parameter set, but as far as I recall, this sets how the boundary conditions are computed, which in turn will give the “correct” (as the model evaluates them) form of the div(fu,fv) terms along the boundary. Does that make sense?</div><div><br></div><div>Martin<br><div><br><blockquote type="cite"><div>On 5. Sep 2023, at 08:03, Paola Cessi <pcessi@ucsd.edu> wrote:</div><br class="Apple-interchange-newline"><div>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">Dear colleagues,<div><br></div><div>I am trying to calculate the circulation balance over a semi-enclosed area, by taking the vorticity balance and integrating over said area. One of the terms in the vorticity balance is div (fu, fv), which, when integrated over the area, should give the integral over the boundary of (fu,fv) \cdot n, where n is the normal to the boundary. Over the closed portion of the boundary, i.e. on the solid walls, (fu,fv) \cdot n should be zero. However, because the vorticity points are not on the velocity boundary faces, these terms do not vanish. The residual left from this contribution dominates the balance, so it is not a small error.</div><div><br></div><div>I have tried different versions of the vorticity definitions, i.e. the volume flux form as in (equation 13) of Adcroft et al <a href="https://journals.ametsoc.org/view/journals/mwre/132/12/mwr2823.1.xml">https://journals.ametsoc.org/view/journals/mwre/132/12/mwr2823.1.xml</a>, or the simple derivatives and neither of them work. </div><div><br></div><div>Does anyone have any words of wisdom or references that could be helpful?</div><div><br></div><div>Many thanks,<br><div>
<div dir="auto" style="caret-color: rgb(0, 0, 0); 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; overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><div style="font-variant-caps: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;"><span style="font-style: normal;">Paola Cessi<br></span><br><span style="font-family: Helvetica; font-style: normal; font-weight: normal; font-size: 12px;">--------------------------------------------------------------------------------------------<br>Paola Cessi </span><span style="font-family: Helvetica; font-style: normal; font-weight: normal; font-size: 12px;">Tel: +1 858 534 0622 </span><span style="font-family: Helvetica; font-style: normal; font-weight: normal; font-size: 12px;"><br>Scripps Institution of Oceanography </span><span style="font-family: Helvetica; font-style: normal; font-weight: normal; font-size: 12px;"> Fax: +1 858 534 8045</span><span style="font-family: Helvetica; font-style: normal; font-weight: normal; font-size: 12px;"><br>9500 Gilman Drive #0213 </span><span style="font-family: Helvetica; font-style: normal; font-weight: normal; font-size: 12px;"><a href="mailto:pcessi@ucsd.edu">e-mail: pcessi@ucsd.edu</a> </span><span style="font-family: Helvetica; font-style: normal; font-weight: normal; font-size: 12px;"><br>La Jolla, CA 92093-0213 </span><span style="font-family: Helvetica; font-style: normal; font-weight: normal; font-size: 12px;">Web: <a href="http://pordlabs.ucsd.edu/pcessi">http://pordlabs.ucsd.edu/pcessi</a></span><span style="font-family: Helvetica; font-style: normal; font-weight: normal; font-size: 12px;"><br>USA <br><br></span></div></div></div></div></div></div>_______________________________________________<br>MITgcm-support mailing list<br>MITgcm-support@mitgcm.org<br>http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br></div></blockquote></div><br></div></body></html>