! Where appropriate, parameters use usually given in MKS units. ! This particular file is for the example in ice_ocean_SIS2/Baltic_OM4_025. ! This MOM_input file typically contains only the non-default values that are needed to reproduce this example. ! A full list of parameters for this example can be found in the corresponding MOM_parameter_doc.all file ! which is generated by the model at run-time. ! === module MOM_domains === REENTRANT_X = False ! [Boolean] default = True ! If true, the domain is zonally reentrant. !SYMMETRIC_MEMORY_ = True !STATIC_MEMORY_ = False NIGLOBAL = @[NIGLOBAL] ! ! 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 = @[NJGLOBAL] ! ! 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. With ! STATIC_MEMORY_ this is set as NIHALO_ in MOM_memory.h at compile time; without ! STATIC_MEMORY_ the default is NIHALO_ in MOM_memory.h (if defined) or 2. !NJHALO = 4 ! default = 4 ! The number of halo points on each side in the y-direction. With ! STATIC_MEMORY_ this is set as NJHALO_ in MOM_memory.h at compile time; without ! STATIC_MEMORY_ the default is NJHALO_ in MOM_memory.h (if defined) or 2. ! LAYOUT = 4, 4 ! ! The processor layout that was actually used. ! IO_LAYOUT = 1, 1 ! default = 1 ! The processor layout to be used, or 0,0 to automatically set the io_layout to ! be the same as the layout. ! === module MOM === 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_REGRIDDING = True ! [Boolean] default = False ! If True, use the ALE algorithm (regridding/remapping). If False, use the ! layered isopycnal algorithm. THICKNESSDIFFUSE = True ! [Boolean] default = False ! If true, interface heights are diffused with a coefficient of KHTH. THICKNESSDIFFUSE_FIRST = True ! [Boolean] default = False ! If true, do thickness diffusion before dynamics. This is only used if ! THICKNESSDIFFUSE is true. DT = 360.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 = 360.0 ! [s] default = 3600.0 ! The thermodynamic and tracer advection time step. Ideally DT_THERM should be ! an integer multiple of DT and less than the forcing or coupling time-step, ! unless THERMO_SPANS_COUPLING is true, in which case DT_THERM can be an integer ! multiple of the coupling timestep. By default DT_THERM is set to DT. THERMO_SPANS_COUPLING = True ! [Boolean] default = False ! If true, the MOM will take thermodynamic and tracer 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. USE_PSURF_IN_EOS = False ! [Boolean] default = False ! If true, always include the surface pressure contributions in equation of ! state calculations. FRAZIL = True ! [Boolean] default = False ! If true, water freezes if it gets too cold, and the accumulated heat deficit ! is returned in the surface state. FRAZIL is only used if ! ENABLE_THERMODYNAMICS is true. DO_GEOTHERMAL = True ! [Boolean] default = False ! If true, apply geothermal heating. BOUND_SALINITY = True ! [Boolean] default = False ! If true, limit salinity to being positive. (The sea-ice model may ask for more ! salt than is available and drive the salinity negative otherwise.) MIN_SALINITY = 0.01 ! [PPT] default = 0.01 ! The minimum value of salinity when BOUND_SALINITY=True. The default is 0.01 ! for backward compatibility but ideally should be 0. C_P = 3925.0 !3992.0 ! [J kg-1 K-1] default = 3991.86795711963 ! The heat capacity of sea water, approximated as a constant. This is only used ! if ENABLE_THERMODYNAMICS is true. The default value is from the TEOS-10 ! definition of conservative temperature. CHECK_BAD_SURFACE_VALS = True ! [Boolean] default = False ! If true, check the surface state for ridiculous values. BAD_VAL_SSH_MAX = 50.0 ! [m] default = 20.0 ! The value of SSH above which a bad value message is triggered, if ! CHECK_BAD_SURFACE_VALS is true. BAD_VAL_SSS_MAX = 75.0 ! [PPT] default = 45.0 ! The value of SSS above which a bad value message is triggered, if ! CHECK_BAD_SURFACE_VALS is true. BAD_VAL_SST_MAX = 55.0 ! [deg C] default = 45.0 ! The value of SST above which a bad value message is triggered, if ! CHECK_BAD_SURFACE_VALS is true. BAD_VAL_SST_MIN = -3.0 ! [deg C] default = -2.1 ! The value of SST below which a bad value message is triggered, if ! CHECK_BAD_SURFACE_VALS is true. DEFAULT_2018_ANSWERS = False !True ! [Boolean] default = True ! This sets the default value for the various _2018_ANSWERS parameters. WRITE_GEOM = 0 ! 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. SAVE_INITIAL_CONDS = False ! [Boolean] default = False ! If true, write the initial conditions to a file given by IC_OUTPUT_FILE. IC_OUTPUT_FILE = "MOM_IC" ! default = "MOM_IC" ! The file into which to write the initial conditions. ! === module MOM_hor_index === ! Sets the horizontal array index types. ! === module MOM_fixed_initialization === INPUTDIR = "INPUT" ! default = "." ! The directory in which input files are found. ! === module MOM_grid_init === GRID_CONFIG = "mosaic" ! ! A character string that determines the method for defining the horizontal ! grid. Current options are: ! mosaic - read the grid from a mosaic (supergrid) ! file set by GRID_FILE. ! cartesian - use a (flat) Cartesian grid. ! spherical - use a simple spherical grid. ! mercator - use a Mercator spherical grid. GRID_FILE = "ocean_hgrid.nc" ! ! Name of the file from which to read horizontal grid data. TOPO_CONFIG = "file" ! ! This specifies how bathymetry is specified: ! file - read bathymetric information from the file ! specified by (TOPO_FILE). ! flat - flat bottom set to MAXIMUM_DEPTH. ! bowl - an analytically specified bowl-shaped basin ! ranging between MAXIMUM_DEPTH and MINIMUM_DEPTH. ! spoon - a similar shape to 'bowl', but with an vertical ! wall at the southern face. ! halfpipe - a zonally uniform channel with a half-sine ! profile in the meridional direction. ! benchmark - use the benchmark test case topography. ! Neverland - use the Neverland 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. TOPO_FILE = "ocean_topog.nc" ! default = "topog.nc" ! The file from which the bathymetry is read. MAXIMUM_DEPTH = 6500.0 ! [m] ! The maximum depth of the ocean. MINIMUM_DEPTH = 5.0 !9.5 ! [m] default = 0.0 ! If MASKING_DEPTH is unspecified, then anything shallower than MINIMUM_DEPTH is ! assumed to be land and all fluxes are masked out. If MASKING_DEPTH is ! specified, then all depths shallower than MINIMUM_DEPTH but deeper than ! MASKING_DEPTH are rounded to MINIMUM_DEPTH. ! === module MOM_open_boundary === ! Controls where open boundaries are located, what kind of boundary condition to impose, and what data to apply, ! if any. ! khs_testMASKING_DEPTH = 1.0 ! [m] default = -9999.0 MASKING_DEPTH = 5.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 negative. 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. !CHANNEL_LIST_FILE = "MOM_channels_global_025" ! default = "MOM_channel_list" ! The file from which the list of narrowed channels is read. OBC_NUMBER_OF_SEGMENTS = 4 ! default = 0 ! The number of open boundary segments. OBC_FREESLIP_VORTICITY = True ! [Boolean] default = True ! If true, sets the normal gradient of tangential velocity to zero in the ! relative vorticity on open boundaries. This cannot be true if another ! OBC_XXX_VORTICITY option is True. OBC_COMPUTED_VORTICITY = False ! [Boolean] default = False ! If true, uses the external values of tangential velocity in the relative ! vorticity on open boundaries. This cannot be true if another OBC_XXX_VORTICITY ! option is True. OBC_ZERO_BIHARMONIC = False ! [Boolean] default = False ! If true, zeros the Laplacian of flow on open boundaries in the biharmonic ! viscosity term. OBC_FREESLIP_STRAIN = False ! [Boolean] default = True OBC_COMPUTED_STRAIN = True ! [Boolean] default = False OBC_SEGMENT_001 = "J=N,I=N:0,FLATHER,ORLANSKI,NUDGED,ORLANSKI_TAN,NUDGED_TAN" ! north OBC_SEGMENT_002 = "J=0,I=0:N,FLATHER,ORLANSKI,NUDGED,ORLANSKI_TAN,NUDGED_TAN" ! south OBC_SEGMENT_003 = "I=N,J=0:N,FLATHER,ORLANSKI,NUDGED,ORLANSKI_TAN,NUDGED_TAN" ! east OBC_SEGMENT_004 = "I=0,J=N:0,FLATHER,ORLANSKI,NUDGED,ORLANSKI_TAN,NUDGED_TAN" ! west ! Documentation needs to be dynamic????? OBC_SEGMENT_001_VELOCITY_NUDGING_TIMESCALES = 3, 360.0 !3, 360.0 !0.3, 360.0 !3, 3600.0 !0.3, 360.0 ! [days] default = 0.0 ! Timescales in days for nudging along a segment, for inflow, then outflow. ! Setting both to zero should behave like SIMPLE obcs for the baroclinic ! velocities. ! Documentation needs to be dynamic????? OBC_SEGMENT_002_VELOCITY_NUDGING_TIMESCALES = 3, 360.0 !3, 360.0 !0.3, 360.0 !3, 3600.0 !0.3, 360.0 ! [days] default = 0.0 ! Timescales in days for nudging along a segment, for inflow, then outflow. ! Setting both to zero should behave like SIMPLE obcs for the baroclinic ! velocities. OBC_SEGMENT_003_VELOCITY_NUDGING_TIMESCALES = 3, 360.0 !3, 360.0 !0.3, 360.0 !3, 3600.0 !0.3, 360.0 !3, 3600.0 ! [days] default = 0.0 ! Timescales in days for nudging along a segment, for inflow, then outflow. ! Setting both to zero should behave like SIMPLE obcs for the baroclinic ! velocities. ! Documentation needs to be dynamic????? OBC_SEGMENT_004_VELOCITY_NUDGING_TIMESCALES = 3, 360.0 !3, 360.0 !0.3, 360.0 !3, 3600.0 !0.3, 360.0 !3, 3600.0 ! [days] default = 0.0 ! Timescales in days for nudging along a segment, for inflow, then outflow. ! Setting both to zero should behave like SIMPLE obcs for the baroclinic ! velocities. ! Documentation needs to be dynamic????? BRUSHCUTTER_MODE = True ! [Boolean] default = False ! If true, read external OBC data on the supergrid. OBC_SEGMENT_001_DATA = "U=file:ocean_uv_obc_north.nc(u),V=file:ocean_uv_obc_north.nc(v),SSH=file:ocean_ssh_obc_north.nc(ssh),TEMP=file:ocean_ts_obc_north.nc(temp),SALT=file:ocean_ts_obc_north.nc(salt)" ! OBC_SEGMENT_002_DATA = "U=file:ocean_uv_obc_south.nc(u),V=file:ocean_uv_obc_south.nc(v),SSH=file:ocean_ssh_obc_south.nc(ssh),TEMP=file:ocean_ts_obc_south.nc(temp),SALT=file:ocean_ts_obc_south.nc(salt)" ! OBC_SEGMENT_003_DATA = "U=file:ocean_uv_obc_east.nc(u),V=file:ocean_uv_obc_east.nc(v),SSH=file:ocean_ssh_obc_east.nc(ssh),TEMP=file:ocean_ts_obc_east.nc(temp),SALT=file:ocean_ts_obc_east.nc(salt)" ! OBC_SEGMENT_004_DATA = "U=file:ocean_uv_obc_west.nc(u),V=file:ocean_uv_obc_west.nc(v),SSH=file:ocean_ssh_obc_west.nc(ssh),TEMP=file:ocean_ts_obc_west.nc(temp),SALT=file:ocean_ts_obc_west.nc(salt)" ! OBC_TRACER_RESERVOIR_LENGTH_SCALE_OUT = 3000.0 !0.0 !3000.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 = 3.0E+04 !0.0 !3.0E+04 ! [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. RAMP_OBCS = False ! default = False OBC_RAMP_TIMESCALE = 2.0 ! [days] default = 1.0 OBC_RADIATION_MAX = 1.0 ! [nondim] default = 1.0 ! === module MOM_verticalGrid === ! Parameters providing information about the vertical grid. NK = 55 ! [nondim] ! The number of model layers. ! === module MOM_tracer_registry === ! === module MOM_EOS === DTFREEZE_DP = -7.75E-08 ! [deg C Pa-1] default = 0.0 ! When TFREEZE_FORM=LINEAR, this is the derivative of the freezing potential ! temperature with pressure. ! === module MOM_restart === PARALLEL_RESTARTFILES = False ! [Boolean] default = False ! If true, each processor writes its own restart file, otherwise a single ! restart file is generated RESTART_CHECKSUMS_REQUIRED = False ! === module MOM_tracer_flow_control === USE_IDEAL_AGE_TRACER = True ! [Boolean] default = False ! If true, use the ideal_age_example tracer package. ! === module ideal_age_example === ! === module MOM_coord_initialization === COORD_CONFIG = "gprime" !"file" ! 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. GINT = 0.0098 ! [m s-2] !COORD_FILE = "layer_coord.nc" ! ! The file from which the coordinate densities are read. REMAP_UV_USING_OLD_ALG = True ! [Boolean] default = True ! If true, uses the old remapping-via-a-delta-z method for remapping u and v. If ! false, uses the new method that remaps between grids described by an old and ! new thickness. REGRIDDING_COORDINATE_MODE = "Z*" !"HYCOM1" ! default = "LAYER" ! Coordinate mode for vertical regridding. Choose among the following ! possibilities: LAYER - Isopycnal or stacked shallow water layers ! ZSTAR, Z* - stretched geopotential z* ! SIGMA_SHELF_ZSTAR - stretched geopotential z* ignoring shelf ! SIGMA - terrain following coordinates ! RHO - continuous isopycnal ! HYCOM1 - HyCOM-like hybrid coordinate ! SLIGHT - stretched coordinates above continuous isopycnal ! ADAPTIVE - optimize for smooth neutral density surfaces BOUNDARY_EXTRAPOLATION = True ! [Boolean] default = False ! When defined, a proper high-order reconstruction scheme is used within ! boundary cells rather than PCM. E.g., if PPM is used for remapping, a PPM ! reconstruction will also be used within boundary cells. ALE_COORDINATE_CONFIG = "FILE:ocean_vgrid.nc,dz" !"HYBRID:hycom1_75_800m.nc,sigma2,FNC1:2,4000,4.5,.01" ! default = "UNIFORM" ! Determines how to specify the coordinate resolution. Valid options are: ! PARAM - use the vector-parameter ALE_RESOLUTION ! UNIFORM[:N] - uniformly distributed ! FILE:string - read from a file. The string specifies ! the filename and variable name, separated ! by a comma or space, e.g. FILE:lev.nc,dz ! or FILE:lev.nc,interfaces=zw ! WOA09[:N] - the WOA09 vertical grid (approximately) ! FNC1:string - FNC1:dz_min,H_total,power,precision ! HYBRID:string - read from a file. The string specifies ! the filename and two variable names, separated ! by a comma or space, for sigma-2 and dz. e.g. ! HYBRID:vgrid.nc,sigma2,dz !ALE_RESOLUTION = 7*2.0, 2*2.01, 2.02, 2.03, 2.05, 2.08, 2.11, 2.15, 2.21, 2.2800000000000002, 2.37, 2.48, 2.61, 2.77, 2.95, 3.17, 3.4299999999999997, 3.74, 4.09, 4.49, 4.95, 5.48, 6.07, 6.74, 7.5, 8.34, 9.280000000000001, 10.33, 11.49, 12.77, 14.19, 15.74, 17.450000000000003, 19.31, 21.35, 23.56, 25.97, 28.580000000000002, 31.41, 34.47, 37.77, 41.32, 45.14, 49.25, 53.65, 58.370000000000005, 63.42, 68.81, 74.56, 80.68, 87.21000000000001, 94.14, 101.51, 109.33, 117.62, 126.4, 135.68, 145.5, 155.87, 166.81, 178.35, 190.51, 203.31, 216.78, 230.93, 245.8, 261.42, 277.83 ! [m] ! The distribution of vertical resolution for the target ! grid used for Eulerian-like coordinates. For example, ! in z-coordinate mode, the parameter is a list of level ! thicknesses (in m). In sigma-coordinate mode, the list ! is of non-dimensional fractions of the water column. !TARGET_DENSITIES = 1010.0, 1014.3034, 1017.8088, 1020.843, 1023.5566, 1025.813, 1027.0275, 1027.9114, 1028.6422, 1029.2795, 1029.852, 1030.3762, 1030.8626, 1031.3183, 1031.7486, 1032.1572, 1032.5471, 1032.9207, 1033.2798, 1033.6261, 1033.9608, 1034.2519, 1034.4817, 1034.6774, 1034.8508, 1035.0082, 1035.1533, 1035.2886, 1035.4159, 1035.5364, 1035.6511, 1035.7608, 1035.8661, 1035.9675, 1036.0645, 1036.1554, 1036.2411, 1036.3223, 1036.3998, 1036.4739, 1036.5451, 1036.6137, 1036.68, 1036.7441, 1036.8062, 1036.8526, 1036.8874, 1036.9164, 1036.9418, 1036.9647, 1036.9857, 1037.0052, 1037.0236, 1037.0409, 1037.0574, 1037.0738, 1037.0902, 1037.1066, 1037.123, 1037.1394, 1037.1558, 1037.1722, 1037.1887, 1037.206, 1037.2241, 1037.2435, 1037.2642, 1037.2866, 1037.3112, 1037.3389, 1037.3713, 1037.4118, 1037.475, 1037.6332, 1037.8104, 1038.0 ! [m] ! HYBRID target densities for interfaces REGRID_COMPRESSIBILITY_FRACTION = 0.01 ! [nondim] default = 0.0 ! When interpolating potential density profiles we can add some artificial ! compressibility solely to make homogeneous regions appear stratified. MAXIMUM_INT_DEPTH_CONFIG = "FNC1:5,8000.0,1.0,.01" ! default = "NONE" ! Determines how to specify the maximum interface depths. ! Valid options are: ! NONE - there are no maximum interface depths ! PARAM - use the vector-parameter MAXIMUM_INTERFACE_DEPTHS ! 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,Z ! FNC1:string - FNC1:dz_min,H_total,power,precision !MAXIMUM_INT_DEPTHS = 0.0, 5.0, 12.75, 23.25, 36.49, 52.480000000000004, 71.22, 92.71000000000001, 116.94000000000001, 143.92000000000002, 173.65, 206.13, 241.36, 279.33000000000004, 320.05000000000007, 363.5200000000001, 409.7400000000001, 458.7000000000001, 510.4100000000001, 564.8700000000001, 622.0800000000002, 682.0300000000002, 744.7300000000002, 810.1800000000003, 878.3800000000003, 949.3300000000004, 1023.0200000000004, 1099.4600000000005, 1178.6500000000005, 1260.5900000000006, 1345.2700000000007, 1432.7000000000007, 1522.8800000000008, 1615.8100000000009, 1711.490000000001, 1809.910000000001, 1911.080000000001, 2015.0000000000011, 2121.670000000001, 2231.080000000001, 2343.2400000000007, 2458.1500000000005, 2575.8100000000004, 2696.2200000000003, 2819.3700000000003, 2945.2700000000004, 3073.9200000000005, 3205.3200000000006, 3339.4600000000005, 3476.3500000000004, 3615.9900000000002, 3758.38, 3903.52, 4051.4, 4202.03, 4355.41, 4511.54, 4670.41, 4832.03, 4996.4, 5163.5199999999995, 5333.379999999999, 5505.989999999999, 5681.3499999999985, 5859.459999999998, 6040.319999999998, 6223.919999999998, 6410.269999999999, 6599.369999999999, 6791.219999999999, 6985.8099999999995, 7183.15, 7383.24, 7586.08, 7791.67, 8000.0 ! The list of maximum depths for each interface. MAX_LAYER_THICKNESS_CONFIG = "FNC1:400,31000.0,0.1,.01" ! default = "NONE" ! Determines how to specify the maximum layer thicknesses. ! Valid options are: ! NONE - there are no maximum layer thicknesses ! PARAM - use the vector-parameter MAX_LAYER_THICKNESS ! 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,Z ! FNC1:string - FNC1:dz_min,H_total,power,precision !MAX_LAYER_THICKNESS = 400.0, 409.63, 410.32, 410.75, 411.07, 411.32, 411.52, 411.7, 411.86, 412.0, 412.13, 412.24, 412.35, 412.45, 412.54, 412.63, 412.71, 412.79, 412.86, 412.93, 413.0, 413.06, 413.12, 413.18, 413.24, 413.29, 413.34, 413.39, 413.44, 413.49, 413.54, 413.58, 413.62, 413.67, 413.71, 413.75, 413.78, 413.82, 413.86, 413.9, 413.93, 413.97, 414.0, 414.03, 414.06, 414.1, 414.13, 414.16, 414.19, 414.22, 414.24, 414.27, 414.3, 414.33, 414.35, 414.38, 414.41, 414.43, 414.46, 414.48, 414.51, 414.53, 414.55, 414.58, 414.6, 414.62, 414.65, 414.67, 414.69, 414.71, 414.73, 414.75, 414.77, 414.79, 414.83 ! [m] ! The list of maximum thickness for each layer. REMAPPING_SCHEME = "PPM_H4" ! default = "PLM" ! This sets the reconstruction scheme used for vertical remapping for all ! variables. It can be one of the following schemes: PCM (1st-order ! accurate) ! PLM (2nd-order accurate) ! PPM_H4 (3rd-order accurate) ! PPM_IH4 (3rd-order accurate) ! PQM_IH4IH3 (4th-order accurate) ! PQM_IH6IH5 (5th-order accurate) ! === module MOM_grid === ! Parameters providing information about the lateral grid. ! === module MOM_state_initialization === TS_CONFIG = "file" ! ! A string that determines how the initial tempertures and salinities are ! specified for a new run: ! file - read velocities from the file specified ! by (TS_FILE). ! fit - find the temperatures that are consistent with ! the layer densities and salinity S_REF. ! TS_profile - use temperature and salinity profiles ! (read from TS_FILE) to set layer densities. ! benchmark - use the benchmark test case T & S. ! linear - linear in logical layer space. ! DOME2D - 2D DOME initialization. ! ISOMIP - ISOMIP initialization. ! adjustment2d - 2d lock exchange T/S ICs. ! sloshing - sloshing mode T/S ICs. ! seamount - no motion test with seamount ICs. ! dumbbell - sloshing channel ICs. ! rossby_front - a mixed layer front in thermal wind balance. ! SCM_CVMix_tests - used in the SCM CVMix tests. ! USER - call a user modified routine. TS_FILE = "ocean_ts_ic.nc" ! ! The initial condition file for temperature. INIT_LAYERS_FROM_Z_FILE = True ! [Boolean] default = False ! If true, initialize the layer thicknesses, temperatures, and salinities from a ! Z-space file on a latitude-longitude grid. ! === module MOM_initialize_layers_from_Z === TEMP_SALT_Z_INIT_FILE = "ocean_ts_ic.nc" ! default = "temp_salt_z.nc" ! The name of the z-space input file used to initialize temperatures (T) and ! salinities (S). If T and S are not in the same file, TEMP_Z_INIT_FILE and ! SALT_Z_INIT_FILE must be set. !TEMP_Z_INIT_FILE = "" ! default = "" ! The name of the z-space input file used to initialize temperatures, only. !SALT_Z_INIT_FILE = "" ! default = "" ! The name of the z-space input file used to initialize temperatures, only. Z_INIT_FILE_PTEMP_VAR = "Temp" ! default = "ptemp" ! The name of the potential temperature variable in TEMP_Z_INIT_FILE. Z_INIT_FILE_SALT_VAR = "Salt" ! default = "salt" ! The name of the salinity variable in SALT_Z_INIT_FILE. Z_INIT_ALE_REMAPPING = True ! [Boolean] default = False ! If True, then remap straight to model coordinate from file. Z_INIT_REMAP_OLD_ALG = True ! [Boolean] default = True ! If false, uses the preferred remapping algorithm for initialization. If true, ! use an older, less robust algorithm for remapping. Z_INIT_REMAP_GENERAL = True ! [Boolean] default = False ! If false, only initializes to z* coordinates. ! If true, allows initialization directly to general coordinates. Z_INIT_REMAP_FULL_COLUMN = True ! [Boolean] default = True ! If false, only reconstructs profiles for valid data points. ! If true, inserts vanished layers below the valid data. DEPRESS_INITIAL_SURFACE = True ! [Boolean] default = False ! If true, depress the initial surface to avoid huge tsunamis when a large ! surface pressure is applied. SURFACE_HEIGHT_IC_FILE = "ocean_ssh_ic.nc" ! ! The initial condition file for the surface height. SURFACE_HEIGHT_IC_VAR = "ave_ssh" ! default = "SSH" ! The initial condition variable for the surface height VELOCITY_CONFIG = "file" VELOCITY_FILE = "ocean_uv_ic.nc" ! === module MOM_diag_mediator === NUM_DIAG_COORDS = 2 ! default = 1 ! The number of diagnostic vertical coordinates to use. For each coordinate, an ! entry in DIAG_COORDS must be provided. DIAG_COORDS = "z Z ZSTAR", "rho2 RHO2 RHO" ! ! 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_COORD_DEF_RHO2 = "RFNC1:35,999.5,1028,1028.5,8.,1038.,0.0078125" ! default = "WOA09" ! Determines how to specify the coordinate resolution. Valid options are: ! PARAM - use the vector-parameter DIAG_COORD_RES_RHO2 ! UNIFORM[:N] - uniformly distributed ! FILE:string - read from a file. The string specifies ! the filename and variable name, separated ! by a comma or space, e.g. FILE:lev.nc,dz ! or FILE:lev.nc,interfaces=zw ! WOA09[:N] - the WOA09 vertical grid (approximately) ! FNC1:string - FNC1:dz_min,H_total,power,precision ! HYBRID:string - read from a file. The string specifies ! the filename and two variable names, separated ! by a comma or space, for sigma-2 and dz. e.g. ! HYBRID:vgrid.nc,sigma2,dz DIAG_COORD_DEF_Z = "FILE:ocean_vgrid.nc,interfaces=zw" ! default = "WOA09" ! Determines how to specify the coordinate ! resolution. Valid options are: ! PARAM - use the vector-parameter DIAG_COORD_RES_RHO2 ! UNIFORM[:N] - uniformly distributed ! FILE:string - read from a file. The string specifies ! the filename and variable name, separated ! by a comma or space, e.g. FILE:lev.nc,dz ! or FILE:lev.nc,interfaces=zw ! WOA09[:N] - the WOA09 vertical grid (approximately) ! FNC1:string - FNC1:dz_min,H_total,power,precision ! HYBRID:string - read from a file. The string specifies ! the filename and two variable names, separated ! by a comma or space, for sigma-2 and dz. e.g. ! HYBRID:vgrid.nc,sigma2,dz ! === module MOM_MEKE === USE_MEKE = True ! [Boolean] default = False ! If true, turns on the MEKE scheme which calculates a sub-grid mesoscale eddy ! kinetic energy budget. MEKE_GMCOEFF = 1.0 ! [nondim] default = -1.0 ! The efficiency of the conversion of potential energy into MEKE by the ! thickness mixing parameterization. If MEKE_GMCOEFF is negative, this ! conversion is not used or calculated. MEKE_BGSRC = 1.0E-13 ! [W kg-1] default = 0.0 ! A background energy source for MEKE. MEKE_KHMEKE_FAC = 1.0 ! [nondim] default = 0.0 ! A factor that maps MEKE%Kh to Kh for MEKE itself. MEKE_ALPHA_RHINES = 0.15 ! [nondim] default = 0.05 ! If positive, is a coefficient weighting the Rhines scale in the expression for ! mixing length used in MEKE-derived diffusivity. MEKE_ALPHA_EADY = 0.15 ! [nondim] default = 0.05 ! If positive, is a coefficient weighting the Eady length scale in the ! expression for mixing length used in MEKE-derived diffusivity. ! === module MOM_lateral_mixing_coeffs === USE_VARIABLE_MIXING = True ! [Boolean] default = False ! If true, the variable mixing code will be called. This allows diagnostics to ! be created even if the scheme is not used. If KHTR_SLOPE_CFF>0 or ! KhTh_Slope_Cff>0, this is set to true regardless of what is in the parameter ! file. RESOLN_SCALED_KH = True ! [Boolean] default = False ! If true, the Laplacian lateral viscosity is scaled away when the first ! baroclinic deformation radius is well resolved. RESOLN_SCALED_KHTH = True ! [Boolean] default = False ! If true, the interface depth diffusivity is scaled away when the first ! baroclinic deformation radius is well resolved. KHTR_SLOPE_CFF = 0.25 ! [nondim] default = 0.0 ! The nondimensional coefficient in the Visbeck formula for the epipycnal tracer ! diffusivity USE_STORED_SLOPES = True ! [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. INTERPOLATE_RES_FN = False ! [Boolean] default = True ! If true, interpolate the resolution function to the velocity points from the ! thickness points; otherwise interpolate the wave speed and calculate the ! resolution function independently at each point. GILL_EQUATORIAL_LD = True ! [Boolean] default = False ! If true, uses Gill's definition of the baroclinic equatorial deformation ! radius, otherwise, if false, use Pedlosky's definition. These definitions ! differ by a factor of 2 in front of the beta term in the denominator. Gill's ! is the more appropriate definition. ! === module MOM_set_visc === CHANNEL_DRAG = True ! [Boolean] default = False ! If true, the bottom drag is exerted directly on each layer proportional to the ! fraction of the bottom it overlies. PRANDTL_TURB = 1.25 ! [nondim] default = 1.0 ! The turbulent Prandtl number applied to shear instability. HBBL = 10.0 ! [m] ! The thickness of a bottom boundary layer with a viscosity of KVBBL if ! BOTTOMDRAGLAW is not defined, or the thickness over which near-bottom ! velocities are averaged for the drag law if BOTTOMDRAGLAW is defined but ! LINEAR_DRAG is not. DRAG_BG_VEL = 0.01 ! [m s-1] default = 0.0 ! DRAG_BG_VEL is either the assumed bottom velocity (with LINEAR_DRAG) or an ! unresolved velocity that is combined with the resolved velocity to estimate ! the velocity magnitude. DRAG_BG_VEL is only used when BOTTOMDRAGLAW is ! defined. BBL_USE_EOS = True ! [Boolean] default = False ! If true, use the equation of state in determining the properties of the bottom ! boundary layer. Otherwise use the layer target potential densities. BBL_THICK_MIN = 0.1 ! [m] default = 0.0 ! The minimum bottom boundary layer thickness that can be used with ! BOTTOMDRAGLAW. This might be Kv/(cdrag*drag_bg_vel) to give Kv as the minimum ! near-bottom viscosity. KV = 1.0E-06 ! [m2 s-1] ! The background kinematic viscosity in the interior. The molecular value, ~1e-6 ! m2 s-1, may be used. KV_BBL_MIN = 0.0 ! [m2 s-1] default = 1.0E-04 ! The minimum viscosities in the bottom boundary layer. KV_TBL_MIN = 0.0 ! [m2 s-1] default = 1.0E-04 ! The minimum viscosities in the top boundary layer. ! === module MOM_thickness_diffuse === KHTH_MAX_CFL = 0.1 ! [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. USE_GM_WORK_BUG = True ! [Boolean] default = True ! If true, compute the top-layer work tendency on the u-grid with the incorrect ! sign, for legacy reproducibility. ! === module MOM_continuity === ! === module MOM_continuity_PPM === ETA_TOLERANCE = 1.0E-06 ! [m] default = 3.75E-09 ! The tolerance for the differences between the barotropic and baroclinic ! estimates of the sea surface height due to the fluxes through each face. The ! total tolerance for SSH is 4 times this value. The default is ! 0.5*NK*ANGSTROM, and this should not be set less than about ! 10^-15*MAXIMUM_DEPTH. ETA_TOLERANCE_AUX = 0.001 ! [m] default = 1.0E-06 ! The tolerance for free-surface height discrepancies between the barotropic ! solution and the sum of the layer thicknesses when calculating the auxiliary ! corrected velocities. By default, this is the same as ETA_TOLERANCE, but can ! be made larger for efficiency. ! === module MOM_CoriolisAdv === CORIOLIS_SCHEME = "SADOURNY75_ENSTRO" ! 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 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. ! === module MOM_PressureForce === ! === module MOM_PressureForce_AFV === MASS_WEIGHT_IN_PRESSURE_GRADIENT = True ! [Boolean] default = False ! If true, use mass weighting when interpolating T/S for integrals near the ! bathymetry in AFV pressure gradient calculations. ! === module MOM_hor_visc === LAPLACIAN = True ! [Boolean] default = False ! If true, use a Laplacian horizontal viscosity. AH_VEL_SCALE = 0.005 ! [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. SMAGORINSKY_AH = False ! [Boolean] default = False ! If true, use a biharmonic Smagorinsky nonlinear eddy viscosity. SMAG_BI_CONST = 0.06 ! [nondim] default = 0.0 ! The nondimensional biharmonic Smagorinsky constant, typically 0.015 - 0.06. USE_LAND_MASK_FOR_HVISC = False ! [Boolean] default = False ! If true, use 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. Default is False in order to maintain answers with ! legacy experiments but should be changed to True for new experiments. ! === module MOM_vert_friction === HMIX_FIXED = 0.5 ! [m] ! The prescribed depth over which the near-surface viscosity and diffusivity are ! elevated when the bulk mixed layer is not used. ! khs_testMAXVEL = 6.0 ! [m s-1] default = 3.0E+08 MAXVEL = 3.0E+08 ! [m s-1] default = 3.0E+08 ! The maximum velocity allowed before the velocity components are truncated. ! === module MOM_PointAccel === U_TRUNC_FILE = "U_velocity_truncations" ! default = "" ! The absolute path to a file into which the accelerations leading to zonal ! velocity truncations are written. Undefine this for efficiency if this ! diagnostic is not needed. V_TRUNC_FILE = "V_velocity_truncations" ! default = "" ! The absolute path to a file into which the accelerations leading to meridional ! velocity truncations are written. Undefine this for efficiency if this ! diagnostic is not needed. ! === module MOM_barotropic === BOUND_BT_CORRECTION = True ! [Boolean] default = False ! If true, the corrective pseudo mass-fluxes into the barotropic solver are ! limited to values that require less than maxCFL_BT_cont to be accommodated. 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. DYNAMIC_SURFACE_PRESSURE = True ! [Boolean] default = False ! If true, add a dynamic pressure due to a viscous ice shelf, for instance. DT_BT_FILTER = 360 ! [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 ! take to be longer than 2*dt. Set this to 0 to apply no filtering. 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.9 ! [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 = True ! [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. ! === module MOM_mixed_layer_restrat === MIXEDLAYER_RESTRAT = True ! [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. ! FOX_KEMPER_ML_RESTRAT_COEF = 1.0 ! [nondim] default = 0.0 ! A nondimensional coefficient that is proportional to the ratio of the ! deformation radius to the dominant lengthscale of the submesoscale mixed layer ! instabilities, times the minimum of the ratio of the mesoscale eddy kinetic ! energy to the large-scale geostrophic kinetic energy or 1 plus the square of ! the grid spacing over the deformation radius, as detailed by Fox-Kemper et al. ! (2010) MLE_FRONT_LENGTH = 500.0 ! [m] default = 0.0 ! If non-zero, is the frontal-length scale used to calculate the upscaling of ! buoyancy gradients that is otherwise represented by the parameter ! FOX_KEMPER_ML_RESTRAT_COEF. If MLE_FRONT_LENGTH is non-zero, it is recommended ! to set FOX_KEMPER_ML_RESTRAT_COEF=1.0. MLE_USE_PBL_MLD = True ! [Boolean] default = False ! If true, the MLE parameterization will use the mixed-layer depth provided by ! the active PBL parameterization. If false, MLE will estimate a MLD based on a ! density difference with the surface using the parameter MLE_DENSITY_DIFF. MLE_MLD_DECAY_TIME = 2.592E+06 ! [s] default = 0.0 ! The time-scale for a running-mean filter applied to the mixed-layer depth used ! in the MLE restratification parameterization. When the MLD deepens below the ! current running-mean the running-mean is instantaneously set to the current ! MLD. ! === module MOM_diabatic_driver === ! The following parameters are used for diabatic processes. ENERGETICS_SFC_PBL = True ! [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 = False ! [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. ! === module MOM_CVMix_KPP === ! This is the MOM wrapper to CVMix:KPP ! See http://cvmix.github.io/ ! === module MOM_tidal_mixing === ! Vertical Tidal Mixing Parameterization INT_TIDE_DISSIPATION = True ! [Boolean] default = False ! If true, use an internal tidal dissipation scheme to drive diapycnal mixing, ! along the lines of St. Laurent et al. (2002) and Simmons et al. (2004). INT_TIDE_PROFILE = "POLZIN_09" ! default = "STLAURENT_02" ! INT_TIDE_PROFILE selects the vertical profile of energy dissipation with ! INT_TIDE_DISSIPATION. Valid values are: ! STLAURENT_02 - Use the St. Laurent et al exponential ! decay profile. ! POLZIN_09 - Use the Polzin WKB-stretched algebraic ! decay profile. INT_TIDE_DECAY_SCALE = 300.3003003003003 ! [m] default = 500.0 ! The decay scale away from the bottom for tidal TKE with the new coding when ! INT_TIDE_DISSIPATION is used. KAPPA_ITIDES = 6.28319E-04 ! [m-1] default = 6.283185307179586E-04 ! A topographic wavenumber used with INT_TIDE_DISSIPATION. The default is 2pi/10 ! km, as in St.Laurent et al. 2002. KAPPA_H2_FACTOR = 0.84 ! [nondim] default = 1.0 ! A scaling factor for the roughness amplitude with INT_TIDE_DISSIPATION. TKE_ITIDE_MAX = 0.1 ! [W m-2] default = 1000.0 ! The maximum internal tide energy source available to mix above the bottom ! boundary layer with INT_TIDE_DISSIPATION. READ_TIDEAMP = True ! [Boolean] default = False ! If true, read a file (given by TIDEAMP_FILE) containing the tidal amplitude ! with INT_TIDE_DISSIPATION. TIDEAMP_FILE = "ocean_tidal_amplitude.nc" ! default = "tideamp.nc" ! The path to the file containing the spatially varying tidal amplitudes with ! INT_TIDE_DISSIPATION. H2_FILE = "ocean_topog.nc" ! ! The path to the file containing the sub-grid-scale topographic roughness ! amplitude with INT_TIDE_DISSIPATION. ! === module MOM_CVMix_conv === ! Parameterization of enhanced mixing due to convection via CVMix ! === module MOM_geothermal === GEOTHERMAL_SCALE = 1.0 ! [W m-2 or various] default = 0.0 ! The constant geothermal heat flux, a rescaling factor for the heat flux read ! from GEOTHERMAL_FILE, or 0 to disable the geothermal heating. GEOTHERMAL_FILE = "ocean_geothermal.nc" ! default = "" ! The file from which the geothermal heating is to be read, or blank to use a ! constant heating rate. GEOTHERMAL_VARNAME = "geothermal_hf" ! default = "geo_heat" ! The name of the geothermal heating variable in GEOTHERMAL_FILE. ! === module MOM_set_diffusivity === BBL_MIXING_AS_MAX = False ! [Boolean] default = True ! If true, take the maximum of the diffusivity from the BBL mixing and the other ! diffusivities. Otherwise, diffusivity from the BBL_mixing is simply added. USE_LOTW_BBL_DIFFUSIVITY = True ! [Boolean] default = False ! If true, uses a simple, imprecise but non-coordinate dependent, model of BBL ! mixing diffusivity based on Law of the Wall. Otherwise, uses the original BBL ! scheme. SIMPLE_TKE_TO_KD = True ! [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.5E-05 ! [m2 s-1] ! The background diapycnal diffusivity of density in the interior. Zero or the ! molecular value, ~1e-7 m2 s-1, may be used. KD_MIN = 2.0E-06 ! [m2 s-1] default = 1.5E-07 ! The minimum diapycnal diffusivity. HENYEY_IGW_BACKGROUND = True ! [Boolean] default = False ! If true, use a latitude-dependent scaling for the near surface background ! diffusivity, as described in Harrison & Hallberg, JPO 2008. KD_MAX = 0.1 ! [m2 s-1] default = -1.0 ! The maximum permitted increment for the diapycnal diffusivity from TKE-based ! parameterizations, or a negative value for no limit. ! === 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. MAX_RINO_IT = 25 ! [nondim] default = 50 ! The maximum number of iterations that may be used to estimate the Richardson ! number driven mixing. KAPPA_SHEAR_ITER_BUG = True ! [Boolean] default = True ! 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 = True ! [Boolean] default = True ! 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. ! === module MOM_CVMix_shear === ! Parameterization of shear-driven turbulence via CVMix (various options) ! === module MOM_CVMix_ddiff === ! Parameterization of mixing due to double diffusion processes via CVMix ! === module MOM_diabatic_aux === ! The following parameters are used for auxiliary diabatic processes. PRESSURE_DEPENDENT_FRAZIL = True ! [Boolean] default = False ! If true, use a pressure dependent freezing temperature when making frazil. The ! default is false, which will be faster but is inappropriate with ice-shelf ! cavities. VAR_PEN_SW = True ! [Boolean] default = False ! If true, use one of the CHL_A schemes specified by OPACITY_SCHEME to determine ! the e-folding depth of incoming short wave radiation. CHL_FILE = "ocean_chla.nc" ! ! CHL_FILE is the file containing chl_a concentrations in the variable CHL_A. It ! is used when VAR_PEN_SW and CHL_FROM_FILE are true. EPBL_VEL_SCALE_FACTOR = 0.5 ! [nondim] default = 1.0 ! An overall nondimensional scaling factor for wT. Making this larger increases ! the PBL diffusivity. ! === module MOM_CVMix_KPP === ! This is the MOM wrapper to CVMix:KPP ! See http://cvmix.github.io/ USE_KPP = True ! [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. KPP% RI_CRIT = 0.20 ! [nondim] default = 0.3 ! Critical bulk Richardson number used to define depth of the surface Ocean ! Boundary Layer (OBL). USE_KPP_LT_K = False ! [Boolean] default = False ! Flag for Langmuir turbulence enhancement of turbulentmixing coefficient. STOKES_MIXING = False ! [Boolean] default = False ! Flag for Langmuir turbulence enhancement of turbulentmixing coefficient. USE_KPP_LT_VT2 = False ! [Boolean] default = False ! Flag for Langmuir turbulence enhancement of Vt2in Bulk Richardson Number. %KPP ! === module MOM_regularize_layers === ! === module MOM_opacity === PEN_SW_NBANDS = 3 ! default = 1 ! The number of bands of penetrating shortwave radiation. ! === module MOM_tracer_advect === TRACER_ADVECTION_SCHEME = "PPM:H3" ! default = "PLM" ! The horizontal transport scheme for tracers: ! PLM - Piecewise Linear Method ! PPM:H3 - Piecewise Parabolic Method (Huyhn 3rd order) ! PPM - Piecewise Parabolic Method (Colella-Woodward) ! === module MOM_tracer_hor_diff === CHECK_DIFFUSIVE_CFL = True ! [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. ! === module MOM_neutral_diffusion === ! This module implements neutral diffusion of tracers ! === module MOM_lateral_boundary_diffusion === ! This module implements lateral diffusion of tracers near boundaries ! === module MOM_sum_output === MAXTRUNC = 100000 ! [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. ENERGYSAVEDAYS = 0.25 ! [days] default = 1.0 ! The interval in units of TIMEUNIT between saves of the energies of the run and ! other globally summed diagnostics. ! === module ocean_model_init === ! === module MOM_surface_forcing === BUOY_CONFIG = "data_override" ! default = "zero" ! The character string that indicates how buoyancy forcing is specified. Valid ! options include (file), (zero), (linear), (USER), (BFB) and (NONE). ARCHAIC_OMIP_FORCING_FILE = False ! [Boolean] default = True ! If true, use the forcing variable decomposition from the old German OMIP ! prescription that predated CORE. If false, use the variable groupings ! available from MOM output diagnostics of forcing variables. MAX_P_SURF = 0.0 ! [Pa] default = -1.0 ! The maximum surface pressure that can be exerted by the atmosphere and ! floating sea-ice or ice shelves. This is needed because the FMS coupling ! structure does not limit the water that can be frozen out of the ocean and the ! ice-ocean heat fluxes are treated explicitly. No limit is applied if a ! negative value is used. CD_TIDES = 0.0018 ! [nondim] default = 1.0E-04 ! The drag coefficient that applies to the tides. WIND_CONFIG = "data_override" ! default = "zero" ! The character string that indicates how wind forcing is specified. Valid ! options include (file), (2gyre), (1gyre), (gyres), (zero), and (USER). RESTOREBUOY = True ! [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. !FLUXCONST = 0.5 ! [m day-1] default = 0.0 ! The constant that relates the restoring surface fluxes to the relative surface ! anomalies (akin to a piston velocity). Note the non-MKS units. !FLUXCONST_T = 0.5 ! [m day-1] default = 0.5 ! The constant that relates the restoring surface temperature flux to the ! relative surface anomaly (akin to a piston velocity). Note the non-MKS units. !FLUXCONST_S = 0.5 ! [m day-1] default = 0.5 ! The constant that relates the restoring surface salinity flux to the relative ! surface anomaly (akin to a piston velocity). Note the non-MKS units. GUST_CONST = 0.0 ! [Pa] default = 0.0 ! The background gustiness in the winds. FIX_USTAR_GUSTLESS_BUG = True ! [Boolean] default = True ! If true correct a bug in the time-averaging of the gustless wind friction ! velocity ! These parameters are used in CORE mode but should not be used in ! the coupled model (CM4). RESTORE_SALINITY = False ! [Boolean] default = False ! If true, the coupled driver will add a globally-balanced ! fresh-water flux that drives sea-surface salinity ! toward specified values. ADJUST_NET_FRESH_WATER_TO_ZERO = True ! [Boolean] default = False ! If true, adjusts the net fresh-water forcing seen ! by the ocean (including restoring) to zero. FLUXCONST = 0.1667 !0.500 !0.1667 ! [m day-1] ! The constant that relates the restoring surface fluxes ! to the relative surface anomalies (akin to a piston ! velocity). Note the non-MKS units. SALT_RESTORE_FILE = "ocean_salt_restore.nc" !"salt_restore_correct.nc" ! default = "salt_restore.nc" ! A file in which to find the surface salinity to use for restoring. SALT_RESTORE_VARIABLE = "SALT" SRESTORE_AS_SFLUX = True ! [Boolean] default = False ! If true, the restoring of salinity is applied as a salt ! flux instead of as a freshwater flux. MAX_DELTA_SRESTORE = 5.0 ! [PSU or g kg-1] default = 999.0 ! The maximum salinity difference used in restoring terms. USE_NET_FW_ADJUSTMENT_SIGN_BUG = False ! [Boolean] default = True ! If true, use the wrong sign for the adjustment to ! the net fresh-water. ! === module MOM_restart === ! === module MOM_file_parser ===