[MITgcm-support] how to change the vertical mixing strength

979583702 979583702 at qq.com
Fri Dec 20 09:33:35 EST 2019


Hi everyone,


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.


In addition, I didn't switch on the KPP package in the data.pkg. (I comment out the line #useKPP = .TRUE. in data.pkg)


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?


Thank you so much.


Best,


Fei






Here is my data file:
# ====================
# | 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  = .TRUE.,
 no_slip_bottom = .TRUE.,
#
 viscAr=0.5E-4,
#
 viscAh=1.E0,
 viscAhGrid=2.E-2,
# viscAh=2.0e4,
#
 diffKhT=1.E1,
 diffKhS=1.E1,
# change the vertical mixing parameters from 1.E-5 to 5.E-5
 diffKrT=1.E-5,
 diffKrS=1.E-5,
#
### diffKrBL79surf=0.1E-4,
### diffKrBL79deep=1.0E-4,
 bottomDragQuadratic = 0.001,
#when using ggl90
 ivdc_kappa=10.,
 implicitDiffusion=.TRUE.,
 implicitViscosity=.TRUE.,
 useRealFreshWaterFlux=.TRUE.,
# balanceThetaClimRelax=.TRUE.,
 balanceSaltClimRelax=.TRUE.,
# balanceEmPmR=.TRUE.,
# balanceQnet=.TRUE.,
 allowFreezing=.FALSE.,
### hFacInf=0.2,
### hFacSup=2.0,
 hFacMin=.2,
 hFacMinDr=5.,
 select_rStar=2,
 nonlinFreeSurf=4,
 gravity=9.81,
 rhonil=1029.,
 rhoConst=1029.,
 rhoConstFresh=1000.,
 convertFW2Salt=-1.,
 eosType='JMD95Z',
 implicitFreeSurface=.TRUE.,
 exactConserv=.TRUE.,
 useSingleCpuIO=.TRUE.,
 tempAdvScheme=30,
 saltAdvScheme=30,
 tempVertAdvScheme=3,
 saltVertAdvScheme=3,
 tempImplVertAdv=.TRUE.,
 saltImplVertAdv=.TRUE.,
 staggerTimeStep=.TRUE.,
 vectorInvariantMomentum=.TRUE.,
#when using the cd scheme:
# useCDscheme=.TRUE.,
 useJamartWetPoints=.TRUE.,
 readBinaryPrec=32,
 writeBinaryPrec=32,
 debugLevel=1,
 /


# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=300,
#cg2dTargetResWunit=1.E-12,
 /


# Time stepping parameters
 &PARM03
 nIter0=0,
#2 lev2 for testing:
#nTimeSteps=8,
 nTimeSteps=350000,
#
 forcing_In_AB=.FALSE.,
 momDissip_In_AB=.FALSE.,
#when using the cd scheme:
# epsAB_CD = 0.25,
# tauCD=172800.0,
 deltaTmom   =3600.,
 deltaTtracer=3600.,
 deltaTfreesurf=3600.,
 deltaTClock =3600.,
#when using ab2:
# abEps = 0.1,
#when using ab3:
 doAB_onGtGs=.FALSE.,
 alph_AB=0.5,
 beta_AB=0.281105,
#
 pChkptFreq  =10368000.0,
 chkptFreq   =10368000.0,
# taveFreq    =2635200.0,
# dumpFreq    =2635200.0,
# monitorFreq = 7200.0,
 monitorFreq = 2592000.0,
 dumpInitAndLast = .TRUE.,
 adjDumpFreq = 31536000.0,
#adjDumpFreq = 604800.0,
 adjMonitorFreq = 864000.0,
 pickupStrictlyMatch=.FALSE.,
 /


# Gridding parameters
 &PARM04
 usingCurvilinearGrid=.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,
 /


# Input datasets
 &PARM05
 diffKrFile='total_diffkr_r009bit11.bin',
 adTapeDir='tapes',
#bathyFile      ='bathy_eccollc_90x50.bin',
 bathyFile      ='bathy_eccollc_90x50_min2pts.bin',
 hydrogThetaFile='T_OWPv1_M_eccollc_90x50.bin',
 hydrogSaltFile ='S_OWPv1_M_eccollc_90x50.bin',
 viscA4Dfile    ='fenty_biharmonic_visc_v11.bin',
 viscA4Zfile    ='fenty_biharmonic_visc_v11.bin',
 geothermalFile='geothermalFlux.bin',
#
 /
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20191220/220361fb/attachment.html>


More information about the MITgcm-support mailing list