R/chem_edgar.R

In eixport: Export Emissions to Atmospheric Models

#' Agregate EDGAR emissions NetCDF files into a RasterStack by
#'
#' @description The Emissions Database for Global Atmospheric Research (EDGAR) is a
#' project from the Joint Research Centre. This function reads the NetCDF and
#' merge/aggregate into diferent chemical mechanisms
#'
#'
#' @param path Character; path to the NetCDF files from EDGAR. The directory
#' **must** have one file for each of the following pollutants:
#' "voc" from 1 to 25, "co", "nox", "nmvoc","so2", "nh3",
#' "pm10", "pm2.5", "bc" and "oc"
#' @param chem Character; chemical mechanism: "edgar", "radm", "radmsorg",
#' "cbmz_mosaic", "cptec", "ecb05_opt1",
#' "neu_cb05" (thanks to Daniel Schuch) and
#' "ufpr_cbmz" (thanks to Leila Martins).
#'
#' - When chem is "edgar" units are: "g km-2 h-1"
#' - Other mechanisms: gases "mol km-2 h-1" and aerosols: "ug m-2 s-1"
#' @param merge Logical; in the case that tehre are more than one NetCDF per pollutant,
#' merge = TRUE will merge them with sum. Default is FALSE.
#' @param k, Numeric; Value to factorize each pollutant.
#' @param verbose Logical to print more information
#' @return RasterStack
#' @importFrom raster stack raster rotate
#' @importFrom ncdf4 nc_open nc_close
#' @note Molecular weights were obtained from
#'
#' Development of  Improved Chemical Speciation Database for Processing Emissions of
#' Volatile Organic Compounds for Air Quality Models
#' https://intra.engr.ucr.edu/~carter/emitdb/
#'
#' Some mappings were obtained from:
#'
#' Carter, W. P. (2015). Development of a database for chemical
#' mechanism assignments for volatile organic emissions.
#' Journal of the Air & Waste Management Association, 65(10), 1171-1184.
#'
#' Lopez-Norena, Ana and Fernandez, Rafael & Puliafito, SALVADOR. (2019).
#' ESPECIACION DE INVENTARIOS DE EMISIONES DE AEROSOLES Y COMPUESTOS ORGANICOS
#' VOLATILES PARA EL MODELO WRF-CHEM, APLICADO A LOS ESQUEMAS RADM-2,
#' CBM-Z Y MOZART-4.
#'
#'
#' @export
#' @examples
#' \dontrun{
#' # Not run
#' # Downloading EDGAR data ####
#' get_edgar(
#'   dataset = "v432_VOC_spec",
#'   destpath = "V50_432_AP/TOT/",
#'   sector = c("TOTALS"),
#'   type = "nc",
#'   year = 2012
#' )
#'
#' get_edgar(
#'   dataset = "v50_AP",
#'   destpath = "V50_432_AP/TOT",
#'   sector = c("TOTALS"),
#'   type = "nc",
#'   year = 2014
#' )
#'
#' get_edgar(
#'   dataset = "v432_VOC_spec",
#'   destpath = "V50_432_AP/TRO/",
#'   sector = c("TRO"),
#'   type = "nc",
#'   year = 2012, ask = F
#' )
#'
#' get_edgar(
#'   dataset = "v50_AP",
#'   destpath = "V50_432_AP/TRO",
#'   sector = c("TRO_RES", "TRO_noRES"),
#'   type = "nc",
#'   year = 2014
#' )
#'
#' totals <- list.files(
#'   path = "V50_432_AP/TOT/",
#'   full.names = TRUE,
#'   pattern = ".zip"
#' )
#' lapply(totals, unzip, exdir = "V50_432_AP/TOT//")
#'
#'
#' tros <- list.files(
#'   path = "V50_432_AP/TRO",
#'   full.names = TRUE,
#'   pattern = ".zip"
#' )
#' lapply(tros, unzip, exdir = "V50_432_AP/TRO/")
#' edgar_chem("V50_432_AP/TOT", "radm")
#' }
chem_edgar <- function(path,
                       chem,
                       merge = FALSE,
                       k = rep(1, 34),
                       verbose = TRUE) {
  dte <- sysdata$dte                         # nocov start
  emis_opt <- sysdata$emis_opt               # to heavy to check cod-cov

  if (length(k) < 34) stop("k must have 34 elements")

  # no covr, becuase it would take too much time
  ncs <- list.files(
    path = path,
    full.names = TRUE,
    pattern = ".nc"
  )
  # just checking units again
  a <- ncdf4::nc_open(ncs[1])
  cat("units: ", a$var[[1]]$units, "\n")
  ncdf4::nc_close(a)

  # unidades  kg m-2 s-1

  lncs <- list(
    voc1 = grep("voc1_", ncs, value = TRUE),
    voc2 = grep("voc2_", ncs, value = TRUE),
    voc3 = grep("voc3_", ncs, value = TRUE),
    voc4 = grep("voc4_", ncs, value = TRUE),
    voc5 = grep("voc5_", ncs, value = TRUE),
    voc6 = grep("voc6_", ncs, value = TRUE),
    voc7 = grep("voc7_", ncs, value = TRUE),
    voc8 = grep("voc8_", ncs, value = TRUE),
    voc9 = grep("voc9_", ncs, value = TRUE),
    voc10 = grep("voc10_", ncs, value = TRUE),
    voc11 = grep("voc11_", ncs, value = TRUE),
    voc12 = grep("voc12_", ncs, value = TRUE),
    voc13 = grep("voc13_", ncs, value = TRUE),
    voc14 = grep("voc14_", ncs, value = TRUE),
    voc15 = grep("voc15_", ncs, value = TRUE),
    voc16 = grep("voc16_", ncs, value = TRUE),
    voc17 = grep("voc17_", ncs, value = TRUE),
    voc18 = grep("voc18_", ncs, value = TRUE),
    voc19 = grep("voc19_", ncs, value = TRUE),
    voc20 = grep("voc20_", ncs, value = TRUE),
    vpc21 = grep("voc21_", ncs, value = TRUE),
    voc22 = grep("voc22_", ncs, value = TRUE),
    voc23 = grep("voc23_", ncs, value = TRUE),
    voc24 = grep("voc24_", ncs, value = TRUE),
    voc25 = grep("voc25_", ncs, value = TRUE),
    co = grep("CO_", ncs, value = TRUE),
    nox = grep("NOx_", ncs, value = TRUE),
    nmvoc = grep("NMVOC_", ncs, value = TRUE),
    so2 = grep("SO2_", ncs, value = TRUE),
    nh3 = grep("NH3_", ncs, value = TRUE),
    pm10 = grep("PM10_", ncs, value = TRUE),
    pm25 = grep("PM2.5_", ncs, value = TRUE),
    bc = grep("BC_", ncs, value = TRUE),
    oc = grep("_OC_", ncs, value = TRUE)
  )

  if(verbose) cat("Detecting the following files:\n")

  if(verbose) print(as.data.frame(cbind(lncs)))

  la <- unique(unlist(lapply(lncs, length)))

  if (length(la) > 1 & !merge) stop("There should be 1 NetCDF per pollutantt")

  fr <- function(x) raster::rotate(raster::raster(x)) * 1000 * 3600 * 1000 * 1000

  # EDGAR ####

  cat("EDGAR data: ")
  NCS_EDGAR <- data.frame(
    GEIA_id = c(
      paste0("voc", 1:25),
      c(
        "co", "nox", "nmvoc", "so2", "nh3",
        "pm10", "pm2.5", "bc", "oc"
      )
    )
  )


  dte <- cbind(dte, NCS_EDGAR)

  lp <- lapply(1:length(lncs), function(i) {
    cat(dte$GEIA_id[i], " ")

    if (merge & length(lncs[[i]]) > 1) {
      cat("merging ", lncs[[i]])
      cat("\n")
      a <- do.call("+", lapply(lncs[[i]], fr)) * k[i]
    } else {
      a <- fr(lncs[[i]]) * k[i]
    }
    a
  })
  # lp <- lapply(1:nrow(dte), function(i){
  #   cat(dte$GEIA_id[i], " ")
  #   raster::rotate(
  #     raster::raster(
  #       dte$ncs[i]))*1000*3600*1000*1000
  # })

  # ug/m2/s
  bp <- raster::stack(lp)
  bp@history <- list("units: g km-2 h-1")
  names(bp) <- NCS_EDGAR$GEIA_id

  if (missing(chem)) {
    mech <- c("edgar", "radm", "radmsorg", "cbmz_mosaic", "cptec", "ecb05_opt1")
    choice <- utils::menu(mech, title = "Choose:")
    chem <- mech[choice]
  }
  # chemical mechanisms
  if (chem == "edgar") {
    return(bp)
  }

  # radm ####
  if (chem == "radm") {
    E_ISO <- bp$voc10 / dte[dte$GEIA_id == "voc10", ]$g_mol # emis_opt$eradm[1]
    E_SO2 <- bp$so2 / (64 + 16 * 2) # emis_opt$eradm[2]
    E_NO <- bp$nox * 0.9 / (14 + 16) # emis_opt$eradm[3]
    E_CO <- bp$co / (12 + 18) # emis_opt$eradm[4]
    E_ETH <- bp$voc7 / dte[dte$GEIA_id == "voc7", ]$g_mol # emis_opt$eradm[5]
    E_HC3 <- bp$voc1 / dte[dte$GEIA_id == "voc1", ]$g_mol +
      bp$voc2 / dte[dte$GEIA_id == "voc2", ]$g_mol # emis_opt$eradm[6]
    E_HC5 <- bp$voc3 / dte[dte$GEIA_id == "voc3", ]$g_mol +
      bp$voc4 / dte[dte$GEIA_id == "voc4", ]$g_mol # emis_opt$eradm[7]
    E_HC8 <- bp$voc5 / dte[dte$GEIA_id == "voc5", ]$g_mol +
      bp$voc6 / dte[dte$GEIA_id == "voc6", ]$g_mol # emis_opt$eradm[8]
    E_XYL <- bp$voc13 / dte[dte$GEIA_id == "voc13", ]$g_mol +
      bp$voc14 / dte[dte$GEIA_id == "voc14", ]$g_mol # emis_opt$eradm[9]
    E_OL2 <- bp$voc9 / dte[dte$GEIA_id == "voc9", ]$g_mol +
      bp$voc11 / dte[dte$GEIA_id == "voc11", ]$g_mol # emis_opt$eradm[10]
    E_OLT <- bp$voc8 / dte[dte$GEIA_id == "voc8", ]$g_mol # emis_opt$eradm[11]
    E_OLI <- bp$voc12 / dte[dte$GEIA_id == "voc12", ]$g_mol # emis_opt$eradm[12]
    E_TOL <- bp$voc15 / dte[dte$GEIA_id == "voc15", ]$g_mol +
      bp$voc16 / dte[dte$GEIA_id == "voc16", ]$g_mol +
      bp$voc17 / dte[dte$GEIA_id == "voc17", ]$g_mol # emis_opt$eradm[13]
    E_CSL <- bp$voc18 / dte[dte$GEIA_id == "voc18", ]$g_mol +
      bp$voc19 / dte[dte$GEIA_id == "voc19", ]$g_mol # emis_opt$eradm[14]
    E_HCHO <- bp$voc21 / dte[dte$GEIA_id == "voc21", ]$g_mol # emis_opt$eradm[15]
    E_ALD <- bp$voc22 / dte[dte$GEIA_id == "voc22", ]$g_mol # emis_opt$eradm[16]
    E_KET <- bp$voc23 / dte[dte$GEIA_id == "voc23", ]$g_mol # emis_opt$eradm[17]
    E_ORA2 <- bp$voc24 / dte[dte$GEIA_id == "voc24", ]$g_mol # emis_opt$eradm[18]
    E_NH3 <- bp$nh3 / (14 + 3) # emis_opt$eradm[19]

    RADM <- raster::brick(list(
      E_ISO, E_SO2, E_NO, E_CO, E_ETH, E_HC3, E_HC5, E_HC8, E_XYL,
      E_OL2, E_OLT, E_OLI, E_TOL, E_CSL, E_HCHO, E_ALD,
      E_KET, E_ORA2, E_NH3
    ))
    names(RADM) <- emis_opt$eradm
    cat("units of gases: mol/km2/h\n")
    return(RADM)
    # radmsorg ####
  } else if (chem == "radmsorg") {
    E_ISO <- bp$voc10 / dte[dte$GEIA_id == "voc10", ]$g_mol # emis_opt$eradmsorg[1]
    E_SO2 <- bp$so2 / (64 + 16 * 2) # emis_opt$eradmsorg[2]
    E_NO <- bp$nox * 0.9 / (14 + 16) # emis_opt$eradmsorg[3]
    E_NO2 <- bp$nox * 0.1 / (14 + 16 * 1) # emis_opt$eradmsorg[4]
    E_CO <- bp$co / (12 + 18) # emis_opt$eradmsorg[5]
    E_CH4 <- bp$nmvoc * (1 / 100) / (12 + 4) # emis_opt$eradmsorg[6]
    E_ETH <- bp$voc7 / dte[dte$GEIA_id == "voc7", ]$g_mol # emis_opt$eradmsorg[7]
    E_HC3 <- bp$voc1 / dte[dte$GEIA_id == "voc1", ]$g_mol +
      bp$voc2 / dte[dte$GEIA_id == "voc2", ]$g_mol # emis_opt$eradmsorg[8]
    E_HC5 <- bp$voc3 / dte[dte$GEIA_id == "voc3", ]$g_mol +
      bp$voc4 / dte[dte$GEIA_id == "voc4", ]$g_mol # emis_opt$eradmsorg[9]
    dte[dte$MECH == "E_HC8", ]$GEIA_id
    E_HC8 <- bp$voc5 / dte[dte$GEIA_id == "voc5", ]$g_mol +
      bp$voc6 / dte[dte$GEIA_id == "voc6", ]$g_mol # emis_opt$eradmsorg[10]
    E_XYL <- bp$voc13 / dte[dte$GEIA_id == "voc13", ]$g_mol +
      bp$voc14 / dte[dte$GEIA_id == "voc14", ]$g_mol # emis_opt$eradmsorg[11]
    E_OL2 <- bp$voc9 / dte[dte$GEIA_id == "voc9", ]$g_mol +
      bp$voc11 / dte[dte$GEIA_id == "voc11", ]$g_mol # emis_opt$eradmsorg[12]
    E_OLT <- bp$voc8 / dte[dte$GEIA_id == "voc8", ]$g_mol # emis_opt$eradmsorg[13]
    E_OLI <- bp$voc12 / dte[dte$GEIA_id == "voc12", ]$g_mol # emis_opt$eradmsorg[14]
    E_TOL <- bp$voc15 / dte[dte$GEIA_id == "voc15", ]$g_mol +
      bp$voc16 / dte[dte$GEIA_id == "voc16", ]$g_mol +
      bp$voc17 / dte[dte$GEIA_id == "voc17", ]$g_mol # emis_opt$eradmsorg[15]
    E_CSL <- bp$voc18 / dte[dte$GEIA_id == "voc18", ]$g_mol +
      bp$voc19 / dte[dte$GEIA_id == "voc19", ]$g_mol # emis_opt$eradmsorg[16]
    E_HCHO <- bp$voc21 / dte[dte$GEIA_id == "voc21", ]$g_mol # emis_opt$eradmsorg[17]
    E_ALD <- bp$voc22 / dte[dte$GEIA_id == "voc22", ]$g_mol # emis_opt$eradmsorg[18]
    E_KET <- bp$voc23 / dte[dte$GEIA_id == "voc23", ]$g_mol # emis_opt$eradmsorg[19]
    E_ORA2 <- bp$voc24 / dte[dte$GEIA_id == "voc24", ]$g_mol # emis_opt$eradmsorg[20]
    E_NH3 <- bp$nh3 / (14 + 3) # emis_opt$eradmsorg[21]
    E_PM25I <- bp$pm2.5 * 0.2 / 3600 # emis_opt$eradmsorg[22]
    E_PM25J <- bp$pm2.5 * 0.8 / 3600 # emis_opt$eradmsorg[23]
    E_PM_10 <- bp$pm10 / 3600 # emis_opt$eradmsorg[24]
    E_ECI <- bp$bc * 0.2 / 3600 # emis_opt$eradmsorg[25]
    E_ECJ <- bp$bc * 0.8 / 3600 # emis_opt$eradmsorg[26]
    E_ORGI <- bp$oc * 0.2 / 3600 # emis_opt$eradmsorg[27]
    E_ORGJ <- bp$oc * 0.8 / 3600 # emis_opt$eradmsorg[28]
    E_SO4I <- bp$pm2.5 * 0.0077 / 3600 # emis_opt$eradmsorg[29]
    E_SO4J <- bp$pm2.5 * 0.0623 / 3600 # emis_opt$eradmsorg[30]
    E_NO3I <- bp$pm2.5 * 0.00247 / 3600 # emis_opt$eradmsorg[31]
    E_NO3J <- bp$pm2.5 * 0.01053 / 3600 # emis_opt$eradmsorg[32]
    E_NAAJ <- bp$pm2.5 * 0 / 3600 # emis_opt$eradmsorg[33]
    E_NAAI <- bp$pm2.5 * 0 / 3600 # emis_opt$eradmsorg[34]
    E_ORGI_A <- bp$pm2.5 * 0 / 3600 # emis_opt$eradmsorg[35]
    E_ORGJ_A <- bp$pm2.5 * 0 / 3600 # emis_opt$eradmsorg[36]
    E_ORGI_BB <- bp$pm2.5 * 0 / 3600 # emis_opt$eradmsorg[37]
    E_ORGJ_BB <- bp$pm2.5 * 0 / 3600 # emis_opt$eradmsorg[38]
    E_HCL <- bp$voc20 * 0.45 / dte[dte$GEIA_id == "voc20", ]$g_mol[1] # emis_opt$eradmsorg[39]
    E_CLI <- bp$voc20 * 0.02 / dte[dte$GEIA_id == "voc20", ]$g_mol[1] # emis_opt$eradmsorg[40]
    E_CLJ <- bp$voc20 * 0.08 / dte[dte$GEIA_id == "voc20", ]$g_mol[1] # emis_opt$eradmsorg[41]
    E_CH3CL <- bp$voc20 * 0.45 / dte[dte$GEIA_id == "voc20", ]$g_mol[1] # emis_opt$eradmsorg[42]

    RADMSORG <- raster::brick(list(
      E_ISO, E_SO2, E_NO, E_NO2, E_CO, E_CH4, E_ETH, E_HC3, E_HC5,
      E_HC8, E_XYL, E_OL2, E_OLT, E_OLI, E_TOL, E_CSL, E_HCHO, E_ALD,
      E_KET, E_ORA2, E_NH3, E_PM25I, E_PM25J, E_PM_10, E_ECI, E_ECJ,
      E_ORGI, E_ORGJ, E_SO4I, E_SO4J, E_NO3I, E_NO3J, E_NAAJ, E_NAAI,
      E_ORGI_A, E_ORGJ_A, E_ORGI_BB,
      E_ORGJ_BB, E_HCL, E_CLI, E_CLJ, E_CH3CL
    ))
    names(RADMSORG) <- emis_opt$eradmsorg
    cat("units of gases: mol/km2/h\n")
    cat("units of aerosols: ug/m2/s\n")
    return(RADMSORG)
    # cbmz_mosaic ####
  } else if (chem == "cbmz_mosaic") {
    E_ISO <- bp$voc10 / dte[dte$GEIA_id == "voc10", ]$g_mol # emis_opt$ecbmz_mosaic[1]

    E_SO2 <- bp$so2 / (64 + 16 * 2) # emis_opt$ecbmz_mosaic[2]

    E_NO <- bp$nox * 0.9 / (14 + 16) # emis_opt$ecbmz_mosaic[3]

    E_CO <- bp$co / (12 + 18) # emis_opt$ecbmz_mosaic[4]

    E_ETH <- bp$voc7 / dte[dte$GEIA_id == "voc7", ]$g_mol # emis_opt$ecbmz_mosaic[5]

    E_HC3 <- bp$voc1 / dte[dte$GEIA_id == "voc1", ]$g_mol +
      bp$voc2 / dte[dte$GEIA_id == "voc2", ]$g_mol # emis_opt$ecbmz_mosaic[6]

    E_HC5 <- bp$voc3 / dte[dte$GEIA_id == "voc3", ]$g_mol +
      bp$voc4 / dte[dte$GEIA_id == "voc4", ]$g_mol # emis_opt$ecbmz_mosaic[7]

    E_HC8 <- bp$voc5 / dte[dte$GEIA_id == "voc5", ]$g_mol +
      bp$voc6 / dte[dte$GEIA_id == "voc6", ]$g_mol # emis_opt$ecbmz_mosaic[8]

    E_XYL <- bp$voc13 / dte[dte$GEIA_id == "voc13", ]$g_mol +
      bp$voc14 / dte[dte$GEIA_id == "voc14", ]$g_mol # emis_opt$ecbmz_mosaic[9]

    E_OL2 <- bp$voc9 / dte[dte$GEIA_id == "voc9", ]$g_mol +
      bp$voc11 / dte[dte$GEIA_id == "voc11", ]$g_mol # emis_opt$ecbmz_mosaic[10]

    E_OLT <- bp$voc8 / dte[dte$GEIA_id == "voc8", ]$g_mol # emis_opt$ecbmz_mosaic[11]

    E_OLI <- bp$voc12 / dte[dte$GEIA_id == "voc12", ]$g_mol # emis_opt$ecbmz_mosaic[12]

    E_TOL <- bp$voc15 / dte[dte$GEIA_id == "voc15", ]$g_mol +
      bp$voc16 / dte[dte$GEIA_id == "voc16", ]$g_mol +
      bp$voc17 / dte[dte$GEIA_id == "voc17", ]$g_mol # emis_opt$ecbmz_mosaic[13]

    E_CSL <- bp$voc18 / dte[dte$GEIA_id == "voc18", ]$g_mol +
      bp$voc19 / dte[dte$GEIA_id == "voc19", ]$g_mol # emis_opt$ecbmz_mosaic[14]

    E_HCHO <- bp$voc21 / dte[dte$GEIA_id == "voc21", ]$g_mol # emis_opt$ecbmz_mosaic[15]

    E_ALD <- bp$voc22 / dte[dte$GEIA_id == "voc22", ]$g_mol # emis_opt$ecbmz_mosaic[16]

    E_KET <- bp$voc23 / dte[dte$GEIA_id == "voc23", ]$g_mol # emis_opt$ecbmz_mosaic[17]

    E_ORA2 <- bp$voc24 / dte[dte$GEIA_id == "voc24", ]$g_mol # emis_opt$ecbmz_mosaic[18]

    E_NH3 <- bp$nh3 / (14 + 3) # emis_opt$ecbmz_mosaic[19]

    E_NO2 <- bp$nox * 0.1 / (14 + 16 * 1) # emis_opt$ecbmz_mosaic[20]

    E_CH3OH <- bp$nmvoc * (1 / 100) / (12 + 3 + 16 + 1) # emis_opt$ecbmz_mosaic[21]

    E_C2H5OH <- bp$nmvoc * (5 / 100) / (12 * 2 + 5 + 16 + 1) # emis_opt$ecbmz_mosaic[22]

    E_PM25I <- bp$pm2.5 * 0.2 / 3600 # emis_opt$ecbmz_mosaic[23]

    E_PM25J <- bp$pm2.5 * 0.8 / 3600 # emis_opt$ecbmz_mosaic[24]

    E_ECI <- bp$bc * 0.2 / 3600 # emis_opt$ecbmz_mosaic[25]

    E_ECJ <- bp$bc * 0.8 / 3600 # emis_opt$ecbmz_mosaic[26]

    E_ORGI <- bp$oc * 0.2 / 3600 # emis_opt$ecbmz_mosaic[27]

    E_ORGJ <- bp$oc * 0.8 / 3600 # emis_opt$ecbmz_mosaic[28]

    E_SO4I <- bp$pm2.5 * 0.0077 / 3600 # emis_opt$ecbmz_mosaic[29]

    E_SO4J <- bp$pm2.5 * 0.0623 / 3600 # emis_opt$ecbmz_mosaic[30]

    E_NO3I <- bp$pm2.5 * 0.00247 / 3600 # emis_opt$ecbmz_mosaic[31]

    E_NO3J <- bp$pm2.5 * 0.01053 / 3600 # emis_opt$ecbmz_mosaic[32]

    E_SO4C <- bp$pm2.5 * 0 / 3600 # emis_opt$ecbmz_mosaic[33]

    E_NO3C <- bp$pm2.5 * 0 / 3600 # emis_opt$ecbmz_mosaic[34]

    E_ORGC <- bp$pm2.5 * 0 / 3600 # emis_opt$ecbmz_mosaic[35]

    E_ECC <- bp$pm2.5 * 0 / 3600 # emis_opt$ecbmz_mosaic[36]

    CBMZ_MOSAIC <- raster::brick(list(
      E_ISO, E_SO2, E_NO, E_CO, E_ETH, E_HC3, E_HC5, E_HC8, E_XYL, E_OL2,
      E_OLT, E_OLI, E_TOL, E_CSL, E_HCHO, E_ALD, E_KET, E_ORA2, E_NH3, E_NO2,
      E_CH3OH, E_C2H5OH, E_PM25I, E_PM25J, E_ECI, E_ECJ, E_ORGI, E_ORGJ, E_SO4I, E_SO4J,
      E_NO3I, E_NO3J, E_SO4C, E_NO3C, E_ORGC, E_ECC
    ))
    names(CBMZ_MOSAIC) <- emis_opt$ecbmz_mosaic
    cat("units of gases: mol/km2/h\n")
    cat("units of aerosols: ug/m2/s\n")
    return(CBMZ_MOSAIC)
    # cptec ####
  } else if (chem == "cptec") {
    E_ISO <- bp$voc10 / dte[dte$GEIA_id == "voc10", ]$g_mol # emis_opt$ecptec[1]

    E_SO2 <- bp$so2 / (64 + 16 * 2) # emis_opt$ecptec[2]

    E_NO <- bp$nox * 0.9 / (14 + 16) # emis_opt$ecptec[3]

    E_NO2 <- bp$nox * 0.1 / (14 + 16 * 2) # emis_opt$ecptec[4]

    E_CO <- bp$co / (12 + 18) # emis_opt$ecptec[5]

    E_ETH <- bp$voc7 / dte[dte$GEIA_id == "voc7", ]$g_mol # emis_opt$ecptec[6]

    E_HC3 <- bp$voc1 / dte[dte$GEIA_id == "voc1", ]$g_mol +
      bp$voc2 / dte[dte$GEIA_id == "voc2", ]$g_mol # emis_opt$ecptec[7]

    E_HC5 <- bp$voc3 / dte[dte$GEIA_id == "voc3", ]$g_mol +
      bp$voc4 / dte[dte$GEIA_id == "voc4", ]$g_mol # emis_opt$ecptec[8]

    E_HC8 <- bp$voc5 / dte[dte$GEIA_id == "voc5", ]$g_mol +
      bp$voc6 / dte[dte$GEIA_id == "voc6", ]$g_mol # emis_opt$ecptec[9]

    E_XYL <- bp$voc13 / dte[dte$GEIA_id == "voc13", ]$g_mol +
      bp$voc14 / dte[dte$GEIA_id == "voc14", ]$g_mol # emis_opt$ecptec[10]

    E_OL2 <- bp$voc9 / dte[dte$GEIA_id == "voc9", ]$g_mol +
      bp$voc11 / dte[dte$GEIA_id == "voc11", ]$g_mol # emis_opt$ecptec[11]

    E_OLT <- bp$voc8 / dte[dte$GEIA_id == "voc8", ]$g_mol # emis_opt$ecptec[12]

    E_OLI <- bp$voc12 / dte[dte$GEIA_id == "voc12", ]$g_mol # emis_opt$ecptec[13]

    E_TOL <- bp$voc15 / dte[dte$GEIA_id == "voc15", ]$g_mol +
      bp$voc16 / dte[dte$GEIA_id == "voc16", ]$g_mol +
      bp$voc17 / dte[dte$GEIA_id == "voc17", ]$g_mol # emis_opt$ecptec[14]

    E_CSL <- bp$voc18 / dte[dte$GEIA_id == "voc18", ]$g_mol +
      bp$voc19 / dte[dte$GEIA_id == "voc19", ]$g_mol # emis_opt$ecptec[15]

    E_HCHO <- bp$voc21 / dte[dte$GEIA_id == "voc21", ]$g_mol # emis_opt$ecptec[16]

    E_ALD <- bp$voc22 / dte[dte$GEIA_id == "voc22", ]$g_mol # emis_opt$ecptec[17]

    E_KET <- bp$voc23 / dte[dte$GEIA_id == "voc23", ]$g_mol # emis_opt$ecptec[18]

    E_ORA2 <- bp$voc24 / dte[dte$GEIA_id == "voc24", ]$g_mol # emis_opt$ecptec[19]

    E_NH3 <- bp$nh3 / (14 + 3) # emis_opt$ecptec[20]

    E_PM_25 <- bp$pm2.5 / 3600 # emis_opt$ecptec[21]

    E_PM_10 <- bp$pm10 / 3600 # emis_opt$ecptec[22]

    E_OC <- bp$oc / 3600 # emis_opt$ecptec[23]

    E_SULF <- bp$pm2.5 * 0.27 / 3600 # emis_opt$ecptec[24]

    E_BC <- bp$bc / 3600 # emis_opt$ecptec[25]

    CPTEC <- raster::brick(list(
      E_ISO, E_SO2, E_NO, E_NO2, E_CO, E_ETH, E_HC3, E_HC5, E_HC8, E_XYL, E_OL2,
      E_OLT, E_OLI, E_TOL, E_CSL, E_HCHO, E_ALD, E_KET, E_ORA2, E_NH3, E_PM_25, E_PM_10,
      E_OC, E_SULF, E_BC
    ))
    names(CPTEC) <- emis_opt$ecptec
    cat("units of gases: mol/km2/h\n")
    cat("units of aerosols: ug/m2/s\n")
    return(CPTEC)
    # ecb05_opt1 ####
  } else if (chem == "ecb05_opt1") {
    E_NO2 <- bp$nox * 0.1 / (14 + 16 * 2) # emis_opt$ecb05_opt1[1]

    E_XYL <- bp$voc13 / dte[dte$GEIA_id == "voc13", ]$g_mol +
      bp$voc14 / dte[dte$GEIA_id == "voc14", ]$g_mol # emis_opt$ecb05_opt1[2]

    E_TOL <- bp$voc15 / dte[dte$GEIA_id == "voc15", ]$g_mol +
      bp$voc16 / dte[dte$GEIA_id == "voc16", ]$g_mol +
      bp$voc17 / dte[dte$GEIA_id == "voc17", ]$g_mol # emis_opt$ecb05_opt1[3]

    E_TERP <- bp$voc10 / dte[dte$GEIA_id == "voc11", ]$g_mol # emis_opt$ecb05_opt1[4]

    E_SO2 <- bp$so2 / (64 + 16 * 2) # emis_opt$ecb05_opt1[5]

    E_ORA2 <- bp$voc24 / dte[dte$GEIA_id == "voc24", ]$g_mol # emis_opt$ecb05_opt1[6]

    E_OLT <- bp$voc8 / dte[dte$GEIA_id == "voc8", ]$g_mol # emis_opt$ecb05_opt1[7]

    E_OLI <- bp$voc12 / dte[dte$GEIA_id == "voc12", ]$g_mol # emis_opt$ecb05_opt1[8]

    E_OL2 <- bp$voc9 / dte[dte$GEIA_id == "voc9", ]$g_mol +
      bp$voc11 / dte[dte$GEIA_id == "voc11", ]$g_mol # emis_opt$ecb05_opt1[9]

    E_NO <- bp$nox * 0.9 / (14 + 16) # emis_opt$ecb05_opt1[10]

    E_NH3 <- bp$nh3 / (14 + 3) # emis_opt$ecb05_opt1[11]

    E_ISO <- bp$voc10 / dte[dte$GEIA_id == "voc10", ]$g_mol # emis_opt$ecb05_opt1[12]

    E_HCL <- bp$voc20 * 1 / dte[dte$GEIA_id == "voc20", ]$g_mol[1] # emis_opt$ecb05_opt1[13]

    E_HCHO <- bp$voc21 / dte[dte$GEIA_id == "voc21", ]$g_mol # emis_opt$ecb05_opt1[14]

    E_ETH <- bp$voc7 / dte[dte$GEIA_id == "voc7", ]$g_mol # emis_opt$ecb05_opt1[15]

    E_CSL <- bp$voc18 / dte[dte$GEIA_id == "voc18", ]$g_mol +
      bp$voc19 / dte[dte$GEIA_id == "voc19", ]$g_mol # emis_opt$ecb05_opt1[16]

    E_CO <- bp$co / (12 + 16) # emis_opt$ecb05_opt1[17]

    E_CH3OH <- bp$nmvoc * (1 / 100) / (12 + 3 + 16 + 1) # emis_opt$ecb05_opt1[18]

    E_C2H5OH <- bp$nmvoc * (5 / 100) / (12 * 2 + 5 + 16 + 1) # emis_opt$ecb05_opt1[19]

    # Acetaldehyde?
    E_ALD <- bp$voc22 / dte[dte$GEIA_id == "voc22", ]$g_mol * 0.9 # emis_opt$ecb05_opt1[20]

    # C3+ Aldehydes
    E_ALDX <- bp$voc22 / dte[dte$GEIA_id == "voc22", ]$g_mol * 0.1 # emis_opt$ecb05_opt1[21]

    E_HC3 <- bp$voc1 / dte[dte$GEIA_id == "voc1", ]$g_mol +
      bp$voc2 / dte[dte$GEIA_id == "voc2", ]$g_mol # emis_opt$ecb05_opt1[22]

    E_HC5 <- bp$voc3 / dte[dte$GEIA_id == "voc3", ]$g_mol +
      bp$voc4 / dte[dte$GEIA_id == "voc4", ]$g_mol # emis_opt$ecb05_opt1[23]

    E_HC8 <- bp$voc5 / dte[dte$GEIA_id == "voc5", ]$g_mol +
      bp$voc6 / dte[dte$GEIA_id == "voc6", ]$g_mol # emis_opt$ecb05_opt1[24]

    E_KET <- bp$voc23 / dte[dte$GEIA_id == "voc23", ]$g_mol # emis_opt$ecb05_opt1[25]

    E_PM25I <- bp$pm2.5 * 0.2 / 3600 # not read by mechanism   # emis_opt$ecb05_opt1[26]

    E_PM25J <- bp$pm2.5 * (0.8 + 0.2) / 3600 # emis_opt$ecb05_opt1[27]

    E_ECI <- bp$bc * 0.2 / 3600 # not read by mechanism   # emis_opt$ecb05_opt1[28]

    E_ECJ <- bp$bc * (0.8 + 0.2) / 3600 # emis_opt$ecb05_opt1[29]

    E_ORGI <- bp$oc * 0.2 / 3600 # not read by mechanism   # emis_opt$ecb05_opt1[30]

    E_ORGJ <- bp$oc * (0.8 + 0.2) / 3600 # emis_opt$ecb05_opt1[31]

    E_SO4I <- bp$pm2.5 * 0.0077 / 3600 # not read by mechanism  # emis_opt$ecb05_opt1[32]

    E_SO4J <- bp$pm2.5 * (0.0623 + 0.0077) / 3600 # emis_opt$ecb05_opt1[33]

    E_NO3I <- bp$pm2.5 * 0.00247 / 3600 # not read by mechanism  # emis_opt$ecb05_opt1[34]

    E_NO3J <- bp$pm2.5 * (0.01053 + 0.00247) / 3600 # emis_opt$ecb05_opt1[35]

    E_SO4C <- bp$pm2.5 * 0 / 3600 # emis_opt$ecb05_opt1[36]

    E_NO3C <- bp$pm2.5 * 0 / 3600 # emis_opt$ecb05_opt1[37]

    E_ORGC <- bp$pm2.5 * 0 / 3600 # emis_opt$ecb05_opt1[38]

    E_ECC <- bp$pm2.5 * 0 / 3600 # emis_opt$ecb05_opt1[39]

    E_PM10 <- bp$pm10 * 1 / 3600 # emis_opt$ecb05_opt1[40]

    CB05 <- raster::brick(list(
      E_NO2, E_XYL, E_TOL, E_TERP, E_SO2, E_ORA2, E_OLT, E_OLI, E_OL2, E_NO,
      E_NH3, E_ISO, E_HCL, E_HCHO, E_ETH, E_CSL, E_CO, E_CH3OH, E_C2H5OH, E_ALD,
      E_ALDX, E_HC3, E_HC5, E_HC8, E_KET, E_PM25I, E_PM25J, E_ECI, E_ECJ, E_ORGI,
      E_ORGJ, E_SO4I, E_SO4J, E_NO3I, E_NO3J, E_SO4C, E_NO3C, E_ORGC, E_ECC, E_PM10
    ))
    names(CB05) <- emis_opt$ecb05_opt1
    cat("units of gases: mol/km2/h\n")
    cat("units of aerosols: ug/m2/s\n")
    return(CB05)
    # neu_ecb05 ####
  } else if (chem %in% c("neu_ecb05", "neu_cb05")) {
    cb05 <- sysdata$edgar_cb05

    # emis_opt$ecb05_opt2[1]
    E_ACET <- 0 * bp$voc9

    # emis_opt$ecb05_opt2[2]
    E_PAR <- cb05[cb05$GEIA_id == "voc3", ]$M * bp$voc3  / cb05[cb05$GEIA_id == "voc3", ]$g_mol +
             cb05[cb05$GEIA_id == "voc4", ]$M * bp$voc4  / cb05[cb05$GEIA_id == "voc4", ]$g_mol +
             cb05[cb05$GEIA_id == "voc5", ]$M * bp$voc5  / cb05[cb05$GEIA_id == "voc5", ]$g_mol +
             cb05[cb05$GEIA_id == "voc6", ]$M * bp$voc6  / cb05[cb05$GEIA_id == "voc6", ]$g_mol +
             cb05[cb05$GEIA_id == "voc9", ]$M * bp$voc9  / cb05[cb05$GEIA_id == "voc9", ]$g_mol +
             cb05[cb05$GEIA_id == "voc13",]$M * bp$voc13 / cb05[cb05$GEIA_id == "voc13",]$g_mol +
             cb05[cb05$GEIA_id == "voc18",]$M * bp$voc18 / cb05[cb05$GEIA_id == "voc18",]$g_mol +
             cb05[cb05$GEIA_id == "voc19",]$M * bp$voc19 / cb05[cb05$GEIA_id == "voc19",]$g_mol

    # emis_opt$ecb05_opt2[3]
    E_ALK3 <- bp$voc9 * 0

    # emis_opt$ecb05_opt2[4]
    E_ALK4 <- bp$voc9 * 0

    # emis_opt$ecb05_opt2[5]
    E_ALK5 <- bp$voc9 * 0

    # emis_opt$ecb05_opt2[6]
    E_TOL <- bp$voc14 / cb05[cb05$group == "E_TOL", ]$g_mol

    # emis_opt$ecb05_opt2[7]
    E_XYL <- bp$voc15 / cb05[cb05$GEIA_id == "voc15", ]$g_mol +
             bp$voc16 / cb05[cb05$GEIA_id == "voc16", ]$g_mol +
             bp$voc17 / cb05[cb05$GEIA_id == "voc17", ]$g_mol

    # emis_opt$ecb05_opt2[8]
    E_BALD <- bp$voc13 * 0

    # emis_opt$ecb05_opt2[9]
    E_ALD2 <- bp$voc22 / cb05[cb05$group == "E_ALD2", ]$g_mol

    # emis_opt$ecb05_opt2[10]
    E_CCOOH <- bp$voc22 * 0

    # emis_opt$ecb05_opt2[11]
    E_CO <- bp$co / (12 + 16)

    # emis_opt$ecb05_opt2[12]
    E_CRES <- 0 * bp$voc18

    # emis_opt$ecb05_opt2[13]
    E_ETH <- bp$voc7 / cb05[cb05$group == "E_ETH", ]$g_mol

    # emis_opt$ecb05_opt2[14]
    E_ETHA <- bp$voc2 / cb05[cb05$group == "E_ETHA", ]$g_mol

    # emis_opt$ecb05_opt2[15]
    E_GLY <- bp$voc7 * 0

    # emis_opt$ecb05_opt2[16]
    E_FORM <- bp$voc21 / cb05[cb05$group == "E_FORM", ]$g_mol

    # emis_opt$ecb05_opt2[17]
    E_HCOOH <- bp$voc21 * 0

    # emis_opt$ecb05_opt2[18]
    E_IPROD <- bp$voc21 * 0

    # emis_opt$ecb05_opt2[19]
    E_ISOP <- bp$voc10 / cb05[cb05$group == "E_ISOP", ]$g_mol

    # emis_opt$ecb05_opt2[20]
    E_MACR <- bp$voc10 * 0

    # emis_opt$ecb05_opt2[21]
    E_MEK <- bp$voc10 * 0

    # emis_opt$ecb05_opt2[22]
    E_MEOH <- cb05[cb05$group == "E_MEOH", ]$M * bp$voc1 / cb05[cb05$group == "E_MEOH", ]$g_mol

    # emis_opt$ecb05_opt2[23]
    E_MEO2 <- bp$voc10 * 0

    # emis_opt$ecb05_opt2[24]
    E_ETOH <- cb05[cb05$group == "E_ETOH", ]$M * bp$voc1 / cb05[cb05$group == "E_ETOH", ]$g_mol

    # emis_opt$ecb05_opt2[25]
    E_MGLY <- bp$nmvoc * 0

    # emis_opt$ecb05_opt2[26]
    E_NH3 <- bp$nh3 / (14 + 3)

    # emis_opt$ecb05_opt2[27]
    E_HCL <- bp$voc20 * 0
    # Chlorinated hydrocarbons	CH3Cl

    # emis_opt$ecb05_opt2[28]
    E_NO <- bp$nox * cb05[cb05$group == "E_NO", ]$M / (14 + 16)

    # emis_opt$ecb05_opt2[29]
    E_NO2 <- bp$nox * cb05[cb05$group == "E_NO2", ]$M / (14 + 16 * 2)

    # emis_opt$ecb05_opt2[30]
    E_IOLE <- cb05[cb05$group == "E_IOLE", ]$M * bp$voc12 / cb05[cb05$group == "E_IOLE", ]$g_mol

    # emis_opt$ecb05_opt2[31]
    E_OLE <- bp$voc8 / cb05[cb05$group == "E_OLE", ]$g_mol

    # emis_opt$ecb05_opt2[32]
    E_PHEN <- bp$voc8 * 0

    # emis_opt$ecb05_opt2[33]
    E_PROD2 <- bp$voc8 * 0

    # emis_opt$ecb05_opt2[34]
    E_ALDX <- cb05[cb05$group == "E_ALDX", ]$M * bp$voc12 * bp$voc12 / cb05[cb05$group == "E_ALDX", ]$g_mol

    # emis_opt$ecb05_opt2[35]
    E_SO2 <- bp$so2 / (64 + 16 * 2)

    # emis_opt$ecb05_opt2[36]
    E_PSULF <- bp$pm2.5 * 0.27 / 3600

    # emis_opt$ecb05_opt2[37]
    E_TERP <- bp$voc11 / cb05[cb05$group == "E_TERP", ]$g_mol

    # emis_opt$ecb05_opt2[38]
    E_PM25I <- bp$pm2.5 * 0.0 # not read by mechanism

    # emis_opt$ecb05_opt2[39]
    E_PM25J <- bp$pm2.5 * (0.8 + 0.2) / 3600

    # emis_opt$ecb05_opt2[40]
    E_ECI <- bp$bc * 0.0  # not read by mechanism

    # emis_opt$ecb05_opt2[41]
    E_ECJ <- bp$bc * (0.8 + 0.2) / 3600

    # emis_opt$ecb05_opt2[42]
    E_ORGI <- bp$oc * 0.0 # not read by mechanism

    # emis_opt$ecb05_opt2[43]
    E_ORGJ <- bp$oc * (0.8 + 0.2) / 3600

    # emis_opt$ecb05_opt2[44]
    E_SO4I <- bp$pm2.5 * 0.0 # not read by mechanism

    # emis_opt$ecb05_opt2[45]
    E_SO4J <- bp$pm2.5 * (0.0623 + 0.0077) / 3600

    # emis_opt$ecb05_opt2[46]
    E_NO3I <- bp$pm2.5 * 0.0 # not read by mechanism

    # emis_opt$ecb05_opt2[47]
    E_NO3J <- bp$pm2.5 * (0.01053 + 0.00247) / 3600

    # emis_opt$ecb05_opt2[48]
    E_SO4C <- bp$pm2.5 * 0

    # emis_opt$ecb05_opt2[49]
    E_NO3C <- bp$pm2.5 * 0

    # emis_opt$ecb05_opt2[50]
    E_ORGC <- bp$pm2.5 * 0

    # emis_opt$ecb05_opt2[51]
    E_ECC <- bp$pm2.5 * 0

    # emis_opt$ecb05_opt2[52]
    E_PM10 <- (bp$pm10 - bp$pm2.5) * 1 / 3600


    emis_opt$ecb05_opt2

    CB05 <- raster::brick(list(
      E_ACET, E_PAR, E_ALK3, E_ALK4, E_ALK5, E_TOL,
      E_XYL, E_BALD, E_ALD2, E_CCOOH, E_CO,
      E_CRES, E_ETH, E_ETHA, E_GLY, E_FORM,
      E_HCOOH, E_IPROD, E_ISOP, E_MACR, E_MEK, E_MEOH,
      E_MEO2, E_ETOH, E_MGLY, E_NH3, E_HCL,
      E_NO, E_NO2, E_IOLE, E_OLE, E_PHEN, E_PROD2,
      E_ALDX, E_SO2, E_PSULF, E_TERP, E_PM25I,
      E_PM25J, E_ECI, E_ECJ, E_ORGI, E_ORGJ,
      E_SO4I, E_SO4J, E_NO3I, E_NO3J,
      E_SO4C, E_NO3C, E_ORGC, E_ECC, E_PM10
    ))
    names(CB05) <- emis_opt$ecb05_opt2
    cat("units of gases: mol/km2/h\n")
    cat("units of aerosols: ug/m2/s\n")
    return(CB05)
  } else if (chem == "utfpr_cbmz") {
    # emis_opt$ecbmz_mosaic[1]
    E_ISO <- bp$voc10 / dte[dte$GEIA_id == "voc10", ]$g_mol +
      bp$voc11 / dte[dte$GEIA_id == "voc11", ]$g_mol

    # emis_opt$ecbmz_mosaic[2]
    E_SO2 <- bp$so2 / (64 + 16 * 2)

    # emis_opt$ecbmz_mosaic[3]
    E_NO <- bp$nox * 0.9 / (14 + 16)

    # emis_opt$ecbmz_mosaic[4]
    E_CO <- bp$co / (12 + 18)

    # emis_opt$ecbmz_mosaic[5]
    E_ETH <- bp$voc2 / dte[dte$GEIA_id == "voc2", ]$g_mol

    # emis_opt$ecbmz_mosaic[6]
    E_HC3 <- bp$voc3 / dte[dte$GEIA_id == "voc3", ]$g_mol +
      bp$voc9 / dte[dte$GEIA_id == "voc9", ]$g_mol

    # emis_opt$ecbmz_mosaic[7]
    E_HC5 <- bp$voc5 / dte[dte$GEIA_id == "voc5", ]$g_mol +
      bp$voc4 / dte[dte$GEIA_id == "voc4", ]$g_mol

    # emis_opt$ecbmz_mosaic[8]
    E_HC8 <- bp$voc6 / dte[dte$GEIA_id == "voc6", ]$g_mol

    # emis_opt$ecbmz_mosaic[9]
    E_XYL <- bp$voc16 / dte[dte$GEIA_id == "voc16", ]$g_mol +
      bp$voc15 / dte[dte$GEIA_id == "voc15", ]$g_mol +
      bp$voc17 / dte[dte$GEIA_id == "voc17", ]$g_mol

    # emis_opt$ecbmz_mosaic[10]
    E_OL2 <- bp$voc7 / dte[dte$GEIA_id == "voc7", ]$g_mol

    # emis_opt$ecbmz_mosaic[11]
    E_OLT <- bp$voc8 / dte[dte$GEIA_id == "voc8", ]$g_mol

    # emis_opt$ecbmz_mosaic[12]
    E_OLI <- bp$voc12 / dte[dte$GEIA_id == "voc12", ]$g_mol

    # emis_opt$ecbmz_mosaic[13]
    E_TOL <- bp$voc13 / dte[dte$GEIA_id == "voc13", ]$g_mol +
      bp$voc14 / dte[dte$GEIA_id == "voc14", ]$g_mol

    # emis_opt$ecbmz_mosaic[14]
    E_CSL <- bp$voc18 * 0

    # emis_opt$ecbmz_mosaic[15]
    E_HCHO <- bp$voc21 / dte[dte$GEIA_id == "voc21", ]$g_mol

    # emis_opt$ecbmz_mosaic[16]
    E_ALD <- bp$voc22 / dte[dte$GEIA_id == "voc22", ]$g_mol

    # emis_opt$ecbmz_mosaic[17]
    E_KET <- bp$voc23 / dte[dte$GEIA_id == "voc23", ]$g_mol +
      bp$voc18 / dte[dte$GEIA_id == "voc18", ]$g_mol +
      bp$voc19 / dte[dte$GEIA_id == "voc19", ]$g_mol

    # emis_opt$ecbmz_mosaic[18]
    E_ORA2 <- bp$voc24 / dte[dte$GEIA_id == "voc24", ]$g_mol

    # emis_opt$ecbmz_mosaic[19]
    E_NH3 <- bp$nh3 / (14 + 3)

    # emis_opt$ecbmz_mosaic[20]
    E_NO2 <- bp$nox * 0.1 / (14 + 16 * 1)

    # emis_opt$ecbmz_mosaic[21]
    E_CH3OH <- bp$voc1 * 0.2 / (12 + 3 + 16 + 1)

    # emis_opt$ecbmz_mosaic[22]
    E_C2H5OH <- bp$voc1 * 0.8 / (12 * 2 + 5 + 16 + 1)

    # emis_opt$ecbmz_mosaic[23]
    E_PM25I <- bp$pm2.5 * 0.25 / 3600

    # emis_opt$ecbmz_mosaic[24]
    E_PM25J <- bp$pm2.5 * 0.75 / 3600

    # emis_opt$ecbmz_mosaic[25]
    E_ECI <- bp$bc * 0.94 / 3600

    # emis_opt$ecbmz_mosaic[26]
    E_ECJ <- bp$bc * 0.06 / 3600

    # emis_opt$ecbmz_mosaic[27]
    E_ORGI <- bp$oc * 0.190 / 3600

    # emis_opt$ecbmz_mosaic[28]
    E_ORGJ <- bp$oc * 0.81 / 3600

    # emis_opt$ecbmz_mosaic[29]
    E_SO4I <- bp$pm2.5 * 0.070 * 0.136 / 3600

    # emis_opt$ecbmz_mosaic[30]
    E_SO4J <- bp$pm2.5 * 0.070 * 0.864 / 3600

    # emis_opt$ecbmz_mosaic[31]
    E_NO3I <- bp$pm2.5 * 0.016 * 0.230 / 3600

    # emis_opt$ecbmz_mosaic[32]
    E_NO3J <- bp$pm2.5 * 0.016 * 0.770 / 3600

    # emis_opt$ecbmz_mosaic[33]
    E_SO4C <- bp$pm2.5 * 0 / 3600

    # emis_opt$ecbmz_mosaic[34]
    E_NO3C <- bp$pm2.5 * 0 / 3600

    # emis_opt$ecbmz_mosaic[35]
    E_ORGC <- bp$pm2.5 * 0 / 3600

    # emis_opt$ecbmz_mosaic[36]
    E_ECC <- bp$pm2.5 * 0 / 3600

    CBMZ_MOSAIC <- raster::brick(list(
      E_ISO, E_SO2, E_NO, E_CO, E_ETH, E_HC3, E_HC5, E_HC8, E_XYL, E_OL2,
      E_OLT, E_OLI, E_TOL, E_CSL, E_HCHO, E_ALD, E_KET, E_ORA2, E_NH3, E_NO2,
      E_CH3OH, E_C2H5OH, E_PM25I, E_PM25J, E_ECI, E_ECJ, E_ORGI, E_ORGJ, E_SO4I, E_SO4J,
      E_NO3I, E_NO3J, E_SO4C, E_NO3C, E_ORGC, E_ECC
    ))
    names(CBMZ_MOSAIC) <- emis_opt$ecbmz_mosaic
    cat("units of gases: mol/km2/h\n")
    cat("units of aerosols: ug/m2/s\n")
    return(CBMZ_MOSAIC)
  } else {
    stop("other mechanisms not implemented")
  }                                               # nocov end
}