[MITgcm-support] a question about the currents

[979583702]

Hi everyone,


Recently, I ran a global ocean model configuration which is similar with the global_with_exf but with 1 grid*1 grid resolution. Strangely, I found the currents  mainly located on the equator (e.g., strong westward currents over eastern equatorial Pacific), and there are no western boundary currents and  ACC. 


Attached are the currents at the time of 40 years.


Here is my data file:


# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef= 13.72 , 12.82 , 11.57 , 10.35 , 9.38 , 
       8.62 , 7.99 , 7.45 , 6.91 , 6.45 , 
       5.99 , 5.63 , 5.26 , 4.97 , 4.68 , 
       4.44 , 4.20 , 4.06 , 3.85 , 3.68 , 
       3.51 , 3.22 , 2.99 , 2.80 , 2.62 , 
       2.47 , 1.93 , 1.61 , 1.27 , 1.23 ,
 sRef = 30*34.,     
 viscAr=1.E-4,
# viscAh=2.E2,
 viscAh=1.E0,
 viscAhGrid=2.E-2,
 diffKhT=0.,
 diffKrT=2.E-5,
 diffKhS=0.,
 diffKrS=2.E-5,
 rhonil=1035.,
 rhoConstFresh=1000.,
 eosType = 'JMD95Z',
 ivdc_kappa=100.,
 implicitDiffusion=.TRUE.,
 allowFreezing=.TRUE.,
 exactConserv=.TRUE.,
 useRealFreshWaterFlux=.TRUE.,
 useCDscheme=.TRUE.,
# turn on looped cells
 hFacMin=.05,
 hFacMindr=50.,
 useSingleCpuIO=.TRUE.,
# set precision of data files
 readBinaryPrec=32,

 &


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


# Time stepping parameters
 &PARM03
 nIter0 =      0,
 nTimeSteps = 864000,
# 100 years of integration will yield a reasonable flow field
# startTime  =          0.,
# endTime    = 3110400000.,
 deltaTmom = 300.0,
 tauCD =     321428.,
 deltaTtracer= 1800.0,
 deltaTClock = 3600.0,
# if you are using a version later than checkpoint45d on the main branch
# you can uncomment the following line and increase the time step
# deltaTtracer and deltaTClock to 172800.0 as well to speed up the
# asynchronous time stepping
# deltaTfreesurf = 172800.0,
 abEps = 0.1,
 pChkptFreq= 31104000.,
# dumpFreq=   311040000.,
# dumpFreq=   864000.,
# taveFreq=   311040000.,
 taveFreq=   2592000.,
 monitorFreq=720.,
 &


# Gridding parameters
 &PARM04
 usingSphericalPolarGrid=.TRUE.,
 delR= 10., 15., 20., 30., 40., 
      50., 50., 60., 60., 60.,
      60., 60., 60., 60., 60.,
      60., 60., 60., 70., 70.,
      70., 100., 140., 200., 250.,
      320., 420., 550., 650., 800.,
 ygOrigin=-80.,
 xgOrigin=0.,
 dySpacing=1.,
 dxSpacing=1.,
 &


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

Thank you so much.
Best,
Fei
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20200107/6c5bb168/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: currents.png
Type: application/octet-stream
Size: 95039 bytes
Desc: not available
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20200107/6c5bb168/attachment-0001.obj>