# ==================== # | Model parameters | # ==================== # # Continuous equation parameters # # tRef - Reference vertical potential temperature (deg C) # sRef - Reference vertical salinity (PSU) # viscAh - Horizontal eddy viscosity coefficient (m^2/s) # viscAz - Vertical eddy viscosity coefficient (m^2/s) # diffKhT - Horizontal temperature diffusivity (m^2/s) # diffKzT - Vertical temperature diffusivity (m^2/s) # diffKhS - Horizontal salt diffusivity (m^2/s) # diffKzS - Vertical salt diffusivity (m^2/s) # f0 - Reference coriolis parameter, # south edge of f on beta plane (1/s) # beta - df/dy (s^-1.m^-1) # tAlpha - Linear EOS thermal expansion coefficient (1/oC) # sBeta - Linear EOS haline contraction coefficient (1/ppt) # gravity - Acceleration due to gravity (m/s^2) # gBaro - Accel. due to gravity used in barotropic equation (m/s^2) # rigidLid - Set to true to use rigid lid # implicitFreeSurface - Set to true to use implicit free surface # eosType - Flag for linear or polynomial equation of state # momAdvection - On/Off flag for momentum self transport # momViscosity - On/Off flag for momentum mixing # &PARM01 tRef= 24.0 , sRef= 34.65, no_slip_sides=.FALSE., no_slip_bottom=.TRUE., # -=-=- Orig Values -=-=- # viscAz=1.93e-5, # viscAh=5.E4, # diffKhT=0.0, # diffKzT=1.46e-5, # diffKhS=0.0, # diffKzS=1.46e-5, # -=-=- IF 150305 -=-=- viscAz=1.E-3, viscAh=1.E4, diffKhT=1.E2, diffKzT=1.E-5, diffKhS=1.E2, diffKzS=1.E-5, beta=1.E-11, tAlpha=2.E-4, sBeta =7.4E-4, rigidLid=.FALSE., implicitFreeSurface=.TRUE., eosType='JMD95Z', readBinaryPrec=32, writeBinaryPrec=32, saltStepping=.TRUE., tempStepping=.TRUE., momStepping=.TRUE., implicitDiffusion=.TRUE., implicitViscosity=.TRUE., globalFiles=.TRUE., allowFreezing=.FALSE., HeatCapacity_Cp = 3986.D0 gravity = 9.8156, gBaro = 9.8156, rhoNil = 1027.D0 useCDscheme=.FALSE., useNHMTerms=.TRUE., & # Elliptic solver parameters # # cg2dMaxIters - Maximum number of 2d solver iterations # cg2dTargetResidual - Solver target residual # &PARM02 cg2dMaxIters=500, cg2dTargetResidual=1.D-8 & # Time stepping parameters # # startTime - Integration starting time (s) # endTime - Integration ending time (s) # tauCD - CD scheme coupling timescale (s) # deltaTMom - Timestep for momemtum equations (s) # deltaTtracer - Tracer timestep (s) # deltaTClock - Timestep used as model "clock" (s) # abEps - Adams-Bashforth stabilising factor # pChkPtFreq - Frequency of permanent check pointing (s) # chkPtFreq - Frequency of rolling check pointing (s) # dumpFreq - Frequency at which model state is stored (s) # tauThetaClimRelax - Relaxation to climatology time scale (s) # tauSaltClimRelax - Relaxation to climatology time scale (s) # &PARM03 # tauCD=172800., # ============================================= # S I M U L A T I O N T I M E + D E L T A T # ============================================= niter0=0, # one month # nTimeSteps= 732, # two months ... # nTimeSteps=1464, # nTimeSteps = 168, nTimeSteps = 72, # nTimeSteps = 720, # one year # nTimeSteps=8784, # two years # nTimeSteps=17568, # six years # nTimeSteps=52704, # nine years # nTimeSteps=79056, # ten years # nTimeSteps = 87840, # -=-=- IF 150305 -=-=- deltaT = 3600.0, cAdjFreq=0., abEps=0.1, # ============================================= # C H P K O U T P U T F R E Q U E N C Y # ============================================= # -=-=- IF 180305 -=-=- # checkpoint once per year pChkptFreq= 31622400.0, # checkpoint once per month... # pChkptFreq = 2635200., # NEVER # chkptFreq= 0., # ============================================= # T A V E + D U M P O U T P U T F R E Q # ============================================= # NEVER dumpFreq = 0., # EVERY TIME STEP # dumpFreq = 3600., # DAILY # dumpFreq = 86400., # ONCE A WEEK # dumpFreq = 604800., # ---T I M E---A V E R A G E--- # NEVER # tavefreq = 0., # DAILY tavefreq = 864000., # ONCE A WEEK # tavefreq = 604800., # EVERY 10 HOURS # tavefreq = 36000., # ABOUT ONCE A MONTH # tavefreq = 2635200., # annually # tavefreq = 31622400., # ============================================= # R E L A X A T I O N T I M E S # ============================================= # -=-=- Orig Value -=-=- # tauSaltClimRelax = 4142330., # -=-=- IF 150305 -=-=- tauSaltClimRelax = 0., tauThetaClimRelax = 0., # ============================================= # M O N I T O R O U T P U T F R E Q # ============================================= # original value (output constantly) # monitorFreq=72., # quarterly... # monitorFreq=7905600., # semi-annually monitorFreq=15811200., # ============================================= # outputTypesInclusive=.TRUE., forcing_In_AB = .FALSE. & # Gridding parameters # # usingSphericalPolarGrid - On/Off flag for spherical polar coordinates # usingCartesianGrid - On/Off flag for selecting cartesian coordinates # delX - Zonal grid spacing (degrees) # delY - Meridional grid spacing (degrees) # delZ - Vertical grid spacing (m) # phiMin - Southern boundary latitude (degrees) # &PARM04 usingCartesianGrid=.FALSE., usingSphericalPolarGrid=.TRUE., # -=-=- IF 150305 -=-=- delX=41*1.E0, delY=33*1.E0, delZ= 10., phimin=45.5, thetamin=279.5, rSphere = 6371.D3 & # Input datasets # # bathyFile - File containing bathymetry # hydrogThetaFile - File containing initial potential temperature data # hydrogSaltFile - File containing initial salinity data # zonalWindFile - File containing zonal wind data # meridWindFile - File containing meridional wind data # thetaClimFile - File containing theta climatology used for relaxation # saltClimFile - File containing salt climatology used for relaxation # &PARM05 # -=-=- Orig Values -=-=- # -=-=- IF 150305 -=-=- bathyFile = 'LABSEA_bathy_fl_closed_w.bin', hydrogThetaFile = 'LABSEA_LEVITUS_1x1_ptmp_NEW_corK', hydrogSaltFile = 'LABSEA_LEVITUS_1x1_salt_NEW_corK', &