R/readREMIND_11Regi.R

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

#' Read REMIND region dependent data
#'
#' Read-in an csv files that contains regional data
#'
#' @param subtype Name of the regional data, e.g. "p4", "biomass", "ch4waste",
#' "tradecost", "pe2se", "xpres_tax", "deltacapoffset", "capacityFactorGlobal",
#' "capacityFactorRules", "residuesShare", "taxConvergence", "maxFeSubsidy",
#' "maxPeSubsidy", "propFeSubsidy", "fossilExtractionCoeff", "uraniumExtractionCoeff",
#' "RLDCCoefficientsLoB", "RLDCCoefficientsPeak", "earlyRetirementAdjFactor"
#' @return magpie object of region dependent data
#' @author original: not defined, capacity factor, tax, fossil and RLDC changes: Renato Rodrigues
#' @seealso \code{\link{readSource}}
#' @examples
#' \dontrun{
#' a <- readSource(type = "REMIND_11Regi", subtype = "capacityFactorGlobal")
#' }
readREMIND_11Regi <- function(subtype) {

  if (subtype == "p4") {
    x <- read.csv("EconometricEmissionParameter_p4.csv", sep = ";", row.names = 1)
    x <- as.magpie(x)
  } else if (subtype == "biomass") {
    x <- read_excel("biomass.xlsx")
    x <- as.magpie(x, datacol = 4)
  } else if (subtype == "ch4waste") {
    x <- read.csv("emimac0_ch4waste.csv", sep = ";", row.names = 1)
    x <- as.magpie(x)
  } else if (subtype == "tradecost") {
    x <- read.csv("LueckenDiss_TradeCost.csv", sep = ";", row.names = 1)
    x <- as.magpie(x)
  } else if (subtype == "pe2se") {
    x <- read.csv("tax_pe2se_sub.csv", sep = ";")
    x <- as.magpie(x, datacol = 2)
  } else if (subtype == "xpres_tax") {
    x <- read.csv("p21_tau_xpres_tax.csv", sep = ";")
    x <- as.magpie(x, datacol = 2)
    getYears(x) <- "y2005"
    getNames(x) <- "peoil"
  } else if (subtype == "deltacapoffset") {
    x <- read.csv("p_adj_deltacapoffset.csv", sep = ";")
    x <- as.magpie(x, datacol = 2)
  } else if (subtype == "capacityFactorGlobal") {
    x <- read.csv("f_cf-global.csv", sep = ";")
    x[x$Tech == "rockgrind", "Tech"] <- "weathering"
    x <- as.magpie(x, datacol = 2)
  } else if (subtype == "capacityFactorRules") {
    x <- read.csv("f_cf-rules.csv", sep = ";")
    x <- as.magpie(x, datacol = 4)
  } else if (subtype == "storageFactor") {
    x <- read.csv("storageFactor.csv", sep = ";")
    x <- as.magpie(x, datacol = 2)
  } else if (subtype == "residuesShare") {
    x <- read.csv("residuesShare.csv", row.names = 1)
    x <- as.magpie(x, datacol = 4)
  } else if (subtype == "shareIndFE") {
    x <- read.csv("shareIndustyFE.csv", sep = ";", skip = 3)
    x <- as.magpie(x, datacol = 3)
  } else if (subtype == "taxConvergence") {
    x <- read.csv("tax_convergence.csv", sep = ";")
    x <- as.magpie(x, datacol = 4)
  } else if (subtype == "maxFeSubsidy") {
    x <- read.csv("max_FE_subsidy.csv", sep = ";")
    x <- as.magpie(x, datacol = 4)
  } else if (subtype == "maxPeSubsidy") {
    x <- read.csv("max_PE_subsidy.csv", sep = ";")
    x <- as.magpie(x, datacol = 4)
  } else if (subtype == "propFeSubsidy") {
    x <- read.csv("prop_FE_subsidy.csv", sep = ";")
    x <- as.magpie(x, datacol = 4)
  } else if (subtype == "gridFactor") {
    x <- read.csv("homogenous_regions_for grids.csv", sep = ";")
    x$X <- NULL
    x <- as.magpie(x, datacol = 2)
  } else if (subtype == "AP_starting_values") {
    x <- read.csv("f11_emiAPexsolve.cs4r", sep = ",", skip = 1, header = FALSE)
    x <- as.magpie(x, datacol = 6)
  } else if (subtype == "ccs") {
    x <- read.csv("p_dataccs.csv", sep = ";")
    x <- as.magpie(x, spatial = 1, datacol = 2)
  } else if (subtype == "ffPolyRent") {
    x <- read.csv("ffPolyRent.csv", sep = ";")
    x <- as.magpie(x, spatial = 1, datacol = 5)
  } else if (subtype == "ffPolyCumEx") {
    x <- read.csv("ffPolyCumEx.csv", sep = ";")
    x <- as.magpie(x, spatial = 1, datacol = 5)
  } else if (subtype == "fossilExtractionCoeff") {
    x <- read.csv("fossil_extraction_cost_eq_coefficients.csv", sep = ";")
    # removing the X string added to the column names because how the read.table call inside
    # the read.csv function converts column name numbers to valid variable strings (by using check.names)
    colnames(x) <- gsub("^X", "",  colnames(x))
    x <- as.magpie(x, spatial = 1, temporal = 0, datacol = 3)
    # JPN SSP5 gas extraction zero-order coeff was originally negative.
    x["JPN", , "highGas.0"] <- 0
  } else if (subtype == "uraniumExtractionCoeff") {
    x <- read.csv("uranium_extraction_cost_eq_coefficients.csv", sep = ";")
    x <- as.magpie(x, spatial = 1, temporal = 0, datacol = 3)
  } else if (subtype == "RLDCCoefficientsLoB") {
    x <- read.csv("RLDC_Coefficients_LoB.csv", sep = ";")
    x <- as.magpie(x, spatial = 1, temporal = 0, datacol = 3)
  } else if (subtype == "RLDCCoefficientsPeak") {
    x <- read.csv("RLDC_Coefficients_Peak.csv", sep = ";")
    x <- as.magpie(x, spatial = 1, temporal = 0, datacol = 3)
  } else if (subtype == "earlyRetirementAdjFactor") {
    y <- read.csv("earlyRetirementAdjFactor.csv", sep = ";", skip = 5)
    x <- as.magpie(y, spatial = 1, temporal = 0, datacol = 2)
    x <- setNames(x, colnames(y)[-1])
  } else if (subtype == "nashWeight") {
    x <- read.csv("nash_weights.csv", sep = ";")
    x <- as.magpie(x, spatial = 1, datacol = 2)
  } else {
    stop("Not a valid subtype!")
  }
  return(x)
}