hiep.duc@environment.nsw.
New member
I want to study the effect of hydrogen emission from motor vehicle fleet using H2 fuel. I add H2 emission into the inventory of vehicle emission consisting of NOx, NH3, Aldehyde, ethene,... in the wrfchemi_d0<n> files. The E_H2 as a new emission species in addition to other species emission such as E_CO, E_NO,.... I want to use the MOZART+GOCART mechanism as it has hydrogen as part of the chemical reaction chain. So I use chem_opt = 112 and emiss_opt = 8.
To ingest E_H2 emission into WRF-Chem,
(1) I add to the registry.chem as follows
state real e_h2 i+jf emis_ant 1 Z i5r "E_H2" "H2 EMISSIONS" "mol km^-2 hr^-1"
And in package mozcem section of the registry.chem, I add e_h2 corresponding to emiss_opt = 8
package mozcem emiss_opt==8 - emis_ant:e_h2,e_co,e_no,e_no2,e_bigalk,e_bigene,e_c2h4,e_c2h5oh,e_c2h6,e_c3h6,e_c3h8,e_ch2o,e_ch3cho,e_ch3coch3,e_ch3oh,e_mek,e_so2,e_toluene,e_nh3,e_isop,e_c10h16,e_pm_10,e_pm_25,e_bc,e_oc,e_sulf
(2) Then In chem/module_emissions_anthropogenics.F, I add the following
if( config_flags%chem_opt == MOZCART_KPP .or. config_flags%chem_opt == T1_MOZCART_KPP ) then
if( p_h2 >= param_first_scalar ) then
if (k == kts .and. j == jts) then
write(0,*) 'p_h2=', p_h2
write(0,*) 'p_e_h2=', p_e_h2
write(0,*) 'MAX E_CO = ', maxval(emis_ant
,:,:,p_e_co))
write(0,*) 'MAX E_H2 = ', maxval(emis_ant,:,:,p_e_h2))
endif
chem(its:ite,k,j,p_h2) = chem(its:ite,k,j,p_h2) + emis_ant(its:ite,k,j,p_e_h2) * conv_rho(its:ite)
endif
I then compile wrf_chem, and run the simulation
But no ingestion of H2 and no concentration of H2 in the wrf output.
Any tip is much appreciated.
Best
Hiep Nguyen
To ingest E_H2 emission into WRF-Chem,
(1) I add to the registry.chem as follows
state real e_h2 i+jf emis_ant 1 Z i5r "E_H2" "H2 EMISSIONS" "mol km^-2 hr^-1"
And in package mozcem section of the registry.chem, I add e_h2 corresponding to emiss_opt = 8
package mozcem emiss_opt==8 - emis_ant:e_h2,e_co,e_no,e_no2,e_bigalk,e_bigene,e_c2h4,e_c2h5oh,e_c2h6,e_c3h6,e_c3h8,e_ch2o,e_ch3cho,e_ch3coch3,e_ch3oh,e_mek,e_so2,e_toluene,e_nh3,e_isop,e_c10h16,e_pm_10,e_pm_25,e_bc,e_oc,e_sulf
(2) Then In chem/module_emissions_anthropogenics.F, I add the following
if( config_flags%chem_opt == MOZCART_KPP .or. config_flags%chem_opt == T1_MOZCART_KPP ) then
if( p_h2 >= param_first_scalar ) then
if (k == kts .and. j == jts) then
write(0,*) 'p_h2=', p_h2
write(0,*) 'p_e_h2=', p_e_h2
write(0,*) 'MAX E_CO = ', maxval(emis_ant
,:,:,p_e_co))write(0,*) 'MAX E_H2 = ', maxval(emis_ant
,:,:,p_e_h2))endif
chem(its:ite,k,j,p_h2) = chem(its:ite,k,j,p_h2) + emis_ant(its:ite,k,j,p_e_h2) * conv_rho(its:ite)
endif
I then compile wrf_chem, and run the simulation
But no ingestion of H2 and no concentration of H2 in the wrf output.
Any tip is much appreciated.
Best
Hiep Nguyen