Skip to content

MOM6 PR test: gfdl-to-main-2025-09-25 on MC 100km RYF+WOMBATlite #932

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Closed
dougiesquire wants to merge 2 commits into from
Closed

MOM6 PR test: gfdl-to-main-2025-09-25 on MC 100km RYF+WOMBATlite #932

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
fd20032
Use access-om3/pr164-1
dougiesquire Nov 28, 2025
b110b64
Set parameters whose defaults changed
dougiesquire Nov 28, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
113 changes: 76 additions & 37 deletions MOM_input
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -143,7 +143,7 @@ NK = 75 ! [nondim]
! The number of model layers.

! === module MOM_EOS ===
EQN_OF_STATE = "ROQUET_RHO" ! default = "WRIGHT"
EQN_OF_STATE = "ROQUET_RHO" ! default = "WRIGHT_FULL"
! 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
Expand Down Expand Up @@ -193,8 +193,23 @@ ALE_COORDINATE_CONFIG = "FILE:ocean_vgrid.nc,interfaces=zeta" ! default = "UNIFO
! FNC1:string - FNC1:dz_min,H_total,power,precision
! HYBRID:string - read from a file. The string specifies
! the filename and two variable names, separated
! by a comma or space, for sigma-2 and dz. e.g.
! HYBRID:vgrid.nc,sigma2,dz
! by a comma or space, for sigma-2 and dz.
! e.g. HYBRID:vgrid.nc,sigma2,dz
! HYBRID_3D:string - read from a file. The string specifies
! the filename and two 3D variable names, separated
! by a comma or space, for sigma-2 and dz. The
! latter can be FNC1:string which is used everywhere.
! e.g. HYBRID_3D:vgrid.nc,sigma2,dz
! HYBRID_MAP:string - read from a file. The string specifies
! the filename and three variable names, separated
! by a comma or space, for map, sigma-2 and dz.
! Map is a spatial index array with, maxval(map)=N,
! and the others are 2D arrays containing N profiles.
! Map typically contains integer values, but it can
! contain real values, I+w, which imply using
! the weighted sum of profiles I and I+1.
! Dz can be FNC1:string which is used everywhere.
! e.g. HYBRID_MAP:vgrid.nc,map,sigma2,dz
!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
REMAPPING_SCHEME = "PPM_H4" ! default = "PLM"
! This sets the reconstruction scheme used for vertical remapping for all
Expand Down Expand Up @@ -256,13 +271,23 @@ DIAG_COORD_DEF_Z = "FILE:ocean_vgrid.nc,interfaces=zeta" ! default = "WOA09"
! FNC1:string - FNC1:dz_min,H_total,power,precision
! HYBRID:string - read from a file. The string specifies
! the filename and two variable names, separated
! by a comma or space, for sigma-2 and dz. e.g.
! HYBRID:vgrid.nc,sigma2,dz
REGRIDDING_ANSWER_DATE = 99991231 ! default = 20181231
! The vintage of the expressions and order of arithmetic to use for regridding.
! 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.
! by a comma or space, for sigma-2 and dz.
! e.g. HYBRID:vgrid.nc,sigma2,dz
! HYBRID_3D:string - read from a file. The string specifies
! the filename and two 3D variable names, separated
! by a comma or space, for sigma-2 and dz. The
! latter can be FNC1:string which is used everywhere.
! e.g. HYBRID_3D:vgrid.nc,sigma2,dz
! HYBRID_MAP:string - read from a file. The string specifies
! the filename and three variable names, separated
! by a comma or space, for map, sigma-2 and dz.
! Map is a spatial index array with, maxval(map)=N,
! and the others are 2D arrays containing N profiles.
! Map typically contains integer values, but it can
! contain real values, I+w, which imply using
! the weighted sum of profiles I and I+1.
! Dz can be FNC1:string which is used everywhere.
! e.g. HYBRID_MAP:vgrid.nc,map,sigma2,dz
DIAG_COORD_DEF_RHO2 = "RFNC1:76,999.5,1020.,1034.1,3.1,1041.,0.002" ! default = "WOA09"
! Determines how to specify the coordinate resolution. Valid options are:
! PARAM - use the vector-parameter DIAG_COORD_RES_RHO2
Expand All @@ -277,8 +302,23 @@ DIAG_COORD_DEF_RHO2 = "RFNC1:76,999.5,1020.,1034.1,3.1,1041.,0.002" ! default =
! FNC1:string - FNC1:dz_min,H_total,power,precision
! HYBRID:string - read from a file. The string specifies
! the filename and two variable names, separated
! by a comma or space, for sigma-2 and dz. e.g.
! HYBRID:vgrid.nc,sigma2,dz
! by a comma or space, for sigma-2 and dz.
! e.g. HYBRID:vgrid.nc,sigma2,dz
! HYBRID_3D:string - read from a file. The string specifies
! the filename and two 3D variable names, separated
! by a comma or space, for sigma-2 and dz. The
! latter can be FNC1:string which is used everywhere.
! e.g. HYBRID_3D:vgrid.nc,sigma2,dz
! HYBRID_MAP:string - read from a file. The string specifies
! the filename and three variable names, separated
! by a comma or space, for map, sigma-2 and dz.
! Map is a spatial index array with, maxval(map)=N,
! and the others are 2D arrays containing N profiles.
! Map typically contains integer values, but it can
! contain real values, I+w, which imply using
! the weighted sum of profiles I and I+1.
! Dz can be FNC1:string which is used everywhere.
! e.g. HYBRID_MAP:vgrid.nc,map,sigma2,dz

! === module MOM_MEKE ===
USE_MEKE = True ! [Boolean] default = False
Expand Down Expand Up @@ -370,11 +410,6 @@ FGNV_FILTER_SCALE = 0.1 ! [nondim] default = 1.0
! streamfunction formulation.

! === module MOM_dynamics_split_RK2 ===
VISC_REM_BUG = False ! [Boolean] default = True
! If true, visc_rem_[uv] in split mode is incorrectly calculated or accounted
! for in two places. This parameter controls the defaults of two individual
! flags, VISC_REM_TIMESTEP_BUG in MOM_dynamics_split_RK2(b) and
! VISC_REM_BT_WEIGHT_BUG in MOM_barotropic.

! === module MOM_continuity_PPM ===
ETA_TOLERANCE = 1.0E-06 ! [m] default = 3.75E-09
Expand All @@ -397,6 +432,9 @@ CORIOLIS_SCHEME = "SADOURNY75_ENSTRO" ! default = "SADOURNY75_ENERGY"
! ARAKAWA_LAMB81 - Arakawa & Lamb, 1981; En. + Enst.
! ARAKAWA_LAMB_BLEND - A blend of Arakawa & Lamb with
! Arakawa & Hsu and Sadourny energy
! WENOVI5TH_PV_ENSTRO - 5th-order WENO PV enstrophy
! WENOVI3RD_PV_ENSTRO - 3rd-order WENO PV enstrophy
! WENOVI7TH_PV_ENSTRO - 7th-order WENO PV enstrophy
BOUND_CORIOLIS = True ! [Boolean] default = False
! If true, the Coriolis terms at u-points are bounded by the four estimates of
! (f+rv)v from the four neighboring v-points, and similarly at v-points. This
Expand All @@ -409,6 +447,10 @@ 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 = False ! [Boolean] default = True
! 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.

! === module MOM_Zanna_Bolton ===

Expand All @@ -428,26 +470,13 @@ AH_VEL_SCALE = 0.01 ! [m s-1] default = 0.0
! scaled viscosity, the Smagorinsky and Leith viscosities, and AH.
SMAGORINSKY_AH = True ! [Boolean] default = False
! If true, use a biharmonic Smagorinsky nonlinear eddy viscosity.
BACKSCATTER_UNDERBOUND = False ! [Boolean] default = True
! If true, the bounds on the biharmonic viscosity are allowed to increase where
! the Laplacian viscosity is negative (due to backscatter parameterizations)
! beyond the largest timestep-dependent stable values of biharmonic viscosity
! when no Laplacian viscosity is applied. The default is true for historical
! reasons, but this option probably should not be used because it can contribute
! to numerical instabilities.
SMAG_BI_CONST = 0.06 ! [nondim] default = 0.0
! The nondimensional biharmonic Smagorinsky constant, typically 0.015 - 0.06.
FRICTWORK_BUG = False ! [Boolean] default = True
! If true, retain an answer-changing bug in calculating the FrictWork, which
! cancels the h in thickness flux and the h at velocity point. This isnot
! recommended.

! === module MOM_vert_friction ===
HMIX_FIXED = 0.5 ! [m]
! The prescribed depth over which the near-surface viscosity and diffusivity are
! elevated when the bulk mixed layer is not used.
MAXVEL = 6.0 ! [m s-1] default = 3.0E+08
! The maximum velocity allowed before the velocity components are truncated.
CFL_TRUNCATE_RAMP_TIME = 7200.0 ! [s] default = 0.0
! The time over which the CFL truncation value is ramped up at the beginning of
! the run.
Expand Down Expand Up @@ -564,11 +593,6 @@ USE_LOTW_BBL_DIFFUSIVITY = True ! [Boolean] default = False
! If true, uses a simple, imprecise but non-coordinate dependent, model of BBL
! mixing diffusivity based on Law of the Wall. Otherwise, uses the original BBL
! scheme.
LOTW_BBL_ANSWER_DATE = 99991231 ! default = 20190101
! The vintage of the order of arithmetic and expressions in the LOTW_BBL
! 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
Expand Down Expand Up @@ -648,10 +672,10 @@ EPBL_MSTAR_SCHEME = "OM4" ! default = "CONSTANT"
! REICHL_H18 - Use the scheme documented in Reichl & Hallberg, 2018.
MSTAR_CAP = 1.25 ! [nondim] default = -1.0
! If this value is positive, it sets the maximum value of mstar allowed in ePBL.
! (This is not used if EPBL_MSTAR_SCHEME = CONSTANT).
! (This is not used if EPBL_mstar_scheme = CONSTANT).
MSTAR2_COEF1 = 0.29 ! [nondim] default = 0.3
! Coefficient in computing mstar when rotation and stabilizing effects are both
! important (used if EPBL_MSTAR_SCHEME = OM4).
! important (used if EPBL_mstar_scheme = OM4).
MSTAR2_COEF2 = 0.152 ! [nondim] default = 0.085
! Coefficient in computing mstar when only rotation limits the total mixing
! (used if EPBL_MSTAR_SCHEME = OM4)
Expand Down Expand Up @@ -721,8 +745,15 @@ MAX_TR_DIFFUSION_CFL = 2.0 ! [nondim] default = -1.0
! This module implements neutral diffusion of tracers
USE_NEUTRAL_DIFFUSION = True ! [Boolean] default = False
! If true, enables the neutral diffusion module.
NDIFF_ANSWER_DATE = 20240101 ! default = 99991231
! 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_sum_output ===
MAXVEL = 6.0 ! [m s-1] default = 3.0E+08
! The maximum velocity allowed before the velocity components are truncated.

! === module ocean_stochastics_init ===

Expand Down Expand Up @@ -769,3 +800,11 @@ USE_RIGID_SEA_ICE = True ! [Boolean] default = False
SEA_ICE_RIGID_MASS = 100.0 ! [kg m-2] default = 1000.0
! The mass of sea-ice per unit area at which the sea-ice starts to exhibit
! rigidity
WAVE_INTERFACE_ANSWER_DATE = 20221231 ! default = 99991231
! 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,
! while higher values use updated and more robust forms of the same expressions:
! < 20230101 - Original answers for wave interface routines
! >= 20230101 - More robust expressions for Update_Stokes_Drift
! >= 20230102 - More robust expressions for get_StokesSL_LiFoxKemper
! >= 20230103 - More robust expressions for ust_2_u10_coare3p5
4 changes: 2 additions & 2 deletions config.yaml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -16,9 +16,9 @@ sync:
#Model software version
modules:
use:
- /g/data/vk83/modules
- /g/data/vk83/prerelease/modules
load:
- access-om3/2025.08.001
- access-om3/pr164-1

queue: normal
ncpus: 240
Expand Down
Loading
Loading