[Wrf-users] E_PM_10 variable not being read from wrfchemi for emiss_opt=4

[Tabish Ansari]

Hi

I'm running WRF-Chem v3.7.1 with chem_opt=10 (CBMZ-MOSAIC with aq.
reactions). I'm specifying emiss_opt=4 for this chemical mechanism. It
turns out that the model is not reading my PM10 emissions (which actually
stores PMcoarse). Ideally, the PMcoarse primary emissions should go into
oin_a07 and oin_a08 variables (see lines of code below) in the wrfout file
but this is not happening.

Lines of code in module_mosaic_addemis.F:

aem_oin = fr_aem_mosaic_f(n)*(emis_ant(i,k,j,p_e_pm25j)+ebu(i,k,j,p_ebu_pm25))
  &

                + *fr_aem_mosaic_c(n)***emis_ant(i,k,j,p_e_pm_10**) *  &

                + fr_aem_sorgam_i(n)*(emis_ant(i,k,j,p_e_pm25i)+0.25*
ebu(i,k,j,p_ebu_pm25))   &

                + fr_aem_sorgam_j(n)*(emis_ant(i,k,j,p_e_pm25j)+0.75*
ebu(i,k,j,p_ebu_pm25))  &

                + fr_tno_mosaic_f(n)*emis_ant(i,k,j,p_e_oin_25) &

            + fr_tno_mosaic_c(n)*emis_ant(i,k,j,p_e_oin_10)


When I looked in the registry.chem file, it turns out that e_pm_10 variable
isn't listed for emiss_opt=4. Here's how it looked like:

package   ecbmz_mosaic    emiss_opt==4                   -
emis_ant:e_iso,e_so2,e_no,e_co,e_eth,e_hc3,e_hc5,e_hc8,e_
xyl,e_ol2,e_olt,e_oli,e_tol,e_csl,e_hcho,e_ald,e_ket,e_ora2,e_
nh3,e_no2,e_ch3oh,e_c2h5oh,e_pm25i,e_pm25j,e_eci,e_ecj,e_
orgi,e_orgj,e_so4i,e_so4j,e_no3i,e_no3j,e_so4c,e_no3c,e_orgc,e_ecc

So I manually added e_pm_10 in the above list in the end and recompiled the
model, however the problem still persists.

Please let me know if this is a standard bug with WRF-Chem v3.7.1 or also
persists in other versions; and how to resolve it.

Thanks

Tabish


Tabish U Ansari
PhD student, Lancaster Environment Center
Lancaster Univeristy
Bailrigg, Lancaster,
LA1 4YW, United Kingdom
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/wrf-users/attachments/20170124/ccfc14cf/attachment.html