<div dir="ltr">
<div>Dear Team,</div><div><br></div><div>I am having slight confusion 
about setting the horizontal eddy viscosity in my model. Here, I am only
 considering Laplacian horizontal viscosity (i.e. not considering Leith,
 Biharmonic viscosity, etc). I have a 360x780 grid, with 3000m x 3000m 
spatial resolution, length of mesh ~20 deg, width of mesh ~10 deg and 
timestep of 120 sec.</div><div><br></div><div>Referring to the manual, I
 find two criterias to be satisfied i.e., CFL Constraint on Viscosity 
(Section 2.19.1.7) and Numerical Stability Criteria (Section
 4.3.2.1, eq 4.26). Currently, I am specifying viscAh=5000 m2/s. <br><div><img src="https://mail.google.com/mail/u/0?ui=2&ik=e397a587cd&attid=0.1&permmsgid=msg-a:r5298866202815878677&th=174e917745dffb75&view=fimg&sz=s0-l75-ft&attbid=ANGjdJ9reRgvTG-piwXdubfUz6sOvGI2lqYn2O4C4DFJOMDFzwm4ygP_U5iNiGe_uQPG_XX0pAzkOVEeffpykgnMDXHPIcaX0SZ3WGd_z-UcNGT6I_8gYBcvR2KT0JY&disp=emb&realattid=ii_kfs5alcb0" alt="image.png" style="margin-right: 0px;" class="gmail-CToWUd" width="87" height="46">The upper limit for Ah comes to be 18750 > 5000. Also,</div><div><div><img src="https://mail.google.com/mail/u/0?ui=2&ik=e397a587cd&attid=0.2&permmsgid=msg-a:r5298866202815878677&th=174e917745dffb75&view=fimg&sz=s0-l75-ft&attbid=ANGjdJ-aM9dFME_XVC_RTaTS6mxnAB8TRPHYdPGaLxzle2mO4tx8FbD0xjpPeMXoACq_kvWF8-dZWy2okz7rgAzEUPF8arJerTXjzZlj5JjjRYzBQIk51O3MhbwiGCY&disp=emb&realattid=ii_kfs5ep841" alt="image.png" style="margin-right: 0px;" class="gmail-CToWUd" width="182" height="42">i.e.
 the stability for horizontal viscosity. For my setup, it gave a value 
of 0.53 < 0.6 (solving using Adams-Bashforth). Now, when I use 
viscAh=6000, Sl becomes 0.64 > 0.6, but my model simulation went 
through without crashing. I am trying to understand how the model still 
ran.</div><div><br></div><div>Also, I believe using 
viscAhgrid would be a better option. But the manual doesn't mention how 
viscAhgrid is being computed/ what's the formula (range is between 0 and 1). I have set 
viscAhGridMax=0.5 in data 

</div><div><br></div><div>I am trying to address the following questions:<br></div><div>1. How to compute viscAhgrid? (or do I have to use Sl value as viscAhgrid?)<br></div><div>2.
 Can both viscAz and Implicitviscocity=.TRUE. be set at the same time 
(since KPP needs implicitviscocity=.TRUE). I have kept viscAz=0.001<br></div><div><span style="font-size:12pt;line-height:107%;font-family:"Times New Roman","serif""></span>



</div><div>3. Is there any way to check if the run has been stable by using the U 
and V momentum tendencies from dissipation and advection term in the 
output? Or any other variable in the output (S,T,U,V,W,eta)? <br></div><div>Note: I am specifically interested in using the Laplacian eddy viscosity. <br></div><div>
 I am also attaching my data file. Kindly assist how I can improve the setup further. <br></div><div><br></div><div>Regards</div><div>Kunal</div></div></div>

</div>