R/freight_transport_so_data.R

In russellhz/gcamchartr: GCAM Charting

#' freight_transport_so_data
#'
#' Process transportation service output query - filtering out only trn_freight_road
#' @param query GCAM query containing transportation service output data of one or multiple scenarios
#' @param scenario GCAM scenarios to include in processed data
#' @keywords freight transportation service output
#' @import dplyr tidyr
#' @export
#' @examples
#' freight_transport_so_data("queryA.csv", c("Reference1,date=2017-9-6T13:43:53-07:00", "Reference2,date=2017-9-6T13:43:53-07:00"))

freight_transport_so_data <- function(query, scenarios, query_dir = QUERY_FOLDER){
  tech_order <- c("Liquids", "Hybrid Liquids", "NG", "BEV", "FCEV")

  query_title <- query_id(query_dir) %>%
    filter(file == query) %>%
    select(title) %>%
    as.character

  FSO <- read_query(paste0(query_dir,query), skip = 1) %>%
    filter(scenario != query_title, scenario != "scenario",
           scenario %in% scenarios,
           sector == "trn_freight_road") %>%
    gather(year, value, `1990`:`2100`) %>%
    mutate(year = as.integer(year)) %>%
    filter(year >= 2010) %>%
    group_by(scenario, region, technology, Units, year) %>%
    summarise(value = sum(value)) %>%
    ungroup() %>%
    mutate(value = value/1000,
           Units = "billion-ton-km",
           technology = factor(technology, levels = tech_order),
           scenario = if_else(grepl(",date", scenario),
                              substr(scenario, 1, regexpr(",date", scenario)[1]-1),
                              scenario))
  attr(FSO, "query") <- paste("Freight", query_title)
  attr(FSO, "colors") <- TRANSPORT_COLORS
  attr(FSO, "fill") <- "technology"
  attr(FSO, "default_plot") <- "bar"
  return(FSO)
}