<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">I agree with Dimitris, and usually use 7 or 77, but note that then you will have spurious energy dissipation, as in \nu nabla^2 u < dE/dt, where as the second order scheme does a better job of preserving energy (you can always tune that away by increasing the lateral viscosity, but…)<div class=""><br class=""></div><div class="">After trying to get this all to work consistently, I think the best approach I generally take is to do an external budget (either energy or tracer) and assume the the residual relative to the dissipation terms in numerical dissipation. I believe that it is generally true that the advection schemes tend to develop numerical instability where you would expect physical instability anyway, so if you track it with an external budget you should be reflecting the model’s reality. </div><div class=""><br class=""></div><div class="">This is all assuming you are not trying to do a DNS. If you are doing a DNS, you should increase the viscosity and diffusivities to the point where this stop being a problem. But from your description it doesn’t sound like you are trying to do a DNS.</div><div class=""><br class=""></div><div class="">Cheers, Jody<br class=""><div><br class=""><blockquote type="cite" class=""><div class="">On 2 Feb 2021, at 10:49, Dimitris Menemenlis <<a href="mailto:menemenlis@jpl.nasa.gov" class="">menemenlis@jpl.nasa.gov</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><meta http-equiv="Content-Type" content="text/html; charset=utf-8" class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Dear Xiaozhou, did anybody respond to your question?<div class=""><br class=""><div class="">It is possible to reduce (but not completely eliminate) spurious mixing by choosing a different advection scheme.<div class="">I am partial to:</div><div class=""><br class=""></div><div class=""><div class="">C ENUM_OS7MP :: 7th Order One Step method with Monotonicity Preserving Limiter</div><div class=""> INTEGER ENUM_OS7MP</div><div class=""> PARAMETER(ENUM_OS7MP=7)</div><div class=""><br class=""></div><div class="">which you set in the “data” runtime parameter file using: </div><div class=""><div class=""> tempAdvScheme=7,</div><div class=""> saltAdvScheme=7,</div></div><div class=""><br class=""></div><div class="">Prather advection:</div><div class=""><br class=""></div><div class=""><div class="">C ENUM_SOM_PRATHER :: 2nd Order-Moment Advection Scheme, Prather, 1986</div><div class=""> INTEGER ENUM_SOM_PRATHER</div><div class=""> PARAMETER(ENUM_SOM_PRATHER=80)</div><div class=""><br class=""></div><div class="">C ENUM_SOM_LIMITER :: 2nd Order-Moment Advection Scheme, Prather Limiter</div><div class=""> INTEGER ENUM_SOM_LIMITER</div><div class=""> PARAMETER(ENUM_SOM_LIMITER=81)</div></div><div class=""><br class=""></div><div class="">is supposed to be even less diffusive but more expensive to compute.</div><div class=""><br class=""></div><div class="">There is a nice discussion on advection schemes in the MITgcm user manual:</div><div class=""><a href="http://mitgcm.org/sealion/online_documents/node71.html" class="">http://mitgcm.org/sealion/online_documents/node71.html</a></div><div class=""><a href="http://mitgcm.org/sealion/online_documents/node76.html" class="">http://mitgcm.org/sealion/online_documents/node76.html</a></div><div class=""><a href="http://mitgcm.org/sealion/online_documents/node81.html" class="">http://mitgcm.org/sealion/online_documents/node81.html</a></div><div class=""><br class=""></div><div class="">Cheers, Dimitris</div><div class=""><br class=""></div><div class=""><br class=""><blockquote type="cite" class=""><div class="">On Jan 29, 2021, at 9:58 AM, Ruan Xiaozhou <<a href="mailto:saberruan@hotmail.com" class="">saberruan@hotmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class="">
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
<div dir="auto" style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Dear MITgcm users,
<div class=""><br class="">
</div>
<div class="">Hope this message finds you well. I’ve been running simple 2D simulations with internal tides driven by oscillatory mean flows over rough bathymetry. When diagnosing the volume-integrated tracer variance budget, although the l.h.s. (time-tendency)
and r.h.s (advection + diffusion) terms balance exactly, I found some extra destruction of variance coming from the volume integral of the advection term using nonlinear advection schemes (scheme 33) which is as large as the diffusion term. Theoretically this
volume integral of variance advection should vanish... Switching to a centered 2nd order scheme (scheme 2) kills this contribution but the total variance destruction rate becomes about twice as large which, according to the output, can be explained by the
noisy tracer field and thus a larger gradient term.</div>
<div class=""><br class="">
</div>
<div class="">I was wondering if anyone had experience beating down this *spurious* contribution from the advection term? I noticed this contribution not only in the variance budget but also other 2nd order budgets involving a volume integral.</div>
<div class=""><br class="">
</div>
<div class="">Cheers,</div>
<div class="">Xiaozhou</div>
<div class="">
<div class="">
<div class="" style="text-align: -webkit-auto; line-height: normal; orphans: 2; widows: 2; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
--<br class="">
Xiaozhou Ruan</div>
<div class="" style="text-align: -webkit-auto; line-height: normal; orphans: 2; widows: 2; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
Postdoctoral researcher</div>
<div class="" style="text-align: -webkit-auto; line-height: normal; orphans: 2; widows: 2; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
Department of Earth, Atmospheric and Planetary Sciences</div>
<div class="" style="text-align: -webkit-auto; line-height: normal; orphans: 2; widows: 2; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
Massachusetts Institute of Technology<br class="">
54-1622<br class="">
Cambridge, MA 02139</div>
</div>
<div class=""><br class="">
</div>
<div class="">
<div class="">
<div class="" style="text-align: -webkit-auto; line-height: normal; orphans: 2; widows: 2; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
email: <a href="mailto:xiaozhour@caltech.edu" class="">xruan@mit.edu</a></div>
</div>
<div class="" style="text-align: -webkit-auto; line-height: normal; orphans: 2; widows: 2; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
web: <a href="https://urldefense.us/v3/__http://www.its.caltech.edu/*xruan__;fg!!PvBDto6Hs4WbVuu7!YRuAOFjbX-Ud-ObyNQx3QCA_mHxkdrR9q1gbhLpGtf0UsoeYre8nvFKMppLmQdKyrxkxhE2anEU$" style="orphans: auto; widows: auto;" class="">http://www.mit.edu/~xruan</a></div>
</div>
<div class=""><br class="">
</div>
<br class="">
</div>
</div>
</div>
_______________________________________________<br class="">MITgcm-support mailing list<br class=""><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class=""><a href="https://urldefense.us/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!PvBDto6Hs4WbVuu7!YRuAOFjbX-Ud-ObyNQx3QCA_mHxkdrR9q1gbhLpGtf0UsoeYre8nvFKMppLmQdKyrxkxXSPnbl0$" class="">https://urldefense.us/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!PvBDto6Hs4WbVuu7!YRuAOFjbX-Ud-ObyNQx3QCA_mHxkdrR9q1gbhLpGtf0UsoeYre8nvFKMppLmQdKyrxkxXSPnbl0$</a> <br class=""></div></blockquote></div><br class=""></div></div></div></div>_______________________________________________<br class="">MITgcm-support mailing list<br class=""><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br class=""></div></blockquote></div><br class=""><div class="">
<div style="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; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px; color: rgb(0, 0, 0); font-family: "Lucida Sans Typewriter"; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px;"><div style="word-wrap: break-word;" class=""><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px; color: rgb(0, 0, 0); font-family: "Lucida Sans Typewriter"; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px;"><div class="">--</div><div class="">Jody Klymak </div><div class=""><a href="http://ocean-physics.seos.uvic.ca/~jklymak/" class="">http://ocean-physics.seos.uvic.ca/~jklymak/</a></div><div class=""><br class="khtml-block-placeholder"></div><div class=""><br class="khtml-block-placeholder"></div><br class="Apple-interchange-newline"></span></div></span></div><br class="Apple-interchange-newline"><br class="Apple-interchange-newline">
</div>
<br class=""></div></body></html>