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