[MITgcm-support] non-hydrostatic momentum budget

Andrea Cimatoribus andrea.cimatoribus at epfl.ch
Thu Jul 6 11:30:15 EDT 2017


Dear all,
I am trying to compute the budget of momentum in my non-hydrostatic 
configuration. First of all, thanks to all those who contributed in the 
past on this topic, it certainly saved me huge amounts of time and 
head-scratching.

I put here the recipe from JMC for reference:
[du/dt] = - gravity * ( ETAN(i) - ETAN(i-1) ) / DXC
         + Um_dPHdx
         + Um_Advec   (+ Um_Cori : but only if using CD-Scheme )
         + Um_Diss    (+ Implicit vertical viscosity tendency)
         + Um_Ext
         + AB_gU

The problem:
the recipe above works beautifully (relative error <<1e-3) in the 
hydrostatic case, with only some issues near the boundary, I guess 
because of my cheap ETA derivatives. However, in the nonhydrostatic 
case, the above budget gives relative errors around 10% in many places 
(10% of the "truth" diagnosed via TOTUTEND/TOTVTEND).

What term am I missing? Could it be that TOTUTEND does not have all the 
terms in the NH case? I tried looking into mom_vecinv.F, but I could not 
reach any conclusion. I also tried playing with ADVrE_Um/ADVrE_Vm fluxes 
(treating it as the implicit vertical viscosity term), but that didn't 
help either (I thought the vertical advection was missing).

Thank you for your insights,
Andrea


PS: I attach the data file and a figure for the curious ones, it shows 
the relative error of the diagnosed acceleration (magnitude of vector) 
with respect to the one directly produced via TOTUTEND/TOTVTEND. There 
seems to be some correlation (not perfect, though) between the relative 
error and the advective term.

-- 
Andrea Cimatoribus
postdoctoral researcher
EPFL ENAC IIE ECOL
https://people.epfl.ch/andrea.cimatoribus
-------------- next part --------------
# Model parameters
# Continuous equation parameters
 &PARM01
 tRef=5.4,
 sRef=0.03,
# Horizontal viscosity / diffusivity
# see http://mitgcm.org/pipermail/mitgcm-support/2010-January/006439.html
 viscAh=1.0E-6,
 diffKhT=2.0E-7,
 viscC2Smag=0.3,
 viscA4Grid=0.01,
 usejamartwetpoints=.TRUE.,
 vectorInvariantMomentum=.TRUE.,
 highOrderVorticity=.TRUE.,
###
# Vertical viscosity / diffusivity
 viscAz=5.0E-5,
 diffKzT=2.0E-7,
###
 f0=1.05e-4,
 beta=0.E-11,
 selectCoriMap=0,
 useEnergyConservingCoriolis=.TRUE.,
#tAlpha=2.E-4,
#sBeta =0.,
 gravity=9.81,
# free surface
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 exactConserv=.TRUE.,
 nonHydrostatic=.TRUE.,
 hFacMin=0.2,
 hFacInf=0.2,
 hFacSup=2.0,
 nonlinFreeSurf=0,
 select_rStar=0,
 useRealFreshWaterFlux=.FALSE.,
# for avoiding dumping of surface waves
# unconditionally stable
# http://mitgcm.org/sealion/online_documents/node40.html
 implicSurfPress=0.5,
 implicDiv2DFlow=0.5,
###
 implicitViscosity=.TRUE.,
 implicitDiffusion=.TRUE.,
# drag
 allowFreezing=.FALSE.,
 no_slip_sides=.TRUE.,
 no_slip_bottom=.TRUE.,
 bottomDragQuadratic=2.5E-3,
#sideDragFactor=2,
 eosType='MDJWF',
 readBinaryPrec=64,
# SOM, see Hill 2012
 tempAdvScheme=81,
# do not use salt for now
 saltStepping=.FALSE.,
 staggerTimeStep=.TRUE.,
# IO configuration
 useSingleCpuIO=.TRUE.,
 globalFiles=.FALSE.,
 &

# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=1000,
 cg2dTargetResidual=1.E-8,
 cg3dMaxIters=1000,
 cg3dTargetResidual=1.E-8,
 &

# Time stepping parameters
 &PARM03
 nIter0=0,
 endTime=2592000.0,
 deltaT=6.0,
 abEps=0.1,
 pChkptFreq=864000.0,
 chkptFreq=0.0,
 dumpFreq=0,
 monitorSelect=1,
 monitorFreq=43200.,
 &

# Gridding parameters
 &PARM04
 usingCartesianGrid=.TRUE.,
 delX=560*113.0,
 delY=196*113.0,
 delZ= 0.30, 0.54, 0.78, 1.02, 1.26, 1.50, 1.74, 1.98, 2.22, 2.46,
       2.70, 2.94, 3.18, 3.42, 3.66, 3.90, 4.14, 4.38, 4.62, 4.86,
       5.10, 5.34, 5.58, 5.82, 6.06, 6.30, 6.54, 6.78, 7.02, 7.26,
       7.50, 7.74, 7.98, 8.22, 8.46, 8.70, 8.94, 9.18, 9.42, 9.66,
       9.90,10.14,10.38,10.62,10.86,11.10,11.34,11.58,11.82,12.06,
 &

# Input datasets
 &PARM05
 bathyFile='lemanbathy.bin',
 hydrogThetaFile='temp_2016-06-05T1200.bin',
 hydrogSaltFile='',
 uVelInitFile='',
 vVelInitFile='',
 zonalWindFile=,
 meridWindFile=,
 &
-------------- next part --------------
A non-text attachment was scrubbed...
Name: momentum_budget_imbalance_ref.png
Type: image/png
Size: 319678 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20170706/6f32f5d9/attachment-0001.png>


More information about the MITgcm-support mailing list