@@ -242,6 +242,9 @@ USE_TRIPOLAR_GEOLONB_BUG = False ! [Boolean] default = False
242242 ! along the tripolar fold to be off by 360 degrees.
243243RAD_EARTH = 6.378E+06 ! [m] default = 6.378E+06
244244 ! The radius of the Earth.
245+ READ_MEAN_SEA_LEVEL = False ! [Boolean] default = False
246+ ! If true, use a 2D map for time mean sea level, which is used to calculate time
247+ ! mean ocean total thickness.
245248TOPO_CONFIG = "file" !
246249 ! This specifies how bathymetry is specified:
247250 ! file - read bathymetric information from the file
@@ -377,6 +380,12 @@ OBC_TRACER_RESERVOIR_LENGTH_SCALE_OUT = 9000.0 ! [m] default = 0.0
377380OBC_TRACER_RESERVOIR_LENGTH_SCALE_IN = 9000.0 ! [m] default = 0.0
378381 ! An effective length scale for restoring the tracer concentration at the
379382 ! boundaries to values from the interior when the flow is entering the domain.
383+ OBC_THICKNESS_RESERVOIR_LENGTH_SCALE_OUT = 0.0 ! [m] default = 0.0
384+ ! An effective length scale for restoring the layer thickness at the boundaries
385+ ! to externally imposed values when the flow is exiting the domain.
386+ OBC_THICKNESS_RESERVOIR_LENGTH_SCALE_IN = 0.0 ! [m] default = 0.0
387+ ! An effective length scale for restoring the layer thickness at the boundaries
388+ ! to values from the interior when the flow is entering the domain.
380389OBC_REMAPPING_SCHEME = "PPM_H4" ! default = "PPM_H4"
381390 ! This sets the reconstruction scheme used for OBC vertical remapping for all
382391 ! variables. It can be one of the following schemes:
@@ -914,6 +923,9 @@ MIXING_COEFS_OBC_BUG = True ! [Boolean] default = True
914923 ! If false, use only interior data for thickness weighting in lateral mixing
915924 ! coefficient calculations and to calculate stratification and other fields at
916925 ! open boundary condition faces.
926+ RESOLN_FUNCTION_OBC_BUG = True ! [Boolean] default = True
927+ ! If false, use only interior data for calculating the resolution functions at
928+ ! open boundary condition faces and vertices.
917929KD_SMOOTH = 1.0E-06 ! [m2 s-1] default = 1.0E-06
918930 ! A diapycnal diffusivity that is used to interpolate more sensible values of T
919931 ! & S into thin layers.
@@ -937,6 +949,13 @@ DRAG_AS_BODY_FORCE = False ! [Boolean] default = False
937949CHANNEL_DRAG = True ! [Boolean] default = False
938950 ! If true, the bottom drag is exerted directly on each layer proportional to the
939951 ! fraction of the bottom it overlies.
952+ CHANNEL_DRAG_SHELFBREAK_DEPTH = 0.0 ! [m] default = 0.0
953+ ! When CHANNEL_DRAG is true, the bathymetric depth interpolated to the vorticity
954+ ! point is a combination of the harmonic mean of the adjacent velocity point
955+ ! depths below this depth and the arithmetic mean of the depths above it, to
956+ ! roughly mimic a continental shelf break profile. Setting this to exceed
957+ ! MAXIMUM_DEPTH leads to linear interpolation of the topography between velocity
958+ ! points.
940959LINEAR_DRAG = False ! [Boolean] default = False
941960 ! If LINEAR_DRAG and BOTTOMDRAGLAW are defined the drag law is
942961 ! cdrag*DRAG_BG_VEL*u.
@@ -951,6 +970,9 @@ HBBL = 10.0 ! [m]
951970CDRAG = 0.003 ! [nondim] default = 0.003
952971 ! CDRAG is the drag coefficient relating the magnitude of the velocity field to
953972 ! the bottom stress. CDRAG is only used if BOTTOMDRAGLAW is defined.
973+ CDRAG_MAP = False ! [Boolean] default = False
974+ ! If true, apply a spatially varying scaling factor to CDRAG, specified by
975+ ! CDRAG_VAR in CDRAG_FILE.
954976BBL_USE_TIDAL_BG = False ! [Boolean] default = False
955977 ! Flag to use the tidal RMS amplitude in place of constant background velocity
956978 ! for computing u* in the BBL. This flag is only used when BOTTOMDRAGLAW is true
@@ -1056,6 +1078,8 @@ TIDES = True ! [Boolean] default = False
10561078 ! If true, apply tidal momentum forcing.
10571079CALCULATE_SAL = True ! [Boolean] default = True
10581080 ! If true, calculate self-attraction and loading.
1081+ USE_HA = False ! [Boolean] default = False
1082+ ! If true, perform inline harmonic analysis.
10591083BE = 0.6 ! [nondim] default = 0.6
10601084 ! If SPLIT is true, BE determines the relative weighting of a 2nd-order
10611085 ! Runga-Kutta baroclinic time stepping scheme (0.5) and a backward Euler scheme
@@ -1136,6 +1160,9 @@ CONT_PPM_USE_VISC_REM_MAX = True ! [Boolean] default = True
11361160CONT_PPM_MARGINAL_FACE_AREAS = True ! [Boolean] default = True
11371161 ! If true, use the marginal face areas from the continuity solver for use as the
11381162 ! weights in the barotropic solver. Otherwise use the transport averaged areas.
1163+ CONT_USE_H_MARG_MIN = False ! [Boolean] default = False
1164+ ! If true, the marginal thickness used and returned from continuity is bounded
1165+ ! from below by a sub-roundoff value. Otherwise the minimum is 0.
11391166
11401167! === module MOM_CoriolisAdv ===
11411168NOSLIP = False ! [Boolean] default = False
@@ -1296,23 +1323,6 @@ TIDE_ADD_NODAL = True ! [Boolean] default = False
12961323TIDE_NODAL_REF_DATE = 1993, 7, 2 ! default = 0, 0, 0
12971324 ! Fixed reference date to use for nodal modulation of astronomical tidal
12981325 ! forcing.
1299- HA_START_TIME = 0.0 ! [days] default = 0.0
1300- ! Start time of harmonic analysis, in units of days after the start of the
1301- ! current run segment. Must be smaller than HA_END_TIME, otherwise harmonic
1302- ! analysis will not be performed. If negative, |HA_START_TIME| determines the
1303- ! length of harmonic analysis, and harmonic analysis will start |HA_START_TIME|
1304- ! days before HA_END_TIME, or at the beginning of the run segment, whichever
1305- ! occurs later.
1306- HA_END_TIME = 0.0 ! [days] default = 0.0
1307- ! End time of harmonic analysis, in units of days after the start of the current
1308- ! run segment. Must be positive and smaller than the length of the currnet run
1309- ! segment, otherwise harmonic analysis will not be performed.
1310- HA_SSH = False ! [Boolean] default = False
1311- ! If true, perform harmonic analysis of sea serface height.
1312- HA_UBT = False ! [Boolean] default = False
1313- ! If true, perform harmonic analysis of zonal barotropic velocity.
1314- HA_VBT = False ! [Boolean] default = False
1315- ! If true, perform harmonic analysis of meridional barotropic velocity.
13161326
13171327! === module MOM_PressureForce ===
13181328ANALYTIC_FV_PGF = True ! [Boolean] default = True
@@ -1381,6 +1391,7 @@ PRESSURE_RECONSTRUCTION_SCHEME = 1 ! default = 1
13811391 ! 0: PCM or no reconstruction.
13821392 ! 1: PLM reconstruction.
13831393 ! 2: PPM reconstruction.
1394+ ! 3: PLM with least squares slope.
13841395BOUNDARY_EXTRAPOLATION_PRESSURE = True ! [Boolean] default = True
13851396 ! If true, the reconstruction of T & S for pressure in boundary cells is
13861397 ! extrapolated, rather than using PCM in these cells. If true, the same order
0 commit comments