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

Sequence of input processing for bioemiss + BC + IC

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.

HJORQUER

New member
Hello everybody

I am setting up a nested (4 domain) wrf-chem 4.2.2 run.

I use the following sequence of processes for getting the additional inputs:

a) create wrfbdy_d01, wrfinp_d01-d04 and wrflowinp_d01-d04 using real.exe

b) use bioemiss to generate wrfbiochemi_d01-d04

c) use real to update wrfbdy_d01 and wrfinp_d01-d04

d) use mozbc to generate BC + ICs for D01-D04

e) launch wrf.exe to generate the results

The user's guide is not totally clear (to me) regarding the namelist setting for this modeling case. The following is an excerpt from the namelist.input file:

&time_control
run_days = 0,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2021, 2021, 2021, 2021,
start_month = 03, 03, 03, 03,
start_day = 03, 03, 03, 03,
start_hour = 00, 00, 00, 00,
start_minute = 00, 00, 00, 00,
start_second = 00, 00, 00, 00,
end_year = 2021, 2021, 2021, 2021,
end_month = 03, 03, 03, 03,
end_day = 13, 13, 13, 13,
end_hour = 00, 00, 00, 00,
end_minute = 00, 00, 00, 00,
end_second = 00, 00, 00, 00,
interval_seconds = 21600
input_from_file = .true.,.true.,.true., .true.,
fine_input_stream = 0, 0, 0, 0, 2
history_interval = 180, 60, 60, 60,
frames_per_outfile = 8, 24, 24, 24,
restart = .false.,
restart_interval = 14400,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
io_form_auxinput2 = 2,
io_form_auxinput4 = 0,
io_form_auxinput5 = 2,
io_form_auxinput6 = 2,
io_form_auxinput7 = 0,
io_form_auxinput12 = 0,
auxinput4_interval_m = 360, 360, 360, 360,
auxinput5_interval_m = 14400, 14400, 14400, 14400,
auxinput7_interval_m = 60, 60, 60, 60,
auxinput4_inname = "wrflowinp_d<domain>",
auxinput4_interval = 360, 360, 360, 360,
auxinput5_inname = "wrfchemi_d<domain>_<date>",
auxinput6_inname = "wrfbiochemi_d<domain>",
frames_per_auxinput5 = 1, 1, 1, 1,
frames_per_auxinput6 = 1, 1, 1, 1,
debug_level = 0,
force_use_old_data = .true.,
/

&chem
kemit = 15,
chem_opt = 2, 2, 2, 2,
bioemdt = 30, 30, 30, 30,
photdt = 30, 30, 30, 30,
chemdt = 0, 0, 0, 0,
io_style_emissions = 2,
emiss_inpt_opt = 1, 1, 1, 1,
emiss_opt = 3, 3, 3, 3,
emi_inname = "wrfchemi_d<domain>_<date>",
phot_opt = 1, 1, 1, 1,
gas_drydep_opt = 1, 1, 1, 1,
aer_drydep_opt = 1, 1, 1, 1,
bio_emiss_opt = 3, 3, 3, 3,
ne_area = 500,
opt_pars_out = 1,
dust_opt = 0,
dmsemis_opt = 0,
seas_opt = 2,
gas_bc_opt = 1, 1, 1, 1,
gas_ic_opt = 1, 1, 1, 1,
aer_bc_opt = 1, 1, 1, 1,
aer_ic_opt = 1, 1, 1, 1,
gaschem_onoff = 1, 1, 1, 1,
aerchem_onoff = 1, 1, 1, 1,
wetscav_onoff = 0, 0, 0, 0,
cldchem_onoff = 0, 0, 0, 0,
vertmix_onoff = 1, 1, 1, 1,
chem_conv_tr = 0, 0, 0, 0,
biomass_burn_opt = 0, 0, 0, 0,
aer_ra_feedback = 0, 0, 0, 0,
chem_in_opt = 1, 1, 1, 1,
input_chem_inname = "wrfinput_d<domain>",
have_bcs_chem = .true., .false., .false., .false.,
/

Am I setting this up fine? Or am I making a mistake?

Regards
 
Top