MOM_parameter_doc.all

! This file was written by the model and records all non-layout or debugging parameters used at run-time.

! === module MOM ===
SPLIT = True                    !   [Boolean] default = True
                                ! Use the split time stepping if true.
SPLIT_RK2B = False              !   [Boolean] default = False
                                ! If true, use a version of the split explicit time stepping scheme that
                                ! exchanges velocities with step_MOM that have the average barotropic phase over
                                ! a baroclinic timestep rather than the instantaneous barotropic phase.
CALC_RHO_FOR_SEA_LEVEL = False  !   [Boolean] default = False
                                ! If true, the in-situ density is used to calculate the effective sea level that
                                ! is returned to the coupler. If false, the Boussinesq parameter RHO_0 is used.
ENABLE_THERMODYNAMICS = False   !   [Boolean] default = True
                                ! If true, Temperature and salinity are used as state variables.
USE_EOS = False                 !   [Boolean] default = False
                                ! If true,  density is calculated from temperature and salinity with an equation
                                ! of state.  If USE_EOS is true, ENABLE_THERMODYNAMICS must be true as well.
DIABATIC_FIRST = False          !   [Boolean] default = False
                                ! If true, apply diabatic and thermodynamic processes, including buoyancy
                                ! forcing and mass gain or loss, before stepping the dynamics forward.
USE_CONTEMP_ABSSAL = False      !   [Boolean] default = False
                                ! If true, the prognostics T&S are the conservative temperature and absolute
                                ! salinity. Care should be taken to convert them to potential temperature and
                                ! practical salinity before exchanging them with the coupler and/or reporting
                                ! T&S diagnostics.
ADIABATIC = False               !   [Boolean] default = False
                                ! There are no diapycnal mass fluxes if ADIABATIC is true.  This assumes that KD
                                ! = 0.0 and that there is no buoyancy forcing, but makes the model faster by
                                ! eliminating subroutine calls.
DO_DYNAMICS = True              !   [Boolean] default = True
                                ! If False, skips the dynamics calls that update u & v, as well as the gravity
                                ! wave adjustment to h. This may be a fragile feature, but can be useful during
                                ! development
OFFLINE_TRACER_MODE = False     !   [Boolean] default = False
                                ! If true, barotropic and baroclinic dynamics, thermodynamics are all bypassed
                                ! with all the fields necessary to integrate the tracer advection and diffusion
                                ! equation are read in from files stored from a previous integration of the
                                ! prognostic model. NOTE: This option only used in the ocean_solo_driver.
USE_REGRIDDING = False          !   [Boolean] default = False
                                ! If True, use the ALE algorithm (regridding/remapping). If False, use the
                                ! layered isopycnal algorithm.
BULKMIXEDLAYER = False          !   [Boolean] default = False
                                ! If true, use a Kraus-Turner-like bulk mixed layer with transitional buffer
                                ! layers.  Layers 1 through NKML+NKBL have variable densities. There must be at
                                ! least NKML+NKBL+1 layers if BULKMIXEDLAYER is true. BULKMIXEDLAYER can not be
                                ! used with USE_REGRIDDING. The default is influenced by ENABLE_THERMODYNAMICS.
USE_POROUS_BARRIER = False      !   [Boolean] default = False
                                ! If true, use porous barrier to constrain the widths and face areas at the
                                ! edges of the grid cells.
BATHYMETRY_AT_VEL = False       !   [Boolean] default = False
                                ! If true, there are separate values for the basin depths at velocity points.
                                ! Otherwise the effects of topography are entirely determined from thickness
                                ! points.
ENABLE_BUGS_BY_DEFAULT = True   !   [Boolean] default = True
                                ! If true, the defaults for certain recently added bug-fix flags are set to
                                ! recreate the bugs so that the code can be moved forward without changing
                                ! answers for existing configurations.  The defaults for groups of bug-fix flags
                                ! are periodically changed to correct the bugs, at which point this parameter
                                ! will no longer be used to set their default.  Setting this to false means that
                                ! bugs are only used if they are actively selected, but it also means that
                                ! answers may change when code is updated due to newly found bugs.
DT = 120.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 = 120.0                !   [s] default = 120.0
                                ! The thermodynamic time step. Ideally DT_THERM should be an integer multiple of
                                ! DT and of DT_TRACER_ADVECT and less than the forcing or coupling time-step.
                                ! However, if THERMO_SPANS_COUPLING is true, DT_THERM can be an integer multiple
                                ! of the coupling timestep. By default DT_THERM is set to DT.
THERMO_SPANS_COUPLING = False   !   [Boolean] default = False
                                ! If true, the MOM will take thermodynamic timesteps that can be longer than the
                                ! coupling timestep. The actual thermodynamic timestep that is used in this case
                                ! is the largest integer multiple of the coupling timestep that is less than or
                                ! equal to DT_THERM.
DT_TRACER_ADVECT = 120.0        !   [s] default = 120.0
                                ! The tracer advection time step. Ideally DT_TRACER_ADVECT should be an integer
                                ! multiple of DT, less than DT_THERM, and less than the forcing or coupling
                                ! time-step. However, if TRADV_SPANS_COUPLING is true, DT_TRACER_ADVECT can be
                                ! longer than the coupling timestep. By default DT_TRACER_ADVECT is set to
                                ! DT_THERM.
TRADV_SPANS_COUPLING = False    !   [Boolean] default = False
                                ! If true, the MOM will take tracer advection timesteps that can be longer than
                                ! the coupling timestep. The actual tracer advection timestep that is used in
                                ! this case is the largest integer multiple of the coupling timestep that is
                                ! less than or equal to DT_TRACER_ADVECT.
THICKNESSDIFFUSE = False        !   [Boolean] default = False
                                ! If true, isopycnal surfaces are diffused with a Laplacian coefficient of KHTH.
APPLY_INTERFACE_FILTER = False  !   [Boolean] default = False
                                ! If true, model interface heights are subjected to a grid-scale dependent
                                ! spatial smoothing, often with biharmonic filter.
HMIX_SFC_PROP = 1.0             !   [m] default = 1.0
                                ! If BULKMIXEDLAYER is false, HMIX_SFC_PROP is the depth over which to average
                                ! to find surface properties like SST and SSS or density (but not surface
                                ! velocities).
HMIX_UV_SFC_PROP = 0.0          !   [m] default = 0.0
                                ! If BULKMIXEDLAYER is false, HMIX_UV_SFC_PROP is the depth over which to
                                ! average to find surface flow properties, SSU, SSV. A non-positive value
                                ! indicates no averaging.
HFREEZE = -1.0                  !   [m] default = -1.0
                                ! If HFREEZE > 0, melt potential will be computed. The actual depth over which
                                ! melt potential is computed will be min(HFREEZE, OBLD), where OBLD is the
                                ! boundary layer depth. If HFREEZE <= 0 (default), melt potential will not be
                                ! computed.
INTERPOLATE_P_SURF = False      !   [Boolean] default = False
                                ! If true, linearly interpolate the surface pressure over the coupling time
                                ! step, using the specified value at the end of the step.
DTBT_RESET_PERIOD = -1.0        !   [s] default = 120.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 0, DTBT will be set every dynamics time step. The default
                                ! is set by DT_THERM.  This is only used if SPLIT is true.
FIRST_DIRECTION = 0             ! default = 0
                                ! An integer that indicates which direction goes first in parts of the code that
                                ! use directionally split updates, with even numbers (or 0) used for x- first
                                ! and odd numbers used for y-first.
ALTERNATE_FIRST_DIRECTION = False !   [Boolean] default = False
                                ! If true, after every dynamic timestep alternate whether the x- or y- direction
                                ! updates occur first in directionally split parts of the calculation. If this
                                ! is true, FIRST_DIRECTION applies at the start of a new run or if the next
                                ! first direction can not be found in the restart file.
CHECK_BAD_SURFACE_VALS = False  !   [Boolean] default = False
                                ! If true, check the surface state for ridiculous values.
DEFAULT_ANSWER_DATE = 99991231  ! default = 99991231
                                ! This sets the default value for the various _ANSWER_DATE parameters.
SURFACE_ANSWER_DATE = 99991231  ! default = 99991231
                                ! The vintage of the expressions for the surface properties.  Values below
                                ! 20190101 recover the answers from the end of 2018, while higher values use
                                ! updated and more robust forms of the same expressions.
SAVE_INITIAL_CONDS = True       !   [Boolean] default = False
                                ! If true, write the initial conditions to a file given by IC_OUTPUT_FILE.
IC_OUTPUT_FILE = "GOLD_IC"      ! default = "MOM_IC"
                                ! The file into which to write the initial conditions.
WRITE_GEOM = 1                  ! default = 1
                                ! If =0, never write the geometry and vertical grid files. If =1, write the
                                ! geometry and vertical grid files only for a new simulation. If =2, always
                                ! write the geometry and vertical grid files. Other values are invalid.
GEOM_FILE = "ocean_geometry"    ! default = "ocean_geometry"
                                ! The file into which to write the ocean geometry.
USE_DBCLIENT = False            !   [Boolean] default = False
                                ! If true, initialize a client to a remote database that can be used for online
                                ! analysis and machine-learning inference.
ICE_SHELF = False               !   [Boolean] default = False
                                ! If true, enables the ice shelf model.
USE_PARTICLES = False           !   [Boolean] default = False
                                ! If true, use the particles package.
ENSEMBLE_OCEAN = False          !   [Boolean] default = False
                                ! If False, The model is being run in serial mode as a single realization. If
                                ! True, The current model realization is part of a larger ensemble and at the
                                ! end of step MOM, we will perform a gather of the ensemble members for
                                ! statistical evaluation and/or data assimilation.
HOMOGENIZE_FORCINGS = False     !   [Boolean] default = False
                                ! If True, homogenize the forces and fluxes.

! === module MOM_domains ===
REENTRANT_X = False             !   [Boolean] default = True
                                ! If true, the domain is zonally reentrant.
REENTRANT_Y = False             !   [Boolean] default = False
                                ! If true, the domain is meridionally reentrant.
TRIPOLAR_N = False              !   [Boolean] default = False
                                ! Use tripolar connectivity at the northern edge of the domain.  With
                                ! TRIPOLAR_N, NIGLOBAL must be even.
NIGLOBAL = 25                   !
                                ! 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 = 25                   !
                                ! 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.
NIHALO = 4                      ! default = 4
                                ! The number of halo points on each side in the x-direction.  How this is set
                                ! varies with the calling component and static or dynamic memory configuration.
NJHALO = 4                      ! default = 4
                                ! The number of halo points on each side in the y-direction.  How this is set
                                ! varies with the calling component and static or dynamic memory configuration.

! === module MOM_hor_index ===
! Sets the horizontal array index types.

! === module MOM_grid ===
! Parameters providing information about the lateral grid.
REFERENCE_HEIGHT = 0.0          !   [m] default = 0.0
                                ! A reference value for geometric height fields, such as bathyT.

! === module MOM_grid_init ===
GRID_CONFIG = "cartesian"       !
                                ! 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.
AXIS_UNITS = "k"                ! default = "degrees"
                                ! The units for the Cartesian axes. Valid entries are:
                                !     degrees - degrees of latitude and longitude
                                !     m or meter(s) - meters
                                !     k or km or kilometer(s) - kilometers
SOUTHLAT = 0.0                  !   [km]
                                ! The southern latitude of the domain or the equivalent starting value for the
                                ! y-axis.
LENLAT = 100.0                  !   [km]
                                ! The latitudinal or y-direction length of the domain.
WESTLON = 0.0                   !   [km] default = 0.0
                                ! The western longitude of the domain or the equivalent starting value for the
                                ! x-axis.
LENLON = 100.0                  !   [km]
                                ! The longitudinal or x-direction length of the domain.
RAD_EARTH = 6.378E+06           !   [m] default = 6.378E+06
                                ! The radius of the Earth.
READ_MEAN_SEA_LEVEL = False     !   [Boolean] default = False
                                ! If true, use a 2D map for time mean sea level, which is used to calculate time
                                ! mean ocean total thickness.
TOPO_CONFIG = "flat"            !
                                ! 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.
                                !     bbuilder - build topography from list of functions.
                                !     benchmark - use the benchmark test case topography.
                                !     Neverworld - use the Neverworld test case topography.
                                !     DOME - use a slope and channel configuration for the
                                !       DOME sill-overflow test case.
                                !     ISOMIP - use a slope and channel configuration for the
                                !       ISOMIP test case.
                                !     DOME2D - use a shelf and slope configuration for the
                                !       DOME2D gravity current/overflow test case.
                                !     Kelvin - flat but with rotated land mask.
                                !     seamount - Gaussian bump for spontaneous motion test case.
                                !     dumbbell - Sloshing channel with reservoirs on both ends.
                                !     shelfwave - exponential slope for shelfwave test case.
                                !     Phillips - ACC-like idealized topography used in the Phillips config.
                                !     dense - Denmark Strait-like dense water formation and overflow.
                                !     USER - call a user modified routine.
MINIMUM_DEPTH = 1.0             !   [m] default = 0.0
                                ! The minimum depth of the ocean.
MAXIMUM_DEPTH = 600.0           !   [m]
                                ! The maximum depth of the ocean.

! === module MOM_open_boundary ===
! Controls where open boundaries are located, what kind of boundary condition to impose, and what data to apply,
! if any.
OBC_NUMBER_OF_SEGMENTS = 4      ! default = 0
                                ! The number of open boundary segments.
OBC_VORTICITY_CONFIG = "freeslip" ! default = "freeslip"
                                ! Configuration for relative vorticity in momentum advection at open boundaries.
                                ! Options are:
                                !     none      - No adjustment.
                                !     zero      - Sets relative vorticity to zero.
                                !     freeslip  - Sets the normal gradient of tangential velocity to zero.
                                !     computed  - Computes the normal gradient of tangential velocity using
                                !                 external values of tangential velocity.
                                !     specified - Uses the external values of the normal gradient of
                                !                 tangential velocity.
OBC_STRAIN_CONFIG = "freeslip"  ! default = "freeslip"
                                ! Configuration for strain in horizontal viscosity at open boundaries.  Options
                                ! are:
                                !     none      - No adjustment.
                                !     zero      - Sets strain to zero.
                                !     freeslip  - Sets the normal gradient of tangential velocity to zero.
                                !     computed  - Computes the normal gradient of tangential velocity using
                                !                 external values of tangential velocity.
                                !     specified - Uses the external values of the normal gradient of
                                !                 tangential velocity.
OBC_ZERO_BIHARMONIC = True      !   [Boolean] default = False
                                ! If true, zeros the Laplacian of flow on open boundaries in the biharmonic
                                ! viscosity term.
MASK_OUTSIDE_OBCS = False       !   [Boolean] default = False
                                ! If true, set the areas outside open boundaries to be land.
RAMP_OBCS = False               !   [Boolean] default = False
                                ! If true, ramps from zero to the external values over time, with a ramping
                                ! timescale given by RAMP_TIMESCALE. Ramping SSH only so far.
OBC_RAMP_TIMESCALE = 1.0        !   [days] default = 1.0
                                ! If RAMP_OBCS is true, this sets the ramping timescale.
OBC_TIDE_N_CONSTITUENTS = 0     ! default = 0
                                ! Number of tidal constituents being added to the open boundary.
EXTERIOR_OBC_BUG = True         !   [Boolean] default = True
                                ! If true, recover a bug in barotropic solver and other routines when boundary
                                ! contitions interior to the domain are used.
OBC_HOR_INDEXING_BUG = True     !   [Boolean] default = True
                                ! If true, recover set of a horizontal indexing bugs in the OBC code.
OBC_TEMP_SALT_NEEDED_BUG = True !   [Boolean] default = True
                                ! If true, recover a bug that OBC temperature and salinity can be ignored even
                                ! if they are registered tracers in the rest of the model.
OBC_SEGMENT_001 = "J=N,I=N:0,FLATHER,ORLANSKI" !
                                ! Documentation needs to be dynamic?????
OBC_SEGMENT_002 = "J=0,I=0:N,FLATHER,ORLANSKI" !
                                ! Documentation needs to be dynamic?????
OBC_SEGMENT_003 = "I=N,J=0:N,FLATHER,ORLANSKI" !
                                ! Documentation needs to be dynamic?????
OBC_SEGMENT_004 = "I=0,J=N:0,FLATHER,ORLANSKI" !
                                ! Documentation needs to be dynamic?????
OBC_RADIATION_MAX = 10.0        !   [nondim] default = 1.0
                                ! The maximum magnitude of the baroclinic radiation velocity (or speed of
                                ! characteristics), in gridpoints per timestep.  This is only used if one of the
                                ! open boundary segments is using Orlanski.
OBC_RAD_VEL_WT = 0.3            !   [nondim] default = 0.3
                                ! The relative weighting for the baroclinic radiation velocities (or speed of
                                ! characteristics) at the new time level (1) or the running mean (0) for
                                ! velocities. Valid values range from 0 to 1. This is only used if one of the
                                ! open boundary segments is using Orlanski.
OBC_TRACER_RESERVOIR_LENGTH_SCALE_OUT = 0.0 !   [m] default = 0.0
                                ! An effective length scale for restoring the tracer concentration at the
                                ! boundaries to externally imposed values when the flow is exiting the domain.
OBC_TRACER_RESERVOIR_LENGTH_SCALE_IN = 0.0 !   [m] default = 0.0
                                ! An effective length scale for restoring the tracer concentration at the
                                ! boundaries to values from the interior when the flow is entering the domain.
OBC_THICKNESS_RESERVOIR_LENGTH_SCALE_OUT = 0.0 !   [m] default = 0.0
                                ! An effective length scale for restoring the layer thickness at the boundaries
                                ! to externally imposed values when the flow is exiting the domain.
OBC_THICKNESS_RESERVOIR_LENGTH_SCALE_IN = 0.0 !   [m] default = 0.0
                                ! An effective length scale for restoring the layer thickness at the boundaries
                                ! to values from the interior when the flow is entering the domain.
OBC_REMAPPING_SCHEME = "PLM"    ! default = "PLM"
                                ! This sets the reconstruction scheme used for OBC vertical remapping for all
                                ! variables. It can be one of the following schemes:
                                ! PCM         (1st-order accurate)
                                ! PLM         (2nd-order accurate)
                                ! PLM_HYBGEN  (2nd-order accurate)
                                ! PPM_H4      (3rd-order accurate)
                                ! PPM_IH4     (3rd-order accurate)
                                ! PPM_HYBGEN  (3rd-order accurate)
                                ! WENO_HYBGEN (3rd-order accurate)
                                ! PQM_IH4IH3  (4th-order accurate)
                                ! PQM_IH6IH5  (5th-order accurate)
BRUSHCUTTER_MODE = False        !   [Boolean] default = False
                                ! If true, read external OBC data on the supergrid.
REMAPPING_ANSWER_DATE = 99991231 ! default = 99991231
                                ! The vintage of the expressions and order of arithmetic to use for remapping.
                                ! Values below 20190101 result in the use of older, less accurate expressions
                                ! that were in use at the end of 2018.  Higher values result in the use of more
                                ! robust and accurate forms of mathematically equivalent expressions.
OBC_REMAPPING_USE_OM4_SUBCELLS = True !   [Boolean] default = True
                                ! If true, use the OM4 remapping-via-subcells algorithm for neutral diffusion.
                                ! See REMAPPING_USE_OM4_SUBCELLS for more details. We recommend setting this
                                ! option to false.
OBC_PROJECTION_BUG = True       !   [Boolean] default = True
                                ! If false, use only interior ocean points at OBCs to specify several
                                ! calculations at OBC points, and it avoids applying a land mask at the bay-like
                                ! intersection of orthogonal OBC segments.  Otherwise the calculation of terms
                                ! like the potential vorticity used in the barotropic solver relies on
                                ! bathymetry or other fields being projected outward across OBCs.  This option
                                ! changes answers for some configurations that use OBCs.
MASKING_DEPTH = -9999.0         !   [m] default = -9999.0
                                ! The depth below which to mask points as land points, for which all fluxes are
                                ! zeroed out. MASKING_DEPTH is ignored if it has the special default value.

! === module MOM_fixed_initialization ===
CHANNEL_CONFIG = "none"         ! 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.
SUBGRID_TOPO_AT_VEL = False     !   [Boolean] default = False
                                ! If true, use variables from TOPO_AT_VEL_FILE as parameters for porous barrier.
ROTATION = "beta"               ! default = "2omegasinlat"
                                ! This specifies how the Coriolis parameter is specified:
                                !     2omegasinlat - Use twice the planetary rotation rate
                                !       times the sine of latitude.
                                !     betaplane - Use a beta-plane or f-plane.
                                !     USER - call a user modified routine.
F_0 = 0.0                       !   [s-1] default = 0.0
                                ! The reference value of the Coriolis parameter with the betaplane option.
BETA = 0.0                      !   [m-1 s-1] default = 0.0
                                ! The northward gradient of the Coriolis parameter with the betaplane option.
BETA_LAT_REF = 0.0              !   [kilometers] default = 0.0
                                ! The reference latitude (origin) of the beta-plane
GRID_ROTATION_ANGLE_BUGS = False !   [Boolean] default = False
                                ! If true, use an older algorithm to calculate the sine and cosines needed
                                ! rotate between grid-oriented directions and true north and east.  Differences
                                ! arise at the tripolar fold.

! === module MOM_verticalGrid ===
! Parameters providing information about the vertical grid.
G_EARTH = 9.8                   !   [m s-2] default = 9.8
                                ! The gravitational acceleration of the Earth.
RHO_0 = 1031.0                  !   [kg m-3] default = 1035.0
                                ! The mean ocean density used with BOUSSINESQ true to calculate accelerations
                                ! and the mass for conservation properties, or with BOUSSINESQ false to convert
                                ! some parameters from vertical units of m to kg m-2.
BOUSSINESQ = True               !   [Boolean] default = True
                                ! If true, make the Boussinesq approximation.
ANGSTROM = 1.0E-10              !   [m] default = 1.0E-10
                                ! The minimum layer thickness, usually one-Angstrom.
H_TO_M = 1.0                    !   [m H-1] default = 1.0
                                ! A constant that translates the model's internal units of thickness into m.
NK = 10                         !   [nondim]
                                ! The number of model layers.

! === module MOM_tracer_registry ===

! === module MOM_restart ===
PARALLEL_RESTARTFILES = False   !   [Boolean] default = False
                                ! If true, the IO layout is used to group processors that write to the same
                                ! restart file or each processor writes its own (numbered) restart file. If
                                ! false, a single restart file is generated combining output from all PEs.
RESTARTFILE = "MOM.res"         ! default = "MOM.res"
                                ! The name-root of the restart file.
MAX_FIELDS = 100                ! default = 100
                                ! The maximum number of restart fields that can be used.
RESTART_CHECKSUMS_REQUIRED = True !   [Boolean] default = True
                                ! If true, require the restart checksums to match and error out otherwise. Users
                                ! may want to avoid this comparison if for example the restarts are made from a
                                ! run with a different mask_table than the current run, in which case the
                                ! checksums will not match and cause crash.
RESTART_SYMMETRIC_CHECKSUMS = False !   [Boolean] default = False
                                ! If true, do the restart checksums on all the edge points for a non-reentrant
                                ! grid.  This requires that SYMMETRIC_MEMORY_ is defined at compile time.
RESTART_UNSIGNED_ZEROS = False  !   [Boolean] default = False
                                ! If true, convert any negative zeros that would be written to the restart file
                                ! into ordinary unsigned zeros.  This does not change answers, but it can be
                                ! helpful in comparing restart files after grid rotation, for example.
USE_FILTER = False              !   [Boolean] default = False
                                ! If true, use streaming band-pass filters to detect the instantaneous tidal
                                ! signals in the simulation.

! === module MOM_tracer_flow_control ===
USE_USER_TRACER_EXAMPLE = False !   [Boolean] default = False
                                ! If true, use the USER_tracer_example tracer package.
USE_DOME_TRACER = True          !   [Boolean] default = False
                                ! If true, use the DOME_tracer tracer package.
USE_ISOMIP_TRACER = False       !   [Boolean] default = False
                                ! If true, use the ISOMIP_tracer tracer package.
USE_RGC_TRACER = False          !   [Boolean] default = False
                                ! If true, use the RGC_tracer tracer package.
USE_IDEAL_AGE_TRACER = False    !   [Boolean] default = False
                                ! If true, use the ideal_age_example tracer package.
USE_MARBL_TRACERS = False       !   [Boolean] default = False
                                ! If true, use the MARBL tracer package.
USE_REGIONAL_DYES = False       !   [Boolean] default = False
                                ! If true, use the regional_dyes tracer package.
USE_OIL_TRACER = False          !   [Boolean] default = False
                                ! If true, use the oil_tracer tracer package.
USE_ADVECTION_TEST_TRACER = False !   [Boolean] default = False
                                ! If true, use the advection_test_tracer tracer package.
USE_OCMIP2_CFC = False          !   [Boolean] default = False
                                ! If true, use the MOM_OCMIP2_CFC tracer package.
USE_CFC_CAP = False             !   [Boolean] default = False
                                ! If true, use the MOM_CFC_cap tracer package.
USE_generic_tracer = False      !   [Boolean] default = False
                                ! If true and _USE_GENERIC_TRACER is defined as a preprocessor macro, use the
                                ! MOM_generic_tracer packages.
USE_PSEUDO_SALT_TRACER = False  !   [Boolean] default = False
                                ! If true, use the pseudo salt tracer, typically run as a diagnostic.
USE_BOUNDARY_IMPULSE_TRACER = False !   [Boolean] default = False
                                ! If true, use the boundary impulse tracer.
USE_DYED_OBC_TRACER = False     !   [Boolean] default = False
                                ! If true, use the dyed_obc_tracer tracer package.
USE_NW2_TRACERS = False         !   [Boolean] default = False
                                ! If true, use the NeverWorld2 tracers.

! === module DOME_tracer ===
DOME_TRACER_IC_FILE = ""        ! default = ""
                                ! The name of a file from which to read the initial conditions for the DOME
                                ! tracers, or blank to initialize them internally.
DOME_TRACER_STRIPE_WIDTH = 50.0 !   [km] default = 50.0
                                ! The meridional width of the vertical stripes in the initial condition for the
                                ! DOME tracers.
DOME_TRACER_STRIPE_LAT = 350.0  !   [km] default = 350.0
                                ! The southern latitude of the first vertical stripe in the initial condition
                                ! for the DOME tracers.
DOME_TRACER_SHEET_SPACING = 600.0 !   [m] default = 600.0
                                ! The vertical spacing between successive horizontal sheets of tracer in the
                                ! initial conditions for the DOME tracers, and twice the thickness of these
                                ! tracer sheets.
SPONGE = False                  !   [Boolean] default = False
                                ! If true, sponges may be applied anywhere in the domain. The exact location and
                                ! properties of those sponges are specified from MOM_initialization.F90.

! === module MOM_boundary_update ===
OBC_VALUE_UPDATE_BUG = True     !   [Boolean] default = True
                                ! If true, recover a bug that OBC segment data does not update if all segments
                                ! use 'value' and none uses 'file'.
USE_FILE_OBC = False            !   [Boolean] default = False
                                ! If true, use external files for the open boundary.
USE_TIDAL_BAY_OBC = False       !   [Boolean] default = False
                                ! If true, use the tidal_bay open boundary.
USE_KELVIN_WAVE_OBC = False     !   [Boolean] default = False
                                ! If true, use the Kelvin wave open boundary.
USE_SHELFWAVE_OBC = False       !   [Boolean] default = False
                                ! If true, use the shelfwave open boundary.
USE_DYED_CHANNEL_OBC = False    !   [Boolean] default = False
                                ! If true, use the dyed channel open boundary.
OBC_SEGMENT_001_DATA = "U=value:0.0,V=value:0.0,SSH=value:0.0" !
                                ! OBC segment docs
OBC_SEGMENT_002_DATA = "U=value:0.0,V=value:0.0,SSH=value:0.0" !
                                ! OBC segment docs
OBC_SEGMENT_003_DATA = "U=value:0.0,V=value:0.0,SSH=value:0.0" !
                                ! OBC segment docs
OBC_SEGMENT_004_DATA = "U=value:0.0,V=value:0.0,SSH=value:0.0" !
                                ! OBC segment docs

! === module MOM_coord_initialization ===
COORD_CONFIG = "layer_ref"      ! default = "none"
                                ! This specifies how layers are to be defined:
                                !     ALE or none - used to avoid defining layers in ALE mode
                                !     file - read coordinate information from the file
                                !       specified by (COORD_FILE).
                                !     BFB - Custom coords for buoyancy-forced basin case
                                !       based on SST_S, T_BOT and DRHO_DT.
                                !     linear - linear based on interfaces not layers
                                !     layer_ref - linear based on layer densities
                                !     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.
GFS = 9.8                       !   [m s-2] default = 9.8
                                ! The reduced gravity at the free surface.
LIGHTEST_DENSITY = 1030.0       !   [kg m-3] default = 1031.0
                                ! The reference potential density used for layer 1.
DENSITY_RANGE = 2.0             !   [kg m-3] default = 2.0
                                ! The range of reference potential densities in the layers.

! === module MOM_state_initialization ===
INPUTDIR = "INPUT"              ! default = "."
                                ! The directory in which input files are found.
FATAL_INCONSISTENT_RESTART_TIME = False !   [Boolean] default = False
                                ! If true and a time_in value is provided to MOM_initialize_state, verify that
                                ! the time read from a restart file is the same as time_in, and issue a fatal
                                ! error if it is not.  Otherwise, simply set the time to time_in if present.
INIT_LAYERS_FROM_Z_FILE = False !   [Boolean] default = False
                                ! If true, initialize the layer thicknesses, temperatures, and salinities from a
                                ! Z-space file on a latitude-longitude grid.
THICKNESS_CONFIG = "circle_obcs" ! default = "uniform"
                                ! A string that determines how the initial layer thicknesses are specified for a
                                ! new run:
                                !     file - read interface heights from the file specified
                                !       by (THICKNESS_FILE).
                                !     thickness_file - read thicknesses from the file specified
                                !       by (THICKNESS_FILE).
                                !     mass_file - read thicknesses in units of mass per unit area from the file
                                !       specified by (THICKNESS_FILE).
                                !     coord - determined by ALE coordinate.
                                !     uniform - uniform thickness layers evenly distributed
                                !       between the surface and MAXIMUM_DEPTH.
                                !     list - read a list of positive interface depths.
                                !     param - use thicknesses from parameter THICKNESS_INIT_VALUES.
                                !     DOME - use a slope and channel configuration for the
                                !       DOME sill-overflow test case.
                                !     ISOMIP - use a configuration for the
                                !       ISOMIP test case.
                                !     benchmark - use the benchmark test case thicknesses.
                                !     Neverworld - use the Neverworld test case thicknesses.
                                !     search - search a density profile for the interface
                                !       densities. This is not yet implemented.
                                !     circle_obcs - the circle_obcs test case is used.
                                !     DOME2D - 2D version of DOME initialization.
                                !     adjustment2d - 2D lock exchange thickness ICs.
                                !     sloshing - sloshing gravity thickness ICs.
                                !     seamount - no motion test with seamount ICs.
                                !     dumbbell - sloshing channel ICs.
                                !     soliton - Equatorial Rossby soliton.
                                !     rossby_front - a mixed layer front in thermal wind balance.
                                !     USER - call a user modified routine.

! === module circle_obcs_initialization ===
DISK_RADIUS = 24.0              !   [km]
                                ! The radius of the initially elevated disk in the circle_obcs test case.
DISK_X_OFFSET = 0.0             !   [km] default = 0.0
                                ! The x-offset of the initially elevated disk in the circle_obcs test case.
DISK_IC_AMPLITUDE = 5.0         !   [m] default = 5.0
                                ! Initial amplitude of interface height displacements in the circle_obcs test
                                ! case.
DEPRESS_INITIAL_SURFACE = False !   [Boolean] default = False
                                ! If true,  depress the initial surface to avoid huge tsunamis when a large
                                ! surface pressure is applied.
TRIM_IC_FOR_P_SURF = False      !   [Boolean] default = False
                                ! If true, cuts way the top of the column for initial conditions at the depth
                                ! where the hydrostatic pressure matches the imposed surface pressure which is
                                ! read from file.
VELOCITY_CONFIG = "zero"        ! default = "zero"
                                ! A string that determines how the initial velocities are specified for a new
                                ! run:
                                !     file - read velocities from the file specified
                                !       by (VELOCITY_FILE).
                                !     zero - the fluid is initially at rest.
                                !     uniform - the flow is uniform (determined by
                                !       parameters INITIAL_U_CONST and INITIAL_V_CONST).
                                !     rossby_front - a mixed layer front in thermal wind balance.
                                !     soliton - Equatorial Rossby soliton.
                                !     USER - call a user modified routine.
ODA_INCUPD = False              !   [Boolean] default = False
                                ! If true, oda incremental updates will be applied everywhere in the domain.
OBC_RESERVOIR_INIT_BUG = True   !   [Boolean] default = True
                                ! If true, set the OBC tracer reservoirs at the startup of a new run from the
                                ! interior tracer concentrations regardless of properties that may be explicitly
                                ! specified for the reservoir concentrations.
OBC_USER_CONFIG = "none"        ! default = "none"
                                ! A string that sets how the user code is invoked to set open boundary data:
                                !    DOME - specified inflow on northern boundary
                                !    dyed_channel - supercritical with dye on the inflow boundary
                                !    dyed_obcs - circle_obcs with dyes on the open boundaries
                                !    Kelvin - barotropic Kelvin wave forcing on the western boundary
                                !    shelfwave - Flather with shelf wave forcing on western boundary
                                !    supercritical - now only needed here for the allocations
                                !    tidal_bay - Flather with tidal forcing on eastern boundary
                                !    USER - user specified

! === module MOM_diag_mediator ===
NUM_DIAG_COORDS = 1             ! default = 1
                                ! The number of diagnostic vertical coordinates to use. For each coordinate, an
                                ! entry in DIAG_COORDS must be provided.
DIAG_REMAPPING_USE_OM4_SUBCELLS = True !   [Boolean] default = True
                                ! If true, use the OM4 remapping-via-subcells algorithm for diagnostics. See
                                ! REMAPPING_USE_OM4_SUBCELLS for details. We recommend setting this option to
                                ! false.
USE_INDEX_DIAGNOSTIC_AXES = False !   [Boolean] default = False
                                ! If true, use a grid index coordinate convention for diagnostic axes.
DIAG_COORDS = "z Z ZSTAR"       ! 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           !   [various] 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 text file
                                ! 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 = "UNIFORM"    ! default = "UNIFORM"
                                ! 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)
                                !  WOA09INT[:N] - layers spanned by the WOA09 depths
                                !  WOA23INT[:N] - layers spanned by the WOA23 depths
                                !  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
                                !  HYBRID_3D:string - read from a file. The string specifies
                                !                the filename and two 3D variable names, separated
                                !                by a comma or space, for sigma-2 and dz. The
                                !                latter can be FNC1:string which is used everywhere.
                                !                e.g. HYBRID_3D:vgrid.nc,sigma2,dz
                                !  HYBRID_MAP:string - read from a file. The string specifies
                                !                the filename and three variable names, separated
                                !                by a comma or space, for map, sigma-2 and dz.
                                !                Map is a spatial index array with, maxval(map)=N,
                                !                and the others are 2D arrays containing N profiles.
                                !                Map typically contains integer values, but it can
                                !                contain real values, I+w, which imply using
                                !                the weighted sum of profiles I and I+1.
                                !                Dz can be FNC1:string which is used everywhere.
                                !                e.g. HYBRID_MAP:vgrid.nc,map,sigma2,dz

! === 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.

! === module MOM_lateral_mixing_coeffs ===
USE_VARIABLE_MIXING = False     !   [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.
USE_VISBECK = False             !   [Boolean] default = False
                                ! If true, use the Visbeck et al. (1997) formulation for
                                ! thickness diffusivity.
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.
DEPTH_SCALED_KHTH = False       !   [Boolean] default = False
                                ! If true, KHTH is scaled away when the depth is shallower than a reference
                                ! depth: KHTH = MIN(1,H/H0)**N * KHTH, where H0 is a reference depth, controlled
                                ! via DEPTH_SCALED_KHTH_H0, and the exponent (N) is controlled via
                                ! DEPTH_SCALED_KHTH_EXP.
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.
RESOLN_SCALED_KHTR = False      !   [Boolean] default = False
                                ! If true, the epipycnal tracer diffusivity is scaled away when the first
                                ! baroclinic deformation radius is well resolved.
RESOLN_USE_EBT = False          !   [Boolean] default = False
                                ! If true, uses the equivalent barotropic wave speed instead of first baroclinic
                                ! wave for calculating the resolution function.
BACKSCAT_EBT_POWER = 0.0        !   [nondim] default = 0.0
                                ! Power to raise EBT vertical structure to when backscatter has vertical
                                ! structure.
BS_USE_SQG_STRUCT = False       !   [Boolean] default = False
                                ! If true, the SQG vertical structure is used for backscatter on the condition
                                ! that BS_EBT_power=0
KHTH_USE_EBT_STRUCT = False     !   [Boolean] default = False
                                ! If true, uses the equivalent barotropic structure as the vertical structure of
                                ! thickness diffusivity.
KHTH_USE_SQG_STRUCT = False     !   [Boolean] default = False
                                ! If true, uses the surface quasigeostrophic structure as the vertical structure
                                ! of thickness diffusivity.
KHTR_USE_EBT_STRUCT = False     !   [Boolean] default = False
                                ! If true, uses the equivalent barotropic structure as the vertical structure of
                                ! tracer diffusivity.
KHTR_USE_SQG_STRUCT = False     !   [Boolean] default = False
                                ! If true, uses the surface quasigeostrophic structure as the vertical structure
                                ! of tracer diffusivity.
KD_GL90_USE_EBT_STRUCT = False  !   [Boolean] default = False
                                ! If true, uses the equivalent barotropic structure as the vertical structure of
                                ! diffusivity in the GL90 scheme.
KD_GL90_USE_SQG_STRUCT = False  !   [Boolean] default = False
                                ! If true, uses the equivalent barotropic structure as the vertical structure of
                                ! diffusivity in the GL90 scheme.
KHTH_SLOPE_CFF = 0.0            !   [nondim] default = 0.0
                                ! The nondimensional coefficient in the Visbeck formula for the interface depth
                                ! diffusivity
KHTR_SLOPE_CFF = 0.0            !   [nondim] default = 0.0
                                ! The nondimensional coefficient in the Visbeck formula for the epipycnal tracer
                                ! diffusivity
USE_STORED_SLOPES = False       !   [Boolean] default = False
                                ! If true, the isopycnal slopes are calculated once and stored for re-use. This
                                ! uses more memory but avoids calling the equation of state more times than
                                ! should be necessary.
VERY_SMALL_FREQUENCY = 1.0E-17  !   [s-1] default = 1.0E-17
                                ! A miniscule frequency that is used to avoid division by 0.  The default value
                                ! is roughly (pi / (the age of the universe)).
USE_STANLEY_ISO = False         !   [Boolean] default = False
                                ! If true, turn on Stanley SGS T variance parameterization in isopycnal slope
                                ! code.
MIXING_COEFS_OBC_BUG = True     !   [Boolean] default = True
                                ! If false, use only interior data for thickness weighting in lateral mixing
                                ! coefficient calculations and to calculate stratification and other fields at
                                ! open boundary condition faces.
RESOLN_FUNCTION_OBC_BUG = True  !   [Boolean] default = True
                                ! If false, use only interior data for calculating the resolution functions at
                                ! open boundary condition faces and vertices.
KD_SMOOTH = 1.0E-06             !   [m2 s-1] default = 1.0E-06
                                ! A diapycnal diffusivity that is used to interpolate more sensible values of T
                                ! & S into thin layers.
USE_QG_LEITH_GM = False         !   [Boolean] default = False
                                ! If true, use the QG Leith viscosity as the GM coefficient.

! === module MOM_set_visc ===
SET_VISC_ANSWER_DATE = 99991231 ! default = 99991231
                                ! The vintage of the order of arithmetic and expressions in the set viscosity
                                ! calculations.  Values below 20190101 recover the answers from the end of 2018,
                                ! while higher values use updated and more robust forms of the same expressions.
BOTTOMDRAGLAW = True            !   [Boolean] default = True
                                ! If true, the bottom stress is calculated with a drag law of the form
                                ! c_drag*|u|*u. The velocity magnitude may be an assumed value or it may be
                                ! based on the actual velocity in the bottommost HBBL, depending on LINEAR_DRAG.
DRAG_AS_BODY_FORCE = False      !   [Boolean] default = False
                                ! If true, the bottom stress is imposed as an explicit body force applied over a
                                ! fixed distance from the bottom, rather than as an implicit calculation based
                                ! on an enhanced near-bottom viscosity. The thickness of the bottom boundary
                                ! layer is HBBL.
CHANNEL_DRAG = False            !   [Boolean] default = False
                                ! If true, the bottom drag is exerted directly on each layer proportional to the
                                ! fraction of the bottom it overlies.
LINEAR_DRAG = False             !   [Boolean] default = False
                                ! If LINEAR_DRAG and BOTTOMDRAGLAW are defined the drag law is
                                ! cdrag*DRAG_BG_VEL*u.
DYNAMIC_VISCOUS_ML = False      !   [Boolean] default = False
                                ! If true, use a bulk Richardson number criterion to determine the mixed layer
                                ! thickness for viscosity.
OMEGA = 7.2921E-05              !   [s-1] default = 7.2921E-05
                                ! The rotation rate of the earth.
HBBL = 10.0                     !   [m]
                                ! The thickness of a bottom boundary layer with a viscosity increased by
                                ! KV_EXTRA_BBL 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.
CDRAG = 0.002                   !   [nondim] default = 0.003
                                ! CDRAG is the drag coefficient relating the magnitude of the velocity field to
                                ! the bottom stress. CDRAG is only used if BOTTOMDRAGLAW is defined.
CDRAG_MAP = False               !   [Boolean] default = False
                                ! If true, apply a spatially varying scaling factor to CDRAG, specified by
                                ! CDRAG_VAR in CDRAG_FILE.
BBL_USE_TIDAL_BG = False        !   [Boolean] default = False
                                ! Flag to use the tidal RMS amplitude in place of constant background velocity
                                ! for computing u* in the BBL. This flag is only used when BOTTOMDRAGLAW is true
                                ! and LINEAR_DRAG is false.
DRAG_BG_VEL = 0.05              !   [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_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.
HTBL_SHELF_MIN = 0.1            !   [m] default = 0.1
                                ! The minimum top boundary layer thickness that can be used with BOTTOMDRAGLAW.
                                ! This might be Kv/(cdrag*drag_bg_vel) to give Kv as the minimum near-top
                                ! viscosity.
HTBL_SHELF = 10.0               !   [m] default = 10.0
                                ! The thickness over which near-surface velocities are averaged for the drag law
                                ! under an ice shelf.  By default this is the same as HBBL
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.
KV_BBL_MIN = 1.0E-04            !   [m2 s-1] default = 1.0E-04
                                ! The minimum viscosities in the bottom boundary layer.
KV_TBL_MIN = 1.0E-04            !   [m2 s-1] default = 1.0E-04
                                ! The minimum viscosities in the top boundary layer.
CORRECT_BBL_BOUNDS = False      !   [Boolean] default = False
                                ! If true, uses the correct bounds on the BBL thickness and viscosity so that
                                ! the bottom layer feels the intended drag.

! === module MOM_thickness_diffuse ===
KHTH = 0.0                      !   [m2 s-1] default = 0.0
                                ! The background horizontal thickness diffusivity.
READ_KHTH = False               !   [Boolean] default = False
                                ! If true, read a file (given by KHTH_FILE) containing the spatially varying
                                ! horizontal isopycnal height diffusivity.
KHTH_MIN = 0.0                  !   [m2 s-1] default = 0.0
                                ! The minimum horizontal thickness diffusivity.
KHTH_MAX = 0.0                  !   [m2 s-1] default = 0.0
                                ! The maximum horizontal thickness diffusivity.
KHTH_MAX_CFL = 0.8              !   [nondimensional] default = 0.8
                                ! The maximum value of the local diffusive CFL ratio that is permitted for the
                                ! thickness diffusivity. 1.0 is the marginally unstable value in a pure layered
                                ! model, but much smaller numbers (e.g. 0.1) seem to work better for ALE-based
                                ! models.
KH_ETA_CONST = 0.0              !   [m2 s-1] default = 0.0
                                ! The background horizontal diffusivity of the interface heights (without
                                ! considering the layer density structure).  If diffusive CFL limits are
                                ! encountered, the diffusivities of the isopycnals and the interfaces heights
                                ! are scaled back proportionately.
KH_ETA_VEL_SCALE = 0.0          !   [m s-1] default = 0.0
                                ! A velocity scale that is multiplied by the grid spacing to give a contribution
                                ! to the horizontal diffusivity of the interface heights (without considering
                                ! the layer density structure).
DETANGLE_INTERFACES = False     !   [Boolean] default = False
                                ! If defined add 3-d structured enhanced interface height diffusivities to
                                ! horizontally smooth jagged layers.
KHTH_SLOPE_MAX = 0.01           !   [nondim] default = 0.01
                                ! A slope beyond which the calculated isopycnal slope is not reliable and is
                                ! scaled away.
KHTH_USE_FGNV_STREAMFUNCTION = False !   [Boolean] default = False
                                ! If true, use the streamfunction formulation of Ferrari et al., 2010, which
                                ! effectively emphasizes graver vertical modes by smoothing in the vertical.
USE_STANLEY_GM = False          !   [Boolean] default = False
                                ! If true, turn on Stanley SGS T variance parameterization in GM code.
MEKE_GM_SRC_ALT = False         !   [Boolean] default = False
                                ! If true, use the GM energy conversion form S^2*N^2*kappa rather than the
                                ! streamfunction for the GM source term.
MEKE_GEOMETRIC = False          !   [Boolean] default = False
                                ! If true, uses the GM coefficient formulation from the GEOMETRIC framework
                                ! (Marshall et al., 2012).
USE_GME = False                 !   [Boolean] default = False
                                ! If true, use the GM+E backscatter scheme in association with the Gent and
                                ! McWilliams parameterization.
USE_GM_WORK_BUG = False         !   [Boolean] default = False
                                ! If true, compute the top-layer work tendency on the u-grid with the incorrect
                                ! sign, for legacy reproducibility.

! === module MOM_dynamics_split_RK2 ===
TIDES = False                   !   [Boolean] default = False
                                ! If true, apply tidal momentum forcing.
CALCULATE_SAL = False           !   [Boolean] default = False
                                ! If true, calculate self-attraction and loading.
USE_HA = False                  !   [Boolean] default = False
                                ! If true, perform inline harmonic analysis.
BE = 0.6                        !   [nondim] default = 0.6
                                ! If SPLIT is true, BE determines the relative weighting of a 2nd-order
                                ! Runga-Kutta baroclinic time stepping scheme (0.5) and a backward Euler scheme
                                ! (1) that is used for the Coriolis and inertial terms.  BE may be from 0.5 to
                                ! 1, but instability may occur near 0.5. BE is also applicable if SPLIT is false
                                ! and USE_RK2 is true.
BEGW = 0.0                      !   [nondim] default = 0.0
                                ! If SPLIT is true, BEGW is a number from 0 to 1 that controls the extent to
                                ! which the treatment of gravity waves is forward-backward (0) or simulated
                                ! backward Euler (1).  0 is almost always used. If SPLIT is false and USE_RK2 is
                                ! true, BEGW can be between 0 and 0.5 to damp gravity waves.
SET_DTBT_USE_BT_CONT = False    !   [Boolean] default = False
                                ! If true, use BT_CONT to calculate DTBT if possible.
SPLIT_BOTTOM_STRESS = False     !   [Boolean] default = False
                                ! If true, provide the bottom stress calculated by the vertical viscosity to the
                                ! barotropic solver.
BT_USE_LAYER_FLUXES = True      !   [Boolean] default = True
                                ! If true, use the summed layered fluxes plus an adjustment due to the change in
                                ! the barotropic velocity in the barotropic continuity equation.
BT_ADJ_CORR_MASS_SRC = True     !   [Boolean] default = True
                                ! If true, recalculates the barotropic mass source after predictor step. This
                                ! should make little difference in the deep ocean but appears to help for
                                ! vanished layers. If false, uses the same mass source as from the predictor
                                ! step.
STORE_CORIOLIS_ACCEL = True     !   [Boolean] default = True
                                ! If true, calculate the Coriolis accelerations at the end of each timestep for
                                ! use in the predictor step of the next split RK2 timestep.
FPMIX = False                   !   [Boolean] default = False
                                ! If true, add non-local momentum flux increments and diffuse down the Eulerian
                                ! gradient.
VISC_REM_BUG = True             !   [Boolean] default = False
                                ! If true, visc_rem_[uv] in split mode is incorrectly calculated or accounted
                                ! for in two places. This parameter controls the defaults of two individual
                                ! flags, VISC_REM_TIMESTEP_BUG in MOM_dynamics_split_RK2(b) and
                                ! VISC_REM_BT_WEIGHT_BUG in MOM_barotropic.
VISC_REM_TIMESTEP_BUG = True    !   [Boolean] default = True
                                ! If true, recover a bug that uses dt_pred rather than dt in vertvisc_remnant()
                                ! at the end of predictor stage for the following continuity() and btstep()
                                ! calls in the corrector step. Default of this flag is set by VISC_REM_BUG

! === module MOM_continuity_PPM ===
MONOTONIC_CONTINUITY = False    !   [Boolean] default = False
                                ! If true, CONTINUITY_PPM uses the Colella and Woodward monotonic limiter.  The
                                ! default (false) is to use a simple positive definite limiter.
SIMPLE_2ND_PPM_CONTINUITY = False !   [Boolean] default = False
                                ! If true, CONTINUITY_PPM uses a simple 2nd order (arithmetic mean)
                                ! interpolation of the edge values. This may give better PV conservation
                                ! properties. While it formally reduces the accuracy of the continuity solver
                                ! itself in the strongly advective limit, it does not reduce the overall order
                                ! of accuracy of the dynamic core.
UPWIND_1ST_CONTINUITY = False   !   [Boolean] default = False
                                ! If true, CONTINUITY_PPM becomes a 1st-order upwind continuity solver.  This
                                ! scheme is highly diffusive but may be useful for debugging or in single-column
                                ! mode where its minimal stencil is useful.
ETA_TOLERANCE = 5.0E-10         !   [m] default = 5.0E-10
                                ! 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 than about
                                ! 10^-15*MAXIMUM_DEPTH.
VELOCITY_TOLERANCE = 3.0E+08    !   [m s-1] default = 3.0E+08
                                ! The tolerance for barotropic velocity discrepancies between the barotropic
                                ! solution and  the sum of the layer thicknesses.
CONT_PPM_AGGRESS_ADJUST = False !   [Boolean] default = False
                                ! If true, allow the adjusted velocities to have a relative CFL change up to
                                ! 0.5.
CONT_PPM_VOLUME_BASED_CFL = False !   [Boolean] default = False
                                ! If true, use the ratio of the open face lengths to the tracer cell areas when
                                ! estimating CFL numbers.  The default is set by CONT_PPM_AGGRESS_ADJUST.
CONTINUITY_CFL_LIMIT = 0.5      !   [nondim] default = 0.5
                                ! The maximum CFL of the adjusted velocities.
CONT_PPM_BETTER_ITER = True     !   [Boolean] default = True
                                ! If true, stop corrective iterations using a velocity based criterion and only
                                ! stop if the iteration is better than all predecessors.
CONT_PPM_USE_VISC_REM_MAX = True !   [Boolean] default = True
                                ! If true, use more appropriate limiting bounds for corrections in strongly
                                ! viscous columns.
CONT_PPM_MARGINAL_FACE_AREAS = True !   [Boolean] default = True
                                ! If true, use the marginal face areas from the continuity solver for use as the
                                ! weights in the barotropic solver. Otherwise use the transport averaged areas.
CONT_USE_H_MARG_MIN = False     !   [Boolean] default = False
                                ! If true, the marginal thickness used and returned from continuity is bounded
                                ! from below by a sub-roundoff value. Otherwise the minimum is 0.

! === module MOM_CoriolisAdv ===
NOSLIP = False                  !   [Boolean] default = False
                                ! If true, no slip boundary conditions are used; otherwise free slip boundary
                                ! conditions are assumed. The implementation of the free slip BCs on a C-grid is
                                ! much cleaner than the no slip BCs. The use of free slip BCs is strongly
                                ! encouraged, and no slip BCs are not used with the biharmonic viscosity.
CORIOLIS_EN_DIS = False         !   [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.
CORIOLIS_SCHEME = "SADOURNY75_ENERGY" ! default = "SADOURNY75_ENERGY"
                                ! CORIOLIS_SCHEME selects the discretization for the Coriolis terms. Valid
                                ! values are:
                                !    SADOURNY75_ENERGY - Sadourny, 1975; energy cons.
                                !    ARAKAWA_HSU90     - Arakawa & Hsu, 1990
                                !    SADOURNY75_ENSTRO - Sadourny, 1975; enstrophy cons.
                                !    ARAKAWA_LAMB81    - Arakawa & Lamb, 1981; En. + Enst.
                                !    ARAKAWA_LAMB_BLEND - A blend of Arakawa & Lamb with
                                !                         Arakawa & Hsu and Sadourny energy
                                !    WENOVI5TH_PV_ENSTRO   - 5th-order WENO PV enstrophy
                                !    WENOVI3RD_PV_ENSTRO   - 3rd-order WENO PV enstrophy
                                !    WENOVI7TH_PV_ENSTRO  - 7th-order WENO PV enstrophy
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 would have no effect on the SADOURNY Coriolis scheme if it were
                                ! possible to use centered difference thickness fluxes.
KE_SCHEME = "KE_ARAKAWA"        ! default = "KE_ARAKAWA"
                                ! KE_SCHEME selects the discretization for acceleration due to the kinetic
                                ! energy gradient. Valid values are:
                                !    KE_ARAKAWA, KE_SIMPLE_GUDONOV, KE_GUDONOV, KE_UP3
PV_ADV_SCHEME = "PV_ADV_CENTERED" ! default = "PV_ADV_CENTERED"
                                ! PV_ADV_SCHEME selects the discretization for PV advection. Valid values are:
                                !    PV_ADV_CENTERED - centered (aka Sadourny, 75)
                                !    PV_ADV_UPWIND1  - upwind, first order

! === module MOM_PressureForce ===
ANALYTIC_FV_PGF = True          !   [Boolean] default = True
                                ! If true the pressure gradient forces are calculated with a finite volume form
                                ! that analytically integrates the equations of state in pressure to avoid any
                                ! possibility of numerical thermobaric instability, as described in Adcroft et
                                ! al., O. Mod. (2008).

! === module MOM_PressureForce_FV ===
RHO_PGF_REF = 1031.0            !   [kg m-3] default = 1031.0
                                ! The reference density that is subtracted off when calculating pressure
                                ! gradient forces.  Its inverse is subtracted off of specific volumes when in
                                ! non-Boussinesq mode.  The default is RHO_0.
RHO_PGF_REF_BUG = True          !   [Boolean] default = True
                                ! If true, recover a bug that RHO_0 (the mean seawater density in Boussinesq
                                ! mode) and RHO_PGF_REF (the subtracted reference density in finite volume
                                ! pressure gradient forces) are incorrectly interchanged in several instances in
                                ! Boussinesq mode.
SSH_IN_EOS_PRESSURE_FOR_PGF = False !   [Boolean] default = False
                                ! If true, include contributions from the sea surface height in the height-based
                                ! pressure used in the equation of state calculations for the Boussinesq
                                ! pressure gradient forces, including adjustments for atmospheric or sea-ice
                                ! pressure.
MASS_WEIGHT_IN_PRESSURE_GRADIENT = False !   [Boolean] default = False
                                ! If true, use mass weighting when interpolating T/S for integrals near the
                                ! bathymetry in FV pressure gradient calculations.
MASS_WEIGHT_IN_PRESSURE_GRADIENT_TOP = False !   [Boolean] default = False
                                ! If true and MASS_WEIGHT_IN_PRESSURE_GRADIENT is true, use mass weighting when
                                ! interpolating T/S for integrals near the top of the water column in FV
                                ! pressure gradient calculations.
MASS_WEIGHT_IN_PGF_VANISHED_ONLY = False !   [Boolean] default = False
                                ! If true, use mass weighting when interpolating T/S for integrals only if one
                                ! side is vanished according to RESET_INTXPA_H_NONVANISHED.
USE_INACCURATE_PGF_RHO_ANOM = False !   [Boolean] default = False
                                ! If true, use a form of the PGF that uses the reference density in an
                                ! inaccurate way. This is not recommended.
RECONSTRUCT_FOR_PRESSURE = False !   [Boolean] default = False
                                ! If True, use vertical reconstruction of T & S within the integrals of the FV
                                ! pressure gradient calculation. If False, use the constant-by-layer algorithm.
                                ! The default is set by USE_REGRIDDING.
PRESSURE_RECONSTRUCTION_SCHEME = 1 ! default = 1
                                ! Order of vertical reconstruction of T/S to use in the integrals within the FV
                                ! pressure gradient calculation.
                                !  0: PCM or no reconstruction.
                                !  1: PLM reconstruction.
                                !  2: PPM reconstruction.
                                !  3: PLM with least squares slope.
BOUNDARY_EXTRAPOLATION_PRESSURE = True !   [Boolean] default = True
                                ! If true, the reconstruction of T & S for pressure in boundary cells is
                                ! extrapolated, rather than using PCM in these cells. If true, the same order
                                ! polynomial is used as is used for the interior cells.
USE_STANLEY_PGF = False         !   [Boolean] default = False
                                ! If true, turn on Stanley SGS T variance parameterization in PGF code.

! === module MOM_Zanna_Bolton ===
USE_ZB2020 = False              !   [Boolean] default = False
                                ! If true, turns on Zanna-Bolton-2020 (ZB) subgrid momentum parameterization of
                                ! mesoscale eddies.

! === module MOM_hor_visc ===
USE_CIRCULATION_IN_HORVISC = False !   [Boolean] default = False
                                ! Use circulation theorem to compute vorticity in horvisc module (for ZB20 or
                                ! Leith)
USE_CONT_THICKNESS = False      !   [Boolean] default = False
                                ! If true, use thickness at velocity points from continuity solver. This option
                                ! currently only works with split mode.
LAPLACIAN = True                !   [Boolean] default = False
                                ! If true, use a Laplacian horizontal viscosity.
KH = 25.0                       !   [m2 s-1] default = 0.0
                                ! The background Laplacian horizontal viscosity.
KH_BG_MIN = 0.0                 !   [m2 s-1] default = 0.0
                                ! The minimum value allowed for Laplacian horizontal viscosity, KH.
KH_VEL_SCALE = 0.003            !   [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 and Leith viscosities, and KH.
KH_SIN_LAT = 0.0                !   [m2 s-1] default = 0.0
                                ! The amplitude of a latitudinally-dependent background viscosity of the form
                                ! KH_SIN_LAT*(SIN(LAT)**KH_PWR_OF_SINE).
SMAGORINSKY_KH = True           !   [Boolean] default = False
                                ! If true, use a Smagorinsky nonlinear eddy viscosity.
SMAG_LAP_CONST = 0.15           !   [nondim] default = 0.0
                                ! The nondimensional Laplacian Smagorinsky constant, often 0.15.
LEITH_KH = False                !   [Boolean] default = False
                                ! If true, use a Leith nonlinear eddy viscosity.
BOUND_KH = True                 !   [Boolean] default = True
                                ! If true, the Laplacian coefficient is locally limited to be stable.
EY24_EBT_BS = False             !   [Boolean] default = False
                                ! If true, use the backscatter scheme (EBT mode with kill switch) developed by
                                ! Yankovsky et al. (2024).
ANISOTROPIC_VISCOSITY = False   !   [Boolean] default = False
                                ! If true, allow anistropic viscosity in the Laplacian horizontal viscosity.
ADD_LES_VISCOSITY = False       !   [Boolean] default = False
                                ! If true, adds the viscosity from Smagorinsky and Leith to the background
                                ! viscosity instead of taking the maximum.
BIHARMONIC = True               !   [Boolean] default = True
                                ! If true, use a biharmonic horizontal viscosity. BIHARMONIC may be used with
                                ! LAPLACIAN.
AH = 0.0                        !   [m4 s-1] default = 0.0
                                ! The background biharmonic horizontal viscosity.
AH_VEL_SCALE = 0.003            !   [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 and Leith viscosities, and AH.
AH_TIME_SCALE = 0.0             !   [s] default = 0.0
                                ! A time scale whose inverse is multiplied by the fourth power of the grid
                                ! spacing to calculate biharmonic viscosity. The final viscosity is the largest
                                ! of all viscosity formulations in use. 0.0 means that it's not used.
SMAGORINSKY_AH = True           !   [Boolean] default = False
                                ! If true, use a biharmonic Smagorinsky nonlinear eddy viscosity.
LEITH_AH = False                !   [Boolean] default = False
                                ! If true, use a biharmonic Leith nonlinear eddy viscosity.
USE_LEITHY = False              !   [Boolean] default = False
                                ! If true, use a biharmonic Leith nonlinear eddy viscosity together with a
                                ! harmonic backscatter.
BOUND_AH = True                 !   [Boolean] default = True
                                ! If true, the biharmonic coefficient is locally limited to be stable.
RE_AH = 0.0                     !   [nondim] default = 0.0
                                ! If nonzero, the biharmonic coefficient is scaled so that the biharmonic
                                ! Reynolds number is equal to this.
BACKSCATTER_UNDERBOUND = False  !   [Boolean] default = False
                                ! If true, the bounds on the biharmonic viscosity are allowed to increase where
                                ! the Laplacian viscosity is negative (due to backscatter parameterizations)
                                ! beyond the largest timestep-dependent stable values of biharmonic viscosity
                                ! when no Laplacian viscosity is applied.  The default is true for historical
                                ! reasons, but this option probably should not be used because it can contribute
                                ! to numerical instabilities.
SMAG_BI_CONST = 0.06            !   [nondim] default = 0.0
                                ! The nondimensional biharmonic Smagorinsky constant, typically 0.015 - 0.06.
BOUND_CORIOLIS_BIHARM = True    !   [Boolean] default = True
                                ! If true use a viscosity that increases with the square of the velocity shears,
                                ! so that the resulting viscous drag is of comparable magnitude to the Coriolis
                                ! terms when the velocity differences between adjacent grid points is
                                ! 0.5*BOUND_CORIOLIS_VEL.  The default is the value of BOUND_CORIOLIS (or
                                ! false).
BOUND_CORIOLIS_VEL = 10.0       !   [m s-1] default = 10.0
                                ! The velocity scale at which BOUND_CORIOLIS_BIHARM causes the biharmonic drag
                                ! to have comparable magnitude to the Coriolis acceleration.  The default is set
                                ! by MAXVEL.
USE_LAND_MASK_FOR_HVISC = True  !   [Boolean] default = True
                                ! If true, use the land mask for the computation of thicknesses at velocity
                                ! locations. This eliminates the dependence on arbitrary values over land or
                                ! outside of the domain.
HORVISC_BOUND_COEF = 0.8        !   [nondim] default = 0.8
                                ! The nondimensional coefficient of the ratio of the viscosity bounds to the
                                ! theoretical maximum for stability without considering other terms.
USE_KH_BG_2D = False            !   [Boolean] default = False
                                ! If true, read a file containing 2-d background harmonic viscosities. The final
                                ! viscosity is the maximum of the other terms and this background value.
FRICTWORK_BUG = False           !   [Boolean] default = False
                                ! If true, retain an answer-changing bug in calculating the FrictWork, which
                                ! cancels the h in thickness flux and the h at velocity point. This is not
                                ! recommended.
OBC_SPECIFIED_STRAIN_BUG = True !   [Boolean] default = True
                                ! If true, recover a bug that specified shear strain option at open boundaries
                                ! cannot be applied.

! === module MOM_vert_friction ===
VERT_FRICTION_ANSWER_DATE = 99991231 ! default = 99991231
                                ! The vintage of the order of arithmetic and expressions in the viscous
                                ! calculations.  Values below 20190101 recover the answers from the end of 2018,
                                ! while higher values use expressions that do not use an arbitrary hard-coded
                                ! maximum viscous coupling coefficient between layers.  Values below 20230601
                                ! recover a form of the viscosity within the mixed layer that breaks up the
                                ! magnitude of the wind stress in some non-Boussinesq cases.
DIRECT_STRESS = True            !   [Boolean] default = False
                                ! If true, the wind stress is distributed over the topmost HMIX_STRESS of fluid
                                ! (like in HYCOM), and an added mixed layer viscosity or a physically based
                                ! boundary layer turbulence parameterization is not needed for stability.
FIXED_DEPTH_LOTW_ML = False     !   [Boolean] default = False
                                ! If true, use a Law-of-the-wall prescription for the mixed layer viscosity
                                ! within a boundary layer that is the lesser of HMIX_FIXED and the total depth
                                ! of the ocean in a column.
LOTW_VISCOUS_ML_FLOOR = False   !   [Boolean] default = False
                                ! If true, use a Law-of-the-wall prescription to set a lower bound on the
                                ! viscous coupling between layers within the surface boundary layer, based the
                                ! distance of interfaces from the surface.  This only acts when there are large
                                ! changes in the thicknesses of successive layers or when the viscosity is set
                                ! externally and the wind stress has subsequently increased.
VON_KARMAN_CONST = 0.41         !   [nondim] default = 0.41
                                ! The value the von Karman constant as used for mixed layer viscosity.
HARMONIC_VISC = True            !   [Boolean] default = False
                                ! If true, use the harmonic mean thicknesses for calculating the vertical
                                ! viscosity.
HARMONIC_BL_SCALE = 0.0         !   [nondim] default = 0.0
                                ! A scale to determine when water is in the boundary layers based solely on
                                ! harmonic mean thicknesses for the purpose of determining the extent to which
                                ! the thicknesses used in the viscosities are upwinded.
HMIX_FIXED = 20.0               !   [m]
                                ! The prescribed depth over which the near-surface viscosity and diffusivity are
                                ! elevated when the bulk mixed layer is not used.
HMIX_STRESS = 20.0              !   [m] default = 20.0
                                ! The depth over which the wind stress is applied if DIRECT_STRESS is true.
USE_GL90_IN_SSW = False         !   [Boolean] default = False
                                ! If true, use simpler method to calculate 1/N^2 in GL90 vertical viscosity
                                ! coefficient. This method is valid in stacked shallow water mode.
KV_ML_INVZ2 = 0.0               !   [m2 s-1] default = 0.0
                                ! An extra kinematic viscosity in a mixed layer of thickness HMIX_FIXED, with
                                ! the actual viscosity scaling as 1/(z*HMIX_FIXED)^2, where z is the distance
                                ! from the surface, to allow for finite wind stresses to be transmitted through
                                ! infinitesimally thin surface layers.  This is an older option for numerical
                                ! convenience without a strong physical basis, and its use is now discouraged.
CFL_TRUNCATE = 0.5              !   [nondim] default = 0.5
                                ! The value of the CFL number that will cause velocity components to be
                                ! truncated; instability can occur past 0.5.
CFL_REPORT = 0.5                !   [nondim] default = 0.5
                                ! The value of the CFL number that causes accelerations to be reported; the
                                ! default is CFL_TRUNCATE.
CFL_TRUNCATE_RAMP_TIME = 0.0    !   [s] default = 0.0
                                ! The time over which the CFL truncation value is ramped up at the beginning of
                                ! the run.
CFL_TRUNCATE_START = 0.0        !   [nondim] default = 0.0
                                ! The start value of the truncation CFL number used when ramping up CFL_TRUNC.
STOKES_MIXING_COMBINED = False  !   [Boolean] default = False
                                ! Flag to use Stokes drift Mixing via the Lagrangian  current (Eulerian plus
                                ! Stokes drift).  Still needs work and testing, so not recommended for use.
VEL_UNDERFLOW = 0.0             !   [m s-1] default = 0.0
                                ! A negligibly small velocity magnitude below which velocity components are set
                                ! to 0.  A reasonable value might be 1e-30 m/s, which is less than an Angstrom
                                ! divided by the age of the universe.

! === module MOM_barotropic ===
USE_BT_CONT_TYPE = True         !   [Boolean] default = True
                                ! If true, use a structure with elements that describe effective face areas from
                                ! the summed continuity solver as a function the barotropic flow in coupling
                                ! between the barotropic and baroclinic flow.  This is only used if SPLIT is
                                ! true.
INTEGRAL_BT_CONTINUITY = False  !   [Boolean] default = False
                                ! If true, use the time-integrated velocity over the barotropic steps to
                                ! determine the integrated transports used to update the continuity equation.
                                ! Otherwise the transports are the sum of the transports based on a series of
                                ! instantaneous velocities and the BT_CONT_TYPE for transports.  This is only
                                ! valid if USE_BT_CONT_TYPE = True.
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 maxCFL_BT_cont to be accommodated.
BT_CONT_CORR_BOUNDS = True      !   [Boolean] default = True
                                ! If true, and BOUND_BT_CORRECTION is true, use the BT_cont_type variables to
                                ! set limits determined by MAXCFL_BT_CONT on the CFL number of the velocities
                                ! that are likely to be driven by the corrective mass fluxes.
ADJUST_BT_CONT = False          !   [Boolean] default = False
                                ! If true, adjust the curve fit to the BT_cont type that is used by the
                                ! barotropic solver to match the transport about which the flow is being
                                ! linearized.
GRADUAL_BT_ICS = False          !   [Boolean] default = False
                                ! If true, adjust the initial conditions for the barotropic solver to the values
                                ! from the layered solution over a whole timestep instead of instantly. This is
                                ! a decent approximation to the inclusion of sum(u dh_dt) while also correcting
                                ! for truncation errors.
BT_USE_VISC_REM_U_UH0 = False   !   [Boolean] default = False
                                ! If true, use the viscous remnants when estimating the barotropic velocities
                                ! that were used to calculate uh0 and vh0.  False is probably the better choice.
NONLIN_BT_CONT_UPDATE_PERIOD = 1 ! default = 1
                                ! If NONLINEAR_BT_CONTINUITY is true, this is the number of barotropic time
                                ! steps between updates to the face areas, or 0 to update only before the
                                ! barotropic stepping.
BT_PROJECT_VELOCITY = True      !   [Boolean] default = False
                                ! If true, step the barotropic velocity first and project out the velocity
                                ! tendency 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_NONLIN_STRESS = False        !   [Boolean] default = False
                                ! If true, use the full depth of the ocean at the start of the barotropic step
                                ! when calculating the surface stress contribution to the barotropic
                                ! acclerations.  Otherwise use the depth based on bathyT.
DYNAMIC_SURFACE_PRESSURE = False !   [Boolean] default = False
                                ! If true, add a dynamic pressure due to a viscous ice shelf, for instance.
BT_CORIOLIS_SCALE = 1.0         !   [nondim] default = 1.0
                                ! A factor by which the barotropic Coriolis anomaly terms are scaled.
BAROTROPIC_ANSWER_DATE = 99991231 ! default = 99991231
                                ! The vintage of the expressions in the barotropic solver. Values below 20190101
                                ! recover the answers from the end of 2018, while higher values use more
                                ! efficient or general expressions.
VISC_REM_BT_WEIGHT_BUG = True   !   [Boolean] default = True
                                ! If true, recover a bug in barotropic solver that uses an unnormalized weight
                                ! function for vertical averages of baroclinic velocity and forcing. Default of
                                ! this flag is set by VISC_REM_BUG.
TIDAL_SAL_FLATHER = False       !   [Boolean] default = False
                                ! If true, then apply adjustments to the external gravity wave speed used with
                                ! the Flather OBC routine consistent with the barotropic solver. This applies to
                                ! cases with  tidal forcing using the scalar self-attraction approximation. The
                                ! default is currently False in order to retain previous answers but should be
                                ! set to True for new experiments
SADOURNY = True                 !   [Boolean] default = True
                                ! If true, the Coriolis terms are discretized with the Sadourny (1975) energy
                                ! conserving scheme, otherwise the Arakawa & Hsu scheme is used.  If the
                                ! internal deformation radius is not resolved, the Sadourny scheme should
                                ! probably be used.
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 (the default) - 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
BT_STRONG_DRAG = False          !   [Boolean] default = False
                                ! If true, use a stronger estimate of the retarding effects of strong bottom
                                ! drag, by making it implicit with the barotropic time-step instead of implicit
                                ! with the baroclinic time-step and dividing by the number of barotropic steps.
BT_LINEAR_WAVE_DRAG = False     !   [Boolean] default = False
                                ! If true, apply a linear drag to the barotropic velocities, using rates set by
                                ! lin_drag_u & _v divided by the depth of the ocean.  This was introduced to
                                ! facilitate tide modeling.
BT_LINEAR_FREQ_DRAG = False     !   [Boolean] default = False
                                ! If true, apply frequency-dependent drag to the tidal velocities. The streaming
                                ! band-pass filter must be turned on.
CLIP_BT_VELOCITY = False        !   [Boolean] default = False
                                ! If true, limit any velocity components that exceed CFL_TRUNCATE.  This should
                                ! only be used as a desperate debugging measure.
MAXCFL_BT_CONT = 0.25           !   [nondim] default = 0.25
                                ! The maximum permitted CFL number associated with the barotropic accelerations
                                ! from the summed velocities times the time-derivatives of thicknesses.
DT_BT_FILTER = -0.25            !   [sec or nondim] default = -0.25
                                ! A time-scale over which the barotropic mode solutions are filtered, in seconds
                                ! if positive, or as a fraction of DT if negative. When used this can never be
                                ! taken to be longer than 2*dt.  Set this to 0 to apply no filtering.
G_BT_EXTRA = 0.0                !   [nondim] default = 0.0
                                ! A nondimensional factor by which gtot is enhanced.
SSH_EXTRA = 10.0                !   [m] default = 10.0
                                ! An estimate of how much higher SSH might get, for use in calculating the safe
                                ! external wave speed. The default is the minimum of 10 m or 5% of
                                ! MAXIMUM_DEPTH.
LINEARIZED_BT_CORIOLIS = True   !   [Boolean] default = True
                                ! If true use the bottom depth instead of the total water column thickness in
                                ! the barotropic Coriolis term calculations.
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.
BT_USE_OLD_CORIOLIS_BRACKET_BUG = False !   [Boolean] default = False
                                ! If True, use an order of operations that is not bitwise rotationally symmetric
                                ! in the meridional Coriolis term of the barotropic solver.
MASK_COASTAL_PRESSURE_FORCE = False !   [Boolean] default = False
                                ! If true, use the land masks to zero out the diagnosed barotropic pressure
                                ! gradient accelerations at coastal or land points.  This changes diagnostics
                                ! and improves the reproducibility of certain debugging checksums, but it does
                                ! not alter the solutions themselves.

! === module MOM_mixed_layer_restrat ===
MIXEDLAYER_RESTRAT = False      !   [Boolean] default = False
                                ! If true, a density-gradient dependent re-stratifying flow is imposed in the
                                ! mixed layer. Can be used in ALE mode without restriction but in layer mode can
                                ! only be used if BULKMIXEDLAYER is true.

! === module MOM_diagnostics ===
DIAG_EBT_MONO_N2_COLUMN_FRACTION = 0.0 !   [nondim] default = 0.0
                                ! The lower fraction of water column over which N2 is limited as monotonic for
                                ! the purposes of calculating the equivalent barotropic wave speed.
DIAG_EBT_MONO_N2_DEPTH = -1.0   !   [m] default = -1.0
                                ! The depth below which N2 is limited as monotonic for the purposes of
                                ! calculating the equivalent barotropic wave speed.
INTERNAL_WAVE_SPEED_TOL = 0.001 !   [nondim] default = 0.001
                                ! The fractional tolerance for finding the wave speeds.
INTERNAL_WAVE_SPEED_MIN = 0.0   !   [m s-1] default = 0.0
                                ! A floor in the first mode speed below which 0 used instead.
INTERNAL_WAVE_SPEED_BETTER_EST = True !   [Boolean] default = True
                                ! If true, use a more robust estimate of the first mode wave speed as the
                                ! starting point for iterations.
INTWAVE_REMAPPING_USE_OM4_SUBCELLS = True !   [Boolean] default = True
                                ! If true, use the OM4 remapping-via-subcells algorithm for calculating EBT
                                ! structure. See REMAPPING_USE_OM4_SUBCELLS for details. We recommend setting
                                ! this option to false.

! === module MOM_diabatic_driver ===
! The following parameters are used for diabatic processes.
USE_LEGACY_DIABATIC_DRIVER = True !   [Boolean] default = True
                                ! If true, use a legacy version of the diabatic subroutine. This is temporary
                                ! and is needed to avoid change in answers.
ENERGETICS_SFC_PBL = False      !   [Boolean] default = False
                                ! If true, use an implied energetics planetary boundary layer scheme to
                                ! determine the diffusivity and viscosity in the surface boundary layer.
EPBL_IS_ADDITIVE = True         !   [Boolean] default = True
                                ! If true, the diffusivity from ePBL is added to all other diffusivities.
                                ! Otherwise, the larger of kappa-shear and ePBL diffusivities are used.
INTERNAL_TIDES = False          !   [Boolean] default = False
                                ! If true, use the code that advances a separate set of equations for the
                                ! internal tide energy density.
MASSLESS_MATCH_TARGETS = True   !   [Boolean] default = True
                                ! If true, the temperature and salinity of massless layers are kept consistent
                                ! with their target densities. Otherwise the properties of massless layers
                                ! evolve diffusively to match massive neighboring layers.
AGGREGATE_FW_FORCING = True     !   [Boolean] default = True
                                ! If true, the net incoming and outgoing fresh water fluxes are combined and
                                ! applied as either incoming or outgoing depending on the sign of the net. If
                                ! false, the net incoming fresh water flux is added to the model and thereafter
                                ! the net outgoing is removed from the topmost non-vanished layers of the
                                ! updated state.
MIX_BOUNDARY_TRACERS = True     !   [Boolean] default = True
                                ! If true, mix the passive tracers in massless layers at the bottom into the
                                ! interior as though a diffusivity of KD_MIN_TR were operating.
KD_MIN_TR = 1.0E-05             !   [m2 s-1] default = 1.0E-05
                                ! A minimal diffusivity that should always be applied to tracers, especially in
                                ! massless layers near the bottom. The default is 0.1*KD.
KD_BBL_TR = 0.0                 !   [m2 s-1] default = 0.0
                                ! A bottom boundary layer tracer diffusivity that will allow for explicitly
                                ! specified bottom fluxes. The entrainment at the bottom is at least
                                ! sqrt(Kd_BBL_tr*dt) over the same distance.
TRACER_TRIDIAG = False          !   [Boolean] default = False
                                ! If true, use the passive tracer tridiagonal solver for T and S
MINIMUM_FORCING_DEPTH = 0.001   !   [m] default = 0.001
                                ! The smallest depth over which forcing can be applied. This only takes effect
                                ! when near-surface layers become thin relative to this scale, in which case the
                                ! forcing tendencies scaled down by distributing the forcing over this depth
                                ! scale.
EVAP_CFL_LIMIT = 0.8            !   [nondim] default = 0.8
                                ! The largest fraction of a layer than can be lost to forcing (e.g. evaporation,
                                ! sea-ice formation) in one time-step. The unused mass loss is passed down
                                ! through the column.
DIAG_MLD_DENSITY_DIFF = 0.1     !   [kg/m3] default = 0.1
                                ! The density difference used to determine a diagnostic mixed layer depth,
                                ! MLD_user, following the definition of Levitus 1982. The MLD is the depth at
                                ! which the density is larger than the surface density by the specified amount.
DIAG_DEPTH_SUBML_N2 = 50.0      !   [m] default = 50.0
                                ! The distance over which to calculate a diagnostic of the stratification at the
                                ! base of the mixed layer.

! === module MOM_CVMix_KPP ===
! This is the MOM wrapper to CVMix:KPP
! See http://cvmix.github.io/
USE_KPP = False                 !   [Boolean] default = False
                                ! If true, turns on the [CVMix] KPP scheme of Large et al., 1994, to calculate
                                ! diffusivities and non-local transport in the OBL.

! === module MOM_CVMix_conv ===
! Parameterization of enhanced mixing due to convection via CVMix
USE_CVMix_CONVECTION = False    !   [Boolean] default = False
                                ! If true, turns on the enhanced mixing due to convection via CVMix. This scheme
                                ! increases diapycnal diffs./viscs. at statically unstable interfaces. Relevant
                                ! parameters are contained in the CVMix_CONVECTION% parameter block.

! === 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 = 1.095445115010332E-05
                                ! The tolerance with which to solve for entrainment values.

! === module MOM_set_diffusivity ===
FLUX_RI_MAX = 0.2               !   [nondim] default = 0.2
                                ! The flux Richardson number where the stratification is large enough that N2 >
                                ! omega2.  The full expression for the Flux Richardson number is usually
                                ! FLUX_RI_MAX*N2/(N2+OMEGA2).
SET_DIFF_ANSWER_DATE = 99991231 ! default = 99991231
                                ! The vintage of the order of arithmetic and expressions in the set diffusivity
                                ! calculations.  Values below 20190101 recover the answers from the end of 2018,
                                ! while higher values use updated and more robust forms of the same expressions.
                                ! Values above 20250301 also use less confusing expressions to set the
                                ! bottom-drag generated diffusivity when USE_LOTW_BBL_DIFFUSIVITY is false.

! === module MOM_tidal_mixing ===
! Vertical Tidal Mixing Parameterization
USE_CVMix_TIDAL = False         !   [Boolean] default = False
                                ! If true, turns on tidal mixing via CVMix
INT_TIDE_DISSIPATION = False    !   [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).
ML_RADIATION = False            !   [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.
BBL_EFFIC = 0.2                 !   [nondim] default = 0.2
                                ! The efficiency with which the energy extracted by bottom drag drives BBL
                                ! diffusion.  This is only used if BOTTOMDRAGLAW is true.
BBL_MIXING_MAX_DECAY = 200.0    !   [m] default = 200.0
                                ! The maximum decay scale for the BBL diffusion, or 0 to allow the mixing to
                                ! penetrate as far as stratification and rotation permit.  The default for now
                                ! is 200 m. This is only used if BOTTOMDRAGLAW is true.
BBL_MIXING_AS_MAX = True        !   [Boolean] default = True
                                ! If true, take the maximum of the diffusivity from the BBL mixing and the other
                                ! diffusivities. Otherwise, diffusivity from the BBL_mixing is simply added.
USE_LOTW_BBL_DIFFUSIVITY = False !   [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.
DRAG_DIFFUSIVITY_ANSWER_DATE = 20250101 ! default = 99991231
                                ! The vintage of the order of arithmetic in the drag diffusivity calculations.
                                ! Values above 20250301 use less confusing expressions to set the bottom-drag
                                ! generated diffusivity when USE_LOTW_BBL_DIFFUSIVITY is false.
DZ_BBL_AVG_MIN = 0.0            !   [m] default = 0.0
                                ! A minimal distance over which to average to determine the average bottom
                                ! boundary layer density.
SIMPLE_TKE_TO_KD = False        !   [Boolean] default = False
                                ! If true, uses a simple estimate of Kd/TKE that will work for arbitrary
                                ! vertical coordinates. If false, calculates Kd/TKE and bounds based on exact
                                ! energetics for an isopycnal layer-formulation.

! === module MOM_bkgnd_mixing ===
! Adding static vertical background mixing coefficients
KD = 1.0E-04                    !   [m2 s-1] default = 0.0
                                ! The background diapycnal diffusivity of density in the interior. Zero or the
                                ! molecular value, ~1e-7 m2 s-1, may be used.
KD_MIN = 1.0E-06                !   [m2 s-1] default = 1.0E-06
                                ! The minimum diapycnal diffusivity.
KD_ML_TOT = 1.0E-04             !   [m2 s-1] default = 1.0E-04
                                ! The total diapcynal diffusivity in the surface mixed layer when there is not a
                                ! physically based parameterization of mixing in the mixed layer, such as bulk
                                ! mixed layer or KPP or ePBL.
BRYAN_LEWIS_DIFFUSIVITY = False !   [Boolean] default = False
                                ! If true, use a Bryan & Lewis (JGR 1979) like tanh profile of background
                                ! diapycnal diffusivity with depth. This is done via CVMix.
HORIZ_VARYING_BACKGROUND = False !   [Boolean] default = False
                                ! If true, apply vertically uniform, latitude-dependent background diffusivity,
                                ! as described in Danabasoglu et al., 2012
PRANDTL_BKGND = 1.0             !   [nondim] default = 1.0
                                ! Turbulent Prandtl number used to convert vertical background diffusivities
                                ! into viscosities.
HENYEY_IGW_BACKGROUND = False   !   [Boolean] default = False
                                ! If true, use a latitude-dependent scaling for the near surface background
                                ! diffusivity, as described in Harrison & Hallberg, JPO 2008.
KD_TANH_LAT_FN = False          !   [Boolean] default = False
                                ! If true, use a tanh dependence of Kd_sfc on latitude, like CM2.1/CM2M.  There
                                ! is no physical justification for this form, and it can not be used with
                                ! HENYEY_IGW_BACKGROUND.
KD_MAX = -1.0                   !   [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.
KD_ADD = 0.0                    !   [m2 s-1] default = 0.0
                                ! A uniform diapycnal diffusivity that is added everywhere without any filtering
                                ! or scaling.
USER_CHANGE_DIFFUSIVITY = False !   [Boolean] default = False
                                ! If true, call user-defined code to change the diffusivity.
DISSIPATION_MIN = 0.0           !   [W m-3] default = 0.0
                                ! The minimum dissipation by which to determine a lower bound of Kd (a floor).
DISSIPATION_N0 = 0.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.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
DISSIPATION_KD_MIN = 0.0        !   [m2 s-1] default = 0.0
                                ! The minimum vertical diffusivity applied as a floor.
DOUBLE_DIFFUSION = False        !   [Boolean] default = False
                                ! If true, increase diffusivities for temperature or salinity based on the
                                ! double-diffusive parameterization described in Large et al. (1994).

! === module MOM_kappa_shear ===
! Parameterization of shear-driven turbulence following Jackson, Hallberg and Legg, JPO 2008
USE_JACKSON_PARAM = True        !   [Boolean] default = False
                                ! If true, use the Jackson-Hallberg-Legg (JPO 2008) shear mixing
                                ! parameterization.
VERTEX_SHEAR = False            !   [Boolean] default = False
                                ! If true, do the calculations of the shear-driven mixing at the cell vertices
                                ! (i.e., the vorticity points).
VERTEX_SHEAR_THICKNESS_MEAN = False !   [Boolean] default = False
                                ! If true, apply thickness weighting to horizontal averagings of diffusivity to
                                ! tracer points in the kappa shear solver.
RINO_CRIT = 0.25                !   [nondim] default = 0.25
                                ! The critical Richardson number for shear mixing.
SHEARMIX_RATE = 0.089           !   [nondim] default = 0.089
                                ! A nondimensional rate scale for shear-driven entrainment. Jackson et al find
                                ! values in the range of 0.085-0.089.
MAX_RINO_IT = 25                !   [nondim] default = 50
                                ! The maximum number of iterations that may be used to estimate the Richardson
                                ! number driven mixing.
KD_KAPPA_SHEAR_0 = 1.0E-04      !   [m2 s-1] default = 1.0E-04
                                ! The background diffusivity that is used to smooth the density and shear
                                ! profiles before solving for the diffusivities.  The default is the greater of
                                ! KD and 1e-7 m2 s-1.
KD_SEED_KAPPA_SHEAR = 1.0       !   [m2 s-1] default = 1.0
                                ! A moderately large seed value of diapycnal diffusivity that is used as a
                                ! starting turbulent diffusivity in the iterations to find an energetically
                                ! constrained solution for the shear-driven diffusivity.
KD_TRUNC_KAPPA_SHEAR = 1.0E-06  !   [m2 s-1] default = 1.0E-06
                                ! The value of shear-driven diffusivity that is considered negligible and is
                                ! rounded down to 0. The default is 1% of KD_KAPPA_SHEAR_0.
FRI_CURVATURE = -0.97           !   [nondim] default = -0.97
                                ! The nondimensional curvature of the function of the Richardson number in the
                                ! kappa source term in the Jackson et al. scheme.
TKE_N_DECAY_CONST = 0.24        !   [nondim] default = 0.24
                                ! The coefficient for the decay of TKE due to stratification (i.e. proportional
                                ! to N*tke). The values found by Jackson et al. are 0.24-0.28.
TKE_SHEAR_DECAY_CONST = 0.14    !   [nondim] default = 0.14
                                ! The coefficient for the decay of TKE due to shear (i.e. proportional to
                                ! |S|*tke). The values found by Jackson et al. are 0.14-0.12.
KAPPA_BUOY_SCALE_COEF = 0.82    !   [nondim] default = 0.82
                                ! The coefficient for the buoyancy length scale in the kappa equation.  The
                                ! values found by Jackson et al. are in the range of 0.81-0.86.
KAPPA_N_OVER_S_SCALE_COEF2 = 0.0 !   [nondim] default = 0.0
                                ! The square of the ratio of the coefficients of the buoyancy and shear scales
                                ! in the diffusivity equation, Set this to 0 (the default) to eliminate the
                                ! shear scale. This is only used if USE_JACKSON_PARAM is true.
LZ_RESCALE = 1.0                !   [nondim] default = 1.0
                                ! A coefficient to rescale the distance to the nearest solid boundary. This
                                ! adjustment is to account for regions where 3 dimensional turbulence prevents
                                ! the growth of shear instabilities [nondim].
KAPPA_SHEAR_TOL_ERR = 0.1       !   [nondim] default = 0.1
                                ! The fractional error in kappa that is tolerated. Iteration stops when changes
                                ! between subsequent iterations are smaller than this everywhere in a column.
                                ! The peak diffusivities usually converge most rapidly, and have much smaller
                                ! errors than this.
TKE_BACKGROUND = 0.0            !   [m2 s-2] default = 0.0
                                ! A background level of TKE used in the first iteration of the kappa equation.
                                ! TKE_BACKGROUND could be 0.
KAPPA_SHEAR_ELIM_MASSLESS = True !   [Boolean] default = True
                                ! If true, massless layers are merged with neighboring massive layers in this
                                ! calculation.  The default is true and I can think of no good reason why it
                                ! should be false. This is only used if USE_JACKSON_PARAM is true.
MAX_KAPPA_SHEAR_IT = 13         ! default = 13
                                ! The maximum number of iterations that may be used to estimate the
                                ! time-averaged diffusivity.
PRANDTL_TURB = 1.0              !   [nondim] default = 1.0
                                ! The turbulent Prandtl number applied to shear instability.
KAPPA_SHEAR_MAX_KAP_SRC_CHG = 10.0 !   [nondim] default = 10.0
                                ! The maximum permitted increase in the kappa source within an iteration
                                ! relative to the local source; this must be greater than 1.  The lower limit
                                ! for the permitted fractional decrease is (1 - 0.5/kappa_src_max_chg).  These
                                ! limits could perhaps be made dynamic with an improved iterative solver.
KAPPA_SHEAR_ITER_BUG = False    !   [Boolean] default = False
                                ! If true, use an older, dimensionally inconsistent estimate of the derivative
                                ! of diffusivity with energy in the Newton's method iteration.  The bug causes
                                ! undercorrections when dz > 1 m.
KAPPA_SHEAR_ALL_LAYER_TKE_BUG = False !   [Boolean] default = False
                                ! If true, report back the latest estimate of TKE instead of the time average
                                ! TKE when there is mass in all layers.  Otherwise always report the time
                                ! averaged TKE, as is currently done when there are some massless layers.
USE_RESTRICTIVE_TOLERANCE_CHECK = False !   [Boolean] default = False
                                ! If true, uses the more restrictive tolerance check to determine if a timestep
                                ! is acceptable for the KS_it outer iteration loop.  False uses the original
                                ! less restrictive check.

! === module MOM_CVMix_shear ===
! Parameterization of shear-driven turbulence via CVMix (various options)
USE_LMD94 = False               !   [Boolean] default = False
                                ! If true, use the Large-McWilliams-Doney (JGR 1994) shear mixing
                                ! parameterization.
USE_PP81 = False                !   [Boolean] default = False
                                ! If true, use the Pacanowski and Philander (JPO 1981) shear mixing
                                ! parameterization.

! === module MOM_CVMix_ddiff ===
! Parameterization of mixing due to double diffusion processes via CVMix
USE_CVMIX_DDIFF = False         !   [Boolean] default = False
                                ! If true, turns on double diffusive processes via CVMix. Note that double
                                ! diffusive processes on viscosity are ignored in CVMix, see
                                ! http://cvmix.github.io/ for justification.

! === module MOM_diabatic_aux ===
! The following parameters are used for auxiliary diabatic processes.
DO_BRINE_PLUME = False          !   [Boolean] default = False
                                ! If true, use a brine plume parameterization from Nguyen et al., 2009.

! === module MOM_regularize_layers ===
REGULARIZE_SURFACE_LAYERS = False !   [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_tracer_advect ===
TRACER_ADVECTION_SCHEME = "PLM" ! default = "PLM"
                                ! The horizontal transport scheme for tracers:
                                !   PLM    - Piecewise Linear Method
                                !   PPM:H3 - Piecewise Parabolic Method (Huyhn 3rd order)
                                !   PPM    - Piecewise Parabolic Method (Colella-Woodward)

! === module MOM_tracer_hor_diff ===
KHTR = 0.0                      !   [m2 s-1] default = 0.0
                                ! The background along-isopycnal tracer diffusivity.
KHTR_MIN = 0.0                  !   [m2 s-1] default = 0.0
                                ! The minimum along-isopycnal tracer diffusivity.
KHTR_MAX = 0.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 passivity is the ratio between along isopycnal mixing of tracers to
                                ! thickness mixing. A non-zero value enables this parameterization.
KHTR_PASSIVITY_MIN = 0.5        !   [nondim] default = 0.5
                                ! The minimum passivity which is the ratio between along isopycnal mixing of
                                ! tracers to thickness mixing.
DIFFUSE_ML_TO_INTERIOR = False  !   [Boolean] default = False
                                ! If true, enable epipycnal mixing between the surface boundary layer and the
                                ! interior.
CHECK_DIFFUSIVE_CFL = False     !   [Boolean] default = False
                                ! If true, use enough iterations the diffusion to ensure that the diffusive
                                ! equivalent of the CFL limit is not violated.  If false, always use the greater
                                ! of 1 or MAX_TR_DIFFUSION_CFL iteration.
MAX_TR_DIFFUSION_CFL = -1.0     !   [nondim] default = -1.0
                                ! If positive, locally limit the along-isopycnal tracer diffusivity to keep the
                                ! diffusive CFL locally at or below this value.  The number of diffusive
                                ! iterations is often this value or the next greater integer.
RECALC_NEUTRAL_SURF = False     !   [Boolean] default = False
                                ! If true, then recalculate the neutral surfaces if the
                                ! diffusive CFL is exceeded. If false, assume that the
                                ! positions of the surfaces do not change

! === module MOM_neutral_diffusion ===
! This module implements neutral diffusion of tracers
USE_NEUTRAL_DIFFUSION = False   !   [Boolean] default = False
                                ! If true, enables the neutral diffusion module.

! === module MOM_hor_bnd_diffusion ===
! This module implements horizontal diffusion of tracers near boundaries
USE_HORIZONTAL_BOUNDARY_DIFFUSION = False !   [Boolean] default = False
                                ! If true, enables the horizontal boundary tracer's diffusion module.
OBSOLETE_DIAGNOSTIC_IS_FATAL = True !   [Boolean] default = True
                                ! If an obsolete diagnostic variable appears in the diag_table, cause a FATAL
                                ! error rather than issue a WARNING.

! === module MOM_sum_output ===
CALCULATE_APE = True            !   [Boolean] default = True
                                ! If true, calculate the available potential energy of the interfaces.  Setting
                                ! this to false reduces the memory footprint of high-PE-count models
                                ! dramatically.
WRITE_STOCKS = True             !   [Boolean] default = True
                                ! If true, write the integrated tracer amounts to stdout when the energy files
                                ! are written.
MAXTRUNC = 10                   !   [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.
MAX_ENERGY = 0.0                !   [m2 s-2] default = 0.0
                                ! The maximum permitted average energy per unit mass; the model will be stopped
                                ! if there is more energy than this.  If zero or negative, this is set to
                                ! 10*MAXVEL^2.
MAXVEL = 10.0                   !   [m s-1] default = 3.0E+08
                                ! The maximum velocity allowed before the velocity components are truncated.
ENERGYFILE = "ocean.stats"      ! default = "ocean.stats"
                                ! The file to use to write the energies and globally summed diagnostics.
DATE_STAMPED_STDOUT = False     !   [Boolean] default = True
                                ! If true, use dates (not times) in messages to stdout
ISO_DATE_STAMPED_STDOUT = False !   [Boolean] default = False
                                ! If true, use ISO formatted dates in messages to stdout
TIMEUNIT = 3600.0               !   [s] default = 8.64E+04
                                ! The time unit in seconds a number of input fields
READ_DEPTH_LIST = False         !   [Boolean] default = False
                                ! Read the depth list from a file if it exists or create that file otherwise.
DEPTH_LIST_MIN_INC = 1.0E-10    !   [m] default = 1.0E-10
                                ! The minimum increment between the depths of the entries in the depth-list
                                ! file.
ENERGYSAVEDAYS = 0.5            !   [hours] default = 24.0
                                ! The interval in units of TIMEUNIT between saves of the energies of the run and
                                ! other globally summed diagnostics.
ENERGYSAVEDAYS_GEOMETRIC = 0.0  !   [hours] default = 0.0
                                ! The starting interval in units of TIMEUNIT for the first call to save the
                                ! energies of the run and other globally summed diagnostics. The interval
                                ! increases by a factor of 2. after each call to write_energy.

! === module ocean_stochastics_init ===
DO_SPPT = False                 !   [Boolean] default = False
                                ! If true, then stochastically perturb the thermodynamic tendencies of T,S, and
                                ! h.  Amplitude and correlations are controlled by the nam_stoch namelist in the
                                ! UFS model only.
DO_SKEB = False                 !   [Boolean] default = False
                                ! If true, then stochastically perturb the currents using the stochastic kinetic
                                ! energy backscatter scheme.
PERT_EPBL = False               !   [Boolean] default = False
                                ! If true, then stochastically perturb the kinetic energy production and
                                ! dissipation terms.  Amplitude and correlations are controlled by the nam_stoch
                                ! namelist in the UFS model only.

! === module MOM_surface_forcing ===
VARIABLE_WINDS = False          !   [Boolean] default = True
                                ! If true, the winds vary in time after the initialization.
VARIABLE_BUOYFORCE = False      !   [Boolean] default = True
                                ! If true, the buoyancy forcing varies in time after the initialization of the
                                ! model.
BUOY_CONFIG = "zero"            ! default = "zero"
                                ! The character string that indicates how buoyancy forcing is specified.  Valid
                                ! options include (file), (data_override), (zero), (const), (linear), (MESO),
                                ! (SCM_CVmix_tests), (BFB), (dumbbell), (USER) and (NONE).
WIND_CONFIG = "zero"            ! default = "zero"
                                ! The character string that indicates how wind forcing is specified.  Valid
                                ! options include (file), (data_override), (2gyre), (1gyre), (gyres), (zero),
                                ! (const), (Neverworld), (scurves), (ideal_hurr), (SCM_CVmix_tests) and (USER).
RESTOREBUOY = False             !   [Boolean] default = False
                                ! If true, the buoyancy fluxes drive the model back toward some specified
                                ! surface state with a rate given by FLUXCONST.
LATENT_HEAT_FUSION = 3.34E+05   !   [J/kg] default = 3.34E+05
                                ! The latent heat of fusion.
LATENT_HEAT_VAPORIZATION = 2.5E+06 !   [J/kg] default = 2.5E+06
                                ! The latent heat of fusion.
GUST_CONST = 0.0                !   [Pa] default = 0.0
                                ! The background gustiness in the winds.
USTAR_GUSTLESS_BUG = False      !   [Boolean] default = False
                                ! If true include a bug in the time-averaging of the gustless wind friction
                                ! velocity
READ_GUST_2D = False            !   [Boolean] default = False
                                ! If true, use a 2-dimensional gustiness supplied from an input file

! === module MOM_restart ===
USE_WAVES = False               !   [Boolean] default = False
                                ! If true, enables surface wave modules.

! === module MOM_main (MOM_driver) ===
DT_FORCING = 600.0              !   [s] default = 120.0
                                ! The time step for changing forcing, coupling with other components, or
                                ! potentially writing certain diagnostics. The default value is given by DT.
DAYMAX = 2.0                    !   [hours]
                                ! The final time of the whole simulation, in units of TIMEUNIT seconds.  This
                                ! also sets the potential end time of the present run segment if the end time is
                                ! not set via ocean_solo_nml in input.nml.
SINGLE_STEPPING_CALL = True     !   [Boolean] default = True
                                ! If true, advance the state of MOM with a single step including both dynamics
                                ! and thermodynamics.  If false the two phases are advanced with separate calls.
RESTART_CONTROL = 3             ! default = 1
                                ! An integer whose bits encode which restart files are written. Add 2 (bit 1)
                                ! for a time-stamped file, and odd (bit 0) for a non-time-stamped file. A
                                ! non-time-stamped restart file is saved at the end of the run segment for any
                                ! non-negative value.
RESTINT = 10.0                  !   [hours] default = 0.0
                                ! The interval between saves of the restart file in units of TIMEUNIT.  Use 0
                                ! (the default) to not save incremental restart files at all.
WRITE_CPU_STEPS = 1000          ! default = 1000
                                ! The number of coupled timesteps between writing the cpu time. If this is not
                                ! positive, do not check cpu time, and the segment run-length can not be set via
                                ! an elapsed CPU time.

! === module MOM_write_cputime ===
MAXCPU = -1.0                   !   [wall-clock seconds] default = -1.0
                                ! The maximum amount of cpu time per processor for which MOM should run before
                                ! saving a restart file and quitting with a return value that indicates that a
                                ! further run is required to complete the simulation. If automatic restarts are
                                ! not desired, use a negative value for MAXCPU.  MAXCPU has units of wall-clock
                                ! seconds, so the actual CPU time used is larger by a factor of the number of
                                ! processors used.
CPU_TIME_FILE = "CPU_stats"     ! default = "CPU_stats"
                                ! The file into which CPU time is written.

! === module MOM_file_parser ===
SEND_LOG_TO_STDOUT = False      !   [Boolean] default = False
                                ! If true, all log messages are also sent to stdout.
DOCUMENT_FILE = "MOM_parameter_doc" ! default = "MOM_parameter_doc"
                                ! The basename for files where run-time parameters, their settings, units and
                                ! defaults are documented. Blank will disable all parameter documentation.
COMPLETE_DOCUMENTATION = True   !   [Boolean] default = True
                                ! If true, all run-time parameters are documented in MOM_parameter_doc.all .
MINIMAL_DOCUMENTATION = True    !   [Boolean] default = True
                                ! If true, non-default run-time parameters are documented in
                                ! MOM_parameter_doc.short .