<div dir="ltr"><div>Dear All,</div><div><br></div><div>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.</div><div><br></div><div>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></div><div><br></div><div>I am attaching a comparison plot of model surface currents with observation. Can anyone suggest to me why this is happening?</div><div><br></div><div>Here is my data file:</div><div># ====================<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> &</div><div><br></div><div>Thanks,</div><div>Kunal<br></div></div>