<div dir="ltr">Hello, dear colleagues!<br><br>I am configuring the MITgcm for the Kara Sea region. With the<b> curvilinear grid</b>, the model's resolution is: <b>dx=500-1200 m, dz=2 m</b> in upper layers, and up to<b> 20 m</b> in the deepest. The minimum model depth is set equal to <b>5 m</b>. I use the <b>pure r</b> vertical coordinate (not <b>r*</b>) for now. Because of the curvilinear grid, I use <b>vectorinvariant</b> package. Other packages used are: <b>seaice, my82, exf, cal, obcs</b>.<br><br>Having carefully read (once again) the newest version of the online-manual and browsed the code, I have configured almost all necessary model parameters, but still some questions remain. They are qiute specialized, so I could not find the needed information in puplished papers. But for those experienced with MITgcm, they might be easy to answer. I will appreciate any advise!<br><br>1) The newly added feature with partial cells (#undef <b>EXCLUDE_PCELL_MIX_CODE in CPP_OPTIONS</b>). The parameter <b>pCellMix_select = 0</b> by default, which means that the enhanced mixing is OFF in bottom and surface layers. The parameters intended for the enhanced mixing in the inner layers are also OFF by default: <b>interViscAr_pCell</b> and <b>interDiffKr_pCell = FALSE</b>. Are this option and the corresponding run-time parameters really important for such configuration (Kara Sea)? Any experience with them?<br><br>2) I use the following viscocity and diffusivity values, intending to use the <b>Smagorinsky viscocity</b> for lateral turbulent exchange:<br> <b>viscAh = 1.0,<br> viscC2smag = 3.0,<br> viscAr = 1.0e-5,<br> diffKhT = 1.0,<br> diffKhS= 1.0,<br> diffKrT = 1.0e-5,<br> diffKrS = 1.0e-5</b>,<br>But for diffusivity there is no <b>Smagorinsky</b> analogue. What approach can be the most suitable in this case?<br><br>3) In the paper "<u>ECCO version 4 - an integrated framework... " by Forget et al., 2015</u>, the authors decided to switch OFF the <b>C-D sheme</b> and to switch ON the <b>Implicit vertical advection</b> option. But they investigated the global scales, not high-res. From your experience, can their advise be applyed to the high-res ocean simulations?<br><br>4) The <b>MY82</b> package does not allow using the convective adjustment (<b>cAdjFreq</b>) and <b>ivdc_kappa</b>. Does it mean that the MY in MITgcm works fine with vertical instability cases? Actually, it should, as follows from J. Mellor's papers, but this question arose when I was investigating the <b>GGL90</b> code, where I found the line "<i>convective adjustment might be needed even with ggl90</i>"...<br><br>Once again, I will appreciate any advise!<div><br></div><div>PS: Sorry for large wall-of-text here, but I tried to be concrete. :-)</div></div>