<div dir='auto'>Hi,<div dir="auto"><br></div><div dir="auto"><span style="color: rgb(209, 213, 219); font-family: söhne, ui-sans-serif, "system-ui", -apple-system, "segoe ui", roboto, ubuntu, cantarell, "noto sans", sans-serif, "helvetica neue", arial, "apple color emoji", "segoe ui emoji", "segoe ui symbol", "noto color emoji"; font-size: 16px; white-space-collapse: preserve; background-color: rgb(68, 70, 84);">I wanted to inform you that by implementing your suggestions, I was able to successfully resolve the problem. Thank you so much for your valuable input. </span><br></div><div dir="auto"><span style="color: rgb(209, 213, 219); font-family: söhne, ui-sans-serif, "system-ui", -apple-system, "segoe ui", roboto, ubuntu, cantarell, "noto sans", sans-serif, "helvetica neue", arial, "apple color emoji", "segoe ui emoji", "segoe ui symbol", "noto color emoji"; font-size: 16px; white-space-collapse: preserve; background-color: rgb(68, 70, 84);"><br></span></div><div dir="auto"><span style="color: rgb(209, 213, 219); font-family: söhne, ui-sans-serif, "system-ui", -apple-system, "segoe ui", roboto, ubuntu, cantarell, "noto sans", sans-serif, "helvetica neue", arial, "apple color emoji", "segoe ui emoji", "segoe ui symbol", "noto color emoji"; font-size: 16px; white-space-collapse: preserve; background-color: rgb(68, 70, 84);">Best regards,</span></div><div dir="auto"><span style="color: rgb(209, 213, 219); font-family: söhne, ui-sans-serif, "system-ui", -apple-system, "segoe ui", roboto, ubuntu, cantarell, "noto sans", sans-serif, "helvetica neue", arial, "apple color emoji", "segoe ui emoji", "segoe ui symbol", "noto color emoji"; font-size: 16px; white-space-collapse: preserve; background-color: rgb(68, 70, 84);"><br></span></div><div dir="auto"><span style="color: rgb(209, 213, 219); font-family: söhne, ui-sans-serif, "system-ui", -apple-system, "segoe ui", roboto, ubuntu, cantarell, "noto sans", sans-serif, "helvetica neue", arial, "apple color emoji", "segoe ui emoji", "segoe ui symbol", "noto color emoji"; font-size: 16px; white-space-collapse: preserve; background-color: rgb(68, 70, 84);">Pragnya Makar </span></div></div><div class="gmail_extra"><br><div class="gmail_quote">On 15 Jun 2023 21:30, mitgcm-support-request@mitgcm.org wrote:<br type="attribution" /><blockquote class="quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><p dir="ltr">Send MITgcm-support mailing list submissions to <br>
mitgcm-support@mitgcm.org <br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit <br>
http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support <br>
or, via email, send a message with subject or body 'help' to <br>
mitgcm-support-request@mitgcm.org <br>
<br>
You can reach the person managing the list at <br>
mitgcm-support-owner@mitgcm.org <br>
<br>
When replying, please edit your Subject line so it is more specific <br>
than "Re: Contents of MITgcm-support digest..." <br>
<br>
<br>
Today's Topics: <br>
<br>
   1. Re: [EXTERNAL] Re: Query regarding horizontal strips in <br>
      currents (Menemenlis, Dimitris (US 329B)) (Martin Losch) <br>
<br>
<br>
---------------------------------------------------------------------- <br>
<br>
Message: 1 <br>
Date: Thu, 15 Jun 2023 12:12:13 +0200 <br>
From: Martin Losch <Martin.Losch@awi.de> <br>
To: MITgcm Support <mitgcm-support@mitgcm.org> <br>
Subject: Re: [MITgcm-support] [EXTERNAL] Re: Query regarding <br>
horizontal strips in currents (Menemenlis, Dimitris (US 329B)) <br>
Message-ID: <A9CE3394-7180-41DA-B3F1-21F1486659B3@awi.de> <br>
Content-Type: text/plain; charset="utf-8" <br>
<br>
Hi Pragnya, <br>
<br>
I have these comments. <br>
1. with your grid resolution of 3.7km and deltaT=90. the dimensionless viscosity parameter viscAhGrid that you would need to get viscAh = 0.7 is 0.7/ (0.25*3700**2 / 90) = 1.84e-05 (see pkg/mom_common/mom_calc_visc.F). I would try something larger, e.g. viscAhGrid = 0.01 or 0.001. I am quite confident that the noise will away with this and you can use hFacMin = 0.1 again. In general I usually use the parameter viscAhGrid instead of vischAh. Also for this high resolution you may want to consider biharmonic viscosity to stabilise the flow, e.g. viscA4Grid = 0.01 ( or smaller ) and then you can use viscAh or viscAhGrid to add some ?physical? viscosity (as small as you like). <br>
<br>
2. vertical grid: as a rule of thumb the vertical grid spacing should not vary abruptly (to preserve order of approximation). I normally use vertical grid spacings that do not exceed the condition dz(k+1)/dz(k) < 1.4. You have biggers jumps of factor 2 and 5. That may also be related to your observation that increasing the hFacMin helps. <br>
<br>
3. Not relevant for your stability, but the combination of no_slip_bottom and a drag parameter is not very physical in my view. The no-slip BC already introduces a drag so that you now add drag at the bottom twice. <br>
<br>
<br>
Martin <br>
> On 15. Jun 2023, at 08:30, Pragnya Makar <asz218003@cas.iitd.ac.in> wrote: <br>
>  <br>
> Hi, <br>
> Thank you for your feedback. Well, I have increased the hFacMin value from 0.1 to 0.7; which has significantly reduced the instability in the current simulation, as shown in the attached image. I will also investigate the lateral boundary conditions as you suggested. Additionally, I will examine the viscosity aspect.Also, please find the attached data file. <br>
>  <br>
> Regards, <br>
> Pragnya Makar <br>
<br>
<br>
<br>
------------------------------ <br>
<br>
Subject: Digest Footer <br>
<br>
_______________________________________________ <br>
MITgcm-support mailing list <br>
MITgcm-support@mitgcm.org <br>
http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support <br>
<br>
<br>
------------------------------ <br>
<br>
End of MITgcm-support Digest, Vol 240, Issue 8 <br>
********************************************** <br>
CAUTION: This email originated from outside of IIT Delhi. Do not click links or open attachments unless you recognize the sender and know the content is safe. <br>
<br>
</p>
</blockquote></div><br></div>