convert prices in calcEU_ReferenceScenario

pik-piam · Jul 17, 2024 · 6f5dfe7 · 6f5dfe7
1 parent 2d3370c
commit 6f5dfe7
Showing 1 changed file with 38 additions and 9 deletions.
diff --git 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,