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\_unit} ~ interp\_unsafe\left( \mathtt{ACET\_data}, 1 + \frac{ - \mathtt{ACET\_tstart} + t + \mathtt{t\_ref}}{\mathtt{ACET\_tstep}}, 1 + \frac{ - \mathtt{ACET\_s1start} + \mathtt{lon}}{\mathtt{ACET\_s1step}}, 1 + \frac{ - \mathtt{ACET\_s2start} + \mathtt{lat}}{\mathtt{ACET\_s2step}}, \mathtt{ACET\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_ACET} \right) \\ \mathtt{ACROLEIN}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{ACROLEIN\_unit} ~ interp\_unsafe\left( \mathtt{ACROLEIN\_data}, 1 + \frac{ - \mathtt{ACROLEIN\_tstart} + t + \mathtt{t\_ref}}{\mathtt{ACROLEIN\_tstep}}, 1 + \frac{ - \mathtt{ACROLEIN\_s1start} + \mathtt{lon}}{\mathtt{ACROLEIN\_s1step}}, 1 + \frac{ - \mathtt{ACROLEIN\_s2start} + \mathtt{lat}}{\mathtt{ACROLEIN\_s2step}}, \mathtt{ACROLEIN\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_ACROLEIN} \right) \\ \mathtt{ALD2}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{ALD2\_unit} ~ interp\_unsafe\left( \mathtt{ALD2\_data}, 1 + \frac{ - \mathtt{ALD2\_tstart} + t + \mathtt{t\_ref}}{\mathtt{ALD2\_tstep}}, 1 + \frac{ - \mathtt{ALD2\_s1start} + \mathtt{lon}}{\mathtt{ALD2\_s1step}}, 1 + \frac{ - \mathtt{ALD2\_s2start} + \mathtt{lat}}{\mathtt{ALD2\_s2step}}, \mathtt{ALD2\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_ALD2} \right) \\ \mathtt{ALD2\_PRIMARY}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{ALD2\_PRIMARY\_unit} ~ interp\_unsafe\left( \mathtt{ALD2\_PRIMARY\_data}, 1 + \frac{ - \mathtt{ALD2\_PRIMARY\_tstart} + t + \mathtt{t\_ref}}{\mathtt{ALD2\_PRIMARY\_tstep}}, 1 + \frac{ - \mathtt{ALD2\_PRIMARY\_s1start} + \mathtt{lon}}{\mathtt{ALD2\_PRIMARY\_s1step}}, 1 + \frac{ - \mathtt{ALD2\_PRIMARY\_s2start} + \mathtt{lat}}{\mathtt{ALD2\_PRIMARY\_s2step}}, \mathtt{ALD2\_PRIMARY\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_ALD2\_PRIMARY} \right) \\ \mathtt{ALDX}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{ALDX\_unit} ~ interp\_unsafe\left( \mathtt{ALDX\_data}, 1 + \frac{ - \mathtt{ALDX\_tstart} + t + \mathtt{t\_ref}}{\mathtt{ALDX\_tstep}}, 1 + \frac{ - \mathtt{ALDX\_s1start} + \mathtt{lon}}{\mathtt{ALDX\_s1step}}, 1 + \frac{ - \mathtt{ALDX\_s2start} + \mathtt{lat}}{\mathtt{ALDX\_s2step}}, \mathtt{ALDX\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_ALDX} \right) \\ \mathtt{APIN}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{APIN\_unit} ~ interp\_unsafe\left( \mathtt{APIN\_data}, 1 + \frac{ - \mathtt{APIN\_tstart} + t + \mathtt{t\_ref}}{\mathtt{APIN\_tstep}}, 1 + \frac{ - \mathtt{APIN\_s1start} + \mathtt{lon}}{\mathtt{APIN\_s1step}}, 1 + \frac{ - \mathtt{APIN\_s2start} + \mathtt{lat}}{\mathtt{APIN\_s2step}}, \mathtt{APIN\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_APIN} \right) \\ \mathtt{BENZ}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{BENZ\_unit} ~ interp\_unsafe\left( \mathtt{BENZ\_data}, 1 + \frac{ - \mathtt{BENZ\_tstart} + t + \mathtt{t\_ref}}{\mathtt{BENZ\_tstep}}, 1 + \frac{ - \mathtt{BENZ\_s1start} + \mathtt{lon}}{\mathtt{BENZ\_s1step}}, 1 + \frac{ - \mathtt{BENZ\_s2start} + \mathtt{lat}}{\mathtt{BENZ\_s2step}}, \mathtt{BENZ\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_BENZ} \right) \\ \mathtt{BPIN}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{BPIN\_unit} ~ interp\_unsafe\left( \mathtt{BPIN\_data}, 1 + \frac{ - \mathtt{BPIN\_tstart} + t + \mathtt{t\_ref}}{\mathtt{BPIN\_tstep}}, 1 + \frac{ - \mathtt{BPIN\_s1start} + \mathtt{lon}}{\mathtt{BPIN\_s1step}}, 1 + \frac{ - \mathtt{BPIN\_s2start} + \mathtt{lat}}{\mathtt{BPIN\_s2step}}, \mathtt{BPIN\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_BPIN} \right) \\ \mathtt{BUTADIENE13}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{BUTADIENE13\_unit} ~ interp\_unsafe\left( \mathtt{BUTADIENE13\_data}, 1 + \frac{ - \mathtt{BUTADIENE13\_tstart} + t + \mathtt{t\_ref}}{\mathtt{BUTADIENE13\_tstep}}, 1 + \frac{ - \mathtt{BUTADIENE13\_s1start} + \mathtt{lon}}{\mathtt{BUTADIENE13\_s1step}}, 1 + \frac{ - \mathtt{BUTADIENE13\_s2start} + \mathtt{lat}}{\mathtt{BUTADIENE13\_s2step}}, \mathtt{BUTADIENE13\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_BUTADIENE13} \right) \\ \mathtt{CH4}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{CH4\_unit} ~ interp\_unsafe\left( \mathtt{CH4\_data}, 1 + \frac{ - \mathtt{CH4\_tstart} + t + \mathtt{t\_ref}}{\mathtt{CH4\_tstep}}, 1 + \frac{ - \mathtt{CH4\_s1start} + \mathtt{lon}}{\mathtt{CH4\_s1step}}, 1 + \frac{ - \mathtt{CH4\_s2start} + \mathtt{lat}}{\mathtt{CH4\_s2step}}, \mathtt{CH4\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_CH4} \right) \\ \mathtt{CH4\_INV}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{CH4\_INV\_unit} ~ interp\_unsafe\left( \mathtt{CH4\_INV\_data}, 1 + \frac{ - \mathtt{CH4\_INV\_tstart} + t + \mathtt{t\_ref}}{\mathtt{CH4\_INV\_tstep}}, 1 + \frac{ - \mathtt{CH4\_INV\_s1start} + \mathtt{lon}}{\mathtt{CH4\_INV\_s1step}}, 1 + \frac{ - \mathtt{CH4\_INV\_s2start} + \mathtt{lat}}{\mathtt{CH4\_INV\_s2step}}, \mathtt{CH4\_INV\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_CH4\_INV} \right) \\ \mathtt{CL2}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{CL2\_unit} ~ interp\_unsafe\left( \mathtt{CL2\_data}, 1 + \frac{ - \mathtt{CL2\_tstart} + t + \mathtt{t\_ref}}{\mathtt{CL2\_tstep}}, 1 + \frac{ - \mathtt{CL2\_s1start} + \mathtt{lon}}{\mathtt{CL2\_s1step}}, 1 + \frac{ - \mathtt{CL2\_s2start} + \mathtt{lat}}{\mathtt{CL2\_s2step}}, \mathtt{CL2\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_CL2} \right) \\ \mathtt{CO}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{CO\_unit} ~ interp\_unsafe\left( \mathtt{CO\_data}, 1 + \frac{ - \mathtt{CO\_tstart} + t + \mathtt{t\_ref}}{\mathtt{CO\_tstep}}, 1 + \frac{ - \mathtt{CO\_s1start} + \mathtt{lon}}{\mathtt{CO\_s1step}}, 1 + \frac{ - \mathtt{CO\_s2start} + \mathtt{lat}}{\mathtt{CO\_s2step}}, \mathtt{CO\_extrap} \right) ~ diurnal\_itp\left( t + \mathtt{t\_ref}, \mathtt{lon} \right) ~ dayofweek\_itp\_CO\left( t + \mathtt{t\_ref}, \mathtt{lon} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_CO} \right) \\ \mathtt{CO2\_INV}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{CO2\_INV\_unit} ~ interp\_unsafe\left( \mathtt{CO2\_INV\_data}, 1 + \frac{ - \mathtt{CO2\_INV\_tstart} + t + \mathtt{t\_ref}}{\mathtt{CO2\_INV\_tstep}}, 1 + \frac{ - \mathtt{CO2\_INV\_s1start} + \mathtt{lon}}{\mathtt{CO2\_INV\_s1step}}, 1 + \frac{ - \mathtt{CO2\_INV\_s2start} + \mathtt{lat}}{\mathtt{CO2\_INV\_s2step}}, \mathtt{CO2\_INV\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_CO2\_INV} \right) \\ \mathtt{ETH}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{ETH\_unit} ~ interp\_unsafe\left( \mathtt{ETH\_data}, 1 + \frac{ - \mathtt{ETH\_tstart} + t + \mathtt{t\_ref}}{\mathtt{ETH\_tstep}}, 1 + \frac{ - \mathtt{ETH\_s1start} + \mathtt{lon}}{\mathtt{ETH\_s1step}}, 1 + \frac{ - \mathtt{ETH\_s2start} + \mathtt{lat}}{\mathtt{ETH\_s2step}}, \mathtt{ETH\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_ETH} \right) \\ \mathtt{ETHA}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{ETHA\_unit} ~ interp\_unsafe\left( \mathtt{ETHA\_data}, 1 + \frac{ - \mathtt{ETHA\_tstart} + t + \mathtt{t\_ref}}{\mathtt{ETHA\_tstep}}, 1 + \frac{ - \mathtt{ETHA\_s1start} + \mathtt{lon}}{\mathtt{ETHA\_s1step}}, 1 + \frac{ - \mathtt{ETHA\_s2start} + \mathtt{lat}}{\mathtt{ETHA\_s2step}}, \mathtt{ETHA\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_ETHA} \right) \\ \mathtt{ETHY}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{ETHY\_unit} ~ interp\_unsafe\left( \mathtt{ETHY\_data}, 1 + \frac{ - \mathtt{ETHY\_tstart} + t + \mathtt{t\_ref}}{\mathtt{ETHY\_tstep}}, 1 + \frac{ - \mathtt{ETHY\_s1start} + \mathtt{lon}}{\mathtt{ETHY\_s1step}}, 1 + \frac{ - \mathtt{ETHY\_s2start} + \mathtt{lat}}{\mathtt{ETHY\_s2step}}, \mathtt{ETHY\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_ETHY} \right) \\ \mathtt{ETOH}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{ETOH\_unit} ~ interp\_unsafe\left( \mathtt{ETOH\_data}, 1 + \frac{ - \mathtt{ETOH\_tstart} + t + \mathtt{t\_ref}}{\mathtt{ETOH\_tstep}}, 1 + \frac{ - \mathtt{ETOH\_s1start} + \mathtt{lon}}{\mathtt{ETOH\_s1step}}, 1 + \frac{ - \mathtt{ETOH\_s2start} + \mathtt{lat}}{\mathtt{ETOH\_s2step}}, \mathtt{ETOH\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_ETOH} \right) \\ \mathtt{FORM}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{FORM\_unit} ~ interp\_unsafe\left( \mathtt{FORM\_data}, 1 + \frac{ - \mathtt{FORM\_tstart} + t + \mathtt{t\_ref}}{\mathtt{FORM\_tstep}}, 1 + \frac{ - \mathtt{FORM\_s1start} + \mathtt{lon}}{\mathtt{FORM\_s1step}}, 1 + \frac{ - \mathtt{FORM\_s2start} + \mathtt{lat}}{\mathtt{FORM\_s2step}}, \mathtt{FORM\_extrap} \right) ~ diurnal\_itp\left( t + \mathtt{t\_ref}, \mathtt{lon} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_FORM} \right) \\ \mathtt{FORM\_PRIMARY}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{FORM\_PRIMARY\_unit} ~ interp\_unsafe\left( \mathtt{FORM\_PRIMARY\_data}, 1 + \frac{ - \mathtt{FORM\_PRIMARY\_tstart} + t + \mathtt{t\_ref}}{\mathtt{FORM\_PRIMARY\_tstep}}, 1 + \frac{ - \mathtt{FORM\_PRIMARY\_s1start} + \mathtt{lon}}{\mathtt{FORM\_PRIMARY\_s1step}}, 1 + \frac{ - \mathtt{FORM\_PRIMARY\_s2start} + \mathtt{lat}}{\mathtt{FORM\_PRIMARY\_s2step}}, \mathtt{FORM\_PRIMARY\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_FORM\_PRIMARY} \right) \\ \mathtt{HCL}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{HCL\_unit} ~ interp\_unsafe\left( \mathtt{HCL\_data}, 1 + \frac{ - \mathtt{HCL\_tstart} + t + \mathtt{t\_ref}}{\mathtt{HCL\_tstep}}, 1 + \frac{ - \mathtt{HCL\_s1start} + \mathtt{lon}}{\mathtt{HCL\_s1step}}, 1 + \frac{ - \mathtt{HCL\_s2start} + \mathtt{lat}}{\mathtt{HCL\_s2step}}, \mathtt{HCL\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_HCL} \right) \\ \mathtt{HONO}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{HONO\_unit} ~ interp\_unsafe\left( \mathtt{HONO\_data}, 1 + \frac{ - \mathtt{HONO\_tstart} + t + \mathtt{t\_ref}}{\mathtt{HONO\_tstep}}, 1 + \frac{ - \mathtt{HONO\_s1start} + \mathtt{lon}}{\mathtt{HONO\_s1step}}, 1 + \frac{ - \mathtt{HONO\_s2start} + \mathtt{lat}}{\mathtt{HONO\_s2step}}, \mathtt{HONO\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_HONO} \right) \\ \mathtt{IOLE}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{IOLE\_unit} ~ interp\_unsafe\left( \mathtt{IOLE\_data}, 1 + \frac{ - \mathtt{IOLE\_tstart} + t + \mathtt{t\_ref}}{\mathtt{IOLE\_tstep}}, 1 + \frac{ - \mathtt{IOLE\_s1start} + \mathtt{lon}}{\mathtt{IOLE\_s1step}}, 1 + \frac{ - \mathtt{IOLE\_s2start} + \mathtt{lat}}{\mathtt{IOLE\_s2step}}, \mathtt{IOLE\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_IOLE} \right) \\ \mathtt{ISOP}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{ISOP\_unit} ~ diurnal\_itp\_ISOP\left( t + \mathtt{t\_ref}, \mathtt{lon} \right) ~ interp\_unsafe\left( \mathtt{ISOP\_data}, 1 + \frac{ - \mathtt{ISOP\_tstart} + t + \mathtt{t\_ref}}{\mathtt{ISOP\_tstep}}, 1 + \frac{ - \mathtt{ISOP\_s1start} + \mathtt{lon}}{\mathtt{ISOP\_s1step}}, 1 + \frac{ - \mathtt{ISOP\_s2start} + \mathtt{lat}}{\mathtt{ISOP\_s2step}}, \mathtt{ISOP\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_ISOP} \right) \\ \mathtt{KET}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{KET\_unit} ~ interp\_unsafe\left( \mathtt{KET\_data}, 1 + \frac{ - \mathtt{KET\_tstart} + t + \mathtt{t\_ref}}{\mathtt{KET\_tstep}}, 1 + \frac{ - \mathtt{KET\_s1start} + \mathtt{lon}}{\mathtt{KET\_s1step}}, 1 + \frac{ - \mathtt{KET\_s2start} + \mathtt{lat}}{\mathtt{KET\_s2step}}, \mathtt{KET\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_KET} \right) \\ \mathtt{MEOH}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{MEOH\_unit} ~ interp\_unsafe\left( \mathtt{MEOH\_data}, 1 + \frac{ - \mathtt{MEOH\_tstart} + t + \mathtt{t\_ref}}{\mathtt{MEOH\_tstep}}, 1 + \frac{ - \mathtt{MEOH\_s1start} + \mathtt{lon}}{\mathtt{MEOH\_s1step}}, 1 + \frac{ - \mathtt{MEOH\_s2start} + \mathtt{lat}}{\mathtt{MEOH\_s2step}}, \mathtt{MEOH\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_MEOH} \right) \\ \mathtt{N2O\_INV}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{N2O\_INV\_unit} ~ interp\_unsafe\left( \mathtt{N2O\_INV\_data}, 1 + \frac{ - \mathtt{N2O\_INV\_tstart} + t + \mathtt{t\_ref}}{\mathtt{N2O\_INV\_tstep}}, 1 + \frac{ - \mathtt{N2O\_INV\_s1start} + \mathtt{lon}}{\mathtt{N2O\_INV\_s1step}}, 1 + \frac{ - \mathtt{N2O\_INV\_s2start} + \mathtt{lat}}{\mathtt{N2O\_INV\_s2step}}, \mathtt{N2O\_INV\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_N2O\_INV} \right) \\ \mathtt{NAPH}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{NAPH\_unit} ~ interp\_unsafe\left( \mathtt{NAPH\_data}, 1 + \frac{ - \mathtt{NAPH\_tstart} + t + \mathtt{t\_ref}}{\mathtt{NAPH\_tstep}}, 1 + \frac{ - \mathtt{NAPH\_s1start} + \mathtt{lon}}{\mathtt{NAPH\_s1step}}, 1 + \frac{ - \mathtt{NAPH\_s2start} + \mathtt{lat}}{\mathtt{NAPH\_s2step}}, \mathtt{NAPH\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_NAPH} \right) \\ \mathtt{NH3}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{NH3\_unit} ~ interp\_unsafe\left( \mathtt{NH3\_data}, 1 + \frac{ - \mathtt{NH3\_tstart} + t + \mathtt{t\_ref}}{\mathtt{NH3\_tstep}}, 1 + \frac{ - \mathtt{NH3\_s1start} + \mathtt{lon}}{\mathtt{NH3\_s1step}}, 1 + \frac{ - \mathtt{NH3\_s2start} + \mathtt{lat}}{\mathtt{NH3\_s2step}}, \mathtt{NH3\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_NH3} \right) \\ \mathtt{NH3\_FERT}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{NH3\_FERT\_unit} ~ interp\_unsafe\left( \mathtt{NH3\_FERT\_data}, 1 + \frac{ - \mathtt{NH3\_FERT\_tstart} + t + \mathtt{t\_ref}}{\mathtt{NH3\_FERT\_tstep}}, 1 + \frac{ - \mathtt{NH3\_FERT\_s1start} + \mathtt{lon}}{\mathtt{NH3\_FERT\_s1step}}, 1 + \frac{ - \mathtt{NH3\_FERT\_s2start} + \mathtt{lat}}{\mathtt{NH3\_FERT\_s2step}}, \mathtt{NH3\_FERT\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_NH3\_FERT} \right) \\ \mathtt{NO}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{NO\_unit} ~ dayofweek\_itp\_NOx\left( t + \mathtt{t\_ref}, \mathtt{lon} \right) ~ diurnal\_itp\_NOx\left( t + \mathtt{t\_ref}, \mathtt{lon} \right) ~ interp\_unsafe\left( \mathtt{NO\_data}, 1 + \frac{ - \mathtt{NO\_tstart} + t + \mathtt{t\_ref}}{\mathtt{NO\_tstep}}, 1 + \frac{ - \mathtt{NO\_s1start} + \mathtt{lon}}{\mathtt{NO\_s1step}}, 1 + \frac{ - \mathtt{NO\_s2start} + \mathtt{lat}}{\mathtt{NO\_s2step}}, \mathtt{NO\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_NO} \right) \\ \mathtt{NO2}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{NO2\_unit} ~ dayofweek\_itp\_NOx\left( t + \mathtt{t\_ref}, \mathtt{lon} \right) ~ diurnal\_itp\_NOx\left( t + \mathtt{t\_ref}, \mathtt{lon} \right) ~ interp\_unsafe\left( \mathtt{NO2\_data}, 1 + \frac{ - \mathtt{NO2\_tstart} + t + \mathtt{t\_ref}}{\mathtt{NO2\_tstep}}, 1 + \frac{ - \mathtt{NO2\_s1start} + \mathtt{lon}}{\mathtt{NO2\_s1step}}, 1 + \frac{ - \mathtt{NO2\_s2start} + \mathtt{lat}}{\mathtt{NO2\_s2step}}, \mathtt{NO2\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_NO2} \right) \\ \mathtt{NOX}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{NOX\_unit} ~ interp\_unsafe\left( \mathtt{NOX\_data}, 1 + \frac{ - \mathtt{NOX\_tstart} + t + \mathtt{t\_ref}}{\mathtt{NOX\_tstep}}, 1 + \frac{ - \mathtt{NOX\_s1start} + \mathtt{lon}}{\mathtt{NOX\_s1step}}, 1 + \frac{ - \mathtt{NOX\_s2start} + \mathtt{lat}}{\mathtt{NOX\_s2step}}, \mathtt{NOX\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_NOX} \right) \\ \mathtt{NR}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{NR\_unit} ~ interp\_unsafe\left( \mathtt{NR\_data}, 1 + \frac{ - \mathtt{NR\_tstart} + t + \mathtt{t\_ref}}{\mathtt{NR\_tstep}}, 1 + \frac{ - \mathtt{NR\_s1start} + \mathtt{lon}}{\mathtt{NR\_s1step}}, 1 + \frac{ - \mathtt{NR\_s2start} + \mathtt{lat}}{\mathtt{NR\_s2step}}, \mathtt{NR\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_NR} \right) \\ \mathtt{NVOL}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{NVOL\_unit} ~ interp\_unsafe\left( \mathtt{NVOL\_data}, 1 + \frac{ - \mathtt{NVOL\_tstart} + t + \mathtt{t\_ref}}{\mathtt{NVOL\_tstep}}, 1 + \frac{ - \mathtt{NVOL\_s1start} + \mathtt{lon}}{\mathtt{NVOL\_s1step}}, 1 + \frac{ - \mathtt{NVOL\_s2start} + \mathtt{lat}}{\mathtt{NVOL\_s2step}}, \mathtt{NVOL\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_NVOL} \right) \\ \mathtt{OLE}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{OLE\_unit} ~ interp\_unsafe\left( \mathtt{OLE\_data}, 1 + \frac{ - \mathtt{OLE\_tstart} + t + \mathtt{t\_ref}}{\mathtt{OLE\_tstep}}, 1 + \frac{ - \mathtt{OLE\_s1start} + \mathtt{lon}}{\mathtt{OLE\_s1step}}, 1 + \frac{ - \mathtt{OLE\_s2start} + \mathtt{lat}}{\mathtt{OLE\_s2step}}, \mathtt{OLE\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_OLE} \right) \\ \mathtt{PAL}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PAL\_unit} ~ interp\_unsafe\left( \mathtt{PAL\_data}, 1 + \frac{ - \mathtt{PAL\_tstart} + t + \mathtt{t\_ref}}{\mathtt{PAL\_tstep}}, 1 + \frac{ - \mathtt{PAL\_s1start} + \mathtt{lon}}{\mathtt{PAL\_s1step}}, 1 + \frac{ - \mathtt{PAL\_s2start} + \mathtt{lat}}{\mathtt{PAL\_s2step}}, \mathtt{PAL\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_PAL} \right) \\ \mathtt{PAR}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PAR\_unit} ~ interp\_unsafe\left( \mathtt{PAR\_data}, 1 + \frac{ - \mathtt{PAR\_tstart} + t + \mathtt{t\_ref}}{\mathtt{PAR\_tstep}}, 1 + \frac{ - \mathtt{PAR\_s1start} + \mathtt{lon}}{\mathtt{PAR\_s1step}}, 1 + \frac{ - \mathtt{PAR\_s2start} + \mathtt{lat}}{\mathtt{PAR\_s2step}}, \mathtt{PAR\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_PAR} \right) \\ \mathtt{PCA}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PCA\_unit} ~ interp\_unsafe\left( \mathtt{PCA\_data}, 1 + \frac{ - \mathtt{PCA\_tstart} + t + \mathtt{t\_ref}}{\mathtt{PCA\_tstep}}, 1 + \frac{ - \mathtt{PCA\_s1start} + \mathtt{lon}}{\mathtt{PCA\_s1step}}, 1 + \frac{ - \mathtt{PCA\_s2start} + \mathtt{lat}}{\mathtt{PCA\_s2step}}, \mathtt{PCA\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_PCA} \right) \\ \mathtt{PCL}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PCL\_unit} ~ interp\_unsafe\left( \mathtt{PCL\_data}, 1 + \frac{ - \mathtt{PCL\_tstart} + t + \mathtt{t\_ref}}{\mathtt{PCL\_tstep}}, 1 + \frac{ - \mathtt{PCL\_s1start} + \mathtt{lon}}{\mathtt{PCL\_s1step}}, 1 + \frac{ - \mathtt{PCL\_s2start} + \mathtt{lat}}{\mathtt{PCL\_s2step}}, \mathtt{PCL\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_PCL} \right) \\ \mathtt{PEC}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PEC\_unit} ~ interp\_unsafe\left( \mathtt{PEC\_data}, 1 + \frac{ - \mathtt{PEC\_tstart} + t + \mathtt{t\_ref}}{\mathtt{PEC\_tstep}}, 1 + \frac{ - \mathtt{PEC\_s1start} + \mathtt{lon}}{\mathtt{PEC\_s1step}}, 1 + \frac{ - \mathtt{PEC\_s2start} + \mathtt{lat}}{\mathtt{PEC\_s2step}}, \mathtt{PEC\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_PEC} \right) \\ \mathtt{PFE}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PFE\_unit} ~ interp\_unsafe\left( \mathtt{PFE\_data}, 1 + \frac{ - \mathtt{PFE\_tstart} + t + \mathtt{t\_ref}}{\mathtt{PFE\_tstep}}, 1 + \frac{ - \mathtt{PFE\_s1start} + \mathtt{lon}}{\mathtt{PFE\_s1step}}, 1 + \frac{ - \mathtt{PFE\_s2start} + \mathtt{lat}}{\mathtt{PFE\_s2step}}, \mathtt{PFE\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_PFE} \right) \\ \mathtt{PH2O}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PH2O\_unit} ~ interp\_unsafe\left( \mathtt{PH2O\_data}, 1 + \frac{ - \mathtt{PH2O\_tstart} + t + \mathtt{t\_ref}}{\mathtt{PH2O\_tstep}}, 1 + \frac{ - \mathtt{PH2O\_s1start} + \mathtt{lon}}{\mathtt{PH2O\_s1step}}, 1 + \frac{ - \mathtt{PH2O\_s2start} + \mathtt{lat}}{\mathtt{PH2O\_s2step}}, \mathtt{PH2O\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_PH2O} \right) \\ \mathtt{PK}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PK\_unit} ~ interp\_unsafe\left( \mathtt{PK\_data}, 1 + \frac{ - \mathtt{PK\_tstart} + t + \mathtt{t\_ref}}{\mathtt{PK\_tstep}}, 1 + \frac{ - \mathtt{PK\_s1start} + \mathtt{lon}}{\mathtt{PK\_s1step}}, 1 + \frac{ - \mathtt{PK\_s2start} + \mathtt{lat}}{\mathtt{PK\_s2step}}, \mathtt{PK\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_PK} \right) \\ \mathtt{PM2\_5}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PM2\_5\_unit} ~ interp\_unsafe\left( \mathtt{PM2\_5\_data}, 1 + \frac{ - \mathtt{PM2\_5\_tstart} + t + \mathtt{t\_ref}}{\mathtt{PM2\_5\_tstep}}, 1 + \frac{ - \mathtt{PM2\_5\_s1start} + \mathtt{lon}}{\mathtt{PM2\_5\_s1step}}, 1 + \frac{ - \mathtt{PM2\_5\_s2start} + \mathtt{lat}}{\mathtt{PM2\_5\_s2step}}, \mathtt{PM2\_5\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_PM2\_5} \right) \\ \mathtt{PMC}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PMC\_unit} ~ interp\_unsafe\left( \mathtt{PMC\_data}, 1 + \frac{ - \mathtt{PMC\_tstart} + t + \mathtt{t\_ref}}{\mathtt{PMC\_tstep}}, 1 + \frac{ - \mathtt{PMC\_s1start} + \mathtt{lon}}{\mathtt{PMC\_s1step}}, 1 + \frac{ - \mathtt{PMC\_s2start} + \mathtt{lat}}{\mathtt{PMC\_s2step}}, \mathtt{PMC\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_PMC} \right) \\ \mathtt{PMG}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PMG\_unit} ~ interp\_unsafe\left( \mathtt{PMG\_data}, 1 + \frac{ - \mathtt{PMG\_tstart} + t + \mathtt{t\_ref}}{\mathtt{PMG\_tstep}}, 1 + \frac{ - \mathtt{PMG\_s1start} + \mathtt{lon}}{\mathtt{PMG\_s1step}}, 1 + \frac{ - \mathtt{PMG\_s2start} + \mathtt{lat}}{\mathtt{PMG\_s2step}}, \mathtt{PMG\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_PMG} \right) \\ \mathtt{PMN}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PMN\_unit} ~ interp\_unsafe\left( \mathtt{PMN\_data}, 1 + \frac{ - \mathtt{PMN\_tstart} + t + \mathtt{t\_ref}}{\mathtt{PMN\_tstep}}, 1 + \frac{ - \mathtt{PMN\_s1start} + \mathtt{lon}}{\mathtt{PMN\_s1step}}, 1 + \frac{ - \mathtt{PMN\_s2start} + \mathtt{lat}}{\mathtt{PMN\_s2step}}, \mathtt{PMN\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_PMN} \right) \\ \mathtt{PMOTHR}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PMOTHR\_unit} ~ interp\_unsafe\left( \mathtt{PMOTHR\_data}, 1 + \frac{ - \mathtt{PMOTHR\_tstart} + t + \mathtt{t\_ref}}{\mathtt{PMOTHR\_tstep}}, 1 + \frac{ - \mathtt{PMOTHR\_s1start} + \mathtt{lon}}{\mathtt{PMOTHR\_s1step}}, 1 + \frac{ - \mathtt{PMOTHR\_s2start} + \mathtt{lat}}{\mathtt{PMOTHR\_s2step}}, \mathtt{PMOTHR\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_PMOTHR} \right) \\ \mathtt{PNA}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PNA\_unit} ~ interp\_unsafe\left( \mathtt{PNA\_data}, 1 + \frac{ - \mathtt{PNA\_tstart} + t + \mathtt{t\_ref}}{\mathtt{PNA\_tstep}}, 1 + \frac{ - \mathtt{PNA\_s1start} + \mathtt{lon}}{\mathtt{PNA\_s1step}}, 1 + \frac{ - \mathtt{PNA\_s2start} + \mathtt{lat}}{\mathtt{PNA\_s2step}}, \mathtt{PNA\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_PNA} \right) \\ \mathtt{PNCOM}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PNCOM\_unit} ~ interp\_unsafe\left( \mathtt{PNCOM\_data}, 1 + \frac{ - \mathtt{PNCOM\_tstart} + t + \mathtt{t\_ref}}{\mathtt{PNCOM\_tstep}}, 1 + \frac{ - \mathtt{PNCOM\_s1start} + \mathtt{lon}}{\mathtt{PNCOM\_s1step}}, 1 + \frac{ - \mathtt{PNCOM\_s2start} + \mathtt{lat}}{\mathtt{PNCOM\_s2step}}, \mathtt{PNCOM\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_PNCOM} \right) \\ \mathtt{PNH4}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PNH4\_unit} ~ interp\_unsafe\left( \mathtt{PNH4\_data}, 1 + \frac{ - \mathtt{PNH4\_tstart} + t + \mathtt{t\_ref}}{\mathtt{PNH4\_tstep}}, 1 + \frac{ - \mathtt{PNH4\_s1start} + \mathtt{lon}}{\mathtt{PNH4\_s1step}}, 1 + \frac{ - \mathtt{PNH4\_s2start} + \mathtt{lat}}{\mathtt{PNH4\_s2step}}, \mathtt{PNH4\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_PNH4} \right) \\ \mathtt{PNO3}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PNO3\_unit} ~ interp\_unsafe\left( \mathtt{PNO3\_data}, 1 + \frac{ - \mathtt{PNO3\_tstart} + t + \mathtt{t\_ref}}{\mathtt{PNO3\_tstep}}, 1 + \frac{ - \mathtt{PNO3\_s1start} + \mathtt{lon}}{\mathtt{PNO3\_s1step}}, 1 + \frac{ - \mathtt{PNO3\_s2start} + \mathtt{lat}}{\mathtt{PNO3\_s2step}}, \mathtt{PNO3\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_PNO3} \right) \\ \mathtt{POC}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{POC\_unit} ~ interp\_unsafe\left( \mathtt{POC\_data}, 1 + \frac{ - \mathtt{POC\_tstart} + t + \mathtt{t\_ref}}{\mathtt{POC\_tstep}}, 1 + \frac{ - \mathtt{POC\_s1start} + \mathtt{lon}}{\mathtt{POC\_s1step}}, 1 + \frac{ - \mathtt{POC\_s2start} + \mathtt{lat}}{\mathtt{POC\_s2step}}, \mathtt{POC\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_POC} \right) \\ \mathtt{PRPA}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PRPA\_unit} ~ interp\_unsafe\left( \mathtt{PRPA\_data}, 1 + \frac{ - \mathtt{PRPA\_tstart} + t + \mathtt{t\_ref}}{\mathtt{PRPA\_tstep}}, 1 + \frac{ - \mathtt{PRPA\_s1start} + \mathtt{lon}}{\mathtt{PRPA\_s1step}}, 1 + \frac{ - \mathtt{PRPA\_s2start} + \mathtt{lat}}{\mathtt{PRPA\_s2step}}, \mathtt{PRPA\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_PRPA} \right) \\ \mathtt{PSI}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PSI\_unit} ~ interp\_unsafe\left( \mathtt{PSI\_data}, 1 + \frac{ - \mathtt{PSI\_tstart} + t + \mathtt{t\_ref}}{\mathtt{PSI\_tstep}}, 1 + \frac{ - \mathtt{PSI\_s1start} + \mathtt{lon}}{\mathtt{PSI\_s1step}}, 1 + \frac{ - \mathtt{PSI\_s2start} + \mathtt{lat}}{\mathtt{PSI\_s2step}}, \mathtt{PSI\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_PSI} \right) \\ \mathtt{PSO4}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PSO4\_unit} ~ interp\_unsafe\left( \mathtt{PSO4\_data}, 1 + \frac{ - \mathtt{PSO4\_tstart} + t + \mathtt{t\_ref}}{\mathtt{PSO4\_tstep}}, 1 + \frac{ - \mathtt{PSO4\_s1start} + \mathtt{lon}}{\mathtt{PSO4\_s1step}}, 1 + \frac{ - \mathtt{PSO4\_s2start} + \mathtt{lat}}{\mathtt{PSO4\_s2step}}, \mathtt{PSO4\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_PSO4} \right) \\ \mathtt{PTI}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{PTI\_unit} ~ interp\_unsafe\left( \mathtt{PTI\_data}, 1 + \frac{ - \mathtt{PTI\_tstart} + t + \mathtt{t\_ref}}{\mathtt{PTI\_tstep}}, 1 + \frac{ - \mathtt{PTI\_s1start} + \mathtt{lon}}{\mathtt{PTI\_s1step}}, 1 + \frac{ - \mathtt{PTI\_s2start} + \mathtt{lat}}{\mathtt{PTI\_s2step}}, \mathtt{PTI\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_PTI} \right) \\ \mathtt{SESQ}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{SESQ\_unit} ~ interp\_unsafe\left( \mathtt{SESQ\_data}, 1 + \frac{ - \mathtt{SESQ\_tstart} + t + \mathtt{t\_ref}}{\mathtt{SESQ\_tstep}}, 1 + \frac{ - \mathtt{SESQ\_s1start} + \mathtt{lon}}{\mathtt{SESQ\_s1step}}, 1 + \frac{ - \mathtt{SESQ\_s2start} + \mathtt{lat}}{\mathtt{SESQ\_s2step}}, \mathtt{SESQ\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_SESQ} \right) \\ \mathtt{SO2}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{SO2\_unit} ~ interp\_unsafe\left( \mathtt{SO2\_data}, 1 + \frac{ - \mathtt{SO2\_tstart} + t + \mathtt{t\_ref}}{\mathtt{SO2\_tstep}}, 1 + \frac{ - \mathtt{SO2\_s1start} + \mathtt{lon}}{\mathtt{SO2\_s1step}}, 1 + \frac{ - \mathtt{SO2\_s2start} + \mathtt{lat}}{\mathtt{SO2\_s2step}}, \mathtt{SO2\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_SO2} \right) \\ \mathtt{SOAALK}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{SOAALK\_unit} ~ interp\_unsafe\left( \mathtt{SOAALK\_data}, 1 + \frac{ - \mathtt{SOAALK\_tstart} + t + \mathtt{t\_ref}}{\mathtt{SOAALK\_tstep}}, 1 + \frac{ - \mathtt{SOAALK\_s1start} + \mathtt{lon}}{\mathtt{SOAALK\_s1step}}, 1 + \frac{ - \mathtt{SOAALK\_s2start} + \mathtt{lat}}{\mathtt{SOAALK\_s2step}}, \mathtt{SOAALK\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_SOAALK} \right) \\ \mathtt{SULF}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{SULF\_unit} ~ interp\_unsafe\left( \mathtt{SULF\_data}, 1 + \frac{ - \mathtt{SULF\_tstart} + t + \mathtt{t\_ref}}{\mathtt{SULF\_tstep}}, 1 + \frac{ - \mathtt{SULF\_s1start} + \mathtt{lon}}{\mathtt{SULF\_s1step}}, 1 + \frac{ - \mathtt{SULF\_s2start} + \mathtt{lat}}{\mathtt{SULF\_s2step}}, \mathtt{SULF\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_SULF} \right) \\ \mathtt{TERP}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{TERP\_unit} ~ interp\_unsafe\left( \mathtt{TERP\_data}, 1 + \frac{ - \mathtt{TERP\_tstart} + t + \mathtt{t\_ref}}{\mathtt{TERP\_tstep}}, 1 + \frac{ - \mathtt{TERP\_s1start} + \mathtt{lon}}{\mathtt{TERP\_s1step}}, 1 + \frac{ - \mathtt{TERP\_s2start} + \mathtt{lat}}{\mathtt{TERP\_s2step}}, \mathtt{TERP\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_TERP} \right) \\ \mathtt{TOL}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{TOL\_unit} ~ interp\_unsafe\left( \mathtt{TOL\_data}, 1 + \frac{ - \mathtt{TOL\_tstart} + t + \mathtt{t\_ref}}{\mathtt{TOL\_tstep}}, 1 + \frac{ - \mathtt{TOL\_s1start} + \mathtt{lon}}{\mathtt{TOL\_s1step}}, 1 + \frac{ - \mathtt{TOL\_s2start} + \mathtt{lat}}{\mathtt{TOL\_s2step}}, \mathtt{TOL\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_TOL} \right) \\ \mathtt{UNK}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{UNK\_unit} ~ interp\_unsafe\left( \mathtt{UNK\_data}, 1 + \frac{ - \mathtt{UNK\_tstart} + t + \mathtt{t\_ref}}{\mathtt{UNK\_tstep}}, 1 + \frac{ - \mathtt{UNK\_s1start} + \mathtt{lon}}{\mathtt{UNK\_s1step}}, 1 + \frac{ - \mathtt{UNK\_s2start} + \mathtt{lat}}{\mathtt{UNK\_s2step}}, \mathtt{UNK\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_UNK} \right) \\ \mathtt{UNR}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{UNR\_unit} ~ interp\_unsafe\left( \mathtt{UNR\_data}, 1 + \frac{ - \mathtt{UNR\_tstart} + t + \mathtt{t\_ref}}{\mathtt{UNR\_tstep}}, 1 + \frac{ - \mathtt{UNR\_s1start} + \mathtt{lon}}{\mathtt{UNR\_s1step}}, 1 + \frac{ - \mathtt{UNR\_s2start} + \mathtt{lat}}{\mathtt{UNR\_s2step}}, \mathtt{UNR\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_UNR} \right) \\ \mathtt{VOC\_BEIS}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{VOC\_BEIS\_unit} ~ interp\_unsafe\left( \mathtt{VOC\_BEIS\_data}, 1 + \frac{ - \mathtt{VOC\_BEIS\_tstart} + t + \mathtt{t\_ref}}{\mathtt{VOC\_BEIS\_tstep}}, 1 + \frac{ - \mathtt{VOC\_BEIS\_s1start} + \mathtt{lon}}{\mathtt{VOC\_BEIS\_s1step}}, 1 + \frac{ - \mathtt{VOC\_BEIS\_s2start} + \mathtt{lat}}{\mathtt{VOC\_BEIS\_s2step}}, \mathtt{VOC\_BEIS\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_VOC\_BEIS} \right) \\ \mathtt{VOC\_INV}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{VOC\_INV\_unit} ~ interp\_unsafe\left( \mathtt{VOC\_INV\_data}, 1 + \frac{ - \mathtt{VOC\_INV\_tstart} + t + \mathtt{t\_ref}}{\mathtt{VOC\_INV\_tstep}}, 1 + \frac{ - \mathtt{VOC\_INV\_s1start} + \mathtt{lon}}{\mathtt{VOC\_INV\_s1step}}, 1 + \frac{ - \mathtt{VOC\_INV\_s2start} + \mathtt{lat}}{\mathtt{VOC\_INV\_s2step}}, \mathtt{VOC\_INV\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_VOC\_INV} \right) \\ \mathtt{XYLMN}\left( t \right) &= ifelse\left( \mathtt{lev} < 2, \frac{\mathtt{XYLMN\_unit} ~ interp\_unsafe\left( \mathtt{XYLMN\_data}, 1 + \frac{ - \mathtt{XYLMN\_tstart} + t + \mathtt{t\_ref}}{\mathtt{XYLMN\_tstep}}, 1 + \frac{ - \mathtt{XYLMN\_s1start} + \mathtt{lon}}{\mathtt{XYLMN\_s1step}}, 1 + \frac{ - \mathtt{XYLMN\_s2start} + \mathtt{lat}}{\mathtt{XYLMN\_s2step}}, \mathtt{XYLMN\_extrap} \right)}{\mathtt{g0\_100} ~ delp\_dry\_surface\_itp\left( \mathtt{lon}, \mathtt{lat} \right) ~ \mathtt{{\Delta}z}}, \mathtt{zero\_XYLMN} \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))