[MITgcm-support] a question about the currents

[Matthew Mazloff]

Hi Fei

How does your wind stress look? That is the first thing I would check.

Matt


> On Jan 6, 2020, at 7:21 PM, 979583702 <979583702 at qq.com> wrote:
> 
> 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
> <currents.png>_______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support