gcam-korea-release: Build the GCAM-Korea Input Datasets

Documented in module_emissions_L212.unmgd_nonco2

#' module_emissions_L212.unmgd_nonco2
#'
#' Outputs unmanaged land emissions and emissions coefficients for forest fires, deforestation, and grassland fires.
#'
#' @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{L212.AgSupplySector}, \code{L212.AgSupplySubsector}, \code{L212.ItemName},
#' \code{L212.GRASSEmissions}, \code{L212.FORESTEmissions_FF}, \code{L212.FORESTEmissions_D}, \code{L212.GRASSEmissionsFactors_BCOC},
#' \code{L212.FORESTEmissionsFactors_BCOC_FF}, \code{L212.FORESTEmissionsFactors_BCOC_D}, \code{L212.FORESTEmissionsFactors_future},
#' \code{L212.ItemName_prot}, \code{L212.GRASSEmissions_prot}, \code{L212.GRASSEmissions_noprot}, \code{L212.FORESTEmissions_FF_prot},
#' \code{L212.FORESTEmissions_FF_noprot}, \code{L212.FORESTEmissions_D_prot}, \code{L212.FORESTEmissions_D_noprot}, \code{L212.GRASSEmissionsFactors_BCOC_prot},
#' \code{L212.GRASSEmissionsFactors_BCOC_noprot}, \code{L212.FORESTEmissionsFactors_BCOC_FF_prot}, \code{L212.FORESTEmissionsFactors_BCOC_FF_noprot},
#' \code{L212.FORESTEmissionsFactors_BCOC_D_prot}, \code{L212.FORESTEmissionsFactors_BCOC_D_noprot}. The corresponding file in the
#' original data system was \code{L212.unmgd_nonco2.R} (emissions level2).
#' @details Outputs unmanaged land emissions and emissions coefficients for forest fires, deforestation, and grassland fires.
#' @importFrom assertthat assert_that
#' @importFrom dplyr filter mutate select
#' @importFrom tidyr gather spread
#' @author RLH August 2017
module_emissions_L212.unmgd_nonco2 <- function(command, ...) {
  if(command == driver.DECLARE_INPUTS) {
    return(c(FILE = "common/GCAM_region_names",
             FILE = "water/basin_to_country_mapping",
             FILE = "emissions/A_regions",
             "L124.nonco2_tg_R_grass_Y_GLU",
             "L124.nonco2_tg_R_forest_Y_GLU",
             "L124.deforest_coefs",
             "L125.bcoc_tgbkm2_R_grass_2000",
             "L125.bcoc_tgbkm2_R_forest_2000",
             "L125.deforest_coefs_bcoc"))
  } else if(command == driver.DECLARE_OUTPUTS) {
    return(c("L212.AgSupplySector",
             "L212.AgSupplySubsector",
             "L212.ItemName",
             "L212.GRASSEmissions",
             "L212.FORESTEmissions_FF",
             "L212.FORESTEmissions_D",
             "L212.GRASSEmissionsFactors_BCOC",
             "L212.FORESTEmissionsFactors_BCOC_FF",
             "L212.FORESTEmissionsFactors_BCOC_D",
             "L212.FORESTEmissionsFactors_future",
             "L212.ItemName_prot",
             "L212.GRASSEmissions_prot",
             "L212.GRASSEmissions_noprot",
             "L212.FORESTEmissions_FF_prot",
             "L212.FORESTEmissions_FF_noprot",
             "L212.FORESTEmissions_D_prot",
             "L212.FORESTEmissions_D_noprot",
             "L212.GRASSEmissionsFactors_BCOC_prot",
             "L212.GRASSEmissionsFactors_BCOC_noprot",
             "L212.FORESTEmissionsFactors_BCOC_FF_prot",
             "L212.FORESTEmissionsFactors_BCOC_FF_noprot",
             "L212.FORESTEmissionsFactors_BCOC_D_prot",
             "L212.FORESTEmissionsFactors_BCOC_D_noprot",
             "L212.FORESTEmissionsFactors_future_prot"))
  } else if(command == driver.MAKE) {

    # Silencing package checks
    GCAM_region_ID <- region <- AgSupplySubsector <- GLU <- itemName <- AgSupplySector <-
      UnmanagedLandTechnology <- year <- value <- Non.CO2 <- em_factor <- technology <-
      emiss.coef <- input.emissions <- NULL

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

    # Load required inputs
    GCAM_region_names <- get_data(all_data, "common/GCAM_region_names")
    basin_to_country_mapping <- get_data(all_data, "water/basin_to_country_mapping")
    A_regions <- get_data(all_data, "emissions/A_regions")
    L124.nonco2_tg_R_grass_Y_GLU <- get_data(all_data, "L124.nonco2_tg_R_grass_Y_GLU") %>%
      replace_GLU(basin_to_country_mapping)
    L124.nonco2_tg_R_forest_Y_GLU <- get_data(all_data, "L124.nonco2_tg_R_forest_Y_GLU") %>%
      replace_GLU(basin_to_country_mapping)
    L124.deforest_coefs <- get_data(all_data, "L124.deforest_coefs")
    L125.bcoc_tgbkm2_R_grass_2000 <- get_data(all_data, "L125.bcoc_tgbkm2_R_grass_2000")
    L125.bcoc_tgbkm2_R_forest_2000 <- get_data(all_data, "L125.bcoc_tgbkm2_R_forest_2000")
    L125.deforest_coefs_bcoc <- get_data(all_data, "L125.deforest_coefs_bcoc")

    # ===================================================
    # Unmanaged land sector info
    # L212.AgSupplySector: Sector info for unmanaged land technology
    # NOTE: only making unmanaged land sectors in regions that have aglu related emissions
    # Also note, not including the BC/OC tables here, which are from a different data source, and have no GLU detail
    # If regions/LTs have default BC/OC emissions factors but don't have any other gases (from EDGAR), they are dropped
    L212.regions <- unique(c(L124.nonco2_tg_R_grass_Y_GLU$GCAM_region_ID, L124.nonco2_tg_R_forest_Y_GLU$GCAM_region_ID))
    L212.regions <- GCAM_region_names %>%
      filter(GCAM_region_ID %in% L212.regions) %>%
      select(region)

    L212.AgSupplySector <- L212.regions %>%
      mutate(AgSupplySector = "UnmanagedLand",
             output.unit = "none",
             input.unit = "thous km2",
             price.unit = "none",
             calPrice = -1,
             market = region,
             logit.year.fillout = min(MODEL_BASE_YEARS),
             logit.exponent = -3,
             # Not in old-data-system, using because of removal of get_logit_fn_tables
             logit.type = NA)

    # L212.ItemName: Land item to relate emissions to
    # Note: only making unmanaged land technologies in the regions/glus with the corresponding land use type available
    L212.ItemName_R_LT_GLU <- bind_rows(unique(L124.nonco2_tg_R_grass_Y_GLU[c("GCAM_region_ID", "Land_Type", "GLU")]),
                                     unique(L124.nonco2_tg_R_forest_Y_GLU[c("GCAM_region_ID", "Land_Type", "GLU")]))

    L212.ItemName <- tibble(AgSupplySector = "UnmanagedLand",
                            AgSupplySubsector = c("ForestFire", "Deforest", "GrasslandFires"),
                            itemName = c("UnmanagedForest", "UnmanagedForest", "Grassland")) %>%
      # Add in region and GLU
      left_join(L212.ItemName_R_LT_GLU, by = c("itemName" = "Land_Type")) %>%
      left_join_error_no_match(GCAM_region_names, by = "GCAM_region_ID") %>%
      mutate(AgSupplySubsector = paste(AgSupplySubsector, GLU, sep = "_"),
             UnmanagedLandTechnology = AgSupplySubsector,
             itemName = paste(itemName, GLU, sep = "_")) %>%
      # Repeat for base years
      repeat_add_columns(tibble(year = MODEL_BASE_YEARS)) %>%
      select(region, AgSupplySector, AgSupplySubsector, UnmanagedLandTechnology, year, itemName)

    # Grassland emissions
    # L212.GrassEmissions: Grassland fire emissions in all regions
    L212.GRASSEmissions <- L124.nonco2_tg_R_grass_Y_GLU %>%
      filter(year %in% emissions.MODEL_BASE_YEARS) %>%
      mutate(value = round(value, emissions.DIGITS_EMISSIONS),
             AgSupplySector = "UnmanagedLand",
             AgSupplySubsector = paste("GrasslandFires", GLU, sep = "_"),
             UnmanagedLandTechnology = AgSupplySubsector) %>%
      left_join_error_no_match(GCAM_region_names, by = "GCAM_region_ID") %>%
      rename_SO2(A_regions) %>%
      select(region, AgSupplySector, AgSupplySubsector, UnmanagedLandTechnology, year, Non.CO2, input.emissions = value)

    # Grassland emissions factors for BC/OC
    # L212.GrassEmissions: Grassland fire emissions factors for BC/OC in all regions
    L212.GRASSEmissionsFactors_BCOC <- L212.GRASSEmissions %>%
      select(region, AgSupplySector, AgSupplySubsector, UnmanagedLandTechnology, year) %>%
      distinct() %>%
      repeat_add_columns(tibble(Non.CO2 = unique(L125.bcoc_tgbkm2_R_grass_2000$Non.CO2))) %>%
      # Add in GCAM_region_ID for joining
      left_join_error_no_match(GCAM_region_names, by = "region") %>%
      # Join in emission factor values
      left_join_error_no_match(L125.bcoc_tgbkm2_R_grass_2000, by = c("Non.CO2", "GCAM_region_ID")) %>%
      mutate(em_factor = round(em_factor, emissions.DIGITS_EMISSIONS)) %>%
      select(region, AgSupplySector, AgSupplySubsector, UnmanagedLandTechnology, year, Non.CO2, emiss.coef = em_factor)

    # L212.Forest: Forest fire and deforestation emissions in all regions
    # Will split up by technology in final product creation
    L212.FOREST <- L124.nonco2_tg_R_forest_Y_GLU %>%
      filter(year %in% emissions.MODEL_BASE_YEARS) %>%
      mutate(value = round(value, emissions.DIGITS_EMISSIONS),
             AgSupplySector = "UnmanagedLand",
             AgSupplySubsector = paste(technology, GLU, sep = "_"),
             UnmanagedLandTechnology = AgSupplySubsector) %>%
      left_join_error_no_match(GCAM_region_names, by = "GCAM_region_ID") %>%
      rename_SO2(A_regions) %>%
      select(region, AgSupplySector, AgSupplySubsector, UnmanagedLandTechnology, technology, year, Non.CO2, input.emissions = value)

    # L212.FORESTEmissionsFactors_BCOC: Forest emissions factors for BC/OC in all regions
    # Will split up by technology in final product creation
    L212.FORESTEmissionsFactors_BCOC <- L212.FOREST %>%
      select(region, AgSupplySector, AgSupplySubsector, UnmanagedLandTechnology, technology, year) %>%
      distinct() %>%
      repeat_add_columns(tibble(Non.CO2 = unique(L125.bcoc_tgbkm2_R_forest_2000$Non.CO2))) %>%
      # Add in GCAM_region_ID for joining
      left_join_error_no_match(GCAM_region_names, by = "region") %>%
      # Join in emission factor values
      left_join_error_no_match(L125.bcoc_tgbkm2_R_forest_2000, by = c("Non.CO2", "technology", "GCAM_region_ID")) %>%
      mutate(em_factor = round(em_factor, emissions.DIGITS_EMISSIONS)) %>%
      select(region, AgSupplySector, AgSupplySubsector, UnmanagedLandTechnology, year, Non.CO2, emiss.coef = em_factor, technology)

    # Reading in default emissions factors for deforestation in future periods.
    # Because the model's driver for deforestation is not necessarily linked to the estimated quantity of emissions,
    # the deforestation-related emissions in the base years may be zero (regions with net afforestation), or very high
    # (regions with slightly net deforestation). This method just uses a global average for future non-co2 emissions
    # from deforestation
    L212.default_coefs <- bind_rows(select(L124.deforest_coefs, Non.CO2, emiss.coef),
                                    select(L125.deforest_coefs_bcoc, Non.CO2, emiss.coef))

    L212.FORESTEmissionsFactors_future <- L212.FOREST %>%
      filter(technology == "Deforest") %>%
      select(region, AgSupplySector, AgSupplySubsector, UnmanagedLandTechnology) %>%
      distinct() %>%
      # Take minimum model year greater than emissions model base years
      mutate(year = as.integer(min(MODEL_YEARS[MODEL_YEARS > max(emissions.MODEL_BASE_YEARS)]))) %>%
      repeat_add_columns(L212.default_coefs) %>%
      mutate(emiss.coef = round(emiss.coef, emissions.DIGITS_EMISSIONS)) %>%
      rename_SO2(A_regions) %>%
      select(region, AgSupplySector, AgSupplySubsector, UnmanagedLandTechnology, year, Non.CO2, emiss.coef)

    # We do not actually care about the logit here but we need a value to avoid errors
    L212.AgSupplySubsector <- L212.ItemName %>%
      select(region, AgSupplySector, AgSupplySubsector) %>%
      mutate(logit.year.fillout = min(MODEL_BASE_YEARS),
             logit.exponent = -3,
             logit.type = NA)

    # Writing files with modified emissions for protected lands.
    # Note: the protected lands input file is intended to be read after the normal one,
    # so rather than copy all of the duplicate information (e.g., supplysector and subsector, itemname),
    # these files only read in the new technologies and modified emissions levels

    # Land item: protected land use types now prefixed with protected; non-protected keeps the same name so no need to write out
    # Land item: copy for prot and noprot, prefix the LTs with "Protected" in the prot file
    L212.ItemName_prot <- L212.ItemName %>%
      mutate(UnmanagedLandTechnology = paste0("Protected", UnmanagedLandTechnology),
             itemName = paste0("Protected", itemName))

    # Emissions: names in protected files are prefixed, and emissions in both are multiplied by protected land fraction
    # Grassland emissions
    L212.GRASSEmissions_prot <- L212.GRASSEmissions %>%
      mutate(UnmanagedLandTechnology = paste0("Protected", UnmanagedLandTechnology),
             input.emissions = input.emissions * aglu.PROTECT_LAND_FRACT)

    L212.GRASSEmissions_noprot <- L212.GRASSEmissions %>%
      mutate(input.emissions = input.emissions * (1 - aglu.PROTECT_LAND_FRACT))

    # Forest emissions - fires and deforest
    L212.FORESTEmissions_prot <- L212.FOREST %>%
      mutate(UnmanagedLandTechnology = paste0("Protected", UnmanagedLandTechnology),
             input.emissions = input.emissions * aglu.PROTECT_LAND_FRACT)

    L212.FORESTEmissions_noprot <- L212.FOREST %>%
      mutate(input.emissions = input.emissions * (1 - aglu.PROTECT_LAND_FRACT))

    # Emissions factors: names in protected files are prefixed, and factors are left unchanged
    L212.GRASSEmissionsFactors_BCOC_prot <- L212.GRASSEmissionsFactors_BCOC %>%
      mutate(UnmanagedLandTechnology = paste0("Protected", UnmanagedLandTechnology))

    L212.FORESTEmissionsFactors_BCOC_prot <- L212.FORESTEmissionsFactors_BCOC %>%
      mutate(UnmanagedLandTechnology = paste0("Protected", UnmanagedLandTechnology))

    # Create future emissions factors for protected lands
    L212.FORESTEmissionsFactors_future_prot <- L212.FORESTEmissionsFactors_future %>%
      mutate(UnmanagedLandTechnology = paste0("Protected", UnmanagedLandTechnology))

    # ===================================================
    # Produce outputs
    L212.AgSupplySector %>%
      add_title("Logit Exponents for Unmanaged Land Sector") %>%
      add_units("Unitless") %>%
      add_comments("Hard-coded constant values") %>%
      add_legacy_name("L212.AgSupplySector") %>%
      add_precursors("common/GCAM_region_names", "L124.nonco2_tg_R_grass_Y_GLU",
                     "L124.nonco2_tg_R_forest_Y_GLU") ->
      L212.AgSupplySector

    L212.ItemName %>%
      add_title("Mapping File for Unmanaged Land") %>%
      add_units("NA") %>%
      add_comments("Mapping categories from L124.nonco2_tg_R_grass_Y_GLU, L124.nonco2_tg_R_forest_Y_GLU, and constants") %>%
      add_legacy_name("L212.ItemName") %>%
      add_precursors("common/GCAM_region_names", "water/basin_to_country_mapping",
                     "L124.nonco2_tg_R_grass_Y_GLU", "L124.nonco2_tg_R_forest_Y_GLU") ->
      L212.ItemName

    L212.AgSupplySubsector %>%
      add_title("Logit Exponents for Unmanaged Land Subsectors") %>%
      add_units("Unitless") %>%
      add_comments("Constant logit exponent added to L212.ItemName") %>%
      add_legacy_name("L212.AgSupplySubsector") %>%
      same_precursors_as(L212.ItemName) ->
      L212.AgSupplySubsector

    L212.GRASSEmissions %>%
      add_title("Historical Emissions-Grassland Fires") %>%
      add_units("Tg/yr") %>%
      add_comments("Supply sectors and region names added to L124.nonco2_tg_R_grass_Y_GLU") %>%
       add_legacy_name("L212.GRASSEmissions") %>%
      add_precursors("common/GCAM_region_names", "emissions/A_regions",
                     "L124.nonco2_tg_R_grass_Y_GLU", "water/basin_to_country_mapping") ->
      L212.GRASSEmissions

    L212.FOREST %>%
      filter(technology == "ForestFire") %>%
      select(-technology) %>%
      add_title("Historical Emissions-Forest Fires") %>%
      add_units("Tg/yr") %>%
      add_comments("Supply sectors and region names added to L124.nonco2_tg_R_forest_Y_GLU") %>%
      add_legacy_name("L212.FORESTEmissions_FF") %>%
      add_precursors("common/GCAM_region_names", "emissions/A_regions",
                     "L124.nonco2_tg_R_forest_Y_GLU", "water/basin_to_country_mapping") ->
      L212.FORESTEmissions_FF

    L212.FOREST %>%
      filter(technology == "Deforest") %>%
      select(-technology) %>%
      add_title("Historical Emissions-Deforestation") %>%
      add_units("Tg/yr") %>%
      add_comments("Supply sectors and region names added to L124.nonco2_tg_R_forest_Y_GLU") %>%
      add_legacy_name("L212.FORESTEmissions_D") %>%
      same_precursors_as(L212.FORESTEmissions_FF) ->
      L212.FORESTEmissions_D

    L212.GRASSEmissionsFactors_BCOC %>%
      add_title("Historical BC/OC Emissions Coefficients-Grassland Fires") %>%
      add_units("Tg/bm2") %>%
      add_comments("comments describing how data generated") %>%
      add_comments("can be multiple lines") %>%
      add_legacy_name("L212.GRASSEmissionsFactors_BCOC") %>%
      add_precursors("common/GCAM_region_names", "L125.bcoc_tgbkm2_R_grass_2000",
                     "L124.nonco2_tg_R_grass_Y_GLU", "water/basin_to_country_mapping") ->
      L212.GRASSEmissionsFactors_BCOC

    L212.FORESTEmissionsFactors_BCOC %>%
      filter(technology == "ForestFire") %>%
      select(-technology) %>%
      add_title("Historical BC/OC Emissions Coefficients-Forest Fires") %>%
      add_units("Tg/bm2") %>%
      add_comments("L125.bcoc_tgbkm2_R_grass_2000 emissions factors added to L212.GRASSEmissions") %>%
      add_comments("can be multiple lines") %>%
      add_legacy_name("L212.FORESTEmissionsFactors_BCOC_FF") %>%
      add_precursors("common/GCAM_region_names", "L125.bcoc_tgbkm2_R_forest_2000",
                     "L124.nonco2_tg_R_forest_Y_GLU", "water/basin_to_country_mapping") ->
      L212.FORESTEmissionsFactors_BCOC_FF

    L212.FORESTEmissionsFactors_BCOC %>%
      filter(technology == "Deforest") %>%
      select(-technology) %>%
      add_title("Historical BC/OC Emissions Coefficients-Deforestation") %>%
      add_units("Tg/bm2") %>%
      add_comments("L125.bcoc_tgbkm2_R_forest_2000 emissions factors") %>%
      add_legacy_name("L212.FORESTEmissionsFactors_BCOC_D") %>%
      same_precursors_as(L212.FORESTEmissionsFactors_BCOC_FF) ->
      L212.FORESTEmissionsFactors_BCOC_D

    L212.FORESTEmissionsFactors_future %>%
      add_title("Future BC/OC Emissions Coefficients-Deforestation") %>%
      add_units("kg/m2") %>%
      add_comments("L212.default_coefs values added to first future year") %>%
      add_legacy_name("L212.FORESTEmissionsFactors_future") %>%
      add_precursors("common/GCAM_region_names", "emissions/A_regions",
                     "L124.nonco2_tg_R_forest_Y_GLU", "water/basin_to_country_mapping",
                     "L124.deforest_coefs", "L125.deforest_coefs_bcoc") ->
      L212.FORESTEmissionsFactors_future

    L212.FORESTEmissionsFactors_future_prot %>%
      add_title("Future BC/OC Emissions Coefficients-Deforestation for Protected Lands") %>%
      add_units("kg/m2") %>%
      add_comments("L212.default_coefs values added to first future year") %>%
      add_legacy_name("L212.FORESTEmissionsFactors_future") %>%
      add_precursors("common/GCAM_region_names", "emissions/A_regions",
                     "L124.nonco2_tg_R_forest_Y_GLU", "water/basin_to_country_mapping",
                     "L124.deforest_coefs", "L125.deforest_coefs_bcoc") ->
      L212.FORESTEmissionsFactors_future_prot

    L212.ItemName_prot %>%
      add_title("Mapping File for Protected Unmanaged Land") %>%
      add_units("NA") %>%
      add_comments("Same info as L212.ItemName") %>%
      add_legacy_name("L212.ItemName_prot") %>%
      same_precursors_as(L212.ItemName) ->
      L212.ItemName_prot

    L212.GRASSEmissions_prot %>%
      add_title("Historical Protected Emissions-Grassland Fires") %>%
      add_units("Tg/yr") %>%
      add_comments("L212.GRASSEmissions values multiplied by protected land fraction") %>%
      add_legacy_name("L212.GRASSEmissions_prot") %>%
      same_precursors_as(L212.GRASSEmissions) ->
      L212.GRASSEmissions_prot

    L212.GRASSEmissions_noprot %>%
      add_title("Historical Unprotected Emissions-Grassland Fires") %>%
      add_units("Tg/yr") %>%
      add_comments("L212.GRASSEmissions values multiplied by 1 minus protected land fraction") %>%
      add_legacy_name("L212.GRASSEmissions_noprot") %>%
      same_precursors_as(L212.GRASSEmissions) ->
      L212.GRASSEmissions_noprot

    L212.FORESTEmissions_prot %>%
      filter(technology == "ForestFire") %>%
      select(-technology) %>%
      add_title("Historical Protected Emissions-Forest Fires") %>%
      add_units("Tg/yr") %>%
      add_comments("L212.FORESTEmissions_FF values multiplied by protected land fraction") %>%
      add_legacy_name("L212.FORESTEmissions_FF_prot") %>%
      same_precursors_as(L212.FORESTEmissions_FF) ->
      L212.FORESTEmissions_FF_prot

    L212.FORESTEmissions_noprot %>%
      filter(technology == "ForestFire") %>%
      select(-technology) %>%
      add_title("Historical Unprotected Emissions-Forest Fires") %>%
      add_units("Tg/yr") %>%
      add_comments("L212.FORESTEmissions_FF values multiplied by 1 minus protected land fraction") %>%
      add_legacy_name("L212.FORESTEmissions_FF_noprot") %>%
      same_precursors_as(L212.FORESTEmissions_FF) ->
      L212.FORESTEmissions_FF_noprot

    L212.FORESTEmissions_prot %>%
      filter(technology == "Deforest") %>%
      select(-technology) %>%
      add_title("Historical Protected Emissions-Deforestation") %>%
      add_units("Tg/yr") %>%
      add_comments("L212.FORESTEmissions_D values multiplied by protected land fraction") %>%
      add_legacy_name("L212.FORESTEmissions_D_prot") %>%
      same_precursors_as(L212.FORESTEmissions_D) ->
      L212.FORESTEmissions_D_prot

    L212.FORESTEmissions_noprot %>%
      filter(technology == "Deforest") %>%
      select(-technology) %>%
      add_title("Historical Unprotected Emissions-Deforestation") %>%
      add_units("Tg/yr") %>%
      add_comments("L212.FORESTEmissions_D values multiplied by 1 minus protected land fraction") %>%
      add_legacy_name("L212.FORESTEmissions_D_noprot") %>%
      same_precursors_as(L212.FORESTEmissions_D) ->
      L212.FORESTEmissions_D_noprot

    L212.GRASSEmissionsFactors_BCOC_prot %>%
      add_title("BC/OC Protected Emissions Factors-Grassland Fires") %>%
      add_units("Tg/bm2") %>%
      add_comments("Same values as L212.GRASSEmissionsFactors_BCOC") %>%
      add_legacy_name("L212.GRASSEmissionsFactors_BCOC_prot") %>%
      same_precursors_as(L212.GRASSEmissionsFactors_BCOC) ->
      L212.GRASSEmissionsFactors_BCOC_prot

    L212.GRASSEmissionsFactors_BCOC %>%
      # Useless mutate to add different attributes
      mutate(region = region) %>%
      add_title("BC/OC Unprotected Emissions Factors-Grassland Fires") %>%
      add_units("Tg/bm2") %>%
      add_comments("Same values as L212.GRASSEmissionsFactors_BCOC") %>%
      add_legacy_name("L212.GRASSEmissionsFactors_BCOC_noprot") %>%
      same_precursors_as(L212.GRASSEmissionsFactors_BCOC) ->
      L212.GRASSEmissionsFactors_BCOC_noprot

    L212.FORESTEmissionsFactors_BCOC_prot %>%
      filter(technology == "ForestFire") %>%
      select(-technology) %>%
      add_title("BC/OC Protected Emissions Factors-Forest Fires") %>%
      add_units("Tg/bm2") %>%
      add_comments("Same values as L212.FORESTEmissionsFactors_BCOC_FF") %>%
      add_legacy_name("L212.FORESTEmissionsFactors_BCOC_FF_prot") %>%
      same_precursors_as(L212.FORESTEmissionsFactors_BCOC_FF) ->
      L212.FORESTEmissionsFactors_BCOC_FF_prot

    L212.FORESTEmissionsFactors_BCOC_FF %>%
      # Useless mutate to add different attributes
      mutate(region = region) %>%
      add_title("BC/OC Unprotected Emissions Factors-Forest Fires") %>%
      add_units("Tg/bm2") %>%
      add_comments("Same values as L212.FORESTEmissionsFactors_BCOC_FF") %>%
      add_legacy_name("L212.FORESTEmissionsFactors_BCOC_FF_noprot") %>%
      same_precursors_as(L212.FORESTEmissionsFactors_BCOC_FF) ->
      L212.FORESTEmissionsFactors_BCOC_FF_noprot

    L212.FORESTEmissionsFactors_BCOC_prot %>%
      filter(technology == "Deforest") %>%
      select(-technology) %>%
      add_title("BC/OC Protected Emissions Factors-Deforestation") %>%
      add_units("Tg/bm2") %>%
      add_comments("Same values as L212.FORESTEmissionsFactors_BCOC_D") %>%
      add_legacy_name("L212.FORESTEmissionsFactors_BCOC_D_prot") %>%
      same_precursors_as(L212.FORESTEmissionsFactors_BCOC_D) ->
      L212.FORESTEmissionsFactors_BCOC_D_prot

    L212.FORESTEmissionsFactors_BCOC_D %>%
      # Useless mutate to add different attributes
      mutate(region = region) %>%
      add_title("BC/OC Unprotected Emissions Factors-Deforestation") %>%
      add_units("Tg/bm2") %>%
      add_comments("Same values as L212.FORESTEmissionsFactors_BCOC_D") %>%
      add_legacy_name("L212.FORESTEmissionsFactors_BCOC_D_noprot") %>%
      same_precursors_as(L212.FORESTEmissionsFactors_BCOC_D) ->
      L212.FORESTEmissionsFactors_BCOC_D_noprot

    return_data(L212.AgSupplySector, L212.AgSupplySubsector, L212.ItemName, L212.GRASSEmissions,
                L212.FORESTEmissions_FF, L212.FORESTEmissions_D, L212.GRASSEmissionsFactors_BCOC,
                L212.FORESTEmissionsFactors_BCOC_FF, L212.FORESTEmissionsFactors_BCOC_D,
                L212.FORESTEmissionsFactors_future, L212.ItemName_prot, L212.GRASSEmissions_prot,
                L212.GRASSEmissions_noprot, L212.FORESTEmissions_FF_prot, L212.FORESTEmissions_FF_noprot,
                L212.FORESTEmissions_D_prot, L212.FORESTEmissions_D_noprot, L212.GRASSEmissionsFactors_BCOC_prot,
                L212.GRASSEmissionsFactors_BCOC_noprot, L212.FORESTEmissionsFactors_BCOC_FF_prot,
                L212.FORESTEmissionsFactors_BCOC_FF_noprot, L212.FORESTEmissionsFactors_BCOC_D_prot,
                L212.FORESTEmissionsFactors_BCOC_D_noprot, L212.FORESTEmissionsFactors_future_prot)
  } else {
    stop("Unknown command")
  }
}