<div dir="ltr"><div>Thank you Dimitris. This helps.</div><div><br></div><div>Kunal<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Oct 12, 2020 at 8:32 PM Dimitris Menemenlis <<a href="mailto:menemenlis@jpl.nasa.gov" target="_blank">menemenlis@jpl.nasa.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div>You can use your sensitivity experiments to select an optimized set of horizontal and vertical viscosities, and other empirical model parameters, for your model configuration as is done here: <a href="https://journals.ametsoc.org/mwr/article/133/5/1224/67525/Using-Green-s-Functions-to-Calibrate-an-Ocean" target="_blank">https://journals.ametsoc.org/mwr/article/133/5/1224/67525/Using-Green-s-Functions-to-Calibrate-an-Ocean</a><div><br><div><blockquote type="cite"><div>On Oct 11, 2020, at 8:36 PM, kunal madkaiker <<a href="mailto:kunal.madkaiker02@gmail.com" target="_blank">kunal.madkaiker02@gmail.com</a>> wrote:</div><br><div><div dir="ltr"><div>Dear All,</div><div><br></div><div>I need some help in setting up my model viscosities. Few model studies done in Arabian Sea and N. Indian ocean suggest that horizontal viscosity can be in the range of 5.E3 to 5.E4 m2/s whereas vertical viscosity can be set as 1.E-3 m2/s in this region. (My model setup has delT=120s, delX=delY=3000m)<br></div><div><br></div><div>After going through the formulae in mom_calc_visc.F, I found 'viscAh=5.E3' computes to be 'viscAhgrid=0.26' (for consideration, selecting viscAhgrid=0.2), which satisfies stability and CFL criteria (Manual, Section 4.3.2.1, eq. 4.24, 4.26). Also, using viscAz=0.03, which satisfies eq. 4.28 for min delZ of 5m. I find this setup to be most stable (KE oscillations are the least) <br></div><div><br></div><div>Now, manual suggests to setup viscAz=0.003, which is 1 order less than my computed value. Also, considering Munk eq. 4.27, ratio of (Ah/Mw^3) computes to be in order of 4x10^-13. So, for Mw=3000 m (delx/dely), Ah becomes 0.0127 m2/s, which is again an order less.</div><div><br></div><div>Can anyone shed some light on this? Incase second combination is more robust, I feel viscAhgrid value will be too small (Ah itself is in order of 1.E-2)<br></div><div></div><div><br></div><div>After doing some preliminary experiments, I observed the below mentioned points:</div><div></div><div>1. Increasing vertical viscosity definitely subdues oscillations in the KE, but same may not be said for horizontal viscosity.<br>2. Reducing vertical viscosity is making the surface currents much stronger.</div><div><br></div><div>Kindly assist.</div><div><br></div><div>Regards</div><div>Kunal<br></div><div><span style="font-size:12pt;line-height:107%;font-family:"Times New Roman","serif""><span></span></span>





</div></div>
_______________________________________________<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.us/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!PvBDto6Hs4WbVuu7!cSU37kr2e5tA0BLBscou8LPfdvZLMuGnjpavof13QouQFAiWDFK25nYik94t2aOxnffB0bhRHH0$" target="_blank">https://urldefense.us/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!PvBDto6Hs4WbVuu7!cSU37kr2e5tA0BLBscou8LPfdvZLMuGnjpavof13QouQFAiWDFK25nYik94t2aOxnffB0bhRHH0$</a> <br></div></blockquote></div><br></div></div>_______________________________________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
</blockquote></div>