<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
Hi Navad, here is a recently-published paper that may be relevant to some of the questions that you ask: <a href="https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2022GL099614" class="">https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2022GL099614</a>
<div class=""><br class="">
</div>
<div class="">But note that the above study pertains to the representation of internal waves in an open ocean region, as opposed to a narrow, elongated gulf.</div>
<div class=""><br class="">
</div>
<div class="">
<div class="">Dimitris Menemenlis </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="">
https://urldefense.us/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!PvBDto6Hs4WbVuu7!Z9iqcR5gWjF-J_Axv-Z5IwjHE3NIjlnukR1MoYYCJsMmMBhhzukH4-GMMAi7fj3nI9TXJBvuY6M$
<br class="">
</div>
</blockquote>
</div>
<br class="">
</div>
</body>
</html>