2016 US EPA National Emissions Inventory (NEI) data

We have a data loader for CMAQ-formatted monthly US National Emissions Inventory data for year 2016,NEI2016MonthlyEmis.

Download Configuration

Because there is an issue with the EPA's FTP server that we download the data from you may need to set the following environment variable before using it:

In Julia:

ENV["JULIA_NO_VERIFY_HOSTS"] = "gaftp.epa.gov"

or in a bash shell:

export JULIA_NO_VERIFY_HOSTS=gaftp.epa.gov

Equations

This is what its equation system looks like:

using EarthSciData, EarthSciMLBase
using ModelingToolkit, DynamicQuantities, DataFrames
using ModelingToolkit: t
using DynamicQuantities: dimension
using Dates

domain = DomainInfo(
    DateTime(2016, 5, 1), DateTime(2016, 5, 2);
    lonrange = deg2rad(-115):deg2rad(2.5):deg2rad(-68.75),
    latrange = deg2rad(25):deg2rad(2):deg2rad(53.7),
    levrange = 1:10,
    u_proto = zeros(Float32, 1, 1, 1, 1)
)

emis = NEI2016MonthlyEmis("mrggrid_withbeis_withrwc", domain)

\[ \begin{align} \mathtt{ACET}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{ACET\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{ACROLEIN}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{ACROLEIN\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{ALD2}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{ALD2\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{ALD2\_PRIMARY}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{ALD2\_PRIMARY\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{ALDX}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{ALDX\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{APIN}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{APIN\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{BENZ}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{BENZ\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{BPIN}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{BPIN\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{BUTADIENE13}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{BUTADIENE13\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{CH4}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{CH4\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{CH4\_INV}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{CH4\_INV\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{CL2}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{CL2\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{CO}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{CO\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{CO2\_INV}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{CO2\_INV\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{ETH}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{ETH\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{ETHA}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{ETHA\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{ETHY}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{ETHY\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{ETOH}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{ETOH\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{FORM}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{FORM\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{FORM\_PRIMARY}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{FORM\_PRIMARY\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{HCL}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{HCL\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{HONO}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{HONO\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{IOLE}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{IOLE\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{ISOP}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{ISOP\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{KET}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{KET\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{MEOH}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{MEOH\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{N2O\_INV}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{N2O\_INV\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{NAPH}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{NAPH\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{NH3}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{NH3\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{NH3\_FERT}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{NH3\_FERT\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{NO}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{NO\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{NO2}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{NO2\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{NOX}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{NOX\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{NR}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{NR\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{NVOL}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{NVOL\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{OLE}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{OLE\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{PAL}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PAL\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{PAR}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PAR\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{PCA}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PCA\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{PCL}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PCL\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{PEC}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PEC\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{PFE}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PFE\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{PH2O}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PH2O\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{PK}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PK\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{PM2\_5}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PM2\_5\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{PMC}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PMC\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{PMG}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PMG\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{PMN}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PMN\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{PMOTHR}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PMOTHR\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{PNA}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PNA\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{PNCOM}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PNCOM\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{PNH4}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PNH4\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{PNO3}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PNO3\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{POC}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{POC\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{PRPA}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PRPA\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{PSI}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PSI\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{PSO4}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PSO4\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{PTI}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PTI\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{SESQ}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{SESQ\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{SO2}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{SO2\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{SOAALK}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{SOAALK\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{SULF}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{SULF\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{TERP}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{TERP\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{TOL}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{TOL\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{UNK}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{UNK\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{UNR}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{UNR\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{VOC\_BEIS}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{VOC\_BEIS\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{VOC\_INV}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{VOC\_INV\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \\ \mathtt{XYLMN}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{XYLMN\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right)}{\mathtt{{\Delta}z}}, \mathtt{zero\_emis} \right) \end{align} \]

Variables

Here are the variables in tabular format:

function table(vars)
    DataFrame(
        :Name => [string(Symbolics.tosymbol(v, escape = false)) for v in vars],
        :Units => [dimension(ModelingToolkit.get_unit(v)) for v in vars],
        :Description => [ModelingToolkit.getdescription(v) for v in vars]
    )
end
table(unknowns(emis))

Parameters

Finally, here are the parameters in tabular format:

table(parameters(emis))