Hi,
I am working with WRF-Chem v4.6.1 using the MOZMEM chemical mechanism and have set emiss_opt ==10 to read point source emissions.
I have a point source emission inventory that includes both CH4 (methane) and benzene. The model correctly reads e_benzene emission and generate benzene output. However mozmem does not correctly read e_ch4 emission input and shows zero values in the wrfout files.
To address this, I modified the following:
Added e_ch4 to the Registry/registry.chem file in the emiss_opt==10 list.
Updated chem/module_emissions_anthropogenics.F after the p_apin line (line 138) to include handling for e_ch4.
if( p_ch4 >= param_first_scalar ) then
chem(its:ite,k,j,p_ch4) = chem(its:ite,k,j,p_ch4) + emis_ant(its:ite,k,j,p_e_ch4)*conv_rho(its:ite)
endif
Despite these changes, methane still does not appear in the output, while benzene is processed and output correctly. Are there additional modifications required to fully register e_ch4 in MOZMEM (emiss_opt == 10)?
Is there a specific location in module_emissions_anthropogenics.F or elsewhere where e_ch4 must be mapped for MOZMEM to ensure e_ch4 is included in the emissions correctly?
The namelist that I used for the chemistry is included:
kemit = 1,
chem_opt = 201, 201, 201, 201, 201,
!chem_opt = 0, 0, 0, 0, 0,
bioemdt = 30, 30, 30, 30, 30,
photdt = 30, 30, 30, 30, 30,
chemdt = 10, 10, 10, 10, 10,
emi_inname = "methane_d<domain>.nc"
io_style_emissions = 2,
mozart_ph_diag = 1
emiss_inpt_opt = 102, 102, 102, 102, 102,
emiss_opt = 10, 10, 10, 10, 10,
chem_in_opt = 0, 0, 0, 0, 0,
phot_opt = 1, 1, 1, 1, 1,
aircraft_emiss_opt = 0, 0, 0, 0, 0,
gas_drydep_opt = 1, 1, 1, 1, 1,
aer_drydep_opt = 1, 1, 1, 1, 1,
bio_emiss_opt = 3, 3, 3, 3, 3,
dust_opt = 1,
dmsemis_opt = 1,
seas_opt = 1,
gas_bc_opt = 1, 1, 1, 1, 1,
gas_ic_opt = 1, 1, 1, 1, 1,
aer_bc_opt = 1, 1, 1, 1, 1,
aer_ic_opt = 1, 1, 1, 1, 1,
gaschem_onoff = 1, 1, 1, 1, 1,
aerchem_onoff = 1, 1, 1, 1, 1,
wetscav_onoff = 0, 0, 0, 0, 0,
cldchem_onoff = 0, 0, 0, 0, 0,
vertmix_onoff = 1, 1, 1, 1, 1,
chem_conv_tr = 0, 0, 0, 0, 0,
biomass_burn_opt = 4, 4, 4, 4, 4,
plumerisefire_frq = 30, 30, 30, 30, 30,
scale_fire_emiss = .true
aer_ra_feedback = 1, 1, 1, 1, 1,
aer_op_opt = 1, 1, 1, 1, 1,
opt_pars_out = 1,
ne_area = 210,
!cld_ra_feedback = 0, 0, 0, 0, 0,
have_bcs_chem = .false.,.false., .false., .false.,.false.,
Thank you,
Sunny
I am working with WRF-Chem v4.6.1 using the MOZMEM chemical mechanism and have set emiss_opt ==10 to read point source emissions.
I have a point source emission inventory that includes both CH4 (methane) and benzene. The model correctly reads e_benzene emission and generate benzene output. However mozmem does not correctly read e_ch4 emission input and shows zero values in the wrfout files.
To address this, I modified the following:
Added e_ch4 to the Registry/registry.chem file in the emiss_opt==10 list.
Updated chem/module_emissions_anthropogenics.F after the p_apin line (line 138) to include handling for e_ch4.
if( p_ch4 >= param_first_scalar ) then
chem(its:ite,k,j,p_ch4) = chem(its:ite,k,j,p_ch4) + emis_ant(its:ite,k,j,p_e_ch4)*conv_rho(its:ite)
endif
Despite these changes, methane still does not appear in the output, while benzene is processed and output correctly. Are there additional modifications required to fully register e_ch4 in MOZMEM (emiss_opt == 10)?
Is there a specific location in module_emissions_anthropogenics.F or elsewhere where e_ch4 must be mapped for MOZMEM to ensure e_ch4 is included in the emissions correctly?
The namelist that I used for the chemistry is included:
kemit = 1,
chem_opt = 201, 201, 201, 201, 201,
!chem_opt = 0, 0, 0, 0, 0,
bioemdt = 30, 30, 30, 30, 30,
photdt = 30, 30, 30, 30, 30,
chemdt = 10, 10, 10, 10, 10,
emi_inname = "methane_d<domain>.nc"
io_style_emissions = 2,
mozart_ph_diag = 1
emiss_inpt_opt = 102, 102, 102, 102, 102,
emiss_opt = 10, 10, 10, 10, 10,
chem_in_opt = 0, 0, 0, 0, 0,
phot_opt = 1, 1, 1, 1, 1,
aircraft_emiss_opt = 0, 0, 0, 0, 0,
gas_drydep_opt = 1, 1, 1, 1, 1,
aer_drydep_opt = 1, 1, 1, 1, 1,
bio_emiss_opt = 3, 3, 3, 3, 3,
dust_opt = 1,
dmsemis_opt = 1,
seas_opt = 1,
gas_bc_opt = 1, 1, 1, 1, 1,
gas_ic_opt = 1, 1, 1, 1, 1,
aer_bc_opt = 1, 1, 1, 1, 1,
aer_ic_opt = 1, 1, 1, 1, 1,
gaschem_onoff = 1, 1, 1, 1, 1,
aerchem_onoff = 1, 1, 1, 1, 1,
wetscav_onoff = 0, 0, 0, 0, 0,
cldchem_onoff = 0, 0, 0, 0, 0,
vertmix_onoff = 1, 1, 1, 1, 1,
chem_conv_tr = 0, 0, 0, 0, 0,
biomass_burn_opt = 4, 4, 4, 4, 4,
plumerisefire_frq = 30, 30, 30, 30, 30,
scale_fire_emiss = .true
aer_ra_feedback = 1, 1, 1, 1, 1,
aer_op_opt = 1, 1, 1, 1, 1,
opt_pars_out = 1,
ne_area = 210,
!cld_ra_feedback = 0, 0, 0, 0, 0,
have_bcs_chem = .false.,.false., .false., .false.,.false.,
Thank you,
Sunny