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

Merged
mgduda merged 9 commits into from
Feb 13, 2026
Merged

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
9 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
c34e74d
Remove etp, etm, ewp, and ewm from the 'invariant' stream in the atmo…
mgduda Feb 11, 2026
0d434b5
Remove old epssm variables and commented-out code from mpas_atm_time_…
mgduda Feb 11, 2026
a18b64f
Report use of config_epssm in the &nhyd_model namlist group with log …
mgduda Feb 12, 2026
6f9441c
Deactivate epssm log writes in atm_compute_damping_coefs using MPAS_D…
mgduda Feb 12, 2026
dfc451e
Copy epssm arrays etp, etm, ewp, and ewm to the device when using Ope…
mgduda Feb 12, 2026
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
35 changes: 34 additions & 1 deletion src/core_atmosphere/Registry.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -226,7 +226,7 @@
description="Mix full $\theta$ and $u$ fields, or mix perturbation from intitial state"
possible_values=".true. or .false."/>

<nml_option name="config_epssm" type="real" default_value="0.1"
<nml_option name="config_epssm" type="real" default_value="0.0" in_defaults="false"
units="-"
description="Off-centering parameter for the vertically implicit acoustic integration"
possible_values="Positive real values"/>
Expand Down 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 @@ -1484,6 +1505,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