[MITgcm-support] Tendency of Salinity

[Mattia Almansi]

Hi all,

I'm trying to reproduce TOTSTEND using the ADV, DF, KPP and SFLUX terms.
I'm following the procedure explained by Abhisek and Jean-Michel here:
http://mitgcm.org/download/daily_snapshot/MITgcm/doc/
Heat_Salt_Budget_MITgcm.pdf

When I plot the difference between my estimation and TOTSTEND at one
specific subsurface layer, it's clear that the error is huge where there is
shallow water (the error is proportional to 1/Depth).
I'm pretty sure that this is due to the non-linear free surface +
r*-coordinate.
I tried to substitute DRF with DRFstar [DRFstar = DRF * H /(H + eta),
right?] , but it doesn't help.

Does anyone know how to fix this?

Thanks,
Mattia

ps: Here is my data file:

 &PARM01
 viscAz =1.E-5,
 diffK4T=0.0,
 diffKzT=1.E-5,
 diffK4S=0.0,
 diffKzS=1.E-5,
 viscAhGridMax=0.25,
 viscC2Leith=0.0,
 viscAhGrid=0.0,
 useAnisotropicViscAGridMax=.TRUE.,
 viscA4GridMax=0.25,
 viscC4Leith=1.0,
 viscA4Grid=0.0,
 no_slip_sides=.TRUE.,
 no_slip_bottom=.TRUE.,
 bottomDragQuadratic = 0.001,
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 useRealFreshWaterFlux=.TRUE.,
#--- z* - begin
 nonlinFreeSurf    = 4,
 exactConserv      = .TRUE.,
 select_rStar      = 2,
 hFacInf           = 0.01,
 hFacSup           = 2.5,
#--- z* - end
#--- - debugMode in eedata
 debugLevel          = 3,
 eosType='JMD95Z',
 readBinaryPrec=64,
 writeBinaryPrec=64,
 useSingleCpuIo=.TRUE.,
 globalFiles=.TRUE.,
 hFacMin=0.1,
 hFacMinDz=5.,
 allowFreezing=.FALSE.,
 HeatCapacity_Cp = 3986.D0,
 gravity         = 9.8156,
 gBaro           = 9.8156,
 rhoNil          = 1027.D0,
 implicitDiffusion=.TRUE.,
 implicitViscosity=.TRUE.,
 vectorInvariantMomentum=.TRUE.,
 staggerTimeStep=.TRUE.,
 tempAdvScheme=33,
 saltAdvScheme=33,
 MultiDimAdvection=.TRUE.,
 &



-- 
Mattia Almansi,
Ph.D. Student,
Johns Hopkins University,
Department of Earth and Planetary Sciences,
315 Olin Hall - 3400 N. Charles Street, Baltimore, MD 21218

e-mail: mattia.almansi at jhu.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20170421/7d35b24d/attachment.html>