<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=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@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:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:"\@DengXian";
panose-1:2 1 6 0 3 1 1 1 1 1;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
span.EmailStyle18
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body lang="EN-US" link="blue" vlink="purple" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal">Hi Matt, <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thanks for the explanation. I track down the code you list to mom_calc_relvort3.F. I guess the terms in the vorticity equation can be coded for diagnostic purpose in a similar way. Or, the vorticity budget can be calculated post-run, since
the velocities, viscosities, momHDiv, and momVort3 can be diagnosed. <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Cheers,<o:p></o:p></p>
<p class="MsoNormal">Yilang<o:p></o:p></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Times New Roman",serif;color:black"><o:p> </o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal"><b><span style="font-size:12.0pt;color:black">From: </span></b><span style="font-size:12.0pt;color:black">MITgcm-support <mitgcm-support-bounces@mitgcm.org> on behalf of Matthew Mazloff <mmazloff@ucsd.edu><br>
<b>Reply-To: </b>"mitgcm-support@mitgcm.org" <mitgcm-support@mitgcm.org><br>
<b>Date: </b>Wednesday, January 25, 2023 at 14:27<br>
<b>To: </b>"mitgcm-support@mitgcm.org" <mitgcm-support@mitgcm.org><br>
<b>Cc: </b>Yilang Xu <yxu@whoi.edu><br>
<b>Subject: </b>Re: [MITgcm-support] Diagnose terms in vorticity equation<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Hello<o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Once you have closed the momentum budget you can close the vorticity budget by taking the curl.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">When we do this we use the stencil:<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> vort3(I,J)=<o:p></o:p></p>
</div>
<p class="MsoNormal"> & recip_rAz(I,J,bi,bj)*(<br>
& ( vFld(I,J)*dyC(I,J,bi,bj)<br>
& -vFld(I-1,J)*dyC(I-1,J,bi,bj) )<br>
& -( uFld(I,J)*dxC(I,J,bi,bj)<br>
& -uFld(I,J-1)*dxC(I,J-1,bi,bj) )<br>
& )<br>
<br>
As this makes pxy - pyx = 0 <o:p></o:p></p>
<div>
<p class="MsoNormal">and otherwise it doesn't always (e.g. dx is often a function of y).<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">-Matt<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">On Wed, Jan 25, 2023 at 6:36 AM Yilang Xu <<a href="mailto:yxu@whoi.edu">yxu@whoi.edu</a>> wrote:<o:p></o:p></p>
</div>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-right:0in">
<div>
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Dear MITgcm community,
<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">I am wondering if anyone has diagnosed the terms in vorticity equation in MITgcm. I see that there are diagnostics for momentum budget, and I am curious if there are similar things
for vorticity equation. <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Thank you very much.
<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Best,<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Yilang<o:p></o:p></p>
</div>
</div>
<p class="MsoNormal">_______________________________________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
<a href="https://urldefense.com/v3/__http:/mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!Hn9HembHHugP5pPzonRkldVlrpJux6OKounalValIwp9nJMXH2HsglZDt7FRQ9sfBvCc2LrSEldRTw$" target="_blank">https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!Hn9HembHHugP5pPzonRkldVlrpJux6OKounalValIwp9nJMXH2HsglZDt7FRQ9sfBvCc2LrSEldRTw$</a>
<o:p></o:p></p>
</div>
</blockquote>
</div>
</div>
</body>
</html>