emission: Emissions in the format for atmospheric models

Description Usage Arguments Format Note See Also Examples

View source: R/emission.R

Description

Combine area sources and total emissions to model output

Usage

1
2
emission(total, pol, area, grid, inventory = NULL, mm = 1,
  aerosol = F, plot = F, verbose = T)

Arguments

total

list of total emission

pol

pollutant name

area

list of area sources or matrix with a spatial distribution

grid

grid information

inventory

a inventory raster from read

mm

pollutant molar mass

aerosol

TRUE for aerosols and FALSE (defoult) for gazes

plot

TRUE for plot the final emissions

verbose

display additional information

Format

matrix of emission

Note

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

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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
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,"CO",list(SP = SP, RJ = RJ),grid,mm=28)

EmissV documentation built on May 2, 2019, 7:16 a.m.