eixport (version 0.4.0)

to_wrf: Combine total/spatial/temporal/split and write emission to file

Description

Function to expand, split and write emissions. The input is expanded into time by profile and split between variables with diferent weights.

Usage

to_wrf(POL, file = file.choose(), name = NA, total = NA, norm = F,
  profile = 1, weights = 1, verbose = T)

Arguments

POL

matrix or array of emissions of spatial weights

file

emission file name

name

species to be write

total

total of emited species (modifier)

norm

if the spatial weights need to be normalized (modifier)

profile

temporal profile to expand the emissions (modifier)

weights

weight of each species (modifier)

verbose

display additional information

See Also

wrf_create, wrf_get,wrf_profile and wrf_plot

Examples

# NOT RUN {
dir.create(file.path(tempdir(), "EMISS"))
wrf_create(wrfinput_dir = system.file("extdata", package = "eixport"),
          wrfchemi_dir = file.path(tempdir(), "EMISS"),
          frames_per_auxinput5 = 24)

# get the name of created file
files <- list.files(path = file.path(tempdir(), "EMISS"),
                   pattern = "wrfchemi",
                   full.names = TRUE)

data(Lights)

perfil <- c(0.010760058, 0.005280596, 0.002883553, 0.002666932,
            0.005781312, 0.018412838, 0.051900411, 0.077834636,
            0.067919758, 0.060831614, 0.055852868, 0.052468599,
            0.050938043, 0.051921718, 0.052756244, 0.052820165,
            0.058388406, 0.072855890, 0.075267137, 0.063246412,
            0.042713523, 0.029108975, 0.022091855, 0.015298458)

plot(perfil, ty = "l", col= "purple", xlab = "Hour", main = "Time profile",
    ylab = "Weight", axes = FALSE, xlim = c(0, 24))
axis(2)
axis(1, at = c(0, 6, 12, 18, 24),
    labels = c("00:00","06:00","12:00","18:00","00:00"))

to_wrf(Lights, files[1], total = 1521983, profile = perfil, name = "E_CO")
# }