Hi all,
I am running WRF-Chem 4.0.3 on 9km, 3km and 1km domain (two-way nesting) but surface values of aerosol concentrations (partly also gases) are in area of inner domains clearly 2-3 times lower in comparison to a run on 9km domain only. It is not problem with emissions but it seems there is a higher vertical mixing in 3km and 1km domains, pollutants are distributed more equally in PBL, while in 9km only run, they are concentrated more near surface. Any idea, tips to solute this discrepancy?
Thanks, Jan
My namelist:
&time_control
run_days = 15,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2015, 2015, 2015,
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 = 2015, 2015, 2015,
end_month = 02, 02, 02,
end_day = 25, 25, 25,
end_hour = 00, 00, 00,
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 = 20000, 20000, 20000,
restart = .false.,
restart_interval = 21600,
auxinput12_inname = 'wrf_chem_input',
io_form_auxinput5 = 2,
auxinput5_inname = 'wrfchemi_d<domain>_<date>',
io_form_auxinput5 = 2,
auxinput5_interval_m = 60,
frames_per_auxinput5 = 24,
auxinput6_inname = 'wrfbiochemi_d<domain>',
auxinput8_interval_m = 1440,
io_form_auxinput2 = 2,
io_form_auxinput6 = 2,
io_form_auxinput8 = 0,
io_form_auxinput4 = 2
auxinput4_inname = 'wrflowinp_d<domain>'
auxinput4_interval = 360, 360, 360,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
debug_level = 0
output_diagnostics = 0
auxhist3_outname = 'wrfxtrm_d01_2015-02-10_00:00:00',
io_form_auxhist3 = 2,
auxhist3_interval = 1440, 1440, 1440,
frames_per_auxhist3 = 1000, 1000, 1000,
force_use_old_data = .true.,
/
&domains
time_step = 50,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
e_we = 190, 166, 106,
e_sn = 166, 148, 106,
e_vert = 40, 40, 40,
p_top_requested = 5000,
num_metgrid_levels = 38,
num_metgrid_soil_levels = 4,
dx = 9000, 3000, 1000,
dy = 9000, 3000, 1000,
grid_id = 1, 2, 3,
parent_id = 0, 1, 2,
i_parent_start = 1, 68, 66,
j_parent_start = 1, 59, 57,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 0
/
&physics
mp_physics = 2, 2, 2,
ra_lw_physics = 4, 4, 4,
ra_sw_physics = 4, 4, 4,
radt = 5, 5, 5,
swint_opt = 1,
sf_sfclay_physics = 2, 2, 2,
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 8, 8, 8,
bldt = 0, 0, 0,
cu_physics = 5, 0, 0,
cudt = 5, 5, 5,
cu_diag = 1, 0, 0,
cu_rad_feedback = .true.,.false.,.false.
isfflx = 1,
ifsnow = 0,
icloud = 1,
surface_input_source = 1,
num_land_cat = 24,
num_soil_layers = 4,
sst_update = 1,
tmn_update = 1,
lagday = 150,
sst_skin = 1,
bucket_mm = 10000.0,
bucket_J = -1,
sf_urban_physics = 1, 1, 1,
/
&fdda
/
&dynamics
w_damping = 0,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12,
base_temp = 290.
damp_opt = 0,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
non_hydrostatic = .true., .true., .true.,
moist_adv_opt = 1, 1, 1,
scalar_adv_opt = 1, 1, 1,
chem_adv_opt = 1, 1, 1,
tke_adv_opt = 1, 1, 1,
gwd_opt = 1,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,.false.,
spec_exp = 0.33,
nested = .false., .true., .true.,
/
&grib2
/
&chem
chem_opt = 2, 2, 2,
chemdt = 8, 4, 1,
input_chem_inname = 'wrfchemi_d<domain>_<date>',
emiss_opt = 3, 3, 3,
emiss_inpt_opt = 102,102, 102,
phot_opt = 1, 1, 1,
photdt = 12, 4, 1,
chem_conv_tr = 1, 0, 0,
io_style_emissions = 2,
chem_in_opt = 1,1,1
have_bcs_chem = .true.,.false.,.false.,
have_bcs_upper = .true.,.true.,.true.
fixed_ubc_inname = 'ubvals_b40.20th.track1_1996-2005.nc',
trop_lev_inname = 'clim_p_trop.nc',
bio_emiss_opt = 3,3,3
bioemdt = 12, 4, 1,
ne_area = 296,
biomass_burn_opt = 0, 0, 0,
plumerisefire_frq = 30, 30, 30,
scale_fire_emiss = .false.,.false.,.false.
kemit = 9,
wetscav_onoff = -1, -1, -1,
cldchem_onoff = 0, 0, 0,
gas_drydep_opt = 1, 1, 1,
aer_drydep_opt = 1, 1, 1,
dust_opt = 0,
seas_opt = 0,
dmsemis_opt = 0,
aer_op_opt = 1, 1, 1,
opt_pars_out = 1,
diagnostic_dep = 1, 1, 1,
aer_ra_feedback = 1, 1, 1,
lnox_opt = 1, 2, 2,
N_IC = 125., 500., 500.,
N_CG = 125., 500., 500.,
/
&control
domains = 3,
start_lai_mnth = 1,
end_lai_mnth = 12,
wrf_dir = './',
megan_dir = '/net/meop-nas2.priv/volume2/data2/ucci/WRFchem/ICBC_chem/MEGAN'
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
I am running WRF-Chem 4.0.3 on 9km, 3km and 1km domain (two-way nesting) but surface values of aerosol concentrations (partly also gases) are in area of inner domains clearly 2-3 times lower in comparison to a run on 9km domain only. It is not problem with emissions but it seems there is a higher vertical mixing in 3km and 1km domains, pollutants are distributed more equally in PBL, while in 9km only run, they are concentrated more near surface. Any idea, tips to solute this discrepancy?
Thanks, Jan
My namelist:
&time_control
run_days = 15,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2015, 2015, 2015,
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 = 2015, 2015, 2015,
end_month = 02, 02, 02,
end_day = 25, 25, 25,
end_hour = 00, 00, 00,
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 = 20000, 20000, 20000,
restart = .false.,
restart_interval = 21600,
auxinput12_inname = 'wrf_chem_input',
io_form_auxinput5 = 2,
auxinput5_inname = 'wrfchemi_d<domain>_<date>',
io_form_auxinput5 = 2,
auxinput5_interval_m = 60,
frames_per_auxinput5 = 24,
auxinput6_inname = 'wrfbiochemi_d<domain>',
auxinput8_interval_m = 1440,
io_form_auxinput2 = 2,
io_form_auxinput6 = 2,
io_form_auxinput8 = 0,
io_form_auxinput4 = 2
auxinput4_inname = 'wrflowinp_d<domain>'
auxinput4_interval = 360, 360, 360,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
debug_level = 0
output_diagnostics = 0
auxhist3_outname = 'wrfxtrm_d01_2015-02-10_00:00:00',
io_form_auxhist3 = 2,
auxhist3_interval = 1440, 1440, 1440,
frames_per_auxhist3 = 1000, 1000, 1000,
force_use_old_data = .true.,
/
&domains
time_step = 50,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
e_we = 190, 166, 106,
e_sn = 166, 148, 106,
e_vert = 40, 40, 40,
p_top_requested = 5000,
num_metgrid_levels = 38,
num_metgrid_soil_levels = 4,
dx = 9000, 3000, 1000,
dy = 9000, 3000, 1000,
grid_id = 1, 2, 3,
parent_id = 0, 1, 2,
i_parent_start = 1, 68, 66,
j_parent_start = 1, 59, 57,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 0
/
&physics
mp_physics = 2, 2, 2,
ra_lw_physics = 4, 4, 4,
ra_sw_physics = 4, 4, 4,
radt = 5, 5, 5,
swint_opt = 1,
sf_sfclay_physics = 2, 2, 2,
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 8, 8, 8,
bldt = 0, 0, 0,
cu_physics = 5, 0, 0,
cudt = 5, 5, 5,
cu_diag = 1, 0, 0,
cu_rad_feedback = .true.,.false.,.false.
isfflx = 1,
ifsnow = 0,
icloud = 1,
surface_input_source = 1,
num_land_cat = 24,
num_soil_layers = 4,
sst_update = 1,
tmn_update = 1,
lagday = 150,
sst_skin = 1,
bucket_mm = 10000.0,
bucket_J = -1,
sf_urban_physics = 1, 1, 1,
/
&fdda
/
&dynamics
w_damping = 0,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12,
base_temp = 290.
damp_opt = 0,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
non_hydrostatic = .true., .true., .true.,
moist_adv_opt = 1, 1, 1,
scalar_adv_opt = 1, 1, 1,
chem_adv_opt = 1, 1, 1,
tke_adv_opt = 1, 1, 1,
gwd_opt = 1,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,.false.,
spec_exp = 0.33,
nested = .false., .true., .true.,
/
&grib2
/
&chem
chem_opt = 2, 2, 2,
chemdt = 8, 4, 1,
input_chem_inname = 'wrfchemi_d<domain>_<date>',
emiss_opt = 3, 3, 3,
emiss_inpt_opt = 102,102, 102,
phot_opt = 1, 1, 1,
photdt = 12, 4, 1,
chem_conv_tr = 1, 0, 0,
io_style_emissions = 2,
chem_in_opt = 1,1,1
have_bcs_chem = .true.,.false.,.false.,
have_bcs_upper = .true.,.true.,.true.
fixed_ubc_inname = 'ubvals_b40.20th.track1_1996-2005.nc',
trop_lev_inname = 'clim_p_trop.nc',
bio_emiss_opt = 3,3,3
bioemdt = 12, 4, 1,
ne_area = 296,
biomass_burn_opt = 0, 0, 0,
plumerisefire_frq = 30, 30, 30,
scale_fire_emiss = .false.,.false.,.false.
kemit = 9,
wetscav_onoff = -1, -1, -1,
cldchem_onoff = 0, 0, 0,
gas_drydep_opt = 1, 1, 1,
aer_drydep_opt = 1, 1, 1,
dust_opt = 0,
seas_opt = 0,
dmsemis_opt = 0,
aer_op_opt = 1, 1, 1,
opt_pars_out = 1,
diagnostic_dep = 1, 1, 1,
aer_ra_feedback = 1, 1, 1,
lnox_opt = 1, 2, 2,
N_IC = 125., 500., 500.,
N_CG = 125., 500., 500.,
/
&control
domains = 3,
start_lai_mnth = 1,
end_lai_mnth = 12,
wrf_dir = './',
megan_dir = '/net/meop-nas2.priv/volume2/data2/ucci/WRFchem/ICBC_chem/MEGAN'
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/