Skip to content

cam6_3_028: Science Updates and bug fixes for SE dycore. #401

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
jtruesdal merged 14 commits into from
Aug 19, 2021
Merged

cam6_3_028: Science Updates and bug fixes for SE dycore. #401

Show file tree
Hide file tree
Changes from 10 commits
Commits
Show all changes
14 commits
Select commit Hold shift + click to select a range
f06ba04
Full-physics physics-dynamics coupling with MPAS dynamical core
jtruesdal Jul 8, 2021
9d9e8f6
Changelog testing update, se_hypervis_subcycle_sponge mod for WACCM-X…
jtruesdal Jul 12, 2021
bba0979
updates to ChangeLog (Had wrong description)
jtruesdal Jul 15, 2021
e639045
move del4 sponge layer parameters to namelist
PeterHjortLauritzen Jul 19, 2021
ecdb4ca
bug in previous commit
PeterHjortLauritzen Jul 20, 2021
6f759fc
Generalize sponge layer setup and fix hard coded layer settings.
jtruesdal Jul 20, 2021
ee0edee
updates requested as part of PR#401
jtruesdal Jul 21, 2021
ef544d0
PR request to rename nu_s/nu_s_lev to nu_t/nu_t_lev for clarity
jtruesdal Jul 21, 2021
b5a6e0d
bring science updates for CESM2.2 from cam6_3_025 to cam_development …
jtruesdal Jul 21, 2021
3db5d65
ChangeLog merge with cam6_3_027
jtruesdal Jul 21, 2021
5c88161
PR requests
jtruesdal Jul 23, 2021
97f989b
i0 format suggestion by Steve
jtruesdal Aug 19, 2021
a8dcb53
update Changelog with climate validation runs
jtruesdal Aug 19, 2021
de2e10f
update Changelog for tag
jtruesdal Aug 19, 2021
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
6 changes: 3 additions & 3 deletions bld/build-namelist
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3619,6 +3619,9 @@ if ($dyn =~ /se/) {
se_nu_div
se_nu_p
se_nu_top
se_sponge_del4_nu_fac
se_sponge_del4_nu_div_fac
se_sponge_del4_lev
se_qsplit
se_rsplit
se_statediag_numtrac
Expand All @@ -3634,9 +3637,6 @@ if ($dyn =~ /se/) {
se_kmin_jet
se_kmax_jet
se_phys_dyn_cp
se_raytau0
se_raykrange
se_rayk0
se_molecular_diff
);

Expand Down
57 changes: 24 additions & 33 deletions bld/namelist_files/namelist_defaults_cam.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -649,11 +649,11 @@
<pgwv waccmx="1">18</pgwv>

<!-- gravity wave settings -->
<effgw_beres_dp>0.1D0</effgw_beres_dp>
<effgw_beres_dp>0.7D0</effgw_beres_dp>
<effgw_beres_dp hgrid="0.9x1.25" >0.4D0</effgw_beres_dp>
<effgw_beres_dp waccm_phys="0" >0.55D0</effgw_beres_dp>
<effgw_beres_dp waccm_phys="1" hgrid="0.9x1.25" >0.5D0</effgw_beres_dp>
<effgw_beres_dp waccm_phys="1" hgrid="ne30np4" nlev="70" >0.325D0</effgw_beres_dp>
<effgw_beres_dp waccm_phys="1" hgrid="ne30np4" nlev="70" >0.5D0</effgw_beres_dp>
<effgw_beres_dp waccm_phys="1" hgrid="ne30np4" nlev="110" >0.7D0</effgw_beres_dp>
<effgw_beres_dp chem="trop_strat_mam4_vbs" hgrid="0.9x1.25" >0.5D0</effgw_beres_dp>
<effgw_beres_dp chem="trop_strat_mam4_vbsext" hgrid="0.9x1.25">0.5D0</effgw_beres_dp>
Expand Down Expand Up @@ -1996,6 +1996,11 @@
<cld_macmic_num_steps microphys="mg2" clubb_sgs="1" dtime="600" > 1 </cld_macmic_num_steps>
<cld_macmic_num_steps microphys="mg2" clubb_sgs="1" dtime="300" > 1 </cld_macmic_num_steps>
<cld_macmic_num_steps microphys="mg3" clubb_sgs="1"> 3 </cld_macmic_num_steps>
<!-- hi-res and var-res support for ne120, ARCTIC, ARCTCIGRIS, CONUS -->
<cld_macmic_num_steps microphys="mg2" clubb_sgs="1" dtime="450" > 2 </cld_macmic_num_steps>
<cld_macmic_num_steps microphys="mg2" clubb_sgs="1" dtime="225" > 1 </cld_macmic_num_steps>
<cld_macmic_num_steps microphys="mg3" clubb_sgs="1" dtime="450" > 2 </cld_macmic_num_steps>
<cld_macmic_num_steps microphys="mg3" clubb_sgs="1" dtime="225" > 1 </cld_macmic_num_steps>

<!-- Ice nucleation -->
<nucleate_ice_subgrid >1.0D0</nucleate_ice_subgrid>
Expand Down Expand Up @@ -2658,70 +2663,57 @@
<!-- Defaults for SE -->
<!-- ================================================================== -->

<se_fine_ne > 0 </se_fine_ne>
<se_fine_ne hgrid="ne0np4CONUS.ne30x8" > 240 </se_fine_ne>

<se_ftype> 2 </se_ftype>

<se_hypervis_dynamic_ref_state> .false. </se_hypervis_dynamic_ref_state>
<se_lcp_moist > .true. </se_lcp_moist>

<se_large_Courant_incr > .true. </se_large_Courant_incr>

<se_hypervis_power> 0 </se_hypervis_power>

<se_hypervis_scaling se_refined_mesh="1" hypervis_type="tensor" >3.0D0 </se_hypervis_scaling>
<se_hypervis_scaling se_refined_mesh="1" hypervis_type="tensor" >3.22D0 </se_hypervis_scaling>

<se_hypervis_subcycle > 3 </se_hypervis_subcycle>
<se_hypervis_subcycle waccm_phys="1" > 1 </se_hypervis_subcycle>
<se_hypervis_subcycle waccm_phys="1" > 2 </se_hypervis_subcycle>
<se_hypervis_subcycle hgrid="ne16np4" > 4 </se_hypervis_subcycle>
<se_hypervis_subcycle se_refined_mesh="1" > 2 </se_hypervis_subcycle>
<se_hypervis_subcycle se_refined_mesh="1" > 3 </se_hypervis_subcycle>
<se_hypervis_subcycle hgrid="ne0np4CONUS.ne30x8" waccm_phys="1" > 1 </se_hypervis_subcycle>

<se_hypervis_subcycle_sponge > 1 </se_hypervis_subcycle_sponge>
<se_hypervis_subcycle_sponge waccmx="1" > 2 </se_hypervis_subcycle_sponge>
<se_hypervis_subcycle_sponge waccm_phys="1" > 1 </se_hypervis_subcycle_sponge>
<se_hypervis_subcycle_sponge waccmx="1" > 5 </se_hypervis_subcycle_sponge>
<se_hypervis_subcycle_sponge hgrid="ne120np4" > 4 </se_hypervis_subcycle_sponge>
<se_hypervis_subcycle_sponge se_refined_mesh="1" > 4 </se_hypervis_subcycle_sponge>
<se_hypervis_subcycle_sponge hgrid="ne0np4CONUS.ne30x8" waccm_phys="1"> 7 </se_hypervis_subcycle_sponge>
<se_hypervis_subcycle_sponge se_refined_mesh="1" > 2 </se_hypervis_subcycle_sponge>
<se_hypervis_subcycle_sponge hgrid="ne0np4CONUS.ne30x8" waccm_phys="1"> 4 </se_hypervis_subcycle_sponge>

<se_hypervis_subcycle_q>1 </se_hypervis_subcycle_q>
<se_hypervis_subcycle_q hgrid="ne16np4">2 </se_hypervis_subcycle_q>

<se_limiter_option> 8 </se_limiter_option>

<se_fine_ne> 0 </se_fine_ne>
<se_hypervis_power> 0 </se_hypervis_power>

<se_max_hypervis_courant > 1.0e99 </se_max_hypervis_courant>
<se_max_hypervis_courant se_refined_mesh="1" hypervis_type="scalar" > 1.9 </se_max_hypervis_courant>

<se_phys_dyn_cp>1</se_phys_dyn_cp>

<se_nu>-1</se_nu>
<se_nu se_refined_mesh="1" hypervis_type="scalar" > 1.0e13 </se_nu>
<se_nu waccmx="1">5.e15</se_nu>

<se_nu_div> -1 </se_nu_div>
<se_nu_div se_refined_mesh="1" hypervis_type="scalar" > 1.5625e13 </se_nu_div>
<se_nu_div waccmx="1"> 10.e15 </se_nu_div>

<se_nu_p>-1 </se_nu_p>
<se_nu_p se_refined_mesh="1" hypervis_type="scalar" > 1.5625e13 </se_nu_p>

<se_nu_top > 5.0e5 </se_nu_top>
<se_nu_top se_refined_mesh="1"> 2.0e5 </se_nu_top>
<se_nu_top > 1.25e5 </se_nu_top>
<se_nu_top waccmx="1" > 1.0e6 </se_nu_top>
<!-- WACCM uses molecular/thermal diffusion in sponge -->
<se_nu_top waccm_phys="1" > 0.0 </se_nu_top>

<se_raytau0 > 0.0 </se_raytau0>
<se_raykrange > 0.5 </se_raykrange>
<se_rayk0 > 2 </se_rayk0>
<!-- suggested value for WACCM se_raytau0=0.02 -->
<se_raytau0 waccm_phys="1" > 0.0 </se_raytau0>
<se_raykrange waccm_phys="1" > 3 </se_raykrange>
<se_rayk0 waccm_phys="1" > 2 </se_rayk0>

<se_molecular_diff > 0.0 </se_molecular_diff>
<se_molecular_diff waccm_phys="1" waccmx="0">100.0 </se_molecular_diff>
<se_molecular_diff waccmx="1" > 1.0 </se_molecular_diff>

<se_sponge_del4_nu_fac > -1 </se_sponge_del4_nu_fac>
<se_sponge_del4_nu_div_fac> -1 </se_sponge_del4_nu_div_fac>
<se_sponge_del4_lev > -1 </se_sponge_del4_lev>

<se_qsplit > 1 </se_qsplit>

Expand All @@ -2738,14 +2730,13 @@
<se_nsplit hgrid="ne16np4" waccm_phys="1" waccmx="1" > 5 </se_nsplit>
<se_nsplit hgrid="ne16np4" waccm_phys="1" > 3 </se_nsplit>
<se_nsplit hgrid="ne30np4" waccm_phys="1" > 5 </se_nsplit>
<se_nsplit hgrid="ne30np4" waccm_phys="1" nlev="110"> 4 </se_nsplit>
<se_nsplit hgrid="ne0np4CONUS.ne30x8" waccm_phys="1" > 10 </se_nsplit>

<se_nsplit hgrid="ne0np4CONUS.ne30x8" waccm_phys="1"> 5 </se_nsplit>
<se_nsplit hgrid="ne0np4TESTONLY.ne5x4"> 7 </se_nsplit>

<se_rsplit > 3 </se_rsplit>
<se_rsplit waccm_phys="1" > 4 </se_rsplit>
<se_rsplit hgrid="ne0np4CONUS.ne30x8" waccm_phys="1"> 5 </se_rsplit>
<se_rsplit hgrid="ne0np4CONUS.ne30x8" waccm_phys="1"> 3 </se_rsplit>

<se_fvm_supercycling>-1</se_fvm_supercycling>
<se_fvm_supercycling_jet>-1</se_fvm_supercycling_jet>
Expand Down
60 changes: 38 additions & 22 deletions bld/namelist_files/namelist_definition.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -7210,28 +7210,6 @@ Number of hyperviscosity subcycles per dynamics timestep in sponge del2 diffusio
Default: Set by build-namelist
</entry>


<entry id="se_rayk0" type="integer" category="se"
group="dyn_se_inparm" valid_values="" >
Variable to specify the vertical index at which the
Rayleigh friction term is centered (the peak value).
Default: 2
</entry>

<entry id="se_raykrange" type="real" category="se"
group="dyn_se_inparm" valid_values="" >
Rayleigh friction parameter to determine the width of the profile. If set
to 0 then a width is chosen by the algorithm (see rayleigh_friction.F90).
Default: 0.5.
</entry>

<entry id="se_raytau0" type="real" category="se"
group="dyn_se_inparm" valid_values="" >
Rayleigh friction parameter to determine the approximate value of the decay
time (days) at model top. If 0.0 then no Rayleigh friction is applied.
Default: 0.
</entry>

<entry id="se_molecular_diff" type="real" category="se"
group="dyn_se_inparm" valid_values="" >
Used by SE dycore to apply sponge layer diffusion to u, v, and T for
Expand Down Expand Up @@ -7330,6 +7308,44 @@ Second-order viscosity applied only near the model top [m^2/s].
Default: Set by build-namelist.
</entry>

<entry id="se_sponge_del4_nu_fac" type="real" category="se"
group="dyn_se_inparm" valid_values="" >
Hyperviscosity coefficient se_nu [m^4/s] for u,v, T is increased to
se_nu_p*se_sponge_del4_nu_fac following a hyperbolic tangent function
centered around pressure at vertical index se_sponge_del4_lev:

0.5_r8*(1.0_r8+tanh(2.0_r8*log(pmid(se_sponge_del4_lev)/press)))

where press is pressure

If &lt; 0, se_sponge_del4_nu_fac is automatically set based on model top location.
Default: Set by build-namelist.
</entry>

<entry id="se_sponge_del4_nu_div_fac" type="integer" category="se"
group="dyn_se_inparm" valid_values="" >
Divergence damping hyperviscosity coefficient se_nu_div [m^4/s] for u,v is increased to
se_nu_p*se_sponge_del4_nu_fac following a hyperbolic tangent function
Copy link
Collaborator

@nusbaume nusbaume Jul 22, 2021

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I believe it should be se_nu_p*se_sponge_del4_nu_div_fac here.

Copy link
Collaborator Author

@jtruesdal jtruesdal Jul 23, 2021

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

corrected

centered around pressure at vertical index se_sponge_del4_lev:

0.5_r8*(1.0_r8+tanh(2.0_r8*log(pmid(se_sponge_del4_lev)/press)))

where press is pressure

If &lt; 0, se_sponge_del4_nu_div_fac is automatically set based on model top location.
Default: Set by build-namelist.
</entry>

<entry id="se_sponge_del4_lev" type="real" category="se"
group="dyn_se_inparm" valid_values="" >
Level index around which increased del4 damping is centered.

See se_sponge_del4_nu_fac and se_sponge_del4_nu_div_fac

If &lt; 0, se_sponge_del4_lev is automatically set based on model top location.
Default: Set by build-namelist.
</entry>

<entry id="se_hypervis_dynamic_ref_state" type="logical" category="se"
group="dyn_se_inparm" valid_values="" >
Hyperscosity for T and dp is applied to (T-Tref) and (dp-dp_ref) where
Expand Down
4 changes: 0 additions & 4 deletions cime_config/testdefs/testmods_dirs/cam/outfrq9s_wcm_ne30/user_nl_cam
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -4,9 +4,5 @@ nhtfrq=9,9,9,9,9,9,9,9,9
inithist='ENDOFRUN'

se_nsplit = 30
se_rayk0= 10
se_raykrange= 5
se_raytau0=0.002

state_debug_checks = .true.

Loading