[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