/* This input file provides the adjustable run-time parameters for version 6 of
the Modular Ocean Model (MOM6), a numerical ocean model developed at NOAA-GFDL.
Where appropriate, parameters use usually given in MKS units.
This particular file is for the example in SIS2.
This MOM_input file typically contains only the non-default values that are
needed to reproduce this example. A full list of parameters for this example
can be found in the corresponding MOM_parameter_doc.all file which is
generated by the model at run-time. */
TRIPOLAR_N = True ! [Boolean] default = False
! Use tripolar connectivity at the northern edge of the
! domain. With TRIPOLAR_N, NIGLOBAL must be even.
!SYMMETRIC_MEMORY_ = False ! [Boolean]
! If defined, the velocity point data domain includes
! every face of the thickness points. In other words,
! some arrays are larger than others, depending on where
! they are on the staggered grid. Also, the starting
! index of the velocity-point arrays is usually 0, not 1.
! This can only be set at compile time.
!STATIC_MEMORY_ = False ! [Boolean]
! If STATIC_MEMORY_ is defined, the principle variables
! will have sizes that are statically determined at
! compile time. Otherwise the sizes are not determined
! until run time. The STATIC option is substantially
! faster, but does not allow the PE count to be changed
! at run time. This can only be set at compile time.
NIHALO = 4 ! default = 2
! The number of halo points on each side in the
! x-direction. With STATIC_MEMORY_ this is set as NIHALO_
! in MOM_memory.h at compile time; without STATIC_MEMORY_
! the default is NIHALO_ in MOM_memory.h (if defined) or 2.
NJHALO = 4 ! default = 2
! The number of halo points on each side in the
! y-direction. With STATIC_MEMORY_ this is set as NJHALO_
! in MOM_memory.h at compile time; without STATIC_MEMORY_
! the default is NJHALO_ in MOM_memory.h (if defined) or 2.
NIGLOBAL = 120 !
! The total number of thickness grid points in the
! x-direction in the physical domain. With STATIC_MEMORY_
! this is set in MOM_memory.h at compile time.
NJGLOBAL = 84 !
! The total number of thickness grid points in the
! y-direction in the physical domain. With STATIC_MEMORY_
! this is set in MOM_memory.h at compile time.
!NIPROC = 10 !
! The number of processors in the x-direction. With
! STATIC_MEMORY_ this is set in MOM_memory.h at compile time.
!NJPROC = 6 !
! The number of processors in the x-direction. With
! STATIC_MEMORY_ this is set in MOM_memory.h at compile time.
!LAYOUT = 10, 6 !
! The processor layout that was acutally used.
!IO_LAYOUT = 1, 1 ! default = 0
! The processor layout to be used, or 0,0 to automatically
! set the io_layout to be the same as the layout.
! === module MOM_grid ===
! Parameters providing information about the vertical grid.
NK = 28 ! [nondim]
! The number of model layers.
! === module MOM_diag_mediator ===
NUM_DIAG_COORDS = 2 ! default = 1
! The number of diagnostic vertical coordinates to use.
! For each coordinate, an entry in DIAG_COORDS must be provided.
DIAG_COORDS = "z Z ZSTAR, rho RHO RHO" ! default = "z Z ZSTAR"
! A list of string tuples associating diag_table modules to
! a coordinate definition used for diagnostics. Each string
! is of the form "MODULE_SUFFIX PARAMETER_SUFFIX COORDINATE_NAME".
DIAG_MISVAL = 1.0E+20 ! [not defined] default = 1.0E+20
! Set the default missing value to use for diagnostics.
DIAG_AS_CHKSUM = False ! [Boolean] default = False
! Instead of writing diagnostics to the diag manager, write
! a textfile containing the checksum (bitcount) of the array.
AVAILABLE_DIAGS_FILE = "available_diags.000000" ! default = "available_diags.000000"
! A file into which to write a list of all available
! ocean diagnostics that can be included in a diag_table.
DIAG_COORD_DEF_Z = "WOA09" ! default = "WOA09"
! Determines how to specify the coordinate
! resolution. Valid options are:
! PARAM - use the vector-parameter DIAG_COORD_RES_Z
! UNIFORM[:N] - uniformly distributed
! FILE:string - read from a file. The string specifies
! the filename and variable name, separated
! by a comma or space, e.g. FILE:lev.nc,dz
! or FILE:lev.nc,interfaces=zw
! WOA09[:N] - the WOA09 vertical grid (approximately)
! FNC1:string - FNC1:dz_min,H_total,power,precision
! HYBRID:string - read from a file. The string specifies
! the filename and two variable names, separated
! by a comma or space, for sigma-2 and dz. e.g.
! HYBRID:vgrid.nc,sigma2,dz
DIAG_COORD_DEF_RHO = "RFNC1:35,999.5,1028,1028.5,8.,1038.,0.0078125" ! Definition of the density vertical coordinate
! === module MOM_verticalGrid ===
! Parameters providing information about the vertical grid.
! === module MOM ===
DIABATIC_FIRST = True ! [Boolean] default = False
! If true, apply diabatic and thermodynamic processes,
! including buoyancy forcing and mass gain or loss,
! before stepping the dynamics forward.
THICKNESSDIFFUSE = True ! [Boolean] default = False
! If true, interfaces or isopycnal surfaces are diffused,
! depending on the value of FULL_THICKNESSDIFFUSE.
THICKNESSDIFFUSE_FIRST = True ! [Boolean] default = False
! If true, do thickness diffusion before dynamics.
! This is only used if THICKNESSDIFFUSE is true.
MIXEDLAYER_RESTRAT = True ! [Boolean] default = False
! If true, a density-gradient dependent re-stratifying
! flow is imposed in the mixed layer.
! This is only used if BULKMIXEDLAYER is true.
DT = 7200.0 ! [s]
! The (baroclinic) dynamics time step. The time-step that
! is actually used will be an integer fraction of the
! forcing time-step (DT_FORCING in ocean-only mode or the
! coupling timestep in coupled mode.)
DT_THERM = 7200.0 ! [s] default = 3600.0
! The thermodynamic and tracer advection time step.
! Ideally DT_THERM should be an integer multiple of DT
! and less than the forcing or coupling time-step.
! By default DT_THERM is set to DT.
MIN_Z_DIAG_INTERVAL = 2.16E+04 ! [s] default = 0.0
! The minimum amount of time in seconds between
! calculations of depth-space diagnostics. Making this
! larger than DT_THERM reduces the performance penalty
! of regridding to depth online.
DTBT_RESET_PERIOD = 0.0 ! [s] default = 7200.0
! The period between recalculations of DTBT (if DTBT <= 0).
! If DTBT_RESET_PERIOD is negative, DTBT is set based
! only on information available at initialization. If
! dynamic, DTBT will be set at least every forcing time
! step, and if 0, every dynamics time step. The default is
! set by DT_THERM. This is only used if SPLIT is true.
FRAZIL = True ! [Boolean] default = False
! If true, water freezes if it gets too cold, and the
! the accumulated heat deficit is returned in the
! surface state. FRAZIL is only used if
! ENABLE_THERMODYNAMICS is true.
DO_GEOTHERMAL = True ! [Boolean] default = False
! If true, apply geothermal heating.
BOUND_SALINITY = True ! [Boolean] default = False
! If true, limit salinity to being positive. (The sea-ice
! model may ask for more salt than is available and
! drive the salinity negative otherwise.)
C_P = 3925.0 ! [J kg-1 K-1] default = 3991.86795711963
! The heat capacity of sea water, approximated as a
! constant. This is only used if ENABLE_THERMODYNAMICS is
! true. The default value is from the TEOS-10 definition
! of conservative temperature.
CHECK_BAD_SURFACE_VALS = True ! [Boolean] default = False
! If true, check the surface state for ridiculous values.
NKML = 2 ! [nondim] default = 2
! The number of sublayers within the mixed layer if
! BULKMIXEDLAYER is true.
NKBL = 2 ! [nondim] default = 2
! The number of layers that are used as variable density
! buffer layers if BULKMIXEDLAYER is true.
SAVE_INITIAL_CONDS = True ! [Boolean] default = False
! If true, write the initial conditions to a file given
! by IC_OUTPUT_FILE.
IC_OUTPUT_FILE = "MOM_IC" ! default = "MOM_IC"
! The file into which to write the initial conditions.
! === module MOM_tracer_registry ===
! === module MOM_tracer_flow_control ===
USE_IDEAL_AGE_TRACER = True ! [Boolean] default = False
! If true, use the ideal_age_example tracer package.
USE_OCMIP2_CFC = False ! [Boolean] default = False
! If true, use the MOM_OCMIP2_CFC tracer package.
! === module ideal_age_example ===
! === module MOM_OCMIP2_CFC ===
INPUTDIR = "INPUT" ! default = "."
! The directory in which input files are found.
COORD_CONFIG = "file" !
! This specifies how layers are to be defined:
! file - read coordinate information from the file
! specified by (COORD_FILE).
! linear - linear based on interfaces not layesrs.
! ts_ref - use reference temperature and salinity
! ts_range - use range of temperature and salinity
! (T_REF and S_REF) to determine surface density
! and GINT calculate internal densities.
! gprime - use reference density (RHO_0) for surface
! density and GINT calculate internal densities.
! ts_profile - use temperature and salinity profiles
! (read from COORD_FILE) to set layer densities.
! USER - call a user modified routine.
COORD_FILE = "Layer_coord28.nc" !
! The file from which the coordinate densities are read.
! === module MOM_grid_init ===
GRID_CONFIG = "mosaic" !
! A character string that determines the method for
! defining the horizontal grid. Current options are:
! mosaic - read the grid from a mosaic (supergrid)
! file set by GRID_FILE.
! cartesian - use a (flat) Cartesian grid.
! spherical - use a simple spherical grid.
! mercator - use a Mercator spherical grid.
GRID_FILE = "ocean_hgrid.nc" !
! Name of the file from which to read horizontal grid data.
TOPO_CONFIG = "file" !
! This specifies how bathymetry is specified:
! file - read bathymetric information from the file
! specified by (TOPO_FILE).
! flat - flat bottom set to MAXIMUM_DEPTH.
! bowl - an analytically specified bowl-shaped basin
! ranging between MAXIMUM_DEPTH and MINIMUM_DEPTH.
! spoon - a similar shape to 'bowl', but with an vertical
! wall at the southern face.
! halfpipe - a zonally uniform channel with a half-sine
! profile in the meridional direction.
! benchmark - use the benchmark test case topography.
! DOME - use a slope and channel configuration for the
! DOME sill-overflow test case.
! DOME2D - use a shelf and slope configuration for the
! DOME2D gravity current/overflow test case.
! seamount - Gaussian bump for spontaneous motion test case.
! Phillips - ACC-like idealized topography used in the Phillips config.
! USER - call a user modified routine.
MAXIMUM_DEPTH = 6000.0 ! [m]
! The maximum depth of the ocean.
MINIMUM_DEPTH = 30 ! [m] default = 0.0
! If MASKING_DEPTH is unspecified, then anything shallower than
! MINIMUM_DEPTH is assumed to be land and all fluxes are masked out.
! If MASKING_DEPTH is specified, then all depths shallower than
! MINIMUM_DEPTH but depper than MASKING_DEPTH are rounded to MINIMUM_DEPTH.
CHANNEL_CONFIG = "list" ! default = "none"
! A parameter that determines which set of channels are
! restricted to specific widths. Options are:
! none - All channels have the grid width.
! global_1deg - Sets 16 specific channels appropriate
! for a 1-degree model, as used in CM2G.
! list - Read the channel locations and widths from a
! text file, like MOM_channel_list in the MOM_SIS
! test case.
! file - Read open face widths everywhere from a
! NetCDF file on the model grid.
CHANNEL_LIST_360_LON_CHECK = True ! [Boolean] default = True
! If true, the channel configuration list works for any
! longitudes in the range of -360 to 360.
REMAPPING_SCHEME = "PPM_H4" ! default = "PLM"
! This sets the reconstruction scheme used
! for vertical remapping for all variables.
! It can be one of the following schemes:
! PCM (1st-order accurate)
! PLM (2nd-order accurate)
! PPM_H4 (3rd-order accurate)
! PPM_IH4 (3rd-order accurate)
! PQM_IH4IH3 (4th-order accurate)
! PQM_IH6IH5 (5th-order accurate)
INIT_LAYERS_FROM_Z_FILE = True ! [Boolean] default = False
! If true, intialize the layer thicknesses, temperatures,
! and salnities from a Z-space file on a latitude-
! longitude grid.
TEMP_SALT_Z_INIT_FILE = "WOA05_pottemp_salt.nc" ! default = "temp_salt_z.nc"
! The name of the z-space input file used to initialize
! the layer thicknesses, temperatures and salinities.
Z_INIT_FILE_PTEMP_VAR = "PTEMP" ! default = "ptemp"
! The name of the potential temperature variable in
! TEMP_SALT_Z_INIT_FILE.
Z_INIT_FILE_SALT_VAR = "SALT" ! default = "salt"
! The name of the salinity variable in
! TEMP_SALT_Z_INIT_FILE.
! === module MOM_MEKE ===
USE_MEKE = False ! [Boolean] default = False
! If true, turns on the MEKE scheme which calculates
! a sub-grid mesoscale eddy kinetic energy budget.
MEKE_GMCOEFF = 1.0 ! [nondim] default = -1.0
! The efficiency of the conversion of potential energy
! into MEKE by the thickness mixing parameterization.
! If MEKE_GMCOEFF is negative, this conversion is not
! used or calculated.
! === module MOM_lateral_mixing_coeffs ===
USE_VARIABLE_MIXING = True ! [Boolean] default = False
! If true, the variable mixing code will be called. This
! allows diagnostics to be created even if the scheme is
! not used. If KHTR_SLOPE_CFF>0 or KhTh_Slope_Cff>0,
! this is set to true regardless of what is in the
! parameter file.
RESOLN_SCALED_KH = False ! [Boolean] default = False
! If true, the Laplacian lateral viscosity is scaled away
! when the first baroclinic deformation radius is well
! resolved.
RESOLN_SCALED_KHTH = False ! [Boolean] default = False
! If true, the interface depth diffusivity is scaled away
! when the first baroclinic deformation radius is well
! resolved.
KHTH_SLOPE_CFF = 0.25 ! [nondim] default = 0.0
! The nondimensional coefficient in the Visbeck formula
! for the interface depth diffusivity
KHTR_SLOPE_CFF = 0.25 ! [nondim] default = 0.0
! The nondimensional coefficient in the Visbeck formula
! for the epipycnal tracer diffusivity
VARMIX_KTOP = 2 ! [nondim] default = 2
! The layer number at which to start vertical integration
! of S*N for purposes of finding the Eady growth rate.
VISBECK_L_SCALE = 3.0E+04 ! [m] default = 0.0
! The fixed length scale in the Visbeck formula.
ETA_TOLERANCE = 1.0E-06 ! [m] default = 3.15E-09
! The tolerance for the differences between the
! barotropic and baroclinic estimates of the sea surface
! height due to the fluxes through each face. The total
! tolerance for SSH is 4 times this value. The default
! is 0.5*NK*ANGSTROM, and this should not be set less x
! than about 10^-15*MAXIMUM_DEPTH.
VELOCITY_TOLERANCE = 1.0E-04 ! [m s-1] default = 3.0E+08
! The tolerance for barotropic velocity discrepancies
! between the barotropic solution and the sum of the
! layer thicknesses.
CORIOLIS_EN_DIS = True ! [Boolean] default = False
! If true, two estimates of the thickness fluxes are used
! to estimate the Coriolis term, and the one that
! dissipates energy relative to the other one is used.
BOUND_CORIOLIS = True ! [Boolean] default = False
! If true, the Coriolis terms at u-points are bounded by
! the four estimates of (f+rv)v from the four neighboring
! v-points, and similarly at v-points. This option is
! always effectively false with CORIOLIS_EN_DIS defined and
! CORIOLIS_SCHEME set to SADOURNY75_ENERGY.
! === module MOM_hor_visc ===
LAPLACIAN = True ! [Boolean] default = False
! If true, use a Laplacian horizontal viscosity.
KH_VEL_SCALE = 0.01 ! [m s-1] default = 0.0
! The velocity scale which is multiplied by the grid
! spacing to calculate the Laplacian viscosity.
! The final viscosity is the largest of this scaled
! viscosity, the Smagorinsky viscosity and KH.
AH_VEL_SCALE = 0.05 ! [m s-1] default = 0.0
! The velocity scale which is multiplied by the cube of
! the grid spacing to calculate the biharmonic viscosity.
! The final viscosity is the largest of this scaled
! viscosity, the Smagorinsky viscosity and AH.
SMAGORINSKY_AH = True ! [Boolean] default = False
! If true, use a biharmonic Smagorinsky nonlinear eddy
! viscosity.
SMAG_BI_CONST = 0.06 ! [nondim] default = 0.0
! The nondimensional biharmonic Smagorinsky constant,
! typically 0.015 - 0.06.
! === module MOM_vert_friction ===
CHANNEL_DRAG = True ! [Boolean] default = False
! If true, the bottom drag is exerted directly on each
! layer proportional to the fraction of the bottom it
! overlies.
PRANDTL_TURB = 1.0 ! [nondim] default = 1.0
! The turbulent Prandtl number applied to shear
! instability.
DYNAMIC_VISCOUS_ML = True ! [Boolean] default = False
! If true, use a bulk Richardson number criterion to
! determine the mixed layer thickness for viscosity.
U_TRUNC_FILE = "U_velocity_truncations" ! default = ""
! The absolute path to a file into which the accelerations
! leading to zonal velocity truncations are written.
! Undefine this for efficiency if this diagnostic is not
! needed.
V_TRUNC_FILE = "V_velocity_truncations" ! default = ""
! The absolute path to a file into which the accelerations
! leading to meridional velocity truncations are written.
! Undefine this for efficiency if this diagnostic is not
! needed.
KV = 1.0E-04 ! [m2 s-1]
! The background kinematic viscosity in the interior.
! The molecular value, ~1e-6 m2 s-1, may be used.
HBBL = 10.0 ! [m]
! The thickness of a bottom boundary layer with a
! viscosity of KVBBL if BOTTOMDRAGLAW is not defined, or
! the thickness over which near-bottom velocities are
! averaged for the drag law if BOTTOMDRAGLAW is defined
! but LINEAR_DRAG is not.
MAXVEL = 6.0 ! [m s-1] default = 3.0E+08
! The maximum velocity allowed before the velocity
! components are truncated.
! === module MOM_set_visc ===
USE_JACKSON_PARAM = True ! [Boolean] default = False
! If true, use the Jackson-Hallberg-Legg (JPO 2008)
! shear mixing parameterization.
ML_OMEGA_FRAC = 1.0 ! [nondim] default = 0.0
! When setting the decay scale for turbulence, use this
! fraction of the absolute rotation rate blended with the
! local value of f, as sqrt((1-of)*f^2 + of*4*omega^2).
DRAG_BG_VEL = 0.1 ! [m s-1] default = 0.0
! DRAG_BG_VEL is either the assumed bottom velocity (with
! LINEAR_DRAG) or an unresolved velocity that is
! combined with the resolved velocity to estimate the
! velocity magnitude. DRAG_BG_VEL is only used when
! BOTTOMDRAGLAW is defined.
BBL_USE_EOS = True ! [Boolean] default = False
! If true, use the equation of state in determining the
! properties of the bottom boundary layer. Otherwise use
! the layer target potential densities.
BBL_THICK_MIN = 0.1 ! [m] default = 0.0
! The minimum bottom boundary layer thickness that can be
! used with BOTTOMDRAGLAW. This might be
! Kv / (cdrag * drag_bg_vel) to give Kv as the minimum
! near-bottom viscosity.
! === module MOM_barotropic ===
BOUND_BT_CORRECTION = True ! [Boolean] default = False
! If true, the corrective pseudo mass-fluxes into the
! barotropic solver are limited to values that require
! less than 0.1*MAXVEL to be accommodated.
!BT x-halo = 0 !
! The barotropic x-halo size that is actually used.
!BT y-halo = 0 !
! The barotropic y-halo size that is actually used.
BT_PROJECT_VELOCITY = True ! [Boolean] default = False
! If true, step the barotropic velocity first and project
! out the velocity tendancy by 1+BEBT when calculating the
! transport. The default (false) is to use a predictor
! continuity step to find the pressure field, and then
! to do a corrector continuity step using a weighted
! average of the old and new velocities, with weights
! of (1-BEBT) and BEBT.
BT_THICK_SCHEME = "FROM_BT_CONT" ! default = "FROM_BT_CONT"
! A string describing the scheme that is used to set the
! open face areas used for barotropic transport and the
! relative weights of the accelerations. Valid values are:
! ARITHMETIC - arithmetic mean layer thicknesses
! HARMONIC - harmonic mean layer thicknesses
! HYBRID - use arithmetic means for layers above the
! shallowest bottom, the harmonic mean for layers
! below, and a weighted average for layers that
! straddle that depth
! FROM_BT_CONT - use the average thicknesses kept
! in the h_u and h_v fields of the BT_cont_type
BEBT = 0.2 ! [nondim] default = 0.1
! BEBT determines whether the barotropic time stepping
! uses the forward-backward time-stepping scheme or a
! backward Euler scheme. BEBT is valid in the range from
! 0 (for a forward-backward treatment of nonrotating
! gravity waves) to 1 (for a backward Euler treatment).
! In practice, BEBT must be greater than about 0.05.
DTBT = -0.95 ! [s or nondim] default = -0.98
! The barotropic time step, in s. DTBT is only used with
! the split explicit time stepping. To set the time step
! automatically based the maximum stable value use 0, or
! a negative value gives the fraction of the stable value.
! Setting DTBT to 0 is the same as setting it to -0.98.
! The value of DTBT that will actually be used is an
! integer fraction of DT, rounding down.
! === module MOM_thickness_diffuse ===
KHTH = 600.0 ! [m2 s-1] default = 0.0
! The background horizontal thickness diffusivity.
KHTH_MAX = 900.0 ! [m2 s-1] default = 0.0
! The maximum horizontal thickness diffusivity.
! === module MOM_mixed_layer_restrat ===
FOX_KEMPER_ML_RESTRAT_COEF = 20.0 ! [nondim] default = 0.0
! A nondimensional coefficient that is proportional to
! the ratio of the deformation radius to the dominant
! lengthscale of the submesoscale mixed layer
! instabilities, times the minimum of the ratio of the
! mesoscale eddy kinetic energy to the large-scale
! geostrophic kinetic energy or 1 plus the square of the
! grid spacing over the deformation radius, as detailed
! by Fox-Kemper et al. (2010)
Z_OUTPUT_GRID_FILE = "analysis_vgrid_lev35.v1.nc" ! default = ""
! The file that specifies the vertical grid for
! depth-space diagnostics, or blank to disable
! depth-space output.
!NK_ZSPACE (from file) = 50 ! [nondim]
! The number of depth-space levels. This is determined
! from the size of the variable zw in the output grid file.
! === module MOM_diabatic_driver ===
! The following parameters are used for diabatic processes.
BBL_MIXING_AS_MAX = False ! [Boolean] default = True
! If true, take the maximum of the diffusivity from the
! BBL mixing and the other diffusivities. Otherwise,
! diffusiviy from the BBL_mixing is simply added.
USE_LOTW_BBL_DIFFUSIVITY = True ! [Boolean] default = False
! If true, uses a simple, imprecise but non-coordinate dependent, model
! of BBL mixing diffusivity based on Law of the Wall. Otherwise, uses
! the original BBL scheme.
KD_MAX = 0.1 ! [m2 s-1] default = -1.0
! The maximum permitted increment for the diapycnal
! diffusivity from TKE-based parameterizations, or a
! negative value for no limit.
HENYEY_IGW_BACKGROUND = True ! [Boolean] default = False
! If true, use a latitude-dependent scaling for the near
! surface background diffusivity, as described in
! Harrison & Hallberg, JPO 2008.
N2_FLOOR_IOMEGA2 = 0.0 ! [nondim] default = 1.0
! The floor applied to N2(k) scaled by Omega^2:
! If =0., N2(k) is simply positive definite.
! If =1., N2(k) > Omega^2 everywhere.
KD = 2.0E-05 ! [m2 s-1]
! The background diapycnal diffusivity of density in the
! interior. Zero or the molecular value, ~1e-7 m2 s-1,
! may be used.
ML_RADIATION = True ! [Boolean] default = False
! If true, allow a fraction of TKE available from wind
! work to penetrate below the base of the mixed layer
! with a vertical decay scale determined by the minimum
! of: (1) The depth of the mixed layer, (2) an Ekman
! length scale.
ML_RAD_COEFF = 0.1 ! [nondim] default = 0.2
! The coefficient which scales MSTAR*USTAR^3 to obtain
! the energy available for mixing below the base of the
! mixed layer. This is only used if ML_RADIATION is true.
MSTAR = 0.3 ! [units=nondim] default = 1.2
! The ratio of the friction velocity cubed to the TKE
! input to the mixed layer.
TKE_DECAY = 10.0 ! [nondim] default = 2.5
! The ratio of the natural Ekman depth to the TKE decay scale.
KD_MIN = 2.0E-06 ! [m2 s-1] default = 2.0E-07
! The minimum diapycnal diffusivity.
INT_TIDE_DISSIPATION = True ! [Boolean] default = False
! If true, use an internal tidal dissipation scheme to
! drive diapycnal mixing, along the lines of St. Laurent
! et al. (2002) and Simmons et al. (2004).
DISSIPATION_N0 = 0 ! [W m-3] default = 0.0
! The intercept when N=0 of the N-dependent expression
! used to set a minimum dissipation by which to determine
! a lower bound of Kd (a floor): A in eps_min = A + B*N.
DISSIPATION_N1 = 0 ! [J m-3] default = 0.0
! The coefficient multiplying N, following Gargett, used to
! set a minimum dissipation by which to determine a lower
! bound of Kd (a floor): B in eps_min = A + B*N
INT_TIDE_DECAY_SCALE = 300.3003003003003 ! [m] default = 0.0
! The decay scale away from the bottom for tidal TKE with
! the new coding when INT_TIDE_DISSIPATION is used.
KAPPA_ITIDES = 6.28319E-04 ! [m-1] default = 6.283185307179586E-04
! A topographic wavenumber used with INT_TIDE_DISSIPATION.
! The default is 2pi/10 km, as in St.Laurent et al. 2002.
KAPPA_H2_FACTOR = 0.75 ! [nondim] default = 1.0
! A scaling factor for the roughness amplitude with nINT_TIDE_DISSIPATION.
TKE_ITIDE_MAX = 0.1 ! [W m-2] default = 1000.0
! The maximum internal tide energy source availble to mix
! above the bottom boundary layer with INT_TIDE_DISSIPATION.
READ_TIDEAMP = True ! [Boolean] default = False
! If true, read a file (given by TIDEAMP_FILE) containing
! the tidal amplitude with INT_TIDE_DISSIPATION.
H2_FILE = "sgs_h2_MOM6.nc" !
! The path to the file containing the sub-grid-scale
! topographic roughness amplitude with INT_TIDE_DISSIPATION.
! === module MOM_KPP ===
! This is the MOM wrapper to CVmix:KPP
! See http://code.google.com/p/cvmix/
KPP%
%KPP
! === module MOM_diffConvection ===
! This module implements enhanced diffusivity as a
! function of static stability, N^2.
CONVECTION%
%CONVECTION
! === module MOM_entrain_diffusive ===
MAX_ENT_IT = 20 ! default = 5
! The maximum number of iterations that may be used to
! calculate the interior diapycnal entrainment.
TOLERANCE_ENT = 1.0E-05 ! [m] default = 2.683281572999748E-05
! The tolerance with which to solve for entrainment values.
! === module MOM_geothermal ===
GEOTHERMAL_SCALE = 0.001 ! [W m-2 or various] default = 0.0
! The constant geothermal heat flux, a rescaling
! factor for the heat flux read from GEOTHERMAL_FILE, or
! 0 to disable the geothermal heating.
GEOTHERMAL_FILE = "geo_hflux_newgrid.nc" ! default = ""
! The file from which the geothermal heating is to be
! read, or blank to use a constant heating rate.
! === module MOM_kappa_shear ===
MAX_RINO_IT = 25 ! [nondim] default = 50
! The maximum number of iterations that may be used to
! estimate the Richardson number driven mixing.
! === module MOM_mixed_layer ===
BULK_RI_ML = 0.05 ! [nondim]
! The efficiency with which mean kinetic energy released
! by mechanically forced entrainment of the mixed layer
! is converted to turbulent kinetic energy.
! ABSORB_ALL_SW = True ! [Boolean] default = False
! If true, all shortwave radiation is absorbed by the
! ocean, instead of passing through to the bottom mud.
HMIX_MIN = 2.0 ! [m] default = 0.0
! The minimum mixed layer depth if the mixed layer depth
! is determined dynamically.
! LIMIT_BUFFER_DETRAIN = True ! [Boolean] default = False
! If true, limit the detrainment from the buffer layers
! to not be too different from the neighbors.
DEPTH_LIMIT_FLUXES = 0.1 ! [m] default = 0.2
! The surface fluxes are scaled away when the total ocean
! depth is less than DEPTH_LIMIT_FLUXES.
ML_RESORT = True ! [Boolean] default = False
! If true, resort the topmost layers by potential density
! before the mixed layer calculations.
ML_PRESORT_NK_CONV_ADJ = 4 ! [nondim] default = 0
! Convectively mix the first ML_PRESORT_NK_CONV_ADJ
! layers before sorting when ML_RESORT is true.
CORRECT_ABSORPTION_DEPTH = True ! [Boolean] default = False
! If true, the depth at which penetrating shortwave
! radiation is absorbed is corrected by moving some of
! the heating upward in the water column.
DO_RIVERMIX = True ! [Boolean] default = False
! If true, apply additional mixing whereever there is
! runoff, so that it is mixed down to RIVERMIX_DEPTH,
! if the ocean is that deep.
RIVERMIX_DEPTH = 40.0 ! [m] default = 0.0
! The depth to which rivers are mixed if DO_RIVERMIX is
! defined.
! === module MOM_regularize_layers ===
REGULARIZE_SURFACE_LAYERS = True ! [Boolean] default = False
! If defined, vertically restructure the near-surface
! layers when they have too much lateral variations to
! allow for sensible lateral barotropic transports.
! === module MOM_opacity ===
VAR_PEN_SW = False ! [Boolean] default = False
! If true, use one of the CHL_A schemes specified by
! OPACITY_SCHEME to determine the e-folding depth of
! incoming short wave radiation.
CHL_FILE = "seawifs_1998-2006_GOLD_smoothed_2X.nc" !
! CHL_FILE is the file containing chl_a concentrations in
! the variable CHL_A. It is used when VAR_PEN_SW and
! CHL_FROM_FILE are true.
PEN_SW_NBANDS = 1 ! default = 1
! The number of bands of penetrating shortwave radiation.
! === module MOM_tracer_advect ===
TRACER_ADVECTION_SCHEME = "PPM:H3" ! default = "PLM"
! The horizontal transport scheme for tracers:
! PLM - Piecewise Linear Method
! PPM:H3 - Piecewise Parabolic Method (Huyhn 3rd order)
! === module MOM_tracer_hor_diff ===
KHTR = 600.0 ! [m2 s-1] default = 0.0
! The background along-isopycnal tracer diffusivity.
KHTR_MIN = 50.0 ! [m2 s-1] default = 0.0
! The minimum along-isopycnal tracer diffusivity.
KHTR_MAX = 900.0 ! [m2 s-1] default = 0.0
! The maximum along-isopycnal tracer diffusivity.
KHTR_PASSIVITY_COEFF = 0.0 ! [nondim] default = 0.0
! The coefficient that scales deformation radius over
! grid-spacing in passivity, where passiviity is the ratio
! between along isopycnal mxiing of tracers to thickness mixing.
! A non-zero value enables this parameterization.
DIFFUSE_ML_TO_INTERIOR = False ! [Boolean] default = False
! If true, enable epipycnal mixing between the surface
! boundary layer and the interior.
ML_KHTR_SCALE = 0.0 ! [nondim] default = 1.0
! With Diffuse_ML_interior, the ratio of the truly
! horizontal diffusivity in the mixed layer to the
! epipycnal diffusivity. The valid range is 0 to 1.
! === module ocean_model_init ===
ENERGYSAVEDAYS = 1 ! [days] default = 1.0
! The interval in units of TIMEUNIT between saves of the
! energies of the run and other globally summed diagnostics.
RESTORE_SALINITY = True ! [Boolean] default = False
! If true, the coupled driver will add a globally-balanced
! fresh-water flux that drives sea-surface salinity
! toward specified values.
OCEAN_SURFACE_STAGGER = "B" ! default = "B"
! A case-insensitive character string to indicate the
! staggering of the surface velocity field that is
! returned to the coupler. Valid values include
! 'A', 'B', or 'C'.
! === module MOM_surface_forcing ===
MAX_P_SURF = 7.0E+04 ! [Pa] default = -1.0
! The maximum surface pressure that can be exerted by the
! atmosphere and floating sea-ice or ice shelves. This is
! needed because the FMS coupling structure does not
! limit the water that can be frozen out of the ocean and
! the ice-ocean heat fluxes are treated explicitly. No
! limit is applied if a negative value is used.
ADJUST_NET_FRESH_WATER_TO_ZERO = True ! [Boolean] default = False
! If true, adjusts the net fresh-water forcing seen
! by the ocean (including restoring) to zero.
USE_LIMITED_PATM_SSH = False ! [Boolean] default = True
! If true, return the the sea surface height with the
! correction for the atmospheric (and sea-ice) pressure
! limited by max_p_surf instead of the full atmospheric
! pressure.
WIND_STAGGER = "B" ! default = "C"
! A case-insensitive character string to indicate the
! staggering of the input wind stress field. Valid
! values are 'A', 'B', or 'C'.
FLUXCONST = 0.5 ! [m day-1]
! The constant that relates the restoring surface fluxes
! to the relative surface anomalies (akin to a piston
! velocity). Note the non-MKS units.
CD_TIDES = 0.0025 ! [nondim] default = 1.0E-04
! The drag coefficient that applies to the tides.
READ_GUST_2D = False ! [Boolean] default = False
! If true, use a 2-dimensional gustiness supplied from
! an input file
GUST_2D_FILE = "gustiness_qscat.nc" !
! The file in which the wind gustiness is found in
! variable gustiness.
! === module MOM_sum_output ===
MAXTRUNC = 5000 ! [truncations save_interval-1] default = 0
! The run will be stopped, and the day set to a very
! large value if the velocity is truncated more than
! MAXTRUNC times between energy saves. Set MAXTRUNC to 0
! to stop if there is any truncation of velocities.