Skip to content
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

User/mjh/matts dumbbell defaults #12

Open
wants to merge 9 commits into
base: dev/esmg
Choose a base branch
from
Open

User/mjh/matts dumbbell defaults #12

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
9 commits
Select commit Hold shift + click to select a range
5e0bc35
MOM_input cleanup
MJHarrison-GFDL May 13, 2018
7862db1
Fixes for create_bcs.py
MJHarrison-GFDL May 13, 2018
e05939e
Matt's suggested defaults for regional model
MJHarrison-GFDL May 13, 2018
5095fc1
Matt's logfiles
MJHarrison-GFDL May 13, 2018
012e975
update src pointers
MJHarrison-GFDL May 13, 2018
71fb26d
dumbbell diag table
MJHarrison-GFDL May 13, 2018
98ac455
fix create_bcs.py
MJHarrison-GFDL May 13, 2018
02d860a
new answers for proposed reference case
MJHarrison-GFDL May 13, 2018
bcf7b98
Merge branch 'user/mjh/matts_dumbbell_defaults' into dev/esmg
MJHarrison-GFDL May 13, 2018
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
12 changes: 6 additions & 6 deletions ocean_only/dumbbell/common/MOM_input
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -250,7 +250,7 @@ ETA_TOLERANCE = 1.0E-12 ! [m] default = 1.0E-09
! tolerance for SSH is 4 times this value. The default
! is 0.5*NK*ANGSTROM, and this should not be set less x
! than about 10^-15*MAXIMUM_DEPTH.
CORIOLIS_SCHEME = "SADOURNY75_ENERGY" ! default = "SADOURNY75_ENERGY"
CORIOLIS_SCHEME = "SADOURNY75_ENSTRO" ! default = "SADOURNY75_ENERGY"
! CORIOLIS_SCHEME selects the discretization for the
! Coriolis terms. Valid values are:
! SADOURNY75_ENERGY - Sadourny, 1975; energy cons.
Expand All @@ -271,17 +271,17 @@ BOUND_CORIOLIS = True ! [Boolean] default = False
! CORIOLIS_SCHEME set to SADOURNY75_ENERGY.

! === module MOM_hor_visc ===
LAPLACIAN = False ! [Boolean] default = False
LAPLACIAN = True ! [Boolean] default = False
! If true, use a Laplacian horizontal viscosity.
BIHARMONIC = True ! [Boolean] default = True
! If true, se a biharmonic horizontal viscosity.
! BIHARMONIC may be used with LAPLACIAN.
KH_VEL_SCALE = 0.003 ! [m s-1] default = 0.0
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 viscosity and KH.
SMAG_BI_CONST = 0.03 ! [nondim] default = 0.0
SMAG_BI_CONST = 0.03 ! [nondim] default = 0.0
! The nondimensional biharmonic Smagorinsky constant,
! typically 0.015 - 0.06.
BOUND_CORIOLIS_VEL = 10.0 ! [m s-1] default = 1.0
Expand Down Expand Up @@ -459,11 +459,11 @@ MAXCPU = 2.88E+04 ! [wall-clock seconds] default = -1.0
! factor of the number of processors used.

! === module MOM_main (MOM_driver) ===
DT_FORCING = 600.0 ! [s] default = 900.0
DT_FORCING = 300.0 ! [s] default = 900.0
! The time step for changing forcing, coupling with other
! components, or potentially writing certain diagnostics.
! The default value is given by DT.
DAYMAX = 5.0 ! [days]
DAYMAX = 10.0 ! [days]
! The final time of the whole simulation, in units of
! TIMEUNIT seconds. This also sets the potential end
! time of the present run segment if the end time is
Expand Down
9 changes: 5 additions & 4 deletions ocean_only/dumbbell/common/diag_table
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,19 +1,20 @@
"Dumbbell experiment"
1 1 1 0 0 0
"prog", 10,"minutes",1,"days","Time"
"prog", 5,"minutes",1,"days","Time"

#This is the field section of the diag_table.

# Prognostic Ocean fields:
#=========================

"ocean_model","u","u","prog","all",.false.,"none",2
"ocean_model","v","v","prog","all",.false.,"none",2
"ocean_model","ssh","ssh","prog","all",.false.,"none",2
"ocean_model","uav","u","prog","all",.false.,"none",2
"ocean_model","vav","v","prog","all",.false.,"none",2
"ocean_model","h","h","prog","all",.false.,"none",1
"ocean_model","p_surf","p_surf","prog","all",.false.,"none",1
"ocean_model","e","e","prog","all",.false.,"none",2
"ocean_model","SSH","SSH","prog","all",.false.,"none",2
"ocean_model","RV","RV","prog","all",.false.,"none",2
"ocean_model","PV","PV","prog","all",.false.,"none",2
"ocean_model","salt","salt","prog","all",.false.,"none",2


Expand Down
2 changes: 1 addition & 1 deletion ocean_only/dumbbell/z/MOM_override
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -7,4 +7,4 @@ REMAPPING_SCHEME = "PPM_IH4"
SPONGE = True
SPONGE_CONFIG = "DUMBBELL"
DUMBBELL_SPONGE_TIME_SCALE = 8.64e4
! #override DAYMAX = 30.

45 changes: 35 additions & 10 deletions ocean_only/dumbbell/z/MOM_parameter_doc.all
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -54,12 +54,12 @@ 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.
DT = 600.0 ! [s]
DT = 300.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 = 600.0 ! [s] default = 600.0
DT_THERM = 300.0 ! [s] default = 300.0
! The thermodynamic and tracer advection time step.
! Ideally DT_THERM should be an integer multiple of DT
! and less than the forcing or coupling time-step, unless
Expand Down Expand Up @@ -839,7 +839,7 @@ CORIOLIS_EN_DIS = True ! [Boolean] default = False
! If true, two estimates of the thickness fluxes are used
! to estimate the Coriolis term, and the one that
! dissipates energy relative to the other one is used.
CORIOLIS_SCHEME = "SADOURNY75_ENERGY" ! default = "SADOURNY75_ENERGY"
CORIOLIS_SCHEME = "SADOURNY75_ENSTRO" ! default = "SADOURNY75_ENERGY"
! CORIOLIS_SCHEME selects the discretization for the
! Coriolis terms. Valid values are:
! SADOURNY75_ENERGY - Sadourny, 1975; energy cons.
Expand All @@ -851,9 +851,9 @@ CORIOLIS_SCHEME = "SADOURNY75_ENERGY" ! default = "SADOURNY75_ENERGY"
BOUND_CORIOLIS = True ! [Boolean] default = False
! If true, the Coriolis terms at u-points are bounded by
! the four estimates of (f+rv)v from the four neighboring
! v-points, and similarly at v-points. This option is
! always effectively false with CORIOLIS_EN_DIS defined and
! CORIOLIS_SCHEME set to SADOURNY75_ENERGY.
! v-points, and similarly at v-points. This option would
! have no effect on the SADOURNY Coriolis scheme if it
! were possible to use centered difference thickness fluxes.
KE_SCHEME = "KE_ARAKAWA" ! default = "KE_ARAKAWA"
! KE_SCHEME selects the discretization for acceleration
! due to the kinetic energy gradient. Valid values are:
Expand Down Expand Up @@ -894,8 +894,33 @@ BOUNDARY_EXTRAPOLATION_PRESSURE = True ! [Boolean] default = True
! used as is used for the interior cells.

! === module MOM_hor_visc ===
LAPLACIAN = False ! [Boolean] default = False
LAPLACIAN = True ! [Boolean] default = False
! If true, use a Laplacian horizontal viscosity.
KH = 0.0 ! [m2 s-1] default = 0.0
! The background Laplacian horizontal viscosity.
KH_BG_MIN = 0.0 ! [m2 s-1] default = 0.0
! The minimum value allowed for Laplacian horizontal viscosity, KH.
KH_VEL_SCALE = 0.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 = 0.0 ! [m2 s-1] default = 0.0
! The amplitude of a latidutinally-dependent background
! viscosity of the form KH_SIN_LAT*(SIN(LAT)**KH_PWR_OF_SINE).
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.
MODIFIED_LEITH = False ! [Boolean] default = False
! If true, add a term to Leith viscosity which is
! proportional to the gradient of divergence.
BOUND_KH = True ! [Boolean] default = True
! If true, the Laplacian coefficient is locally limited
! to be stable.
BETTER_BOUND_KH = True ! [Boolean] default = True
! If true, the Laplacian coefficient is locally limited
! to be stable with a better bounding than just BOUND_KH.
BIHARMONIC = True ! [Boolean] default = True
! If true, use a biharmonic horizontal viscosity.
! BIHARMONIC may be used with LAPLACIAN.
Expand All @@ -918,7 +943,7 @@ BOUND_AH = True ! [Boolean] default = True
BETTER_BOUND_AH = True ! [Boolean] default = True
! If true, the biharmonic coefficient is locally limited
! to be stable with a better bounding than just BOUND_AH.
SMAG_BI_CONST = 0.015 ! [nondim] default = 0.0
SMAG_BI_CONST = 0.03 ! [nondim] default = 0.0
! The nondimensional biharmonic Smagorinsky constant,
! typically 0.015 - 0.06.
BOUND_CORIOLIS_BIHARM = True ! [Boolean] default = True
Expand Down Expand Up @@ -1578,11 +1603,11 @@ LA_DEPTH_RATIO = 0.04 ! [nondim] default = 0.04
! Lanmguir number calculation, where La = sqrt(ust/Stokes).

! === module MOM_main (MOM_driver) ===
DT_FORCING = 600.0 ! [s] default = 600.0
DT_FORCING = 300.0 ! [s] default = 300.0
! The time step for changing forcing, coupling with other
! components, or potentially writing certain diagnostics.
! The default value is given by DT.
DAYMAX = 30.0 ! [days]
DAYMAX = 10.0 ! [days]
! The final time of the whole simulation, in units of
! TIMEUNIT seconds. This also sets the potential end
! time of the present run segment if the end time is
Expand Down
28 changes: 22 additions & 6 deletions ocean_only/dumbbell/z/MOM_parameter_doc.short
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -4,7 +4,7 @@
USE_REGRIDDING = True ! [Boolean] default = False
! If True, use the ALE algorithm (regridding/remapping).
! If False, use the layered isopycnal algorithm.
DT = 600.0 ! [s]
DT = 300.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
Expand Down Expand Up @@ -322,22 +322,38 @@ CORIOLIS_EN_DIS = True ! [Boolean] default = False
! If true, two estimates of the thickness fluxes are used
! to estimate the Coriolis term, and the one that
! dissipates energy relative to the other one is used.
CORIOLIS_SCHEME = "SADOURNY75_ENSTRO" ! default = "SADOURNY75_ENERGY"
! CORIOLIS_SCHEME selects the discretization for the
! Coriolis terms. Valid values are:
! SADOURNY75_ENERGY - Sadourny, 1975; energy cons.
! ARAKAWA_HSU90 - Arakawa & Hsu, 1990
! SADOURNY75_ENSTRO - Sadourny, 1975; enstrophy cons.
! ARAKAWA_LAMB81 - Arakawa & Lamb, 1981; En. + Enst.
! ARAKAWA_LAMB_BLEND - A blend of Arakawa & Lamb with
! Arakawa & Hsu and Sadourny energy
BOUND_CORIOLIS = True ! [Boolean] default = False
! If true, the Coriolis terms at u-points are bounded by
! the four estimates of (f+rv)v from the four neighboring
! v-points, and similarly at v-points. This option is
! always effectively false with CORIOLIS_EN_DIS defined and
! CORIOLIS_SCHEME set to SADOURNY75_ENERGY.
! v-points, and similarly at v-points. This option would
! have no effect on the SADOURNY Coriolis scheme if it
! were possible to use centered difference thickness fluxes.

! === module MOM_PressureForce ===

! === module MOM_PressureForce_AFV ===

! === module MOM_hor_visc ===
LAPLACIAN = True ! [Boolean] default = False
! If true, use a Laplacian horizontal viscosity.
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.
SMAGORINSKY_AH = True ! [Boolean] default = False
! If true, use a biharmonic Smagorinsky nonlinear eddy
! viscosity.
SMAG_BI_CONST = 0.015 ! [nondim] default = 0.0
SMAG_BI_CONST = 0.03 ! [nondim] default = 0.0
! The nondimensional biharmonic Smagorinsky constant,
! typically 0.015 - 0.06.

Expand Down Expand Up @@ -466,7 +482,7 @@ WIND_CONFIG = "zero" !
! === module MOM_restart ===

! === module MOM_main (MOM_driver) ===
DAYMAX = 30.0 ! [days]
DAYMAX = 10.0 ! [days]
! The final time of the whole simulation, in units of
! TIMEUNIT seconds. This also sets the potential end
! time of the present run segment if the end time is
Expand Down
Loading