Hi everyone, I’m facing a tricky issue with my WRF-Chem simulation and would really appreciate any help or insights from those who have experience with anthro_emis tool and EDGAR inventories.
I’m using WRF-Chem v4.0 and the anthro_emis preprocessor (from UCAR ACOM: WRF-Chem Tools for the Community | Atmospheric Chemistry Observations & Modeling) to process the EDGAR-HTAP 2010 emission inventory. I’ve generated hourly emission files (e.g., wrfchemi_d01_2020-02-08_20:00:00) and the model runs smoothly with 24 cores—no errors reported at all.
However, the simulated concentrations of PM2.5, NH3, and other species are 3–4 orders of magnitude lower than observed values, almost approaching zero. After some troubleshooting, I suspect the problem lies in two main areas:
1. Unit conversion or species mapping errors in anthro_emis: Due to version mismatches between WRF-Chem, EDGAR inventory, and the chemical mechanism (MOZART-MOSAIC), there might be an order-of-magnitude discrepancy during the inventory processing.
2. Incorrect namelist.input settings: The model may not be reading the offline emission files properly, and instead only uses extremely low initial background concentrations, which decay to near-zero over integration.
From past forum posts, I’ve learned that parameters like auxinput5_interval_m, chem_opt, io_form_auxinput5, and emiss_inpt_opt can significantly affect emission input. I’ve done some test adjustments:
- When I changed auxinput5_interval_m from 60 to 100000, the simulated concentrations became more reasonable—but I know this is not a physically correct setting.
- Previously, the model would crash at the first hour of the second day because it couldn’t recognize the next day’s emission files. After adding force_use_old_data=.true., the model runs without crashing, but the concentrations remain extremely low.
I’ve attached my MOZART-MOSAIC.inp file (used in anthro_emis) and my namelist.input (which I’ve heavily modified during troubleshooting, so it’s a bit of a mess now).
If anyone has experience with this combination (anthro_emis + EDGAR-HTAP 2010 + WRF-Chem v4.0 + MOZART-MOSAIC), I would be incredibly grateful for your help with:
- Identifying potential errors in my namelist.input or MOZART-MOSAIC.inp settings;
- Sharing a working set of anthro_emis MOZART-MOSAIC.inp and WRF-Chem namelist.input files for reference.
Thank you so much in advance for your time and assistance!
Best regards,
Donglei Fu
Email: 2025092@ynau.edu.cn
namelist.input:
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 = 100000, 100000, 100000,
! 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,
force_use_old_data=.true.
&chem
kemit = 1,
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,
phot_opt = 1, 1, 1,
gas_drydep_opt = 0, 0, 0,
aer_drydep_opt = 0, 0, 0,
bio_emiss_opt = 0, 0, 0,
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 = 0, 0, 0,
have_bcs_chem = .false., .false., .false.,
MOZART-MOSAIC.inp:
&CONTROL
anthro_dir = '/home/lulu/models/Emis/EDGAR-HTAP/MOZART_MOSAIC'
wrf_dir = '/home/lulu/models/WRF3.9/WRFCHEM/WRF/test/em_real'
src_file_prefix = 'EDGAR_HTAP_emi_'
src_file_suffix = '_2010.0.1x0.1.nc'
src_names = 'CO(28)','NOx(30)','SO2(64)','NH3(17)','BC(12)','OC(12)','PM2.5(1)','PM10(1)',
'BIGALK(72)','BIGENE(56)','C2H4(28)','C2H5OH(46)','C2H6(30)','CH2O(30)',
'CH3CHO(44)','CH3COCH3(58)','CH3OH(32)','MEK(72)','TOLUENE(92)','C3H6(42)','C3H8(44)',
'BENZENE(78)','XYLENE(106)'
sub_categories = 'emis_tot'
cat_var_prefix = ' '
serial_output = .true.
start_output_time = '2020-01-01_00:00:00'
stop_output_time = '2020-01-10_00:00:00'
output_interval = 3600
data_yrs_offset = 10,
domains =2
emissions_zdim_stag = 1
emis_map = 'CO->CO','NO->NOx','SO2->SO2','NH3->NH3','BIGALK->BIGALK','BIGENE->BIGENE',
'C2H4->C2H4','C2H5OH->C2H5OH','C2H6->C2H6','CH2O->CH2O','CH3CHO->CH3CHO',
'CH3COCH3->CH3COCH3','CH3OH->CH3OH','MEK->MEK','TOLUENE->TOLUENE',
'C3H6->C3H6','C3H8->C3H8','NO2->0.0*NOx','ISOP->0.0*CO','C10H16->0.0*CO',
'SULF->0.0*SO2','C2H2->0.00561790*CO','BENZENE->BENZENE','XYLENE->XYLENE',
'GLY->0.0*CO','MACR->0.0*CO','MGLY->0.0*CO','MVK->0.0*CO',
'HCOOH->0.0*CO','HONO->0.0*CO',
'ECI(a)->0.2*BC','ECJ(a)->0.8*BC','ORGI(a)->0.2*OC','ORGJ(a)->0.8*OC','PM25I(a)->0.2*PM2.5',
'PM25J(a)->0.8*PM2.5','PM_10(a)->PM10','SO4I(a)->0.0*PM10','SO4J(a)->0.0*PM10','NO3I(a)->0.0*PM10',
'NO3J(a)->0.0*PM10','NH4I(a)->0.0*PM10','NH4J(a)->0.0*PM10','NAI(a)->0.0*PM10','NAJ(a)->0.0*PM10',
'CLI(a)->0.0*PM10','CLJ(a)->0.0*PM10','CO_A->CO','CO_BB->0.0*CO','ORGI_A(a)->0.0*PM10',
'ORGI_BB(a)->0.0*PM10','ORGJ_A(a)->0.0*PM10','ORGJ_BB(a)->0.0*PM10','VOCA->0.0*PM10','VOCBB->0.0*PM10'
I’m using WRF-Chem v4.0 and the anthro_emis preprocessor (from UCAR ACOM: WRF-Chem Tools for the Community | Atmospheric Chemistry Observations & Modeling) to process the EDGAR-HTAP 2010 emission inventory. I’ve generated hourly emission files (e.g., wrfchemi_d01_2020-02-08_20:00:00) and the model runs smoothly with 24 cores—no errors reported at all.
However, the simulated concentrations of PM2.5, NH3, and other species are 3–4 orders of magnitude lower than observed values, almost approaching zero. After some troubleshooting, I suspect the problem lies in two main areas:
1. Unit conversion or species mapping errors in anthro_emis: Due to version mismatches between WRF-Chem, EDGAR inventory, and the chemical mechanism (MOZART-MOSAIC), there might be an order-of-magnitude discrepancy during the inventory processing.
2. Incorrect namelist.input settings: The model may not be reading the offline emission files properly, and instead only uses extremely low initial background concentrations, which decay to near-zero over integration.
From past forum posts, I’ve learned that parameters like auxinput5_interval_m, chem_opt, io_form_auxinput5, and emiss_inpt_opt can significantly affect emission input. I’ve done some test adjustments:
- When I changed auxinput5_interval_m from 60 to 100000, the simulated concentrations became more reasonable—but I know this is not a physically correct setting.
- Previously, the model would crash at the first hour of the second day because it couldn’t recognize the next day’s emission files. After adding force_use_old_data=.true., the model runs without crashing, but the concentrations remain extremely low.
I’ve attached my MOZART-MOSAIC.inp file (used in anthro_emis) and my namelist.input (which I’ve heavily modified during troubleshooting, so it’s a bit of a mess now).
If anyone has experience with this combination (anthro_emis + EDGAR-HTAP 2010 + WRF-Chem v4.0 + MOZART-MOSAIC), I would be incredibly grateful for your help with:
- Identifying potential errors in my namelist.input or MOZART-MOSAIC.inp settings;
- Sharing a working set of anthro_emis MOZART-MOSAIC.inp and WRF-Chem namelist.input files for reference.
Thank you so much in advance for your time and assistance!
Best regards,
Donglei Fu
Email: 2025092@ynau.edu.cn
namelist.input:
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 = 100000, 100000, 100000,
! 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,
force_use_old_data=.true.
&chem
kemit = 1,
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,
phot_opt = 1, 1, 1,
gas_drydep_opt = 0, 0, 0,
aer_drydep_opt = 0, 0, 0,
bio_emiss_opt = 0, 0, 0,
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 = 0, 0, 0,
have_bcs_chem = .false., .false., .false.,
MOZART-MOSAIC.inp:
&CONTROL
anthro_dir = '/home/lulu/models/Emis/EDGAR-HTAP/MOZART_MOSAIC'
wrf_dir = '/home/lulu/models/WRF3.9/WRFCHEM/WRF/test/em_real'
src_file_prefix = 'EDGAR_HTAP_emi_'
src_file_suffix = '_2010.0.1x0.1.nc'
src_names = 'CO(28)','NOx(30)','SO2(64)','NH3(17)','BC(12)','OC(12)','PM2.5(1)','PM10(1)',
'BIGALK(72)','BIGENE(56)','C2H4(28)','C2H5OH(46)','C2H6(30)','CH2O(30)',
'CH3CHO(44)','CH3COCH3(58)','CH3OH(32)','MEK(72)','TOLUENE(92)','C3H6(42)','C3H8(44)',
'BENZENE(78)','XYLENE(106)'
sub_categories = 'emis_tot'
cat_var_prefix = ' '
serial_output = .true.
start_output_time = '2020-01-01_00:00:00'
stop_output_time = '2020-01-10_00:00:00'
output_interval = 3600
data_yrs_offset = 10,
domains =2
emissions_zdim_stag = 1
emis_map = 'CO->CO','NO->NOx','SO2->SO2','NH3->NH3','BIGALK->BIGALK','BIGENE->BIGENE',
'C2H4->C2H4','C2H5OH->C2H5OH','C2H6->C2H6','CH2O->CH2O','CH3CHO->CH3CHO',
'CH3COCH3->CH3COCH3','CH3OH->CH3OH','MEK->MEK','TOLUENE->TOLUENE',
'C3H6->C3H6','C3H8->C3H8','NO2->0.0*NOx','ISOP->0.0*CO','C10H16->0.0*CO',
'SULF->0.0*SO2','C2H2->0.00561790*CO','BENZENE->BENZENE','XYLENE->XYLENE',
'GLY->0.0*CO','MACR->0.0*CO','MGLY->0.0*CO','MVK->0.0*CO',
'HCOOH->0.0*CO','HONO->0.0*CO',
'ECI(a)->0.2*BC','ECJ(a)->0.8*BC','ORGI(a)->0.2*OC','ORGJ(a)->0.8*OC','PM25I(a)->0.2*PM2.5',
'PM25J(a)->0.8*PM2.5','PM_10(a)->PM10','SO4I(a)->0.0*PM10','SO4J(a)->0.0*PM10','NO3I(a)->0.0*PM10',
'NO3J(a)->0.0*PM10','NH4I(a)->0.0*PM10','NH4J(a)->0.0*PM10','NAI(a)->0.0*PM10','NAJ(a)->0.0*PM10',
'CLI(a)->0.0*PM10','CLJ(a)->0.0*PM10','CO_A->CO','CO_BB->0.0*CO','ORGI_A(a)->0.0*PM10',
'ORGI_BB(a)->0.0*PM10','ORGJ_A(a)->0.0*PM10','ORGJ_BB(a)->0.0*PM10','VOCA->0.0*PM10','VOCBB->0.0*PM10'