emission: Emissions in the format for atmospheric models

View source: R/emission.R

emissionR Documentation

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 Dec. 24, 2024, 3:15 a.m.