# ==================== # | Model parameters | # ==================== # to get a definition of the parameters # see /model/src/ini_parms.F, /model/src/set_defaults.F # see also model/src/modeldata_example.F &PARM01 # # Continuous equation parameters # # tRef - Reference vertical potential temperature (deg C) # sRef - Reference vertical salinity (PSU) # tRef and sRef are overridden if an initialization is read in later # 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 # eosType poly3 requires file: POLY3.COEFFS # GMkbackground - background value of GM/Redi coefficient # momAdvection - On/Off flag for momentum self transport # momViscosity - On/Off flag for momentum mixing # ##### tRef = 48*20., ##### sRef = 50*33., no_slip_sides=.FALSE., no_slip_bottom=.TRUE., viscAh = 2.E-2, viscAhGrid = 2.E-2, viscAr = 1.E-5, #- This is the standard Leith Viscosity, prop to grad(curl v_h) viscC4leith = 1.85, # -This is an additional Leith Viscosity, prop to grad(div.v_h) # It traps instabilities that have no signal in grad(curl v_h) viscC4leithD = 1.85, # This is harmonic Smagorinsky Coefficient # 0.2-0.9 = value from 3-d turbulence theory; Smagorinsky (1993) # 2.2-4 = stable oceanic values (Griffies and Hallberg, 2000) # 1. smag is approximately the same as 1.85 leith+leithd for this grid spacing viscC4smag = 1., #- This is the max Viscosity used, be it Smag or Leith, and it scales # with grid size and timestep #viscAhGridMax = 1., viscA4GridMax = 1., #- This uses the full grad(vort) and grad(div) useFullLeith=.TRUE., # Vertical diffusion diffKrT=1.E-6, diffKrS=1.E-6, # Biharmonic diffusion of Heat/Salt laterally diffK4T=0.1E7, diffK4S=0.1E7, # Bryan and Lewis 1979 diffusion # diffKrBL79surf=0.1E-4, # diffKrBL79deep=1.0E-4, implicitDiffusion=.TRUE., implicitViscosity=.TRUE., implicitFreeSurface=.TRUE., tempAdvScheme=30, saltAdvScheme=30, nonlinFreeSurf=0, useRealFreshWaterFlux=.FALSE., convertFW2Salt=-1, linFSConserveTr=.true. exactConserv=.FALSE., staggerTimeStep=.TRUE., vectorInvariantMomentum=.TRUE., useEnergyConservingCoriolis=.true., bottomDragQuadratic= 0.002, gravity=9.81, rhonil=1029., rhoConstFresh=1000., eosType='JMD95Z', useSingleCpuIO=.TRUE., allowFreezing=.false., hFacMin=.3, hFacMinDr=30., readBinaryPrec=32, debugLevel=1, / &PARM02 # # Elliptic solver parameters # # cg2dMaxIters - Maximum number of 2d solver iterations # cg2dTargetResidual - Solver target residual # cg2dMaxIters=1000, #FY cg2dMaxIters=5000, cg2dTargetResWunit=1.E-9, #FY cg2dTargetResWunit=1.E-11, # Elliptic solver parameters # cg2dMaxIters=1000, # cg2dTargetResidual=1.E-11, # cg3dMaxIters=40, # cg3dTargetResidual=1.E-9, / &PARM03 # Time stepping parameters # # nIter0 - Integration starting step (s) # nTimeSteps - Integration time in steps (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) # cAdjFreq - Convective adjustment period (s) # # nIter0=10483200, # nTimeSteps=350400, # nIter0=7714560, # nTimeSteps=350400, nIter0=0, nTimeSteps=9000000, deltaT =45., deltaTmom = 45., deltaTtracer = 45., deltaTClock = 45., pickupStrictlyMatch=.FALSE., cAdjFreq=0., abEps = 0.01, forcing_In_AB=.FALSE., pChkptFreq=7862400, chkptFreq=7862400, taveFreq=2628000, dumpFreq=86400, monitorFreq=2592000., # 2 months restoring timescale for temperature # tauThetaClimRelax= 5184000., # 3 months restoring timescale for salinity # climsssTauRelax=7776000., # tauSaltClimRelax = 2592000., # periodicExternalForcing=.TRUE., externForcingPeriod=10800., # externForcingCycle=31104000., # tauThetaClimRelax= 2592000., # tauSaltClimRelax= 7776000., # tauThetaClimRelax= 2592000., # tauSaltClimRelax= 2592000., / &PARM04 # 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) # usingCartesianGrid = .FALSE., usingSphericalPolarGrid= .TRUE., delR= 1.02,1.07, 1.137, 1.21, 1.31, 1.42,1.56, 1.73, 1.94, 2.19, 2.50, 2.88, 3.34, 3.90, 4.59, 5.42, 6.43,7.67, 9.16, 10.98, 13.18, 15.85, 19.06, 22.93, 27.58, 33.14, 39.75, 47.59, 56.83, 67.62, 80.13, 94.49, 110.79, 129.07, 149.26, 171.22, 194.71, 219.37, 244.78, 270.44, 295.85, 320.51, 343.99, 365.96, 386.15, 404.42, 420.73, 435.09, xgOrigin=105.5771, ygOrigin=15.0979, delX=500*0.03, delY=404*0.03, / &PARM05 # Input datasets # # bathyFile - File containing bathymetry # hydrogThetaFile - File containing initial potential temperature data # hydrogSaltFile - File containing initial potential salinity data # uVelInitFile - File containing initial zonal velocity data # vVelInitFile - File containing initial meridional velocity data # bathyFile='/mnt/disk_B/user/MITgcm/verification/case_scs/input/InputFiles/mybathy.bin', hydrogThetaFile='/mnt/disk_B/user/MITgcm/verification/case_scs/input/InputFiles/T0.bin', hydrogSaltFile='/mnt/disk_B/user/MITgcm/verification/case_scs/input/InputFiles/S0.bin', uVelInitFile='/mnt/disk_B/user/MITgcm/verification/case_scs/input/InputFiles/U0.bin', vVelInitFile='/mnt/disk_B/user/MITgcm/verification/case_scs/input/InputFiles/V0.bin', /