pointSource {EmissV}

Emissions from point sources

Description

Transform a set of points into a grinded output

Usage

pointSource(emissions, grid, verbose = TRUE)

Arguments

emissions

list of points

grid

grid object with the grid information

verbose

display additional information

Value

a raster

See Also

gridInfo and rasterSource

Examples

d1 <- gridInfo(paste(system.file("extdata", package = "EmissV"),"/wrfinput_d01",sep=""))

p = data.frame(lat      = c(-22,-22,-23.5),
               lon      = c(-46,-48,-47  ),
               z        = c(0  ,  0,  0  ),
               emission = c(666,444,111  ) )

p_emissions <- pointSource(emissions = p, grid = d1)

raster::plot(p_emissions,ylab="Lat", xlab="Lon",
             main = "3 point sources for domain d1")