emission: Emissions in the format for atmospheric models

In Schuch666/EmissV: Tools for Create Emissions for Air Quality Models

Emissions in the format for atmospheric models

Description

Combine area sources and total emissions to model output

Usage

emission(
  inventory = NULL,
  grid,
  mm = 1,
  aerosol = FALSE,
  check = TRUE,
  total,
  pol,
  area,
  plot = FALSE,
  verbose = TRUE
)

Arguments

inventory	a inventory raster from read
grid	grid information
mm	pollutant molar mass
aerosol	TRUE for aerosols and FALSE (defoult) for gazes
check	TRUE (defoult) to check negative and NA values and replace it for zero
total	list of total emission
pol	pollutant name
area	list of area sources or matrix with a spatial distribution
plot	TRUE for plot the final emissions
verbose	display additional information

Format

matrix of emission

Value

a vector of emissions in MOL / mk2 h for gases and ug / m2 s for aerosols.

Note

if Inventory is provided, the firsts tree arguments are not be used by the function.

Is a good practice use the set_units(fe,your_unity), where fe is your emission factory and your_unity is usually g/km on your emission factory

the list of area must be in the same order as defined in vehicles and total emission.

just WRF-Chem is suported by now

See Also

totalEmission and areaSource

Examples

fleet  <- vehicles(example = TRUE)

EmissionFactors <- emissionFactor(example = TRUE)

TOTAL  <- totalEmission(fleet,EmissionFactors,pol = c("CO"),verbose = TRUE)

grid   <- gridInfo(paste0(system.file("extdata", package = "EmissV"),"/wrfinput_d01"))
shape  <- raster::shapefile(paste0(system.file("extdata", package = "EmissV"),"/BR.shp"))
raster <- raster::raster(paste0(system.file("extdata", package = "EmissV"),"/dmsp.tiff"))

SP     <- areaSource(shape[22,1],raster,grid,name = "SP")
RJ     <- areaSource(shape[17,1],raster,grid,name = "RJ")

e_CO   <- emission(total = TOTAL,
                   pol = "CO",
                   area = list(SP = SP, RJ = RJ),
                   grid = grid,
                   mm = 28)

Schuch666/EmissV documentation built on July 19, 2023, 10:36 a.m.