R/calcEconometricEmiParameter.R

In pik-piam/mrremind: MadRat REMIND Input Data Package

#' Calculate baseline emissions of waste
#'
#' Provides REMIND data for CO2 parameters to calculate baseline emissions of
#' waste from population and investment.
#'
#'
#' @return REMIND data forCO2 parameters to calculate baseline emissions of
#' waste from population and investment and corresonding weights (population)
#' as a list of two MAgPIE objects
#' @author Lavinia Baumstark
#' @seealso \code{\link{calcOutput}}, \code{\link{readSource}}
#' @examples
#' \dontrun{
#' calcOutput("calcEconometricEmiParameter")
#' }
#' @importFrom utils read.csv
#' @importFrom magclass getNames<- getYears<-


calcEconometricEmiParameter <- function() {

  # calculate parameter p3
  edgar <- readSource("EDGAR", subtype = "co2") * 12 / 44 * 1e-6
  pop   <- calcOutput("Population", aggregate = FALSE)[, 2005, "pop_SSP2"] / 1000

  p3Country <- dimSums(edgar[, , c("2A1", "2A2")], dim = 3) / pop
  mapping    <- toolGetMapping(type = "regional", name = "regionmappingEconometricEmi.csv", where = "mappingfolder")
  p3         <- toolAggregate(p3Country[unique(mapping$RepCountry), , ], mapping, weight = NULL)
  p3[is.na(p3)] <- 0
  getNames(p3) <- "co2cement_process.p3"

  # calculate parameter p4
  p4 <- readSource("REMIND_11Regi", subtype = "p4")
  getNames(p4) <- "co2cement_process.p4"
  getYears(p4) <- getYears(p3)

  # combine all parameters
  x <- mbind(p3, p4)
  getYears(x) <- NULL
  return(list(x = x,
              weight = pop,
              unit = "p3 - GtC/Cap, p4 - $US2005/Cap",
              description = "CO2 parameters to calculate baseline emissions of waste from population and investment",
              note = c("more to be added by Jessi")))
}