R/zchunk_L254.transportation_korea.R
In rohmin9122/gcam-korea-release: Build the GCAM-Korea Input Datasets
Defines functions
module_gcam.korea_L254.transportation_korea
#' module_gcam.korea_L254.transportation_korea
#'
#' Generates gcam-korea model inputs for transportation sector by states.
#'
#' @param command API command to execute
#' @param ... other optional parameters, depending on command
#' @return Depends on \code{command}: either a vector of required inputs,
#' a vector of output names, or (if \code{command} is "MAKE") all
#' the generated outputs: \code{L254.DeleteSupplysector_koreaTrn}, \code{L254.DeleteFinalDemand_koreaTrn},
#' \code{L254.Supplysector_trn_korea}, \code{L254.FinalEnergyKeyword_trn_korea}, \code{L254.tranSubsectorLogit_korea},
#' \code{L254.tranSubsectorShrwtFllt_korea}, \code{L254.tranSubsectorInterp_korea}, \code{L254.tranSubsectorSpeed_korea},
#' \code{L254.tranSubsectorSpeed_passthru_korea}, \code{L254.tranSubsectorSpeed_noVOTT_korea},
#' \code{L254.tranSubsectorSpeed_nonmotor_korea}, \code{L254.tranSubsectorVOTT_korea}, \code{L254.tranSubsectorFuelPref_korea},
#' \code{L254.StubTranTech_korea}, \code{L254.StubTranTech_passthru_korea}, \code{L254.StubTranTech_nonmotor_korea},
#' \code{L254.StubTranTechLoadFactor_korea}, \code{L254.StubTranTechCost_korea}, \code{L254.StubTranTechCoef_korea},
#' \code{L254.PerCapitaBased_trn_korea}, \code{L254.PriceElasticity_trn_korea}, \code{L254.IncomeElasticity_trn_korea},
#' \code{L254.StubTranTechCalInput_korea}, \code{L254.StubTranTechProd_nonmotor_korea}, \code{L254.StubTranTechCalInput_passthru_korea},
#' \code{L254.BaseService_trn_korea}. The corresponding file in the
#' original data system was \code{L254.transportation_korea.R} (gcam-korea level2).
#' @details This chunk generates input files for transportation sector with generic information for supplysector,
#' subsector and technologies, as well as calibrated inputs and outputs by the US states.
#' @note The transportation structure is heavily nested. The GCAM structure of sector/subsector/technology only
#' allows two levels of nesting within any sector, but a technology of one sector (e.g., trn_pass) can consume the
#' output of another "sector" (e.g., trn_pass_road) that is really just used to represent lower nesting levels of
#' that first, or parent, sector. In the transportation sector, each lower-level nesting "sector" is named by
#' appending a string to the parent sector. So, \code{trn_pass} contains \code{trn_pass_road} which has
#' \code{trn_pass_road_LDV} which has \code{trn_pass_road_LDV_4W}. Each of the links between any two of those sectors
#' is done with a pass-through technology within the parent sector that consumes the output of the child sector.
#' The technology is called a "pass-through" because it (generally) only consumes the output of the child "sector"
#' without making any changes to it. There's an additional complication in the transportation sector, that the
#' pass-through technologies are normal, standard GCAM technologies, not "tranTechnologies" which have different
#' parameters read in, and perform a bunch of hard-wired unit conversions between inputs and outputs.
#' @importFrom assertthat assert_that
#' @importFrom dplyr filter mutate select
#' @importFrom tidyr gather spread
#' @author S. Jeon
module_gcam.korea_L254.transportation_korea <- function(command, ...) {
if(command == driver.DECLARE_INPUTS) {
return(c(FILE = "energy/mappings/UCD_techs",
FILE = "energy/A54.globaltech_nonmotor",
FILE = "energy/A54.globaltech_passthru",
FILE = "energy/A54.sector",
FILE = "gcam-korea/states_subregions",
FILE = "gcam-korea/trn_road_car_intensity",
FILE = "gcam-korea/trn_road_car_loadfactor",
FILE = "gcam-korea/kor_trn_inc_elasticity_for_realitycheck",
"L254.Supplysector_trn",
"L254.FinalEnergyKeyword_trn",
"L254.tranSubsectorLogit",
"L254.tranSubsectorShrwtFllt",
"L254.tranSubsectorInterp",
"L254.tranSubsectorSpeed",
"L254.tranSubsectorSpeed_passthru",
"L254.tranSubsectorSpeed_noVOTT",
"L254.tranSubsectorSpeed_nonmotor",
"L254.tranSubsectorVOTT",
"L254.tranSubsectorFuelPref",
"L254.StubTranTech",
"L254.StubTech_passthru",
"L254.StubTech_nonmotor",
"L254.StubTranTechLoadFactor",
"L254.StubTranTechCost",
"L254.StubTranTechCoef",
"L254.PerCapitaBased_trn",
"L254.PriceElasticity_trn",
"L254.IncomeElasticity_trn",
"L154.in_EJ_korea_trn_m_sz_tech_F",
"L154.out_mpkm_korea_trn_nonmotor_Yh"))
} else if(command == driver.DECLARE_OUTPUTS) {
return(c("L254.DeleteSupplysector_koreaTrn",
"L254.DeleteFinalDemand_koreaTrn",
"L254.Supplysector_trn_korea",
"L254.FinalEnergyKeyword_trn_korea",
"L254.tranSubsectorLogit_trn_korea",
"L254.tranSubsectorShrwtFllt_trn_korea",
"L254.tranSubsectorInterp_trn_korea",
"L254.tranSubsectorSpeed_trn_korea",
"L254.tranSubsectorSpeed_passthru_trn_korea",
"L254.tranSubsectorSpeed_noVOTT_trn_korea",
"L254.tranSubsectorSpeed_nonmotor_trn_korea",
"L254.tranSubsectorVOTT_trn_korea",
"L254.tranSubsectorFuelPref_trn_korea",
"L254.StubTranTech_trn_korea",
"L254.StubTranTech_passthru_trn_korea",
"L254.StubTranTech_nonmotor_trn_korea",
"L254.StubTranTechLoadFactor_trn_korea",
"L254.StubTranTechCost_trn_korea",
"L254.StubTranTechCoef_trn_korea",
"L254.PerCapitaBased_trn_korea",
"L254.PriceElasticity_trn_korea",
"L254.IncomeElasticity_trn_korea",
"L254.StubTranTechCalInput_trn_korea",
"L254.StubTranTechProd_nonmotor_korea",
"L254.StubTranTechCalInput_passthru_trn_korea",
"L254.BaseService_trn_korea"))
} else if(command == driver.MAKE) {
all_data <- list(...)[[1]]
# Silence package notes
region <- supplysector <- technology <- minicam.energy.input <- year <- value <- coefficient <-
UCD_fuel <- UCD_sector <- UCD_technology <- calibrated.value <- coefficient <- loadFactor <-
size.class <- state <- tranSubsector <- tranTechnology <- calibrated.value <- stub.technology <-
output <- output_agg <- output_cum <- share.weight <- subs.share.weight <- share.weight.year <-
tech.share.weight <- calOutputValue <- energy.final.demand <- base.service <- . <- NULL
# Load required inputs
UCD_techs <- get_data(all_data, "energy/mappings/UCD_techs") # Mapping file of transportation technology from the UC Davis report (Mishra et al. 2013)
A54.globaltech_nonmotor <- get_data(all_data, "energy/A54.globaltech_nonmotor")
A54.globaltech_passthru <- get_data(all_data, "energy/A54.globaltech_passthru")
A54.sector <- get_data(all_data, "energy/A54.sector")
states_subregions <- get_data(all_data, "gcam-korea/states_subregions")
trn_road_car_intensity <- get_data(all_data, "gcam-korea/trn_road_car_intensity")
trn_road_car_loadfactor <- get_data(all_data, "gcam-korea/trn_road_car_loadfactor")
kor_trn_inc_elasticity_for_realitycheck <- get_data(all_data, "gcam-korea/kor_trn_inc_elasticity_for_realitycheck")
L254.Supplysector_trn <- get_data(all_data, "L254.Supplysector_trn")
L254.FinalEnergyKeyword_trn <- get_data(all_data, "L254.FinalEnergyKeyword_trn")
L254.tranSubsectorLogit <- get_data(all_data, "L254.tranSubsectorLogit")
L254.tranSubsectorShrwtFllt <- get_data(all_data, "L254.tranSubsectorShrwtFllt")
L254.tranSubsectorInterp <- get_data(all_data, "L254.tranSubsectorInterp")
L254.tranSubsectorSpeed <- get_data(all_data, "L254.tranSubsectorSpeed")
L254.tranSubsectorSpeed_passthru <- get_data(all_data, "L254.tranSubsectorSpeed_passthru")
L254.tranSubsectorSpeed_noVOTT <- get_data(all_data, "L254.tranSubsectorSpeed_noVOTT")
L254.tranSubsectorSpeed_nonmotor <- get_data(all_data, "L254.tranSubsectorSpeed_nonmotor")
L254.tranSubsectorVOTT <- get_data(all_data, "L254.tranSubsectorVOTT")
L254.tranSubsectorFuelPref <- get_data(all_data, "L254.tranSubsectorFuelPref")
L254.StubTranTech <- get_data(all_data, "L254.StubTranTech")
L254.StubTech_passthru <- get_data(all_data, "L254.StubTech_passthru")
L254.StubTech_nonmotor <- get_data(all_data, "L254.StubTech_nonmotor")
L254.StubTranTechLoadFactor <- get_data(all_data, "L254.StubTranTechLoadFactor")
L254.StubTranTechCost <- get_data(all_data, "L254.StubTranTechCost")
L254.StubTranTechCoef <- get_data(all_data, "L254.StubTranTechCoef")
L254.PerCapitaBased_trn <- get_data(all_data, "L254.PerCapitaBased_trn")
L254.PriceElasticity_trn <- get_data(all_data, "L254.PriceElasticity_trn")
L254.IncomeElasticity_trn <- get_data(all_data, "L254.IncomeElasticity_trn")
L154.in_EJ_korea_trn_m_sz_tech_F <- get_data(all_data, "L154.in_EJ_korea_trn_m_sz_tech_F")
L154.out_mpkm_korea_trn_nonmotor_Yh <- get_data(all_data, "L154.out_mpkm_korea_trn_nonmotor_Yh")
# Need to delete the transportation sector in the korea region (energy-final-demands and supplysectors)
# L254.DeleteSupplysector_koreaTrn: Delete transportation supplysectors of the korea region
L254.Supplysector_trn %>%
mutate(region = region) %>% # strip off attributes like title, etc.
filter(region == gcamkorea.REGION) %>%
select(region, supplysector) ->
L254.DeleteSupplysector_koreaTrn
# L254.DeleteFinalDemand_koreaTrn: Delete energy final demand sectors of the korea region
L254.PerCapitaBased_trn %>%
mutate(region = region) %>% # strip off attributes like title, etc.
filter(region == gcamkorea.REGION) %>%
select(LEVEL2_DATA_NAMES[["EnergyFinalDemand"]]) ->
L254.DeleteFinalDemand_koreaTrn
# Process tables at the USA region level to the states level.
# All tables for which processing is identical are done by a function.
# This applies to the supplysectors, subsectors, and stub tech characteristics of the states.
process_korea_to_states <- function(data) {
state <- region <- grid_region <- subsector <- market.name <-
minicam.energy.input <- NULL # silence package check notes
data_new <- data %>%
filter(region == gcamkorea.REGION) %>%
write_to_all_states(names = c(names(data), "region"), region_states=gcamkorea.STATES)
# Re-set markets from USA to regional markets, if called for in the GCAM-USA assumptions for selected fuels
if(gcamkorea.USE_REGIONAL_FUEL_MARKETS & "market.name" %in% names(data_new)) {
data_new <- data_new %>%
left_join_error_no_match(select(states_subregions, state, grid_region), by = c("region" = "state")) %>%
mutate(market.name = replace(market.name, minicam.energy.input %in% gcamkorea.REGIONAL_FUEL_MARKETS,
grid_region[minicam.energy.input %in% gcamkorea.REGIONAL_FUEL_MARKETS])) %>%
select(-grid_region)
}
# Electricity is always consumed from state markets
if("market.name" %in% names(data_new)) {
data_new <- data_new %>%
mutate(market.name = replace(market.name, minicam.energy.input %in% gcamkorea.ELECT_TD_SECTORS,
region[minicam.energy.input %in% gcamkorea.ELECT_TD_SECTORS]))
}
data_new
}
process_korea_to_states(L254.Supplysector_trn) -> L254.Supplysector_trn_korea
process_korea_to_states(L254.FinalEnergyKeyword_trn) -> L254.FinalEnergyKeyword_trn_korea
process_korea_to_states(L254.tranSubsectorLogit) -> L254.tranSubsectorLogit_trn_korea
process_korea_to_states(L254.tranSubsectorShrwtFllt) -> L254.tranSubsectorShrwtFllt_trn_korea
process_korea_to_states(L254.tranSubsectorInterp) -> L254.tranSubsectorInterp_trn_korea
process_korea_to_states(L254.tranSubsectorSpeed) -> L254.tranSubsectorSpeed_trn_korea
process_korea_to_states(L254.tranSubsectorSpeed_passthru) -> L254.tranSubsectorSpeed_passthru_trn_korea
process_korea_to_states(L254.tranSubsectorSpeed_noVOTT) -> L254.tranSubsectorSpeed_noVOTT_trn_korea
process_korea_to_states(L254.tranSubsectorSpeed_nonmotor) -> L254.tranSubsectorSpeed_nonmotor_trn_korea
process_korea_to_states(L254.tranSubsectorVOTT) -> L254.tranSubsectorVOTT_trn_korea
process_korea_to_states(L254.tranSubsectorFuelPref) -> L254.tranSubsectorFuelPref_trn_korea
process_korea_to_states(L254.StubTranTech) -> L254.StubTranTech_trn_korea
process_korea_to_states(L254.StubTech_passthru) -> L254.StubTranTech_passthru_trn_korea
process_korea_to_states(L254.StubTech_nonmotor) -> L254.StubTranTech_nonmotor_trn_korea
process_korea_to_states(L254.StubTranTechLoadFactor) -> L254.StubTranTechLoadFactor_trn_korea
####### Written by Jeon ####### Start
## Replace default load factor with KOTI loadfactor
trn_road_car_loadfactor %>% rename(tranSubsector = GCAM_mode) -> trn_road_car_loadfactor
year_2020_2100 <- seq(2020, 2100, by = 5)
for ( i in 1:length(year_2020_2100)) {
trn_road_car_loadfactor %>%
filter( year == 2015 ) %>%
mutate( year = year_2020_2100[i]) -> each_data
trn_road_car_loadfactor %>%
bind_rows(each_data) -> trn_road_car_loadfactor
}
L254.StubTranTechLoadFactor_trn_korea %>%
left_join(trn_road_car_loadfactor) %>%
mutate(loadFactor = ifelse(is.na(loadfactor), loadFactor, loadfactor)) %>%
select(region, supplysector, tranSubsector, stub.technology, year, loadFactor) -> L254.StubTranTechLoadFactor_trn_korea
####### Written by Jeon ####### End
process_korea_to_states(L254.StubTranTechCost) -> L254.StubTranTechCost_trn_korea
L254.StubTranTechCoef %>%
mutate(coefficient = round(coefficient, digits = gcamkorea.DIGITS_TRNUSA_DEFAULT)) %>%
process_korea_to_states ->
L254.StubTranTechCoef_trn_korea
####### Written by Jeon ####### Start
## Replace default fuel intensity with fuel intensity calculated from TS and ECS(KEEI)
trn_road_car_intensity %>%
rename(tranSubsector = GCAM_mode) -> trn_road_car_intensity
L254.StubTranTechCoef_trn_korea %>%
left_join(trn_road_car_intensity, by = c("tranSubsector", "stub.technology", "region")) %>%
mutate(coefficient = ifelse(is.na(intensity), coefficient, intensity)) %>%
select(region, supplysector, tranSubsector, stub.technology, year.x, minicam.energy.input, coefficient, market.name) %>%
rename(year = year.x)-> L254.StubTranTechCoef_trn_korea
####### Written by Jeon ####### End
process_korea_to_states(L254.PerCapitaBased_trn) -> L254.PerCapitaBased_trn_korea
process_korea_to_states(L254.PriceElasticity_trn) -> L254.PriceElasticity_trn_korea
process_korea_to_states(L254.IncomeElasticity_trn) -> L254.IncomeElasticity_trn_korea
## Written by Jeon ## Start
# Make historical tredns plausible
L254.IncomeElasticity_trn_korea %>%
left_join(kor_trn_inc_elasticity_for_realitycheck, by = c("region", "energy.final.demand", "year")) %>%
mutate( income.elasticity = ifelse( is.na(income.elasticity.y),
income.elasticity.x,
income.elasticity.y)) %>%
select(region, energy.final.demand, year, income.elasticity) -> L254.IncomeElasticity_trn_korea
####### Written by Jeon ####### End
# Calibration
# L254.StubTranTechCalInput_trn_korea: calibrated energy consumption by all technologies
L154.in_EJ_korea_trn_m_sz_tech_F %>%
filter(year %in% MODEL_BASE_YEARS) %>%
mutate(calibrated.value = round(value, digits = energy.DIGITS_CALOUTPUT),
region = state) %>%
left_join_error_no_match(select(UCD_techs, UCD_sector, mode, size.class, UCD_technology, UCD_fuel,
supplysector, tranSubsector, stub.technology = tranTechnology, minicam.energy.input),
by = c("UCD_sector", "mode", "size.class", "UCD_technology", "UCD_fuel")) %>%
select(LEVEL2_DATA_NAMES[["StubTranTech"]], year, minicam.energy.input, calibrated.value) ->
L254.StubTranTechCalInput_trn_korea
# NOTE: NEED TO WRITE THIS OUT FOR ALL TECHNOLOGIES, NOT JUST THOSE THAT EXIST IN SOME BASE YEARS.
# Model may make up calibration values otherwise.
L254.StubTranTechCoef_trn_korea %>%
filter(year %in% MODEL_BASE_YEARS) %>%
select(names(.)[names(.) %in% LEVEL2_DATA_NAMES[["StubTranTechCalInput"]]]) %>%
left_join(L254.StubTranTechCalInput_trn_korea,
by = c("region", "supplysector", "tranSubsector", "stub.technology", "year", "minicam.energy.input")) %>%
# Set calibration values as zero for technolgies that do not exist in some base years
replace_na(list(calibrated.value = 0)) %>%
mutate(share.weight.year = year,
# Create the needed variables to use the function set_subsector_shrwt
subsector = tranSubsector, calOutputValue = calibrated.value) %>%
set_subsector_shrwt() %>%
mutate(tech.share.weight = if_else(calibrated.value > 0, 1, 0)) %>%
select(LEVEL2_DATA_NAMES[["StubTranTechCalInput"]]) ->
L254.StubTranTechCalInput_trn_korea
# Non-motorized technologies
# L254.StubTranTechProd_nonmotor_korea: service output of non-motorized transportation technologies
#minyoung
L154.out_mpkm_korea_trn_nonmotor_Yh %>%
filter(year %in% MODEL_BASE_YEARS) %>%
mutate(calOutputValue = round(value, digits = energy.DIGITS_MPKM),
region = state, tranSubsector = mode) %>%
left_join_error_no_match(A54.globaltech_nonmotor, by = "tranSubsector") %>%
mutate(stub.technology = technology) %>%
# There is no need to match shareweights to the calOutputValue because no region should ever have a 0 here
select(LEVEL2_DATA_NAMES[["StubTranTech"]], year, calOutputValue) ->
L254.StubTranTechProd_nonmotor_korea
# L254.StubTranTechCalInput_passthru_trn_korea: calibrated input of passthrough technologies
# trn_pass, trn_pass_road, trn_pass_road_LDV, trn_freight
# The transportation structure is heavily nested.
# The GCAM structure of sector/subsector/technology only allows two levels of nesting within any sector,
# but a technology of one sector (e.g., trn_pass) can consume the output of another "sector" (e.g., trn_pass_road)
# that is really just used to represent lower nesting levels of that first, or parent, sector. In the
# transportation sector, each lower-level nesting "sector" is named by appending a string to the parent sector.
# So, trn_pass contains trn_pass_road which has trn_pass_road_LDV which has trn_pass_road_LDV_4W. Each of the links
# between any two of those sectors is done with a pass-through technology within the parent sector that consumes
# the output of the child sector. The technology is called a "pass-through" because it (generally) only consumes
# the output of the child "sector" without making any changes to it. There's an additional complication in the
# transportation sector: the pass-through technologies are normal, standard GCAM technologies, not "tranTechnologies"
# which have different parameters read in, and perform a bunch of hard-wired unit conversions between inputs and outputs.
# First, need to calculate the service output for all tranTechnologies
# calInput * loadFactor * unit_conversion / (coef * unit conversion)
L254.StubTranTechCalInput_trn_korea %>%
left_join_error_no_match(L254.StubTranTechLoadFactor_trn_korea,
by = c("region", "supplysector", "tranSubsector", "stub.technology", "year")) %>%
left_join_error_no_match(L254.StubTranTechCoef_trn_korea,
by = c("region", "supplysector", "tranSubsector", "stub.technology", "year", "minicam.energy.input")) %>%
mutate(output = round(calibrated.value * loadFactor * CONV_EJ_GJ / (coefficient * CONV_BTU_KJ),
digits = gcamkorea.DIGITS_TRNUSA_DEFAULT)) ->
L254.StubTranTechOutput_korea
# The next step is to calculate the aggregated outputs by supplysector
# Outputs of certain supplysectors are inputs for the passthrough technologies
L254.StubTranTechOutput_korea %>%
group_by(region, year, supplysector) %>%
summarise(output_agg = sum(output)) %>%
ungroup() ->
L254.StubTranTechOutput_korea_agg
# Write all possible pass-through technologies to all regions
A54.globaltech_passthru %>%
repeat_add_columns(tibble(year = MODEL_BASE_YEARS)) %>%
write_to_all_states(names = c(names(.), "region"), region_states=gcamkorea.STATES) %>%
select(region, supplysector, tranSubsector, stub.technology = technology, year, minicam.energy.input) %>%
# Subset only the passthrough technologies that are applicable in each region
semi_join(L254.StubTranTech_passthru_trn_korea,
by = c("region", "supplysector", "tranSubsector", "stub.technology")) %>%
# Match in outputs of supplysectors that are inputs for the passthrough technologies
left_join(L254.StubTranTechOutput_korea_agg,
by = c("region", "year", "minicam.energy.input" = "supplysector")) %>%
# Some of the technologies are sub-totals, assign zero value now, will be calculated below
replace_na(list(output_agg = 0)) %>%
# Arrange input sectors so that sub-total sector is behind the subsectors
arrange(desc(minicam.energy.input)) %>%
group_by(region, year) %>%
# Calculate the cumulative for sub-total sector
mutate(output_cum = cumsum(output_agg)) %>%
ungroup() ->
L254.StubTranTechCalInput_passthru_trn_korea_cum
# Prepare a list of the supplysector in the passthrough input table to filter the sub-total sectors
LIST_supplysector <- unique(L254.StubTranTechCalInput_passthru_trn_korea_cum$supplysector)
L254.StubTranTechCalInput_passthru_trn_korea_cum %>%
# Use the cumulative value for sub-total sectors
mutate(calibrated.value = if_else(minicam.energy.input %in% LIST_supplysector,
output_cum, output_agg)) %>%
mutate(share.weight.year = year,
subs.share.weight = if_else(calibrated.value > 0, 1, 0),
tech.share.weight = if_else(calibrated.value > 0, 1, 0)) %>%
select(LEVEL2_DATA_NAMES[["StubTranTechCalInput"]]) ->
L254.StubTranTechCalInput_passthru_trn_korea
# L254.BaseService_trn_korea: base-year service output of transportation final demand
L254.StubTranTechOutput_korea %>%
select(LEVEL2_DATA_NAMES[["StubTranTech"]], year, base.service = output) %>%
bind_rows(L254.StubTranTechProd_nonmotor_korea %>%
select(LEVEL2_DATA_NAMES[["StubTranTech"]], year, base.service = calOutputValue)) %>%
left_join_error_no_match(select(A54.sector, supplysector, energy.final.demand), by = "supplysector") %>%
group_by(region, energy.final.demand, year) %>%
summarise(base.service = sum(base.service)) %>%
ungroup ->
L254.BaseService_trn_korea
# #minyoung
# L254.StubTranTechOutput_korea %>%
# select(LEVEL2_DATA_NAMES[["StubTranTech"]], year, base.service = output) %>%
# left_join_error_no_match(select(A54.sector, supplysector, energy.final.demand), by = "supplysector") %>%
# group_by(region, energy.final.demand, year) %>%
# summarise(base.service = sum(base.service)) %>%
# ungroup ->
# L254.BaseService_trn_korea
# Produce outputs
L254.DeleteSupplysector_koreaTrn %>%
add_title("Delect transportation supply sectors of the full USA region") %>%
add_units("NA") %>%
add_comments("Delect transportation supply sectors of the full USA region") %>%
add_legacy_name("L254.DeleteSupplysector_koreaTrn") %>%
add_precursors("L254.Supplysector_trn") ->
L254.DeleteSupplysector_koreaTrn
L254.DeleteFinalDemand_koreaTrn %>%
add_title("Delete energy final demand sectors of the full USA region") %>%
add_units("NA") %>%
add_comments("Delete energy final demand sectors of the full USA region") %>%
add_legacy_name("L254.DeleteFinalDemand_koreaTrn") %>%
add_precursors("L254.PerCapitaBased_trn") ->
L254.DeleteFinalDemand_koreaTrn
L254.Supplysector_trn_korea %>%
add_title("Supply sector information for transportation sector in the US states") %>%
add_units("Unitless") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.Supplysector_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.Supplysector_trn") ->
L254.Supplysector_trn_korea
L254.FinalEnergyKeyword_trn_korea %>%
add_title("Supply sector final energy keywords for transportation sector in the US states") %>%
add_units("NA") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.FinalEnergyKeyword_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.FinalEnergyKeyword_trn") ->
L254.FinalEnergyKeyword_trn_korea
L254.tranSubsectorLogit_trn_korea %>%
add_title("Subsector logit exponents of transportation sector in the US states") %>%
add_units("Unitless") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.tranSubsectorLogit_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.tranSubsectorLogit") ->
L254.tranSubsectorLogit_trn_korea
L254.tranSubsectorShrwtFllt_trn_korea %>%
add_title("Subsector shareweights of transportation sector in the US states") %>%
add_units("Unitless") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.tranSubsectorShrwtFllt_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.tranSubsectorShrwtFllt") ->
L254.tranSubsectorShrwtFllt_trn_korea
L254.tranSubsectorInterp_trn_korea %>%
add_title("Temporal subsector shareweight interpolation of transportation sector in the US states") %>%
add_units("Unitless") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.tranSubsectorInterp_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.tranSubsectorInterp") ->
L254.tranSubsectorInterp_trn_korea
L254.tranSubsectorSpeed_trn_korea %>%
add_title("Speeds of transportation modes (not including pass-through sectors) in the US states") %>%
add_units("km / hr") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.tranSubsectorSpeed_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.tranSubsectorSpeed") ->
L254.tranSubsectorSpeed_trn_korea
L254.tranSubsectorSpeed_passthru_trn_korea %>%
add_title("Speeds of pass-through transportation subsectors in the US states") %>%
add_units("km / hr") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.tranSubsectorSpeed_passthru_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.tranSubsectorSpeed_passthru") ->
L254.tranSubsectorSpeed_passthru_trn_korea
L254.tranSubsectorSpeed_noVOTT_trn_korea %>%
add_title("Speeds of transportation subsectors whose time value is not considered in the US states") %>%
add_units("km / hr") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.tranSubsectorSpeed_noVOTT_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.tranSubsectorSpeed_noVOTT") ->
L254.tranSubsectorSpeed_noVOTT_trn_korea
L254.tranSubsectorSpeed_nonmotor_trn_korea %>%
add_title("Speeds of non-motorized transportation subsectors in the US states") %>%
add_units("km / hr") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.tranSubsectorSpeed_nonmotor_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.tranSubsectorSpeed_nonmotor") ->
L254.tranSubsectorSpeed_nonmotor_trn_korea
L254.tranSubsectorVOTT_trn_korea %>%
add_title("descriptive title of data") %>%
add_units("units") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.tranSubsectorVOTT_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.tranSubsectorVOTT") ->
L254.tranSubsectorVOTT_trn_korea
L254.tranSubsectorFuelPref_trn_korea %>%
add_title("Subsector (fuel) preferences elasticity that are tied to GDP in the US states") %>%
add_units("Unitless") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_comments("Fuel preferences are unrelated to time value") %>%
add_legacy_name("L254.tranSubsectorFuelPref_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.tranSubsectorFuelPref") ->
L254.tranSubsectorFuelPref_trn_korea
L254.StubTranTech_trn_korea %>%
add_title("Transportation stub technologies in the US states") %>%
add_units("NA") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.StubTranTech_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.StubTranTech") ->
L254.StubTranTech_trn_korea
L254.StubTranTech_passthru_trn_korea %>%
add_title("Transportation stub technologies for passthrough sectors in the US states") %>%
add_units("NA") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name(" L254.StubTranTech_passthru_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.StubTech_passthru") ->
L254.StubTranTech_passthru_trn_korea
L254.StubTranTech_nonmotor_trn_korea %>%
add_title("Transportation stub technologies for non-motorized subsectors in the US states") %>%
add_units("NA") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.StubTranTech_nonmotor_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.StubTech_nonmotor") ->
L254.StubTranTech_nonmotor_trn_korea
L254.StubTranTechLoadFactor_trn_korea %>%
add_title("Load factors of transportation stub technologies in the US states") %>%
add_units("person/vehicle and tonnes/vehicle") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.StubTranTechLoadFactor_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.StubTranTechLoadFactor",
"gcam-korea/trn_road_car_loadfactor") ->
L254.StubTranTechLoadFactor_trn_korea
L254.StubTranTechCost_trn_korea %>%
add_title("Costs of transportation stub technologies in the US states") %>%
add_units("$1990USD / vkm") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.StubTranTechCost_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.StubTranTechCost") ->
L254.StubTranTechCost_trn_korea
L254.StubTranTechCoef_trn_korea %>%
add_title("Coefficients of transportation stub technologies in the US states") %>%
add_units("BTU / vkm") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_comments("Re-set electricity consumed at the state markets") %>%
add_legacy_name("L254.StubTranTechCoef_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.StubTranTechCoef",
"gcam-korea/trn_road_car_intensity") ->
L254.StubTranTechCoef_trn_korea
L254.PerCapitaBased_trn_korea %>%
add_title("Per-capita based flag for transportation final demand in the US states") %>%
add_units("NA") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.PerCapitaBased_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.PerCapitaBased_trn") ->
L254.PerCapitaBased_trn_korea
L254.PriceElasticity_trn_korea %>%
add_title("Price elasticity of transportation final demand in the US states") %>%
add_units("Unitless") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.PriceElasticity_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.PriceElasticity_trn") ->
L254.PriceElasticity_trn_korea
L254.IncomeElasticity_trn_korea %>%
add_title("Income elasticity of transportation final demand in the US states") %>%
add_units("Unitless") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.IncomeElasticity_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.IncomeElasticity_trn",
"gcam-korea/kor_trn_inc_elasticity_for_realitycheck") ->
L254.IncomeElasticity_trn_korea
L254.StubTranTechCalInput_trn_korea %>%
add_title("Calibrated energy consumption by all transportation stub technologies in the US states") %>%
add_units("EJ") %>%
add_comments("Set calibration values for those technologies that do not exist in some base years as zero") %>%
add_legacy_name("L254.StubTranTechCalInput_trn_korea") %>%
same_precursors_as("L254.StubTranTechCoef_trn_korea") %>%
add_precursors("L154.in_EJ_korea_trn_m_sz_tech_F",
"energy/mappings/UCD_techs") ->
L254.StubTranTechCalInput_trn_korea
L254.StubTranTechProd_nonmotor_korea %>%
add_title("Calibrated service output of non-motorized transportation technologies in the US states") %>%
add_units("Million pass-km") %>%
add_comments("Not match shareweights to the calOutputValue because no region should ever have a zero here") %>%
add_legacy_name("L254.StubTranTechProd_nonmotor_korea") %>%
add_precursors("L154.out_mpkm_korea_trn_nonmotor_Yh",
"energy/A54.globaltech_nonmotor") ->
L254.StubTranTechProd_nonmotor_korea
L254.StubTranTechCalInput_passthru_trn_korea %>%
add_title("Calibrated energy consumption of transportation passthrough technologies in the US states") %>%
add_units("EJ") %>%
add_comments("Use outputs of the supplysectors that are inputs for passthrough technologies") %>%
add_comments("Outputs of all motorized technologies are calculated as calInput * loadFactor / coefficient") %>%
add_legacy_name("L254.StubTranTechCalInput_passthru_trn_korea") %>%
same_precursors_as("L254.StubTranTechCalInput_trn_korea") %>%
same_precursors_as("L254.StubTranTechLoadFactor_trn_korea") %>%
same_precursors_as("L254.StubTranTechCoef_trn_korea") %>%
same_precursors_as("L254.StubTranTech_passthru_trn_korea") %>%
add_precursors("energy/A54.globaltech_passthru") ->
L254.StubTranTechCalInput_passthru_trn_korea
L254.BaseService_trn_korea %>%
add_title("Base-year service output of transportation final demand") %>%
add_units("Million pass-km and million ton-km") %>%
add_comments("Service outputs of all motorized technologies are calculated as calInput * loadFactor / coefficient") %>%
add_comments("Combine with service output of non-motorized transportation technologies") %>%
add_legacy_name("L254.BaseService_trn_korea") %>%
same_precursors_as("L254.StubTranTechCalInput_trn_korea") %>%
same_precursors_as("L254.StubTranTechLoadFactor_trn_korea") %>%
same_precursors_as("L254.StubTranTechCoef_trn_korea") %>%
same_precursors_as("L254.StubTranTechProd_nonmotor_korea") %>%
add_precursors("energy/A54.sector") ->
L254.BaseService_trn_korea
return_data(L254.DeleteSupplysector_koreaTrn, L254.DeleteFinalDemand_koreaTrn,
L254.Supplysector_trn_korea,
L254.FinalEnergyKeyword_trn_korea,
L254.tranSubsectorLogit_trn_korea,
L254.tranSubsectorShrwtFllt_trn_korea,
L254.tranSubsectorInterp_trn_korea,
L254.tranSubsectorSpeed_trn_korea,
L254.tranSubsectorSpeed_passthru_trn_korea,
L254.tranSubsectorSpeed_noVOTT_trn_korea,
L254.tranSubsectorSpeed_nonmotor_trn_korea,
L254.tranSubsectorVOTT_trn_korea,
L254.tranSubsectorFuelPref_trn_korea,
L254.StubTranTech_trn_korea,
L254.StubTranTech_passthru_trn_korea,
L254.StubTranTech_nonmotor_trn_korea,
L254.StubTranTechLoadFactor_trn_korea,
L254.StubTranTechCost_trn_korea,
L254.StubTranTechCoef_trn_korea,
L254.PerCapitaBased_trn_korea,
L254.PriceElasticity_trn_korea,
L254.IncomeElasticity_trn_korea,
L254.StubTranTechCalInput_trn_korea,
L254.StubTranTechProd_nonmotor_korea,
L254.StubTranTechCalInput_passthru_trn_korea, L254.BaseService_trn_korea)
} else {
stop("Unknown command")
}
}
rohmin9122/gcam-korea-release documentation built on Nov. 26, 2020, 8:11 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Defines functions module_gcam.korea_L254.transportation_korea
#' module_gcam.korea_L254.transportation_korea
#'
#' Generates gcam-korea model inputs for transportation sector by states.
#'
#' @param command API command to execute
#' @param ... other optional parameters, depending on command
#' @return Depends on \code{command}: either a vector of required inputs,
#' a vector of output names, or (if \code{command} is "MAKE") all
#' the generated outputs: \code{L254.DeleteSupplysector_koreaTrn}, \code{L254.DeleteFinalDemand_koreaTrn},
#' \code{L254.Supplysector_trn_korea}, \code{L254.FinalEnergyKeyword_trn_korea}, \code{L254.tranSubsectorLogit_korea},
#' \code{L254.tranSubsectorShrwtFllt_korea}, \code{L254.tranSubsectorInterp_korea}, \code{L254.tranSubsectorSpeed_korea},
#' \code{L254.tranSubsectorSpeed_passthru_korea}, \code{L254.tranSubsectorSpeed_noVOTT_korea},
#' \code{L254.tranSubsectorSpeed_nonmotor_korea}, \code{L254.tranSubsectorVOTT_korea}, \code{L254.tranSubsectorFuelPref_korea},
#' \code{L254.StubTranTech_korea}, \code{L254.StubTranTech_passthru_korea}, \code{L254.StubTranTech_nonmotor_korea},
#' \code{L254.StubTranTechLoadFactor_korea}, \code{L254.StubTranTechCost_korea}, \code{L254.StubTranTechCoef_korea},
#' \code{L254.PerCapitaBased_trn_korea}, \code{L254.PriceElasticity_trn_korea}, \code{L254.IncomeElasticity_trn_korea},
#' \code{L254.StubTranTechCalInput_korea}, \code{L254.StubTranTechProd_nonmotor_korea}, \code{L254.StubTranTechCalInput_passthru_korea},
#' \code{L254.BaseService_trn_korea}. The corresponding file in the
#' original data system was \code{L254.transportation_korea.R} (gcam-korea level2).
#' @details This chunk generates input files for transportation sector with generic information for supplysector,
#' subsector and technologies, as well as calibrated inputs and outputs by the US states.
#' @note The transportation structure is heavily nested. The GCAM structure of sector/subsector/technology only
#' allows two levels of nesting within any sector, but a technology of one sector (e.g., trn_pass) can consume the
#' output of another "sector" (e.g., trn_pass_road) that is really just used to represent lower nesting levels of
#' that first, or parent, sector. In the transportation sector, each lower-level nesting "sector" is named by
#' appending a string to the parent sector. So, \code{trn_pass} contains \code{trn_pass_road} which has
#' \code{trn_pass_road_LDV} which has \code{trn_pass_road_LDV_4W}. Each of the links between any two of those sectors
#' is done with a pass-through technology within the parent sector that consumes the output of the child sector.
#' The technology is called a "pass-through" because it (generally) only consumes the output of the child "sector"
#' without making any changes to it. There's an additional complication in the transportation sector, that the
#' pass-through technologies are normal, standard GCAM technologies, not "tranTechnologies" which have different
#' parameters read in, and perform a bunch of hard-wired unit conversions between inputs and outputs.
#' @importFrom assertthat assert_that
#' @importFrom dplyr filter mutate select
#' @importFrom tidyr gather spread
#' @author S. Jeon
module_gcam.korea_L254.transportation_korea <- function(command, ...) {
if(command == driver.DECLARE_INPUTS) {
return(c(FILE = "energy/mappings/UCD_techs",
FILE = "energy/A54.globaltech_nonmotor",
FILE = "energy/A54.globaltech_passthru",
FILE = "energy/A54.sector",
FILE = "gcam-korea/states_subregions",
FILE = "gcam-korea/trn_road_car_intensity",
FILE = "gcam-korea/trn_road_car_loadfactor",
FILE = "gcam-korea/kor_trn_inc_elasticity_for_realitycheck",
"L254.Supplysector_trn",
"L254.FinalEnergyKeyword_trn",
"L254.tranSubsectorLogit",
"L254.tranSubsectorShrwtFllt",
"L254.tranSubsectorInterp",
"L254.tranSubsectorSpeed",
"L254.tranSubsectorSpeed_passthru",
"L254.tranSubsectorSpeed_noVOTT",
"L254.tranSubsectorSpeed_nonmotor",
"L254.tranSubsectorVOTT",
"L254.tranSubsectorFuelPref",
"L254.StubTranTech",
"L254.StubTech_passthru",
"L254.StubTech_nonmotor",
"L254.StubTranTechLoadFactor",
"L254.StubTranTechCost",
"L254.StubTranTechCoef",
"L254.PerCapitaBased_trn",
"L254.PriceElasticity_trn",
"L254.IncomeElasticity_trn",
"L154.in_EJ_korea_trn_m_sz_tech_F",
"L154.out_mpkm_korea_trn_nonmotor_Yh"))
} else if(command == driver.DECLARE_OUTPUTS) {
return(c("L254.DeleteSupplysector_koreaTrn",
"L254.DeleteFinalDemand_koreaTrn",
"L254.Supplysector_trn_korea",
"L254.FinalEnergyKeyword_trn_korea",
"L254.tranSubsectorLogit_trn_korea",
"L254.tranSubsectorShrwtFllt_trn_korea",
"L254.tranSubsectorInterp_trn_korea",
"L254.tranSubsectorSpeed_trn_korea",
"L254.tranSubsectorSpeed_passthru_trn_korea",
"L254.tranSubsectorSpeed_noVOTT_trn_korea",
"L254.tranSubsectorSpeed_nonmotor_trn_korea",
"L254.tranSubsectorVOTT_trn_korea",
"L254.tranSubsectorFuelPref_trn_korea",
"L254.StubTranTech_trn_korea",
"L254.StubTranTech_passthru_trn_korea",
"L254.StubTranTech_nonmotor_trn_korea",
"L254.StubTranTechLoadFactor_trn_korea",
"L254.StubTranTechCost_trn_korea",
"L254.StubTranTechCoef_trn_korea",
"L254.PerCapitaBased_trn_korea",
"L254.PriceElasticity_trn_korea",
"L254.IncomeElasticity_trn_korea",
"L254.StubTranTechCalInput_trn_korea",
"L254.StubTranTechProd_nonmotor_korea",
"L254.StubTranTechCalInput_passthru_trn_korea",
"L254.BaseService_trn_korea"))
} else if(command == driver.MAKE) {
all_data <- list(...)[[1]]
# Silence package notes
region <- supplysector <- technology <- minicam.energy.input <- year <- value <- coefficient <-
UCD_fuel <- UCD_sector <- UCD_technology <- calibrated.value <- coefficient <- loadFactor <-
size.class <- state <- tranSubsector <- tranTechnology <- calibrated.value <- stub.technology <-
output <- output_agg <- output_cum <- share.weight <- subs.share.weight <- share.weight.year <-
tech.share.weight <- calOutputValue <- energy.final.demand <- base.service <- . <- NULL
# Load required inputs
UCD_techs <- get_data(all_data, "energy/mappings/UCD_techs") # Mapping file of transportation technology from the UC Davis report (Mishra et al. 2013)
A54.globaltech_nonmotor <- get_data(all_data, "energy/A54.globaltech_nonmotor")
A54.globaltech_passthru <- get_data(all_data, "energy/A54.globaltech_passthru")
A54.sector <- get_data(all_data, "energy/A54.sector")
states_subregions <- get_data(all_data, "gcam-korea/states_subregions")
trn_road_car_intensity <- get_data(all_data, "gcam-korea/trn_road_car_intensity")
trn_road_car_loadfactor <- get_data(all_data, "gcam-korea/trn_road_car_loadfactor")
kor_trn_inc_elasticity_for_realitycheck <- get_data(all_data, "gcam-korea/kor_trn_inc_elasticity_for_realitycheck")
L254.Supplysector_trn <- get_data(all_data, "L254.Supplysector_trn")
L254.FinalEnergyKeyword_trn <- get_data(all_data, "L254.FinalEnergyKeyword_trn")
L254.tranSubsectorLogit <- get_data(all_data, "L254.tranSubsectorLogit")
L254.tranSubsectorShrwtFllt <- get_data(all_data, "L254.tranSubsectorShrwtFllt")
L254.tranSubsectorInterp <- get_data(all_data, "L254.tranSubsectorInterp")
L254.tranSubsectorSpeed <- get_data(all_data, "L254.tranSubsectorSpeed")
L254.tranSubsectorSpeed_passthru <- get_data(all_data, "L254.tranSubsectorSpeed_passthru")
L254.tranSubsectorSpeed_noVOTT <- get_data(all_data, "L254.tranSubsectorSpeed_noVOTT")
L254.tranSubsectorSpeed_nonmotor <- get_data(all_data, "L254.tranSubsectorSpeed_nonmotor")
L254.tranSubsectorVOTT <- get_data(all_data, "L254.tranSubsectorVOTT")
L254.tranSubsectorFuelPref <- get_data(all_data, "L254.tranSubsectorFuelPref")
L254.StubTranTech <- get_data(all_data, "L254.StubTranTech")
L254.StubTech_passthru <- get_data(all_data, "L254.StubTech_passthru")
L254.StubTech_nonmotor <- get_data(all_data, "L254.StubTech_nonmotor")
L254.StubTranTechLoadFactor <- get_data(all_data, "L254.StubTranTechLoadFactor")
L254.StubTranTechCost <- get_data(all_data, "L254.StubTranTechCost")
L254.StubTranTechCoef <- get_data(all_data, "L254.StubTranTechCoef")
L254.PerCapitaBased_trn <- get_data(all_data, "L254.PerCapitaBased_trn")
L254.PriceElasticity_trn <- get_data(all_data, "L254.PriceElasticity_trn")
L254.IncomeElasticity_trn <- get_data(all_data, "L254.IncomeElasticity_trn")
L154.in_EJ_korea_trn_m_sz_tech_F <- get_data(all_data, "L154.in_EJ_korea_trn_m_sz_tech_F")
L154.out_mpkm_korea_trn_nonmotor_Yh <- get_data(all_data, "L154.out_mpkm_korea_trn_nonmotor_Yh")
# Need to delete the transportation sector in the korea region (energy-final-demands and supplysectors)
# L254.DeleteSupplysector_koreaTrn: Delete transportation supplysectors of the korea region
L254.Supplysector_trn %>%
mutate(region = region) %>% # strip off attributes like title, etc.
filter(region == gcamkorea.REGION) %>%
select(region, supplysector) ->
L254.DeleteSupplysector_koreaTrn
# L254.DeleteFinalDemand_koreaTrn: Delete energy final demand sectors of the korea region
L254.PerCapitaBased_trn %>%
mutate(region = region) %>% # strip off attributes like title, etc.
filter(region == gcamkorea.REGION) %>%
select(LEVEL2_DATA_NAMES[["EnergyFinalDemand"]]) ->
L254.DeleteFinalDemand_koreaTrn
# Process tables at the USA region level to the states level.
# All tables for which processing is identical are done by a function.
# This applies to the supplysectors, subsectors, and stub tech characteristics of the states.
process_korea_to_states <- function(data) {
state <- region <- grid_region <- subsector <- market.name <-
minicam.energy.input <- NULL # silence package check notes
data_new <- data %>%
filter(region == gcamkorea.REGION) %>%
write_to_all_states(names = c(names(data), "region"), region_states=gcamkorea.STATES)
# Re-set markets from USA to regional markets, if called for in the GCAM-USA assumptions for selected fuels
if(gcamkorea.USE_REGIONAL_FUEL_MARKETS & "market.name" %in% names(data_new)) {
data_new <- data_new %>%
left_join_error_no_match(select(states_subregions, state, grid_region), by = c("region" = "state")) %>%
mutate(market.name = replace(market.name, minicam.energy.input %in% gcamkorea.REGIONAL_FUEL_MARKETS,
grid_region[minicam.energy.input %in% gcamkorea.REGIONAL_FUEL_MARKETS])) %>%
select(-grid_region)
}
# Electricity is always consumed from state markets
if("market.name" %in% names(data_new)) {
data_new <- data_new %>%
mutate(market.name = replace(market.name, minicam.energy.input %in% gcamkorea.ELECT_TD_SECTORS,
region[minicam.energy.input %in% gcamkorea.ELECT_TD_SECTORS]))
}
data_new
}
process_korea_to_states(L254.Supplysector_trn) -> L254.Supplysector_trn_korea
process_korea_to_states(L254.FinalEnergyKeyword_trn) -> L254.FinalEnergyKeyword_trn_korea
process_korea_to_states(L254.tranSubsectorLogit) -> L254.tranSubsectorLogit_trn_korea
process_korea_to_states(L254.tranSubsectorShrwtFllt) -> L254.tranSubsectorShrwtFllt_trn_korea
process_korea_to_states(L254.tranSubsectorInterp) -> L254.tranSubsectorInterp_trn_korea
process_korea_to_states(L254.tranSubsectorSpeed) -> L254.tranSubsectorSpeed_trn_korea
process_korea_to_states(L254.tranSubsectorSpeed_passthru) -> L254.tranSubsectorSpeed_passthru_trn_korea
process_korea_to_states(L254.tranSubsectorSpeed_noVOTT) -> L254.tranSubsectorSpeed_noVOTT_trn_korea
process_korea_to_states(L254.tranSubsectorSpeed_nonmotor) -> L254.tranSubsectorSpeed_nonmotor_trn_korea
process_korea_to_states(L254.tranSubsectorVOTT) -> L254.tranSubsectorVOTT_trn_korea
process_korea_to_states(L254.tranSubsectorFuelPref) -> L254.tranSubsectorFuelPref_trn_korea
process_korea_to_states(L254.StubTranTech) -> L254.StubTranTech_trn_korea
process_korea_to_states(L254.StubTech_passthru) -> L254.StubTranTech_passthru_trn_korea
process_korea_to_states(L254.StubTech_nonmotor) -> L254.StubTranTech_nonmotor_trn_korea
process_korea_to_states(L254.StubTranTechLoadFactor) -> L254.StubTranTechLoadFactor_trn_korea
####### Written by Jeon ####### Start
## Replace default load factor with KOTI loadfactor
trn_road_car_loadfactor %>% rename(tranSubsector = GCAM_mode) -> trn_road_car_loadfactor
year_2020_2100 <- seq(2020, 2100, by = 5)
for ( i in 1:length(year_2020_2100)) {
trn_road_car_loadfactor %>%
filter( year == 2015 ) %>%
mutate( year = year_2020_2100[i]) -> each_data
trn_road_car_loadfactor %>%
bind_rows(each_data) -> trn_road_car_loadfactor
}
L254.StubTranTechLoadFactor_trn_korea %>%
left_join(trn_road_car_loadfactor) %>%
mutate(loadFactor = ifelse(is.na(loadfactor), loadFactor, loadfactor)) %>%
select(region, supplysector, tranSubsector, stub.technology, year, loadFactor) -> L254.StubTranTechLoadFactor_trn_korea
####### Written by Jeon ####### End
process_korea_to_states(L254.StubTranTechCost) -> L254.StubTranTechCost_trn_korea
L254.StubTranTechCoef %>%
mutate(coefficient = round(coefficient, digits = gcamkorea.DIGITS_TRNUSA_DEFAULT)) %>%
process_korea_to_states ->
L254.StubTranTechCoef_trn_korea
####### Written by Jeon ####### Start
## Replace default fuel intensity with fuel intensity calculated from TS and ECS(KEEI)
trn_road_car_intensity %>%
rename(tranSubsector = GCAM_mode) -> trn_road_car_intensity
L254.StubTranTechCoef_trn_korea %>%
left_join(trn_road_car_intensity, by = c("tranSubsector", "stub.technology", "region")) %>%
mutate(coefficient = ifelse(is.na(intensity), coefficient, intensity)) %>%
select(region, supplysector, tranSubsector, stub.technology, year.x, minicam.energy.input, coefficient, market.name) %>%
rename(year = year.x)-> L254.StubTranTechCoef_trn_korea
####### Written by Jeon ####### End
process_korea_to_states(L254.PerCapitaBased_trn) -> L254.PerCapitaBased_trn_korea
process_korea_to_states(L254.PriceElasticity_trn) -> L254.PriceElasticity_trn_korea
process_korea_to_states(L254.IncomeElasticity_trn) -> L254.IncomeElasticity_trn_korea
## Written by Jeon ## Start
# Make historical tredns plausible
L254.IncomeElasticity_trn_korea %>%
left_join(kor_trn_inc_elasticity_for_realitycheck, by = c("region", "energy.final.demand", "year")) %>%
mutate( income.elasticity = ifelse( is.na(income.elasticity.y),
income.elasticity.x,
income.elasticity.y)) %>%
select(region, energy.final.demand, year, income.elasticity) -> L254.IncomeElasticity_trn_korea
####### Written by Jeon ####### End
# Calibration
# L254.StubTranTechCalInput_trn_korea: calibrated energy consumption by all technologies
L154.in_EJ_korea_trn_m_sz_tech_F %>%
filter(year %in% MODEL_BASE_YEARS) %>%
mutate(calibrated.value = round(value, digits = energy.DIGITS_CALOUTPUT),
region = state) %>%
left_join_error_no_match(select(UCD_techs, UCD_sector, mode, size.class, UCD_technology, UCD_fuel,
supplysector, tranSubsector, stub.technology = tranTechnology, minicam.energy.input),
by = c("UCD_sector", "mode", "size.class", "UCD_technology", "UCD_fuel")) %>%
select(LEVEL2_DATA_NAMES[["StubTranTech"]], year, minicam.energy.input, calibrated.value) ->
L254.StubTranTechCalInput_trn_korea
# NOTE: NEED TO WRITE THIS OUT FOR ALL TECHNOLOGIES, NOT JUST THOSE THAT EXIST IN SOME BASE YEARS.
# Model may make up calibration values otherwise.
L254.StubTranTechCoef_trn_korea %>%
filter(year %in% MODEL_BASE_YEARS) %>%
select(names(.)[names(.) %in% LEVEL2_DATA_NAMES[["StubTranTechCalInput"]]]) %>%
left_join(L254.StubTranTechCalInput_trn_korea,
by = c("region", "supplysector", "tranSubsector", "stub.technology", "year", "minicam.energy.input")) %>%
# Set calibration values as zero for technolgies that do not exist in some base years
replace_na(list(calibrated.value = 0)) %>%
mutate(share.weight.year = year,
# Create the needed variables to use the function set_subsector_shrwt
subsector = tranSubsector, calOutputValue = calibrated.value) %>%
set_subsector_shrwt() %>%
mutate(tech.share.weight = if_else(calibrated.value > 0, 1, 0)) %>%
select(LEVEL2_DATA_NAMES[["StubTranTechCalInput"]]) ->
L254.StubTranTechCalInput_trn_korea
# Non-motorized technologies
# L254.StubTranTechProd_nonmotor_korea: service output of non-motorized transportation technologies
#minyoung
L154.out_mpkm_korea_trn_nonmotor_Yh %>%
filter(year %in% MODEL_BASE_YEARS) %>%
mutate(calOutputValue = round(value, digits = energy.DIGITS_MPKM),
region = state, tranSubsector = mode) %>%
left_join_error_no_match(A54.globaltech_nonmotor, by = "tranSubsector") %>%
mutate(stub.technology = technology) %>%
# There is no need to match shareweights to the calOutputValue because no region should ever have a 0 here
select(LEVEL2_DATA_NAMES[["StubTranTech"]], year, calOutputValue) ->
L254.StubTranTechProd_nonmotor_korea
# L254.StubTranTechCalInput_passthru_trn_korea: calibrated input of passthrough technologies
# trn_pass, trn_pass_road, trn_pass_road_LDV, trn_freight
# The transportation structure is heavily nested.
# The GCAM structure of sector/subsector/technology only allows two levels of nesting within any sector,
# but a technology of one sector (e.g., trn_pass) can consume the output of another "sector" (e.g., trn_pass_road)
# that is really just used to represent lower nesting levels of that first, or parent, sector. In the
# transportation sector, each lower-level nesting "sector" is named by appending a string to the parent sector.
# So, trn_pass contains trn_pass_road which has trn_pass_road_LDV which has trn_pass_road_LDV_4W. Each of the links
# between any two of those sectors is done with a pass-through technology within the parent sector that consumes
# the output of the child sector. The technology is called a "pass-through" because it (generally) only consumes
# the output of the child "sector" without making any changes to it. There's an additional complication in the
# transportation sector: the pass-through technologies are normal, standard GCAM technologies, not "tranTechnologies"
# which have different parameters read in, and perform a bunch of hard-wired unit conversions between inputs and outputs.
# First, need to calculate the service output for all tranTechnologies
# calInput * loadFactor * unit_conversion / (coef * unit conversion)
L254.StubTranTechCalInput_trn_korea %>%
left_join_error_no_match(L254.StubTranTechLoadFactor_trn_korea,
by = c("region", "supplysector", "tranSubsector", "stub.technology", "year")) %>%
left_join_error_no_match(L254.StubTranTechCoef_trn_korea,
by = c("region", "supplysector", "tranSubsector", "stub.technology", "year", "minicam.energy.input")) %>%
mutate(output = round(calibrated.value * loadFactor * CONV_EJ_GJ / (coefficient * CONV_BTU_KJ),
digits = gcamkorea.DIGITS_TRNUSA_DEFAULT)) ->
L254.StubTranTechOutput_korea
# The next step is to calculate the aggregated outputs by supplysector
# Outputs of certain supplysectors are inputs for the passthrough technologies
L254.StubTranTechOutput_korea %>%
group_by(region, year, supplysector) %>%
summarise(output_agg = sum(output)) %>%
ungroup() ->
L254.StubTranTechOutput_korea_agg
# Write all possible pass-through technologies to all regions
A54.globaltech_passthru %>%
repeat_add_columns(tibble(year = MODEL_BASE_YEARS)) %>%
write_to_all_states(names = c(names(.), "region"), region_states=gcamkorea.STATES) %>%
select(region, supplysector, tranSubsector, stub.technology = technology, year, minicam.energy.input) %>%
# Subset only the passthrough technologies that are applicable in each region
semi_join(L254.StubTranTech_passthru_trn_korea,
by = c("region", "supplysector", "tranSubsector", "stub.technology")) %>%
# Match in outputs of supplysectors that are inputs for the passthrough technologies
left_join(L254.StubTranTechOutput_korea_agg,
by = c("region", "year", "minicam.energy.input" = "supplysector")) %>%
# Some of the technologies are sub-totals, assign zero value now, will be calculated below
replace_na(list(output_agg = 0)) %>%
# Arrange input sectors so that sub-total sector is behind the subsectors
arrange(desc(minicam.energy.input)) %>%
group_by(region, year) %>%
# Calculate the cumulative for sub-total sector
mutate(output_cum = cumsum(output_agg)) %>%
ungroup() ->
L254.StubTranTechCalInput_passthru_trn_korea_cum
# Prepare a list of the supplysector in the passthrough input table to filter the sub-total sectors
LIST_supplysector <- unique(L254.StubTranTechCalInput_passthru_trn_korea_cum$supplysector)
L254.StubTranTechCalInput_passthru_trn_korea_cum %>%
# Use the cumulative value for sub-total sectors
mutate(calibrated.value = if_else(minicam.energy.input %in% LIST_supplysector,
output_cum, output_agg)) %>%
mutate(share.weight.year = year,
subs.share.weight = if_else(calibrated.value > 0, 1, 0),
tech.share.weight = if_else(calibrated.value > 0, 1, 0)) %>%
select(LEVEL2_DATA_NAMES[["StubTranTechCalInput"]]) ->
L254.StubTranTechCalInput_passthru_trn_korea
# L254.BaseService_trn_korea: base-year service output of transportation final demand
L254.StubTranTechOutput_korea %>%
select(LEVEL2_DATA_NAMES[["StubTranTech"]], year, base.service = output) %>%
bind_rows(L254.StubTranTechProd_nonmotor_korea %>%
select(LEVEL2_DATA_NAMES[["StubTranTech"]], year, base.service = calOutputValue)) %>%
left_join_error_no_match(select(A54.sector, supplysector, energy.final.demand), by = "supplysector") %>%
group_by(region, energy.final.demand, year) %>%
summarise(base.service = sum(base.service)) %>%
ungroup ->
L254.BaseService_trn_korea
# #minyoung
# L254.StubTranTechOutput_korea %>%
# select(LEVEL2_DATA_NAMES[["StubTranTech"]], year, base.service = output) %>%
# left_join_error_no_match(select(A54.sector, supplysector, energy.final.demand), by = "supplysector") %>%
# group_by(region, energy.final.demand, year) %>%
# summarise(base.service = sum(base.service)) %>%
# ungroup ->
# L254.BaseService_trn_korea
# Produce outputs
L254.DeleteSupplysector_koreaTrn %>%
add_title("Delect transportation supply sectors of the full USA region") %>%
add_units("NA") %>%
add_comments("Delect transportation supply sectors of the full USA region") %>%
add_legacy_name("L254.DeleteSupplysector_koreaTrn") %>%
add_precursors("L254.Supplysector_trn") ->
L254.DeleteSupplysector_koreaTrn
L254.DeleteFinalDemand_koreaTrn %>%
add_title("Delete energy final demand sectors of the full USA region") %>%
add_units("NA") %>%
add_comments("Delete energy final demand sectors of the full USA region") %>%
add_legacy_name("L254.DeleteFinalDemand_koreaTrn") %>%
add_precursors("L254.PerCapitaBased_trn") ->
L254.DeleteFinalDemand_koreaTrn
L254.Supplysector_trn_korea %>%
add_title("Supply sector information for transportation sector in the US states") %>%
add_units("Unitless") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.Supplysector_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.Supplysector_trn") ->
L254.Supplysector_trn_korea
L254.FinalEnergyKeyword_trn_korea %>%
add_title("Supply sector final energy keywords for transportation sector in the US states") %>%
add_units("NA") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.FinalEnergyKeyword_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.FinalEnergyKeyword_trn") ->
L254.FinalEnergyKeyword_trn_korea
L254.tranSubsectorLogit_trn_korea %>%
add_title("Subsector logit exponents of transportation sector in the US states") %>%
add_units("Unitless") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.tranSubsectorLogit_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.tranSubsectorLogit") ->
L254.tranSubsectorLogit_trn_korea
L254.tranSubsectorShrwtFllt_trn_korea %>%
add_title("Subsector shareweights of transportation sector in the US states") %>%
add_units("Unitless") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.tranSubsectorShrwtFllt_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.tranSubsectorShrwtFllt") ->
L254.tranSubsectorShrwtFllt_trn_korea
L254.tranSubsectorInterp_trn_korea %>%
add_title("Temporal subsector shareweight interpolation of transportation sector in the US states") %>%
add_units("Unitless") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.tranSubsectorInterp_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.tranSubsectorInterp") ->
L254.tranSubsectorInterp_trn_korea
L254.tranSubsectorSpeed_trn_korea %>%
add_title("Speeds of transportation modes (not including pass-through sectors) in the US states") %>%
add_units("km / hr") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.tranSubsectorSpeed_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.tranSubsectorSpeed") ->
L254.tranSubsectorSpeed_trn_korea
L254.tranSubsectorSpeed_passthru_trn_korea %>%
add_title("Speeds of pass-through transportation subsectors in the US states") %>%
add_units("km / hr") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.tranSubsectorSpeed_passthru_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.tranSubsectorSpeed_passthru") ->
L254.tranSubsectorSpeed_passthru_trn_korea
L254.tranSubsectorSpeed_noVOTT_trn_korea %>%
add_title("Speeds of transportation subsectors whose time value is not considered in the US states") %>%
add_units("km / hr") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.tranSubsectorSpeed_noVOTT_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.tranSubsectorSpeed_noVOTT") ->
L254.tranSubsectorSpeed_noVOTT_trn_korea
L254.tranSubsectorSpeed_nonmotor_trn_korea %>%
add_title("Speeds of non-motorized transportation subsectors in the US states") %>%
add_units("km / hr") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.tranSubsectorSpeed_nonmotor_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.tranSubsectorSpeed_nonmotor") ->
L254.tranSubsectorSpeed_nonmotor_trn_korea
L254.tranSubsectorVOTT_trn_korea %>%
add_title("descriptive title of data") %>%
add_units("units") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.tranSubsectorVOTT_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.tranSubsectorVOTT") ->
L254.tranSubsectorVOTT_trn_korea
L254.tranSubsectorFuelPref_trn_korea %>%
add_title("Subsector (fuel) preferences elasticity that are tied to GDP in the US states") %>%
add_units("Unitless") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_comments("Fuel preferences are unrelated to time value") %>%
add_legacy_name("L254.tranSubsectorFuelPref_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.tranSubsectorFuelPref") ->
L254.tranSubsectorFuelPref_trn_korea
L254.StubTranTech_trn_korea %>%
add_title("Transportation stub technologies in the US states") %>%
add_units("NA") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.StubTranTech_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.StubTranTech") ->
L254.StubTranTech_trn_korea
L254.StubTranTech_passthru_trn_korea %>%
add_title("Transportation stub technologies for passthrough sectors in the US states") %>%
add_units("NA") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name(" L254.StubTranTech_passthru_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.StubTech_passthru") ->
L254.StubTranTech_passthru_trn_korea
L254.StubTranTech_nonmotor_trn_korea %>%
add_title("Transportation stub technologies for non-motorized subsectors in the US states") %>%
add_units("NA") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.StubTranTech_nonmotor_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.StubTech_nonmotor") ->
L254.StubTranTech_nonmotor_trn_korea
L254.StubTranTechLoadFactor_trn_korea %>%
add_title("Load factors of transportation stub technologies in the US states") %>%
add_units("person/vehicle and tonnes/vehicle") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.StubTranTechLoadFactor_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.StubTranTechLoadFactor",
"gcam-korea/trn_road_car_loadfactor") ->
L254.StubTranTechLoadFactor_trn_korea
L254.StubTranTechCost_trn_korea %>%
add_title("Costs of transportation stub technologies in the US states") %>%
add_units("$1990USD / vkm") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.StubTranTechCost_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.StubTranTechCost") ->
L254.StubTranTechCost_trn_korea
L254.StubTranTechCoef_trn_korea %>%
add_title("Coefficients of transportation stub technologies in the US states") %>%
add_units("BTU / vkm") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_comments("Re-set electricity consumed at the state markets") %>%
add_legacy_name("L254.StubTranTechCoef_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.StubTranTechCoef",
"gcam-korea/trn_road_car_intensity") ->
L254.StubTranTechCoef_trn_korea
L254.PerCapitaBased_trn_korea %>%
add_title("Per-capita based flag for transportation final demand in the US states") %>%
add_units("NA") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.PerCapitaBased_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.PerCapitaBased_trn") ->
L254.PerCapitaBased_trn_korea
L254.PriceElasticity_trn_korea %>%
add_title("Price elasticity of transportation final demand in the US states") %>%
add_units("Unitless") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.PriceElasticity_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.PriceElasticity_trn") ->
L254.PriceElasticity_trn_korea
L254.IncomeElasticity_trn_korea %>%
add_title("Income elasticity of transportation final demand in the US states") %>%
add_units("Unitless") %>%
add_comments("The same USA region values are repeated for each state") %>%
add_legacy_name("L254.IncomeElasticity_trn_korea") %>%
add_precursors("gcam-korea/states_subregions",
"L254.IncomeElasticity_trn",
"gcam-korea/kor_trn_inc_elasticity_for_realitycheck") ->
L254.IncomeElasticity_trn_korea
L254.StubTranTechCalInput_trn_korea %>%
add_title("Calibrated energy consumption by all transportation stub technologies in the US states") %>%
add_units("EJ") %>%
add_comments("Set calibration values for those technologies that do not exist in some base years as zero") %>%
add_legacy_name("L254.StubTranTechCalInput_trn_korea") %>%
same_precursors_as("L254.StubTranTechCoef_trn_korea") %>%
add_precursors("L154.in_EJ_korea_trn_m_sz_tech_F",
"energy/mappings/UCD_techs") ->
L254.StubTranTechCalInput_trn_korea
L254.StubTranTechProd_nonmotor_korea %>%
add_title("Calibrated service output of non-motorized transportation technologies in the US states") %>%
add_units("Million pass-km") %>%
add_comments("Not match shareweights to the calOutputValue because no region should ever have a zero here") %>%
add_legacy_name("L254.StubTranTechProd_nonmotor_korea") %>%
add_precursors("L154.out_mpkm_korea_trn_nonmotor_Yh",
"energy/A54.globaltech_nonmotor") ->
L254.StubTranTechProd_nonmotor_korea
L254.StubTranTechCalInput_passthru_trn_korea %>%
add_title("Calibrated energy consumption of transportation passthrough technologies in the US states") %>%
add_units("EJ") %>%
add_comments("Use outputs of the supplysectors that are inputs for passthrough technologies") %>%
add_comments("Outputs of all motorized technologies are calculated as calInput * loadFactor / coefficient") %>%
add_legacy_name("L254.StubTranTechCalInput_passthru_trn_korea") %>%
same_precursors_as("L254.StubTranTechCalInput_trn_korea") %>%
same_precursors_as("L254.StubTranTechLoadFactor_trn_korea") %>%
same_precursors_as("L254.StubTranTechCoef_trn_korea") %>%
same_precursors_as("L254.StubTranTech_passthru_trn_korea") %>%
add_precursors("energy/A54.globaltech_passthru") ->
L254.StubTranTechCalInput_passthru_trn_korea
L254.BaseService_trn_korea %>%
add_title("Base-year service output of transportation final demand") %>%
add_units("Million pass-km and million ton-km") %>%
add_comments("Service outputs of all motorized technologies are calculated as calInput * loadFactor / coefficient") %>%
add_comments("Combine with service output of non-motorized transportation technologies") %>%
add_legacy_name("L254.BaseService_trn_korea") %>%
same_precursors_as("L254.StubTranTechCalInput_trn_korea") %>%
same_precursors_as("L254.StubTranTechLoadFactor_trn_korea") %>%
same_precursors_as("L254.StubTranTechCoef_trn_korea") %>%
same_precursors_as("L254.StubTranTechProd_nonmotor_korea") %>%
add_precursors("energy/A54.sector") ->
L254.BaseService_trn_korea
return_data(L254.DeleteSupplysector_koreaTrn, L254.DeleteFinalDemand_koreaTrn,
L254.Supplysector_trn_korea,
L254.FinalEnergyKeyword_trn_korea,
L254.tranSubsectorLogit_trn_korea,
L254.tranSubsectorShrwtFllt_trn_korea,
L254.tranSubsectorInterp_trn_korea,
L254.tranSubsectorSpeed_trn_korea,
L254.tranSubsectorSpeed_passthru_trn_korea,
L254.tranSubsectorSpeed_noVOTT_trn_korea,
L254.tranSubsectorSpeed_nonmotor_trn_korea,
L254.tranSubsectorVOTT_trn_korea,
L254.tranSubsectorFuelPref_trn_korea,
L254.StubTranTech_trn_korea,
L254.StubTranTech_passthru_trn_korea,
L254.StubTranTech_nonmotor_trn_korea,
L254.StubTranTechLoadFactor_trn_korea,
L254.StubTranTechCost_trn_korea,
L254.StubTranTechCoef_trn_korea,
L254.PerCapitaBased_trn_korea,
L254.PriceElasticity_trn_korea,
L254.IncomeElasticity_trn_korea,
L254.StubTranTechCalInput_trn_korea,
L254.StubTranTechProd_nonmotor_korea,
L254.StubTranTechCalInput_passthru_trn_korea, L254.BaseService_trn_korea)
} else {
stop("Unknown command")
}
}
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.