[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