I am seeing unusual tiling in the precipitation field simulated by WRF aerosol microphysics scheme of Song and Zhang, 2011, i.e., mp_physics = 40, which is used together with cu_physics = 11. I have 3 domains in my simulations and tiling is happening in all three domains. I am using WRF V4.4. The Thompson aerosol-aware scheme (mp_physics = 28) doesn't show such tiling and the basic Morrison scheme (mp_physics = 10) also does not show unusual tiling. It is basically happening when I select the mp_physics = 40 option. I also ran mp_physics = 40 but with aercu_opt = 0 but still I see the tiled results in the precipitation field. The tiling is actually happening with only the resolved precipitation (RAINNC) field, the convective precipitation (RAINC) looks just fine (please see the attached figure). So definitely it is related to microphysics aerosol code. But I don't know how to deal with this issue. Any suggestions please?
You are using an out of date browser. It may not display this or other websites correctly.
You should upgrade or use an alternative browser.
You should upgrade or use an alternative browser.
Unusual tiling in precipitation field while using WRF aerosol-aware microphysics scheme
- Thread starter psagar
- Start date
-
- Tags
- microphysics wrf physics
Here is my namelist.input file:
&time_control
start_year = 2023, 2023, 2023,
start_month = 02, 02, 02,
start_day = 10, 10, 10,
start_hour = 00, 00, 00,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2023, 2023, 2023,
end_month = 02, 02, 02,
end_day = 28, 28, 28,
end_hour = 12, 12, 12,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 21600
input_from_file = .true.,.true.,.true.,
history_interval = 60, 60, 60,
frames_per_outfile = 1, 1, 1,
restart = .false.,
restart_interval = 11520
write_hist_at_0h_rst = .true.,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
bdy_inname = "wrfbdy_d<domain>_<date>",
iofields_filename = "myoutfields.txt","myoutfields.txt","myoutfields.txt",
auxinput1_inname = 'met_em.d<domain>.<date>'
auxinput4_inname = 'wrflowinp_d<domain>',
!auxinput5_inname = 'wrfchemi_<hr>z_d<domain>',
!auxinput6_inname = 'wrfbiochemi_d01',
!auxinput7_inname = 'wrffirechemi_d<domain>',
!auxinput8_inname = 'wrfchemi_gocart_bg_d<domain>',
!auxinput12_inname = 'wrf_chem_input',
!auxinput13_inname = 'wrfchemv_d<domain>'
!history_outname = '/project/k1090/sagar/WRF_edit_fractions_nest_ssm_jja/wrfchem-3.8.1/test/em_real/large_files/wrfout_d<domain>_<date>',
auxinput4_interval = 360, 360, 360,
!auxinput5_interval_m = 60480, 60480, 60480,
!auxinput7_interval_m = 1440, 1440, 1440,
!auxinput8_interval_m = 1440, 1440, 1440,
io_form_auxinput2 = 2,
io_form_auxinput4 = 2, !sst data input
!io_form_auxinput5 = 2, !anthropogenic-emission input netcdf
!io_form_auxinput6 = 0, !biogenic emissions input netcdf
!io_form_auxinput7 = 0, !biomass burning input netcdf
!io_form_auxinput8 = 0, !gocart background emissions input netcdf
!io_form_auxinput12 = 0, !set to use previous simulation data to initilize chemistry
!io_form_auxinput13 = 0, !volcanic emissions input netcdf
!io_form_auxinput14 = 0, !aircraft emissions input netcdf
!io_form_auxinput15 = 0, !co2 or ghg gases emissions input netcdf only for chem_opt = 16 or 17
!io_form_auxinput = 0,
debug_level = 0,
/
&domains
time_step = 24,
time_step_fract_num = 0,
time_step_fract_den = 1,
!use_adaptive_time_step = .false.,
step_to_output_time = .true.,
target_cfl = 1.28, 1.28, 1.28,
!target_hcfl = 0.9, 0.9, 0.9,
!max_step_increase_pct = 5, 5, 5,
!starting_time_step = 36, 12, 4,
!max_time_step = 54, 18, 6,
!min_time_step = 9, 3, 1,
adaptation_domain = 1,
time_step_fract_num = 0,
time_step_fract_den = 1,
reasonable_time_step_ratio = 20,
max_dom = 3,
e_we = 828, 718, 508,
e_sn = 570, 754, 502,
e_vert = 60, 60, 60,
auto_levels_opt = 2,
dzbot = 40.,
dzstretch_s = 1.3,
dzstretch_u = 1.1,
p_top_requested = 10000,
num_metgrid_levels = 138,
num_metgrid_soil_levels = 4,
dx = 8894.198, 2964.733, 988.2442,
dy = 8894.198, 2964.733, 988.2442,
num_metgrid_soil_levels = 4,
grid_id = 1, 2, 3,
parent_id = 0, 1, 2,
i_parent_start = 1, 504, 272,
j_parent_start = 1, 163, 194,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 0,
smooth_option = 1
zap_close_levels = 50
interp_type = 1
t_extrap_type = 2
force_sfc_in_vinterp = 0
use_levels_below_ground = .true.
use_surface = .true.
lagrange_order = 1
sfcp_to_sfcp = .true.,
/
&physics
mp_physics = 40, 40, 40,
!progn = 1, 1, 1,
ra_lw_physics = 4, 4, 4,
ra_sw_physics = 4, 4, 4,
radt = 9, 3, 1,
sf_sfclay_physics = 2, 2, 2,
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 2, 2, 2,
bldt = 9, 3, 1,
cu_physics = 11, 11, 11,
cugd_avedx = 1,
cu_diag = 0, 0, 0,
cudt = 0, 0, 0,
ishallow = 1,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 0, 0, 0,
mp_zero_out = 2,
mp_zero_out_thresh = 1.e-12,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
cu_rad_feedback = .true.,
sst_update = 1
aercu_opt = 2,
aercu_fct = 1,
no_src_types_cu = 10,
alevsiz_cu = 30,
!aer_opt = 3,
/
&fdda
/
&dynamics
rk_ord = 3,
w_damping = 1,
diff_opt = 1,
km_opt = 4,
base_temp = 290.
damp_opt = 0,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.01, 0.01, 0.01,
diff_6th_opt = 0,
diff_6th_factor = 0.12, 0.12, 0.12,
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
non_hydrostatic = .true., .true., .true.,
moist_adv_opt = 2, 2, 1,
scalar_adv_opt = 2, 2, 1,
chem_adv_opt = 2, 2, 1,
tke_adv_opt = 2, 2, 1,
time_step_sound = 4, 4, 4,
h_mom_adv_order = 5, 5, 5,
v_mom_adv_order = 3, 3, 3,
h_sca_adv_order = 5, 5, 5,
v_sca_adv_order = 3, 3, 3,
/
&bdy_control
multi_bdy_files = .true.,
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false., .false.,
nested = .false., .true., .true.,
polar = .false., .false., .false.
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
&time_control
start_year = 2023, 2023, 2023,
start_month = 02, 02, 02,
start_day = 10, 10, 10,
start_hour = 00, 00, 00,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2023, 2023, 2023,
end_month = 02, 02, 02,
end_day = 28, 28, 28,
end_hour = 12, 12, 12,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 21600
input_from_file = .true.,.true.,.true.,
history_interval = 60, 60, 60,
frames_per_outfile = 1, 1, 1,
restart = .false.,
restart_interval = 11520
write_hist_at_0h_rst = .true.,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
bdy_inname = "wrfbdy_d<domain>_<date>",
iofields_filename = "myoutfields.txt","myoutfields.txt","myoutfields.txt",
auxinput1_inname = 'met_em.d<domain>.<date>'
auxinput4_inname = 'wrflowinp_d<domain>',
!auxinput5_inname = 'wrfchemi_<hr>z_d<domain>',
!auxinput6_inname = 'wrfbiochemi_d01',
!auxinput7_inname = 'wrffirechemi_d<domain>',
!auxinput8_inname = 'wrfchemi_gocart_bg_d<domain>',
!auxinput12_inname = 'wrf_chem_input',
!auxinput13_inname = 'wrfchemv_d<domain>'
!history_outname = '/project/k1090/sagar/WRF_edit_fractions_nest_ssm_jja/wrfchem-3.8.1/test/em_real/large_files/wrfout_d<domain>_<date>',
auxinput4_interval = 360, 360, 360,
!auxinput5_interval_m = 60480, 60480, 60480,
!auxinput7_interval_m = 1440, 1440, 1440,
!auxinput8_interval_m = 1440, 1440, 1440,
io_form_auxinput2 = 2,
io_form_auxinput4 = 2, !sst data input
!io_form_auxinput5 = 2, !anthropogenic-emission input netcdf
!io_form_auxinput6 = 0, !biogenic emissions input netcdf
!io_form_auxinput7 = 0, !biomass burning input netcdf
!io_form_auxinput8 = 0, !gocart background emissions input netcdf
!io_form_auxinput12 = 0, !set to use previous simulation data to initilize chemistry
!io_form_auxinput13 = 0, !volcanic emissions input netcdf
!io_form_auxinput14 = 0, !aircraft emissions input netcdf
!io_form_auxinput15 = 0, !co2 or ghg gases emissions input netcdf only for chem_opt = 16 or 17
!io_form_auxinput = 0,
debug_level = 0,
/
&domains
time_step = 24,
time_step_fract_num = 0,
time_step_fract_den = 1,
!use_adaptive_time_step = .false.,
step_to_output_time = .true.,
target_cfl = 1.28, 1.28, 1.28,
!target_hcfl = 0.9, 0.9, 0.9,
!max_step_increase_pct = 5, 5, 5,
!starting_time_step = 36, 12, 4,
!max_time_step = 54, 18, 6,
!min_time_step = 9, 3, 1,
adaptation_domain = 1,
time_step_fract_num = 0,
time_step_fract_den = 1,
reasonable_time_step_ratio = 20,
max_dom = 3,
e_we = 828, 718, 508,
e_sn = 570, 754, 502,
e_vert = 60, 60, 60,
auto_levels_opt = 2,
dzbot = 40.,
dzstretch_s = 1.3,
dzstretch_u = 1.1,
p_top_requested = 10000,
num_metgrid_levels = 138,
num_metgrid_soil_levels = 4,
dx = 8894.198, 2964.733, 988.2442,
dy = 8894.198, 2964.733, 988.2442,
num_metgrid_soil_levels = 4,
grid_id = 1, 2, 3,
parent_id = 0, 1, 2,
i_parent_start = 1, 504, 272,
j_parent_start = 1, 163, 194,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 0,
smooth_option = 1
zap_close_levels = 50
interp_type = 1
t_extrap_type = 2
force_sfc_in_vinterp = 0
use_levels_below_ground = .true.
use_surface = .true.
lagrange_order = 1
sfcp_to_sfcp = .true.,
/
&physics
mp_physics = 40, 40, 40,
!progn = 1, 1, 1,
ra_lw_physics = 4, 4, 4,
ra_sw_physics = 4, 4, 4,
radt = 9, 3, 1,
sf_sfclay_physics = 2, 2, 2,
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 2, 2, 2,
bldt = 9, 3, 1,
cu_physics = 11, 11, 11,
cugd_avedx = 1,
cu_diag = 0, 0, 0,
cudt = 0, 0, 0,
ishallow = 1,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 0, 0, 0,
mp_zero_out = 2,
mp_zero_out_thresh = 1.e-12,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
cu_rad_feedback = .true.,
sst_update = 1
aercu_opt = 2,
aercu_fct = 1,
no_src_types_cu = 10,
alevsiz_cu = 30,
!aer_opt = 3,
/
&fdda
/
&dynamics
rk_ord = 3,
w_damping = 1,
diff_opt = 1,
km_opt = 4,
base_temp = 290.
damp_opt = 0,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.01, 0.01, 0.01,
diff_6th_opt = 0,
diff_6th_factor = 0.12, 0.12, 0.12,
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
non_hydrostatic = .true., .true., .true.,
moist_adv_opt = 2, 2, 1,
scalar_adv_opt = 2, 2, 1,
chem_adv_opt = 2, 2, 1,
tke_adv_opt = 2, 2, 1,
time_step_sound = 4, 4, 4,
h_mom_adv_order = 5, 5, 5,
v_mom_adv_order = 3, 3, 3,
h_sca_adv_order = 5, 5, 5,
v_sca_adv_order = 3, 3, 3,
/
&bdy_control
multi_bdy_files = .true.,
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false., .false.,
nested = .false., .true., .true.,
polar = .false., .false., .false.
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/