[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