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

A quasi-monotone flux limiter for isopycnal diffusion (Redi) #53

Closed
wants to merge 5 commits into from
Closed

A quasi-monotone flux limiter for isopycnal diffusion (Redi) #53

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
831eda7
Add a quasi-monotone limiter for isopycnal diffusion
dengwirda Apr 28, 2023
cb595d2
Add a quasi-monotone limiter for isopycnal diffusion
dengwirda Apr 28, 2023
02a36bd
Add a quasi-monotone limiter for isopycnal diffusion
dengwirda May 1, 2023
dfe41e3
Fix errors in isopycnal diffusion flux terms
dengwirda May 5, 2023
b8c3c99
initialize redi_flux_scale
mark-petersen Sep 20, 2023
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
8 changes: 8 additions & 0 deletions components/mpas-ocean/src/Registry.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -326,6 +326,14 @@
description="If true, the N2 limiting is applied to the horizontal diffusion term"
possible_values=".true. or .false."
/>
<nml_option name="config_Redi_use_quasi_monotone_limiter" type="logical" default_value=".false."
description="If true, fluxes are reduced to prevent tracers from violating monotonicity. Cross-term fluxes are scaled toward zero to prevent tracers from under/overshooting the min/max values in adjacent cells and layers"
possible_values=".true. or .false."
/>
<nml_option name="config_Redi_quasi_monotone_safety_factor" type="real" default_value="1.0"
description="A safety factor applied to flux scaling when monotonicity is violated. Smaller values scale fluxes toward zero more aggressively."
possible_values="A value between 0 and 1"
/>
<nml_option name="config_Redi_min_layers_diag_terms" type="integer" default_value="6"
description="Redi diagonal terms (2 and 3) are turned off from layer 1 through config_Redi_min_layers_diag_terms-1, and on from config_Redi_min_layers_diag_terms to nVertLevels. The Redi diagonal terms are not guaranteed to produce bounded tracer fields, and in practice produce growing temperatures in a few columns with fewer than 5 vertical cells. Redi is meant for isopycnal mixing in the deep ocean, so not applying Redi diagonal terms in very shallow regions is an acceptable solution."
possible_values="any integer between 0 (all layers on) and nVertLevels (all layers off)"
Expand Down
Loading