diff --git a/CITATION.cff b/CITATION.cff index 95c6b2e00..3a0d60b95 100644 --- a/CITATION.cff +++ b/CITATION.cff @@ -102,6 +102,34 @@ authors: family-names: Yang email: Luwei.Yang@anu.edu.au affiliation: Research School of Earth Sciences, Australian National University + - given-names: Claire + family-names: Yung + email: Claire.Yung@anu.edu.au + affiliation: Research School of Earth Sciences, Australian National University + orcid: 'https://orcid.org/0000-0002-0052-7668' + - given-names: Wilton + family-names: Aguiar + email: wilton.aguiar@anu.edu.au + affiliation: Research School of Earth Sciences, Australian National University + orcid: 'https://orcid.org/0000-0003-2453-9791' + - given-names: Helen + family-names: Macdonald + email: helen.macdonald@anu.edu.au + affiliation: ACCESS-NRI + orcid: 'https://orcid.org/0000-0003-4734-3150' + - given-names: Edward + family-names: Yang + email: edward.yang@anu.edu.au + affiliation: ACCESS-NRI + - given-names: David + family-names: Gwyther + email: d.gwyther@uq.edu.au + affiliation: School of the Environment, University of Queensland + orcid: 'https://orcid.org/0000-0002-7218-2785' + - given-names: Fabio + family-names: Boeira Dias + email: f.boeira_dias@unsw.edu.au + orcid: 'https://orcid.org/0000-0002-2965-2120' identifiers: - description: "The concept DOI of all ACCESS-OM3 configurations." type: doi diff --git a/MOM_input b/MOM_input index ef16bbe15..457072464 100644 --- a/MOM_input +++ b/MOM_input @@ -10,58 +10,59 @@ V_TRUNC_FILE = "V_velocity_truncations" ! default = "" ! velocity truncations are written. Undefine this for efficiency if this ! diagnostic is not needed. -! === module MOM_file_parser === -FATAL_UNUSED_PARAMS = True ! [Boolean] default = False - ! If true, kill the run if there are any unused parameters. - ! This section contains non default layout parameters which are set. ! ------------------------ ! === module MOM_domains === -AUTO_MASKTABLE = True ! [Boolean] default = False - ! Turn on automatic mask table generation to eliminate land blocks. -AUTO_IO_LAYOUT_FAC = 6 ! default = 0 - ! When AUTO_MASKTABLE is enabled, io layout is calculated by performing integer - ! division of the runtime-determined domain layout with this factor. If the - ! factor is set to 0 (default), the io layout is set to 1,1. +MASKTABLE = "mask_table.1490.84x49" ! default = "MOM_mask_table" + ! A text file to specify n_mask, layout and mask_list. This feature masks out + ! processors that contain only land points. The first line of mask_table is the + ! number of regions to be masked out. The second line is the layout of the model + ! and must be consistent with the actual model layout. The following (n_mask) + ! lines give the logical positions of the processors that are masked out. The + ! mask_table can be created by tools like check_mask. The following example of + ! mask_table masks out 2 processors, (1,2) and (3,6), out of the 24 in a 4x6 + ! layout: + ! 2 + ! 4,6 + ! 1,2 + ! 3,6 +LAYOUT = 84, 49 ! + ! The processor layout that was actually used. +IO_LAYOUT = 12, 7 ! default = 1, 1 + ! The processor layout to be used, or 0,0 to automatically set the io_layout to + ! be the same as the layout. ! This section is formatted the same as the MOM_parameter_short.doc output file. ! ------------------------ ! This file was written by the model and records the non-default parameters used at run-time. ! === module MOM === -USE_CONTEMP_ABSSAL = True ! [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. USE_REGRIDDING = True ! [Boolean] default = False ! If True, use the ALE algorithm (regridding/remapping). If False, use the ! layered isopycnal algorithm. -DT = 900.0 ! [s] +DT = 150.0 ! [s] ! The (baroclinic) dynamics time step. The time-step that is actually used will ! be an integer fraction of the forcing time-step (DT_FORCING in ocean-only mode ! or the coupling timestep in coupled mode.) -DT_THERM = 7200.0 ! [s] default = 900.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 = True ! [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. -THICKNESSDIFFUSE = True ! [Boolean] default = False - ! If true, isopycnal surfaces are diffused with a Laplacian coefficient of KHTH. -THICKNESSDIFFUSE_FIRST = True ! [Boolean] default = False - ! If true, do thickness diffusion or interface height smoothing before dynamics. - ! This is only used if THICKNESSDIFFUSE or APPLY_INTERFACE_FILTER is true. +HMIX_SFC_PROP = 1.08 ! [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 = 1.08 ! [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 = 10.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. -DTBT_RESET_PERIOD = 0.0 ! [s] default = 7200.0 +DTBT_RESET_PERIOD = 0.0 ! [s] default = 150.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 @@ -77,20 +78,9 @@ C_P = 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. +USE_PSURF_IN_EOS = False ! [Boolean] default = True + ! If true, always include the surface pressure contributions in equation of + ! state calculations. WRITE_GEOM = 2 ! 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 @@ -99,13 +89,10 @@ GEOM_FILE = "access-om3.mom6.geometry.nc" ! default = "ocean_geometry" ! The file into which to write the ocean geometry. ! === module MOM_domains === -TRIPOLAR_N = True ! [Boolean] default = False - ! Use tripolar connectivity at the northern edge of the domain. With - ! TRIPOLAR_N, NIGLOBAL must be even. -NIGLOBAL = 1440 ! +NIGLOBAL = 4320 ! ! 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 = 1152 ! +NJGLOBAL = 1442 ! ! 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. @@ -156,21 +143,45 @@ TOPO_CONFIG = "file" ! MAXIMUM_DEPTH = 6000.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 = 1 ! default = 0 + ! The number of open boundary segments. +OBC_ZERO_BIHARMONIC = True ! [Boolean] default = False + ! If true, zeros the Laplacian of flow on open boundaries in the biharmonic + ! viscosity term. +OBC_SEGMENT_001 = "J=N,I=N:0,FLATHER,ORLANSKI,NUDGED" ! + ! Documentation needs to be dynamic????? +OBC_SEGMENT_001_VELOCITY_NUDGING_TIMESCALES = 0.3, 360.0 ! [days] + ! 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_TRACER_RESERVOIR_LENGTH_SCALE_OUT = 3.0E+04 ! [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 = 3000.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. +BRUSHCUTTER_MODE = True ! [Boolean] default = False + ! If true, read external OBC data on the supergrid. + ! === module MOM_verticalGrid === ! Parameters providing information about the vertical grid. +ANGSTROM = 1.0E-15 ! [m] default = 1.0E-10 + ! The minimum layer thickness, usually one-Angstrom. NK = 75 ! [nondim] ! The number of model layers. ! === module MOM_EOS === -EQN_OF_STATE = "ROQUET_RHO" ! default = "WRIGHT" +EQN_OF_STATE = "WRIGHT_FULL" ! default = "WRIGHT" ! EQN_OF_STATE determines which ocean equation of state should be used. ! Currently, the valid choices are "LINEAR", "UNESCO", "JACKETT_MCD", "WRIGHT", ! "WRIGHT_REDUCED", "WRIGHT_FULL", "NEMO", "ROQUET_RHO", "ROQUET_SPV" and ! "TEOS10". This is only used if USE_EOS is true. -TFREEZE_FORM = "TEOS_POLY" ! default = "TEOS10" - ! TFREEZE_FORM determines which expression should be used for the freezing - ! point. Currently, the valid choices are "LINEAR", "MILLERO_78", "TEOS_POLY", - ! "TEOS10" +DTFREEZE_DP = -7.75E-08 ! [degC 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 = True ! [Boolean] default = False @@ -184,12 +195,18 @@ USE_IDEAL_AGE_TRACER = True ! [Boolean] default = False ! === module ideal_age_example === +! === module MOM_boundary_update === + ! === module MOM_fixed_initialization === ! === module MOM_grid_init === +! === module MOM_open_boundary === +! Controls where open boundaries are located, what kind of boundary condition to impose, and what data to apply, +! if any. + ! === module MOM_coord_initialization === -REGRIDDING_COORDINATE_MODE = "ZSTAR" ! default = "LAYER" +REGRIDDING_COORDINATE_MODE = "SIGMA_SHELF_ZSTAR" ! default = "LAYER" ! Coordinate mode for vertical regridding. Choose among the following ! possibilities: LAYER - Isopycnal or stacked shallow water layers ! ZSTAR, Z* - stretched geopotential z* @@ -216,6 +233,8 @@ ALE_COORDINATE_CONFIG = "FILE:ocean_vgrid.nc,interfaces=zeta" ! default = "UNIFO ! by a comma or space, for sigma-2 and dz. e.g. ! HYBRID:vgrid.nc,sigma2,dz !ALE_RESOLUTION = 1.0825614929199219, 1.1963462829589844, 1.322089672088623, 1.4610481262207031, 1.614609718322754, 1.784308910369873, 1.9718408584594727, 2.1790781021118164, 2.4080896377563477, 2.661160469055176, 2.940814971923828, 3.249845504760742, 3.591329574584961, 3.968667984008789, 4.385614395141602, 4.846321105957031, 5.355350494384766, 5.917766571044922, 6.539115905761719, 7.225547790527344, 7.983818054199219, 8.821372985839844, 9.746376037597656, 10.767845153808594, 11.895652770996094, 13.140586853027344, 14.514511108398438, 16.030319213867188, 17.7020263671875, 19.544876098632812, 21.575271606445312, 23.810821533203125, 26.270294189453125, 28.973419189453125, 31.94091796875, 35.194000244140625, 38.75390625, 42.641632080078125, 46.876739501953125, 51.476593017578125, 56.45489501953125, 61.82025146484375, 67.5743408203125, 73.70965576171875, 80.207763671875, 87.03759765625, 94.1534423828125, 101.4951171875, 108.9879150390625, 116.5452880859375, 124.0714111328125, 131.4671630859375, 138.6346435546875, 145.484130859375, 151.938232421875, 157.93701171875, 163.439697265625, 168.42431640625, 172.8876953125, 176.842041015625, 180.3125, 183.33154296875, 185.938720703125, 188.175048828125, 190.08251953125, 191.701171875 +MIN_THICKNESS = 1.0E-12 ! [m] default = 0.001 + ! When regridding, this is the minimum layer thickness allowed. 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: @@ -228,6 +247,9 @@ REMAPPING_SCHEME = "PPM_H4" ! default = "PLM" ! WENO_HYBGEN (3rd-order accurate) ! PQM_IH4IH3 (4th-order accurate) ! PQM_IH6IH5 (5th-order accurate) +INIT_BOUNDARY_EXTRAP = True ! [Boolean] default = False + ! If true, values at the interfaces of boundary cells are extrapolated instead + ! of piecewise constant during initialization.Defaults to REMAP_BOUNDARY_EXTRAP. REMAPPING_USE_OM4_SUBCELLS = False ! [Boolean] default = True ! This selects the remapping algorithm used in OM4 that does not use the full ! reconstruction for the top- and lower-most sub-layers, but instead assumes @@ -240,10 +262,12 @@ INIT_LAYERS_FROM_Z_FILE = True ! [Boolean] default = False ! Z-space file on a latitude-longitude grid. ! === module MOM_initialize_layers_from_Z === -TEMP_SALT_Z_INIT_FILE = "woa23_ts_01_mom.nc" ! default = "temp_salt_z.nc" +TEMP_SALT_Z_INIT_FILE = "ocean_temp_salt.res.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. +Z_INIT_FILE_PTEMP_VAR = "temp" ! default = "ptemp" + ! The name of the potential temperature variable in TEMP_Z_INIT_FILE. Z_INIT_ALE_REMAPPING = True ! [Boolean] default = False ! If True, then remap straight to model coordinate from file. Z_INIT_REMAP_GENERAL = True ! [Boolean] default = False @@ -253,6 +277,27 @@ TEMP_SALT_INIT_VERTICAL_REMAP_ONLY = True ! [Boolean] default = False ! If true, initial conditions are on the model horizontal grid. Extrapolation ! over missing ocean values is done using an ICE-9 procedure with vertical ALE ! remapping . +TRIM_IC_FOR_P_SURF = True ! [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. +SURFACE_PRESSURE_FILE = "ice_shelf_thick.nc" ! + ! The initial condition file for the surface pressure exerted by ice. +SURFACE_PRESSURE_VAR = "thick" ! default = "" + ! The initial condition variable for the surface pressure exerted by ice. +SURFACE_PRESSURE_SCALE = 8986.6 ! [file dependent] default = 1.0 + ! A scaling factor to convert SURFACE_PRESSURE_VAR from file + ! SURFACE_PRESSURE_FILE into a surface pressure. +TRIM_IC_Z_TOLERANCE = 1.0E-10 ! [m] default = 1.0E-05 + ! The tolerance with which to find the depth matching the specified surface + ! pressure with TRIM_IC_FOR_P_SURF. +FRAC_DP_AT_POS_NEGATIVE_P_BUGFIX = True ! [Boolean] default = False + ! If true, use bugfix in ice shelf TRIM_IC initialization. Otherwise, pressure + ! input to density EOS is negative. +TRIMMING_USES_REMAPPING = True ! [Boolean] default = False + ! When trimming the column, also remap T and S. +OBC_SEGMENT_001_DATA = "U=file:OBC37S_forcing_access_yr2_4km.nc(u),V=file:OBC37S_forcing_access_yr2_4km.nc(v),SSH=file:OBC37S_forcing_access_yr2_4km.nc(eta_t),TEMP=file:OBC37S_forcing_access_yr2_4km.nc(pot_temp),SALT=file:OBC37S_forcing_access_yr2_4km.nc(salt)" ! + ! OBC segment docs ! === module MOM_diag_mediator === NUM_DIAG_COORDS = 2 ! default = 1 @@ -283,7 +328,7 @@ REGRIDDING_ANSWER_DATE = 99991231 ! default = 20181231 ! 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. -DIAG_COORD_DEF_RHO2 = "RFNC1:76,999.5,1020.,1034.1,3.1,1041.,0.002" ! default = "WOA09" +DIAG_COORD_DEF_RHO2 = "FILE:diag_rho2.nc,interfaces=rho2" ! 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 @@ -300,70 +345,24 @@ DIAG_COORD_DEF_RHO2 = "RFNC1:76,999.5,1020.,1034.1,3.1,1041.,0.002" ! default = ! 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_KHTH_FAC = 0.5 ! [nondim] default = 0.0 - ! A factor that maps MEKE%Kh to KhTh. -MEKE_KHTR_FAC = 0.5 ! [nondim] default = 0.0 - ! A factor that maps MEKE%Kh to KhTr. -MEKE_KHMEKE_FAC = 1.0 ! [nondim] default = 0.0 - ! A factor that maps MEKE%Kh to Kh for MEKE itself. -MEKE_VISCOSITY_COEFF_KU = 1.0 ! [nondim] default = 0.0 - ! If non-zero, is the scaling coefficient in the expression forviscosity used to - ! parameterize harmonic lateral momentum mixing byunresolved eddies represented - ! by MEKE. Can be negative torepresent backscatter from the unresolved eddies. -MEKE_ALPHA_RHINES = 0.15 ! [nondim] default = 0.0 - ! 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.0 - ! 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. -KHTH_USE_EBT_STRUCT = True ! [Boolean] default = False - ! If true, uses the equivalent barotropic structure as the vertical structure of - ! thickness diffusivity. 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. -KH_RES_FN_POWER = 100 ! default = 2 - ! The power of dx/Ld in the Kh resolution function. Any positive integer may be - ! used, although even integers are more efficient to calculate. Setting this - ! greater than 100 results in a step-function being used. ! === 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. +LINEAR_DRAG = True ! [Boolean] default = False + ! If LINEAR_DRAG and BOTTOMDRAGLAW are defined the drag law is + ! cdrag*DRAG_BG_VEL*u. 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. -DRAG_BG_VEL = 0.1 ! [m s-1] default = 0.0 +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 @@ -372,7 +371,7 @@ 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 = 0.0 ! [m2 s-1] +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. @@ -382,12 +381,6 @@ KHTH_MAX_CFL = 0.1 ! [nondimensional] default = 0.8 ! 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. -KHTH_USE_FGNV_STREAMFUNCTION = True ! [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. -FGNV_FILTER_SCALE = 0.1 ! [nondim] default = 1.0 - ! A coefficient scaling the vertical smoothing term in the Ferrari et al., 2010, - ! streamfunction formulation. ! === module MOM_dynamics_split_RK2 === VISC_REM_BUG = False ! [Boolean] default = True @@ -397,7 +390,7 @@ VISC_REM_BUG = False ! [Boolean] default = True ! VISC_REM_BT_WEIGHT_BUG in MOM_barotropic. ! === module MOM_continuity_PPM === -ETA_TOLERANCE = 1.0E-06 ! [m] default = 3.75E-09 +ETA_TOLERANCE = 1.0E-06 ! [m] default = 3.75E-14 ! 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 @@ -429,6 +422,25 @@ BOUND_CORIOLIS = True ! [Boolean] default = False MASS_WEIGHT_IN_PRESSURE_GRADIENT = True ! [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 = True ! [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 = True ! [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. +RESET_INTXPA_INTEGRAL = True ! [Boolean] default = False + ! If true, reset INTXPA to match pressures at first nonvanished cell. Includes + ! pressure correction. +RESET_INTXPA_INTEGRAL_FLATTEST = True ! [Boolean] default = False + ! If true, use flattest interface as reference interface where there is no + ! better choice for RESET_INTXPA_INTEGRAL. Otherwise, use surface interface. +PRESSURE_RECONSTRUCTION_SCHEME = 2 ! 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. ! === module MOM_Zanna_Bolton === @@ -439,9 +451,6 @@ KH_VEL_SCALE = 0.01 ! [m s-1] default = 0.0 ! The velocity scale which is multiplied by the grid spacing to calculate the ! Laplacian viscosity. The final viscosity is the largest of this scaled ! viscosity, the Smagorinsky and Leith viscosities, and KH. -KH_SIN_LAT = 2000.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). AH_VEL_SCALE = 0.01 ! [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 @@ -463,6 +472,10 @@ FRICTWORK_BUG = False ! [Boolean] default = True ! recommended. ! === module MOM_vert_friction === +HARMONIC_BL_SCALE = 1.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 = 0.5 ! [m] ! The prescribed depth over which the near-surface viscosity and diffusivity are ! elevated when the bulk mixed layer is not used. @@ -542,6 +555,14 @@ ENERGETICS_SFC_PBL = True ! [Boolean] default = False 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. +MINIMUM_FORCING_DEPTH = 1.08 ! [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. +FRAZIL_NOT_UNDER_ICESHELF = True ! [Boolean] default = False + ! If true, do not use frazil scheme underneath ice shelves defined by + ! frac_shelf_h greater than 0. ! === module MOM_set_diffusivity === @@ -550,24 +571,6 @@ EPBL_IS_ADDITIVE = False ! [Boolean] default = True 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. @@ -589,10 +592,6 @@ LOTW_BBL_ANSWER_DATE = 99991231 ! default = 20190101 ! calculations. Values below 20240630 recover the original answers, while ! higher values use more accurate expressions. This only applies when ! USE_LOTW_BBL_DIFFUSIVITY is true. -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 @@ -644,14 +643,6 @@ USE_RIVER_HEAT_CONTENT = True ! [Boolean] default = False USE_CALVING_HEAT_CONTENT = True ! [Boolean] default = False ! If true, use the fluxes%calving_Hflx field to set the heat carried by runoff, ! instead of using SST*CP*froz_runoff. -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 = "chl_globcolour_monthly_clim.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. -CHL_VARNAME = "CHL" ! default = "CHL_A" - ! Name of CHL_A variable in CHL_FILE. ! === module MOM_energetic_PBL === ML_OMEGA_FRAC = 0.001 ! [nondim] default = 0.0 @@ -715,8 +706,11 @@ LT_MOD_LAC5 = 0.22 ! [nondim] default = 0.95 ! unstable Obukhov depth. ! === module MOM_opacity === -PEN_SW_NBANDS = 3 ! default = 1 - ! The number of bands of penetrating shortwave radiation. +PEN_SW_SCALE = 15.0 ! [m] default = 0.0 + ! The vertical absorption e-folding depth of the penetrating shortwave + ! radiation. +PEN_SW_FRAC = 0.42 ! [nondim] default = 0.0 + ! The fraction of the shortwave radiation that penetrates below the surface. ! === module MOM_tracer_advect === TRACER_ADVECTION_SCHEME = "PPM:H3" ! default = "PLM" @@ -726,8 +720,6 @@ TRACER_ADVECTION_SCHEME = "PPM:H3" ! default = "PLM" ! PPM - Piecewise Parabolic Method (Colella-Woodward) ! === module MOM_tracer_hor_diff === -KHTR = 50.0 ! [m2 s-1] default = 0.0 - ! The background along-isopycnal tracer diffusivity. 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 @@ -737,13 +729,8 @@ MAX_TR_DIFFUSION_CFL = 2.0 ! [nondim] default = -1.0 ! diffusive CFL locally at or below this value. The number of diffusive ! iterations is often this value or the next greater integer. -! === module MOM_neutral_diffusion === -! This module implements neutral diffusion of tracers -USE_NEUTRAL_DIFFUSION = True ! [Boolean] default = False - ! If true, enables the neutral diffusion module. - ! === module MOM_sum_output === -MAXTRUNC = 10000 ! [truncations save_interval-1] default = 0 +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. @@ -761,15 +748,23 @@ EPS_OMESH = 1.0E-13 ! [degrees] default = 1.0E-04 RESTORE_SALINITY = True ! [Boolean] default = False ! If true, the coupled driver will add a globally-balanced fresh-water flux that ! drives sea-surface salinity toward specified values. +ICE_SHELF = True ! [Boolean] default = False + ! If true, enables the ice shelf model. ! === module MOM_surface_forcing_nuopc === LATENT_HEAT_FUSION = 3.337E+05 ! [J/kg] default = 3.34E+05 ! The latent heat of fusion. LATENT_HEAT_VAPORIZATION = 2.501E+06 ! [J/kg] default = 2.5E+06 ! The latent heat of fusion. -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. +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. +ADJUST_NET_SRESTORE_TO_ZERO = False ! [Boolean] default = True + ! If true, adjusts the salinity restoring seen to zero whether restoring is via + ! a salt flux or virtual precip. WIND_STAGGER = "A" ! default = "C" ! A case-insensitive character string to indicate the staggering of the input ! wind stress field. Valid values are 'A', 'B', or 'C'. @@ -781,12 +776,11 @@ FLUXCONST = 0.11 ! [m day-1] default = 0.0 ! anomalies (akin to a piston velocity). Note the non-MKS units. SALT_RESTORE_FILE = "salt_sfc_restore.nc" ! default = "salt_restore.nc" ! A file in which to find the surface salinity to use for restoring. -SALT_RESTORE_VARIABLE = "asalt" ! default = "salt" - ! The name of the surface salinity variable to read from SALT_RESTORE_FILE for - ! restoring salinity. 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. +MASK_SRESTORE = True ! [Boolean] default = False + ! If true, read a file (salt_restore_mask) containing a mask for SSS restoring. USE_RIGID_SEA_ICE = True ! [Boolean] default = False ! If true, sea-ice is rigid enough to exert a nonhydrostatic pressure that ! resist vertical motion. @@ -795,3 +789,41 @@ SEA_ICE_RIGID_MASS = 100.0 ! [kg m-2] default = 1000.0 ! rigidity ! === module MOM_restart === + +! === module MOM_domains min_halo === + +! === module MOM_grid_init === + +! === module MOM_ice_shelf === +SHELF_THERMO = True ! [Boolean] default = False + ! If true, use a thermodynamically interactive ice shelf. +ICE_SHELF_TEMPERATURE = -20.0 ! [degC] default = -15.0 + ! The temperature at the center of the ice shelf. +COL_THICK_MELT_THRESHOLD = 0.001 ! [m] default = 0.0 + ! The minimum ocean column thickness where melting is allowed. +ICE_SHELF_USTAR_FROM_VEL_BUGFIX = True ! [Boolean] default = False + ! Bug fix for ice-area weighting of squared ocean velocities used to calculate + ! friction velocity under ice shelves +ICE_SHELF_TIDEAMP_SCALING_FACTOR = 0.66 ! [nondim] default = 1.0 + ! Scale TIDEAMP_FILE or UTIDE by number in melt parameterisation friction + ! velocity calculation. + +! === module MOM_EOS === +DENSITY_ICE = 917.0 ! [kg m-3] default = 900.0 + ! A typical density of ice. +USTAR_SHELF_BG = 6.0E-04 ! [m s-1] default = 0.0 + ! The minimum value of ustar under ice shelves. +CDRAG_SHELF = 0.007 ! [nondim] default = 0.003 + ! CDRAG is the drag coefficient relating the magnitude of the velocity field to + ! the surface stress. + +! === module MOM_restart === +ICE_PROFILE_CONFIG = "FILE" ! + ! This specifies how the initial ice profile is specified. Valid values are: + ! CHANNEL, FILE, and USER. +ICE_THICKNESS_FILE = "ice_shelf_thick.nc" ! default = "ice_shelf_h.nc" + ! The file from which the bathymetry is read. +ICE_THICKNESS_VARNAME = "thick" ! default = "h_shelf" + ! The name of the thickness variable in ICE_THICKNESS_FILE. +ICE_AREA_VARNAME = "area" ! default = "area_shelf_h" + ! The name of the area variable in ICE_THICKNESS_FILE. diff --git a/README.md b/README.md index 10b4fd07d..ac2f7f63a 100644 --- a/README.md +++ b/README.md @@ -1,3 +1,7 @@ +## MOM6-CICE6 4km JRA55-do RYF regional pan-Antarctic ACCESS-OM3 configuration with ice shelf cavities + +**WARNING: This configuration is still under development and should not be used for production.** + See [`main` branch README](https://github.com/ACCESS-NRI/access-om3-configs/blob/main/README.md) for usage information. \ No newline at end of file diff --git a/config.yaml b/config.yaml index a1bf4a32e..c3cd40aa4 100644 --- a/config.yaml +++ b/config.yaml @@ -2,9 +2,9 @@ # If submitting to a different project to your default, uncomment line below # and change project code as appropriate; also set shortpath below -# project: x77 +# project: e14 # Force payu to always find, and save, files in this scratch project directory -# shortpath: /scratch/v45 +# shortpath: /scratch/x77 # Using payu sync is recommend to copy data from ephemeral scratch space to # longer term storage @@ -18,35 +18,41 @@ sync: modules: use: - /g/data/vk83/modules + - /g/data/vk83/prerelease/modules load: - - access-om3/2025.08.001 - - model-tools/mppnccombine-fast + - access-om3/pr142-36 + - model-tools/mppnccombine-fast/2025.07.000 queue: normalsr -ncpus: 1976 +ncpus: 2990 jobfs: 10GB -mem: 9500GB -walltime: 6:00:00 -jobname: 25km_jra_ryf +mem: 14700GB +walltime: 4:00:00 +jobname: 4km_panAn_om3_isf model: access-om3 exe: access-om3-MOM6-CICE6 input: - - /g/data/vk83/configurations/inputs/access-om3/cice/grids/global.25km/2025.11.27/kmt.nc - - /g/data/vk83/configurations/inputs/access-om3/mom/grids/mosaic/global.25km/2025.09.02/ocean_hgrid.nc - - /g/data/vk83/configurations/inputs/access-om3/mom/grids/vertical/global.25km/2025.03.12/ocean_vgrid.nc - - /g/data/vk83/configurations/inputs/access-om3/mom/initial_conditions/global.25km/2025.10.24/woa23_ts_01_mom.nc - - /g/data/vk83/configurations/inputs/access-om3/mom/surface_salt_restoring/global.25km/2025.10.24/salt_sfc_restore.nc - - /g/data/vk83/configurations/inputs/access-om3/mom/tidal_external_files/global.25km/2025.09.03/bottom_roughness.nc - - /g/data/vk83/configurations/inputs/access-om3/mom/tidal_external_files/global.25km/2025.09.03/tideamp.nc - - /g/data/vk83/configurations/inputs/access-om3/share/grids/global.25km/2025.11.27/topog.nc - - /g/data/vk83/configurations/inputs/access-om3/mom/chlorophyll/global.25km/2025.11.21/chl_globcolour_monthly_clim.nc - - /g/data/vk83/configurations/inputs/access-om3/share/meshes/global.25km/2025.11.27/access-om3-25km-ESMFmesh.nc - - /g/data/vk83/configurations/inputs/access-om3/share/meshes/global.25km/2025.11.27/access-om3-25km-nomask-ESMFmesh.nc + - /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/kmt.nc + - /g/data/vk83/prerelease/configurations/inputs/access-om3/panan.4km/2026.01.08/ocean_hgrid.nc + - /g/data/vk83/prerelease/configurations/inputs/access-om3/panan.4km/2026.01.08/ocean_vgrid.nc + - /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/ocean_temp_salt.res.nc + - /g/data/vk83/prerelease/configurations/inputs/access-om3/panan.4km/2026.01.08/salt_sfc_restore.nc + - /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/salt_restore_mask.nc + - /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/bottom_roughness.nc + - /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/tideamp.nc + - /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/topog_seaice.nc + - /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/topog.nc + - /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/access-om3-4km-ESMFmesh.nc + - /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/access-om3-4km-nomask-ESMFmesh.nc - /g/data/vk83/prerelease/configurations/inputs/access-om3/share/meshes/share/2026.01.21/JRA55do-datm-ESMFmesh.nc - /g/data/vk83/prerelease/configurations/inputs/access-om3/share/meshes/share/2026.01.21/JRA55do-drof-ESMFmesh.nc - - /g/data/vk83/configurations/inputs/access-om3/cmeps/remap_weights/global.25km/2025.11.27/access-om3-25km-rof-remap-weights.nc + - /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/access-om3-4km-rof-remap-weights.nc - /g/data/vk83/configurations/inputs/JRA-55/RYF/v1-6/data + - /g/data/vk83/prerelease/configurations/inputs/access-om3/panan.4km/2026.01.08/OBC37S_forcing_access_yr2_4km.nc + - /g/data/vk83/prerelease/configurations/inputs/access-om3/panan.4km/2026.01.08/diag_rho2.nc + - /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/ice_shelf_thick.nc + - /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/mask_table.1490.84x49 runlog: true metadata: @@ -63,9 +69,9 @@ platform: collate: restart: true mpi: true - walltime: 4:00:00 - mem: 30GB - ncpus: 4 + walltime: 1:00:00 + mem: 80GB + ncpus: 16 queue: normalsr exe: mppnccombine-fast @@ -76,5 +82,5 @@ postscript: -v PROJECT,SCRIPTS_DIR=/g/data/vk83/apps/om3-scripts -lstorage=${PBS qsub_flags: -W umask=027 -restart_freq: 1YS +restart_freq: 1 payu_minimum_version: 1.2.0 diff --git a/datm.streams.xml b/datm.streams.xml index a0c86274e..de490178a 100644 --- a/datm.streams.xml +++ b/datm.streams.xml @@ -67,7 +67,7 @@ cycle - 1.0 + 2 single bilinear 1900 diff --git a/datm_in b/datm_in index 8dd0b9566..42f36741f 100644 --- a/datm_in +++ b/datm_in @@ -8,10 +8,10 @@ flds_preso3 = .false. flds_wiso = .false. iradsw = 1 - model_maskfile = "./INPUT/access-om3-25km-nomask-ESMFmesh.nc" - model_meshfile = "./INPUT/access-om3-25km-nomask-ESMFmesh.nc" - nx_global = 1440 - ny_global = 1152 + model_maskfile = "./INPUT/access-om3-4km-nomask-ESMFmesh.nc" + model_meshfile = "./INPUT/access-om3-4km-nomask-ESMFmesh.nc" + nx_global = 4320 + ny_global = 1442 restfilm = "null" skip_restart_read = .false. / diff --git a/diag_table b/diag_table index 3daef19e8..324ca6106 100644 --- a/diag_table +++ b/diag_table @@ -43,8 +43,8 @@ ACCESS-OM3 # monthly 3d fields on z -"access-om3.mom6.3d.thkcello+z.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model_z", "thkcello", "thkcello", "access-om3.mom6.3d.thkcello+z.1mon.mean.%4yr", "all", "average", "none", 2 +"access-om3.mom6.3d.h+z.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" +"ocean_model_z", "h", "h", "access-om3.mom6.3d.h+z.1mon.mean.%4yr", "all", "average", "none", 2 "access-om3.mom6.3d.agessc+z.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" "ocean_model_z", "agessc", "agessc", "access-om3.mom6.3d.agessc+z.1mon.mean.%4yr", "all", "average", "none", 2 @@ -70,11 +70,11 @@ ACCESS-OM3 "access-om3.mom6.3d.KE+z.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" "ocean_model_z", "KE", "KE", "access-om3.mom6.3d.KE+z.1mon.mean.%4yr", "all", "average", "none", 2 -"access-om3.mom6.3d.Kd_Itidal_Work+z.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model_z", "Kd_Itidal_Work", "Kd_Itidal_Work", "access-om3.mom6.3d.Kd_Itidal_Work+z.1mon.mean.%4yr", "all", "average", "none", 2 +"access-om3.mom6.3d.e+z.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" +"ocean_model_z", "e", "e", "access-om3.mom6.3d.e+z.1mon.mean.%4yr", "all", "average", "none", 2 -"access-om3.mom6.3d.rsdoabsorb+z.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model_z", "rsdoabsorb", "rsdoabsorb", "access-om3.mom6.3d.rsdoabsorb+z.1mon.mean.%4yr", "all", "average", "none", 2 +"access-om3.mom6.3d.Kd_interface+z.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" +"ocean_model_z", "Kd_interface", "Kd_interface", "access-om3.mom6.3d.Kd_interface+z.1mon.mean.%4yr", "all", "average", "none", 2 # monthly 3d fields on rho2 @@ -91,30 +91,6 @@ ACCESS-OM3 "access-om3.mom6.3d.vmo+rho2.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" "ocean_model_rho2", "vmo", "vmo", "access-om3.mom6.3d.vmo+rho2.1mon.mean.%4yr", "all", "average", "none", 2 -"access-om3.mom6.3d.uhGM+rho2.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model_rho2", "uhGM", "uhGM", "access-om3.mom6.3d.uhGM+rho2.1mon.mean.%4yr", "all", "average", "none", 2 - -"access-om3.mom6.3d.vhGM+rho2.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model_rho2", "vhGM", "vhGM", "access-om3.mom6.3d.vhGM+rho2.1mon.mean.%4yr", "all", "average", "none", 2 - -"access-om3.mom6.3d.GM_sfn_y+rho2.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model_rho2", "GM_sfn_y", "GM_sfn_y", "access-om3.mom6.3d.GM_sfn_y+rho2.1mon.mean.%4yr", "all", "average", "none", 2 - - -# monthly 3d fields on native - -"access-om3.mom6.3d.e.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "e", "e", "access-om3.mom6.3d.e.1mon.mean.%4yr", "all", "average", "none", 2 - -"access-om3.mom6.3d.thkcello.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "thkcello", "thkcello", "access-om3.mom6.3d.thkcello.1mon.mean.%4yr", "all", "average", "none", 2 - -"access-om3.mom6.3d.Kd_interface.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "Kd_interface", "Kd_interface", "access-om3.mom6.3d.Kd_interface.1mon.mean.%4yr", "all", "average", "none", 2 - -"access-om3.mom6.3d.Kd_ePBL.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "Kd_ePBL", "Kd_ePBL", "access-om3.mom6.3d.Kd_ePBL.1mon.mean.%4yr", "all", "average", "none", 2 - # monthly 2d fields @@ -169,41 +145,20 @@ ACCESS-OM3 "access-om3.mom6.2d.pbo.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" "ocean_model", "pbo", "pbo", "access-om3.mom6.2d.pbo.1mon.mean.%4yr", "all", "average", "none", 2 -"access-om3.mom6.2d.Chl_opac.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "Chl_opac", "Chl_opac", "access-om3.mom6.2d.Chl_opac.1mon.mean.%4yr", "all", "average", "none", 2 - -"access-om3.mom6.2d.SW_pen.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "SW_pen", "SW_pen", "access-om3.mom6.2d.SW_pen.1mon.mean.%4yr", "all", "average", "none", 2 - -"access-om3.mom6.2d.T_adx_2d.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "T_adx_2d", "T_adx_2d", "access-om3.mom6.2d.T_adx_2d.1mon.mean.%4yr", "all", "average", "none", 2 - -"access-om3.mom6.2d.T_ady_2d.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "T_ady_2d", "T_ady_2d", "access-om3.mom6.2d.T_ady_2d.1mon.mean.%4yr", "all", "average", "none", 2 - -"access-om3.mom6.2d.T_diffx_2d.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "T_diffx_2d", "T_diffx_2d", "access-om3.mom6.2d.T_diffx_2d.1mon.mean.%4yr", "all", "average", "none", 2 - -"access-om3.mom6.2d.T_diffy_2d.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "T_diffy_2d", "T_diffy_2d", "access-om3.mom6.2d.T_diffy_2d.1mon.mean.%4yr", "all", "average", "none", 2 - -"access-om3.mom6.2d.T_hbd_diffx_2d.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "T_hbd_diffx_2d", "T_hbd_diffx_2d", "access-om3.mom6.2d.T_hbd_diffx_2d.1mon.mean.%4yr", "all", "average", "none", 2 - -"access-om3.mom6.2d.T_hbd_diffy_2d.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "T_hbd_diffy_2d", "T_hbd_diffy_2d", "access-om3.mom6.2d.T_hbd_diffy_2d.1mon.mean.%4yr", "all", "average", "none", 2 - "access-om3.mom6.2d.friver.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" "ocean_model", "friver", "friver", "access-om3.mom6.2d.friver.1mon.mean.%4yr", "all", "average", "none", 2 "access-om3.mom6.2d.ficeberg.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" "ocean_model", "ficeberg", "ficeberg", "access-om3.mom6.2d.ficeberg.1mon.mean.%4yr", "all", "average", "none", 2 -"access-om3.mom6.2d.ePBL_h_ML.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "ePBL_h_ML", "ePBL_h_ML", "access-om3.mom6.2d.ePBL_h_ML.1mon.mean.%4yr", "all", "average", "none", 2 +"access-om3.mom6.2d.fsitherm.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" +"ocean_model", "fsitherm", "fsitherm", "access-om3.mom6.2d.fsitherm.1mon.mean.%4yr", "all", "average", "none", 2 + +"access-om3.mom6.2d.precip.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" +"ocean_model", "precip", "precip", "access-om3.mom6.2d.precip.1mon.mean.%4yr", "all", "average", "none", 2 -"access-om3.mom6.2d.h_ML.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "h_ML", "h_ML", "access-om3.mom6.2d.h_ML.1mon.mean.%4yr", "all", "average", "none", 2 +"access-om3.mom6.2d.frazil_heat_tendency_2d.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" +"ocean_model", "frazil_heat_tendency_2d", "frazil_heat_tendency_2d", "access-om3.mom6.2d.frazil_heat_tendency_2d.1mon.mean.%4yr", "all", "average", "none", 2 # monthly 2d max fields @@ -211,36 +166,6 @@ ACCESS-OM3 "access-om3.mom6.2d.mlotst.1mon.max.%4yr", 1, "months", 1, "days", "time", 1, "years" "ocean_model", "mlotst", "mlotst", "access-om3.mom6.2d.mlotst.1mon.max.%4yr", "all", "max", "none", 2 -"access-om3.mom6.2d.speed.1mon.max.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "speed", "speed", "access-om3.mom6.2d.speed.1mon.max.%4yr", "all", "max", "none", 2 - -"access-om3.mom6.2d.zos.1mon.max.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "zos", "zos", "access-om3.mom6.2d.zos.1mon.max.%4yr", "all", "max", "none", 2 - -"access-om3.mom6.2d.tos.1mon.max.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "tos", "tos", "access-om3.mom6.2d.tos.1mon.max.%4yr", "all", "max", "none", 2 - -"access-om3.mom6.2d.sos.1mon.max.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "sos", "sos", "access-om3.mom6.2d.sos.1mon.max.%4yr", "all", "max", "none", 2 - -"access-om3.mom6.2d.ePBL_h_ML.1mon.max.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "ePBL_h_ML", "ePBL_h_ML", "access-om3.mom6.2d.ePBL_h_ML.1mon.max.%4yr", "all", "max", "none", 2 - -"access-om3.mom6.2d.h_ML.1mon.max.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "h_ML", "h_ML", "access-om3.mom6.2d.h_ML.1mon.max.%4yr", "all", "max", "none", 2 - - -# monthly 2d min fields - -"access-om3.mom6.2d.zos.1mon.min.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "zos", "zos", "access-om3.mom6.2d.zos.1mon.min.%4yr", "all", "min", "none", 2 - -"access-om3.mom6.2d.tos.1mon.min.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "tos", "tos", "access-om3.mom6.2d.tos.1mon.min.%4yr", "all", "min", "none", 2 - -"access-om3.mom6.2d.sos.1mon.min.%4yr", 1, "months", 1, "days", "time", 1, "years" -"ocean_model", "sos", "sos", "access-om3.mom6.2d.sos.1mon.min.%4yr", "all", "min", "none", 2 - # daily 2D fields @@ -259,32 +184,43 @@ ACCESS-OM3 "access-om3.mom6.2d.speed.1day.mean.%4yr", 1, "days", 1, "days", "time", 1, "years" "ocean_model", "speed", "speed", "access-om3.mom6.2d.speed.1day.mean.%4yr", "all", "average", "none", 2 -"access-om3.mom6.2d.evs.1day.mean.%4yr", 1, "days", 1, "days", "time", 1, "years" -"ocean_model", "evs", "evs", "access-om3.mom6.2d.evs.1day.mean.%4yr", "all", "average", "none", 2 -"access-om3.mom6.2d.fsitherm.1day.mean.%4yr", 1, "days", 1, "days", "time", 1, "years" -"ocean_model", "fsitherm", "fsitherm", "access-om3.mom6.2d.fsitherm.1day.mean.%4yr", "all", "average", "none", 2 +# daily 2D ice shelf fields + +"access-om3.mom6.2d.melt.1day.mean.%4yr", 1, "days", 1, "days", "time", 1, "years" +"ice_shelf_model", "melt", "melt", "access-om3.mom6.2d.melt.1day.mean.%4yr", "all", "average", "none", 2 + +"access-om3.mom6.2d.thermal_driving.1day.mean.%4yr", 1, "days", 1, "days", "time", 1, "years" +"ice_shelf_model", "thermal_driving", "thermal_driving", "access-om3.mom6.2d.thermal_driving.1day.mean.%4yr", "all", "average", "none", 2 + +"access-om3.mom6.2d.ustar_shelf.1day.mean.%4yr", 1, "days", 1, "days", "time", 1, "years" +"ice_shelf_model", "ustar_shelf", "ustar_shelf", "access-om3.mom6.2d.ustar_shelf.1day.mean.%4yr", "all", "average", "none", 2 + + +# monthly 2d ice shelf fields + +"access-om3.mom6.2d.melt.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" +"ice_shelf_model", "melt", "melt", "access-om3.mom6.2d.melt.1mon.mean.%4yr", "all", "average", "none", 2 -"access-om3.mom6.2d.precip.1day.mean.%4yr", 1, "days", 1, "days", "time", 1, "years" -"ocean_model", "precip", "precip", "access-om3.mom6.2d.precip.1day.mean.%4yr", "all", "average", "none", 2 +"access-om3.mom6.2d.mass_flux.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" +"ice_shelf_model", "mass_flux", "mass_flux", "access-om3.mom6.2d.mass_flux.1mon.mean.%4yr", "all", "average", "none", 2 -"access-om3.mom6.2d.net_massout.1day.mean.%4yr", 1, "days", 1, "days", "time", 1, "years" -"ocean_model", "net_massout", "net_massout", "access-om3.mom6.2d.net_massout.1day.mean.%4yr", "all", "average", "none", 2 +"access-om3.mom6.2d.thermal_driving.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" +"ice_shelf_model", "thermal_driving", "thermal_driving", "access-om3.mom6.2d.thermal_driving.1mon.mean.%4yr", "all", "average", "none", 2 -"access-om3.mom6.2d.net_massin.1day.mean.%4yr", 1, "days", 1, "days", "time", 1, "years" -"ocean_model", "net_massin", "net_massin", "access-om3.mom6.2d.net_massin.1day.mean.%4yr", "all", "average", "none", 2 +"access-om3.mom6.2d.haline_driving.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" +"ice_shelf_model", "haline_driving", "haline_driving", "access-om3.mom6.2d.haline_driving.1mon.mean.%4yr", "all", "average", "none", 2 -"access-om3.mom6.2d.massout_flux.1day.mean.%4yr", 1, "days", 1, "days", "time", 1, "years" -"ocean_model", "massout_flux", "massout_flux", "access-om3.mom6.2d.massout_flux.1day.mean.%4yr", "all", "average", "none", 2 +"access-om3.mom6.2d.tfreeze.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" +"ice_shelf_model", "tfreeze", "tfreeze", "access-om3.mom6.2d.tfreeze.1mon.mean.%4yr", "all", "average", "none", 2 -"access-om3.mom6.2d.massin_flux.1day.mean.%4yr", 1, "days", 1, "days", "time", 1, "years" -"ocean_model", "massin_flux", "massin_flux", "access-om3.mom6.2d.massin_flux.1day.mean.%4yr", "all", "average", "none", 2 +"access-om3.mom6.2d.ustar_shelf.1mon.mean.%4yr", 1, "months", 1, "days", "time", 1, "years" +"ice_shelf_model", "ustar_shelf", "ustar_shelf", "access-om3.mom6.2d.ustar_shelf.1mon.mean.%4yr", "all", "average", "none", 2 # daily scalar timeseries "access-om3.mom6.scalar.1day.snap.%4yr", 1, "days", 1, "days", "time", 1, "years" -"ocean_model", "masso", "masso", "access-om3.mom6.scalar.1day.snap.%4yr", "all", "none", "none", 1 "ocean_model", "soga", "soga", "access-om3.mom6.scalar.1day.snap.%4yr", "all", "none", "none", 1 "ocean_model", "thetaoga", "thetaoga", "access-om3.mom6.scalar.1day.snap.%4yr", "all", "none", "none", 1 "ocean_model", "tosga", "tosga", "access-om3.mom6.scalar.1day.snap.%4yr", "all", "none", "none", 1 @@ -304,8 +240,3 @@ ACCESS-OM3 "ocean_model", "total_friver", "total_friver", "access-om3.mom6.scalar.1day.snap.%4yr", "all", "none", "none", 1 "ocean_model", "total_net_massout", "total_net_massout", "access-om3.mom6.scalar.1day.snap.%4yr", "all", "none", "none", 1 "ocean_model", "total_net_massin", "total_net_massin", "access-om3.mom6.scalar.1day.snap.%4yr", "all", "none", "none", 1 -"ocean_model", "total_PRCmE", "total_PRCmE", "access-om3.mom6.scalar.1day.snap.%4yr", "all", "none", "none", 1 -"ocean_model", "total_evap", "total_evap", "access-om3.mom6.scalar.1day.snap.%4yr", "all", "none", "none", 1 -"ocean_model", "total_icemelt", "total_icemelt", "access-om3.mom6.scalar.1day.snap.%4yr", "all", "none", "none", 1 -"ocean_model", "total_frunoff", "total_frunoff", "access-om3.mom6.scalar.1day.snap.%4yr", "all", "none", "none", 1 -"ocean_model", "total_lrunoff", "total_lrunoff", "access-om3.mom6.scalar.1day.snap.%4yr", "all", "none", "none", 1 diff --git a/diag_table_source.yaml b/diag_table_source.yaml index 59ff928a9..1af986949 100644 --- a/diag_table_source.yaml +++ b/diag_table_source.yaml @@ -173,7 +173,7 @@ diag_table: module_name: ocean_model_z file_name_diag_coord: "z" fields: - thkcello: # Cell Thickness + h: # Layer thickness agessc: # Sea water age since surface contact rhopot0: # Potential density referenced to surface rhopot2: # Potential density referenced to 2000 dbar @@ -182,8 +182,8 @@ diag_table: uo: # u velocity vo: # v velocity KE: # Layer Kinetic energy per unit mass - Kd_Itidal_Work: # Work done by Internal Tide Diapycnal Mixing - rsdoabsorb: # Convergence of Penetrative Shortwave Flux in Sea Water Layer + e: # Interface height relative to mean sea level + Kd_interface: #Total diapycnal diffusivity at interfaces 'monthly 3d fields on rho2': defaults: # these can be overridden for individual fields below @@ -196,19 +196,6 @@ diag_table: e: # Interface height relative to mean sea level umo: # Ocean Mass X Transport vmo: # Ocean Mass Y Transport - uhGM: # Time Mean Diffusive Zonal Thickness Flux - vhGM: # Time Mean Diffusive Meridional Thickness Flux - GM_sfn_y: # Parameterized Meridional Overturning Streamfunction - - 'monthly 3d fields on native': - defaults: # these can be overridden for individual fields below - file_name_dimension: 3d # descriptor for filename, e.g. 3d, 2d, scalar - reduction_method: mean - fields: - e: # Interface height relative to mean sea level - thkcello: # Cell Thickness - Kd_interface: # Total diapycnal diffusivity at interfaces - Kd_ePBL: # ePBL diapycnal diffusivity at interfaces 'monthly 2d fields': defaults: # these can be overridden for individual fields below @@ -232,18 +219,11 @@ diag_table: umo_2d: # Ocean Mass X Transport Vertical Sum vmo_2d: # Ocean Mass Y Transport Vertical Sum pbo: # Sea Water Pressure at Sea Floor - Chl_opac: # Surface chlorophyll A concentration used to find opacity - SW_pen: # Penetrating shortwave radiation flux into ocean - T_adx_2d: # Vertically Integrated Advective Zonal Flux of Heat - T_ady_2d: # Vertically Integrated Advective Meridional Flux of Heat - T_diffx_2d: # Vertically Integrated Diffusive Zonal Flux of Heat - T_diffy_2d: # Vertically Integrated Diffusive Meridional Flux of Heat - T_hbd_diffx_2d: # Vertically-integrated zonal diffusive flux from the horizontal boundary diffusion scheme for Heat - T_hbd_diffy_2d: # Vertically-integrated meridional diffusive flux from the horizontal boundary diffusion scheme for Heat friver: # Liquid runoff (rivers) into ocean ficeberg: # Frozen runoff (calving) and iceberg melt into ocean - ePBL_h_ML: # Surface boundary layer depth - h_ML: # Surface mixed layer depth based on active turbulence + fsitherm: # water flux to ocean from snow/sea ice melting(> 0) or formation(< 0) + precip: # Liquid+frozen precip into ocean + frazil_heat_tendency_2d: # Frazil heat tendency integrated vertically 'monthly 2d max fields': defaults: # these can be overridden for individual fields below @@ -251,21 +231,6 @@ diag_table: reduction_method: max # mean, snap, rms, pow##, min, max, or diurnal## fields: mlotst: # Mixed layer depth (delta rho = 0.03) MLD_003 - speed: # Sea Surface Speed - zos: # Sea surface height above geoid - tos: # Sea Surface Temperature - sos: # Sea Surface Salinity - ePBL_h_ML: # Surface boundary layer depth - h_ML: # Surface mixed layer depth based on active turbulence - - 'monthly 2d min fields': - defaults: # these can be overridden for individual fields below - file_name_dimension: 2d # descriptor for filename, e.g. 3d, 2d, scalar - reduction_method: min # mean, snap, rms, pow##, min, max, or diurnal## - fields: - zos: # Sea surface height above geoid - tos: # Sea Surface Temperature - sos: # Sea Surface Salinity 'daily 2D fields': defaults: @@ -280,13 +245,32 @@ diag_table: tos: # Sea Surface Temperature sos: # Sea Surface Salinity speed: # Sea Surface Speed - evs: # Evaporation/condensation at ocean surface (evaporation is negative) - fsitherm: # water flux to ocean from snow/sea ice melting(> 0) or formation(< 0) - precip: # Liquid + frozen precipitation into ocean - net_massout: # Net mass leaving the ocean due to evaporation, seaice formation - net_massin: # Net mass entering ocean due to precip, runoff, ice melt - massout_flux: # Net mass flux of freshwater out of the ocean (used in the boundary flux calculation) - massin_flux: # Net mass flux of freshwater into the ocean (used in boundary flux calculation) + + 'daily 2D ice shelf fields': + defaults: + file_name_dimension: 2d # descriptor for filename, e.g. 3d, 2d, scalar + module_name: ice_shelf_model + '': + - output_freq_units + output_freq_units: days # time units for output: years, months, days, hours, minutes, or seconds + reduction_method: mean # mean, snap, rms, pow##, min, max, or diurnal## + fields: + melt: # melt rate + thermal_driving: # thermal driving + ustar_shelf: #friction velocity + + 'monthly 2d ice shelf fields': + defaults: # these can be overridden for individual fields below + file_name_dimension: 2d # descriptor for filename, e.g. 3d, 2d, scalar + reduction_method: mean + module_name: ice_shelf_model + fields: + melt: # melt rate + mass_flux: # mass flux + thermal_driving: # thermal driving + haline_driving: # haline driving + tfreeze: # tfreeze + ustar_shelf: #friction velocity 'daily scalar timeseries': defaults: # these can be overridden for individual fields below @@ -314,7 +298,6 @@ diag_table: reduction_method: snap # mean, snap, rms, pow##, min, max, or diurnal## packing: 1 # double precision fields: - masso: # Mass of liquid ocean soga: # Global Mean Ocean Salinity thetaoga: # Global Mean Ocean Potential Temperature tosga: # Global Area Average Sea Surface Temperature sst_global @@ -333,10 +316,4 @@ diag_table: total_ficeberg: # Area integrated frozen runoff (calving) & iceberg melt into ocean total_frunoff total_friver: # Area integrated liquid runoff into ocean total_lrunoff total_net_massout: # Area integrated mass leaving ocean due to evap and seaice form - total_net_massin: # Area integrated mass entering ocean due to precip, runoff, ice melt - total_PRCmE: # Area integrated net surface water flux (precip+melt+liq runoff+ice calving-evap) - total_evap: # Area integrated evap/condense at ocean surface - total_icemelt: # Area integrated sea ice melt (>0) or form (<0) - total_precip: # Area integrated liquid+frozen precip into ocean - total_frunoff: # Area integrated frozen runoff (calving) & iceberg melt into ocean - total_lrunoff: # Area integrated liquid runoff into ocean + total_net_massin: # Area integrated mass entering ocean due to predip, runoff, ice melt diff --git a/docs/MOM_parameter_doc.all b/docs/MOM_parameter_doc.all index 858628026..86c4099bd 100644 --- a/docs/MOM_parameter_doc.all +++ b/docs/MOM_parameter_doc.all @@ -18,7 +18,7 @@ USE_EOS = True ! [Boolean] default = True 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 = True ! [Boolean] default = False +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 @@ -60,19 +60,11 @@ 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 periodcially 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 = 900.0 ! [s] +DT = 150.0 ! [s] ! The (baroclinic) dynamics time step. The time-step that is actually used will ! be an integer fraction of the forcing time-step (DT_FORCING in ocean-only mode ! or the coupling timestep in coupled mode.) -DT_THERM = 7200.0 ! [s] default = 900.0 +DT_THERM = 150.0 ! [s] default = 150.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 @@ -82,7 +74,7 @@ THERMO_SPANS_COUPLING = True ! [Boolean] default = False ! 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 = 7200.0 ! [s] default = 7200.0 +DT_TRACER_ADVECT = 150.0 ! [s] default = 150.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 @@ -93,19 +85,16 @@ TRADV_SPANS_COUPLING = True ! [Boolean] default = True ! 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 = True ! [Boolean] default = False +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. -THICKNESSDIFFUSE_FIRST = True ! [Boolean] default = False - ! If true, do thickness diffusion or interface height smoothing before dynamics. - ! This is only used if THICKNESSDIFFUSE or APPLY_INTERFACE_FILTER is true. -HMIX_SFC_PROP = 1.0 ! [m] default = 1.0 +HMIX_SFC_PROP = 1.08 ! [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 +HMIX_UV_SFC_PROP = 1.08 ! [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. @@ -117,7 +106,7 @@ HFREEZE = 10.0 ! [m] default = -1.0 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 = 0.0 ! [s] default = 7200.0 +DTBT_RESET_PERIOD = 0.0 ! [s] default = 150.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 @@ -142,7 +131,7 @@ C_P = 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. -USE_PSURF_IN_EOS = True ! [Boolean] default = True +USE_PSURF_IN_EOS = False ! [Boolean] default = True ! If true, always include the surface pressure contributions in equation of ! state calculations. P_REF = 2.0E+07 ! [Pa] default = 2.0E+07 @@ -158,23 +147,8 @@ ALTERNATE_FIRST_DIRECTION = False ! [Boolean] default = False ! 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 = True ! [Boolean] default = False +CHECK_BAD_SURFACE_VALS = False ! [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. -BAD_VAL_COLUMN_THICKNESS = 0.0 ! [m] default = 0.0 - ! The value of column thickness below which a bad value message is triggered, if - ! CHECK_BAD_SURFACE_VALS is true. DEFAULT_ANSWER_DATE = 99991231 ! default = 99991231 ! This sets the default value for the various _ANSWER_DATE parameters. SURFACE_ANSWER_DATE = 99991231 ! default = 99991231 @@ -214,13 +188,13 @@ REENTRANT_X = True ! [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 = True ! [Boolean] default = False +TRIPOLAR_N = False ! [Boolean] default = False ! Use tripolar connectivity at the northern edge of the domain. With ! TRIPOLAR_N, NIGLOBAL must be even. -NIGLOBAL = 1440 ! +NIGLOBAL = 4320 ! ! 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 = 1152 ! +NJGLOBAL = 1442 ! ! 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 @@ -307,8 +281,95 @@ MAXIMUM_DEPTH = 6000.0 ! [m] ! === 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 = 0 ! default = 0 +OBC_NUMBER_OF_SEGMENTS = 1 ! default = 0 ! The number of open boundary segments. +OBC_ZERO_VORTICITY = False ! [Boolean] default = False + ! If true, sets relative vorticity to zero on open boundaries. +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_SPECIFIED_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_STRAIN = False ! [Boolean] default = False + ! If true, sets the strain used in the stress tensor to zero on open boundaries. +OBC_FREESLIP_STRAIN = True ! [Boolean] default = True + ! If true, sets the normal gradient of tangential velocity to zero in the strain + ! use in the stress tensor on open boundaries. This cannot be true if another + ! OBC_XXX_STRAIN option is True. +OBC_COMPUTED_STRAIN = False ! [Boolean] default = False + ! If true, sets the normal gradient of tangential velocity to zero in the strain + ! use in the stress tensor on open boundaries. This cannot be true if another + ! OBC_XXX_STRAIN option is True. +OBC_SPECIFIED_STRAIN = False ! [Boolean] default = False + ! If true, sets the normal gradient of tangential velocity to zero in the strain + ! use in the stress tensor on open boundaries. This cannot be true if another + ! OBC_XXX_STRAIN option is True. +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, witha 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_SEGMENT_001 = "J=N,I=N:0,FLATHER,ORLANSKI,NUDGED" ! + ! Documentation needs to be dynamic????? +OBC_SEGMENT_001_VELOCITY_NUDGING_TIMESCALES = 0.3, 360.0 ! [days] + ! 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_RADIATION_MAX = 1.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 = 3.0E+04 ! [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 = 3000.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_REMAPPING_SCHEME = "PPM_H4" ! default = "PPM_H4" + ! 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 = True ! [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 = False ! [Boolean] default = False + ! 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. CHANNEL_CONFIG = "none" ! default = "none" ! A parameter that determines which set of channels are ! restricted to specific widths. Options are: @@ -345,7 +406,7 @@ RHO_0 = 1035.0 ! [kg m-3] default = 1035.0 ! 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 +ANGSTROM = 1.0E-15 ! [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. @@ -355,25 +416,27 @@ NK = 75 ! [nondim] ! === module MOM_tracer_registry === ! === module MOM_EOS === -EQN_OF_STATE = "ROQUET_RHO" ! default = "WRIGHT" +EQN_OF_STATE = "WRIGHT_FULL" ! default = "WRIGHT" ! EQN_OF_STATE determines which ocean equation of state should be used. ! Currently, the valid choices are "LINEAR", "UNESCO", "JACKETT_MCD", "WRIGHT", ! "WRIGHT_REDUCED", "WRIGHT_FULL", "NEMO", "ROQUET_RHO", "ROQUET_SPV" and ! "TEOS10". This is only used if USE_EOS is true. -EOS_QUADRATURE = True ! [Boolean] default = True +EOS_QUADRATURE = False ! [Boolean] default = False ! If true, always use the generic (quadrature) code code for the integrals of ! density. -TFREEZE_FORM = "TEOS_POLY" ! default = "TEOS10" +TFREEZE_FORM = "LINEAR" ! default = "LINEAR" ! TFREEZE_FORM determines which expression should be used for the freezing ! point. Currently, the valid choices are "LINEAR", "MILLERO_78", "TEOS_POLY", ! "TEOS10" -TFREEZE_S_IS_PRACS = False ! [Boolean] default = False - ! When True, the model will check if the model internal salinity is practical - ! salinity. If the model uses absolute salinity, a conversion will be applied. -TFREEZE_T_IS_POTT = False ! [Boolean] default = False - ! When True, the model will check if the model internal temperature is potential - ! temperature. If the model uses conservative temperature, a conversion will be - ! applied. +TFREEZE_S0_P0 = 0.0 ! [degC] default = 0.0 + ! When TFREEZE_FORM=LINEAR, this is the freezing potential temperature at S=0, + ! P=0. +DTFREEZE_DS = -0.054 ! [degC ppt-1] default = -0.054 + ! When TFREEZE_FORM=LINEAR, this is the derivative of the freezing potential + ! temperature with salinity. +DTFREEZE_DP = -7.75E-08 ! [degC 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 = True ! [Boolean] default = False @@ -463,6 +526,18 @@ TRACERS_MAY_REINIT = False ! [Boolean] default = False ! in the restart files. Otherwise it is a fatal error if the tracers are not ! found in the restart files of a restarted run. +! === module MOM_boundary_update === +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. + ! === module MOM_hor_index === ! Sets the horizontal array index types. @@ -497,7 +572,7 @@ GFS = 9.8 ! [m s-2] default = 9.8 ! The reduced gravity at the free surface. LIGHTEST_DENSITY = 1035.0 ! [kg m-3] default = 1035.0 ! The reference potential density used for layer 1. -REGRIDDING_COORDINATE_MODE = "ZSTAR" ! default = "LAYER" +REGRIDDING_COORDINATE_MODE = "SIGMA_SHELF_ZSTAR" ! default = "LAYER" ! Coordinate mode for vertical regridding. Choose among the following ! possibilities: LAYER - Isopycnal or stacked shallow water layers ! ZSTAR, Z* - stretched geopotential z* @@ -526,7 +601,7 @@ ALE_COORDINATE_CONFIG = "FILE:ocean_vgrid.nc,interfaces=zeta" ! default = "UNIFO ! by a comma or space, for sigma-2 and dz. e.g. ! HYBRID:vgrid.nc,sigma2,dz !ALE_RESOLUTION = 1.0825614929199219, 1.1963462829589844, 1.322089672088623, 1.4610481262207031, 1.614609718322754, 1.784308910369873, 1.9718408584594727, 2.1790781021118164, 2.4080896377563477, 2.661160469055176, 2.940814971923828, 3.249845504760742, 3.591329574584961, 3.968667984008789, 4.385614395141602, 4.846321105957031, 5.355350494384766, 5.917766571044922, 6.539115905761719, 7.225547790527344, 7.983818054199219, 8.821372985839844, 9.746376037597656, 10.767845153808594, 11.895652770996094, 13.140586853027344, 14.514511108398438, 16.030319213867188, 17.7020263671875, 19.544876098632812, 21.575271606445312, 23.810821533203125, 26.270294189453125, 28.973419189453125, 31.94091796875, 35.194000244140625, 38.75390625, 42.641632080078125, 46.876739501953125, 51.476593017578125, 56.45489501953125, 61.82025146484375, 67.5743408203125, 73.70965576171875, 80.207763671875, 87.03759765625, 94.1534423828125, 101.4951171875, 108.9879150390625, 116.5452880859375, 124.0714111328125, 131.4671630859375, 138.6346435546875, 145.484130859375, 151.938232421875, 157.93701171875, 163.439697265625, 168.42431640625, 172.8876953125, 176.842041015625, 180.3125, 183.33154296875, 185.938720703125, 188.175048828125, 190.08251953125, 191.701171875 -MIN_THICKNESS = 0.001 ! [m] default = 0.001 +MIN_THICKNESS = 1.0E-12 ! [m] default = 0.001 ! When regridding, this is the minimum layer thickness allowed. REMAPPING_SCHEME = "PPM_H4" ! default = "PLM" ! This sets the reconstruction scheme used for vertical remapping for all @@ -565,7 +640,7 @@ REMAP_BOUND_INTERMEDIATE_VALUES = False ! [Boolean] default = False REMAP_BOUNDARY_EXTRAP = False ! [Boolean] default = False ! If true, values at the interfaces of boundary cells are extrapolated instead ! of piecewise constant -INIT_BOUNDARY_EXTRAP = False ! [Boolean] default = False +INIT_BOUNDARY_EXTRAP = True ! [Boolean] default = False ! If true, values at the interfaces of boundary cells are extrapolated instead ! of piecewise constant during initialization.Defaults to REMAP_BOUNDARY_EXTRAP. REMAPPING_USE_OM4_SUBCELLS = False ! [Boolean] default = True @@ -573,11 +648,6 @@ REMAPPING_USE_OM4_SUBCELLS = False ! [Boolean] default = True ! reconstruction for the top- and lower-most sub-layers, but instead assumes ! they are always vanished (untrue) and so just uses their edge values. We ! recommend setting this option to false. -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. PARTIAL_CELL_VELOCITY_REMAP = False ! [Boolean] default = False ! If true, use partial cell thicknesses at velocity points that are masked out ! where they extend below the shallower of the neighboring bathymetry for @@ -617,15 +687,15 @@ INIT_LAYERS_FROM_Z_FILE = True ! [Boolean] default = False ! Z-space file on a latitude-longitude grid. ! === module MOM_initialize_layers_from_Z === -TEMP_SALT_Z_INIT_FILE = "woa23_ts_01_mom.nc" ! default = "temp_salt_z.nc" +TEMP_SALT_Z_INIT_FILE = "ocean_temp_salt.res.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 = "woa23_ts_01_mom.nc" ! default = "woa23_ts_01_mom.nc" +TEMP_Z_INIT_FILE = "ocean_temp_salt.res.nc" ! default = "ocean_temp_salt.res.nc" ! The name of the z-space input file used to initialize temperatures, only. -SALT_Z_INIT_FILE = "woa23_ts_01_mom.nc" ! default = "woa23_ts_01_mom.nc" +SALT_Z_INIT_FILE = "ocean_temp_salt.res.nc" ! default = "ocean_temp_salt.res.nc" ! The name of the z-space input file used to initialize temperatures, only. -Z_INIT_FILE_PTEMP_VAR = "ptemp" ! default = "ptemp" +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. @@ -672,10 +742,25 @@ HORIZ_INTERP_TOL_SALIN = 0.001 ! [ppt] default = 0.001 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 +TRIM_IC_FOR_P_SURF = True ! [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. +SURFACE_PRESSURE_FILE = "ice_shelf_thick.nc" ! + ! The initial condition file for the surface pressure exerted by ice. +SURFACE_PRESSURE_VAR = "thick" ! default = "" + ! The initial condition variable for the surface pressure exerted by ice. +SURFACE_PRESSURE_SCALE = 8986.6 ! [file dependent] default = 1.0 + ! A scaling factor to convert SURFACE_PRESSURE_VAR from file + ! SURFACE_PRESSURE_FILE into a surface pressure. +TRIM_IC_Z_TOLERANCE = 1.0E-10 ! [m] default = 1.0E-05 + ! The tolerance with which to find the depth matching the specified surface + ! pressure with TRIM_IC_FOR_P_SURF. +FRAC_DP_AT_POS_NEGATIVE_P_BUGFIX = True ! [Boolean] default = False + ! If true, use bugfix in ice shelf TRIM_IC initialization. Otherwise, pressure + ! input to density EOS is negative. +TRIMMING_USES_REMAPPING = True ! [Boolean] default = False + ! When trimming the column, also remap T and S. REGRID_ACCELERATE_INIT = False ! [Boolean] default = False ! If true, runs REGRID_ACCELERATE_ITERATIONS iterations of the regridding ! algorithm to push the initial grid to be consistent with the initial @@ -696,6 +781,18 @@ ODA_INCUPD = False ! [Boolean] default = False 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 via SPONGE_CONFIG. +OBC_SEGMENT_001_DATA = "U=file:OBC37S_forcing_access_yr2_4km.nc(u),V=file:OBC37S_forcing_access_yr2_4km.nc(v),SSH=file:OBC37S_forcing_access_yr2_4km.nc(eta_t),TEMP=file:OBC37S_forcing_access_yr2_4km.nc(pot_temp),SALT=file:OBC37S_forcing_access_yr2_4km.nc(salt)" ! + ! OBC segment docs +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 = 2 ! default = 1 @@ -756,7 +853,7 @@ REGRIDDING_ANSWER_DATE = 99991231 ! default = 20181231 ! 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. -DIAG_COORD_DEF_RHO2 = "RFNC1:76,999.5,1020.,1034.1,3.1,1041.,0.002" ! default = "WOA09" +DIAG_COORD_DEF_RHO2 = "FILE:diag_rho2.nc,interfaces=rho2" ! 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 @@ -781,156 +878,12 @@ DIAG_COORD_REGRID_COMPRESSIBILITY_FRACTION_RHO2 = 0.0 ! [nondim] default = 0.0 ! compressibility solely to make homogeneous regions appear stratified. ! === module MOM_MEKE === -USE_MEKE = True ! [Boolean] default = False +USE_MEKE = False ! [Boolean] default = False ! If true, turns on the MEKE scheme which calculates a sub-grid mesoscale eddy ! kinetic energy budget. -MEKE_IN_DYNAMICS = True ! [Boolean] default = True - ! If true, step MEKE forward with the dynamicsotherwise with the tracer - ! timestep. -EKE_SOURCE = "prog" ! default = "prog" - ! Determine the where EKE comes from: - ! 'prog': Calculated solving EKE equation - ! 'file': Read in from a file - ! 'dbclient': Retrieved from ML-database -MEKE_DAMPING = 0.0 ! [s-1] default = 0.0 - ! The local depth-independent MEKE dissipation rate. -MEKE_CD_SCALE = 0.0 ! [nondim] default = 0.0 - ! The ratio of the bottom eddy velocity to the column mean eddy velocity, i.e. - ! sqrt(2*MEKE). This should be less than 1 to account for the surface - ! intensification of MEKE. -MEKE_CB = 25.0 ! [nondim] default = 25.0 - ! A coefficient in the expression for the ratio of bottom projected eddy energy - ! and mean column energy (see Jansen et al. 2015). -MEKE_MIN_GAMMA2 = 1.0E-04 ! [nondim] default = 1.0E-04 - ! The minimum allowed value of gamma_b^2. -MEKE_CT = 50.0 ! [nondim] default = 50.0 - ! A coefficient in the expression for the ratio of barotropic eddy energy and - ! mean column energy (see Jansen et al. 2015). -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_GEOMETRIC = False ! [Boolean] default = False - ! If MEKE_GEOMETRIC is true, uses the GM coefficient formulation from the - ! GEOMETRIC framework (Marshall et al., 2012). -MEKE_GEOMETRIC_ALPHA = 0.05 ! [nondim] default = 0.05 - ! The nondimensional coefficient governing the efficiency of the GEOMETRIC - ! thickness diffusion. -MEKE_EQUILIBRIUM_ALT = False ! [Boolean] default = False - ! If true, use an alternative formula for computing the (equilibrium)initial - ! value of MEKE. -MEKE_EQUILIBRIUM_RESTORING = False ! [Boolean] default = False - ! If true, restore MEKE back to its equilibrium value, which is calculated at - ! each time step. -MEKE_FRCOEFF = -1.0 ! [nondim] default = -1.0 - ! The efficiency of the conversion of mean energy into MEKE. If MEKE_FRCOEFF is - ! negative, this conversion is not used or calculated. -MEKE_BHFRCOEFF = -1.0 ! [nondim] default = -1.0 - ! The efficiency of the conversion of mean energy into MEKE by the biharmonic - ! dissipation. If MEKE_bhFRCOEFF is negative, this conversion is not used or - ! calculated. -MEKE_GMECOEFF = -1.0 ! [nondim] default = -1.0 - ! The efficiency of the conversion of MEKE into mean energy by GME. If - ! MEKE_GMECOEFF 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_KH = -1.0 ! [m2 s-1] default = -1.0 - ! A background lateral diffusivity of MEKE. Use a negative value to not apply - ! lateral diffusion to MEKE. -MEKE_K4 = -1.0 ! [m4 s-1] default = -1.0 - ! A lateral bi-harmonic diffusivity of MEKE. Use a negative value to not apply - ! bi-harmonic diffusion to MEKE. -MEKE_DTSCALE = 1.0 ! [nondim] default = 1.0 - ! A scaling factor to accelerate the time evolution of MEKE. -MEKE_POSITIVE = False ! [Boolean] default = False - ! If true, it guarantees that MEKE will always be >= 0. -MEKE_KHCOEFF = 1.0 ! [nondim] default = 1.0 - ! A scaling factor in the expression for eddy diffusivity which is otherwise - ! proportional to the MEKE velocity- scale times an eddy mixing-length. This - ! factor must be >0 for MEKE to contribute to the thickness/ and tracer - ! diffusivity in the rest of the model. -MEKE_USCALE = 0.0 ! [m s-1] default = 0.0 - ! The background velocity that is combined with MEKE to calculate the bottom - ! drag. -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 MEKE GM source term. -MEKE_VISC_DRAG = True ! [Boolean] default = True - ! If true, use the vertvisc_type to calculate the bottom drag acting on MEKE. -MEKE_KHTH_FAC = 0.5 ! [nondim] default = 0.0 - ! A factor that maps MEKE%Kh to KhTh. -MEKE_KHTR_FAC = 0.5 ! [nondim] default = 0.0 - ! A factor that maps MEKE%Kh to KhTr. -MEKE_KHMEKE_FAC = 1.0 ! [nondim] default = 0.0 - ! A factor that maps MEKE%Kh to Kh for MEKE itself. -MEKE_OLD_LSCALE = False ! [Boolean] default = False - ! If true, use the old formula for length scale which is a function of grid - ! spacing and deformation radius. -MEKE_MIN_LSCALE = False ! [Boolean] default = False - ! If true, use a strict minimum of provided length scales rather than harmonic - ! mean. -MEKE_RD_MAX_SCALE = False ! [Boolean] default = False - ! If true, the length scale used by MEKE is the minimum of the deformation - ! radius or grid-spacing. Only used if MEKE_OLD_LSCALE=True -MEKE_VISCOSITY_COEFF_KU = 1.0 ! [nondim] default = 0.0 - ! If non-zero, is the scaling coefficient in the expression forviscosity used to - ! parameterize harmonic lateral momentum mixing byunresolved eddies represented - ! by MEKE. Can be negative torepresent backscatter from the unresolved eddies. -MEKE_VISCOSITY_COEFF_AU = 0.0 ! [nondim] default = 0.0 - ! If non-zero, is the scaling coefficient in the expression forviscosity used to - ! parameterize biharmonic lateral momentum mixing byunresolved eddies - ! represented by MEKE. Can be negative torepresent backscatter from the - ! unresolved eddies. -MEKE_FIXED_MIXING_LENGTH = 0.0 ! [m] default = 0.0 - ! If positive, is a fixed length contribution to the expression for mixing - ! length used in MEKE-derived diffusivity. -MEKE_FIXED_TOTAL_DEPTH = True ! [Boolean] default = True - ! If true, use the nominal bathymetric depth as the estimate of the time-varying - ! ocean depth. Otherwise base the depth on the total ocean massper unit area. -MEKE_ALPHA_DEFORM = 0.0 ! [nondim] default = 0.0 - ! If positive, is a coefficient weighting the deformation scale in the - ! expression for mixing length used in MEKE-derived diffusivity. -MEKE_ALPHA_RHINES = 0.15 ! [nondim] default = 0.0 - ! 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.0 - ! If positive, is a coefficient weighting the Eady length scale in the - ! expression for mixing length used in MEKE-derived diffusivity. -MEKE_ALPHA_FRICT = 0.0 ! [nondim] default = 0.0 - ! If positive, is a coefficient weighting the frictional arrest scale in the - ! expression for mixing length used in MEKE-derived diffusivity. -MEKE_ALPHA_GRID = 0.0 ! [nondim] default = 0.0 - ! If positive, is a coefficient weighting the grid-spacing as a scale in the - ! expression for mixing length used in MEKE-derived diffusivity. -MEKE_COLD_START = False ! [Boolean] default = False - ! If true, initialize EKE to zero. Otherwise a local equilibrium solution is - ! used as an initial condition for EKE. -MEKE_BACKSCAT_RO_C = 0.0 ! [nondim] default = 0.0 - ! The coefficient in the Rossby number function for scaling the biharmonic - ! frictional energy source. Setting to non-zero enables the Rossby number - ! function. -MEKE_BACKSCAT_RO_POW = 0.0 ! [nondim] default = 0.0 - ! The power in the Rossby number function for scaling the biharmonic frictional - ! energy source. -MEKE_ADVECTION_FACTOR = 0.0 ! [nondim] default = 0.0 - ! A scale factor in front of advection of eddy energy. Zero turns advection off. - ! Using unity would be normal but other values could accommodate a mismatch - ! between the advecting barotropic flow and the vertical structure of MEKE. -MEKE_TOPOGRAPHIC_BETA = 0.0 ! [nondim] default = 0.0 - ! A scale factor to determine how much topographic beta is weighed in computing - ! beta in the expression of Rhines scale. Use 1 if full topographic beta effect - ! is considered; use 0 if it's completely ignored. -SQG_USE_MEKE = False ! [Boolean] default = False - ! If true, the eddy scale of MEKE is used for the SQG vertical structure -CDRAG = 0.003 ! [nondim] default = 0.003 - ! CDRAG is the drag coefficient relating the magnitude of the velocity field to - ! the bottom stress. -MEKE_CDRAG = 0.003 ! [nondim] default = 0.003 - ! Drag coefficient relating the magnitude of the velocity field to the bottom - ! stress in MEKE. ! === module MOM_lateral_mixing_coeffs === -USE_VARIABLE_MIXING = True ! [Boolean] default = False +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 @@ -938,7 +891,7 @@ USE_VARIABLE_MIXING = True ! [Boolean] default = False USE_VISBECK = False ! [Boolean] default = False ! If true, use the Visbeck et al. (1997) formulation for ! thickness diffusivity. -RESOLN_SCALED_KH = True ! [Boolean] default = False +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 @@ -946,7 +899,7 @@ DEPTH_SCALED_KHTH = False ! [Boolean] default = False ! 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 = True ! [Boolean] default = False +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 @@ -961,7 +914,10 @@ BACKSCAT_EBT_POWER = 0.0 ! [nondim] default = 0.0 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 = True ! [Boolean] default = False +SQG_EXPO = 1.0 ! [nondim] default = 1.0 + ! Nondimensional exponent coeffecient of the SQG mode that is used for the + ! vertical struture of diffusivities. +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 @@ -985,7 +941,7 @@ KHTH_SLOPE_CFF = 0.0 ! [nondim] default = 0.0 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 +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. @@ -995,14 +951,6 @@ VERY_SMALL_FREQUENCY = 1.0E-17 ! [s-1] default = 1.0E-17 USE_STANLEY_ISO = False ! [Boolean] default = False ! If true, turn on Stanley SGS T variance parameterization in isopycnal slope ! code. -RESOLN_N2_FILTER_DEPTH = -1.0 ! [m] default = -1.0 - ! The depth below which N2 is monotonized to avoid stratification artifacts from - ! altering the equivalent barotropic mode structure. This monotonzization is - ! disabled if this parameter is negative. -VISBECK_MAX_SLOPE = 0.0 ! [nondim] default = 0.0 - ! If non-zero, is an upper bound on slopes used in the Visbeck formula for - ! diffusivity. This does not affect the isopycnal slope calculation used within - ! thickness diffusion. 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. @@ -1012,45 +960,18 @@ USE_SIMPLER_EADY_GROWTH_RATE = False ! [Boolean] default = False VARMIX_KTOP = 2 ! [nondim] default = 2 ! The layer number at which to start vertical integration of S*N for purposes of ! finding the Eady growth rate. +MIN_DZ_FOR_SLOPE_N2 = 1.0 ! [m] default = 1.0 + ! The minimum vertical distance to use in the denominator of the bouyancy + ! frequency used in the slope calculation. +FULL_DEPTH_EADY_GROWTH_RATE = False ! [Boolean] default = False + ! If true, calculate the Eady growth rate based on average slope times + ! stratification that includes contributions from sea-level changes in its + ! denominator, rather than just the nominal depth of the bathymetry. This only + ! applies when using the model interface heights as a proxy for isopycnal + ! slopes. VISBECK_L_SCALE = 0.0 ! [m or nondim] default = 0.0 ! The fixed length scale in the Visbeck formula, or if negative a nondimensional ! scaling factor relating this length scale squared to the cell areas. -KH_RES_SCALE_COEF = 1.0 ! [nondim] default = 1.0 - ! A coefficient that determines how KhTh is scaled away if RESOLN_SCALED_... is - ! true, as F = 1 / (1 + (KH_RES_SCALE_COEF*Rd/dx)^KH_RES_FN_POWER). -KH_RES_FN_POWER = 100 ! default = 2 - ! The power of dx/Ld in the Kh resolution function. Any positive integer may be - ! used, although even integers are more efficient to calculate. Setting this - ! greater than 100 results in a step-function being used. -VISC_RES_SCALE_COEF = 1.0 ! [nondim] default = 1.0 - ! A coefficient that determines how Kh is scaled away if RESOLN_SCALED_... is - ! true, as F = 1 / (1 + (KH_RES_SCALE_COEF*Rd/dx)^KH_RES_FN_POWER). This - ! function affects lateral viscosity, Kh, and not KhTh. -VISC_RES_FN_POWER = 100 ! default = 100 - ! The power of dx/Ld in the Kh resolution function. Any positive integer may be - ! used, although even integers are more efficient to calculate. Setting this - ! greater than 100 results in a step-function being used. This function affects - ! lateral viscosity, Kh, and not KhTh. -INTERPOLATE_RES_FN = False ! [Boolean] default = False - ! 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 = True - ! 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. -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. -EBT_REMAPPING_USE_OM4_SUBCELLS = False ! [Boolean] default = False - ! 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. USE_QG_LEITH_GM = False ! [Boolean] default = False ! If true, use the QG Leith viscosity as the GM coefficient. @@ -1071,7 +992,7 @@ DRAG_AS_BODY_FORCE = False ! [Boolean] default = False 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. -LINEAR_DRAG = False ! [Boolean] default = False +LINEAR_DRAG = True ! [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 @@ -1082,11 +1003,14 @@ HBBL = 10.0 ! [m] ! 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.003 ! [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. 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.1 ! [m s-1] default = 0.0 +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 @@ -1106,12 +1030,12 @@ HTBL_SHELF_MIN = 0.1 ! [m] default = 0.1 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 = 0.0 ! [m2 s-1] +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 = 0.0 +KV_BBL_MIN = 1.0E-06 ! [m2 s-1] default = 1.0E-06 ! The minimum viscosities in the bottom boundary layer. -KV_TBL_MIN = 0.0 ! [m2 s-1] default = 0.0 +KV_TBL_MIN = 1.0E-06 ! [m2 s-1] default = 1.0E-06 ! 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 @@ -1160,23 +1084,17 @@ DETANGLE_INTERFACES = False ! [Boolean] default = False 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 = True ! [Boolean] default = False +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. -FGNV_FILTER_SCALE = 0.1 ! [nondim] default = 1.0 - ! A coefficient scaling the vertical smoothing term in the Ferrari et al., 2010, - ! streamfunction formulation. -FGNV_C_MIN = 0.0 ! [m s-1] default = 0.0 - ! A minium wave speed used in the Ferrari et al., 2010, streamfunction - ! formulation. -FGNV_STRAT_FLOOR = 1.0E-15 ! [nondim] default = 1.0E-15 - ! A floor for Brunt-Vasaila frequency in the Ferrari et al., 2010, - ! streamfunction formulation, expressed as a fraction of planetary rotation, - ! OMEGA. This should be tiny but non-zero to avoid degeneracy. USE_STANLEY_GM = False ! [Boolean] default = False ! If true, turn on Stanley SGS T variance parameterization in GM code. -USE_KH_IN_MEKE = False ! [Boolean] default = False - ! If true, uses the thickness diffusivity calculated here to diffuse MEKE. +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. @@ -1242,7 +1160,7 @@ 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 = 1.0E-06 ! [m] default = 3.75E-09 +ETA_TOLERANCE = 1.0E-06 ! [m] default = 3.75E-14 ! 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 @@ -1327,21 +1245,24 @@ SSH_IN_EOS_PRESSURE_FOR_PGF = False ! [Boolean] default = False MASS_WEIGHT_IN_PRESSURE_GRADIENT = True ! [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 +MASS_WEIGHT_IN_PRESSURE_GRADIENT_TOP = True ! [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 +MASS_WEIGHT_IN_PGF_VANISHED_ONLY = True ! [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. CORRECTION_INTXPA = False ! [Boolean] default = False ! If true, use a correction for surface pressure curvature in intx_pa. -RESET_INTXPA_INTEGRAL = False ! [Boolean] default = False +RESET_INTXPA_INTEGRAL = True ! [Boolean] default = False ! If true, reset INTXPA to match pressures at first nonvanished cell. Includes ! pressure correction. -RESET_INTXPA_INTEGRAL_FLATTEST = False ! [Boolean] default = False +RESET_INTXPA_INTEGRAL_FLATTEST = True ! [Boolean] default = False ! If true, use flattest interface as reference interface where there is no ! better choice for RESET_INTXPA_INTEGRAL. Otherwise, use surface interface. +RESET_INTXPA_H_NONVANISHED = 1.0E-06 ! [m] default = 1.0E-06 + ! A minimal layer thickness that indicates that a layer is thick enough to + ! usefully reestimate the pressure integral across the interface below. 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. @@ -1349,7 +1270,7 @@ RECONSTRUCT_FOR_PRESSURE = True ! [Boolean] default = True ! 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 +PRESSURE_RECONSTRUCTION_SCHEME = 2 ! 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. @@ -1381,19 +1302,13 @@ KH_VEL_SCALE = 0.01 ! [m s-1] default = 0.0 ! The velocity scale which is multiplied by the grid spacing to calculate the ! Laplacian viscosity. The final viscosity is the largest of this scaled ! viscosity, the Smagorinsky and Leith viscosities, and KH. -KH_SIN_LAT = 2000.0 ! [m2 s-1] default = 0.0 +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). -KH_PWR_OF_SINE = 4.0 ! [nondim] default = 4.0 - ! The power used to raise SIN(LAT) when using a latitudinally dependent - ! background viscosity. SMAGORINSKY_KH = False ! [Boolean] default = False ! If true, use a Smagorinsky nonlinear eddy viscosity. LEITH_KH = False ! [Boolean] default = False ! If true, use a Leith nonlinear eddy viscosity. -RES_SCALE_MEKE_VISC = False ! [Boolean] default = False - ! If true, the viscosity contribution from MEKE is scaled by the resolution - ! function. BOUND_KH = True ! [Boolean] default = True ! If true, the Laplacian coefficient is locally limited to be stable. BETTER_BOUND_KH = True ! [Boolean] default = True @@ -1496,7 +1411,7 @@ VON_KARMAN_CONST = 0.41 ! [nondim] default = 0.41 HARMONIC_VISC = False ! [Boolean] default = False ! If true, use the harmonic mean thicknesses for calculating the vertical ! viscosity. -HARMONIC_BL_SCALE = 0.0 ! [nondim] default = 0.0 +HARMONIC_BL_SCALE = 1.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. @@ -1739,6 +1654,13 @@ DIAG_EBT_MONO_N2_COLUMN_FRACTION = 0.0 ! [nondim] default = 0.0 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 = False ! [Boolean] default = False ! If true, use the OM4 remapping-via-subcells algorithm for calculating EBT ! structure. See REMAPPING_USE_OM4_SUBCELLS for details. We recommend setting @@ -1786,7 +1708,7 @@ KD_BBL_TR = 0.0 ! [m2 s-1] default = 0.0 ! 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 +MINIMUM_FORCING_DEPTH = 1.08 ! [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 @@ -1815,6 +1737,9 @@ DIAG_DEPTH_SUBML_N2 = 50.0 ! [m] default = 50.0 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. +FRAZIL_NOT_UNDER_ICESHELF = True ! [Boolean] default = False + ! If true, do not use frazil scheme underneath ice shelves defined by + ! frac_shelf_h greater than 0. ! === module MOM_CVMix_conv === ! Parameterization of enhanced mixing due to convection via CVMix @@ -1846,7 +1771,7 @@ TIDAL_MIXING_ANSWER_DATE = 99991231 ! default = 99991231 ! The vintage of the order of arithmetic and expressions in the tidal mixing ! 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. -INT_TIDE_PROFILE = "POLZIN_09" ! default = "STLAURENT_02" +INT_TIDE_PROFILE = "STLAURENT_02" ! 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 @@ -1862,26 +1787,7 @@ INT_TIDE_LOWMODE_DISSIPATION = False ! [Boolean] default = False ! generated; as with itidal drag on the barotropic tide, use an internal tidal ! dissipation scheme to drive diapycnal mixing, along the lines of St. Laurent ! et al. (2002) and Simmons et al. (2004). -NU_POLZIN = 0.0697 ! [nondim] default = 0.0697 - ! When the Polzin decay profile is used, this is a non-dimensional constant in - ! the expression for the vertical scale of decay for the tidal energy - ! dissipation. -NBOTREF_POLZIN = 9.61E-04 ! [s-1] default = 9.61E-04 - ! When the Polzin decay profile is used, this is the reference value of the - ! buoyancy frequency at the ocean bottom in the Polzin formulation for the - ! vertical scale of decay for the tidal energy dissipation. -POLZIN_DECAY_SCALE_FACTOR = 1.0 ! [nondim] default = 1.0 - ! When the Polzin decay profile is used, this is a scale factor for the vertical - ! scale of decay of the tidal energy dissipation. -POLZIN_SCALE_MAX_FACTOR = 1.0 ! [nondim] default = 1.0 - ! When the Polzin decay profile is used, this is a factor to limit the vertical - ! scale of decay of the tidal energy dissipation to - ! POLZIN_DECAY_SCALE_MAX_FACTOR times the depth of the ocean. -POLZIN_MIN_DECAY_SCALE = 0.0 ! [m] default = 0.0 - ! When the Polzin decay profile is used, this is the minimum vertical decay - ! scale for the vertical profile - ! of internal tide dissipation with the Polzin (2009) formulation -INT_TIDE_DECAY_SCALE = 300.3003003003003 ! [m] default = 500.0 +INT_TIDE_DECAY_SCALE = 500.0 ! [m] default = 500.0 ! The decay scale away from the bottom for tidal TKE with the new coding when ! INT_TIDE_DISSIPATION is used. MU_ITIDES = 0.2 ! [nondim] default = 0.2 @@ -1893,14 +1799,14 @@ GAMMA_ITIDES = 0.3333 ! [nondim] default = 0.3333 MIN_ZBOT_ITIDES = 0.0 ! [m] default = 0.0 ! Turn off internal tidal dissipation when the total ocean depth is less than ! this value. -KAPPA_ITIDES = 6.28319E-04 ! [m-1] default = 6.283185307179586E-04 +KAPPA_ITIDES = 6.283185307179586E-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. UTIDE = 0.0 ! [m s-1] default = 0.0 ! The constant tidal amplitude used with INT_TIDE_DISSIPATION. -KAPPA_H2_FACTOR = 0.84 ! [nondim] default = 1.0 +KAPPA_H2_FACTOR = 1.0 ! [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 +TKE_ITIDE_MAX = 1000.0 ! [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 @@ -1951,7 +1857,7 @@ LOTW_BBL_ANSWER_DATE = 99991231 ! default = 20190101 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 = True ! [Boolean] default = False +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. @@ -2156,16 +2062,9 @@ USE_CALVING_HEAT_CONTENT = True ! [Boolean] default = False ! instead of using SST*CP*froz_runoff. DO_BRINE_PLUME = False ! [Boolean] default = False ! If true, use a brine plume parameterization from Nguyen et al., 2009. -VAR_PEN_SW = True ! [Boolean] default = False +VAR_PEN_SW = False ! [Boolean] default = False ! If true, use one of the CHL_A schemes specified by OPACITY_SCHEME to determine ! the e-folding depth of incoming short wave radiation. -CHL_FROM_FILE = True ! [Boolean] default = True - ! If true, chl_a is read from a file. -CHL_FILE = "chl_globcolour_monthly_clim.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. -CHL_VARNAME = "CHL" ! default = "CHL_A" - ! Name of CHL_A variable in CHL_FILE. ! === module MOM_energetic_PBL === ML_OMEGA_FRAC = 0.001 ! [nondim] default = 0.0 @@ -2301,19 +2200,11 @@ LT_MOD_LAC4 = 0.0 ! [nondim] default = 0.95 LT_MOD_LAC5 = 0.22 ! [nondim] default = 0.95 ! Coefficient for modification of Langmuir number due to ratio of Ekman to ! unstable Obukhov depth. -EPBL_EQD_DIFFUSIVITY_SHAPE = False ! [nondim] default = False - ! Logical flag for activating ML equation for shape function that uses forcing - ! to change its structure. -EPBL_EQD_DIFFUSIVITY_VELOCITY = False ! [nondim] default = False - ! Logical flag for activating ML equation discovery for velocity scale -EPBL_EQD_DIFFUSIVITY_VELOCITY_H = False ! [nondim] default = False - ! Logical flag for activating ML equation discovery for velocity scale with h as - ! input EPBL_OPTIONS_DIFF = 0 ! default = 0 ! If positive, this is a coded integer indicating a pair of settings whose ! differences are diagnosed in a passive diagnostic mode via extra calls to ! ePBL_column. If this is 0 or negative no extra calls occur. -!EPBL_USTAR_MIN = 1.45842E-18 ! [m s-1] +!EPBL_USTAR_MIN = 1.45842E-23 ! [m s-1] ! The (tiny) minimum friction velocity used within the ePBL code, derived from ! OMEGA and ANGSTROM. @@ -2323,15 +2214,17 @@ REGULARIZE_SURFACE_LAYERS = False ! [Boolean] default = False ! much lateral variations to allow for sensible lateral barotropic transports. ! === module MOM_opacity === -OPACITY_SCHEME = "MANIZZA_05" ! default = "MANIZZA_05" - ! This character string specifies how chlorophyll concentrations are translated - ! into opacities. Currently valid options include: - ! MANIZZA_05 - Use Manizza et al., GRL, 2005. - ! MOREL_88 - Use Morel, JGR, 1988. - ! OHLMANN_03 - Use Ohlmann, J Clim, 2003. (only use if dz(1)>2.0m) -BLUE_FRAC_SW = 0.5 ! [nondim] default = 0.5 - ! The fraction of the penetrating shortwave radiation that is in the blue band. -PEN_SW_NBANDS = 3 ! default = 1 +EXP_OPACITY_SCHEME = "SINGLE_EXP" ! default = "SINGLE_EXP" + ! This character string specifies which exponential opacity scheme to utilize. + ! Currently valid options include: + ! SINGLE_EXP - Single Exponent decay. + ! DOUBLE_EXP - Double Exponent decay. +PEN_SW_SCALE = 15.0 ! [m] default = 0.0 + ! The vertical absorption e-folding depth of the penetrating shortwave + ! radiation. +PEN_SW_FRAC = 0.42 ! [nondim] default = 0.0 + ! The fraction of the shortwave radiation that penetrates below the surface. +PEN_SW_NBANDS = 1 ! default = 1 ! The number of bands of penetrating shortwave radiation. OPTICS_ANSWER_DATE = 99991231 ! default = 99991231 ! The vintage of the order of arithmetic and expressions in the optics @@ -2345,18 +2238,6 @@ PEN_SW_FLUX_ABSORB = 2.5E-11 ! [degC m s-1] default = 2.5E-11 PEN_SW_ABSORB_MINTHICK = 1.0 ! [m] default = 1.0 ! A thickness that is used to absorb the remaining penetrating shortwave heat ! flux when it drops below PEN_SW_FLUX_ABSORB. -OPACITY_VALUES_MANIZZA = 0.0232, 0.074, 0.225, 0.037, 2.86, 0.0 ! [m-1] default = 0.0232, 0.074, 0.225, 0.037, 2.86, 0.0 - ! Pairs of opacity coefficients for blue, red and near-infrared radiation with - ! parameterizations following the functional form from Manizza et al., GRL 2005, - ! namely in the form opacity = coef_1 + coef_2 * chl**pow for each band. - ! Although coefficients are set for 3 bands, more or less bands may actually be - ! used, with extra bands following the same properties as band 3. -CHOROPHYLL_POWER_MANIZZA = 0.674, 0.629, 0.0 ! [nondim] default = 0.674, 0.629, 0.0 - ! Powers of chlorophyll for blue, red and near-infrared radiation bands in - ! expressions for opacity of the form opacity = coef_1 + coef_2 * chl**pow. -OPACITY_BAND_WAVELENGTHS = 0.0, 550.0, 700.0, 2800.0 ! [nm] default = 0.0, 550.0, 700.0, 2800.0 - ! The bounding wavelengths for the various bands of shortwave radiation, with - ! defaults that depend on the setting for OPACITY_SCHEME. OPACITY_LAND_VALUE = 10.0 ! [m-1] default = 10.0 ! The value to use for opacity over land. The default is 10 m-1 - a value for ! muddy water. @@ -2373,13 +2254,12 @@ USE_HUYNH_STENCIL_BUG = False ! [Boolean] default = False ! required to reproduce results in legacy simulations. ! === module MOM_tracer_hor_diff === -KHTR = 50.0 ! [m2 s-1] default = 0.0 +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. Set to a positive value to - ! activate. + ! 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 @@ -2405,28 +2285,8 @@ RECALC_NEUTRAL_SURF = False ! [Boolean] default = False ! === module MOM_neutral_diffusion === ! This module implements neutral diffusion of tracers -USE_NEUTRAL_DIFFUSION = True ! [Boolean] default = False +USE_NEUTRAL_DIFFUSION = False ! [Boolean] default = False ! If true, enables the neutral diffusion module. -NDIFF_CONTINUOUS = True ! [Boolean] default = True - ! If true, uses a continuous reconstruction of T and S when finding neutral - ! surfaces along which diffusion will happen. If false, a PPM discontinuous - ! reconstruction of T and S is done which results in a higher order routine but - ! exacts a higher computational cost. -NDIFF_REF_PRES = -1.0 ! [Pa] default = -1.0 - ! The reference pressure (Pa) used for the derivatives of the equation of state. - ! If negative (default), local pressure is used. -NDIFF_INTERIOR_ONLY = False ! [Boolean] default = False - ! If true, only applies neutral diffusion in the ocean interior.That is, the - ! algorithm will exclude the surface and bottomboundary layers. -NDIFF_USE_UNMASKED_TRANSPORT_BUG = False ! [Boolean] default = False - ! If true, use an older form for the accumulation of neutral-diffusion - ! transports that were unmasked, as used prior to Jan 2018. This is not - ! recommended. -NDIFF_ANSWER_DATE = 20240101 ! default = 20240101 - ! The vintage of the order of arithmetic to use for the neutral diffusion. - ! Values of 20240330 or below recover the answers from the original form of the - ! neutral diffusion code, while higher values use mathematically equivalent - ! expressions that recover rotational symmetry. ! === module MOM_hor_bnd_diffusion === ! This module implements horizontal diffusion of tracers near boundaries @@ -2444,7 +2304,7 @@ CALCULATE_APE = True ! [Boolean] default = True WRITE_STOCKS = True ! [Boolean] default = True ! If true, write the integrated tracer amounts to stdout when the energy files ! are written. -MAXTRUNC = 10000 ! [truncations save_interval-1] default = 0 +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. @@ -2456,8 +2316,6 @@ ENERGYFILE = "ocean.stats" ! default = "ocean.stats" ! The file to use to write the energies and globally summed diagnostics. DATE_STAMPED_STDOUT = True ! [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 = 8.64E+04 ! [s] default = 8.64E+04 ! The time unit in seconds a number of input fields READ_DEPTH_LIST = False ! [Boolean] default = False @@ -2509,7 +2367,7 @@ RESTORE_SALINITY = True ! [Boolean] default = False RESTORE_TEMPERATURE = False ! [Boolean] default = False ! If true, the coupled driver will add a heat flux that drives sea-surface ! temperature toward specified values. -ICE_SHELF = False ! [Boolean] default = False +ICE_SHELF = True ! [Boolean] default = False ! If true, enables the ice shelf model. ICEBERGS_APPLY_RIGID_BOUNDARY = False ! [Boolean] default = False ! If true, allows icebergs to change boundary condition felt by ocean @@ -2521,23 +2379,21 @@ LATENT_HEAT_FUSION = 3.337E+05 ! [J/kg] default = 3.34E+05 ! The latent heat of fusion. LATENT_HEAT_VAPORIZATION = 2.501E+06 ! [J/kg] default = 2.5E+06 ! The latent heat of fusion. -MAX_P_SURF = -1.0 ! [Pa] default = -1.0 +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. -ADJUST_NET_SRESTORE_TO_ZERO = True ! [Boolean] default = True +ADJUST_NET_SRESTORE_TO_ZERO = False ! [Boolean] default = True ! If true, adjusts the salinity restoring seen to zero whether restoring is via ! a salt flux or virtual precip. ADJUST_NET_SRESTORE_BY_SCALING = False ! [Boolean] default = False ! If true, adjustments to salt restoring to achieve zero net are made by scaling ! values without moving the zero contour. -ADJUST_NET_FRESH_WATER_TO_ZERO = True ! [Boolean] default = False +ADJUST_NET_FRESH_WATER_TO_ZERO = False ! [Boolean] default = False ! If true, adjusts the net fresh-water forcing seen by the ocean (including ! restoring) to zero. -USE_NET_FW_ADJUSTMENT_SIGN_BUG = False ! [Boolean] default = False - ! If true, use the wrong sign for the adjustment to the net fresh-water. ADJUST_NET_FRESH_WATER_BY_SCALING = False ! [Boolean] default = False ! If true, adjustments to net fresh water to achieve zero net are made by ! scaling values without moving the zero contour. @@ -2562,7 +2418,7 @@ FLUXCONST = 0.11 ! [m day-1] default = 0.0 ! anomalies (akin to a piston velocity). Note the non-MKS units. SALT_RESTORE_FILE = "salt_sfc_restore.nc" ! default = "salt_restore.nc" ! A file in which to find the surface salinity to use for restoring. -SALT_RESTORE_VARIABLE = "asalt" ! default = "salt" +SALT_RESTORE_VARIABLE = "salt" ! default = "salt" ! The name of the surface salinity variable to read from SALT_RESTORE_FILE for ! restoring salinity. SRESTORE_AS_SFLUX = True ! [Boolean] default = False @@ -2578,7 +2434,7 @@ MASK_SRESTORE_MARGINAL_SEAS = False ! [Boolean] default = False ! RESTORE_SALINITY is True. BASIN_FILE = "basin.nc" ! default = "basin.nc" ! A file in which to find the basin masks, in variable 'basin'. -MASK_SRESTORE = False ! [Boolean] default = False +MASK_SRESTORE = True ! [Boolean] default = False ! If true, read a file (salt_restore_mask) containing a mask for SSS restoring. CD_TIDES = 1.0E-04 ! [nondim] default = 1.0E-04 ! The drag coefficient that applies to the tides. @@ -2614,6 +2470,122 @@ LIQUID_RUNOFF_FROM_DATA = False ! [Boolean] default = False ! the component name 'OCN'. ! === module MOM_restart === + +! === module MOM_domains min_halo === + +! === module MOM_grid === +! Parameters providing information about the lateral grid. + +! === module MOM_hor_index === +! Sets the horizontal array index types. + +! === module MOM_grid_init === +MOM_ICESHELF_AVAILABLE_DIAGS_FILE = "MOM_IceShelf.available_diags" ! default = "MOM_IceShelf.available_diags" + ! A file into which to write a list of all available ice shelf diagnostics that + ! can be included in a diag_table. + +! === module MOM_ice_shelf === +DEBUG_IS = False ! [Boolean] default = False + ! If true, write verbose debugging messages for the ice shelf. +DYNAMIC_SHELF_MASS = False ! [Boolean] default = False + ! If true, the ice sheet mass can evolve with time. +DATA_OVERRIDE_MELT = False ! [Boolean] default = False + ! If true, the data override feature is used for melt rate and heat flux from + ! ice shelves. +SHELF_THERMO = True ! [Boolean] default = False + ! If true, use a thermodynamically interactive ice shelf. +SHELF_THREE_EQN = True ! [Boolean] default = True + ! If true, use the three equation expression of consistency to calculate the + ! fluxes at the ice-ocean interface. +SHELF_INSULATOR = False ! [Boolean] default = False + ! If true, the ice shelf is a perfect insulator (no conduction). +MELTING_CUTOFF_DEPTH = 0.0 ! [m] default = 0.0 + ! Depth above which the melt is set to zero (it must be >= 0) Default value + ! won't affect the solution. +CONST_SEA_LEVEL = False ! [Boolean] default = False + ! If true, apply evaporative, heat and salt fluxes in the sponge region. This + ! will avoid a large increase in sea level. This option is needed for some of + ! the ISOMIP+ experiments (Ocean3 and Ocean4). IMPORTANT: it is not currently + ! possible to do prefect restarts using this flag. +CONST_SEA_LEVEL_MISOMIP = False ! [Boolean] default = False + ! If true, constant_sea_level fluxes are applied only over the surface sponge + ! cells from the ISOMIP/MISOMIP configuration +SHELF_3EQ_GAMMA = False ! [Boolean] default = False + ! If true, user specifies a constant nondimensional heat-transfer coefficient + ! (GAMMA_T_3EQ), from which the default salt-transfer coefficient is set as + ! GAMMA_T_3EQ/35. This is used with SHELF_THREE_EQN. +SHELF_S_ROOT = False ! [Boolean] default = False + ! If SHELF_S_ROOT = True, salinity at the ice/ocean interface (Sbdry) is + ! computed from a quadratic equation. Otherwise, the previous interactive method + ! to estimate Sbdry is used. +ICE_SHELF_MASS_FROM_FILE = False ! [Boolean] default = False + ! Read the mass of the ice shelf (every time step) from a file. +C_P_ICE = 2100.0 ! [J kg-1 K-1] default = 2100.0 + ! The heat capacity of ice. +ICE_SHELF_FLUX_FACTOR = 1.0 ! [none] default = 1.0 + ! Non-dimensional factor applied to shelf thermodynamic fluxes. +KV_ICE = 1.0E+10 ! [m2 s-1] default = 1.0E+10 + ! The viscosity of the ice. +ICE_SHELF_SALINITY = 0.0 ! [psu] default = 0.0 + ! The salinity of the ice inside the ice shelf. +ICE_SHELF_TEMPERATURE = -20.0 ! [degC] default = -15.0 + ! The temperature at the center of the ice shelf. +KD_SALT_MOLECULAR = 8.02E-10 ! [m2 s-1] default = 8.02E-10 + ! The molecular diffusivity of salt in sea water at the freezing point. +KD_TEMP_MOLECULAR = 1.41E-07 ! [m2 s-1] default = 1.41E-07 + ! The molecular diffusivity of heat in sea water at the freezing point. +DT_FORCING = 0.0 ! [s] default = 0.0 + ! The time step for changing forcing, coupling with other components, or + ! potentially writing certain diagnostics. The default value is given by DT. +COL_THICK_MELT_THRESHOLD = 0.001 ! [m] default = 0.0 + ! The minimum ocean column thickness where melting is allowed. +ICE_SHELF_LINEAR_SHELF_FRAC = 0.13 ! [nondim] default = 0.13 + ! Ratio of HJ99 stability constant xi_N (ratio of maximum mixing length to + ! planetary boundary layer depth in neutrally stable conditions) to the von + ! Karman constant +ICE_SHELF_VK_CNST = 0.4 ! [nondim] default = 0.4 + ! Von Karman constant. +ICE_SHELF_RC = 0.2 ! [nondim] default = 0.2 + ! Critical flux Richardson number for ice melt +ICE_SHELF_USTAR_FROM_VEL_BUGFIX = True ! [Boolean] default = False + ! Bug fix for ice-area weighting of squared ocean velocities used to calculate + ! friction velocity under ice shelves +ICE_SHELF_BUOYANCY_FLUX_ITT_BUGFIX = True ! [Boolean] default = True + ! Bug fix of buoyancy iteration +ICE_SHELF_SALT_FLUX_ITT_BUGFIX = True ! [Boolean] default = True + ! Bug fix of salt iteration +ICE_SHELF_BUOYANCY_FLUX_ITT_THRESHOLD = 1.0E-04 ! [nondim] default = 1.0E-04 + ! Convergence criterion of Newton's method for ice shelf buoyancy iteration. +ICE_SHELF_TIDEAMP_SCALING_FACTOR = 0.66 ! [nondim] default = 1.0 + ! Scale TIDEAMP_FILE or UTIDE by number in melt parameterisation friction + ! velocity calculation. + +! === module MOM_EOS === +DENSITY_ICE = 917.0 ! [kg m-3] default = 900.0 + ! A typical density of ice. +MIN_THICKNESS_SIMPLE_CALVE = 0.0 ! [m] default = 0.0 + ! Min thickness rule for the very simple calving law +USTAR_SHELF_BG = 6.0E-04 ! [m s-1] default = 0.0 + ! The minimum value of ustar under ice shelves. +CDRAG_SHELF = 0.007 ! [nondim] default = 0.003 + ! CDRAG is the drag coefficient relating the magnitude of the velocity field to + ! the surface stress. +USTAR_SHELF_FROM_VEL = True ! [Boolean] default = True + ! If true, use the surface velocities to set the friction velocity under ice + ! shelves instead of using the previous values of the stresses. + +! === module MOM_restart === +ICE_PROFILE_CONFIG = "FILE" ! + ! This specifies how the initial ice profile is specified. Valid values are: + ! CHANNEL, FILE, and USER. +ICE_THICKNESS_FILE = "ice_shelf_thick.nc" ! default = "ice_shelf_h.nc" + ! The file from which the bathymetry is read. +LEN_SIDE_STRESS = 0.0 ! [axis_units] default = 0.0 + ! position past which shelf sides are stress free. +ICE_THICKNESS_VARNAME = "thick" ! default = "h_shelf" + ! The name of the thickness variable in ICE_THICKNESS_FILE. +ICE_AREA_VARNAME = "area" ! default = "area_shelf_h" + ! The name of the area variable in ICE_THICKNESS_FILE. WAVE_INTERFACE_ANSWER_DATE = 20221231 ! default = 20221231 ! The vintage of the order of arithmetic and expressions in the surface wave ! calculations. Values below 20230101 recover the answers from the end of 2022, diff --git a/docs/MOM_parameter_doc.debugging b/docs/MOM_parameter_doc.debugging index 75ee5e538..e1e65e29c 100644 --- a/docs/MOM_parameter_doc.debugging +++ b/docs/MOM_parameter_doc.debugging @@ -84,8 +84,11 @@ WRITE_TRACER_MIN_MAX = False ! [Boolean] default = False ! If true, write the maximum and minimum values of temperature, salinity and ! some tracer concentrations to stdout when the energy files are written. +! === module MOM_unit_scaling === +! Parameters for doing unit scaling of variables. + ! === module MOM_file_parser === REPORT_UNUSED_PARAMS = True ! [Boolean] default = True ! If true, report any parameter lines that are not used in the run. -FATAL_UNUSED_PARAMS = True ! [Boolean] default = False +FATAL_UNUSED_PARAMS = False ! [Boolean] default = False ! If true, kill the run if there are any unused parameters. diff --git a/docs/MOM_parameter_doc.layout b/docs/MOM_parameter_doc.layout index 6523517fe..97c60ace9 100644 --- a/docs/MOM_parameter_doc.layout +++ b/docs/MOM_parameter_doc.layout @@ -23,9 +23,9 @@ THIN_HALO_UPDATES = True ! [Boolean] default = True ! statically determined at compile time. Otherwise the sizes are not determined ! until run time. The STATIC option is substantially faster, but does not allow ! the PE count to be changed at run time. This can only be set at compile time. -AUTO_MASKTABLE = True ! [Boolean] default = False +AUTO_MASKTABLE = False ! [Boolean] default = False ! Turn on automatic mask table generation to eliminate land blocks. -MASKTABLE = "MOM_auto_mask_table" ! default = "MOM_mask_table" +MASKTABLE = "mask_table.1490.84x49" ! default = "MOM_mask_table" ! A text file to specify n_mask, layout and mask_list. This feature masks out ! processors that contain only land points. The first line of mask_table is the ! number of regions to be masked out. The second line is the layout of the model @@ -38,19 +38,15 @@ MASKTABLE = "MOM_auto_mask_table" ! default = "MOM_mask_table" ! 4,6 ! 1,2 ! 3,6 -NIPROC = 48 ! +NIPROC = 84 ! ! The number of processors in the x-direction. With STATIC_MEMORY_ this is set ! in MOM_memory.h at compile time. -NJPROC = 48 ! +NJPROC = 49 ! ! The number of processors in the y-direction. With STATIC_MEMORY_ this is set ! in MOM_memory.h at compile time. -LAYOUT = 48, 48 ! +LAYOUT = 84, 49 ! ! The processor layout that was actually used. -AUTO_IO_LAYOUT_FAC = 6 ! default = 0 - ! When AUTO_MASKTABLE is enabled, io layout is calculated by performing integer - ! division of the runtime-determined domain layout with this factor. If the - ! factor is set to 0 (default), the io layout is set to 1,1. -IO_LAYOUT = 8, 8 ! +IO_LAYOUT = 12, 7 ! default = 1, 1 ! The processor layout to be used, or 0,0 to automatically set the io_layout to ! be the same as the layout. @@ -67,11 +63,12 @@ BT_USE_WIDE_HALOS = True ! [Boolean] default = True ! efficiency. BTHALO = 0 ! default = 0 ! The minimum halo size for the barotropic solver. -BT_WIDE_HALO_MIN_STENCIL = 0 ! default = 0 - ! The minimum stencil width to use with the wide halo iterations. A nonzero - ! value may be useful for debugging purposes, but at the cost of reducing the - ! efficiency gain from BT_USE_WIDE_HALOS. !BT x-halo = 4 ! ! The barotropic x-halo size that is actually used. !BT y-halo = 4 ! ! The barotropic y-halo size that is actually used. + +! === module MOM_domains min_halo === + +! === module MOM_grid === +! Parameters providing information about the lateral grid. diff --git a/docs/MOM_parameter_doc.short b/docs/MOM_parameter_doc.short index 5f09abe7f..ae820140a 100644 --- a/docs/MOM_parameter_doc.short +++ b/docs/MOM_parameter_doc.short @@ -1,39 +1,32 @@ ! This file was written by the model and records the non-default parameters used at run-time. ! === module MOM === -USE_CONTEMP_ABSSAL = True ! [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. USE_REGRIDDING = True ! [Boolean] default = False ! If True, use the ALE algorithm (regridding/remapping). If False, use the ! layered isopycnal algorithm. -DT = 900.0 ! [s] +DT = 150.0 ! [s] ! The (baroclinic) dynamics time step. The time-step that is actually used will ! be an integer fraction of the forcing time-step (DT_FORCING in ocean-only mode ! or the coupling timestep in coupled mode.) -DT_THERM = 7200.0 ! [s] default = 900.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 = True ! [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. -THICKNESSDIFFUSE = True ! [Boolean] default = False - ! If true, isopycnal surfaces are diffused with a Laplacian coefficient of KHTH. -THICKNESSDIFFUSE_FIRST = True ! [Boolean] default = False - ! If true, do thickness diffusion or interface height smoothing before dynamics. - ! This is only used if THICKNESSDIFFUSE or APPLY_INTERFACE_FILTER is true. +HMIX_SFC_PROP = 1.08 ! [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 = 1.08 ! [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 = 10.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. -DTBT_RESET_PERIOD = 0.0 ! [s] default = 7200.0 +DTBT_RESET_PERIOD = 0.0 ! [s] default = 150.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 @@ -49,20 +42,9 @@ C_P = 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. +USE_PSURF_IN_EOS = False ! [Boolean] default = True + ! If true, always include the surface pressure contributions in equation of + ! state calculations. WRITE_GEOM = 2 ! 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 @@ -71,13 +53,10 @@ GEOM_FILE = "access-om3.mom6.geometry.nc" ! default = "ocean_geometry" ! The file into which to write the ocean geometry. ! === module MOM_domains === -TRIPOLAR_N = True ! [Boolean] default = False - ! Use tripolar connectivity at the northern edge of the domain. With - ! TRIPOLAR_N, NIGLOBAL must be even. -NIGLOBAL = 1440 ! +NIGLOBAL = 4320 ! ! 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 = 1152 ! +NJGLOBAL = 1442 ! ! 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. @@ -128,21 +107,45 @@ TOPO_CONFIG = "file" ! MAXIMUM_DEPTH = 6000.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 = 1 ! default = 0 + ! The number of open boundary segments. +OBC_ZERO_BIHARMONIC = True ! [Boolean] default = False + ! If true, zeros the Laplacian of flow on open boundaries in the biharmonic + ! viscosity term. +OBC_SEGMENT_001 = "J=N,I=N:0,FLATHER,ORLANSKI,NUDGED" ! + ! Documentation needs to be dynamic????? +OBC_SEGMENT_001_VELOCITY_NUDGING_TIMESCALES = 0.3, 360.0 ! [days] + ! 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_TRACER_RESERVOIR_LENGTH_SCALE_OUT = 3.0E+04 ! [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 = 3000.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. +BRUSHCUTTER_MODE = True ! [Boolean] default = False + ! If true, read external OBC data on the supergrid. + ! === module MOM_verticalGrid === ! Parameters providing information about the vertical grid. +ANGSTROM = 1.0E-15 ! [m] default = 1.0E-10 + ! The minimum layer thickness, usually one-Angstrom. NK = 75 ! [nondim] ! The number of model layers. ! === module MOM_EOS === -EQN_OF_STATE = "ROQUET_RHO" ! default = "WRIGHT" +EQN_OF_STATE = "WRIGHT_FULL" ! default = "WRIGHT" ! EQN_OF_STATE determines which ocean equation of state should be used. ! Currently, the valid choices are "LINEAR", "UNESCO", "JACKETT_MCD", "WRIGHT", ! "WRIGHT_REDUCED", "WRIGHT_FULL", "NEMO", "ROQUET_RHO", "ROQUET_SPV" and ! "TEOS10". This is only used if USE_EOS is true. -TFREEZE_FORM = "TEOS_POLY" ! default = "TEOS10" - ! TFREEZE_FORM determines which expression should be used for the freezing - ! point. Currently, the valid choices are "LINEAR", "MILLERO_78", "TEOS_POLY", - ! "TEOS10" +DTFREEZE_DP = -7.75E-08 ! [degC 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 = True ! [Boolean] default = False @@ -156,12 +159,18 @@ USE_IDEAL_AGE_TRACER = True ! [Boolean] default = False ! === module ideal_age_example === +! === module MOM_boundary_update === + ! === module MOM_fixed_initialization === ! === module MOM_grid_init === +! === module MOM_open_boundary === +! Controls where open boundaries are located, what kind of boundary condition to impose, and what data to apply, +! if any. + ! === module MOM_coord_initialization === -REGRIDDING_COORDINATE_MODE = "ZSTAR" ! default = "LAYER" +REGRIDDING_COORDINATE_MODE = "SIGMA_SHELF_ZSTAR" ! default = "LAYER" ! Coordinate mode for vertical regridding. Choose among the following ! possibilities: LAYER - Isopycnal or stacked shallow water layers ! ZSTAR, Z* - stretched geopotential z* @@ -188,6 +197,8 @@ ALE_COORDINATE_CONFIG = "FILE:ocean_vgrid.nc,interfaces=zeta" ! default = "UNIFO ! by a comma or space, for sigma-2 and dz. e.g. ! HYBRID:vgrid.nc,sigma2,dz !ALE_RESOLUTION = 1.0825614929199219, 1.1963462829589844, 1.322089672088623, 1.4610481262207031, 1.614609718322754, 1.784308910369873, 1.9718408584594727, 2.1790781021118164, 2.4080896377563477, 2.661160469055176, 2.940814971923828, 3.249845504760742, 3.591329574584961, 3.968667984008789, 4.385614395141602, 4.846321105957031, 5.355350494384766, 5.917766571044922, 6.539115905761719, 7.225547790527344, 7.983818054199219, 8.821372985839844, 9.746376037597656, 10.767845153808594, 11.895652770996094, 13.140586853027344, 14.514511108398438, 16.030319213867188, 17.7020263671875, 19.544876098632812, 21.575271606445312, 23.810821533203125, 26.270294189453125, 28.973419189453125, 31.94091796875, 35.194000244140625, 38.75390625, 42.641632080078125, 46.876739501953125, 51.476593017578125, 56.45489501953125, 61.82025146484375, 67.5743408203125, 73.70965576171875, 80.207763671875, 87.03759765625, 94.1534423828125, 101.4951171875, 108.9879150390625, 116.5452880859375, 124.0714111328125, 131.4671630859375, 138.6346435546875, 145.484130859375, 151.938232421875, 157.93701171875, 163.439697265625, 168.42431640625, 172.8876953125, 176.842041015625, 180.3125, 183.33154296875, 185.938720703125, 188.175048828125, 190.08251953125, 191.701171875 +MIN_THICKNESS = 1.0E-12 ! [m] default = 0.001 + ! When regridding, this is the minimum layer thickness allowed. 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: @@ -200,6 +211,9 @@ REMAPPING_SCHEME = "PPM_H4" ! default = "PLM" ! WENO_HYBGEN (3rd-order accurate) ! PQM_IH4IH3 (4th-order accurate) ! PQM_IH6IH5 (5th-order accurate) +INIT_BOUNDARY_EXTRAP = True ! [Boolean] default = False + ! If true, values at the interfaces of boundary cells are extrapolated instead + ! of piecewise constant during initialization.Defaults to REMAP_BOUNDARY_EXTRAP. REMAPPING_USE_OM4_SUBCELLS = False ! [Boolean] default = True ! This selects the remapping algorithm used in OM4 that does not use the full ! reconstruction for the top- and lower-most sub-layers, but instead assumes @@ -212,10 +226,12 @@ INIT_LAYERS_FROM_Z_FILE = True ! [Boolean] default = False ! Z-space file on a latitude-longitude grid. ! === module MOM_initialize_layers_from_Z === -TEMP_SALT_Z_INIT_FILE = "woa23_ts_01_mom.nc" ! default = "temp_salt_z.nc" +TEMP_SALT_Z_INIT_FILE = "ocean_temp_salt.res.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. +Z_INIT_FILE_PTEMP_VAR = "temp" ! default = "ptemp" + ! The name of the potential temperature variable in TEMP_Z_INIT_FILE. Z_INIT_ALE_REMAPPING = True ! [Boolean] default = False ! If True, then remap straight to model coordinate from file. Z_INIT_REMAP_GENERAL = True ! [Boolean] default = False @@ -225,6 +241,27 @@ TEMP_SALT_INIT_VERTICAL_REMAP_ONLY = True ! [Boolean] default = False ! If true, initial conditions are on the model horizontal grid. Extrapolation ! over missing ocean values is done using an ICE-9 procedure with vertical ALE ! remapping . +TRIM_IC_FOR_P_SURF = True ! [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. +SURFACE_PRESSURE_FILE = "ice_shelf_thick.nc" ! + ! The initial condition file for the surface pressure exerted by ice. +SURFACE_PRESSURE_VAR = "thick" ! default = "" + ! The initial condition variable for the surface pressure exerted by ice. +SURFACE_PRESSURE_SCALE = 8986.6 ! [file dependent] default = 1.0 + ! A scaling factor to convert SURFACE_PRESSURE_VAR from file + ! SURFACE_PRESSURE_FILE into a surface pressure. +TRIM_IC_Z_TOLERANCE = 1.0E-10 ! [m] default = 1.0E-05 + ! The tolerance with which to find the depth matching the specified surface + ! pressure with TRIM_IC_FOR_P_SURF. +FRAC_DP_AT_POS_NEGATIVE_P_BUGFIX = True ! [Boolean] default = False + ! If true, use bugfix in ice shelf TRIM_IC initialization. Otherwise, pressure + ! input to density EOS is negative. +TRIMMING_USES_REMAPPING = True ! [Boolean] default = False + ! When trimming the column, also remap T and S. +OBC_SEGMENT_001_DATA = "U=file:OBC37S_forcing_access_yr2_4km.nc(u),V=file:OBC37S_forcing_access_yr2_4km.nc(v),SSH=file:OBC37S_forcing_access_yr2_4km.nc(eta_t),TEMP=file:OBC37S_forcing_access_yr2_4km.nc(pot_temp),SALT=file:OBC37S_forcing_access_yr2_4km.nc(salt)" ! + ! OBC segment docs ! === module MOM_diag_mediator === NUM_DIAG_COORDS = 2 ! default = 1 @@ -255,7 +292,7 @@ REGRIDDING_ANSWER_DATE = 99991231 ! default = 20181231 ! 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. -DIAG_COORD_DEF_RHO2 = "RFNC1:76,999.5,1020.,1034.1,3.1,1041.,0.002" ! default = "WOA09" +DIAG_COORD_DEF_RHO2 = "FILE:diag_rho2.nc,interfaces=rho2" ! 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 @@ -272,70 +309,24 @@ DIAG_COORD_DEF_RHO2 = "RFNC1:76,999.5,1020.,1034.1,3.1,1041.,0.002" ! default = ! 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_KHTH_FAC = 0.5 ! [nondim] default = 0.0 - ! A factor that maps MEKE%Kh to KhTh. -MEKE_KHTR_FAC = 0.5 ! [nondim] default = 0.0 - ! A factor that maps MEKE%Kh to KhTr. -MEKE_KHMEKE_FAC = 1.0 ! [nondim] default = 0.0 - ! A factor that maps MEKE%Kh to Kh for MEKE itself. -MEKE_VISCOSITY_COEFF_KU = 1.0 ! [nondim] default = 0.0 - ! If non-zero, is the scaling coefficient in the expression forviscosity used to - ! parameterize harmonic lateral momentum mixing byunresolved eddies represented - ! by MEKE. Can be negative torepresent backscatter from the unresolved eddies. -MEKE_ALPHA_RHINES = 0.15 ! [nondim] default = 0.0 - ! 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.0 - ! 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. -KHTH_USE_EBT_STRUCT = True ! [Boolean] default = False - ! If true, uses the equivalent barotropic structure as the vertical structure of - ! thickness diffusivity. 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. -KH_RES_FN_POWER = 100 ! default = 2 - ! The power of dx/Ld in the Kh resolution function. Any positive integer may be - ! used, although even integers are more efficient to calculate. Setting this - ! greater than 100 results in a step-function being used. ! === 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. +LINEAR_DRAG = True ! [Boolean] default = False + ! If LINEAR_DRAG and BOTTOMDRAGLAW are defined the drag law is + ! cdrag*DRAG_BG_VEL*u. 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. -DRAG_BG_VEL = 0.1 ! [m s-1] default = 0.0 +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 @@ -344,7 +335,7 @@ 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 = 0.0 ! [m2 s-1] +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. @@ -354,12 +345,6 @@ KHTH_MAX_CFL = 0.1 ! [nondimensional] default = 0.8 ! 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. -KHTH_USE_FGNV_STREAMFUNCTION = True ! [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. -FGNV_FILTER_SCALE = 0.1 ! [nondim] default = 1.0 - ! A coefficient scaling the vertical smoothing term in the Ferrari et al., 2010, - ! streamfunction formulation. ! === module MOM_dynamics_split_RK2 === VISC_REM_BUG = False ! [Boolean] default = True @@ -369,7 +354,7 @@ VISC_REM_BUG = False ! [Boolean] default = True ! VISC_REM_BT_WEIGHT_BUG in MOM_barotropic. ! === module MOM_continuity_PPM === -ETA_TOLERANCE = 1.0E-06 ! [m] default = 3.75E-09 +ETA_TOLERANCE = 1.0E-06 ! [m] default = 3.75E-14 ! 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 @@ -401,6 +386,25 @@ BOUND_CORIOLIS = True ! [Boolean] default = False MASS_WEIGHT_IN_PRESSURE_GRADIENT = True ! [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 = True ! [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 = True ! [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. +RESET_INTXPA_INTEGRAL = True ! [Boolean] default = False + ! If true, reset INTXPA to match pressures at first nonvanished cell. Includes + ! pressure correction. +RESET_INTXPA_INTEGRAL_FLATTEST = True ! [Boolean] default = False + ! If true, use flattest interface as reference interface where there is no + ! better choice for RESET_INTXPA_INTEGRAL. Otherwise, use surface interface. +PRESSURE_RECONSTRUCTION_SCHEME = 2 ! 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. ! === module MOM_Zanna_Bolton === @@ -411,9 +415,6 @@ KH_VEL_SCALE = 0.01 ! [m s-1] default = 0.0 ! The velocity scale which is multiplied by the grid spacing to calculate the ! Laplacian viscosity. The final viscosity is the largest of this scaled ! viscosity, the Smagorinsky and Leith viscosities, and KH. -KH_SIN_LAT = 2000.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). AH_VEL_SCALE = 0.01 ! [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 @@ -435,6 +436,10 @@ FRICTWORK_BUG = False ! [Boolean] default = True ! recommended. ! === module MOM_vert_friction === +HARMONIC_BL_SCALE = 1.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 = 0.5 ! [m] ! The prescribed depth over which the near-surface viscosity and diffusivity are ! elevated when the bulk mixed layer is not used. @@ -514,6 +519,14 @@ ENERGETICS_SFC_PBL = True ! [Boolean] default = False 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. +MINIMUM_FORCING_DEPTH = 1.08 ! [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. +FRAZIL_NOT_UNDER_ICESHELF = True ! [Boolean] default = False + ! If true, do not use frazil scheme underneath ice shelves defined by + ! frac_shelf_h greater than 0. ! === module MOM_set_diffusivity === @@ -522,24 +535,6 @@ EPBL_IS_ADDITIVE = False ! [Boolean] default = True 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. @@ -561,10 +556,6 @@ LOTW_BBL_ANSWER_DATE = 99991231 ! default = 20190101 ! calculations. Values below 20240630 recover the original answers, while ! higher values use more accurate expressions. This only applies when ! USE_LOTW_BBL_DIFFUSIVITY is true. -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 @@ -616,14 +607,6 @@ USE_RIVER_HEAT_CONTENT = True ! [Boolean] default = False USE_CALVING_HEAT_CONTENT = True ! [Boolean] default = False ! If true, use the fluxes%calving_Hflx field to set the heat carried by runoff, ! instead of using SST*CP*froz_runoff. -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 = "chl_globcolour_monthly_clim.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. -CHL_VARNAME = "CHL" ! default = "CHL_A" - ! Name of CHL_A variable in CHL_FILE. ! === module MOM_energetic_PBL === ML_OMEGA_FRAC = 0.001 ! [nondim] default = 0.0 @@ -687,8 +670,11 @@ LT_MOD_LAC5 = 0.22 ! [nondim] default = 0.95 ! unstable Obukhov depth. ! === module MOM_opacity === -PEN_SW_NBANDS = 3 ! default = 1 - ! The number of bands of penetrating shortwave radiation. +PEN_SW_SCALE = 15.0 ! [m] default = 0.0 + ! The vertical absorption e-folding depth of the penetrating shortwave + ! radiation. +PEN_SW_FRAC = 0.42 ! [nondim] default = 0.0 + ! The fraction of the shortwave radiation that penetrates below the surface. ! === module MOM_tracer_advect === TRACER_ADVECTION_SCHEME = "PPM:H3" ! default = "PLM" @@ -698,8 +684,6 @@ TRACER_ADVECTION_SCHEME = "PPM:H3" ! default = "PLM" ! PPM - Piecewise Parabolic Method (Colella-Woodward) ! === module MOM_tracer_hor_diff === -KHTR = 50.0 ! [m2 s-1] default = 0.0 - ! The background along-isopycnal tracer diffusivity. 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 @@ -709,13 +693,8 @@ MAX_TR_DIFFUSION_CFL = 2.0 ! [nondim] default = -1.0 ! diffusive CFL locally at or below this value. The number of diffusive ! iterations is often this value or the next greater integer. -! === module MOM_neutral_diffusion === -! This module implements neutral diffusion of tracers -USE_NEUTRAL_DIFFUSION = True ! [Boolean] default = False - ! If true, enables the neutral diffusion module. - ! === module MOM_sum_output === -MAXTRUNC = 10000 ! [truncations save_interval-1] default = 0 +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. @@ -733,15 +712,23 @@ EPS_OMESH = 1.0E-13 ! [degrees] default = 1.0E-04 RESTORE_SALINITY = True ! [Boolean] default = False ! If true, the coupled driver will add a globally-balanced fresh-water flux that ! drives sea-surface salinity toward specified values. +ICE_SHELF = True ! [Boolean] default = False + ! If true, enables the ice shelf model. ! === module MOM_surface_forcing_nuopc === LATENT_HEAT_FUSION = 3.337E+05 ! [J/kg] default = 3.34E+05 ! The latent heat of fusion. LATENT_HEAT_VAPORIZATION = 2.501E+06 ! [J/kg] default = 2.5E+06 ! The latent heat of fusion. -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. +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. +ADJUST_NET_SRESTORE_TO_ZERO = False ! [Boolean] default = True + ! If true, adjusts the salinity restoring seen to zero whether restoring is via + ! a salt flux or virtual precip. WIND_STAGGER = "A" ! default = "C" ! A case-insensitive character string to indicate the staggering of the input ! wind stress field. Valid values are 'A', 'B', or 'C'. @@ -753,12 +740,11 @@ FLUXCONST = 0.11 ! [m day-1] default = 0.0 ! anomalies (akin to a piston velocity). Note the non-MKS units. SALT_RESTORE_FILE = "salt_sfc_restore.nc" ! default = "salt_restore.nc" ! A file in which to find the surface salinity to use for restoring. -SALT_RESTORE_VARIABLE = "asalt" ! default = "salt" - ! The name of the surface salinity variable to read from SALT_RESTORE_FILE for - ! restoring salinity. 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. +MASK_SRESTORE = True ! [Boolean] default = False + ! If true, read a file (salt_restore_mask) containing a mask for SSS restoring. USE_RIGID_SEA_ICE = True ! [Boolean] default = False ! If true, sea-ice is rigid enough to exert a nonhydrostatic pressure that ! resist vertical motion. @@ -767,3 +753,41 @@ SEA_ICE_RIGID_MASS = 100.0 ! [kg m-2] default = 1000.0 ! rigidity ! === module MOM_restart === + +! === module MOM_domains min_halo === + +! === module MOM_grid_init === + +! === module MOM_ice_shelf === +SHELF_THERMO = True ! [Boolean] default = False + ! If true, use a thermodynamically interactive ice shelf. +ICE_SHELF_TEMPERATURE = -20.0 ! [degC] default = -15.0 + ! The temperature at the center of the ice shelf. +COL_THICK_MELT_THRESHOLD = 0.001 ! [m] default = 0.0 + ! The minimum ocean column thickness where melting is allowed. +ICE_SHELF_USTAR_FROM_VEL_BUGFIX = True ! [Boolean] default = False + ! Bug fix for ice-area weighting of squared ocean velocities used to calculate + ! friction velocity under ice shelves +ICE_SHELF_TIDEAMP_SCALING_FACTOR = 0.66 ! [nondim] default = 1.0 + ! Scale TIDEAMP_FILE or UTIDE by number in melt parameterisation friction + ! velocity calculation. + +! === module MOM_EOS === +DENSITY_ICE = 917.0 ! [kg m-3] default = 900.0 + ! A typical density of ice. +USTAR_SHELF_BG = 6.0E-04 ! [m s-1] default = 0.0 + ! The minimum value of ustar under ice shelves. +CDRAG_SHELF = 0.007 ! [nondim] default = 0.003 + ! CDRAG is the drag coefficient relating the magnitude of the velocity field to + ! the surface stress. + +! === module MOM_restart === +ICE_PROFILE_CONFIG = "FILE" ! + ! This specifies how the initial ice profile is specified. Valid values are: + ! CHANNEL, FILE, and USER. +ICE_THICKNESS_FILE = "ice_shelf_thick.nc" ! default = "ice_shelf_h.nc" + ! The file from which the bathymetry is read. +ICE_THICKNESS_VARNAME = "thick" ! default = "h_shelf" + ! The name of the thickness variable in ICE_THICKNESS_FILE. +ICE_AREA_VARNAME = "area" ! default = "area_shelf_h" + ! The name of the area variable in ICE_THICKNESS_FILE. diff --git a/docs/available_diags.000000 b/docs/available_diags.000000 index f64735451..cb5f51bc1 100644 --- a/docs/available_diags.000000 +++ b/docs/available_diags.000000 @@ -168,144 +168,6 @@ ! long_name: heat capacity of sea water ! units: J kg-1 K-1 ! variants: {C_p,cpocean} -"MEKE" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: Mesoscale Eddy Kinetic Energy - ! units: m2 s-2 - ! cell_methods: xh:mean yh:mean area:mean -"MEKE_KH" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: MEKE derived diffusivity - ! units: m2 s-1 - ! cell_methods: xh:mean yh:mean area:mean -"MEKE_KU" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: MEKE derived lateral viscosity - ! units: m2 s-1 - ! cell_methods: xh:mean yh:mean area:mean -"MEKE_AU" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: MEKE derived lateral biharmonic viscosity - ! units: m4 s-1 - ! cell_methods: xh:mean yh:mean area:mean -"MEKE_Ue" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: MEKE derived eddy-velocity scale - ! units: m s-1 - ! cell_methods: xh:mean yh:mean area:mean -"MEKE_Ub" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: MEKE derived bottom eddy-velocity scale - ! units: m s-1 - ! cell_methods: xh:mean yh:mean area:mean -"MEKE_Ut" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: MEKE derived barotropic eddy-velocity scale - ! units: m s-1 - ! cell_methods: xh:mean yh:mean area:mean -"MEKE_src" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: MEKE energy source - ! units: m2 s-3 - ! cell_methods: xh:mean yh:mean area:mean -"MEKE_src_adv" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: MEKE energy source from the horizontal advection of MEKE - ! units: m2 s-3 - ! cell_methods: xh:mean yh:mean area:mean -"MEKE_src_btm_drag" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: MEKE energy source from the bottom drag acting on MEKE - ! units: m2 s-3 - ! cell_methods: xh:mean yh:mean area:mean -"MEKE_src_GM" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: MEKE energy source from the thickness mixing (GM scheme) - ! units: m2 s-3 - ! cell_methods: xh:mean yh:mean area:mean -"MEKE_decay" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: MEKE decay rate - ! units: s-1 - ! cell_methods: xh:mean yh:mean area:mean -"MEKE_GM_src" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: MEKE energy available from thickness mixing - ! units: W m-2 - ! cell_methods: xh:mean yh:mean area:mean -"MEKE_mom_src" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: MEKE energy available from momentum - ! units: W m-2 - ! cell_methods: xh:mean yh:mean area:mean -"MEKE_mom_src_bh" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: MEKE energy available from the biharmonic dissipation of momentum - ! units: W m-2 - ! cell_methods: xh:mean yh:mean area:mean -"MEKE_GME_snk" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: MEKE energy lost to GME backscatter - ! units: W m-2 - ! cell_methods: xh:mean yh:mean area:mean -"MEKE_Le" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: Eddy mixing length used in the MEKE derived eddy diffusivity - ! units: m - ! cell_methods: xh:mean yh:mean area:mean -"MEKE_Lrhines" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: Rhines length scale used in the MEKE derived eddy diffusivity - ! units: m - ! cell_methods: xh:mean yh:mean area:mean -"MEKE_Leady" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: Eady length scale used in the MEKE derived eddy diffusivity - ! units: m - ! cell_methods: xh:mean yh:mean area:mean -"MEKE_gamma_b" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: Ratio of bottom-projected eddy velocity to column-mean eddy velocity - ! units: nondim - ! cell_methods: xh:mean yh:mean area:mean -"MEKE_gamma_t" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: Ratio of barotropic eddy velocity to column-mean eddy velocity - ! units: nondim - ! cell_methods: xh:mean yh:mean area:mean -"KHMEKE_u" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xq, yh - ! long_name: Zonal diffusivity of MEKE - ! units: m2 s-1 - ! cell_methods: xq:point yh:mean -"KHMEKE_v" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yq - ! long_name: Meridional diffusivity of MEKE - ! units: m2 s-1 - ! cell_methods: xh:mean yq:point "SN_u" [Unused] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xq, yh @@ -337,45 +199,6 @@ ! units: nondim ! cell_methods: xh:mean yh:mean zl:mean area:mean ! variants: {sqg_struct,sqg_struct_xyave} -"khth_struct" [Unused] - ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} - ! dimensions: xh, yh, zl - ! long_name: Vertical structure of thickness diffusivity - ! units: nondim - ! cell_methods: xh:mean yh:mean zl:mean area:mean - ! variants: {khth_struct,khth_struct_xyave} -"N2_u" [Unused] - ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} - ! dimensions: xq, yh, zi - ! long_name: Square of Brunt-Vaisala frequency, N^2, at u-points, as used in Visbeck et al. - ! units: s-2 - ! cell_methods: xq:point yh:mean zi:point - ! variants: {N2_u,N2_u_xyave} -"N2_v" [Unused] - ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} - ! dimensions: xh, yq, zi - ! long_name: Square of Brunt-Vaisala frequency, N^2, at v-points, as used in Visbeck et al. - ! units: s-2 - ! cell_methods: xh:mean yq:point zi:point - ! variants: {N2_v,N2_v_xyave} -"S2_u" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xq, yh - ! long_name: Depth average square of slope magnitude, S^2, at u-points, as used in Visbeck et al. - ! units: nondim - ! cell_methods: xq:point yh:mean -"S2_v" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yq - ! long_name: Depth average square of slope magnitude, S^2, at v-points, as used in Visbeck et al. - ! units: nondim - ! cell_methods: xh:mean yq:point -"Res_fn" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: Resolution function for scaling diffusivities - ! units: nondim - ! cell_methods: xh:mean yh:mean area:mean "Rd_dx" [Unused] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xh, yh @@ -432,14 +255,14 @@ ! units: m s-1 ! cell_methods: xh:mean yq:point zl:mean ! variants: {Rayleigh_v,Rayleigh_v_xyave} -"uhGM" [Used] +"uhGM" [Unused] ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xq, yh, zl ! long_name: Time Mean Diffusive Zonal Thickness Flux ! units: kg s-1 ! cell_methods: xq:point yh:sum zl:sum ! variants: {uhGM,uhGM_xyave} -"vhGM" [Used] +"vhGM" [Unused] ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xh, yq, zl ! long_name: Time Mean Diffusive Meridional Thickness Flux @@ -513,7 +336,7 @@ ! units: m3 s-1 ! cell_methods: xq:point yh:mean zi:point ! variants: {GM_sfn_x,GM_sfn_x_xyave} -"GM_sfn_y" [Used] +"GM_sfn_y" [Unused] ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xh, yq, zi ! long_name: Parameterized Meridional Overturning Streamfunction @@ -1712,7 +1535,7 @@ ! standard_name: sea_water_mass_per_unit_area ! cell_methods: xh:mean yh:mean zl:sum area:mean ! variants: {masscello,masscello_xyave} -"masso" [Used] +"masso" [Unused] ! modules: ocean_model ! dimensions: scalar ! long_name: Mass of liquid ocean @@ -1733,20 +1556,6 @@ ! units: m ! cell_methods: xh:mean yh:mean zl:sum area:mean ! variants: {h_pre_sync,h_pre_sync_xyave} -"temp" [Used] (CMOR equivalent is "thetao") - ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} - ! dimensions: xh, yh, zl - ! long_name: Potential Temperature - ! units: degC - ! cell_methods: xh:mean yh:mean zl:mean area:mean - ! variants: {temp,temp_xyave,thetao,thetao_xyave} -"salt" [Used] (CMOR equivalent is "so") - ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} - ! dimensions: xh, yh, zl - ! long_name: Salinity - ! units: psu - ! cell_methods: xh:mean yh:mean zl:mean area:mean - ! variants: {salt,salt_xyave,so,so_xyave} "tob" [Used] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xh, yh @@ -1886,7 +1695,7 @@ ! units: m2 s-2 ! cell_methods: xh:mean yh:mean zl:mean area:mean ! variants: {uv,uv_xyave} -"h" [Unused] +"h" [Used] ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xh, yh, zl ! long_name: Layer Thickness @@ -2851,7 +2660,7 @@ ! units: m2 s-1 ! cell_methods: xh:mean yh:mean zi:point area:mean ! variants: {Kd_interface,Kd_interface_xyave} -"Kd_ePBL" [Used] +"Kd_ePBL" [Unused] ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xh, yh, zi ! long_name: ePBL diapycnal diffusivity at interfaces @@ -2997,7 +2806,7 @@ ! units: W m-2 ! cell_methods: xh:mean yh:mean zl:sum area:mean ! variants: {frazil_heat_tendency,frazil_heat_tendency_xyave} -"frazil_heat_tendency_2d" [Unused] +"frazil_heat_tendency_2d" [Used] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xh, yh ! long_name: Depth integrated heat tendency due to frazil formation @@ -3067,7 +2876,7 @@ ! long_name: Buoyancy frequency squared averaged over the water column ! units: s-2 ! cell_methods: xh:mean yh:mean area:mean -"Kd_Itidal_Work" [Used] +"Kd_Itidal_Work" [Unused] ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xh, yh, zl ! long_name: Work done by Internal Tide Diapycnal Mixing @@ -3214,7 +3023,7 @@ ! long_name: The volume flux added to stop the ocean from drying out and becoming negative in depth ! units: m s-1 ! cell_methods: xh:mean yh:mean area:mean -"rsdoabsorb" [Used] +"rsdoabsorb" [Unused] ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xh, yh, zl ! long_name: Convergence of Penetrative Shortwave Flux in Sea Water Layer @@ -3237,19 +3046,13 @@ ! units: W m-2 ! standard_name: nondownwelling_shortwave_flux_in_sea_water ! cell_methods: xh:mean yh:mean area:mean -"Chl_opac" [Used] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: Surface chlorophyll A concentration used to find opacity - ! units: mg m-3 - ! cell_methods: xh:mean yh:mean area:mean -"ePBL_h_ML" [Used] +"ePBL_h_ML" [Unused] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xh, yh ! long_name: Surface boundary layer depth ! units: m ! cell_methods: xh:mean yh:mean area:mean -"h_ML" [Used] +"h_ML" [Unused] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xh, yh ! long_name: Surface mixed layer depth based on active turbulence @@ -3335,7 +3138,7 @@ ! long_name: Increase in mstar due to Langmuir Turbulence. ! units: nondim ! cell_methods: xh:mean yh:mean area:mean -"SW_pen" [Used] +"SW_pen" [Unused] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xh, yh ! long_name: Penetrating shortwave radiation flux into ocean @@ -3354,20 +3157,6 @@ ! units: m-1 ! cell_methods: xh:mean yh:mean zl:mean area:mean ! variants: {opac_1,opac_1_xyave} -"opac_2" [Unused] - ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} - ! dimensions: xh, yh, zl - ! long_name: Opacity for shortwave radiation in band 2, saved as L^-1 tanh(Opacity * L) for L = 10^-10 m - ! units: m-1 - ! cell_methods: xh:mean yh:mean zl:mean area:mean - ! variants: {opac_2,opac_2_xyave} -"opac_3" [Unused] - ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} - ! dimensions: xh, yh, zl - ! long_name: Opacity for shortwave radiation in band 3, saved as L^-1 tanh(Opacity * L) for L = 10^-10 m - ! units: m-1 - ! cell_methods: xh:mean yh:mean zl:mean area:mean - ! variants: {opac_3,opac_3_xyave} "KHTR_u" [Unused] ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xq, yh, zi @@ -3497,18 +3286,6 @@ ! units: psu2 ! cell_methods: xh:mean yh:mean area:mean ! variants: {SSS_sq,sossq} -"conSST" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: Sea Surface Conservative Temperature - ! units: Celsius - ! cell_methods: xh:mean yh:mean area:mean -"absSSS" [Unused] - ! modules: {ocean_model,ocean_model_d2} - ! dimensions: xh, yh - ! long_name: Sea Surface Absolute Salinity - ! units: g kg-1 - ! cell_methods: xh:mean yh:mean area:mean "frazil" [Unused] (CMOR equivalent is "hfsifrazil") ! modules: {ocean_model,ocean_model_d2} ! dimensions: xh, yh @@ -3619,20 +3396,20 @@ ! units: m s-1 ! cell_methods: xh:mean yh:mean zl:sum area:mean ! variants: {dynamics_h_tendency,dynamics_h_tendency_xyave} -"contemp" [Unused] (CMOR equivalent is "bigthetao") +"temp" [Used] (CMOR equivalent is "thetao") ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xh, yh, zl - ! long_name: Conservative Temperature - ! units: Celsius + ! long_name: Potential Temperature + ! units: degC ! cell_methods: xh:mean yh:mean zl:mean area:mean - ! variants: {contemp,contemp_xyave,bigthetao,bigthetao_xyave} -"contemp_post_horzn" [Unused] + ! variants: {temp,temp_xyave,thetao,thetao_xyave} +"temp_post_horzn" [Unused] ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xh, yh, zl - ! long_name: Conservative Temperature after horizontal transport (advection/diffusion) has occurred - ! units: Celsius + ! long_name: Potential Temperature after horizontal transport (advection/diffusion) has occurred + ! units: degC ! cell_methods: xh:mean yh:mean zl:mean area:mean - ! variants: {contemp_post_horzn,contemp_post_horzn_xyave} + ! variants: {temp_post_horzn,temp_post_horzn_xyave} "T_adx" [Unused] ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xq, yh, zl @@ -3678,52 +3455,52 @@ "T_zint" [Unused] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xh, yh - ! long_name: Thickness-weighted integral of Conservative Temperature - ! units: Celsius m + ! long_name: Thickness-weighted integral of Potential Temperature + ! units: degC m ! cell_methods: xh:mean yh:mean area:mean "T_zint_100m" [Unused] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xh, yh - ! long_name: Thickness-weighted integral of Conservative Temperature over top 100m - ! units: Celsius m + ! long_name: Thickness-weighted integral of Potential Temperature over top 100m + ! units: degC m ! cell_methods: xh:mean yh:mean area:mean "T_SURF" [Unused] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xh, yh - ! long_name: Surface values of Conservative Temperature - ! units: Celsius + ! long_name: Surface values of Potential Temperature + ! units: degC ! cell_methods: xh:mean yh:mean area:mean -"T_adx_2d" [Used] +"T_adx_2d" [Unused] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xq, yh ! long_name: Vertically Integrated Advective Zonal Flux of Heat ! units: W ! cell_methods: xq:point yh:sum -"T_ady_2d" [Used] +"T_ady_2d" [Unused] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xh, yq ! long_name: Vertically Integrated Advective Meridional Flux of Heat ! units: W ! cell_methods: xh:sum yq:point -"T_diffx_2d" [Used] +"T_diffx_2d" [Unused] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xq, yh ! long_name: Vertically Integrated Diffusive Zonal Flux of Heat ! units: W ! cell_methods: xq:point yh:sum -"T_diffy_2d" [Used] +"T_diffy_2d" [Unused] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xh, yq ! long_name: Vertically Integrated Diffusive Meridional Flux of Heat ! units: W ! cell_methods: xh:sum yq:point -"T_hbd_diffx_2d" [Used] +"T_hbd_diffx_2d" [Unused] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xq, yh ! long_name: Vertically-integrated zonal diffusive flux from the horizontal boundary diffusion scheme for Heat ! units: W ! cell_methods: xq:point yh:sum -"T_hbd_diffy_2d" [Used] +"T_hbd_diffy_2d" [Unused] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xh, yq ! long_name: Vertically-integrated meridional diffusive flux from the horizontal boundary diffusion scheme for Heat @@ -3745,8 +3522,8 @@ "T_tendency" [Unused] ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xh, yh, zl - ! long_name: Net time tendency for conservative temperature - ! units: Celsius s-1 + ! long_name: Net time tendency for potential temperature + ! units: degC s-1 ! cell_methods: xh:mean yh:mean zl:mean area:mean ! variants: {T_tendency,T_tendency_xyave} "T_dfxy_cont_tendency" [Unused] (CMOR equivalent is "opottemppmdiff") @@ -3780,14 +3557,14 @@ ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xh, yh, zl ! long_name: Neutral diffusion tracer concentration tendency for T - ! units: Celsius s-1 + ! units: degC s-1 ! cell_methods: xh:mean yh:mean zl:mean area:mean ! variants: {T_dfxy_conc_tendency,T_dfxy_conc_tendency_xyave} "T_hbdxy_conc_tendency" [Unused] ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xh, yh, zl ! long_name: Horizontal diffusion tracer concentration tendency for T - ! units: Celsius s-1 + ! units: degC s-1 ! cell_methods: xh:mean yh:mean zl:mean area:mean ! variants: {T_hbdxy_conc_tendency,T_hbdxy_conc_tendency_xyave} "Th_tendency" [Unused] (CMOR equivalent is "opottemptend") @@ -3807,8 +3584,8 @@ "T_tendency_vert_remap" [Unused] ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xh, yh, zl - ! long_name: Vertical remapping tracer concentration tendency for contemp - ! units: Celsius s-1 + ! long_name: Vertical remapping tracer concentration tendency for temp + ! units: degC s-1 ! cell_methods: xh:mean yh:mean zl:mean area:mean ! variants: {T_tendency_vert_remap,T_tendency_vert_remap_xyave} "Th_tendency_vert_remap" [Unused] @@ -3827,24 +3604,24 @@ "T_vardec" [Unused] ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xh, yh, zl - ! long_name: ALE variance decay for conservative temperature - ! units: Celsius2 s-1 + ! long_name: ALE variance decay for potential temperature + ! units: degC2 s-1 ! cell_methods: xh:mean yh:mean zl:mean area:mean ! variants: {T_vardec,T_vardec_xyave} -"abssalt" [Unused] (CMOR equivalent is "absso") +"salt" [Used] (CMOR equivalent is "so") ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xh, yh, zl - ! long_name: Absolute Salinity - ! units: g kg-1 + ! long_name: Salinity + ! units: psu ! cell_methods: xh:mean yh:mean zl:mean area:mean - ! variants: {abssalt,abssalt_xyave,absso,absso_xyave} -"abssalt_post_horzn" [Unused] + ! variants: {salt,salt_xyave,so,so_xyave} +"salt_post_horzn" [Unused] ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xh, yh, zl - ! long_name: Absolute Salinity after horizontal transport (advection/diffusion) has occurred - ! units: g kg-1 + ! long_name: Salinity after horizontal transport (advection/diffusion) has occurred + ! units: psu ! cell_methods: xh:mean yh:mean zl:mean area:mean - ! variants: {abssalt_post_horzn,abssalt_post_horzn_xyave} + ! variants: {salt_post_horzn,salt_post_horzn_xyave} "S_adx" [Unused] ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xq, yh, zl @@ -3890,20 +3667,20 @@ "S_zint" [Unused] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xh, yh - ! long_name: Thickness-weighted integral of Absolute Salinity - ! units: g kg-1 m + ! long_name: Thickness-weighted integral of Salinity + ! units: psu m ! cell_methods: xh:mean yh:mean area:mean "S_zint_100m" [Unused] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xh, yh - ! long_name: Thickness-weighted integral of Absolute Salinity over top 100m - ! units: g kg-1 m + ! long_name: Thickness-weighted integral of Salinity over top 100m + ! units: psu m ! cell_methods: xh:mean yh:mean area:mean "S_SURF" [Unused] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xh, yh - ! long_name: Surface values of Absolute Salinity - ! units: g kg-1 + ! long_name: Surface values of Salinity + ! units: psu ! cell_methods: xh:mean yh:mean area:mean "S_adx_2d" [Unused] ! modules: {ocean_model,ocean_model_d2} @@ -3957,8 +3734,8 @@ "S_tendency" [Unused] ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xh, yh, zl - ! long_name: Net time tendency for absolute salinity - ! units: g kg-1 s-1 + ! long_name: Net time tendency for salinity + ! units: psu s-1 ! cell_methods: xh:mean yh:mean zl:mean area:mean ! variants: {S_tendency,S_tendency_xyave} "S_dfxy_cont_tendency" [Unused] (CMOR equivalent is "osaltpmdiff") @@ -3992,14 +3769,14 @@ ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xh, yh, zl ! long_name: Neutral diffusion tracer concentration tendency for S - ! units: g kg-1 s-1 + ! units: psu s-1 ! cell_methods: xh:mean yh:mean zl:mean area:mean ! variants: {S_dfxy_conc_tendency,S_dfxy_conc_tendency_xyave} "S_hbdxy_conc_tendency" [Unused] ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xh, yh, zl ! long_name: Horizontal diffusion tracer concentration tendency for S - ! units: g kg-1 s-1 + ! units: psu s-1 ! cell_methods: xh:mean yh:mean zl:mean area:mean ! variants: {S_hbdxy_conc_tendency,S_hbdxy_conc_tendency_xyave} "Sh_tendency" [Unused] (CMOR equivalent is "osalttend") @@ -4019,8 +3796,8 @@ "S_tendency_vert_remap" [Unused] ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xh, yh, zl - ! long_name: Vertical remapping tracer concentration tendency for abssalt - ! units: g kg-1 s-1 + ! long_name: Vertical remapping tracer concentration tendency for salt + ! units: psu s-1 ! cell_methods: xh:mean yh:mean zl:mean area:mean ! variants: {S_tendency_vert_remap,S_tendency_vert_remap_xyave} "Sh_tendency_vert_remap" [Unused] @@ -4039,8 +3816,8 @@ "S_vardec" [Unused] ! modules: {ocean_model,ocean_model_z,ocean_model_rho2,ocean_model_d2,ocean_model_z_d2,ocean_model_rho2_d2} ! dimensions: xh, yh, zl - ! long_name: ALE variance decay for absolute salinity - ! units: (g kg-1)2 s-1 + ! long_name: ALE variance decay for salinity + ! units: psu2 s-1 ! cell_methods: xh:mean yh:mean zl:mean area:mean ! variants: {S_vardec,S_vardec_xyave} "age" [Used] (CMOR equivalent is "agessc") @@ -4444,7 +4221,7 @@ ! standard_name: water_flux_into_sea_water ! cell_methods: xh:mean yh:mean area:mean ! variants: {PRCmE,wfo} -"evap" [Used] (CMOR equivalent is "evs") +"evap" [Unused] (CMOR equivalent is "evs") ! modules: {ocean_model,ocean_model_d2} ! dimensions: xh, yh ! long_name: Evaporation/condensation at ocean surface (evaporation is negative) @@ -4518,45 +4295,45 @@ ! units: kg m-2 s-1 ! standard_name: water_flux_into_sea_water_from_glaciers ! cell_methods: xh:mean yh:mean area:mean -"net_massout" [Used] +"net_massout" [Unused] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xh, yh ! long_name: Net mass leaving the ocean due to evaporation, seaice formation ! units: kg m-2 s-1 ! cell_methods: xh:mean yh:mean area:mean -"net_massin" [Used] +"net_massin" [Unused] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xh, yh ! long_name: Net mass entering ocean due to precip, runoff, ice melt ! units: kg m-2 s-1 ! cell_methods: xh:mean yh:mean area:mean -"massout_flux" [Used] +"massout_flux" [Unused] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xh, yh ! long_name: Net mass flux of freshwater out of the ocean (used in the boundary flux calculation) ! units: kg m-2 ! cell_methods: xh:mean yh:mean area:mean -"massin_flux" [Used] +"massin_flux" [Unused] ! modules: {ocean_model,ocean_model_d2} ! dimensions: xh, yh ! long_name: Net mass flux of freshwater into the ocean (used in boundary flux calculation) ! units: kg m-2 ! cell_methods: xh:mean yh:mean area:mean -"total_PRCmE" [Used] +"total_PRCmE" [Unused] ! modules: ocean_model ! dimensions: scalar ! long_name: Area integrated net surface water flux (precip+melt+liq runoff+ice calving-evap) ! units: kg s-1 ! standard_name: water_flux_into_sea_water_area_integrated ! variants: {total_PRCmE,total_wfo} -"total_evap" [Used] +"total_evap" [Unused] ! modules: ocean_model ! dimensions: scalar ! long_name: Area integrated evap/condense at ocean surface ! units: kg s-1 ! standard_name: water_evaporation_flux_area_integrated ! variants: {total_evap,total_evs} -"total_icemelt" [Used] +"total_icemelt" [Unused] ! modules: ocean_model ! dimensions: scalar ! long_name: Area integrated sea ice melt (>0) or form (<0) @@ -4587,13 +4364,13 @@ ! dimensions: scalar ! long_name: Area integrated virtual liquid precip due to SSS restoring ! units: kg s-1 -"total_frunoff" [Used] +"total_frunoff" [Unused] ! modules: ocean_model ! dimensions: scalar ! long_name: Area integrated frozen runoff (calving) & iceberg melt into ocean ! units: kg s-1 ! variants: {total_frunoff,total_ficeberg} -"total_lrunoff" [Used] +"total_lrunoff" [Unused] ! modules: ocean_model ! dimensions: scalar ! long_name: Area integrated liquid runoff into ocean @@ -4617,7 +4394,7 @@ "total_net_massin" [Used] ! modules: ocean_model ! dimensions: scalar - ! long_name: Area integrated mass entering ocean due to precip, runoff, ice melt + ! long_name: Area integrated mass entering ocean due to predip, runoff, ice melt ! units: kg s-1 "PRCmE_ga" [Unused] ! modules: ocean_model diff --git a/drof.streams.xml b/drof.streams.xml index 895132b06..0e5660cc5 100644 --- a/drof.streams.xml +++ b/drof.streams.xml @@ -5,26 +5,6 @@ 2025-07-04 10:28:08 Created by as2285 (anton-seaice) on 2025-07-04, using https://github.com/ACCESS-NRI/om3-scripts/blob/ffad24e7448649c5797bffe6a4e6184741dd87a5/data_stream_xml_generation/generate_xml_drof.py: /g/data/vk83/apps/om3-scripts/data_stream_xml_generation/generate_xml_drof.py 1900 1900 - - cycle - 1.0 - linear - 0 - single - bilinear - 1900 - 1900 - 1900 - null - ./INPUT/JRA55do-drof-ESMFmesh.nc - null - - ./INPUT/RYF.friver.1990_1991.nc - - - friver Forr_rofl - - cycle 1.0 diff --git a/drof_in b/drof_in index 9fe11d57c..e08e43781 100644 --- a/drof_in +++ b/drof_in @@ -1,9 +1,9 @@ &drof_nml datamode = "copyall" - model_maskfile = "./INPUT/access-om3-25km-nomask-ESMFmesh.nc" - model_meshfile = "./INPUT/access-om3-25km-nomask-ESMFmesh.nc" - nx_global = 1440 - ny_global = 1152 + model_maskfile = "./INPUT/access-om3-4km-nomask-ESMFmesh.nc" + model_meshfile = "./INPUT/access-om3-4km-nomask-ESMFmesh.nc" + nx_global = 4320 + ny_global = 1442 restfilm = "null" skip_restart_read = .false. / diff --git a/ice_in b/ice_in index cb8d79a51..432bd20a5 100644 --- a/ice_in +++ b/ice_in @@ -7,7 +7,7 @@ hist_suffix = ".1day.mean", ".1mon.mean", "x", "x", "x" hist_time_axis = "middle" history_deflate = 1 - history_chunksize = 720, 540 + history_chunksize = 720, 186 history_precision = 8 restart_deflate = 1 ice_ic = 'none' @@ -15,20 +15,21 @@ npt = 35040 pointer_file = './rpointer.ice' print_global = .false. + ndtd = 4 / &grid_nml - bathymetry_file = "./INPUT/topog.nc" + bathymetry_file = "./INPUT/topog_seaice.nc" grid_atm = "A" grid_file = "./INPUT/ocean_hgrid.nc" grid_format = "mom_nc" grid_ice = "B" grid_ocn = "A" - grid_type = "tripole" + grid_type = "regional" grid_outfile = .true. - kcatbound = 0 + kcatbound = 2 kmt_file = "./INPUT/kmt.nc" nblyr = 1 - ncat = 5 + ncat = 7 nfsd = 1 nilyr = 4 nslyr = 1 @@ -68,17 +69,17 @@ ustar_min = 0.0005 / &domain_nml - block_size_x = 30 - block_size_y = 27 - distribution_type = "roundrobin" + block_size_x = 48 + block_size_y = 36 + distribution_type = "sectrobin" distribution_wght = "latitude" maskhalo_bound = .true. maskhalo_dyn = .true. maskhalo_remap = .true. max_blocks = -1 - ns_boundary_type = "tripole" - nx_global = 1440 - ny_global = 1152 + ns_boundary_type = "open" + nx_global = 4320 + ny_global = 1442 processor_shape = "square-ice" / &ice_prescribed_nml diff --git a/input.nml b/input.nml index 15c817e8c..68b070de1 100644 --- a/input.nml +++ b/input.nml @@ -2,7 +2,8 @@ output_directory = './' restart_input_dir = './' restart_output_dir = './' - parameter_filename = 'MOM_input', 'MOM_override' + input_filename = 'n' + parameter_filename = 'MOM_input' / &fms_nml diff --git a/manifests/exe.yaml b/manifests/exe.yaml index 5101ae285..0feb564ae 100644 --- a/manifests/exe.yaml +++ b/manifests/exe.yaml @@ -2,7 +2,7 @@ format: yamanifest version: 1.0 --- work/access-om3-MOM6-CICE6: - fullpath: /g/data/vk83/apps/spack/0.22/release/linux-rocky8-x86_64_v4/oneapi-2025.2.0/access3-2025.08.000-n2qsfzdl7syl4k6m4dx3ydqxarpkmw6l/bin/access-om3-MOM6-CICE6 + fullpath: /g/data/vk83/prerelease/apps/spack/0.22/release/linux-rocky8-x86_64_v4/oneapi-2025.2.0/access3-2025.08.000-vvi3x2ud3qpmwkhecu4okiax5goguqho/bin/access-om3-MOM6-CICE6 hashes: - binhash: 2c767b4b07ae2e86e21b4a50582de8df - md5: cda3625020745cae76fcfc1f800ef4b3 + binhash: ffd7d05f91f51efa0506e471bfb28c06 + md5: a0dba1d7df3aaeced3d2619333c8b2b6 diff --git a/manifests/input.yaml b/manifests/input.yaml index b6bdde012..6bf2af6eb 100644 --- a/manifests/input.yaml +++ b/manifests/input.yaml @@ -11,6 +11,11 @@ work/INPUT/JRA55do-drof-ESMFmesh.nc: hashes: binhash: 9f3ee3ce42e776fcbee8026e0c3880f5 md5: bf62c4e7786ddfbe56d61b3074d5d16d +work/INPUT/OBC37S_forcing_access_yr2_4km.nc: + fullpath: /g/data/vk83/prerelease/configurations/inputs/access-om3/panan.4km/2026.01.08/OBC37S_forcing_access_yr2_4km.nc + hashes: + binhash: 39b50148464d56adc7845d79189c0208 + md5: 487153ac346d933eb3a0f1fefd991a8c work/INPUT/RYF.friver.1990_1991.nc: fullpath: /g/data/vk83/configurations/inputs/JRA-55/RYF/v1-6/data/RYF.friver.1990_1991.nc hashes: @@ -66,63 +71,83 @@ work/INPUT/RYF.vas.1990_1991.nc: hashes: binhash: f353a71a014062723ed2b56353932caf md5: 30b92b1e409cff9a9900692aa961aba4 -work/INPUT/access-om3-25km-ESMFmesh.nc: - fullpath: /g/data/vk83/configurations/inputs/access-om3/share/meshes/global.25km/2025.11.27/access-om3-25km-ESMFmesh.nc +work/INPUT/access-om3-4km-ESMFmesh.nc: + fullpath: /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/access-om3-4km-ESMFmesh.nc hashes: - binhash: 772cf9c667ca5e67a3d2870890e2ea55 - md5: 0ae5ab11f801a02f2537fb7e1ec02674 -work/INPUT/access-om3-25km-nomask-ESMFmesh.nc: - fullpath: /g/data/vk83/configurations/inputs/access-om3/share/meshes/global.25km/2025.11.27/access-om3-25km-nomask-ESMFmesh.nc + binhash: 83ca4dae88dcaa5f405784c536fc4757 + md5: 7dc01fc3a9c875f93690d6c7528b1b70 +work/INPUT/access-om3-4km-nomask-ESMFmesh.nc: + fullpath: /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/access-om3-4km-nomask-ESMFmesh.nc hashes: - binhash: bc20c96c3827517fb82b157b1062a886 - md5: c3cdd8ea948c41e78e55709ebaf4ca91 -work/INPUT/access-om3-25km-rof-remap-weights.nc: - fullpath: /g/data/vk83/configurations/inputs/access-om3/cmeps/remap_weights/global.25km/2025.11.27/access-om3-25km-rof-remap-weights.nc + binhash: 95bde3efc00d398d00c3bb6790f8023b + md5: 28d78cf58801f559effb201dc2bf3ec4 +work/INPUT/access-om3-4km-rof-remap-weights.nc: + fullpath: /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/access-om3-4km-rof-remap-weights.nc hashes: - binhash: 39c2a25bec7d738a530ebd6a5fe603ef - md5: 80067fb4b5b6a46e483e6d346636b3db + binhash: 0b87e16c64b846abcda73a04841c0b80 + md5: 38c29881d380d2911fb31bbc1f987391 work/INPUT/bottom_roughness.nc: - fullpath: /g/data/vk83/configurations/inputs/access-om3/mom/tidal_external_files/global.25km/2025.09.03/bottom_roughness.nc + fullpath: /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/bottom_roughness.nc + hashes: + binhash: 8a947441b2ee9162af7d4a3030bb4153 + md5: 7f82353f4add21c0694a793d4480e50f +work/INPUT/diag_rho2.nc: + fullpath: /g/data/vk83/prerelease/configurations/inputs/access-om3/panan.4km/2026.01.08/diag_rho2.nc hashes: - binhash: 5aee86b7b6d0fcdb232dd9412c37d62f - md5: 99698827fe5c5129dce971d006e4fdfc -work/INPUT/chl_globcolour_monthly_clim.nc: - fullpath: /g/data/vk83/configurations/inputs/access-om3/mom/chlorophyll/global.25km/2025.11.21/chl_globcolour_monthly_clim.nc + binhash: fab458dc2a1c9fbb9ef0b6420aebaa78 + md5: b8abcb148bb3ef71282d6be4531e491e +work/INPUT/ice_shelf_thick.nc: + fullpath: /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/ice_shelf_thick.nc hashes: - binhash: f30a3f69bee2f2174596d238e1de822e - md5: 269f4492ae65f6ef9dd6e2a59f2bff4b + binhash: a92ad5eefa1e3b9adf0992d7cebe7c93 + md5: e3b33e2834ffd0e638f890173e000ce0 work/INPUT/kmt.nc: - fullpath: /g/data/vk83/configurations/inputs/access-om3/cice/grids/global.25km/2025.11.27/kmt.nc + fullpath: /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/kmt.nc + hashes: + binhash: 8d34a488332e54be5d932a997c44d9ab + md5: 551a62a9e522b16126bae843e257313f +work/INPUT/mask_table.1490.84x49: + fullpath: /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/mask_table.1490.84x49 hashes: - binhash: eb989f5ba7c3c3bae2551daefe151698 - md5: 64bc07e338d1a5d4bc2364f49e77fb96 + binhash: ad56cbcee4ee95de2389a555f0fe14d5 + md5: 2d21eca0c746628fc47f1c32ea500e11 work/INPUT/ocean_hgrid.nc: - fullpath: /g/data/vk83/configurations/inputs/access-om3/mom/grids/mosaic/global.25km/2025.09.02/ocean_hgrid.nc + fullpath: /g/data/vk83/prerelease/configurations/inputs/access-om3/panan.4km/2026.01.08/ocean_hgrid.nc hashes: - binhash: a7882fb994abd0e402809e39545dc2c5 - md5: dd9cb78712fedbdd3e5930f7066bcac8 + binhash: 931f526ede0691b223ca4bb42879eeb3 + md5: a1876832e039ca9bf761b8c16a08449e +work/INPUT/ocean_temp_salt.res.nc: + fullpath: /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/ocean_temp_salt.res.nc + hashes: + binhash: 8345a635ab30295f2cd6b3640f591223 + md5: 4bb7af6ab6b226e08760931cb1748fce work/INPUT/ocean_vgrid.nc: - fullpath: /g/data/vk83/configurations/inputs/access-om3/mom/grids/vertical/global.25km/2025.03.12/ocean_vgrid.nc + fullpath: /g/data/vk83/prerelease/configurations/inputs/access-om3/panan.4km/2026.01.08/ocean_vgrid.nc hashes: - binhash: cbc588e29d819ce8898d77760a65e20d + binhash: 55ccd1cf1887ece8c2b5788600f49100 md5: 66f71f2b24bd00b1c220787505e0106e +work/INPUT/salt_restore_mask.nc: + fullpath: /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/salt_restore_mask.nc + hashes: + binhash: 5cb019b840db395e08a95c32001f31d9 + md5: 89e3a69b96d12492c0baab68d3bc9b79 work/INPUT/salt_sfc_restore.nc: - fullpath: /g/data/vk83/configurations/inputs/access-om3/mom/surface_salt_restoring/global.25km/2025.10.24/salt_sfc_restore.nc + fullpath: /g/data/vk83/prerelease/configurations/inputs/access-om3/panan.4km/2026.01.08/salt_sfc_restore.nc hashes: - binhash: 54166d28383dfc561138638479cec1db - md5: 0e59f79a3da9a989eb9ce854277599f8 + binhash: 8b43aba8e111b98d2294c4733ea573ab + md5: 62bcbd4638b849e1fada4f4e08dc228c work/INPUT/tideamp.nc: - fullpath: /g/data/vk83/configurations/inputs/access-om3/mom/tidal_external_files/global.25km/2025.09.03/tideamp.nc + fullpath: /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/tideamp.nc hashes: - binhash: 0b69afaa564866597b1232be5ec00dad - md5: 1aa1db4e7b14f0486d3546b9bc2a035a + binhash: 115f80c9b9f0600c76a4c348a89e7d99 + md5: 12af3ec0ed46cdf425dc72f869567dde work/INPUT/topog.nc: - fullpath: /g/data/vk83/configurations/inputs/access-om3/share/grids/global.25km/2025.11.27/topog.nc + fullpath: /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/topog.nc hashes: - binhash: 6e15c0bea662cdb8f7df8d75af70ad3a - md5: 3a1d6b02a276985b038ddaafac909083 -work/INPUT/woa23_ts_01_mom.nc: - fullpath: /g/data/vk83/configurations/inputs/access-om3/mom/initial_conditions/global.25km/2025.10.24/woa23_ts_01_mom.nc + binhash: 4c89085e6ca52f1c41f7c1a2f056a9e6 + md5: 18ed8e8bd321b8f4375ece4e9afe2f5f +work/INPUT/topog_seaice.nc: + fullpath: /g/data/vk83/prerelease/configurations/inputs/access-om3/panan-isf.4km/2026.01.08/topog_seaice.nc hashes: - binhash: a56e68e159c8007dd75cfa3d2d22cb6c - md5: bd9559e8a3dbefa16080d0340b26be12 + binhash: 22d9f7a24ee41e4f0fc024b5d4511c2a + md5: b6d663421561dd6478744c63216e5b8d diff --git a/nuopc.runconfig b/nuopc.runconfig index 7ec7e3a65..b0294a46f 100644 --- a/nuopc.runconfig +++ b/nuopc.runconfig @@ -27,28 +27,28 @@ DRIVER_attributes:: :: PELAYOUT_attributes:: - atm_ntasks = 291 + atm_ntasks = 364 atm_nthreads = 1 atm_pestride = 1 atm_rootpe = 0 - cpl_ntasks = 291 + cpl_ntasks = 364 cpl_nthreads = 1 cpl_pestride = 1 cpl_rootpe = 0 esmf_logging = ESMF_LOGKIND_NONE - ice_ntasks = 291 + ice_ntasks = 364 ice_nthreads = 1 ice_pestride = 1 ice_rootpe = 0 ninst = 1 - ocn_ntasks = 1685 + ocn_ntasks = 2626 ocn_nthreads = 1 ocn_pestride = 1 - ocn_rootpe = 291 + ocn_rootpe = 364 pio_asyncio_ntasks = 0 pio_asyncio_rootpe = 1 pio_asyncio_stride = 0 - rof_ntasks = 291 + rof_ntasks = 364 rof_nthreads = 1 rof_pestride = 1 rof_rootpe = 0 @@ -98,11 +98,11 @@ ALLCOMP_attributes:: hostname = gadi ice_ncat = 5 mediator_present = true - mesh_atm = ./INPUT/access-om3-25km-ESMFmesh.nc - mesh_ice = ./INPUT/access-om3-25km-ESMFmesh.nc + mesh_atm = ./INPUT/access-om3-4km-nomask-ESMFmesh.nc + mesh_ice = ./INPUT/access-om3-4km-ESMFmesh.nc mesh_lnd = UNSET - mesh_mask = ./INPUT/access-om3-25km-ESMFmesh.nc - mesh_ocn = ./INPUT/access-om3-25km-ESMFmesh.nc + mesh_mask = ./INPUT/access-om3-4km-ESMFmesh.nc + mesh_ocn = ./INPUT/access-om3-4km-ESMFmesh.nc model_version = unknown ocn2glc_coupling = .false. ocn2glc_levels = 1:10:19:26:30:33:35 @@ -119,6 +119,7 @@ ALLCOMP_attributes:: tfreeze_option = linear_salt wav_coupling_to_cice = .false. restart_pointer_append_date = .false. + additional_restart_dir = RESTART :: MED_attributes:: @@ -128,8 +129,8 @@ MED_attributes:: atm2lnd_map = unset atm2ocn_map = unset atm2wav_map = unset - atm_nx = 1440 - atm_ny = 1152 + atm_nx = 4320 + atm_ny = 1442 budget_ann = 1 budget_daily = 0 budget_inst = 0 @@ -234,8 +235,8 @@ MED_attributes:: history_option_wav_inst = never ice2atm_map = unset ice2wav_smapname = unset - ice_nx = 1440 - ice_ny = 1152 + ice_nx = 4320 + ice_ny = 1442 info_debug = 1 lnd2atm_map = unset lnd2rof_map = unset @@ -244,15 +245,15 @@ MED_attributes:: mapuv_with_cart3d = .true. ocn2atm_map = unset ocn2wav_smapname = unset - ocn_nx = 1440 - ocn_ny = 1152 + ocn_nx = 4320 + ocn_ny = 1442 ocn_surface_flux_scheme = 0 rof2lnd_map = unset rof2ocn_fmapname = unset - rof2ocn_ice_rmapname = ./INPUT/access-om3-25km-rof-remap-weights.nc - rof2ocn_liq_rmapname = ./INPUT/access-om3-25km-rof-remap-weights.nc - rof_nx = 1440 - rof_ny = 1152 + rof2ocn_ice_rmapname = ./INPUT/access-om3-4km-rof-remap-weights.nc + rof2ocn_liq_rmapname = ./INPUT/access-om3-4km-rof-remap-weights.nc + rof_nx = 4320 + rof_ny = 1442 wav2ocn_smapname = unset wav_nx = 0 wav_ny = 0 @@ -267,15 +268,15 @@ CLOCK_attributes:: history_ymd = -999 ice_cpl_dt = 99999 #not used lnd_cpl_dt = 99999 #not used - ocn_cpl_dt = 900 #ignored (coupling timestep set by nuopc.runseq) unless stop_option is nsteps + ocn_cpl_dt = 150 #ignored (coupling timestep set by nuopc.runseq) unless stop_option is nsteps restart_n = 1 - restart_option = nyears + restart_option = nmonths restart_ymd = -999 rof_cpl_dt = 99999 #not used start_tod = 0 start_ymd = 19000101 stop_n = 1 - stop_option = nyears + stop_option = nmonths stop_tod = 0 stop_ymd = -999 tprof_n = -999 @@ -302,7 +303,7 @@ OCN_attributes:: ROF_attributes:: Verbosity = off - mesh_rof = ./INPUT/access-om3-25km-nomask-ESMFmesh.nc + mesh_rof = ./INPUT/access-om3-4km-nomask-ESMFmesh.nc :: MED_modelio:: diff --git a/nuopc.runseq b/nuopc.runseq index 4f1c3a71c..c026759b2 100644 --- a/nuopc.runseq +++ b/nuopc.runseq @@ -1,5 +1,5 @@ runSeq:: -@900 +@150 MED med_phases_aofluxes_run MED med_phases_prep_ocn_accum MED med_phases_ocnalb_run diff --git a/testing/checksum/historical-6hr-checksum.json b/testing/checksum/historical-6hr-checksum.json index ff54744e5..ff6bb69e6 100644 --- a/testing/checksum/historical-6hr-checksum.json +++ b/testing/checksum/historical-6hr-checksum.json @@ -1,101 +1,104 @@ { "schema_version": "1-0-0", "output": { - "Salt": [ - "7FB53E6271F127AA" - ], "Temp": [ - "EA900B9C74D1A2EC" + "2EAE1AE4BE468357" + ], + "Salt": [ + "486A28B7627A65DA" ], "h": [ - "599923CF98C8FC84" + "D30C4E9AFA5A63E1" ], "u": [ - "18AAC8E0B9D8B3AB" - ], - "CAu": [ - "E23EBEBB53824FA6" + "E3020B4A47756A4" ], "First_direction": [ "0" ], "ave_ssh": [ - "5B67935AFF6B4C0B" + "4B998CD7D98C596D" ], "frazil": [ - "224A5876C7B1F2D0" + "528FBEB185CC5F68" ], - "p_surf_EOS": [ - "ACB47429463E107" + "hML": [ + "72118BCAED827BEF" ], "sfc": [ - "8690E36A3D932FFC" - ], - "u2": [ - "C6C1EF74A04490B3" + "35E801C1470AA084" ], "v": [ - "A15C8F831E3FEF47" + "5F0E1BA3A7A6ADB" ], - "v2": [ - "742F77B1DAF3BC44" - ], - "CAv": [ - "96BA978BAC3CE46D" - ], - "DTBT": [ - "4032D871A0314B1B" - ], - "MEKE": [ - "A8680CFDEE272F34" + "u2": [ + "D5A6812749A9EF95" ], - "MEKE_Kh": [ - "62F028F6A7DA8441" + "v2": [ + "DAC71C722F259F37" ], - "MEKE_Ku": [ - "197667D8C0685859" + "CAu": [ + "D7D0BA68288FC556" ], - "age": [ - "F170693B459BC917" + "CAv": [ + "FD785F813F7A58D6" ], "diffu": [ - "B488A16C6796CC5B" + "20135250F67AE46B" + ], + "DTBT": [ + "4013AB118F4D8ED4" ], "diffv": [ - "1895EDBECC734BA4" + "4F803E5A0DD775B" ], "ubtav": [ - "36E7F2C1637344E7" + "C047148D25F32C8" ], "vbtav": [ - "BD7BD2391C34BFE0" + "948445BAD64C538D" + ], + "age": [ + "8864E57E2340FBF4" ], "Kd_shear": [ - "E3F43FF114372798" + "AB29A6F5327B0EE" ], "Kv_shear": [ - "C921537BCC7494D7" + "F9CCEF533CBF577D" ], "Kv_shear_Bu": [ - "B1EB0C2FAF0EE796" + "3D4B01B5BCBF5E98" ], "MLD": [ - "4C7107FDD8FC43EC" + "72118BCAED827BEF" ], "MLD_MLE_filtered": [ - "41DF9C26239D3E2E" + "CB085731B2F5074A" ], "MLD_MLE_filtered_slow": [ - "248AD2FE41F12992" + "59FDD0B5EEC74EBE" ], "MLE_Bflux": [ - "3FD53169F90442F3" + "A01F92B5348DF534" ], "SFC_BFLX": [ - "239CBF0347136EE6" + "71A61F881BF4EE70" ], "h_ML": [ - "4C7107FDD8FC43EC" + "72118BCAED827BEF" + ], + "rx_normal": [ + "0" + ], + "ry_normal": [ + "0" + ], + "tres_y_001": [ + "0" + ], + "tres_y_002": [ + "0" ] } } \ No newline at end of file