<div>Hi everyone,</div><div><br></div><div>Recently, I'd like to study the relationship between MOC and ocean vertical mixing. I want to change the strength of vertical mixing, as far as I know, we can change the two parameters:  diffKrT and  diffKrS (the value typically ranges from 1.E-5 to 1.E-4) in the data file. In my experiments, I changed diffKrT and  diffKrS from 1.E-5 to 5.E-5, but I found there is little change in the MOC strength. It seems quite odd.</div><div><br></div><div>In addition, I didn't switch on the KPP package in the data.pkg. (I comment out the line #useKPP = .TRUE. in data.pkg)</div><div><br></div><div>I was wondering if anyone could tell me the reason why MOC changes little when I change the  diffKrT and  diffKrS.  Is there any better way to change the vertical mixing strength (e.g., switch on KPP package and change the vertical mixing strength)?  If there is, how to change?</div><div><br></div><div>Thank you so much.</div><div><br></div><div>Best,</div><div><br></div><div>Fei</div><div><br></div><div><br></div><div><br></div><div>Here is my data file:</div><div><div># ====================</div><div># | Model parameters |</div><div># ====================</div><div>#</div><div># Continuous equation parameters</div><div> &PARM01</div><div> tRef               = 3*23.,3*22.,21.,2*20.,19.,2*18.,17.,2*16.,15.,14.,13.,</div><div>                      12.,11.,2*9.,8.,7.,2*6.,2*5.,3*4.,3*3.,4*2.,12*1.,</div><div> sRef               = 50*34.5,</div><div> no_slip_sides  = .TRUE.,</div><div> no_slip_bottom = .TRUE.,</div><div>#</div><div> viscAr=0.5E-4,</div><div>#</div><div> viscAh=1.E0,</div><div> viscAhGrid=2.E-2,</div><div># viscAh=2.0e4,</div><div>#</div><div> diffKhT=1.E1,</div><div> diffKhS=1.E1,</div><div># change the vertical mixing parameters from 1.E-5 to 5.E-5</div><div> diffKrT=1.E-5,</div><div> diffKrS=1.E-5,</div><div>#</div><div>### diffKrBL79surf=0.1E-4,</div><div>### diffKrBL79deep=1.0E-4,</div><div> bottomDragQuadratic = 0.001,</div><div>#when using ggl90</div><div> ivdc_kappa=10.,</div><div> implicitDiffusion=.TRUE.,</div><div> implicitViscosity=.TRUE.,</div><div> useRealFreshWaterFlux=.TRUE.,</div><div># balanceThetaClimRelax=.TRUE.,</div><div> balanceSaltClimRelax=.TRUE.,</div><div># balanceEmPmR=.TRUE.,</div><div># balanceQnet=.TRUE.,</div><div> allowFreezing=.FALSE.,</div><div>### hFacInf=0.2,</div><div>### hFacSup=2.0,</div><div> hFacMin=.2,</div><div> hFacMinDr=5.,</div><div> select_rStar=2,</div><div> nonlinFreeSurf=4,</div><div> gravity=9.81,</div><div> rhonil=1029.,</div><div> rhoConst=1029.,</div><div> rhoConstFresh=1000.,</div><div> convertFW2Salt=-1.,</div><div> eosType='JMD95Z',</div><div> implicitFreeSurface=.TRUE.,</div><div> exactConserv=.TRUE.,</div><div> useSingleCpuIO=.TRUE.,</div><div> tempAdvScheme=30,</div><div> saltAdvScheme=30,</div><div> tempVertAdvScheme=3,</div><div> saltVertAdvScheme=3,</div><div> tempImplVertAdv=.TRUE.,</div><div> saltImplVertAdv=.TRUE.,</div><div> staggerTimeStep=.TRUE.,</div><div> vectorInvariantMomentum=.TRUE.,</div><div>#when using the cd scheme:</div><div># useCDscheme=.TRUE.,</div><div> useJamartWetPoints=.TRUE.,</div><div> readBinaryPrec=32,</div><div> writeBinaryPrec=32,</div><div> debugLevel=1,</div><div> /</div><div><br></div><div># Elliptic solver parameters</div><div> &PARM02</div><div> cg2dMaxIters=300,</div><div>#cg2dTargetResWunit=1.E-12,</div><div> /</div><div><br></div><div># Time stepping parameters</div><div> &PARM03</div><div> nIter0=0,</div><div>#2 lev2 for testing:</div><div>#nTimeSteps=8,</div><div> nTimeSteps=350000,</div><div>#</div><div> forcing_In_AB=.FALSE.,</div><div> momDissip_In_AB=.FALSE.,</div><div>#when using the cd scheme:</div><div># epsAB_CD = 0.25,</div><div># tauCD=172800.0,</div><div> deltaTmom   =3600.,</div><div> deltaTtracer=3600.,</div><div> deltaTfreesurf=3600.,</div><div> deltaTClock =3600.,</div><div>#when using ab2:</div><div># abEps = 0.1,</div><div>#when using ab3:</div><div> doAB_onGtGs=.FALSE.,</div><div> alph_AB=0.5,</div><div> beta_AB=0.281105,</div><div>#</div><div> pChkptFreq  =10368000.0,</div><div> chkptFreq   =10368000.0,</div><div># taveFreq    =2635200.0,</div><div># dumpFreq    =2635200.0,</div><div># monitorFreq = 7200.0,</div><div> monitorFreq = 2592000.0,</div><div> dumpInitAndLast = .TRUE.,</div><div> adjDumpFreq = 31536000.0,</div><div>#adjDumpFreq = 604800.0,</div><div> adjMonitorFreq = 864000.0,</div><div> pickupStrictlyMatch=.FALSE.,</div><div> /</div><div><br></div><div># Gridding parameters</div><div> &PARM04</div><div> usingCurvilinearGrid=.TRUE.,</div><div> delR = </div><div>     10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.01,</div><div>     10.03, 10.11, 10.32, 10.80, 11.76, 13.42, 16.04, 19.82, 24.85,</div><div>     31.10, 38.42, 46.50, 55.00, 63.50, 71.58, 78.90, 85.15, 90.18,</div><div>     93.96, 96.58, 98.25, 99.25,100.01,101.33,104.56,111.33,122.83,</div><div>     139.09,158.94,180.83,203.55,226.50,249.50,272.50,295.50,318.50,</div><div>     341.50,364.50,387.50,410.50,433.50,456.50,</div><div> /</div><div><br></div><div># Input datasets</div><div> &PARM05</div><div> diffKrFile='total_diffkr_r009bit11.bin',</div><div> adTapeDir='tapes',</div><div>#bathyFile      ='bathy_eccollc_90x50.bin',</div><div> bathyFile      ='bathy_eccollc_90x50_min2pts.bin',</div><div> hydrogThetaFile='T_OWPv1_M_eccollc_90x50.bin',</div><div> hydrogSaltFile ='S_OWPv1_M_eccollc_90x50.bin',</div><div> viscA4Dfile    ='fenty_biharmonic_visc_v11.bin',</div><div> viscA4Zfile    ='fenty_biharmonic_visc_v11.bin',</div><div> geothermalFile='geothermalFlux.bin',</div><div>#</div><div> /</div></div><div><br></div>