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

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

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.

yuzuoing

New member
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.1 KB · Views: 41
Top