How can an integrated emission file be read by WRF-CHEM?

Post Reply
yuzuoing
Posts: 10
Joined: Fri Dec 14, 2018 1:27 am

How can an integrated emission file be read by WRF-CHEM?

Post by yuzuoing » Sat Dec 26, 2020 7:48 am

Hi,

I am using the WRF-CHEM to do runnings. Now I can just create three single emission files contained all simulated periods for the three domains instead of separate hourly data files (e.g. wrfchemi_d01_2010-05-01_08:00:00, wrfchemi_d01_2010-05-01_13:00:00, hundreds of files for a long-term simulation). My namelist.input is provided here.

Two puzzles are concerned:

1) the model runs smoothly without any errors happened while using the single emission file, I am not sure whether the model reads and recognizes the three integrated emission files.

2) interestingly, the model can obtain 10 days simulated result by only employing a 24-h emssion file (an integerated file containing 24-h), is it related to the namelist.input configurations?

Any available responses are grateful!


&time_control
run_days = 9,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2019, 2019, 2019,
start_month = 12, 12, 12,
start_day = 01, 01, 01,
start_hour = 00, 00, 00,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2019, 2019, 2019,
end_month = 12, 12, 12,
end_day = 10, 10, 10,
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 = 360, 360, 180,
frames_per_outfile = 5000, 5000, 5000,
restart = .false.,
write_hist_at_0h_rst = .false.,
restart_interval = 1440,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
io_form_auxinput2 = 2,
auxinput1_inname = "met_em.d<domain>.<date>",
auxinput5_inname = "wrfchemi_d<domain>_<date>",
auxinput6_inname = 'wrfbiochemi_d<domain>',
!auxinput7_inname = 'wrffirechemi_d<domain>',
!auxinput8_inname = 'wrfchemi_gocart_bg_d<domain>',
auxinput12_inname = 'wrf_chem_input',
auxinput5_interval_m = 60, 60, 60,
! auxinput7_interval = 21600, 21600, 21600,
! auxinput8_interval = 21600, 21600, 21600,
io_form_auxinput5 = 2,
io_form_auxinput6 = 2,
io_form_auxinput7 = 0,
io_form_auxinput8 = 0,
io_form_auxinput12 = 0,
io_form_auxinput13 = 0,
io_form_auxinput14 = 0,
io_form_auxinput15 = 0,
!frames_per_auxinput5 = 1,

debug_level = 0,

/

&domains
time_step = 180,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
e_we = 55, 67, 46,
e_sn = 55, 67, 46,
e_vert = 30, 30, 30,
p_top_requested = 5000,
eta_levels = 1.0000,0.9970,0.9930,0.9880,0.9800,0.9700,
0.9600,0.9500,0.9400,0.9300,0.9100,0.8900,
0.8700,0.8400,0.8000,0.7600,0.7200,0.6800,
0.6400,0.6000,0.5500,0.5000,0.4500,0.4000,
0.3250,0.2500,0.1750,0.1000,0.0400,0.0000
num_metgrid_levels = 34,
num_metgrid_soil_levels = 4,
dx = 27000, 9000, 3000,
dy = 27000, 9000, 3000,
grid_id = 1, 2, 3,
parent_id = 1, 1, 2,
i_parent_start = 1, 17, 26,
j_parent_start = 1, 17, 27,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 0
smooth_cg_topo = .true.
/

&physics
mp_physics = 2, 2, 2,
progn = 0, 0, 0,
naer = 1e9
ra_lw_physics = 1, 1, 1,
ra_sw_physics = 1, 1, 1,
radt = 30, 30, 30,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 3, 3, 3,
bl_pbl_physics = 1, 1, 1,
bldt = 0, 0, 0,
usemonalb = .true.,
rdlai2d = .true.,
mosaic_lu = 1
mosaic_soil = 1
cu_physics = 0, 0, 0,
cudt = 0, 0, 0,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 3,
num_soil_layers = 4,
num_land_cat = 24,
sf_urban_physics = 0, 0, 0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
cu_rad_feedback = .false.,
cu_diag = 0
/

&fdda
/

&dynamics
w_damping = 1,
diff_opt = 1, 1, 1,
km_opt = 4, 4, 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 = 2, 2, 2,
scalar_adv_opt = 2, 2, 2,
chem_adv_opt = 1, 1, 1,
/

&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true.,.false.,.false.,
nested = .false., .true., .true.,
/

&grib2
/

&chem
kemit = 6,
chem_opt = 2, 2, 2,
photdt = 30, 30, 30,
chemdt = 2., 0.66666, 0.22222,
!frames_per_emissfile = 1, 1, 1,
io_style_emissions = 2,
emiss_inpt_opt = 1, 1, 1,
emiss_opt = 3, 3, 3,
chem_in_opt = 0, 0, 0,
ne_area = 400,
phot_opt = 1, 1, 1,
gas_drydep_opt = 1, 1, 1,
aer_drydep_opt = 1, 1, 1,
bio_emiss_opt = 2, 2, 2,
dust_opt = 0,
dmsemis_opt = 0,
seas_opt = 0,
gas_bc_opt = 1, 1, 1,
gas_ic_opt = 1, 1, 1,
aer_bc_opt = 1, 1, 1,
aer_ic_opt = 1, 1, 1,
gaschem_onoff = 1, 1, 1,
aerchem_onoff = 1, 1, 1,
wetscav_onoff = 0, 0, 0,
cldchem_onoff = 0, 0, 0,
vertmix_onoff = 1, 1, 1,
chem_conv_tr = 0, 0, 0,
biomass_burn_opt = 0, 0, 0,
plumerisefire_frq = 30, 30, 30,
aer_ra_feedback = 1, 1 1,
have_bcs_chem = .false., .true., .true.,
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
Attachments
namelist.input
(8.13 KiB) Downloaded 8 times

Post Reply

Return to “Emissions”