[MITgcm-support] [EXTERNAL] Coriolis force in non-hydrostatic model
Dimitris Menemenlis
dmenemenlis at gmail.com
Fri Mar 22 16:48:49 EDT 2019
Maybe it would help to get an answer if you sent a list of what you actually changed in your data files?
The before and after. Did you just add explicitly default values?
> On Mar 22, 2019, at 12:10 PM, Павел Лобовиков <plobovikov at gmail.com> wrote:
>
> Kindly reminder :)
>
> чт, 21 мар. 2019 г. в 14:14, Павел Лобовиков <plobovikov at gmail.com <mailto: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
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20190322/b4ade088/attachment.html>
More information about the MITgcm-support
mailing list