R/calcSteelStock.R

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

#' Calculate Steel Stock from Mueller steel stock per capita and WDI population
#'
#' @md
#' @return A [`magpie`][magclass::magclass] object.
#'
#' @author Falk Benke
#'
#' @importFrom dplyr filter select mutate left_join %>%
#' @export

calcSteelStock <- function() {
  steel_stock_per_capita <- readSource(type = "Mueller", subtype = "stocks", convert = TRUE)[, , "Steel stock per-capita|med (t)"]
  population <- calcOutput("PopulationPast", PopulationPast = 'UN_PopDiv', aggregate = FALSE)

  y <- intersect(getItems(steel_stock_per_capita, 2), getItems(population, 2))

  stock <- steel_stock_per_capita[, y, ] * population[, y, ]
  getNames(stock) <- c("Steel stock (million t)")

  return(list(
    x = stock,
    weight = NULL,
    unit = "million t",
    description = "Historical steel stock based on Mueller's medium stock per capita numbers and WDI population"
  ))
}