# ==================== # | 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=.FALSE., no_slip_bottom=.TRUE., viscAh = 1.E2, viscAhGrid = 0.01, viscAr =1.E-5, # diffKhT=0.2E-1, diffKhS=0.2E-1, diffK4T=0., diffK4S=0., diffKrT=1.E-5, diffKrS=1.E-5, diffKrBL79surf=0.1E-4, diffKrBL79deep=1.0E-4, bottomDragQuadratic = 0.002, ### ivdc_kappa=10., implicitDiffusion=.TRUE., implicitViscosity=.TRUE., gravity=9.81, rhonil=1029., rhoConstFresh=1000., eosType='JMD95Z', staggerTimeStep=.TRUE., vectorInvariantMomentum=.TRUE., implicitFreeSurface=.TRUE., exactConserv=.TRUE., useSingleCpuIO=.TRUE., #ph( ### select_rStar=2, hFacInf=0.2, hFacSup=2.0, nonlinFreeSurf=1, useRealFreshWaterFlux=.FALSE., #ph) allowFreezing=.TRUE., hFacMin=.3, hFacMinDr=50., readBinaryPrec=64, #writeBinaryPrec=64, / # Elliptic solver parameters &PARM02 cg2dMaxIters=200, cg2dTargetResWunit=1.E-10, / # Time stepping parameters &PARM03 nIter0=0, # 1 year at 600sec ### nTimeSteps=52560, # 1 year at 1800sec nTimeSteps=17520, # deltaT = 1800., #deltaTmom =1200., #deltaTtracer=86400., #deltaTfreesurf=86400., #deltaTClock =86400., abEps = 0.1, forcing_In_AB=.FALSE., pChkptFreq =2635200.0, chkptFreq =2635200.0, # taveFreq =2635200.0, dumpFreq =2635200.0, monitorFreq =864000., #latBandClimRelax=50., ### tauThetaClimRelax= 2592000., ### tauSaltClimRelax = 2592000., / # Gridding parameters &PARM04 usingCurvilinearGrid=.TRUE., ###horizGridFile='grid_cs32', 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 ='llc90x90x360_bathy_mindepth_noisolated.bin', ### hydrogThetaFile='llc90x90x360_50lev_ptemp.bin', ### hydrogSaltFile ='llc90x90x360_50lev_salt.bin', hydrogThetaFile='DiagTheta_diffus.bin', hydrogSaltFile ='DiagSalt_diffus.bin', /