Skip to content

Epssm(z) generalization of the parameter epssm in the MPAS dycore #1381

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

Open
skamaroc wants to merge 4 commits into
base: master
Choose a base branch
from
+229 −37
Open

Epssm(z) generalization of the parameter epssm in the MPAS dycore #1381

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
d635aed
Variable epssm(zeta) introduction: epssm is an offcentering coeffici…
skamaroc Sep 30, 2025
14ae8f8
Variable epssm(zeta) introduction, step (2): Added code into mpas_atm…
skamaroc Oct 1, 2025
097f9fe
Variable epssm(zeta) introduction, step (3): modified code to set
skamaroc Oct 14, 2025
51d8389
Changed the default values of the variable_epssm configuration variables
skamaroc Oct 14, 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
37 changes: 37 additions & 0 deletions src/core_atmosphere/Registry.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -292,6 +292,27 @@
units="-"
description="Number of layers in which to apply Rayleigh damping on horizontal velocity at top of model; damping linearly ramps to zero by layer number from the top"
possible_values="Positive integer values"/>

<nml_option name="config_epssm_minimum" type="real" default_value="0.1"
units="-"
description="Value of epssm below transition zone"
possible_values="Positive real values between 0 and 1"/>

<nml_option name="config_epssm_maximum" type="real" default_value="0.5"
units="-"
description="Value of epssm above transition zone"
possible_values="Positive real values between 0 and 1"/>

<nml_option name="config_epssm_transition_bottom_z" type="real" default_value="30000."
units="m"
description="Height MSL of bottom of transition zone for epssm"
possible_values="Positive real values"/>

<nml_option name="config_epssm_transition_top_z" type="real" default_value="50000."
units="m"
description="Height MSL of top of transition zone for epssm"
possible_values="Positive real values"/>

</nml_record>

<nml_record name="limited_area" in_defaults="true">
Expand Down Expand Up @@ -495,6 +516,10 @@
<var name="defc_a"/>
<var name="defc_b"/>
<var name="coeffs_reconstruct"/>
<var name="etp"/>
<var name="etm"/>
<var name="ewp"/>
<var name="ewm"/>
#ifdef MPAS_CAM_DYCORE
<var name="cell_gradient_coef_x"/>
<var name="cell_gradient_coef_y"/>
Expand Down Expand Up @@ -1484,6 +1509,18 @@
<var name="rdzw" type="real" dimensions="nVertLevels" units="unitless"
description="Reciprocal dzw"/>

<var name="etp" type="real" dimensions="nVertLevels" units="unitless"
description="etp = (1+epssm(z))/2 at theta point"/>

<var name="etm" type="real" dimensions="nVertLevels" units="unitless"
description="etm = (1-epssm(z))/2 at theta point"/>

<var name="ewp" type="real" dimensions="nVertLevelsP1" units="unitless"
description="ewp = (1+epssm(z))/2 at w point"/>

<var name="ewm" type="real" dimensions="nVertLevelsP1" units="unitless"
description="ewm = (1-epssm(z))/2 at w point"/>

<var name="dzu" type="real" dimensions="nVertLevels" units="unitless"
description="d(zeta) at w levels"/>

Expand Down
Loading