# ==================== # | 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, viscAr= 6.0826e-04, no_slip_sides=.FALSE., no_slip_bottom=.TRUE., diffKhT=0., diffK4T=0., diffKrT=0., diffKhS=0., diffK4S=0., diffKrS=0., tAlpha=2.E-4, sBeta =7.4E-4, gravity=9.81, eosType='JMD95Z', readBinaryPrec=32, writeBinaryPrec=32, useSingleCPUio=.TRUE., hFacMinDr=50., hFacMin=0.3, hFacInf=0.1, hFacSup=5., # z* - begin select_rStar=2, nonlinFreeSurf=4, # * - end implicitDiffusion=.true., implicitViscosity=.true., allowFreezing=.false., viscAh=0., viscAhGrid=0., # Following change took place starting in checkpoint57v_post # Redefine viscC?Leith*. Change the coefficients of Leith viscosity # to be proportional to enstrophy dissipation length. Thus, # viscC?leith?_new=pi*(viscC?leith?_old)^(1/3) # This is now consistent with Smagorinsky and with the manual. # viscC4Leith=1.0, # viscC4Leithd=1.0, # viscC4Leith=3.1416, # viscC4Leithd=3.1416, viscC4Leith=1.5, viscC4Leithd=1.5, viscA4GridMax=0.5, ########################################################## viscC2Leith=0., viscA4=0., useAreaViscLength=.TRUE., sideDragFactor=0., SadournyCoriolis=.FALSE., useAbsVorticity=.FALSE., highOrderVorticity = .TRUE., bottomDragQuadratic = 0.002, tempAdvScheme=33, saltAdvScheme=33, StaggerTimeStep=.true., multiDimAdvection=.true., vectorInvariantMomentum=.TRUE., rigidLid=.FALSE., implicitFreeSurface=.TRUE., exactConserv=.TRUE., debuglevel=-1, # useRealFreshWaterFlux=.TRUE., & # Elliptic solver parameters &PARM02 cg2dMaxIters=100, cg2dTargetResidual=1.E-5, & # Time stepping parameters &PARM03 niter0=216, # endtime for December 31, 2002 endtime=347155200., abEps=0.1, deltaT=1200, cAdjFreq=0., pChkptFreq=31536000., chkptFreq=2592000., # blending: 12/24/06 tauSaltClimRelax = 8722203., # tauSaltClimRelax = 3842380., monitorFreq=86400., forcing_In_AB = .FALSE., dumpInitAndLast=.FALSE., & # Gridding parameters &PARM04 usingCartesianGrid=.FALSE., usingSphericalPolarGrid=.FALSE., 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 bathyFile = 'ETOPO2_510x6x510_ver5e_filled.bin', hydrogThetaFile= 'WGHC_T_510x6x510.bin', hydrogSaltFile = 'WGHC_S_510x6x510.bin', diffKrFile = 'DIFFKR_2_20_1', &