[MITgcm-support] Help me get rid of my noise!

Ryan Abernathey ryan.abernathey at gmail.com
Wed May 7 13:11:59 EDT 2014


Hello MITgcm Support Group,

I am running a 2D model (depth and latitude) in a Southern Ocean
configuration. I am trying to eliminate what appears to be unphysical
grid-scale noise in my W field--columns of alternately very high and low W.
This is a snapshot of W after 5 years of spinup. This pattern is steady in
time.
As you can see, the noise is concentrated mainly (but not exclusively) over
the shelf break. A possibly related issue is that my cg2d solver is
relatively far from convergence:
cg2d: Sum(rhs),rhsMax =  -5.46186194609888E-05  6.64557830841732E-02

My model is 0.5 degree resolution in latitude, meaning that it is
relatively high resolution compared to the tutorials (e.g. ideal_2D_oce),
but of course it has no eddies or any temporal variability at all, since it
is 2D. This puts it in an unknown [at least to me] parameter regime.

I have tried cranking the lateral viscosity (both laplacian and biharmonic)
up as high as they will go. I have enabled both linear and quadratic bottom
drag. I am NOT using the CD scheme; the noise was much worse with it
enabled. I am NOT using KPP but rather GGL90. (I tried with KPP and it was
not very different.)

I would appreciate any suggestions you have about how to get rid of this
noise and achieve smooth fields.

Thanks,
Ryan

####### data file ##########
 &PARM01
 useAreaViscLength=.FALSE.,
 viscAr=1.0E-03,
 viscA4=1e12,
 viscAh=1e5,
 tempAdvScheme=77,
# KrT fields are read in from a file
 diffKrT=0.0E-5,
 diffKhT=0.0,
 diffK4T=0.0,
 staggerTimeStep=.TRUE.,
# initial vertical profiles of T and S
 tRef=44*0.0,
 sRef=44*35.0,
# equation of state
 eosType='LINEAR',
 rhonil=1035.,
 eosType='LINEAR',
 tAlpha=2.000000E-04,
 sBeta=0.000000E+00,
 saltStepping=.FALSE.,
# boundary conditions
 no_slip_sides=.FALSE.,
 no_slip_bottom=.TRUE.,
# needed for points above the U sponge layer (shelf)
 bottomDragLinear=1.E-3,
 bottomDragQuadratic = 0.0012
# additional drag is done with RBCS
# physical parameters
 gravity=9.810000E+00,
# implicit diffusion and convective adjustment
 implicitDiffusion=.TRUE.,
 implicitViscosity=.TRUE.,
 implicitFreeSurface=.TRUE.,
# exact volume conservation
# exactConserv=.TRUE.,
# C-V scheme for Coriolis term
 useCDscheme=.FALSE.,
# partial cells for smooth topography
 hFacMin=5.000000E-02,
# file IO stuff
 readBinaryPrec=64,
 useSingleCpuIO=.TRUE.,
 debugLevel=1,
 &
# elliptic solver parameters
 &PARM02
 cg2dMaxIters=500,
 cg2dTargetResidual=1E-09,
 &
# timestepping parameters
 &PARM03
 nIter0=0000172800,
# nTimeSteps=1000,
# nTimeSteps=12961,
 nTimeSteps=345600,
# deltaT=3600.0,
# deltaT=2400.,
# deltaT=1800.0,
 deltaT=900.0,
 abEps=0.1,
# needed for sea ice
 forcing_In_AB = .FALSE.,
 ChkptFreq=31104000.
 pChkptFreq=311040000.,
 taveFreq=311040000.0,
 dumpFreq=31104000.,
 monitorFreq=2592000.0,
 tauThetaClimRelax=2592000.0,
 dumpInitAndLast=.TRUE.,
 pickupStrictlyMatch=.FALSE.,
 &
# gridding parameters
 &PARM04
 usingSphericalPolarGrid=.TRUE.,
 delX=0.5,
 delY=100*0.5,
 xgOrigin=0.,
 ygOrigin=-80.,
 delR=10., 10., 10., 10., 10., 10., 10., 12., 14.,
  16.,   19.,   22.,   26.,   30.,   36.,   42.,   50.,   60.,
  72.,   85.,  100.,  120.,  140.,  166.,  200.,  200.,  200.,
  200.,  200.,  200.,  200.,  200.,  200.,  200.,  200.,  200.,
  200.,  200.,  200.,  200.,  200.,  200.,  200.,  200.
 &
# Input datasets
 &PARM05
 bathyFile='bathyFile_shelf_2Dchan.bin',
 zonalWindFile='zonalWindFile.bin',
 thetaClimFile='thetaClimFile_WOA98_2Dchan.bin',
 diffKrFile='diffKrFile_nosponge_2Dchan.bin',
 hydrogThetaFile='hydrogThetaFile_WOA98_2Dchan.bin',
 hydrogSaltFile='hydrogSaltFile_WOA98_2Dchan.bin',
 &

########## data.pkg ##############
 &PACKAGES
  useMNC         = .TRUE.,
  useGMREDI      = .TRUE.,
  useGGL90       = .TRUE.,
  useRBCS        = .TRUE.,
 &

######## data.gmredi ##############
 &GM_PARM01
  GM_MNC             = .FALSE,
  GM_AdvForm         = .TRUE.,
  GM_UseBVP          = .TRUE.,
  GM_background_K    = 2000,
  GM_isopycK         = 2000,
  GM_Kmin_horiz      = 50,
  GM_taper_scheme    = 'dm95',
  GM_BVP_ModeNumber  = 1,
  GM_BVP_cMin        = .1,
 &end

########### data.ggl90 ############
 &GGL90_PARM01
 GGL90writeState=.FALSE.,
 GGL90TKEmin = 1.e-7,
 mxlMaxFlag=2,
 GGL90mixingLengthMin=3.,
 &


############ SIZE.h ##############
...
      PARAMETER (
     &           sNx =   1,
     &           sNy = 100,
     &           OLx =   4,
     &           OLy =   4,
     &           nSx =   1,
     &           nSy =   1,
     &           nPx =   1,
     &           nPy =   1,
     &           Nx  = sNx*nSx*nPx,
     &           Ny  = sNy*nSy*nPy,
     &           Nr  =  44)
...
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20140507/96b6a075/attachment-0001.htm>


More information about the MITgcm-support mailing list