<div dir="ltr"><div>Hello Dimitris,</div><div><br></div><div>I am using SODA (T,S,U and V) as my boundary conditions and EXF forcing from ERA5.</div><div><br></div><div>Recently, Martin suggested in a query, not to use leith and smag simultaneously. Guess, I need to change that.</div><div><br></div><div>The eddies are becoming too extreme by half of the year, dismantling the circulation.</div><div>Thanks,</div><div>Kunal<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Mar 14, 2022 at 3:44 AM Menemenlis, Dimitris (US 329B) <<a href="mailto:dimitris.menemenlis@jpl.nasa.gov" target="_blank">dimitris.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">What type of lateral boundary conditions are you using?<br>
<br>
> On Mar 11, 2022, at 5:22 AM, kunal madkaiker <<a href="mailto:kunal.madkaiker02@gmail.com" target="_blank">kunal.madkaiker02@gmail.com</a>> wrote:<br>
> <br>
> Dear All,<br>
> <br>
> I am using MITgcm for regional modelling (660x480 gridpoints, 5km resolution, delZ=5m), over NIO. I simulated the model for 1 year (warm start), with 3-sided BC and daily EXF forcing.<br>
> <br>
> I am observing more turbulence, and it increases as the simulation progresses, especially in the latter half of the year. I am also turning KPP on, and using Leith and Smagorinsky viscosities.<br>
> <br>
> I am attaching a comparison plot of model surface currents with observation. Can anyone suggest to me why this is happening?<br>
> <br>
> Here is my data file:<br>
> # ====================<br>
> # | Model parameters |<br>
> # ====================<br>
> #<br>
> # Continuous equation parameters<br>
>  &PARM01<br>
>  viscAh=100.,<br>
>  viscAz=1.E-5,<br>
>  viscA4=0.,<br>
>  diffKhT=0.1,<br>
>  diffKzT=1.E-5,<br>
>  diffK4T=0.,<br>
>  diffKhS=0.1,<br>
>  diffKzS=1.E-5,<br>
>  diffK4S=0.,<br>
> #<br>
>  useFullLeith=.TRUE.,<br>
>  useStrainTensionVisc=.TRUE.,<br>
> #<br>
>  viscC2leith=0.5,<br>
>  viscC4leith=0.05,<br>
> #<br>
>  viscC2smag=0.5,<br>
>  viscC4Smag=0.05,<br>
> #<br>
>  viscAhGridMax=1.E-2,<br>
>  viscA4GridMax=1.E-9,<br>
> #<br>
>  no_slip_sides=.FALSE.,<br>
>  no_slip_bottom=.TRUE.,<br>
>  tempAdvScheme=33,<br>
>  saltAdvScheme=33,<br>
>  staggerTimeStep=.TRUE.,<br>
> #<br>
>  tAlpha=2.E-4,<br>
>  sBeta =7.E-4,<br>
>  gravity=9.81,<br>
>  f0=0.0000101456,<br>
>  beta=2.2209E-11,<br>
>  rigidLid=.FALSE.,<br>
>  implicitFreeSurface=.TRUE.,<br>
>  implicitDiffusion=.TRUE.,<br>
>  implicitViscosity=.TRUE.,<br>
>  nonHydrostatic=.FALSE.,<br>
>  bottomDragLinear=0.0002,<br>
>  eosType='MDJWF',<br>
>  useRealFreshWaterFlux=.TRUE.,<br>
>  exactConserv=.TRUE.,<br>
>  hFacMin=0.5,<br>
> # hFacMinDz=5.,<br>
>  hFacInf=0.2,<br>
>  hFacSup=2.0,<br>
>  useSingleCpuIO=.TRUE<br>
>  readBinaryPrec=64,<br>
>  writeBinaryPrec=64,<br>
> <br>
>  &<br>
> <br>
> # Elliptic solver parameters<br>
>  &PARM02<br>
>  cg2dMaxIters=1000,<br>
>  cg2dTargetResidual=1.E-13,<br>
>  cg3dMaxIters=400,<br>
>  cg3dTargetResidual=1.E-13,<br>
> <br>
>  &<br>
> <br>
> # Time stepping parameters<br>
>  &PARM03<br>
>  startTime=0.,<br>
>  endTime=31536000.,<br>
>  deltaTmom=180.0,<br>
>  pChkptFreq=31536000., <br>
>  dumpFreq=432000., <br>
> # taveFreq=864000.,<br>
>  &<br>
> <br>
> # Gridding parameters<br>
>  &PARM04<br>
>  usingCurvilinearGrid=.FALSE.,<br>
>  usingCartesianGrid=.TRUE.,<br>
>  usingSphericalPolarGrid=.FALSE.,<br>
>  dxSpacing=5055.960,<br>
>  dySpacing=5100.408,<br>
>  delZ=  5.0 , 5.0 , 5.0 , 5.0 , 5.0, 5.0 , 5.0 , 5.0 , 5.0 , 5.0,<br>
>         5.0 , 5.0 , 5.0 , 5.0 , 5.0, 5.0 , 5.0 , 5.0 , 5.0 , 5.0,<br>
>         5.0 , 5.0 , 5.0 , 5.0 , 5.0, 5.0 , 5.0 , 5.0 , 5.0 , 5.0,<br>
>         10.0 , 10.0 , 10.0 , 10.0 , 10.0 , 10.0 , 10.0 , 10.0 , 10.0 , 10.0,<br>
>         50.0 , 100.0 , 100.0 , 500.0 , 500.0 , 500.0 , 1000.0 , 1000.0 , 1000.0,<br>
> <br>
>  radius_fromHorizGrid=6370.E3,<br>
>  &<br>
> <br>
> # Input datasets<br>
>  &PARM05<br>
>  hydrogThetaFile='temp_dec_init_5km.bin',<br>
>  hydrogSaltFile='salt_dec_init_5km.bin',<br>
>  uVelInitFile='U_dec_init_5km.bin',<br>
>  vVelInitFile='V_dec_init_5km.bin',<br>
>  bathyFile='bathy_BoCo_5km_UPD.bin',<br>
>  checkIniTemp=.false.,<br>
>  checkIniSalt=.false.,<br>
>  &<br>
> <br>
> Thanks,<br>
> Kunal<br>
> <2x4_Curr_seasonal.png>_______________________________________________<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!cOlg4J_PViRrtv36OlmpkcycbEDYcMJdNfroCGUsrk_n7lO9KRSg29T8w-1a5sHzcJ2eQfdM6zA$" rel="noreferrer" target="_blank">https://urldefense.us/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!PvBDto6Hs4WbVuu7!cOlg4J_PViRrtv36OlmpkcycbEDYcMJdNfroCGUsrk_n7lO9KRSg29T8w-1a5sHzcJ2eQfdM6zA$</a> <br>
<br>
_______________________________________________<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>