This README includes instructions on how to adapt the existing 2014v2 emissions modeling 
package for NATA (National Air Toxics Assessment) processing. NATA processing involves 
additional toxic species.

Processing of emissions for NATA requires downloading two additional zip file packages, 
in addition to the entire 2014v2 base package described in 
README_2014_2015_2016_v7.2_platform_package.txt.

For NATA processing, install the following two additional packages in $INSTALL_DIR.
The second package listed is only needed for nonroad sector processing.

ancillary_data/ge_dat_for_NATA_addendum.zip
2014emissions/2014fd_inputs_nonroad_NATA_addendum.zip

Instructions for modifying the scripts in the 2014v2 emissions modeling base package 
for NATA modeling are provided in this README. These instructions concern the 
sector-specific scripts from the base package that are located in 
$INSTALL_DIR/2014fd_cb6_14j/scripts/.

== Inventory Table (INVTABLE) ==

For the onroad and onroad_ca_adj sectors (all scripts):
setenv INVTABLE "$GE_DAT/invtable_MOVES2014_NATA_2014v1_06dec2017_nf_v1.txt"

For the onroad_mex sector:
setenv INVTABLE "$GE_DAT/invtable_for_MOVES_Mexico_10may2017_v2.txt"

For the nonroad sector:
setenv INVTABLE "$GE_DAT/invtable_hapcap_NATA_2016_21nov2017_nf_v10.txt"

For all other sectors:
setenv INVTABLE "$GE_DAT/invtable_hapcap_NATA_2016_04jan2018_nf_v9.txt"

== Speciation profiles (GSPRO / GSPROTMP_A) ==

For the onroad and onroad_ca_adj sectors (all scripts):
setenv GSPRO "$GE_DAT/speciation/gspro_MOVES2014_CB6_NATA_2014v1_06dec2017_nf_v1.txt"

For the onroad_mex sector: 
setenv GSPROTMP_A "$GE_DAT/speciation/gspro_onroad_mex_cmaq_cb6_2014fd_nata_cb6_14j_nf.txt"

For the nonroad sector: 
setenv GSPROTMP_A "$GE_DAT/speciation/gspro_nonroad_cmaq_cb6_2014fd_nata_cb6_14j_14apr2017.txt"

For the onroad_can, othar, othpt, ptfire_othna, and beis sectors, DO NOT CHANGE 
the GSPROTMP_A definition from the base scripts. In other words, these sectors should 
use the same GSPROTMP_A for NATA processing as for standard processing.

For all other sectors:
setenv GSPROTMP_A "$GE_DAT/speciation/gspro_cmaq_cb6_2014fd_nata_cb6_14j_nf.txt"

== Speciation cross-reference (GSREF / GSREFTMP_A) ==

For the onroad and onroad_ca_adj sectors (all scripts), keep the GSREF from the base
scripts. In other words, these sectors should use the same GSREF for NATA processing as 
for standard processing.

For the onroad_mex sector:
setenv GSREFTMP_A "$GE_DAT/speciation/gsref_onroad_mex_cmaq_cb6_2014fd_nata_cb6_14j_17jan2017.txt"

For the othpt sector:
setenv GSREFTMP_A "$GE_DAT/speciation/gsref_othpt_cmaq_cb6_2014fd_nata_cb6_14j_23jan2017.txt"

For the ag and othar sectors: 
setenv GSREFTMP_A "$GE_DAT/speciation/gsref_ag_othar_cmaq_cb6_2014fd_nata_cb6_14j_17jan2017.txt"

For the ptegu sector: 
setenv GSREFTMP_A "$GE_DAT/speciation/gsref_ptegu_cmaq_cb6_2014fd_nata_cb6_14j_nf.txt"

For the cmv_c1c2, cmv_c3, nonroad, and rail sectors:
setenv GSREFTMP_A "$GE_DAT/speciation/gsref_cmv_rail_nonroad_cmaq_cb6_2014fd_nata_cb6_14j_12feb2018.txt"

For all other sectors:
setenv GSREFTMP_A "$GE_DAT/speciation/gsref_cmaq_cb6_2014fd_nata_cb6_14j_nf.txt"

== Other script changes for the onroad and onroad_ca_adj sectors ==

For RPD jobs:
setenv MEPROC "$GE_DAT/onroad/meproc_MOVES2014_RPD_AQ_NATA_06dec2017_nf_v5"
setenv MRCLIST "$GE_DAT/onroad/mrclist_RPD_2014v2_NATA_12jul2017_newpostproc_06dec2017_v0"

For RPH jobs:
setenv MEPROC "$GE_DAT/onroad/meproc_MOVES2014_RPH_AQ_NATA_06dec2017_nf_v4"
setenv MRCLIST "$GE_DAT/onroad/mrclist_RPH_2014v2_NATA_12jul2017_newpostproc_06dec2017_v0"

For RPP jobs:
setenv MEPROC "$GE_DAT/onroad/meproc_MOVES2014_RPP_AQ_NATA_06dec2017_nf_v4"
setenv MRCLIST "$GE_DAT/onroad/mrclist_RPP_2014v2_NATA_12jul2017_newpostproc_06dec2017_v0"

For RPV jobs:
setenv MEPROC "$GE_DAT/onroad/meproc_MOVES2014_RPV_AQ_NATA_06dec2017_nf_v4"
setenv MRCLIST "$GE_DAT/onroad/mrclist_RPV_2014v2_NATA_12jul2017_newpostproc_06dec2017_v0"

NATA processing for SMOKE-MOVES also requires a set of MOVES emission factor tables (EFTABLES)
specific to NATA processing. These NATA-specific EF tables are large and are only furnished upon 
request.

== Volcanic mercury emissions ==

The NATA addendum zip (ge_dat_for_NATA_addendum.zip) includes a file containing ]
gridded volcanic mercury emissions for the 12US2 grid: volcanic_hg/volcanos_12US2.ncf

To include this file in your merged 2-D emissions, make the following script / helper 
file changes.

Add this line to sectorlist_2014fd_29jan2018_v0 (located in 2014fd_cb6_14j/scripts):
"volcanic_hg","none","2014","$VOLCROOT$GRID$EXT","SectBaseYr","N","","Y",""

Add this line to Sector_merge_12US2_2014fd_cb6_14j.csh (located in 2014fd_cb6_14j/scripts/merge):
setenv VOLCROOT "$INSTALL_DIR/volcanic_hg/volcanos_"

== SMKINVEN_FORMULA ==

For NATA processing, you must change the SMKINVEN_FORMULA definition for a few sectors.

For the nonroad sector:
setenv SMKINVEN_FORMULA "EXH__PMC=EXH__PM10-EXH__PM2_5,PMC=PM10-PM2_5,DIESEL_PMC=DIESEL_PM10-DIESEL_PM25"

For the cmv_c1c2, cmv_c3, rail, and ptnonipm sectors:
setenv SMKINVEN_FORMULA "PMC=PM10-PM2_5,DIESEL_PMC=DIESEL_PM10-DIESEL_PM25"

== Additional nonroad inventory ==

To include the extra nonroad inventory from 
2014emissions/2014fd_inputs_nonroad_NATA_addendum.zip
in your processing, which is required for NATA, add this line to 
Monthly_nonroad_12US2_2014fd_cb6_14j.csh:
setenv EMISINV_C "$CASEINPUTS/nonroad/2014NEIv1_nonroad_from_MOVES_forAQ_08dec2017_v9.csv.metals_pahs"