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

Jump to bottom

Add changes for variable epssm and reference sounding option #1131

Open
wants to merge 3 commits into
base: develop
Choose a base branch
from
Open

Add changes for variable epssm and reference sounding option #1131

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
f0b26e6
1) Added code to allow epssm to vary with height according to a speci…
Nov 28, 2023
1914709
Modifications to clean version 8 of MPAS coe:
Dec 6, 2023
24e589d
Revisions to the pull request that contain additions to the clean
Dec 28, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
48 changes: 47 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 @@ -251,7 +251,8 @@
units="-"
description="Coefficient for the divergent component of the Laplacian filter of momentum in the relaxation zone"
possible_values="Positive real values"/>
</nml_record>

</nml_record>

<nml_record name="damping" in_defaults="true">
<nml_option name="config_zd" type="real" default_value="22000.0"
Expand All @@ -264,6 +265,26 @@
description="Maximum w-damping coefficient at model top"
possible_values="0 $\leq$ config_xnutr $\leq$ 1"/>

<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.1"
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="10000."
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="20000."
units="m"
description="Height MSL of top of transition zone for epssm"
possible_values="Positive real values"/>

<nml_option name="config_mpas_cam_coef" type="real" default_value="0.0" in_defaults="false"
units="-"
description="Coefficient for scaling the 2nd-order horizontal mixing in the mpas_cam absorbing layer"
Expand Down Expand Up @@ -471,6 +492,10 @@
<var name="cf3"/>
<var name="zgrid"/>
<var name="rdzw"/>
<var name="etp"/>
<var name="etm"/>
<var name="ewp"/>
<var name="ewm"/>
<var name="dzu"/>
<var name="rdzu"/>
<var name="fzm"/>
Expand All @@ -487,6 +512,7 @@
<var name="deriv_two"/>
<var name="defc_a"/>
<var name="defc_b"/>
<var name="reference_sounding_switch"/>
#ifdef MPAS_CAM_DYCORE
<var name="cell_gradient_coef_x"/>
<var name="cell_gradient_coef_y"/>
Expand Down Expand Up @@ -607,6 +633,10 @@
<var name="cf3"/>
<var name="zgrid"/>
<var name="rdzw"/>
<var name="etp"/>
<var name="etm"/>
<var name="ewp"/>
<var name="ewm"/>
<var name="dzu"/>
<var name="rdzu"/>
<var name="fzm"/>
Expand All @@ -623,6 +653,7 @@
<var name="deriv_two"/>
<var name="defc_a"/>
<var name="defc_b"/>
<var name="reference_sounding_switch"/>
#ifdef MPAS_CAM_DYCORE
<var name="cell_gradient_coef_x"/>
<var name="cell_gradient_coef_y"/>
Expand Down Expand Up @@ -1353,6 +1384,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 Expand Up @@ -1400,6 +1443,9 @@
<var name="qv_init" type="real" dimensions="nVertLevels" units="kg kg^{-1}"
description="qv reference profile"/>

<var name="reference_sounding_switch" type="integer" dimensions="" units="unitless" default_value="1"
description="Switch to designate a reference sounding is used (1) or full variables (0)"/>

<!-- Space needed for advection -->
<var name="deriv_two" type="real" dimensions="FIFTEEN TWO nEdges" units="unitless"
description="weights for cell-centered second derivative, normal to edge, for transport scheme"/>
Expand Down