Merge branch 'master' of https://github.com/pik-piam/remind2 into fix…

…-emissions-inconsistency

.buildlibrary

@@ -1,4 +1,4 @@
-ValidationKey: '223437494'
+ValidationKey: '223680318'
 AcceptedWarnings:
 - 'Warning: package ''.*'' was built under R version'
 - 'Warning: namespace ''.*'' is not available and has been replaced'

CITATION.cff

@@ -2,8 +2,8 @@ cff-version: 1.2.0
 message: If you use this software, please cite it using the metadata from this file.
 type: software
 title: 'remind2: The REMIND R package (2nd generation)'
-version: 1.131.1
-date-released: '2024-02-01'
+version: 1.132.1
+date-released: '2024-02-05'
 abstract: Contains the REMIND-specific routines for data and model output manipulation.
 authors:
 - family-names: Rodrigues

DESCRIPTION

@@ -1,8 +1,8 @@
 Type: Package
 Package: remind2
 Title: The REMIND R package (2nd generation)
-Version: 1.131.1
-Date: 2024-02-01
+Version: 1.132.1
+Date: 2024-02-05
 Authors@R: c(
     person("Renato", "Rodrigues", , "[email protected]", role = c("aut", "cre")),
     person("Lavinia", "Baumstark", role = "aut"),

R/compareScenConf.R

@@ -77,10 +77,15 @@ compareScenConf <- function(fileList = NULL, remindPath = "/p/projects/rd3mod/gi
   settings1 <- readCheckScenarioConfig(fileList[[1]], remindPath = remindPath, fillWithDefault = TRUE, testmode = TRUE)
   settings2 <- readCheckScenarioConfig(fileList[[2]], remindPath = remindPath, fillWithDefault = TRUE, testmode = TRUE)
 
-  # for mapping files
+  # for mapping files use "Variable" if exists, else combine first two columns
   if (is.null(row.names)) {
-    rownames(settings1) <- make.unique(paste0(settings1[, 1], ": ", settings1[, 2]))
-    rownames(settings2) <- make.unique(paste0(settings2[, 1], ": ", settings2[, 2]))
+    if ("Variable" %in% intersect(colnames(settings1), colnames(settings2))) {
+      rownames(settings1) <- make.unique(settings1[, "Variable"])
+      rownames(settings2) <- make.unique(settings2[, "Variable"])
+    } else {
+      rownames(settings1) <- make.unique(paste0(settings1[, 1], ": ", settings1[, 2]))
+      rownames(settings2) <- make.unique(paste0(settings2[, 1], ": ", settings2[, 2]))
+    }
   }
 
   # rename columns and rows in old file to new names after some checks

R/reportCosts.R

@@ -132,7 +132,7 @@ reportCosts <- function(gdx,output=NULL,regionSubsetList=NULL,t=c(seq(2005,2060,
     vm_demFeSector <- readGDX(gdx,name=c("vm_demFeSector"),field="l",restore_zeros=FALSE,format="first_found")
     demFeIndst.m <- readGDX(gdx,name=c("q37_demFeIndst"),types="equations",restore_zeros=FALSE,field="m",format="first_found")
     demFeBuild.m <- readGDX(gdx,name=c("q36_demFeBuild"),types="equations",restore_zeros=FALSE,field="m",format="first_found")
-    balFe.m <- readGDX(gdx,name=c("qm_balFe","q_balFe"),types="equations",restore_zeros=FALSE,field="m",format="first_found") 
+    balFe.m <- readGDX(gdx,name=c("q_balFe","qm_balFe"),types="equations",restore_zeros=FALSE,field="m",format="first_found") 
     #balFeCDR.m <- readGDX(gdx,name=c("q33_balFeCDR"),types="equations",restore_zeros=FALSE,field="m",format="first_found")
     #febalForUe.m <- dimReduce(febalForUe.m, dim_exclude = 2)
     demFeTrans.m <- readGDX(gdx,name=c("q35_demFeTrans"),types="equations",restore_zeros=FALSE,field="m",format="first_found")

R/reportEmi.R

@@ -201,6 +201,9 @@ reportEmi <- function(gdx, output = NULL, regionSubsetList = NULL, t = c(seq(200
   # variable to release captured CO2 when no CCU capacities are standing anymore vent captured CO2
   v_co2capturevalve <- readGDX(gdx, "v_co2capturevalve", field = "l", restore_zeros = F)[, t, ]
 
+  # maximum annual CO2 storage potential assumed
+  max_geolStorage <-  readGDX(gdx, "vm_co2CCS", field = "up", restore_zeros = F)[, t, ]  # CO2 captured per industry subsector
+
   # CO2 captured per industry subsector
   # NOTE: The parameter pm_IndstCO2Captured was calculated without taking into
   # account the different emission factors of energy carriers, so we recalculate
@@ -1626,6 +1629,19 @@ if (!is.null(vm_plasticsCarbon)) {
 
   ### carbon storage ----
 
+  # maximum annual carbon storage and share that is used
+  out <- mbind(out,
+               setNames(dimSums(max_geolStorage, dim = 3, na.rm = T) * GtC_2_MtCO2,
+                        "Carbon Management|Storage|Maximum annual CO2 storage potential (Mt CO2/yr)")
+  )
+  # share of annual storage potential used
+  out <- mbind(out,
+
+             setNames(dimSums(vm_co2CCS, dim = 3, na.rm = T) / dimSums(max_geolStorage, dim = 3, na.rm = T) * 100,
+                      "Carbon Management|Storage|Share of annual potential used (%)")%>%
+            ifelse(is.finite(.), ., 0)
+  )
+
   # calculate carbon storage variables
   out <- mbind(out,
 
@@ -2676,19 +2692,24 @@ if (!is.null(vm_plasticsCarbon)) {
   ## aggregate intensive variables ----
   .regionSubsetList <- c(list('GLO' = getItems(vm_co2CCS, dim = 'all_regi')),
                          regionSubsetList)
-
+  i <- seq_along(.regionSubsetList)
   for (i in seq_along(.regionSubsetList)) {
-    var <- 'Carbon Management|Share of Stored CO2 from Captured CO2 (%)'
-
+    var1 <- "Carbon Management|Share of Stored CO2 from Captured CO2 (%)"
+    var2 <- "Carbon Management|Storage|Share of annual potential used (%)"
     target_region  <- .regionSubsetList[i]
     source_regions <- .regionSubsetList[[i]]
-
-    out[names(target_region),,var] <- (
+    out[names(target_region),,var1] <- (
         dimSums(vm_co2CCS[source_regions,,],     dim = c(1, 3), na.rm = TRUE)
       / dimSums(vm_co2capture[source_regions,,], dim = c(1, 3))
       * 100
       ) %>%
       ifelse(is.finite(.), ., 0)   # set NaN (division by 0) to 0
+
+    out[names(target_region),,var2] <- (
+        dimSums(vm_co2CCS[source_regions,,], dim = c(1,3), na.rm = T) 
+      / dimSums(max_geolStorage[source_regions,,], dim = c(1,3), na.rm = T) 
+      * 100) %>%
+      ifelse(is.finite(.), ., 0)   # set NaN (division by 0) to 0
   }
 
 

R/reportLCOE.R

@@ -1189,8 +1189,8 @@ df.co2price.weighted <- df.pomeg.expand %>%
     sector.mapping <- c("build" = "buildings", "indst" = "industry", "trans" = "transport")
 
 
-    pm_tau_fe_tax <- readGDX(gdx, "pm_tau_fe_tax", restore_zeros = F)
-    pm_tau_fe_sub <- readGDX(gdx, "pm_tau_fe_sub", restore_zeros = F)
+    pm_tau_fe_tax <- readGDX(gdx, c("p21_tau_fe_tax","pm_tau_fe_tax"), restore_zeros = F)
+    pm_tau_fe_sub <- readGDX(gdx, c("p21_tau_fe_sub","pm_tau_fe_sub"), restore_zeros = F)
 
     df.taxrate <- as.quitte(pm_tau_fe_tax * 1.2 / s_twa2mwh * 1e12) %>%
                     rename(taxrate = value)

R/reportPrices.R

@@ -552,8 +552,8 @@ reportPrices <- function(gdx, output=NULL, regionSubsetList=NULL,
     "indst.fehos", "indst.fesos", "indst.feels", "indst.feh2s", "indst.fegas",
     "build.fepet", "indst.fepet"
   )
-  pm_tau_fe_tax <- readGDX(gdx, "pm_tau_fe_tax")[, YearsFrom2005, entyFe2Sector] # [tr USD2005/TWa]
-  pm_tau_fe_sub <- readGDX(gdx, "pm_tau_fe_sub")[, YearsFrom2005, entyFe2Sector] # [tr USD2005/TWa]
+  pm_tau_fe_tax <- readGDX(gdx, c("p21_tau_fe_tax","pm_tau_fe_tax"))[, YearsFrom2005, entyFe2Sector] # [tr USD2005/TWa]
+  pm_tau_fe_sub <- readGDX(gdx, c("p21_tau_fe_sub","pm_tau_fe_sub"))[, YearsFrom2005, entyFe2Sector] # [tr USD2005/TWa]
   price.tax <- (pm_tau_fe_tax + pm_tau_fe_sub) / s_twa2mwh / 3.6 * 1e12 # [USD2005/GJ]
 
   out <- mbind(

R/reportTax.R

@@ -45,8 +45,8 @@ reportTax <- function(gdx,output=NULL,regionSubsetList=NULL,t=c(seq(2005,2060,5)
   out <- NULL
 
   ### FE taxes/subsidies per sector
-  fe_tax  <- readGDX(gdx, name=c("pm_tau_fe_tax"), format="first_found", react = "silent")[,t,] * tdptwyr2dpgj
-  fe_sub  <- readGDX(gdx, name=c("pm_tau_fe_sub"), format="first_found", react = "silent")[,t,] * tdptwyr2dpgj
+  fe_tax  <- readGDX(gdx, name=c("p21_tau_fe_tax","pm_tau_fe_tax"), format="first_found", react = "silent")[,t,] * tdptwyr2dpgj
+  fe_sub  <- readGDX(gdx, name=c("p21_tau_fe_sub","pm_tau_fe_sub"), format="first_found", react = "silent")[,t,] * tdptwyr2dpgj
 
   vm_demFeSector <- readGDX(gdx,name=c("vm_demFeSector"),field="l",format="first_found",restore_zeros=FALSE)[,t,]*TWa_2_EJ
   vm_demFeSector[is.na(vm_demFeSector)] <- 0

README.md

@@ -1,6 +1,6 @@
 # The REMIND R package (2nd generation)
 
-R package **remind2**, version **1.131.1**
+R package **remind2**, version **1.132.1**
 
 [![CRAN status](https://www.r-pkg.org/badges/version/remind2)](https://cran.r-project.org/package=remind2)  [![R build status](https://github.com/pik-piam/remind2/workflows/check/badge.svg)](https://github.com/pik-piam/remind2/actions) [![codecov](https://codecov.io/gh/pik-piam/remind2/branch/master/graph/badge.svg)](https://app.codecov.io/gh/pik-piam/remind2) [![r-universe](https://pik-piam.r-universe.dev/badges/remind2)](https://pik-piam.r-universe.dev/builds)
 
@@ -49,7 +49,7 @@ In case of questions / problems please contact Renato Rodrigues <renato.rodrigue
 
 To cite package **remind2** in publications use:
 
-Rodrigues R, Baumstark L, Benke F, Dietrich J, Dirnaichner A, Duerrwaechter J, Führlich P, Giannousakis A, Hasse R, Hilaire J, Klein D, Koch J, Kowalczyk K, Levesque A, Malik A, Merfort A, Merfort L, Morena-Leiva S, Pehl M, Pietzcker R, Rauner S, Richters O, Rottoli M, Schötz C, Schreyer F, Siala K, Sörgel B, Spahr M, Strefler J, Verpoort P, Weigmann P (2024). _remind2: The REMIND R package (2nd generation)_. R package version 1.131.1, <https://github.com/pik-piam/remind2>.
+Rodrigues R, Baumstark L, Benke F, Dietrich J, Dirnaichner A, Duerrwaechter J, Führlich P, Giannousakis A, Hasse R, Hilaire J, Klein D, Koch J, Kowalczyk K, Levesque A, Malik A, Merfort A, Merfort L, Morena-Leiva S, Pehl M, Pietzcker R, Rauner S, Richters O, Rottoli M, Schötz C, Schreyer F, Siala K, Sörgel B, Spahr M, Strefler J, Verpoort P, Weigmann P (2024). _remind2: The REMIND R package (2nd generation)_. R package version 1.132.1, <URL: https://github.com/pik-piam/remind2>.
 
 A BibTeX entry for LaTeX users is
 
@@ -58,7 +58,7 @@ A BibTeX entry for LaTeX users is
   title = {remind2: The REMIND R package (2nd generation)},
   author = {Renato Rodrigues and Lavinia Baumstark and Falk Benke and Jan Philipp Dietrich and Alois Dirnaichner and Jakob Duerrwaechter and Pascal Führlich and Anastasis Giannousakis and Robin Hasse and Jérome Hilaire and David Klein and Johannes Koch and Katarzyna Kowalczyk and Antoine Levesque and Aman Malik and Anne Merfort and Leon Merfort and Simón Morena-Leiva and Michaja Pehl and Robert Pietzcker and Sebastian Rauner and Oliver Richters and Marianna Rottoli and Christof Schötz and Felix Schreyer and Kais Siala and Björn Sörgel and Mike Spahr and Jessica Strefler and Philipp Verpoort and Pascal Weigmann},
   year = {2024},
-  note = {R package version 1.131.1},
+  note = {R package version 1.132.1},
   url = {https://github.com/pik-piam/remind2},
 }
 ```

inst/markdown/compareScenarios2/cs2_09_carbon_management.Rmd

@@ -16,6 +16,10 @@ showAreaAndBarPlotsPlus(data, "Carbon Management|Carbon Capture|Industry Energy"
 showLinePlots(data, "Carbon Management|Storage")
 showAreaAndBarPlotsPlus(data, "Carbon Management|Storage", scales = "fixed")
 showAreaAndBarPlotsPlus(data, "Carbon Management|Storage|Industry Energy", scales = "fixed")
+showLinePlots(data, "Carbon Management|Storage|Share of annual potential used")
+showAreaAndBarPlots(data,
+                    vars = c("Carbon Management|Storage"),
+                    tot = "Carbon Management|Storage|Maximum annual CO2 storage potential")
 ```
 
 ## Usage