<div dir="ltr">Kindly reminder :)</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">чт, 21 мар. 2019 г. в 14:14, Павел Лобовиков <<a href="mailto:plobovikov@gmail.com">plobovikov@gmail.com</a>>:<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 dir="ltr"><div dir="ltr">Hi everyone!<br><br>I use non-hydrostatic MITgcm model to simulate the propagation of internal waves in three-layer ocean.<br><br>Now I need to check the impact of Coriolis force to internal waves with it long propagation.<br><br>I try to enable Corolis force using useCoriols=.TRUE. and f0=1.176E-4 (this corresponds to 54 latitude), but i don't see any influence on my internal wave.<br><br>What I'm doing wrong?<br><br>My data file looks like this:<br><br><div># ====================</div><div># | Model parameters |</div><div># ====================</div><div>#</div><div># Continuous equation parameters</div><div> &PARM01</div><div> tRef=200*0.,</div><div> sRef=200*0.,</div><div> viscAh=0.E-4,</div><div> viscAr=0.E-5,</div><div> viscA4=0.E-5,</div><div> no_slip_sides=.FALSE.,</div><div> no_slip_bottom=.FALSE.,</div><div> beta=9.99999E-11,</div><div> eosType='LINEAR',</div><div> tAlpha=9.80392E-4,</div><div> sBeta=0.,</div><div> rhoNil=1020,</div><div> rhoConst=1020,</div><div> gravity=9.8066,</div><div> useCoriolis=.TRUE.,</div><div> f0=1.176E-4,</div><div><br></div><div># Free serface formulations.</div><div> rigidLid=.TRUE.,</div><div># Crank-Nicolson numeric technique.</div><div># implicSurfPress=0.5,</div><div># implicDiv2Dflow=0.5,</div><div># implicitNHPress=1.,</div><div> implicitFreeSurface=.FALSE.,</div><div> nonlinFreeSurf=0,</div><div># hFacInf=0.05,</div><div># hFacSup=2.0,</div><div># Nonhydrostatic On/Off parameters.</div><div> exactConserv=.FALSE.,</div><div> nonHydrostatic=.TRUE.,</div><div># Momentum equations.</div><div> momViscosity=.FALSE.,</div><div> momAdvection=.TRUE.,</div><div> momForcing=.TRUE.,</div><div> useCoriolis=.TRUE.,</div><div> momPressureForcing=.TRUE.,</div><div> momStepping=.TRUE.,</div><div> metricTerms=.FALSE.,</div><div> hFacMin=0.001,</div><div># Diffusion.</div><div> implicitDiffusion=.FALSE.,</div><div> diffKhT=0.E-5,</div><div> diffKrT=0.E-7,</div><div> diffKhS=0.E3,</div><div> diffKrS=0.E-5,</div><div># On/off temp and salt parametries.</div><div> tempStepping=.TRUE.,</div><div> tempAdvection=.TRUE.,</div><div> tempImplVertAdv=.FALSE.,</div><div> tempForcing=.FALSE.,</div><div> saltStepping=.FALSE.,</div><div> saltAdvection=.FALSE.,</div><div> saltImplVertAdv=.FALSE.,</div><div> saltForcing=.FALSE.,</div><div><br></div><div># Bottom linear drag coefficient (def=2.E-4)</div><div> bottomDragLinear=0.E-4,</div><div># Bottom quadratic drag coefficient (def=3.E-3)</div><div> bottomDragQuadratic=0.E-3,</div><div> tempAdvScheme=77,</div><div> staggerTimeStep=.TRUE.,</div><div>#- not safe to use globalFiles in multi-processors runs</div><div> globalFiles=.FALSE.,</div><div> useSingleCpuIO=.FALSE.,</div><div> readBinaryPrec=64,</div><div> writeBinaryPrec=64,</div><div> writeStatePrec=64,</div><div> &</div><div><br></div><div># Elliptic solver parameters</div><div> &PARM02</div><div> cg2dMaxIters=1000,</div><div> cg2dTargetResidual=1.E-13,</div><div> cg2dUseMinResSol=1,</div><div> cg3dMaxIters=400,</div><div> cg3dTargetResidual=1.E-13,</div><div> &</div><div><br></div><div># Time stepping parameters</div><div> &PARM03</div><div># External forcing parameters.</div><div># periodicExternalForcing=.TRUE,</div><div># externForcingPeriod=0.025,</div><div># externForcingCycle=1,</div><div># General time stepping parameters.</div><div> nIter0=0,</div><div> nTimeSteps=400000,</div><div> deltaT=0.25,</div><div> abEps=0.1,</div><div> abEps=0.1,</div><div> pChkptFreq=12519,</div><div> chkptFreq=100,</div><div> dumpFreq=100,</div><div> monitorFreq=100,</div><div> monitorSelect=2,</div><div> &</div><div><br></div><div># Gridding parameters</div><div> &PARM04</div><div> usingCartesianGrid=.TRUE.,</div><div> delXFile='delXvar',</div><div> delYFile='delYvar',</div><div> delRFile='delZvar',</div><div> &</div><div><br></div><div># Input datasets</div><div> &PARM05</div><div> hydrogSaltFile='salinity.init',</div><div> checkIniSalt=.FALSE.,</div><div> hydrogThetaFile='temperature.init',</div><div> checkIniTemp=.FALSE.,</div><div> uVelInitFile='uvelocity.init',</div><div> pSurfInitFile='surface.init',</div><div> &</div><div><br><br>Please, help!<br><br>Thanks, Pavel</div></div></div>
</blockquote></div>