Skip to content

Commit

Permalink
UWM changes to namelist files for new clubb external
Browse files Browse the repository at this point in the history
  • Loading branch information
@Katetc
Katetc committed Jan 3, 2024
1 parent 00188f3 commit 1d5007e
Show file tree
Hide file tree
Showing 3 changed files with 107 additions and 0 deletions.
3 changes: 3 additions & 0 deletions bld/build-namelist
View file
Original file line number Diff line number Diff line change
Expand Up @@ -3310,6 +3310,7 @@ if ($clubb_sgs =~ /$TRUE/io) {
add_default($nl, 'clubb_gamma_coefb');
}

add_default($nl, 'clubb_bv_efold');
add_default($nl, 'clubb_C7');
add_default($nl, 'clubb_C7b');
add_default($nl, 'clubb_c_K1');
Expand Down Expand Up @@ -3374,6 +3375,8 @@ if ($clubb_sgs =~ /$TRUE/io) {
add_default($nl, 'clubb_tridiag_solve_method');
add_default($nl, 'clubb_up2_sfc_coef');
add_default($nl, 'clubb_wpxp_L_thresh');
add_default($nl, 'clubb_wpxp_Ri_exp');
add_default($nl, 'clubb_z_displace');

#CLUBB+MF options
add_default($nl, 'do_clubb_mf');
Expand Down
4 changes: 4 additions & 0 deletions bld/namelist_files/namelist_defaults_cam.xml
View file
Original file line number Diff line number Diff line change
Expand Up @@ -2019,6 +2019,7 @@

<!-- CLUBB options -->
<clubb_beta > 2.4 </clubb_beta>
<clubb_bv_efold > 5.0 </clubb_bv_efold>
<clubb_c1 > 1.0 </clubb_c1>
<clubb_c1b > 1.0 </clubb_c1b>
<clubb_c11 > 0.7D0 </clubb_c11>
Expand Down Expand Up @@ -2121,6 +2122,8 @@
<clubb_tridiag_solve_method > 1 </clubb_tridiag_solve_method>
<clubb_up2_sfc_coef > 2.0 </clubb_up2_sfc_coef>
<clubb_wpxp_L_thresh > 60.0 </clubb_wpxp_L_thresh>
<clubb_wpxp_Ri_exp > 0.5 </clubb_wpxp_Ri_exp>
<clubb_z_displace > 25.0 </clubb_z_displace>
<do_hb_above_clubb >.false. </do_hb_above_clubb>
<do_hb_above_clubb phys="cam_dev" >.true. </do_hb_above_clubb>
<do_hb_above_clubb phys="cam6" >.true. </do_hb_above_clubb>
Expand Down Expand Up @@ -2149,6 +2152,7 @@
<clubb_skw_max_mag silhs="1" > 10.0 </clubb_skw_max_mag>
<clubb_up2_sfc_coef silhs="1" > 4.0 </clubb_up2_sfc_coef>
<clubb_C_wp2_splat silhs="1" > 0.0 </clubb_C_wp2_splat>
<clubb_bv_efold silhs="1" > 5.0 </clubb_bv_efold>

<clubb_l_brunt_vaisala_freq_moist silhs="1" > .true. </clubb_l_brunt_vaisala_freq_moist>
<clubb_l_call_pdf_closure_twice silhs="1" > .false. </clubb_l_call_pdf_closure_twice>
Expand Down
100 changes: 100 additions & 0 deletions bld/namelist_files/namelist_definition.xml
View file
Original file line number Diff line number Diff line change
Expand Up @@ -3732,6 +3732,11 @@ air is supersaturated with respect to ice.
Plume widths for theta_l and rt
</entry>

<entry id="clubb_bv_efold" type="real" category="pblrad"
group="clubb_params_nl" valid_values="" >
E-folding parameter for mixed Brunt Vaisala Frequency
</entry>

<entry id="clubb_c1" type="real" category="pblrad"
group="clubb_params_nl" valid_values="" >
Limiting value of C1 when skewness of w (vertical velocity) is small in
Expand Down Expand Up @@ -3971,6 +3976,14 @@ Gaussian PDF, and also decreases the difference between the means of w from
each Gaussian.
</entry>

<entry id="clubb_iiPDF_type" type="integer" category="pblrad"
group="clubb_params_nl" valid_values="" >
Selected option for the two-component normal
(double Gaussian) PDF type to use for the w, rt,
and theta-l (or w, chi, and eta) portion of
CLUBB's multivariate, two-component PDF.
</entry>

<entry id="clubb_ipdf_call_placement" type="integer" category="pblrad"
group="clubb_params_nl" valid_values="1,2,3" >
Option for the placement of the call to CLUBB's PDF closure. The options include: ipdf_pre_advance_fields (1) calls the PDF closure before advancing prognostic fields. ipdf_post_advance_fields (2) calls after advancing prognostic fields, and ipdf_pre_post_advance_fields (3) calls both before and after advancing prognostic fields.
Expand All @@ -3994,13 +4007,34 @@ Flag to uses an alternate equation to calculate the Brunt-Vaisala frequency.
This equation calculates an in-cloud Brunt-Vaisala frequency.
</entry>

<entry id="clubb_l_C2_cloud_frac" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
Flag to use cloud fraction to adjust the value of the
turbulent dissipation coefficient, C2.
</entry>

<entry id="clubb_l_calc_thlp2_rad" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
Include the contribution of radiation to thlp2
</entry>

<entry id="clubb_l_calc_w_corr" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
Calculate the correlations between w and the hydrometeors
</entry>

<entry id="clubb_l_call_pdf_closure_twice" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
Flag to call CLUBB's PDF closure at both thermodynamic and momentum vertical
grid levels. When this flag is turned off, CLUBB's PDF closure is only called
on thermodynamic grid levels.
</entry>

<entry id="clubb_l_const_Nc_in_cloud" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
Use a constant cloud droplet conc. within cloud
</entry>

<entry id="clubb_l_damp_wp2_using_em" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
Flag to use a dissipation formula of -(2/3)*em/tau_zm, as in Bougeault (1981),
Expand All @@ -4022,6 +4056,16 @@ is turned off, Lscale is calculated first, and then dissipation time-scale tau
is calculated as tau = Lscale / sqrt(tke).
</entry>

<entry id="clubb_l_diagnose_correlations" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
Diagnose correlations instead of using fixed ones
</entry>

<entry id="clubb_l_diffuse_rtm_and_thlm" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
Diffuses rtm and thlm
</entry>

<entry id="clubb_l_do_expldiff_rtm_thlm" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
Explicit diffusion on temperature and moisture by CLUBB, in addition to CLUBB's
Expand All @@ -4038,6 +4082,11 @@ Flag to run CLUBB with E3SM settings.
Flag to relax clipping on wpxp in xm_wpxp_clipping_and_stats.
</entry>

<entry id="clubb_l_fix_w_chi_eta_correlations" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
Use a fixed correlation for s and t Mellor(chi/eta)
</entry>

<entry id="clubb_l_godunov_upwind_wpxp_ta" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
This flag determines whether we want to use an upwind differencing approximation
Expand Down Expand Up @@ -4132,6 +4181,11 @@ horizontal winds um and vm. When this flag is turned off, upwp and vpwp are
calculated by down-gradient diffusion.
</entry>

<entry id="clubb_l_prescribed_avg_deltaz" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
used in adj_low_res_nu. If .true., avg_deltaz = deltaz
</entry>

<entry id="clubb_l_rcm_supersat_adj" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
Flag to take any remaining supersaturation after CLUBB PDF call and add it to
Expand All @@ -4141,6 +4195,11 @@ levels and the momentum grid levels and variables are interpolated between the
two grid level types.
</entry>

<entry id="clubb_l_rtm_nudge" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
For rtm nudging
</entry>

<entry id="clubb_l_smooth_Heaviside_tau_wpxp" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
Flag to use smooth Heaviside 'Peskin' in computation of invrs_tau.
Expand All @@ -4153,13 +4212,24 @@ Use the standard discretization for the turbulent advection terms. Setting to
advance_wp2_wp3_module.F90 and in advance_xp2_xpyp_module.F90.
</entry>

<entry id="clubb_l_stability_correct_Kh_N2_zm" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
Divides Kh_N2_zm by a stability factor
</entry>

<entry id="clubb_l_stability_correct_tau_zm" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
Flag to use a stability corrected version of CLUBB's time scale (tau_zm). This
creates a time scale that provides stronger damping at altitudes where
Brunt-Vaisala frequency is large.
</entry>

<entry id="clubb_l_tke_aniso" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
For anisotropic turbulent kinetic energy, i.e.
TKE = 1/2 (u'^2 + v'^2 + w'^2)
</entry>

<entry id="clubb_l_trapezoidal_rule_zm" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
Flag that uses the trapezoidal rule to adjust fields calculated by CLUBB's PDF
Expand All @@ -4176,6 +4246,14 @@ adjacent vertical grid level. The clubb_l_trapezoidal_rule_zt flag applies this
adjustment to PDF fields calculated on thermodynamic vertical grid levels.
</entry>

<entry id="clubb_l_upwind_xm_ma" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
This flag determines whether we want to use an upwind
differencing approximation rather than a centered
differencing for turbulent or mean advection terms. It
affects rtm, thlm, sclrm, um and vm.
</entry>

<entry id="clubb_l_upwind_xpyp_ta" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
Flag to use "upwind" discretization in the turbulent advection term in the
Expand All @@ -4185,6 +4263,11 @@ potential temperature). When this flag is turned off, centered discretization
is used.
</entry>

<entry id="clubb_l_uv_nudge" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
For wind speed nudging.
</entry>

<entry id="clubb_l_use_C11_Richardson" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
Flag to calculate the value of CLUBB's C11 based on Richardson number, where
Expand All @@ -4207,6 +4290,13 @@ levels influence the amount of cloudiness and amount of cloud water in a
grid box.
</entry>

<entry id="clubb_l_use_precip_frac" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
Flag to use precipitation fraction in KK microphysics. The
precipitation fraction is automatically set to 1 when this
flag is turned off.
</entry>

<entry id="clubb_l_use_shear_Richardson" type="logical" category="pblrad"
group="clubb_params_nl" valid_values="" >
Flag to use shear in the calculation of Richardson number.
Expand Down Expand Up @@ -4292,6 +4382,16 @@ clubb_up2_sfc_coef increases the values of up2 and vp2 at the surface.
CLUBB tunable parameter - Lscale threshold: damp C6 and C7 (units: m)
</entry>

<entry id="clubb_wpxp_Ri_exp" type="real" category="pblrad"
group="clubb_params_nl" valid_values="" >
Exponent for Richardson number in calculation of invrs_tau_wpxp term
</entry>

<entry id="clubb_z_displace" type="real" category="pblrad"
group="clubb_params_nl" valid_values="" >
Displacement of log law profile above ground (units: m)
</entry>

<!-- for CLUBB+MF -->
<entry id="do_clubb_mf" type="logical" category="conv"
group="clubb_mf_nl" valid_values="" >
Expand Down

0 comments on commit 1d5007e

Please sign in to comment.