R/calcEmiMac.R

In pik-piam/mrcommons: MadRat commons Input Data Library

#' Calculate baseline emissions for maccs
#'
#' Provides REMIND data for baseline emissions for maccs.
#'
#'
#' @return REMIND data for baseline emissions for maccs and corresonding
#' weights (NULL) as a list of two MAgPIE objects
#' @author Lavinia Baumstark
#' @seealso [calcOutput()], [readSource()]
#' @examples
#' \dontrun{
#' calcOutput("calcEmiMac")
#' }
#' @importFrom magclass getNames<- getYears<-

calcEmiMac <- function() {

  # emissions for the calculation of econometric paramter p1
  co2       <- readSource("EDGAR", subtype = "co2") * 12 / 44 * 1e-6

  n2o       <- readSource("EDGAR", subtype = "n2owaste")  * 28 / 44 * 1e-3
  ch4       <- readSource("EDGAR", subtype = "ch4waste") * 1e-3

  co2cement <- dimSums(co2[, , c("2A1", "2A2")], dim = 3)
  getNames(co2cement) <- "co2cement_process"
  co2luc <- dimSums(co2[, , c("5A", "5D", "5F2")], dim = 3)
  getNames(co2luc) <- "co2luc"
  n2owaste <- dimSums(n2o[, , c("6A", "6B", "6C", "6D")], dim = 3)
  getNames(n2owaste) <- "n2owaste"
  n2otrans <- dimSums(n2o[, , c("1A3a", "1A3b", "1A3c", "1A3d", "1A3e")], dim = 3)
  getNames(n2otrans) <- "n2otrans"
  n2oacid <- dimSums(n2o[, , c("2B")], dim = 3)
  getNames(n2oacid) <- "n2oacid"
  ch4wsts  <- dimSums(ch4[, , c("6B")], dim = 3)
  getNames(ch4wsts) <- "ch4wsts"
  ch4wstl  <- dimSums(ch4[, , c("6A", "6C", "6D")], dim = 3)
  getNames(ch4wstl) <- "ch4wstl"

  # overwritting european countries with eurostat data
  euCountries <- c("ALA", "AUT", "BEL", "BGR", "HRV", "CYP", "CZE", "DNK", "EST", "FRO",
  "FIN", "FRA", "DEU", "GIB", "GRC", "GGY", "HUN", "IRL", "IMN", "ITA", "JEY", "LVA", "LTU",
  "LUX", "MLT", "NLD", "POL", "PRT", "ROU", "SVK", "SVN", "ESP", "SWE", "GBR")
  baselineEurostat <- calcOutput("HistEmissions", subtype = "MAC", aggregate = FALSE)
  n2owaste[euCountries, 2005, "n2owaste"] <- baselineEurostat[euCountries, 2005, "n2owaste"]
  n2otrans[euCountries, 2005, "n2otrans"] <- baselineEurostat[euCountries, 2005, "n2otrans"]
  n2oacid[euCountries, 2005, "n2oacid"] <- dimSums(baselineEurostat[euCountries, 2005, c("n2oadac", "n2onitac")])
  ch4wsts[euCountries, 2005, "ch4wsts"] <- baselineEurostat[euCountries, 2005, "ch4wsts"]
  ch4wstl[euCountries, 2005, "ch4wstl"] <- baselineEurostat[euCountries, 2005, "ch4wstl"]

  # combine all parameters
  x <- mbind(co2cement, co2luc, n2owaste, n2otrans, n2oacid, ch4wsts, ch4wstl)
  getYears(x) <- NULL
  return(list(x = x,
              weight = NULL,
              unit = "GtC, MtCH4, MtN2O",
              description = "emissions in 2005",
              note = c("used to calculate econometric emission parameter p1")))
}