projects/brazil/config/UF/scripts/exhaust.R

In ibarraespinosa/vein: Vehicular Emissions Inventories

year               <- as.numeric(substr(x = getwd(), start = nchar(getwd()) - 6, stop = nchar(getwd()) - 3))

suppressWarnings(file.remove("emi/exhaust.csv"))
if(plot_ef){
  
switch (language,
        "portuguese" = cat("Plotando EF\n"),
        "english" = cat("Plotting EF\n"),
        "spanish" = cat("Plotando EF\n"))

# emission factors
n_PC <- metadata[metadata$family == "PC", ]$vehicles
n_LCV <- metadata[metadata$family == "LCV", ]$vehicles
n_TRUCKS <- metadata[metadata$family == "TRUCKS", ]$vehicles
n_BUS <- metadata[metadata$family == "BUS", ]$vehicles
n_MC <- metadata[metadata$family == "MC", ]$vehicles

ns <- c("PC", "LCV", "TRUCKS", "BUS", "MC")
ln  <- list(n_PC, n_LCV, n_TRUCKS, n_BUS, n_MC)

for(i in seq_along(ns)) {
  
  dl <- lapply(seq_along(pol), function(j){
    data.frame(ef_cetesb(p = pol[j], 
                         veh = ln[[i]],
                         year = year,
                         agemax = 40, 
                         sppm = metadata[metadata$vehicles %in%ln[[i]], ]$sppm,
                         verbose = verbose), 
               pol = pol[j])
  })
  
  dl <- do.call("rbind", dl)
  df <- wide_to_long(df = dl, 
                     column_with_data = ln[[i]], 
                     column_fixed = "pol")
  names(df) <- c("ef", "pol", "veh")
  df$age <- 1:40
  p <- ggplot(df[df$ef > 0, ], 
              aes(x = age, y = ef, colour = veh)) +
    geom_line() +
    facet_wrap(~pol, scales = "free_y") +
    ylim(0, NA) +
    labs(y = "g/km") +
    # scale_y_log10() +
    theme_bw()
  
  png(filename =  paste0("images/EF_", ns[i], ".png"),
      width = 2100, height = 1500, units ="px", pointsize = 12,
      bg = "white",  res = 300)
  print(p)
  dev.off()
}
switch (language,
        "portuguese" = message("\nFiguras em /images\n"),
        "english" = message("\nFigures in /image\n"),
        "spanish" = message("\nFiguras en /images\n"))

}

# Escapamento ####
switch (language,
        "portuguese" = cat("Estimando emissões\n"),
        "english" = cat("Estimating emissions\n"),
        "spanish" = cat("Estimando emisiones\n"))

# Exhaust ####
for(i in seq_along(metadata$vehicles)) {
  if(verbose) {
    cat("\n", metadata$vehicles[i], 
        rep("", max(nchar(metadata$vehicles) + 1) - nchar(metadata$vehicles[i])))
  }
  
  x <- readRDS(paste0("veh/", metadata$vehicles[i], ".rds"))
  
  for(j in seq_along(pol)){
    if(verbose) cat(pol[j], " ")
    
    ef <- ef_cetesb(p = pol[j], 
                    veh = metadata$vehicles[i], 
                    year = year,
                    agemax = ncol(x), 
                    scale = scale, 
                    sppm = metadata$sppm[i],
                    verbose = verbose)
    
    array_x <- emis_hot_td(veh = x, 
                           lkm = mileage[[metadata$vehicles[i]]], 
                           ef = ef, 
                           pro_month = as.numeric(pmonth[fuel == metadata$fuel[i]]$m3),
                           fortran = TRUE,
                           nt = check_nt()*0.9, 
                           verbose = verbose, 
                           params = list(veh = metadata$vehicles[i],
                                         size = metadata$size[i],
                                         fuel = metadata$fuel[i],
                                         pollutant = pol[j],
                                         type_emi = "Exhaust",
                                         subtype_emi = "Exhaust",
                                         baseyear = year))
    
    fwrite(array_x, "emi/exhaust.csv", append = TRUE)
    
  }
}


switch (language,
        "portuguese" = message("\n\nArquivos em: /emi/exhaust.csv:"),
        "english" = message("\n\nFiles in: /emi/exhaust.csv"),
        "spanish" = message("\n\nArchivos en: /emi/exhaust.csv"))


switch (language,
        "portuguese" = message("Limpando..."),
        "english" = message("Cleaning..."),
        "spanish" = message("Limpiando..."))

suppressWarnings(
  rm(i, j, pol, 
     n_PC, n_LCV, n_TRUCKS, n_BUS, n_MC,
     ns, ln, p, df, dl, cores) 
)

ls()   


gc()