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.

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

EAMxx: Change default value of lambda_high in SHOC to 0.08 #6797

Open
wants to merge 2 commits into
base: master
Choose a base branch
from
Open
Show file tree
Hide file tree
Changes from all commits
Commits
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
Original file line number Diff line number Diff line change
Expand Up @@ -827,7 +827,7 @@
<shoc_length_fac > 0.5D0 </shoc_length_fac>
<shoc_c_diag_3rd_mom > 7.0D0 </shoc_c_diag_3rd_mom>
<shoc_lambda_low > 0.001D0 </shoc_lambda_low>
<shoc_lambda_high > 0.04D0 </shoc_lambda_high>
<shoc_lambda_high > 0.08D0 </shoc_lambda_high>
<shoc_lambda_slope > 2.65D0 </shoc_lambda_slope>
<shoc_lambda_thresh > 0.02D0 </shoc_lambda_thresh>
<shoc_Ckh > 0.1D0 </shoc_Ckh>
Expand Down
2 changes: 1 addition & 1 deletion components/eam/src/physics/cam/shoc.F90
Original file line number Diff line number Diff line change
Expand Up @@ -59,7 +59,7 @@ module shoc
real(rtype) :: length_fac = 0.5_rtype ! Length scale factor
real(rtype) :: c_diag_3rd_mom = 7.0_rtype ! w3 factor
real(rtype) :: lambda_low = 0.001_rtype ! lowest value for stability correction
real(rtype) :: lambda_high = 0.04_rtype ! highest value for stability correction
real(rtype) :: lambda_high = 0.08_rtype ! highest value for stability correction
real(rtype) :: lambda_slope = 2.65_rtype ! stability correction slope
real(rtype) :: lambda_thresh = 0.02_rtype ! value to apply stability correction
real(rtype) :: Ckh = 0.1_rtype ! Eddy diffusivity coefficient for heat
Expand Down
2 changes: 1 addition & 1 deletion components/eamxx/cime_config/namelist_defaults_scream.xml
Original file line number Diff line number Diff line change
Expand Up @@ -230,7 +230,7 @@ be lost if SCREAM_HACK_XML is not enabled.
<enable_column_conservation_checks>false</enable_column_conservation_checks>
<check_flux_state_consistency>false</check_flux_state_consistency>
<lambda_low type="real" doc="minimum value of stability correction.">0.001</lambda_low>
<lambda_high type="real" doc="maximum value of stability correction.">0.04</lambda_high>
<lambda_high type="real" doc="maximum value of stability correction.">0.08</lambda_high>
<lambda_slope type="real" doc="slope of change from lambda_low to lambda_high.">2.65</lambda_slope>
<lambda_thresh type="real" doc="stability threshold for which to apply more stability correction.">0.02</lambda_thresh>
<thl2tune type="real" doc="Temperature variance tuning factor">1.0</thl2tune>
Expand Down
2 changes: 1 addition & 1 deletion components/eamxx/docs/old/physics/shoc/shoc_doc.tex
Original file line number Diff line number Diff line change
Expand Up @@ -217,7 +217,7 @@ \subsubsection{Eddy Diffusivities}
\end{equation}
%

Where $\lambda_{min} = 0.001$, $\lambda_{slope}$ = 0.35, and $N_{low}$ = 0.037. Here, $\lambda_{slope}$ is an adjustable parameter. $\lambda_{0}$ has a minimum threshold of 0.001 and a maximum threshold of 0.04.
Where $\lambda_{min} = 0.001$, $\lambda_{slope}$ = 0.35, and $N_{low}$ = 0.037. Here, $\lambda_{slope}$ is an adjustable parameter. $\lambda_{0}$ has a minimum threshold of 0.001 and a maximum threshold of 0.08.

\paragraph{Stable Boundary Layer}

Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -2140,7 +2140,7 @@ void isotropic_ts_host(Int nlev, Int shcol, Real* brunt_int, Real* tke,

// Hard code these runtime options for F90
const Real lambda_low = 0.001;
const Real lambda_high = 0.04;
const Real lambda_high = 0.08;
const Real lambda_slope = 2.65;
const Real lambda_thresh = 0.02;
SHF::isotropic_ts(team, nlev, lambda_low, lambda_high, lambda_slope, lambda_thresh,
Expand Down Expand Up @@ -2950,7 +2950,7 @@ void shoc_tke_host(Int shcol, Int nlev, Int nlevi, Real dtime, Real* wthv_sec, R

// Hardcode for F90 testing
const Real lambda_low = 0.001;
const Real lambda_high = 0.04;
const Real lambda_high = 0.08;
const Real lambda_slope = 2.65;
const Real lambda_thresh = 0.02;
const Real Ckh = 0.1;
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -41,7 +41,7 @@ atmosphere_processes:
shoc:
enable_column_conservation_checks: true
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -46,7 +46,7 @@ atmosphere_processes:
max_total_ni: 740.0e3
shoc:
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -37,7 +37,7 @@ atmosphere_processes:
shoc:
check_flux_state_consistency: true
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -37,7 +37,7 @@ atmosphere_processes:
max_total_ni: 740.0e3
shoc:
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -55,7 +55,7 @@ atmosphere_processes:
shoc:
check_flux_state_consistency: true
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -62,7 +62,7 @@ atmosphere_processes:
shoc:
enable_column_conservation_checks: true
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -43,7 +43,7 @@ atmosphere_processes:
shoc:
enable_column_conservation_checks: true
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -49,7 +49,7 @@ atmosphere_processes:
shoc:
check_flux_state_consistency: true
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -41,7 +41,7 @@ atmosphere_processes:
max_total_ni: 740.0e3
shoc:
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -41,7 +41,7 @@ atmosphere_processes:
do_prescribed_ccn: false
shoc:
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -28,7 +28,7 @@ atmosphere_processes:
do_prescribed_ccn: false
shoc:
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -16,7 +16,7 @@ atmosphere_processes:
max_total_ni: 740.0e3
shoc:
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -21,7 +21,7 @@ atmosphere_processes:
do_prescribed_ccn: false
shoc:
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -24,7 +24,7 @@ atmosphere_processes:
do_prescribed_ccn: false
shoc:
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -24,7 +24,7 @@ atmosphere_processes:
do_prescribed_ccn: false
shoc:
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -21,7 +21,7 @@ atmosphere_processes:
do_prescribed_ccn: false
shoc:
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -21,7 +21,7 @@ atmosphere_processes:
top_level_mam4xx: 6
shoc:
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -18,7 +18,7 @@ atmosphere_processes:
number_of_subcycles: ${MAC_MIC_SUBCYCLES}
shoc:
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -21,7 +21,7 @@ atmosphere_processes:
do_prescribed_ccn: false
shoc:
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -22,7 +22,7 @@ atmosphere_processes:
max_total_ni: 740.0e3
shoc:
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -23,7 +23,7 @@ atmosphere_processes:
source_pressure_file: vertical_remap.nc ## Only used in the case of STATIC_1D_VERTICAL_PROFILE
shoc:
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -24,7 +24,7 @@ atmosphere_processes:
source_pressure_file: vertical_remap.nc ## Only used in the case of STATIC_1D_VERTICAL_PROFILE
shoc:
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -15,7 +15,7 @@ atmosphere_processes:
max_total_ni: 720.0e3
shoc:
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
2 changes: 1 addition & 1 deletion components/eamxx/tests/single-process/shoc/input.yaml
Original file line number Diff line number Diff line change
Expand Up @@ -14,7 +14,7 @@ atmosphere_processes:
number_of_subcycles: ${NUM_SUBCYCLES}
compute_tendencies: [all]
lambda_low: 0.001
lambda_high: 0.04
lambda_high: 0.08
lambda_slope: 2.65
lambda_thresh: 0.02
thl2tune: 1.0
Expand Down
Loading