# ==================== # | 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= 1.77e-5 no_slip_sides=.TRUE., no_slip_bottom=.TRUE., diffKhT=0., diffK4T=0., diffKrT=1.51e-5 diffKhS=0., diffK4S=0., diffKrS=1.51e-5 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., viscA4 = 0., viscA4Grid = 1e-3, viscA4Max = 2e10, viscAh = 0., bottomDragQuadratic = 0.002, tempAdvScheme=33 saltAdvScheme=33 StaggerTimeStep=.true. multiDimAdvection=.true. vectorInvariantMomentum=.TRUE., rigidLid=.FALSE., implicitFreeSurface=.TRUE., exactConserv=.TRUE., debuglevel=0 & # Elliptic solver parameters &PARM02 cg2dMaxIters=100, cg2dTargetResidual=1.E-5, & # Time stepping parameters &PARM03 niter0=216, endTime=31881600, abEps=0.1, deltaT=1200, cAdjFreq=0., # day=369, 04-Jan-1993 pChkptFreq=31881600, chkptFreq=2592000 tavefreq=259200 tauSaltClimRelax = 3842380., monitorFreq=86400 & # 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_18360x85_ver5e_filled.bin', hydrogThetaFile='LEVT01_JAN_18360x85x50.bin', hydrogSaltFile ='LEVS01_JAN_18360x85x50.bin', &