<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=""><div class=""><br class=""></div>Hi Nadav Mantel<div class=""><br class=""></div><div class="">Sorry for the apparent bug in KL10 (tracked <a href="https://github.com/MITgcm/MITgcm/issues/647" class="">https://github.com/MITgcm/MITgcm/issues/647</a>) - there is no fundamental reason it should not work for non-linear equations of state, but apparently has a bug in how density is determined that matters for non-linear equation of state.  Should be straightforward to fix, though I’m </div><div class=""><br class=""></div><div class="">On the other hand, the mixing parameterizations can pretty easily produce large co-efficients if there are unstable regions.  If you are getting these when you don’t expect them, you should check your forcing.  Maybe do a short run with constant viscosity.  You can also limit the largest viscosity in most of them.  </div><div class=""><br class=""></div><div class="">Cheers,   Jody<br class=""><div class=""><br class=""><br class=""><br class="">
</div>
<div><br class=""><blockquote type="cite" class=""><div class="">On Aug 25, 2022, at 2:49 AM, Nadav Mantel <<a href="mailto:nadav.mantel@mail.huji.ac.il" class="">nadav.mantel@mail.huji.ac.il</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class="">Hey everyone,<div class=""><br class=""></div><div class="">This is a topic that comes up every now and then but I would still really appreciate someone's help.</div><div class="">I am running a non-hydrostatic regional model for the Gulf of Aqaba\Eilat, an elongated semi-enclosed basin [~200 km long, ~30 km wide ~1000 m deep] with an entry through a shallow sill in the south. The model runs with a horizontal resolution of 300 m and 32 vertical layers with the  differences in depth of 5,10, 12.5, 15, 17.5, ... m and a time step of 10s. I'm forcing tides through the southern boundary using the obcs package. The barotropic flow crossing the Straits of Tiran (our shallow sill) creates large amplitude (50 m in observations) internal waves along the gulf.</div><div class=""><br class=""></div><div class="">We would like to parameterize the shoaling and general breaking of the internal waves using one of the many packages MITgcm has to offer but we have had some troubles.</div><div class=""><br class=""></div><div class="">We first intended to use the KL10 package, yet we had uncharacteristic seasonal mixing and when checking the viscosity output we had weird results where in some places we got viscosity coefficients of 3 [m^2/s] whereas the rest of the gulf is much smaller. It could be possibly due to non-linear EOS, I saw in the mitgcm bug tracker that there are issues with anything but linear EOS use in the package.</div><div class=""><br class=""></div><div class="">We moved on to the richardson number packages, starting with pp81 as it is the simplest. Our problem was that we consistently got the maximum viscosity (a default of 1) in almost all of the second layer, creating once more uncharacteristic mixing with very large diffusivities (60 m^2/s).</div><div class=""><br class=""></div><div class="">We thought about using KPP as it was the package used in simulating the gulf up until now, but it wasn't designed for this purpose and we are skeptical it would produce good results. </div><div class=""><br class=""></div><div class="">Kind of at a loss at this point, and it would be really amazing to get someone's feedback via zoom call on their experiences with internal wave parameterization using any package for regional models.</div><div class=""><br class=""></div><div class="">Thanks,</div><div class="">Nadav Mantel</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></body></html>