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

[Douglas Lowe]

Hi Tabish,

Which emiss_inpt_opt option are you using? It needs to be set to 101 (emiss_inpt_pnnl_cm) if you want to use e_pm_10.
If you're using 1 (emiss_inpt_default) or 102 (emiss_inpt_pnnl_rs) you'll have to add another line of code around line 269 in module_mosaic_addemiss.F to set fr_aem_mosaic_c to be non-zero (though if you do this you should check to make sure it doesn't lead to double counting in any other of your aerosol emissions!).

cheers,
Doug

________________________________
From: wrfchemuserseu-bounces at lists.science.uu.nl [wrfchemuserseu-bounces at lists.science.uu.nl] on behalf of Tabish Ansari [tabishumaransari at gmail.com]
Sent: 24 January 2017 01:04
To: wrf-users at ucar.edu; wrfchemhelp.gsd at noaa.gov; wrfchemuserseu at lists.science.uu.nl; wrf-chem-discussions at list.woc.noaa.gov
Subject: [Wrfchemuserseu] E_PM_10 variable not being read from wrfchemi for emiss_opt=4

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/20170125/3473a38b/attachment.html