to_wrf.RdFunction to expand, split and write emissions. The input is expanded into time by profile and split between variables with diferent weights.
to_wrf(POL, file = file.choose(), name = NA, total = NA, norm = F, profile = 1, weights = 1, verbose = T)
| 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 |
length(profile) must be the number of times in the emission file (value of frames_per_auxinput5 if wrf_create() was used to create this file).
total is an additional way to calculate or correct the total emissions
sum(profile) = 1 and sum(weights) = 1 to conserve mass
names and weights must have the same length
# 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") # }