[MITgcm-support] Coriolis force in non-hydrostatic model

Павел Лобовиков plobovikov at gmail.com
Thu Mar 21 07:14:09 EDT 2019


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/20190321/a22040e9/attachment.html>


More information about the MITgcm-support mailing list