# ==================== # | Model parameters | # ==================== # # Continuous equation parameters &PARM01 tRef = 18.89, 18.89, 18.89, 18.89, 18.89, 18.87, 18.85, 18.82, 18.80, 18.73, 18.65, 18.57, 18.40, 18.22, 18.00, 17.74, 17.44, 17.12, 16.76, 16.39, 15.98, 15.55, 15.08, 14.59, 14.07, 13.53, 12.99, 12.47, 11.97, 11.49, 11.02, 10.57, 10.12, 9.71, 9.27, 8.88, 8.46, 8.09, 7.71, 7.37, 7.03, 6.72, 6.42, 6.13, 5.86, 5.59, 5.34, 5.09, 4.87, 4.65, 4.45, 4.26, 4.08, 3.91, 3.75, 3.60, 3.47, 3.33, 3.20, 3.08, 2.96, 2.84, 2.73, 2.62, 2.51, 2.42, 2.32, 2.23, 2.14, 2.06, 1.98, 1.90, 1.81, 1.73, 1.65, 1.57, 1.49, 1.41, 1.33, 1.24, 1.15, 1.06, 0.98, 0.94, 0.91, 0.92, 0.98, 0.98, 0.98, 0.98, sRef = 34.84, 34.84, 34.84, 34.84, 34.84, 34.84, 34.85, 34.85, 34.85, 34.86, 34.87, 34.88, 34.89, 34.90, 34.92, 34.94, 34.96, 34.98, 35.00, 35.02, 35.04, 35.06, 35.07, 35.07, 35.07, 35.05, 35.03, 35.01, 34.98, 34.95, 34.92, 34.89, 34.85, 34.82, 34.79, 34.76, 34.73, 34.71, 34.68, 34.66, 34.64, 34.62, 34.61, 34.60, 34.59, 34.59, 34.58, 34.58, 34.59, 34.59, 34.60, 34.60, 34.61, 34.62, 34.63, 34.64, 34.65, 34.66, 34.67, 34.68, 34.69, 34.70, 34.71, 34.71, 34.72, 34.72, 34.73, 34.73, 34.74, 34.74, 34.74, 34.74, 34.75, 34.74, 34.74, 34.74, 34.74, 34.74, 34.74, 34.74, 34.73, 34.73, 34.73, 34.73, 34.73, 34.72, 34.72, 34.72, 34.72, 34.72, viscAr= 5.6614e-04, no_slip_sides = .TRUE., no_slip_bottom = .TRUE., diffKrT=5.44e-7, diffKrS=5.44e-7, rhonil=1027.5, rhoConstFresh=999.8, eosType='JMD95Z', hFacMin=0.3, hFacInf=0.1, hFacSup=5., implicitDiffusion=.TRUE., implicitViscosity=.TRUE., viscC4Leith=2.15, viscC4Leithd=2.15, viscA4GridMax=0.8, useAreaViscLength=.TRUE., highOrderVorticity = .TRUE., bottomDragQuadratic = 0.0021, tempAdvScheme=7, saltAdvScheme=7, StaggerTimeStep=.TRUE., multiDimAdvection=.TRUE., vectorInvariantMomentum=.TRUE., implicitFreeSurface=.TRUE., exactConserv=.TRUE., debuglevel=-1, convertFW2Salt=-1., useRealFreshWaterFlux=.TRUE., # useSingleCPUio=.TRUE., useSingleCPUio=.FALSE., # useSingleCpuInput=.TRUE., # globalFiles=.TRUE., readBinaryPrec=32, writeBinaryPrec=32, temp_EvPrRn=0., implicSurfPress = 0.6, implicDiv2DFlow = 0.6, / # Elliptic solver parameters &PARM02 cg2dMaxIters = 100, cg2dTargetResidual = 1.E-5, / # Time stepping parameters &PARM03 nIter0=0, endtime=63072000., deltaT = 25., abEps = 0.1, forcing_In_AB=.FALSE., momDissip_In_AB=.FALSE., chkptFreq = 43200.0, dumpfreq = 3600., monitorFreq = 360., dumpInitAndLast=.FALSE., / # Gridding parameters &PARM04 usingCurvilinearGrid=.TRUE., delR = 1.00, 1.14, 1.30, 1.49, 1.70, 1.93, 2.20, 2.50, 2.84, 3.21, 3.63, 4.10, 4.61, 5.18, 5.79, 6.47, 7.20, 7.98, 8.83, 9.73, 10.69, 11.70, 12.76, 13.87, 15.03, 16.22, 17.45, 18.70, 19.97, 21.27, 22.56, 23.87, 25.17, 26.46, 27.74, 29.00, 30.24, 31.45, 32.65, 33.82, 34.97, 36.09, 37.20, 38.29, 39.37, 40.45, 41.53, 42.62, 43.73, 44.87, 46.05, 47.28, 48.56, 49.93, 51.38, 52.93, 54.61, 56.42, 58.38, 60.53, 62.87, 65.43, 68.24, 71.33, 74.73, 78.47, 82.61, 87.17, 92.21, 97.79, 103.96, 110.79, 118.35, 126.73, 136.01, 146.30, 157.71, 170.35, 184.37, 199.89, 217.09, 236.13, 257.21, 280.50, 306.24, 334.64, 365.93, 400.38, 438.23, 479.74, / # Input datasets &PARM05 bathyFile ='bathy4320_g5_r4', hydrogThetaFile='THETA_llc2160_10sep2011_4320x56160x90_r4', hydrogSaltFile ='SALT_llc2160_10sep2011_4320x56160x90_r4', uVelInitFile ='UVEL_llc2160_10sep2011_4320x56160x90_r4', vVelInitFile ='VVEL_llc2160_10sep2011_4320x56160x90_r4', pSurfInitFile ='ETAN_llc2160_10sep2011_4320x56160_r4', /