Merge pull request #1833 from laurinks/develop

Adding ScenarioMIP to list of projects for IIASA export
remindmodel · Sep 17, 2024 · 5ea14e2 · 5ea14e2
2 parents a24ea7b + 7059c3d
commit 5ea14e2
Showing 1 changed file with 20 additions and 4 deletions.
diff --git a/scripts/output/export/xlsx_IIASA.R b/scripts/output/export/xlsx_IIASA.R
@@ -21,11 +21,13 @@ model <- paste("REMIND", paste0(strsplit(gms::readDefaultConfig(".")$model_versi
 
 
 removeFromScen <- ""                           # you can use regex such as: "_diff|_expoLinear"
+renameScen <- NULL                             # c(oldname1 = "newname1", …), without the `C_` and `-rem-[0-9]` stuff
 addToScen <- NULL                              # is added at the beginning
 
 # filenames relative to REMIND main directory (or use absolute path) 
 mapping <- NULL                                # file obtained from piamInterfaces, or AR6/SHAPE/NAVIGATE or NULL to get asked
 iiasatemplate <- NULL                          # provided for each project, can be yaml or xlsx with a column 'Variable'
+checkSummation <- TRUE                         # if TRUE, tries to use the one from mapping. Or specify here
 
 # note: you can also pass all these options to output.R, so 'Rscript output.R logFile=mylogfile.txt' works.
 lucode2::readArgs("project")
@@ -41,6 +43,11 @@ projects <- list(
                     mapping = c("AR6", "AR6_NGFS"),
                     iiasatemplate = "https://files.ece.iiasa.ac.at/ngfs-phase-5/ngfs-phase-5-template.xlsx",
                     removeFromScen = "C_|_bIT|_bit|_bIt|_KLW"),
+  ScenarioMIP = list(model = "REMIND-MAgPIE 3.4-4.8",
+                     mapping = "ScenarioMIP",
+                     iiasatemplate = "https://files.ece.iiasa.ac.at/ssp-submission/ssp-submission-template.xlsx",
+                     renameScen = c("SMIPv03-M-SSP2-NPi-def" = "SSP2 - Medium Emissions", "SMIPv03-LOS-SSP2-EcBudg400-def" = "SSP2 - Low Overshoot", "SMIPv03-ML-SSP2-PkPrice200-fromL" = "SSP2 - Medium-Low Emissions","SMIPv03-L-SSP2-PkPrice265-inc6-def" = "SSP2 - Low Emissions", "SMIPv03-VL-SSP2_SDP_MC-PkPrice300-def" = "SSP2 - Very Low Emissions"),
+                     checkSummation = "NAVIGATE"),
   SHAPE      = list(mapping = c("NAVIGATE", "SHAPE")),
   TESTTHAT   = list(mapping = "AR6")
 )
@@ -58,14 +65,15 @@ if (! exists("project")) {
 }
 projectdata <- projects[[project]]
 message("# Overwrite settings with project settings for '", project, "'.")
-varnames <- c("mapping", "iiasatemplate", "addToScen", "removeFromScen", "model", "outputFilename", "logFile")
+varnames <- c("mapping", "iiasatemplate", "addToScen", "removeFromScen", "renameScen",
+              "model", "outputFilename", "logFile", "checkSummation")
 for (p in intersect(varnames, names(projectdata))) {
   assign(p, projectdata[[p]])
 }
 
 # overwrite settings with those specified as command-line arguments
-lucode2::readArgs("outputdirs", "filename_prefix", "outputFilename", "model",
-                  "mapping", "logFile", "removeFromScen", "addToScen", "iiasatemplate")
+lucode2::readArgs("outputdirs", "filename_prefix", "outputFilename", "model", "mapping",
+                  "summationFile", "logFile", "removeFromScen", "addToScen", "iiasatemplate")
 
 if (is.null(mapping)) {
   mapping <- gms::chooseFromList(names(piamInterfaces::mappingNames()), type = "mapping")
@@ -130,6 +138,14 @@ withCallingHandlers({ # piping messages to logFile
     mifdata <- rbind(mifdata, thismifdata)
   }
 
+  mifdata$scenario <- gsub("^C_", "", mifdata$scenario)
+  message("Old names: ", paste(sort(unique(mifdata$scenario)), collapse = ", "))
+  for (i in names(renameScen)) {
+    message("Rename scenario: ", i, " -> ", renameScen[[i]])
+    mifdata$scenario[i == mifdata$scenario] <- renameScen[[i]]
+  }
+  message("New names: ", paste(sort(unique(mifdata$scenario)), collapse = ", "))
+
   message("# ", length(temporarydelete), " variables are in the list to be temporarily deleted, ",
           length(unique(mifdata$variable[mifdata$variable %in% temporarydelete])), " were deleted.")
   write(paste0("  - ", paste(unique(mifdata$variable[mifdata$variable %in% temporarydelete]), collapse = "\n  - "), "\n\n"),
@@ -141,7 +157,7 @@ withCallingHandlers({ # piping messages to logFile
 
   generateIIASASubmission(mifdata, mapping = mapping, model = model,
                           removeFromScen = removeFromScen, addToScen = addToScen,
-                          outputDirectory = outputFolder,
+                          outputDirectory = outputFolder, checkSummation = checkSummation,
                           logFile = logFile, outputFilename = basename(OUTPUT_xlsx),
                           iiasatemplate = iiasatemplate, generatePlots = TRUE)