Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

Lower aerosol concentrations in nested domains

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

jan1203

New member
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,
/
 
Top