[MITgcm-support] Underestimated ETA values for Global Ocean!

Abhisek Chakraborty abhisek.sac at gmail.com
Sat Mar 16 05:31:50 EDT 2013


Dear MITgcmUsers,

I am working with global ocean simulation using MITgcm. The "data" file
configuration is

# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef = 3*23.,3*22.,21.,2*20.,19.,2*18.,17.,2*16.,15.,14.,13.,
                      12.,11.,2*9.,8.,7.,2*6.,2*5.,3*4.,3*3.,4*2.,12*1.,
 sRef = 50*34.5,
#
 no_slip_sides=.FALSE.,
 no_slip_bottom=.TRUE.,
#
 viscAr= 1.77e-5,
 diffKhT=0.,
 diffK4T=0.,
 diffKrT=1.51e-5,
 diffKhS=0.,
 diffK4S=0.,
 diffKrS=1.51e-5,
#
 viscAh=0.,
 viscAhGrid=0.,
 viscC2Leith=0.,
 viscA4=0.,
 useFullLeith=.TRUE.,
 viscC4Leith=2.0,
 viscC4Leithd=2.0,
 useAnisotropicViscAGridMax=.TRUE.,
 bottomDragQuadratic = 0.002,
#
 tAlpha=2.E-4,
 sBeta =7.4E-4,
 rhonil=1035.,
 rhoConstFresh=1000.,
 gravity=9.81,
 eosType='JMD95Z',
 readBinaryPrec=32,
 writeBinaryPrec=32,
 useSingleCPUio=.TRUE.,
 globalFiles=.TRUE.,
 hFacMinDr=50.,
 hFacMin=0.3,
 implicitDiffusion=.true.,
 implicitViscosity=.true.,
 implicitFreeSurface=.TRUE.,
 useRealFreshWaterFlux=.true.,
 exactConserv=.TRUE.,
 allowFreezing=.TRUE.,
 useCDscheme=.TRUE.,
 useNHMTerms=.TRUE.,
#
 tempAdvScheme=33,
 saltAdvScheme=33,
 StaggerTimeStep=.true.,
 multiDimAdvection=.true.,
 vectorInvariantMomentum=.TRUE.,
 rigidLid=.FALSE.,
 debuglevel=-1
&

# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=500,
 cg2dTargetResidual=1.E-13,
 &

# Time stepping parameters
 &PARM03
# nIter0 =     420,
# nTimeSteps = 3180,
# 100 years of integration will yield a reasonable flow field
 startTime  = 0.,
 endTime    = 15552000000.,
 deltaTmom  = 600.0,
 deltaTtracer= 600.0,
 deltaTClock = 86400.,
 deltaTfreesurf= 600.,
 tauCD =     321428.,
 abEps = 0.1,
 pChkptFreq= 3110400000.,
 dumpFreq=   3110400000.,
 taveFreq=   3110400000.,
 monitorFreq=2592000.,
 &

# Gridding parameters
 &PARM04
 usingSphericalPolarGrid=.TRUE.,
 delR=10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.01,
 10.03, 10.11, 10.32, 10.80, 11.76, 13.42, 16.04 , 19.82, 24.85,
 31.10, 38.42, 46.50, 55.00, 63.50, 71.58, 78.90, 85.15, 90.18,
 93.96, 96.58, 98.25, 99.25,100.01,101.33,104.56,111.33,122.83,
 139.09,158.94,180.83,203.55,226.50,249.50,272.50,295.50,318.50,
 341.50,364.50,387.50,410.50,433.50,456.50,
 ygOrigin=-80.0,
 dySpacing=1.,
 dxSpacing=1.,
 &

# Input datasets
 &PARM05
 bathyFile=      '../input/bathymetry.bin',
 hydrogThetaFile='../input/lev_t.bin',
 hydrogSaltFile= '../input/lev_s.bin',
 &

and the atmospheric forcings
(atemp,aqh,uwind,vwind,swdown,lwdown,precip,runof,climsst,climsss) are
specified through EXF package. Apart from these I am using GMRedi and KPP
parametrisation schemes. After spinning up the model for 500 years and then
running the model with monthly NCEP fields, I found that temperature,
salinity, currents are matching well with GODAS (NCEP), but surface
elevation (ETA) is underestimated (though the patterns are matching with
GODAS). I have attached one simulation result ( you can see negative
elevation in most of the places). Can anybody please guide me how to
correct this?

Thanking in advance,
Abhisek
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20130316/a46d9916/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Untitled_eta.gif
Type: image/gif
Size: 317395 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20130316/a46d9916/attachment-0001.gif>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: data.pkg
Type: application/octet-stream
Size: 108 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20130316/a46d9916/attachment-0004.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: data.gmredi
Type: application/octet-stream
Size: 1208 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20130316/a46d9916/attachment-0005.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: data.kpp
Type: application/octet-stream
Size: 182 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20130316/a46d9916/attachment-0006.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: data.exf
Type: application/octet-stream
Size: 9916 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20130316/a46d9916/attachment-0007.obj>


More information about the MITgcm-support mailing list