Daily point emissions being dropped for "grouped" pollutants #131
Description
In the latest version of SMOKE, we are seeing dropped emissions for "grouped" HAPs in sectors with daily point inventories (e.g. fires). By "grouped" HAPs, I mean situations where multiple inventory pollutants are mapped to a single SMOKE pollutant via the INVTABLE, for example:
PAH_880E5 83329 846 Y 1.0 N N tons/yr Acenaphthene
PAH_880E5 86737 883 Y 1.0 N N tons/yr Fluorene
PAH_880E5 90120 105 Y 1.0 N N tons/yr 1-Methylnaphthalene
PAH_880E5 91576 196 Y 1.0 N N tons/yr 2-Methylnaphthalene
(and so on)
In daily point sectors such as ptfire, the emissions from these individual pollutants are not preserved; it appears as through the emissions from one component of the PAH group are overriding the emissions from other components of the group, when SMOKE should instead be summing the emissions together from all pollutants within the PAH group.
This was not an issue in SMOKE 5.0, so a recent code update (within the last couple of years) appears to be the cause.
We can provide a more targeted test package if that is needed, but to replicate, you could use the data from 2021hb package, ptfire-wild sector:
https://gaftp.epa.gov/Air/Emismod/2021/
You would need the "invtable_hapcap" and "gspro_hapcap" files from this package to process the PAHs.
We are seeing that when we run the ptfire-wild sector all the way through SMOKE with SMOKE 5.0 and with SMOKE 5.1, the resulting emissions (psot-Smkmerge) are different for many of the PAH_* pollutants, as well as XYLENES, and that the emissions from SMOKE 5.0 are correct while the emissions from SMOKE 5.1 are too low.