R/zchunk_L222.en_transformation_korea.R

In rohmin9122/gcam-korea-release: Build the GCAM-Korea Input Datasets

#' module_gcam.korea_L222.en_transformation_korea
#'
#' Prepare the assumptions and calibrated outputs for energy transformation supplysectors, subsectors, and technologies specific to subregion's sectors and/or 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{L222.DeleteStubTech_koreaEn}, \code{L222.PassThroughSector_koreaEn}, \code{L222.Tech_koreaEn}, \code{L222.TechShrwt_koreaEn}, \code{L222.TechInterp_koreaEn}, \code{L222.TechShrwt_koreaEn}, \code{L222.TechCoef_koreaEn}, \code{L222.Production_koreaRefining}, \code{L222.SectorLogitTables_USA[[ curr_table ]]$data}, \code{L222.Supplysector_en_korea}, \code{L222.SubsectorShrwtFllt_en_korea}, \code{L222.StubTechProd_refining_korea}, \code{L222.StubTechMarket_en_korea}, \code{L222.CarbonCoef_en_korea}. The corresponding file in the
#' original data system was \code{L222.en_transformation_USA.R} (gcam-korea level2).
#' @details This chunk sets up the subregion's energy transformation technology databases as well as writing out assumptions to all states/sectors/markets for shareweights and logits.
#' Calibrated outputs and I:O coefficients are updated from global values produced by \code{\link{module_energy_L222.en_transformation}}.
#' @importFrom assertthat assert_that
#' @importFrom dplyr filter mutate select
#' @importFrom tidyr gather spread
#' @author S. Jeon
module_gcam.korea_L222.en_transformation_korea <- function(command, ...) {
  if(command == driver.DECLARE_INPUTS) {
    return(c(FILE = "gcam-korea/states_subregions",
             FILE = "energy/calibrated_techs",
             "L222.Supplysector_en",
             "L222.SubsectorLogit_en",
             "L222.StubTech_en",
             "L222.StubTechCoef_refining",
             "L222.GlobalTechInterp_en",
             "L222.GlobalTechCoef_en",
             "L222.GlobalTechCost_en",
             "L222.GlobalTechShrwt_en",
             "L222.GlobalTechCapture_en",
             "L222.GlobalTechSCurve_en",
             "L122.out_EJ_korea_refining_F",
             "L202.CarbonCoef"))
  } else if(command == driver.DECLARE_OUTPUTS) {
    return(c("L222.DeleteStubTech_koreaEn",
             "L222.PassThroughSector_koreaEn",
             "L222.Tech_koreaEn",
             "L222.TechShrwt_koreaEn",
             "L222.TechInterp_koreaEn",
             "L222.TechCoef_koreaEn",
             "L222.Production_koreaRefining",
             "L222.Supplysector_en_korea",
             "L222.SubsectorShrwtFllt_en_korea",
             "L222.StubTechProd_refining_korea",
             "L222.StubTechMarket_en_korea",
             "L222.CarbonCoef_en_korea",
             "L222.GlobalTechSCurve_en_korea",
             "L222.GlobalTechCost_en_korea",
             "L222.SubsectorLogit_en_korea",
             "L222.StubTech_en_korea",
             "L222.StubTechCoef_refining_korea",
             "L222.GlobalTechInterp_en_korea",
             "L222.GlobalTechCoef_en_korea",
             "L222.GlobalTechShrwt_en_korea",
             "L222.GlobalTechCapture_en_korea"))
  } else if(command == driver.MAKE) {

    all_data <- list(...)[[1]]

    # Load required inputs
    states_subregions <- get_data(all_data, "gcam-korea/states_subregions")
    calibrated_techs <- get_data(all_data, "energy/calibrated_techs")
    L222.Supplysector_en <- get_data(all_data, "L222.Supplysector_en")
    L222.SubsectorLogit_en <- get_data(all_data, "L222.SubsectorLogit_en")
    L222.StubTech_en <- get_data(all_data, "L222.StubTech_en")
    L222.StubTechCoef_refining <- get_data(all_data, "L222.StubTechCoef_refining")
    L222.GlobalTechInterp_en <- get_data(all_data, "L222.GlobalTechInterp_en")
    L222.GlobalTechCoef_en <- get_data(all_data, "L222.GlobalTechCoef_en")
    L222.GlobalTechCost_en <- get_data(all_data, "L222.GlobalTechCost_en")
    L222.GlobalTechShrwt_en <- get_data(all_data, "L222.GlobalTechShrwt_en")
    L222.GlobalTechCapture_en <- get_data(all_data, "L222.GlobalTechCapture_en")
    L222.GlobalTechSCurve_en <- get_data(all_data, "L222.GlobalTechSCurve_en")
    L122.out_EJ_korea_refining_F <- get_data(all_data, "L122.out_EJ_korea_refining_F")
    L202.CarbonCoef <- get_data(all_data, "L202.CarbonCoef")

    # silence check package notes
    logit.year.fillout <- year <- from.year <- to.year <- region <- supplysector <- subsector <-
      technology <- sector.name <- subsector.name <- sector <- state <- fuel <- value <- market.name <-
      trash <- calOutputValue <- minicam.energy.input <- supplysector.x <- supplysector.y <-
      calibration <- grid_region <- stub.technology <- key <- share.weight <- NULL

    # Correct some of the inputs
    L222.Supplysector_en %>%
      mutate(logit.year.fillout = as.integer(logit.year.fillout)) -> # was character
      L222.Supplysector_en

    L222.SubsectorLogit_en  %>%
      mutate(logit.year.fillout = as.integer(logit.year.fillout)) -> # was character
      L222.SubsectorLogit_en

    L222.StubTechCoef_refining %>%
      mutate(year = as.integer(year)) -> # was double
      L222.StubTechCoef_refining

    L222.GlobalTechInterp_en %>%
      mutate(from.year = as.integer(from.year), # was character
             to.year = as.integer(to.year)) ->
      L222.GlobalTechInterp_en

    L222.GlobalTechCoef_en %>%
      mutate(year  = as.integer(year)) -> # was character
      L222.GlobalTechCoef_en

    L222.GlobalTechCost_en %>%
      mutate(year = as.integer(year)) -> # was double
      L222.GlobalTechCost_en

    L222.GlobalTechCapture_en %>%
      mutate(year = as.integer(year)) -> # was character
      L222.GlobalTechCapture_en

    L222.GlobalTechSCurve_en %>%
      mutate(year = as.integer(year)) -> # was character
      L222.GlobalTechSCurve_en


    # Some helpful functions:
    #
    # global_energy_to_korea_nonGlobalTech - takes global energy inputs for non global tech
    # from L222.en_transformation.R and processes for use in USA
    global_energy_to_korea_nonGlobalTech <- function(data) {
      data %>%
        filter(region == gcamkorea.REGION,
               supplysector %in% gcamkorea.SECTOR_EN_NAMES) %>%
        write_to_all_states(names = c(names(data), "region"), region_states=gcamkorea.STATES) %>%
        filter((subsector == "oil refining" & region %in% oil_refining_states) |
                 subsector != "oil refining") %>%
        mutate(supplysector = subsector)
    } # global_energy_to_korea_nonGlobalTech

    # global_energy_to_korea_GlobalTech - takes global energy inputs for global tech
    # from L222.en_transformation.R and processes for use in USA
    global_energy_to_korea_GlobalTech <- function(data) {
      data %>%
        filter(sector.name %in% gcamkorea.SECTOR_EN_NAMES) %>%
        mutate(sector.name = subsector.name)
    } # global_energy_to_korea_GlobalTech

    # Oil refining sectors are only created in states where the production is > 0 in the historical period.
    # Collect these states. Other techs are available everywhere
    L122.out_EJ_korea_refining_F %>%
      filter(sector == "oil refining",
             year %in% HISTORICAL_YEARS) %>%
      group_by(state, sector, fuel) %>%
      summarise(value = sum(value)) %>%
      ungroup %>%
      filter(value > 0) %>%
      pull(state) ->
      oil_refining_states

    # L222.DeleteStubTech_koreaEn: remove existing stub technologies in the USA region.
    # The supplysector and subsector structure in the sectors defined in gcamkorea.SECTOR_EN_NAMES are retained
    L222.StubTech_en %>%
      filter(region == gcamkorea.REGION,
             supplysector %in% gcamkorea.SECTOR_EN_NAMES) ->
      L222.DeleteStubTech_koreaEn

    # L222.Tech_koreaEn: Just the technology pass-throughs used to set the proper node name, USA region
    L222.SubsectorLogit_en %>%
      select(region, supplysector, subsector) %>%
      filter(region == gcamkorea.REGION,
             supplysector %in% gcamkorea.SECTOR_EN_NAMES) %>%
      repeat_add_columns(tibble(state = gcamkorea.STATES)) %>%
      filter((subsector == "oil refining" & state %in% oil_refining_states) |
               subsector != "oil refining") %>%
      mutate(technology = paste(state, subsector, sep = gcamkorea.STATE_SUBSECTOR_DELIMITER)) ->
      L222.Tech_koreaEn

    # save some of this information for the PassThroughSector information
    # L222.PassThroughSector_koreaEn: PassThroughSector information to send vintaging info from states to USA.
    L222.Tech_koreaEn %>%
      select(state, subsector, supplysector, region) %>%
      rename(marginal.revenue.market = region,
             region = state,
             pass.through.sector = subsector,
             marginal.revenue.sector = supplysector) ->
      L222.PassThroughSector_koreaEn

    # select only relevant columns for L222.Tech_koreaEn, particularly dropping state
    L222.Tech_koreaEn %>%
      select(one_of(LEVEL2_DATA_NAMES[["Tech"]])) ->
      L222.Tech_koreaEn

    # L222.TechInterp_koreaEn: technology shareweights, USA region
    # Technology interpolation only applies to calibrated technologies.
    # For biomass liquids, allow state shares to shift over time
    # (future techs are different than present techs).
    L222.Tech_koreaEn %>%
      filter(subsector %in% c("oil refining", "biomass liquids")) %>%
      mutate(apply.to = "share-weight",
             from.year = max(MODEL_BASE_YEARS),
             to.year = max(MODEL_YEARS),
             interpolation.function = if_else(subsector == "biomass liquids", "s-curve", "fixed")) ->
      L222.TechInterp_koreaEn

    # L222.TechShrwt_koreaEn: technology shareweights in each year, USA region
    # Default the base year shareweights to 0. This will be over-ridden in calibration
    # Default the future year shareweights to 1.
    L222.Tech_koreaEn %>%
      repeat_add_columns(tibble(year = MODEL_BASE_YEARS)) %>%
      mutate(share.weight = 0) ->
      tmp

    L222.Tech_koreaEn %>%
      repeat_add_columns(tibble(year = MODEL_FUTURE_YEARS)) %>%
      mutate(share.weight = 1) %>%
      bind_rows(tmp) ->
      L222.TechShrwt_koreaEn

    # L222.TechCoef_koreaEn: technology coefficients and market names, USA region
    L222.TechShrwt_koreaEn %>%
      select(one_of(LEVEL2_DATA_NAMES[["TechYr"]])) %>%
      mutate(minicam.energy.input = subsector,
             coefficient = 1,
             market.name = technology) %>%
      separate(market.name, c("market.name", "trash"), extra = "merge", sep = gcamkorea.STATE_SUBSECTOR_DELIMITER) %>%
      select(-trash) ->
      L222.TechCoef_koreaEn

    # L222.Production_koreaRefining: calibrated refinery production in USA (consuming output of states)
    # Aggregated to the supplysector/subsector/technology level
      L122.out_EJ_korea_refining_F %>%
        filter(year %in% MODEL_BASE_YEARS) %>%
        rename(calOutputValue = value) %>%
        mutate(calOutputvalue = round(calOutputValue, gcamkorea.DIGITS_CALOUTPUT ),
               region = gcamkorea.REGION) %>%
        left_join_error_no_match(distinct(select(calibrated_techs, sector, supplysector, subsector)), by = "sector") %>%
        mutate(technology = paste(state, subsector, sep = gcamkorea.STATE_SUBSECTOR_DELIMITER),
               minicam.energy.input = subsector) %>%
        filter((subsector == "oil refining" & state %in% oil_refining_states) |
                 subsector != "oil refining") %>%
        # Aggregate
        group_by(region, supplysector, subsector, technology, minicam.energy.input, year) %>%
        summarise(calOutputValue = sum(calOutputValue)) %>%
        ungroup %>%
        mutate(share.weight.year = year) %>%
        set_subsector_shrwt %>%
        # The following line is equivalent to (but slightly faster than): mutate(tech.share.weight = if_else(calOutputValue == 0, 0, 1)) %>%
        mutate(tech.share.weight = abs(sign(calOutputValue))) %>%
        select(one_of(LEVEL2_DATA_NAMES[["Production"]])) ->
        L222.Production_koreaRefining

      # Process energy files from L222.en_transformation.R for use in the USA,
      # slightly differently processing for global tech vs not inputs
      L222.SubsectorLogit_en_korea      <- global_energy_to_korea_nonGlobalTech(L222.SubsectorLogit_en)
      L222.StubTech_en_korea            <- global_energy_to_korea_nonGlobalTech(L222.StubTech_en)
      L222.StubTechCoef_refining_korea  <- global_energy_to_korea_nonGlobalTech(L222.StubTechCoef_refining)
      L222.GlobalTechInterp_en_korea    <- global_energy_to_korea_GlobalTech(L222.GlobalTechInterp_en)
      L222.GlobalTechCoef_en_korea      <- global_energy_to_korea_GlobalTech(L222.GlobalTechCoef_en)
      L222.GlobalTechCost_en_korea      <- global_energy_to_korea_GlobalTech(L222.GlobalTechCost_en)
      L222.GlobalTechShrwt_en_korea     <- global_energy_to_korea_GlobalTech(L222.GlobalTechShrwt_en)
      L222.GlobalTechCapture_en_korea   <- global_energy_to_korea_GlobalTech(L222.GlobalTechCapture_en)
      L222.GlobalTechSCurve_en_korea    <- global_energy_to_korea_GlobalTech(L222.GlobalTechSCurve_en)


      ### The same processing for optional/currently NULL inputs
      # L222.GlobalTechShutdownProfit_en_korea  <- global_energy_to_korea_GlobalTech(L222.GlobalTechShutdownProfit_en)
      # L222.GlobalTechShutdown_en_korea        <- global_energy_to_korea_GlobalTech(L222.GlobalTechShutdown_en)
      # L222.GlobalTechSCurveProfit_en_korea    <- global_energy_to_korea_GlobalTech(L222.GlobalTechSCurveProfit_en)
      # L222.GlobalTechLifetimeProfit_en_korea  <- global_energy_to_korea_GlobalTech(L222.GlobalTechLifetimeProfit_en)
      # L222.GlobalTechLifetime_en_korea        <- global_energy_to_korea_GlobalTech(L222.GlobalTechLifetime_en)

      # TODO: figure out a better strategy.  We need to have at least one technology be available in the final
      # calibration year so we can get a base cost for the absolute cost logit.  Having a share weight of zero
      # at the subsector is sufficient then to ensure we get no production in the calibration years
      L222.GlobalTechShrwt_en_korea %>%
        mutate(share.weight = if_else(technology == "coal to liquids" & year == max(MODEL_BASE_YEARS), 1.0, share.weight),
               share.weight = if_else(technology == "gas to liquids" & year == max(MODEL_BASE_YEARS), 1.0, share.weight)) ->
        L222.GlobalTechShrwt_en_korea

      # L222.Supplysector_en_korea: Supplysector information, replace name of supplysector with the subsector names
      L222.SubsectorLogit_en_korea %>%
        select(one_of(LEVEL2_DATA_NAMES[["Subsector"]])) %>%
        left_join_error_no_match(distinct(select(L222.SubsectorLogit_en, supplysector, subsector)),
                                 by = "subsector") %>%
        rename(supplysector = supplysector.x,
               old_supplysector = supplysector.y) %>%
        left_join_error_no_match(distinct(select(L222.Supplysector_en, -region)),
                                 by = c("old_supplysector" = "supplysector")) %>%
        select(one_of(LEVEL2_DATA_NAMES[["Supplysector"]])) ->
        L222.Supplysector_en_korea

      # L222.Supplysector_en_korea_logit.type - Note there is no competition here so just use the default logit type
      L222.Supplysector_en_korea %>%
        mutate(logit.type = gcamusa.DEFAULT_LOGIT_TYPE) ->
        L222.Supplysector_en_korea_logit.type

      # L222.SubsectorShrwtFllt_en_korea: Subsector shareweights, there is no competition here, so just fill out with 1s
      # (will be over-ridden by base year calibration where necessary)
      L222.SubsectorLogit_en_korea %>%
        select(one_of(LEVEL2_DATA_NAMES[["Subsector"]])) %>%
        mutate(year = min(MODEL_YEARS),
               share.weight = gcamusa.DEFAULT_SHAREWEIGHT) ->
        L222.SubsectorShrwtFllt_en_korea

      # L222.StubTechProd_refining_korea: calibrated fuel production by state.
      # Only take the tech IDs where the calibration is identified as output.
      #
      # Step 1, process the table of calibrated_techs to only include calibration=output and relevant columns
      calibrated_techs %>%
        filter(calibration == "output") %>%
        select(sector, supplysector, subsector, technology) %>%
        distinct ->
        calibrated_techs_tmp

      # Step 2, process L122.out_EJ_korea_refining_F, joining the processed table of calibrated_techs from step 1,
      # to create L222.StubTechProd_refining_korea. Note the supplysector is the same as the subsector within the states.
      L122.out_EJ_korea_refining_F %>%
        filter(year %in% MODEL_BASE_YEARS) %>%
        rename(region = state,
               calOutputValue = value) %>%
        mutate(calOutputValue = round(calOutputValue, gcamkorea.DIGITS_CALOUTPUT )) %>%
        left_join_error_no_match(calibrated_techs_tmp, by = "sector") %>%
        mutate(supplysector = subsector,
               stub.technology = technology,
               share.weight.year = year) %>%
        set_subsector_shrwt() %>%
        mutate(tech.share.weight = if_else(calOutputValue > 0, 1, 0)) %>%
        select(one_of(LEVEL2_DATA_NAMES[["StubTechProd"]])) %>%
        filter((subsector == "oil refining" & region %in% oil_refining_states) |
                 subsector != "oil refining") ->
        L222.StubTechProd_refining_korea

      # L222.StubTechMarket_en_korea: market names of inputs to state refining sectors
      L222.GlobalTechCoef_en_korea %>%
        select(one_of(LEVEL2_DATA_NAMES[["GlobalTechInput"]])) %>%
        write_to_all_states(names = c(LEVEL2_DATA_NAMES[["GlobalTechInput"]], "region"), region_states=gcamkorea.STATES) %>%
        rename(supplysector = sector.name,
               subsector = subsector.name,
               stub.technology = technology) %>%
        mutate(market.name = gcamkorea.REGION) ->
        L222.StubTechMarket_en_korea

      # If designated, switch fuel market names to the regional markets
      if(gcamkorea.USE_REGIONAL_FUEL_MARKETS) {
        L222.StubTechMarket_en_korea %>%
          select(-market.name) %>%
          filter(minicam.energy.input %in% gcamkorea.REGIONAL_FUEL_MARKETS) %>%
          left_join_error_no_match(distinct(select(states_subregions, state, grid_region)), by = c("region" = "state")) %>%
          rename(market.name = grid_region) ->
          tmp

        L222.StubTechMarket_en_korea %>%
          filter(!(minicam.energy.input %in% gcamkorea.REGIONAL_FUEL_MARKETS)) %>%
          bind_rows(tmp) ->
          L222.StubTechMarket_en_korea
      }

      # Finish L222.StubTechMarket_en_korea by Setting electricity to the state markets
      L222.StubTechMarket_en_korea %>%
        filter(minicam.energy.input %in% gcamkorea.ELECT_TD_SECTORS) %>%
        mutate(market.name = region) ->
        tmp

      # create a key for filtering
      L222.StubTech_en_korea %>%
        select(supplysector, subsector, stub.technology) %>%
        unite(key, supplysector, subsector, stub.technology, sep = "~") %>%
        distinct ->
        L222.StubTech_en_korea_key

      L222.StubTechMarket_en_korea %>%
        filter(!(minicam.energy.input %in% gcamkorea.ELECT_TD_SECTORS)) %>%
        bind_rows(tmp) %>%
        select(one_of(LEVEL2_DATA_NAMES[["StubTechMarket"]])) %>%
        unite(key, supplysector, subsector, stub.technology, sep = "~") %>%
        filter(key %in% L222.StubTech_en_korea_key$key) %>%
        separate(key, c("supplysector", "subsector", "stub.technology"), sep = "~") %>%
        filter((subsector == "oil refining" & region %in% oil_refining_states) |
                 subsector != "oil refining") ->
        L222.StubTechMarket_en_korea

      # L222.CarbonCoef_en_korea: energy carbon coefficients in USA
      #
      # Step 1, process L202.CarbonCoef for joining
      L202.CarbonCoef %>%
        filter(region == gcamkorea.REGION) %>%
        select(-region) %>%
        distinct ->
        L202.CarbonCoef_tmp

      # Step 2, create L222.CarbonCoef_en_korea by joining the table from step 1.
      L222.Supplysector_en_korea %>%
        select(region, supplysector) %>%
        distinct %>%
        left_join_error_no_match(distinct(select(L222.TechShrwt_koreaEn, subsector, supplysector)),
                                 by = c("supplysector" = "subsector")) %>%
        left_join_error_no_match(L202.CarbonCoef_tmp, by =  c("supplysector.y" = "PrimaryFuelCO2Coef.name")) %>%
        select(-supplysector.y) %>%
        rename(PrimaryFuelCO2Coef.name = supplysector) ->
        L222.CarbonCoef_en_korea

    # Produce outputs
    L222.DeleteStubTech_koreaEn %>%
      mutate(region = region) %>%  # strip off attributes so we can re-write title, etc.
      add_title("Removes existing stub technologies in the USA region") %>%
      add_units("NA") %>%
      add_comments("Removes existing stub technologies in the USA region from L222.StubTech_en.") %>%
      add_comments("The supplysector and subsector structure in the sectors defined in gcamkorea.SECTOR_EN_NAMES are retained")  %>%
      add_legacy_name("L222.DeleteStubTech_koreaEn") %>%
      add_precursors("L222.StubTech_en") ->
      L222.DeleteStubTech_koreaEn

    L222.PassThroughSector_koreaEn %>%
      add_title("PassThroughSector information to send vintaging info from states to USA") %>%
      add_units("NA") %>%
      add_comments("state, subsector, supplysector, and region fromj L222.Tech_koreaEn is renamed.") %>%
      add_legacy_name("L222.PassThroughSector_koreaEn") %>%
      same_precursors_as(L222.Tech_koreaEn) ->
      L222.PassThroughSector_koreaEn

    L222.Tech_koreaEn %>%
      add_title("The technology pass-throughs used to set the proper node name, USA region.") %>%
      add_units("units") %>%
      add_comments("USA supplysector and subsector information from L222.SubsectorLogit_en is") %>%
      add_comments("repeated for all US states and updated.") %>%
      add_legacy_name("L222.Tech_koreaEn") %>%
      add_precursors("L222.SubsectorLogit_en") ->
      L222.Tech_koreaEn

    L222.TechShrwt_koreaEn %>%
      add_title("Technology shareweights in each year, USA region") %>%
      add_units("NA") %>%
      add_comments("L222.Tech_koreaEn is repeated for model base year and future years and shareweights of 0 and") %>%
      add_comments("1 are added for each, respectively. Overwritten in calibration.") %>%
      add_legacy_name("L222.TechShrwt_koreaEn") %>%
      same_precursors_as(L222.Tech_koreaEn) ->
      L222.TechShrwt_koreaEn

    L222.TechInterp_koreaEn %>%
      add_title("Technology shareweights, USA region") %>%
      add_units("NA") %>%
      add_comments("Technology interpolation only applies to calibrated technologies.For biomass liquids, ") %>%
      add_comments("allow state shares to shift over time since future techs are different than present techs.") %>%
      add_legacy_name("L222.TechInterp_koreaEn")  %>%
      same_precursors_as(L222.Tech_koreaEn) ->
      L222.TechInterp_koreaEn

    L222.TechCoef_koreaEn %>%
      add_title("Technology coefficients and market names, USA region") %>%
      add_units("units") %>%
      add_comments("Data from L222.TechShrwt_koreaEn is renamed and filled out.") %>%
      add_legacy_name("L222.TechCoef_koreaEn") %>%
      same_precursors_as(L222.TechShrwt_koreaEn) ->
      L222.TechCoef_koreaEn

    L222.Production_koreaRefining %>%
      add_title("Calibrated refinery production in USA (consuming output of states)") %>%
      add_units("NA") %>%
      add_comments("L122.out_EJ_korea_refining_F is aggregated to the supplysector/subsector/technology level.") %>%
      add_legacy_name("L222.Production_koreaRefining") %>%
      add_precursors("energy/calibrated_techs",
                      "L122.out_EJ_korea_refining_F") ->
      L222.Production_koreaRefining

    L222.Supplysector_en_korea_logit.type %>%
      add_title("Supplysector information, replace name of supplysector with the subsector names") %>%
      add_units("Varies") %>%
      add_comments("L222.Supplysector_en and L222.SubsectorLogit_en is repeated and filtered for use in USA states.") %>%
      add_comments("can be multiple lines") %>%
      add_legacy_name("L222.Supplysector_en_korea") %>%
      add_precursors("L222.Supplysector_en",
                     "L222.SubsectorLogit_en") ->
      L222.Supplysector_en_korea

    L222.SubsectorShrwtFllt_en_korea %>%
      add_title("Subsector shareweights for energy in USA") %>%
      add_units("NA") %>%
      add_comments("USA energy subsector shareweights. There is no competition here, so shareweights are defaulted to 1.") %>%
      add_comments("Shareweights will be over-ridden by base year calibration.") %>%
      add_legacy_name("L222.SubsectorShrwtFllt_en_korea") %>%
      same_precursors_as(L222.SubsectorLogit_en_korea) ->
      L222.SubsectorShrwtFllt_en_korea

    L222.StubTechProd_refining_korea %>%
      add_title("USA Calibrated fuel production by state.") %>%
      add_units("varies") %>%
      add_comments("Tech IDs where the calibration is identified as output in the calibrated_techs file are") %>%
      add_comments("are used to adjust data from L122.out_EJ_korea_refining_F.") %>%
      add_legacy_name("L222.StubTechProd_refining_korea") %>%
      add_precursors("energy/calibrated_techs",
                     "L122.out_EJ_korea_refining_F") ->
      L222.StubTechProd_refining_korea

    L222.StubTechMarket_en_korea %>%
      add_title("Market names of inputs to state refining sectors") %>%
      add_units("varies") %>%
      add_comments("Data from L222.GlobalTechCoef_en is adjusted for use in US states, depending") %>%
      add_comments("on whether regional markets are used.") %>%
      add_legacy_name("L222.StubTechMarket_en_korea") %>%
      add_precursors("gcam-korea/states_subregions",
                     "L222.GlobalTechCoef_en") ->
      L222.StubTechMarket_en_korea

    L222.CarbonCoef_en_korea %>%
      add_title("Energy carbon coefficients in USA") %>%
      add_units("varies") %>%
      add_comments("Carbon coefficients from L202.CarbonCoef are updated with USA energy tech shareweights to") %>%
      add_comments("produce energy carbon coefficients in USA.") %>%
      add_legacy_name("L222.CarbonCoef_en_korea") %>%
      add_precursors("L222.SubsectorLogit_en",
                     "L202.CarbonCoef")  ->
      L222.CarbonCoef_en_korea

    L222.SubsectorLogit_en_korea %>%
      add_title("Subsector logit competition info for USA energy states and sectors") %>%
      add_units("NA") %>%
      add_comments("Subsector logit data from L222.SubsectorLogit_en are filtered and repeated") %>%
      add_comments("for USA sectors in each state.") %>%
      add_legacy_name("L222.SubsectorLogit_en_korea") %>%
      add_precursors("L222.SubsectorLogit_en") ->
      L222.SubsectorLogit_en_korea

    L222.StubTech_en_korea %>%
      add_title("Stub technology map for USA energy states and sectors.") %>%
      add_units("NA") %>%
      add_comments("The stub technology table from L222.StubTech_en is filtered and repeated") %>%
      add_comments("for USA energy sectors in each state.") %>%
      add_legacy_name("L222.StubTech_en_korea") %>%
      add_precursors("L222.StubTech_en") ->
      L222.StubTech_en_korea

    L222.StubTechCoef_refining_korea %>%
      add_title("Refining stub tech coefficients for USA energy states and sectors") %>%
      add_units("NA") %>%
      add_comments("Coefficients for refining stub technologies in L222.StubTechCoef_refining are filtered and repeated") %>%
      add_comments("for USA energy sectors in each state.") %>%
      add_legacy_name("L222.StubTechCoef_refining_korea") %>%
      add_precursors("L222.StubTechCoef_refining") ->
      L222.StubTechCoef_refining_korea

    L222.GlobalTechSCurve_en_korea %>%
      add_title("Tech S curve parameters for USA energy sectors.") %>%
      add_units("varies") %>%
      add_comments("S curve parameters from L222.GlobalTechScurve_en are filtered for USA sectors.") %>%
      add_legacy_name("L222.GlobalTechSCurve_en_korea") %>%
      add_precursors("L222.GlobalTechSCurve_en") ->
      L222.GlobalTechSCurve_en_korea

    L222.GlobalTechCost_en_korea %>%
      add_title("Tech costs for USA energy sectors.") %>%
      add_units("varies") %>%
      add_comments("Tech cost data from L222.GlobalTechCost_en are filtered for USA sectors.") %>%
      add_legacy_name("L222.GlobalTechCost_en_korea") %>%
      add_precursors("L222.GlobalTechCost_en") ->
      L222.GlobalTechCost_en_korea

    L222.GlobalTechInterp_en_korea %>%
      add_title("Interpolation function key for USA energy sectors.") %>%
      add_units("units") %>%
      add_comments("Interpolation function key from L222.GlobalTechInterp_en are filtered for USA sectors.") %>%
      add_legacy_name("L222.GlobalTechInterp_en_korea") %>%
      add_precursors("L222.GlobalTechInterp_en") ->
      L222.GlobalTechInterp_en_korea

    L222.GlobalTechCoef_en_korea %>%
      add_title("Technology coefficients for USA energy sectors.") %>%
      add_units("NA") %>%
      add_comments("Global technology coefficients from L222.GlobalTechCoef_en are filtered for USA sectors.") %>%
      add_legacy_name("L222.GlobalTechCoef_en_korea") %>%
      add_precursors("L222.GlobalTechCoef_en") ->
      L222.GlobalTechCoef_en_korea

    L222.GlobalTechShrwt_en_korea %>%
      add_title("Technology shareweights for USA energy sectors") %>%
      add_units("NA") %>%
      add_comments("Shareweights from L222.GlobalTechShrwt_en are filtered for USA energy sectors.") %>%
      add_legacy_name("L222.GlobalTechShrwt_en_korea") %>%
      add_precursors("L222.GlobalTechShrwt_en") ->
      L222.GlobalTechShrwt_en_korea

    L222.GlobalTechCapture_en_korea %>%
      add_title("Carbon capture data for USA energy sectors") %>%
      add_units("NA") %>%
      add_comments("Carbon capture data  from L222.GlobalTechCapture_en are filtered for USA energy sectors.") %>%
      add_legacy_name("L222.GlobalTechCapture_en_korea") %>%
      add_precursors("L222.GlobalTechCapture_en") ->
      L222.GlobalTechCapture_en_korea

    return_data(L222.DeleteStubTech_koreaEn, L222.PassThroughSector_koreaEn, L222.Tech_koreaEn,
                L222.TechShrwt_koreaEn, L222.TechInterp_koreaEn, L222.TechCoef_koreaEn, L222.Production_koreaRefining,
                L222.Supplysector_en_korea, L222.SubsectorShrwtFllt_en_korea, L222.StubTechProd_refining_korea, L222.StubTechMarket_en_korea,
                L222.CarbonCoef_en_korea, L222.GlobalTechSCurve_en_korea,
                L222.GlobalTechCost_en_korea,
                L222.SubsectorLogit_en_korea,
                L222.StubTech_en_korea,
                L222.StubTechCoef_refining_korea,
                L222.GlobalTechInterp_en_korea,
                L222.GlobalTechCoef_en_korea,
                L222.GlobalTechShrwt_en_korea,
                L222.GlobalTechCapture_en_korea)
  } else {
    stop("Unknown command")
  }
}