<div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Hi<br><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">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.<br><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Lines of code in module_mosaic_addemis.F:<br><p class="m_7188156326065247304gmail-m_-5848940748083566413gmail-p1"><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1">aem_oin </span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s2">=</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1"> fr_aem_mosaic_f(n)</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s2">*</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1">(emis_ant(i<wbr>,k,j,p_e_pm25j)</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s2">+</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1">ebu(i,k,j,p_eb<wbr>u_pm25)) </span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s3">&</span></p>
<p class="m_7188156326065247304gmail-m_-5848940748083566413gmail-p1"><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s4"> </span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1"> </span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s2">+</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1"> <b>fr_aem_mosaic_c(n)</b></span><b><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s2">*</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1">emis_ant(i,<wbr>k,j,p_e_</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s3">pm_10</span></b><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1"><b>) </b> </span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s3">&</span></p>
<p class="m_7188156326065247304gmail-m_-5848940748083566413gmail-p1"><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s4"> </span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1"> </span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s2">+</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1"> fr_aem_sorgam_i(n)</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s2">*</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1">(emis_ant(i<wbr>,k,j,p_e_pm25i)</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s2">+</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s5">0.25</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s2">*</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1">ebu(i,k,j<wbr>,p_ebu_pm25)) </span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s3">&</span></p>
<p class="m_7188156326065247304gmail-m_-5848940748083566413gmail-p1"><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s4"> </span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1"> </span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s2">+</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1"> fr_aem_sorgam_j(n)</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s2">*</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1">(emis_ant(i<wbr>,k,j,p_e_pm25j)</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s2">+</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s5">0.75</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s2">*</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1">ebu(i,k,j<wbr>,p_ebu_pm25)) </span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s3">&</span></p>
<p class="m_7188156326065247304gmail-m_-5848940748083566413gmail-p1"><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s4"> </span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1"> </span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s2">+</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1"> fr_tno_mosaic_f(n)</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s2">*</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1">emis_ant(i,<wbr>k,j,p_e_oin_25) </span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s3">&</span></p>
<p class="m_7188156326065247304gmail-m_-5848940748083566413gmail-p1"><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1"> </span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s2">+</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1"> fr_tno_mosaic_c(n)</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s2">*</span><span class="m_7188156326065247304gmail-m_-5848940748083566413gmail-s1">emis_ant(i,<wbr>k,j,p_e_oin_10)</span></p><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">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:<br><br>package ecbmz_mosaic emiss_opt==4 - emis_ant:e_iso,e_so2,e_no,e_<wbr>co,e_eth,e_hc3,e_hc5,e_hc8,e_<wbr>xyl,e_ol2,e_olt,e_oli,e_tol,e_<wbr>csl,e_hcho,e_ald,e_ket,e_ora2,e_<wbr>nh3,e_no2,e_ch3oh,e_c2h5oh,e_<wbr>pm25i,e_pm25j,e_eci,e_ecj,e_<wbr>orgi,e_orgj,e_so4i,e_so4j,e_<wbr>no3i,e_no3j,e_so4c,e_no3c,e_<wbr>orgc,e_ecc<br><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">So I manually added e_pm_10 in the above list in the end and recompiled the model, however the problem still persists.<br><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">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.<br><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Thanks<br><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Tabish<br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br><br clear="all"></div><div><div class="m_7188156326065247304gmail_signature"><div dir="ltr"><div><div><font size="1"><span style="font-family:tahoma,sans-serif">Tabish U Ansari<br></span></font></div><font size="1"><span style="font-family:tahoma,sans-serif">PhD student, Lancaster Environment Center<br></span></font></div><font size="1"><span style="font-family:tahoma,sans-serif">Lancaster Univeristy<br> <span>Bailrigg</span>, <span>Lancaster</span>, <br><span>LA1 4YW</span>, <span>United Kingdom</span></span></font><br></div></div></div>
</div>