R/calcTimePref.R

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

calcTimePref <- function() {
  
  # read data
  x <- readSource("ExpertGuess",subtype="prtp")
  getNames(x) <- NULL
  getYears(x) <- NULL
  
  # read data used for weight
  pop <- calcOutput("Population",aggregate=FALSE)[,2005,"pop_SSP2"]
  gdp <- calcOutput("GDP",aggregate=FALSE)[,2005,"gdp_SSP2"]
  # use half pop and half gdp 
  w <- (pop / dimSums(pop,dim=1)) + setNames(gdp / dimSums(gdp,dim=1),NULL)
 
  
  return(list(x           = x,
              weight      = w,
              unit        = "dimensionless",
              description = "pure rate of time preference"))
}