From 6f5dfe7966b83e60ffff7dabea22152ccd98c178 Mon Sep 17 00:00:00 2001 From: Falk Benke Date: Wed, 17 Jul 2024 17:28:39 +0200 Subject: [PATCH] convert prices in calcEU_ReferenceScenario --- R/calcEU_ReferenceScenario.R | 47 +++++++++++++++++++++++++++++------- 1 file changed, 38 insertions(+), 9 deletions(-) diff --git a/R/calcEU_ReferenceScenario.R b/R/calcEU_ReferenceScenario.R index d06ac50d..01f71253 100644 --- a/R/calcEU_ReferenceScenario.R +++ b/R/calcEU_ReferenceScenario.R @@ -1,18 +1,47 @@ calcEU_ReferenceScenario <- function() { - EU_ReferenceScenario_2016 <- readSource("EU_ReferenceScenario", subtype = "2016") - EU_ReferenceScenario_2016 <- add_dimension(EU_ReferenceScenario_2016, dim = 3.1, add = "model", - nm = "EU_ReferenceScenario_2016") + map <- toolGetMapping(type = "regional", name = "regionmapping_21_EU11.csv", where = "mappingfolder") %>% + dplyr::select("iso3c" = "CountryCode", "region" = "RegionCode") - EU_ReferenceScenario_2020 <- readSource("EU_ReferenceScenario", subtype = "2020") - EU_ReferenceScenario_2020 <- add_columns(EU_ReferenceScenario_2020, "y2000", dim = 2) - EU_ReferenceScenario_2020 <- add_dimension(EU_ReferenceScenario_2020, dim = 3.1, add = "model", - nm = "EU_ReferenceScenario_2020") + euRef2016 <- readSource("EU_ReferenceScenario", subtype = "2016") - x <- mbind(EU_ReferenceScenario_2016, EU_ReferenceScenario_2020) + # convert EUR2013 -> US$2017 + tmp <- GDPuc::convertGDP( + gdp = euRef2016[, , "Price|Secondary Energy|Electricity (EUR2013/GJ)"], + unit_in = "constant 2013 €", + unit_out = "constant 2017 Int$PPP", + with_regions = map, + replace_NAs = "regional_average" + ) + + getNames(tmp) <- "Price|Secondary Energy|Electricity (US$2017/GJ)" + euRef2016 <- euRef2016[, , "Price|Secondary Energy|Electricity (EUR2013/GJ)", invert = TRUE] + euRef2016 <- mbind(euRef2016, tmp) + + euRef2016 <- add_dimension(euRef2016, dim = 3.1, add = "model", nm = "EU_ReferenceScenario_2016") + + euRef2020 <- readSource("EU_ReferenceScenario", subtype = "2020") + + # convert EUR2015 -> US$2017 + tmp <- GDPuc::convertGDP( + gdp = euRef2020[, , "Price|Secondary Energy|Electricity (EUR2015/GJ)"], + unit_in = "constant 2015 €", + unit_out = "constant 2017 Int$PPP", + with_regions = map, + replace_NAs = "regional_average" + ) + + getNames(tmp) <- "Price|Secondary Energy|Electricity (US$2017/GJ)" + euRef2020 <- euRef2020[, , "Price|Secondary Energy|Electricity (EUR2015/GJ)", invert = TRUE] + euRef2020 <- mbind(euRef2020, tmp) + + euRef2020 <- add_columns(euRef2020, "y2000", dim = 2) + euRef2020 <- add_dimension(euRef2020, dim = 3.1, add = "model", nm = "EU_ReferenceScenario_2020") + + x <- mbind(euRef2016, euRef2020) weights <- x weights[, , ] <- NA - weights[, , "Price|Secondary Energy|Electricity (US$2005/GJ)"] <- 1 + weights[, , "Price|Secondary Energy|Electricity (US$2017/GJ)"] <- 1 return(list( x = x, weight = weights, mixed_aggregation = TRUE,