projects/brazil_td_chem/scripts/ressuspensao.R

In ibarraespinosa/vein: Vehicular Emissions Inventories

poly <- pol

# Exhaust ####
for (i in seq_along(metadata$vehicles)) {
   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)) {
      cat(poly[j], " ")

      ef <- ef_cetesb(
         p = poly[j],
         veh = metadata$vehicles[i],
         year = year,
         agemax = ncol(x),
         verbose = verbose
      )

      array_x <- emis_hot_td(
         veh = x,
         lkm = mileage[[metadata$vehicles[i]]],
         ef = ef,
         pro_month = pmonth[[metadata$vehicles[i]]],
         fortran = TRUE,
         nt = check_nt() / 2,
         verbose = verbose,
         params = list(
            veh = metadata$vehicles[i],
            size = metadata$size[i],
            fuel = metadata$fuel[i],
            pollutant = pol[j],
            type_emi = "Paved Roads",
            subtype_emi = "Ressuspenssion",
            baseyear = year
         )
      )

      fwrite(array_x, "emi/exhaust.csv", append = TRUE)
   }
}


suppressWarnings(
   rm(
      ADT,
      PC, LCV, TRUCKS, BUS, MC,
      n_PC, n_LCV, n_TRUCKS, n_BUS, n_MC,
      tf_PC, tf_LCV, tf_TRUCKS, tf_BUS, tf_MC,
      emi, g, lkm,
      metadata, tfs, veh,
      net, nveh, pol,
      verbose,
      vk, vkbus, vkhgv, vklcv, vkmc, vkpc, W, x_DF,
      ef_cetesb2, ef_evaps, i, j, k, mileage, name_file_evap, num_vein,
      q, type_emis, vein_version, x, cores, fuel,
      sL1, sL2, sL3, sL4
   )
)