R/convertNREL.R

In pik-piam/moinput: Model Operations Input Data Library

#' Convert NREL data
#' 
#' Convert NREL data on ISO country level.
#' 
#' 
#' @param x MAgPIE object containing NREL data country-region resolution
#' @return NRELWirsenius data as MAgPIE object aggregated to country level
#' @author Lavinia Baumstark
#' @examples
#' 
#' \dontrun{ a <- convertNREL(x,subtype="onshore")
#' }
#' @importFrom madrat toolCountryFill
convertNREL <- function(x) {
  
  # rename countries with ISO-code
  getRegions(x) <- toolCountry2isocode(getRegions(x))
  
  #-------------------- allocation of aggretations -------------------------------
  # "ANT" -> "SXM", "CUW", "BES"
  # "YUG" -> "SRB", "MNE"  # YUG seems to be only these two, "SVN", "HRV", "MKD", "BIH"
  m <- matrix(c(c("ANT","ANT","ANT",       "YUG","YUG"), #,"YUG","YUG","YUG","YUG"),
                c("SXM","CUW","BES",       "SRB","MNE")),5) #,"SVN","HRV","MKD","BIH")),9)
  w <- calcOutput("Population",aggregate=FALSE)[m[,2],2010,"pop_SSP2"]
  x_split <- toolAggregate(x[c("ANT","YUG"),,],m,weight=w)
  # delete ANT and YUG from x
  x <- x[c("ANT","YUG"),,invert=TRUE]
  x <- mbind(x,x_split)
  
  #--------------------------------------------------------------------------------
  
  # fill all missing countries with 0
  x <- toolCountryFill(x,fill=0)
  
  return(x)
}