2024/10/05 04:02:33 -------------------------------------------------- Latest git log entry -------------------------------------------------- Branch name: HEAD (If no branch name is shown HEAD may be detached) commit 6c01610082f21f5cb13460de2de4d4970351abdc Merge: 771412d05 f518cd935 Author: bmg929 Date: Fri Oct 4 15:39:09 2024 -0500 Merge branch 'master' of github.com:larson-group/clubb A detailed git diff can be found at the end of this file -------------------------------------------------- Tunable Parameters: -------------------------------------------------- Parameter Value --------- ----- C1 = 1.00000000000000000000 C1b = 1.00000000000000000000 C1c = 1.00000000000000000000 C2rt = 2.00000000000000000000 C2thl = 2.00000000000000000000 C2rtthl = 2.00000000000000000000 C4 = 2.00000000000000000000 C_uu_shr = 0.40000000000000002220 C_uu_buoy = 0.29999999999999998890 C6rt = 2.00000000000000000000 C6rtb = 2.00000000000000000000 C6rtc = 1.00000000000000000000 C6thl = 2.00000000000000000000 C6thlb = 2.00000000000000000000 C6thlc = 1.00000000000000000000 C7 = 0.50000000000000000000 C7b = 0.50000000000000000000 C7c = 0.50000000000000000000 C8 = 0.50000000000000000000 C8b = 0.02000000000000000042 C10 = 3.29999999999999982236 C11 = 0.40000000000000002220 C11b = 0.40000000000000002220 C11c = 0.50000000000000000000 C12 = 1.00000000000000000000 C13 = 0.10000000000000000555 C14 = 1.00000000000000000000 C_wp2_pr_dfsn = 0.00000000000000000000 C_wp3_pr_tp = 0.00000000000000000000 C_wp3_pr_turb = 1.00000000000000000000 C_wp3_pr_dfsn = 0.00000000000000000000 C_wp2_splat = 0.00000000000000000000 C6rt_Lscale0 = 14.00000000000000000000 C6thl_Lscale0 = 14.00000000000000000000 C7_Lscale0 = 0.84999999999999997780 wpxp_L_thresh = 60.00000000000000000000 c_K = 0.20000000000000001110 c_K1 = 0.20000000000000001110 nu1 = 20.00000000000000000000 c_K2 = 0.02500000000000000139 nu2 = 1.00000000000000000000 c_K6 = 0.37500000000000000000 nu6 = 5.00000000000000000000 c_K8 = 5.00000000000000000000 nu8 = 20.00000000000000000000 c_K9 = 0.10000000000000000555 nu9 = 10.00000000000000000000 nu10 = 0.00000000000000000000 c_K_hm = 0.75000000000000000000 c_K_hmb = 0.75000000000000000000 K_hm_min_coef = 0.10000000000000000555 nu_hm = 1.50000000000000000000 slope_coef_spread_DG_means_w = 21.00000000000000000000 pdf_component_stdev_factor_w = 1.00000000000000000000 coef_spread_DG_means_rt = 0.80000000000000004441 coef_spread_DG_means_thl = 0.80000000000000004441 gamma_coef = 0.25000000000000000000 gamma_coefb = 0.25000000000000000000 gamma_coefc = 5.00000000000000000000 mu = 0.00100000000000000002 beta = 1.00000000000000000000 lmin_coef = 0.50000000000000000000 omicron = 0.50000000000000000000 zeta_vrnce_rat = 0.00000000000000000000 upsilon_precip_frac_rat = 0.55000000000000004441 lambda0_stability_coef = 0.02999999999999999889 mult_coef = 0.50000000000000000000 taumin = 90.00000000000000000000 taumax = 3600.00000000000000000000 Lscale_mu_coef = 2.00000000000000000000 Lscale_pert_coef = 0.10000000000000000555 alpha_corr = 0.14999999999999999445 Skw_denom_coef = 4.00000000000000000000 c_K10 = 1.00000000000000000000 c_K10h = 1.00000000000000000000 thlp2_rad_coef = 1.00000000000000000000 thlp2_rad_cloud_frac_thresh = 0.10000000000000000555 up2_sfc_coef = 4.00000000000000000000 Skw_max_mag = 10.00000000000000000000 C_invrs_tau_bkgnd = 1.10000000000000008882 C_invrs_tau_sfc = 0.10000000000000000555 C_invrs_tau_shear = 0.14999999999999999445 C_invrs_tau_N2 = 0.40000000000000002220 C_invrs_tau_N2_wp2 = 0.20000000000000001110 C_invrs_tau_N2_xp2 = 0.05000000000000000278 C_invrs_tau_N2_wpxp = 0.00000000000000000000 C_invrs_tau_N2_clear_wp3 = 1.00000000000000000000 C_invrs_tau_wpxp_Ri = 0.34999999999999997780 C_invrs_tau_wpxp_N2_thresh = 0.00033000000000000000 xp3_coef_base = 0.25000000000000000000 xp3_coef_slope = 0.01000000000000000021 altitude_threshold = 100.00000000000000000000 rtp2_clip_coef = 0.50000000000000000000 Cx_min = 0.33000000000000001554 Cx_max = 0.94999999999999995559 Richardson_num_min = 0.25000000000000000000 Richardson_num_max = 400.00000000000000000000 a3_coef_min = 1.00000000000000000000 a_const = 1.80000000000000004441 bv_efold = 5.00000000000000000000 wpxp_Ri_exp = 0.50000000000000000000 z_displace = 25.00000000000000000000 -------------------------------------------------- Preprocessing Directives: -------------------------------------------------- -DNETCDF enabled -DCOAMPS_MICRO enabled -DNR_SP enabled -DTUNER disabled -DSILHS enabled -Dnooverlap enabled -Draoffline enabled -DUSE_BUGSrad_ocast_random disabled -DBYTESWAP_IO disabled -------------------------------------------------- &model_setting -------------------------------------------------- runtype = lba nzmax = 128 grid_type = 2 deltaz = 75.00000000000000 zm_init = 0.000000000000000 zm_top = 28625.00000000000 zt_grid_fname = ../input/grid/deep_convection_128lev_27km_zt_grid.grd zm_grid_fname = day = 23 month = 2 year = 1999 lat_vals = -10.86670000000000 lon_vals = -61.95000000000000 sfc_elevation = 0.000000000000000 time_initial = 41400.00000000000 time_final = 63000.00000000000 dt_main = 60.00000000000000 dt_rad = 60.00000000000000 sfctype = 1 T_sfc = 288.0000000000000 p_sfc = 99129.99988000000 sens_ht = 0.000000000000000 latent_ht = 0.000000000000000 fcor = 3.7599999999999999E-005 T0 = 300.0000000000000 ts_nudge = 7200.000000000000 forcings_file_path = ../input/case_setups/lba_forcings/ l_t_dependent = T l_ignore_forcings = T l_modify_ic_with_cubic_int = F l_modify_bc_for_cnvg_test = F l_input_xpwp_sfc = F saturation_formula = 3 thlm_sponge_damp_settings%l_sponge_damping = F rtm_sponge_damp_settings%l_sponge_damping = F uv_sponge_damp_settings%l_sponge_damping = F wp2_sponge_damp_settings%l_sponge_damping = F wp3_sponge_damp_settings%l_sponge_damping = F up2_vp2_sponge_damp_settings%l_sponge_damping = F thlm_sponge_damp_settings%tau_sponge_damp_min = 60.00000000000000 thlm_sponge_damp_settings%tau_sponge_damp_max = 1800.000000000000 thlm_sponge_damp_settings%sponge_damp_depth = 0.2500000000000000 rtm_sponge_damp_settings%tau_sponge_damp_min = 60.00000000000000 rtm_sponge_damp_settings%tau_sponge_damp_max = 1800.000000000000 rtm_sponge_damp_settings%sponge_damp_depth = 0.2500000000000000 uv_sponge_damp_settings%tau_sponge_damp_min = 60.00000000000000 uv_sponge_damp_settings%tau_sponge_damp_max = 1800.000000000000 uv_sponge_damp_settings%sponge_damp_depth = 0.2500000000000000 wp2_sponge_damp_settings%tau_sponge_damp_min = 60.00000000000000 wp2_sponge_damp_settings%tau_sponge_damp_max = 1800.000000000000 wp2_sponge_damp_settings%sponge_damp_depth = 0.2500000000000000 wp3_sponge_damp_settings%tau_sponge_damp_min = 60.00000000000000 wp3_sponge_damp_settings%tau_sponge_damp_max = 1800.000000000000 wp3_sponge_damp_settings%sponge_damp_depth = 0.2500000000000000 up2_vp2_sponge_damp_settings%tau_sponge_damp_min = 60.00000000000000 up2_vp2_sponge_damp_settings%tau_sponge_damp_max = 1800.000000000000 up2_vp2_sponge_damp_settings%sponge_damp_depth = 0.2500000000000000 l_soil_veg = F l_restart = F l_input_fields = F restart_path_case = restart/lba time_restart = 0.000000000000000 debug_level = 2 sclr_dim = 0 edsclr_dim = 0 iisclr_thl = -2 iisclr_rt = -1 iisclr_CO2 = -1 sclr_tol = -------------------------------------------------- &stats_setting -------------------------------------------------- l_stats = T fname_prefix = lba stats_fmt = netcdf stats_tsamp = 60.00000000000000 stats_tout = 60.00000000000000 l_allow_small_stats_tout = F Constant flags: l_pos_def = F l_hole_fill = T l_gamma_Skw = T l_byteswap_io = F Constant tolerances [units] rt_tol [kg/kg] = 1.0000000000000000E-008 thl_tol [K] = 1.0000000000000000E-002 w_tol [m/s] = 2.0000000000000000E-002 -------------------------------------------------- µphysics_setting -------------------------------------------------- microphys_scheme = morrison l_cloud_sed = F sigma_g = 1.500000000000000 l_graupel = T l_hail = F l_seifert_beheng = F l_predict_Nc = F l_const_Nc_in_cloud = T specify_aerosol = morrison_lognormal l_subgrid_w = T l_arctic_nucl = F l_cloud_edge_activation = F l_fix_pgam = F l_in_cloud_Nc_diff = T l_var_covar_src = T l_upwind_diff_sed = F lh_microphys_type = interactive lh_num_samples = 8 lh_sequence_length = 1 lh_seed = 5489 l_fix_w_chi_eta_correlations = T l_silhs_KK_convergence_adj_mean = F importance_prob_thresh = 1.0000000000000000E-008 host_dx = 1000000.000000000 host_dy = 1000000.000000000 hmp2_ip_on_hmm2_ip_slope%rr = 0.000000000000000 hmp2_ip_on_hmm2_ip_slope%ri = 0.000000000000000 hmp2_ip_on_hmm2_ip_slope%rs = 0.000000000000000 hmp2_ip_on_hmm2_ip_slope%rg = 0.000000000000000 hmp2_ip_on_hmm2_ip_slope%Nr = 0.000000000000000 hmp2_ip_on_hmm2_ip_slope%Ni = 0.000000000000000 hmp2_ip_on_hmm2_ip_slope%Ni = 0.000000000000000 hmp2_ip_on_hmm2_ip_slope%Ng = 0.000000000000000 hmp2_ip_on_hmm2_ip_intrcpt%rr = 5.000000000000000 hmp2_ip_on_hmm2_ip_intrcpt%ri = 1.250000000000000 hmp2_ip_on_hmm2_ip_intrcpt%rs = 1.250000000000000 hmp2_ip_on_hmm2_ip_intrcpt%rg = 1.250000000000000 hmp2_ip_on_hmm2_ip_intrcpt%Nr = 2.500000000000000 hmp2_ip_on_hmm2_ip_intrcpt%Ni = 0.6250000000000000 hmp2_ip_on_hmm2_ip_intrcpt%Ni = 0.6250000000000000 hmp2_ip_on_hmm2_ip_intrcpt%Ng = 0.6250000000000000 Ncnp2_on_Ncnm2 = 3.0000000000000001E-003 C_evap = 0.8600000000000000 r_0 = 2.5000000000000001E-005 microphys_start_time = 0.000000000000000 Nc0_in_cloud = 100000000.0000000 ccnconst = 120.0000000000000 ccnexpnt = 0.4000000059604645 aer_rm1 = 1.1000000021965661E-008 aer_rm2 = 5.9999997858994902E-008 aer_n1 = 125000000.0000000 aer_n2 = 65000000.00000000 aer_sig1 = 1.200000047683716 aer_sig2 = 1.700000047683716 pgam_fixed = 5.000000000000000 precip_frac_calc_type = 2 -------------------------------------------------- &SILHS_setting -------------------------------------------------- cluster_allocation_strategy = 3 l_lh_importance_sampling = T l_Lscale_vert_avg = F l_lh_straight_mc = T l_lh_clustered_sampling = T l_rcm_in_cloud_k_lh_start = T l_random_k_lh_start = F l_max_overlap_in_cloud = T l_lh_instant_var_covar_src = T l_lh_limit_weights = T l_lh_var_frac = F l_lh_normalize_weights = T Correlation array (approximate); in cloud: 1.000 -0.600 0.090 0.090 0.558 0.544 0.212 0.207 0.212 0.207 0.212 0.207 -0.600 1.000 0.027 0.027 0.081 0.093 -0.026 0.087 0.019 0.012 0.000 0.000 0.090 0.027 1.000 0.340 0.223 0.217 0.106 0.156 0.000 0.000 0.000 0.000 0.090 0.027 0.340 1.000 0.000 0.000 0.464 0.320 0.168 0.232 0.000 0.000 0.558 0.081 0.223 0.000 1.000 0.728 0.000 0.000 0.116 0.114 0.226 0.223 0.544 0.093 0.217 0.000 0.728 1.000 0.113 0.111 0.000 0.000 0.222 0.219 0.212 -0.026 0.106 0.464 0.000 0.113 1.000 0.717 0.526 0.522 0.323 0.319 0.207 0.087 0.156 0.320 0.000 0.111 0.717 1.000 0.522 0.516 0.319 0.315 0.212 0.019 0.000 0.168 0.116 0.000 0.526 0.522 1.000 0.717 0.425 0.421 0.207 0.012 0.000 0.232 0.114 0.000 0.522 0.516 0.717 1.000 0.421 0.417 0.212 0.000 0.000 0.000 0.226 0.222 0.323 0.319 0.425 0.421 1.000 0.717 0.207 0.000 0.000 0.000 0.223 0.219 0.319 0.315 0.421 0.417 0.717 1.000 Correlation array (approximate); below cloud: 1.000 0.300 0.090 0.090 0.558 0.544 0.212 0.207 0.212 0.207 0.212 0.207 0.300 1.000 0.027 0.027 0.081 0.093 -0.026 0.087 0.019 0.012 0.000 0.000 0.090 0.027 1.000 0.340 0.223 0.217 0.106 0.156 0.000 0.000 0.000 0.000 0.090 0.027 0.340 1.000 0.000 0.000 0.464 0.320 0.168 0.232 0.000 0.000 0.558 0.081 0.223 0.000 1.000 0.728 0.000 0.000 0.116 0.114 0.226 0.223 0.544 0.093 0.217 0.000 0.728 1.000 0.113 0.111 0.000 0.000 0.222 0.219 0.212 -0.026 0.106 0.464 0.000 0.113 1.000 0.717 0.526 0.522 0.323 0.319 0.207 0.087 0.156 0.320 0.000 0.111 0.717 1.000 0.522 0.516 0.319 0.315 0.212 0.019 0.000 0.168 0.116 0.000 0.526 0.522 1.000 0.717 0.425 0.421 0.207 0.012 0.000 0.232 0.114 0.000 0.522 0.516 0.717 1.000 0.421 0.417 0.212 0.000 0.000 0.000 0.226 0.222 0.323 0.319 0.425 0.421 1.000 0.717 0.207 0.000 0.000 0.000 0.223 0.219 0.319 0.315 0.421 0.417 0.717 1.000 -------------------------------------------------- &radiation_setting -------------------------------------------------- rad_scheme = lba sol_const = 1367.000000000000 alvdr = 0.1000000000000000 alvdf = 0.1000000000000000 alndr = 0.1000000000000000 alndf = 0.1000000000000000 radiation_top = 50000.00000000000 F0 = 100.0000000000000 F1 = 20.00000000000000 kappa = 119.0000000000000 gc = 0.8600000000000000 omega = 0.9965000000000001 slr = 1.000000000000000 l_rad_above_cloud = F l_sw_radiation = F l_fix_cos_solar_zen = F l_use_default_std_atmosphere = T Fs_values = 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 cos_solar_zen_values = -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 cos_solar_zen_times = -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 -999.0000000000000 eff_drop_radius = 1.0000000000000001E-005 -------------------------------------------------- &configurable_clubb_flags_nl -------------------------------------------------- iiPDF_type = 1 ipdf_call_placement = 2 penta_solve_method = 2 tridiag_solve_method = 2 l_use_precip_frac = T l_predict_upwp_vpwp = T l_min_wp2_from_corr_wx = T l_min_xp2_from_corr_wx = T l_C2_cloud_frac = F l_diffuse_rtm_and_thlm = F l_stability_correct_Kh_N2_zm = F l_calc_thlp2_rad = T l_upwind_xpyp_ta = T l_upwind_xm_ma = T l_uv_nudge = T l_rtm_nudge = F l_tke_aniso = T l_vert_avg_closure = F l_trapezoidal_rule_zt = F l_trapezoidal_rule_zm = F l_call_pdf_closure_twice = F l_standard_term_ta = F l_partial_upwind_wp3 = F l_godunov_upwind_wpxp_ta = F l_godunov_upwind_xpyp_ta = F l_use_cloud_cover = F l_diagnose_correlations = F l_calc_w_corr = F l_const_Nc_in_cloud = T l_fix_w_chi_eta_correlations = T l_stability_correct_tau_zm = F l_damp_wp2_using_em = T l_do_expldiff_rtm_thlm = F l_Lscale_plume_centered = F l_diag_Lscale_from_tau = T l_use_C7_Richardson = T l_use_C11_Richardson = F l_use_shear_Richardson = F l_brunt_vaisala_freq_moist = F l_use_thvm_in_bv_freq = F l_rcm_supersat_adj = T l_damp_wp3_Skw_squared = T l_prescribed_avg_deltaz = F l_lmm_stepping = F l_e3sm_config = F l_vary_convect_depth = F l_use_tke_in_wp3_pr_turb_term = T l_use_tke_in_wp2_wp3_K_dfsn = F l_use_wp3_lim_with_smth_Heaviside = F l_smooth_Heaviside_tau_wpxp = T l_modify_limiters_for_cnvg_test = F l_enable_relaxed_clipping = F l_linearize_pbl_winds = F l_mono_flux_lim_thlm = T l_mono_flux_lim_rtm = T l_mono_flux_lim_um = T l_mono_flux_lim_vm = T l_mono_flux_lim_spikefix = T l_host_applies_sfc_fluxes = F -------------------------------------------------- git diff src/ -------------------------------------------------- diff --git a/input/tunable_parameters/configurable_model_flags.in b/input/tunable_parameters/configurable_model_flags.in index 83bcf05bf..f24a35b1b 100644 --- a/input/tunable_parameters/configurable_model_flags.in +++ b/input/tunable_parameters/configurable_model_flags.in @@ -1,8 +1,8 @@ &configurable_clubb_flags_nl iiPDF_type = 1 ipdf_call_placement = 2 -penta_solve_method = 1 -tridiag_solve_method = 1 +penta_solve_method = 2 +tridiag_solve_method = 2 saturation_formula = 3 l_godunov_upwind_wpxp_ta = .false., @@ -60,7 +60,7 @@ cluster_allocation_strategy = 3 l_lh_importance_sampling = .true. l_Lscale_vert_avg = .false. -l_lh_straight_mc = .false. +l_lh_straight_mc = true l_lh_clustered_sampling = .true. l_rcm_in_cloud_k_lh_start = .true. l_random_k_lh_start = .false.