-
Notifications
You must be signed in to change notification settings - Fork 42
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
- Loading branch information
1 parent
5ebd6c3
commit e815803
Showing
7 changed files
with
145 additions
and
10 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -1,8 +1,8 @@ | ||
| Type: Package | ||
| Package: remind2 | ||
| Title: The REMIND R package (2nd generation) | ||
| Version: 1.135.9 | ||
| Date: 2024-03-04 | ||
| Version: 1.135.9.9001 | ||
| Date: 2024-03-08 | ||
| Authors@R: c( | ||
| person("Renato", "Rodrigues", , "[email protected]", role = c("aut", "cre")), | ||
| person("Lavinia", "Baumstark", role = "aut"), | ||
|
|
||
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,135 @@ | ||
| # | (C) 2006-2023 Potsdam Institute for Climate Impact Research (PIK) | ||
| # | authors, and contributors see CITATION.cff file. This file is part | ||
| # | of REMIND and licensed under AGPL-3.0-or-later. Under Section 7 of | ||
| # | AGPL-3.0, you are granted additional permissions described in the | ||
| # | REMIND License Exception, version 1.0 (see LICENSE file). | ||
| # | Contact: [email protected] | ||
|
|
||
| library(magclass) | ||
| library(remind2) | ||
| library(lucode2) | ||
| library(gms) | ||
| library(methods) | ||
| library(edgeTransport) | ||
| library(quitte) | ||
| ############################# BASIC CONFIGURATION ############################# | ||
| gdx_name <- "fulldata.gdx" # name of the gdx | ||
| gdx_ref_name <- "input_ref.gdx" # name of the ref for < cm_startyear | ||
| gdx_refpolicycost_name <- "input_refpolicycost.gdx" # name of the reference gdx (for policy cost calculation) | ||
|
|
||
|
|
||
| if(!exists("source_include")) { | ||
| # Define arguments that can be read from command line | ||
| outputdir <- "." | ||
| readArgs("outputdir", "gdx_name", "gdx_ref_name", "gdx_refpolicycost_name") | ||
| } | ||
|
|
||
| gdx <- file.path(outputdir, gdx_name) | ||
| gdx_ref <- file.path(outputdir, gdx_ref_name) | ||
| gdx_refpolicycost <- file.path(outputdir, gdx_refpolicycost_name) | ||
| if (! file.exists(gdx_ref)) gdx_ref <- NULL | ||
| if (! file.exists(gdx_refpolicycost)) gdx_refpolicycost <- NULL | ||
| scenario <- getScenNames(outputdir) | ||
| ############################################################################### | ||
| # paths of the reporting files | ||
| remind_reporting_file <- file.path(outputdir,paste0("REMIND_generic_", scenario,".mif")) | ||
| magicc_reporting_file <- file.path(outputdir,paste0("REMIND_climate_", scenario, ".mif")) | ||
| LCOE_reporting_file <- file.path(outputdir,paste0("REMIND_LCOE_", scenario, ".csv")) | ||
|
|
||
| remind_policy_reporting_file <- file.path(outputdir,paste0("REMIND_generic_",scenario,"_adjustedPolicyCosts.mif")) | ||
| remind_policy_reporting_file <- remind_policy_reporting_file[file.exists(remind_policy_reporting_file)] | ||
| if (length(remind_policy_reporting_file) > 0) { | ||
| unlink(remind_policy_reporting_file) | ||
| message("\n", paste(basename(remind_policy_reporting_file), collapse = ", "), " deleted.") | ||
| message(paste(basename(remind_reporting_file), collapse = ", "), " will contain policy costs based on input_refpolicycost.gdx.") | ||
| } | ||
|
|
||
| configfile <- file.path(outputdir, "config.Rdata") | ||
| envir <- new.env() | ||
| load(configfile, envir = envir) | ||
|
|
||
| # produce REMIND reporting *.mif based on gdx information | ||
| message("\n### start generation of mif files at ", Sys.time()) | ||
| convGDX2MIF(gdx, gdx_refpolicycost = gdx_refpolicycost, file = remind_reporting_file, | ||
| scenario = scenario, gdx_ref = gdx_ref) | ||
|
|
||
| # MAGICC code not working with REMIND-EU | ||
| # generate MAGICC reporting and append to REMIND reporting | ||
| if (0 == nchar(Sys.getenv('MAGICC_BINARY'))) { | ||
| warning('Can\'t find magicc executable under environment variable MAGICC_BINARY') | ||
| } else { | ||
| message("Generate ", basename(magicc_reporting_file)) | ||
| system(paste("cd ",outputdir ,"/magicc; ", | ||
| "pwd;", | ||
| "sed -f modify_MAGCFG_USER_CFG.sed -i MAGCFG_USER.CFG; ", | ||
| Sys.getenv('MAGICC_BINARY'), '; ', | ||
| "awk -f MAGICC_reporting.awk -v c_expname=\"", scenario, "\"", | ||
| " < climate_reporting_template.txt ", | ||
| " > ","../../../", magicc_reporting_file,"; ", | ||
| "sed -i 's/;glob;/;World;/g' ","../../../", magicc_reporting_file, "; ", | ||
| "cat ", "../../../",magicc_reporting_file, " >> ", "../../../",remind_reporting_file, "; ", | ||
| sep = "")) | ||
| } | ||
|
|
||
|
|
||
| ## generate EDGE-T reporting if it is needed | ||
| ## the reporting is appended to REMIND_generic_<scenario>.MIF | ||
| ## REMIND_generic_<scenario>_withoutPlus.MIF is replaced. | ||
|
|
||
| edgetOutputDir <- file.path(outputdir, "EDGE-T") | ||
| if(file.exists(edgetOutputDir)) { | ||
| if (! file.exists(file.path(edgetOutputDir, "demandF_plot_pkm.RDS"))) { | ||
| message("EDGE-T reporting files are missing, probably because the run was killed.") | ||
| message("Rerunning toolIterativeEDGETransport(reporting = TRUE).") | ||
| savewd <- getwd() | ||
| setwd(outputdir) | ||
| edgeTransport::toolIterativeEDGETransport(reporting = TRUE) | ||
| setwd(savewd) | ||
| } | ||
| message("start generation of EDGE-T reporting") | ||
| EDGET_output <- toolReportEDGET(edgetOutputDir, | ||
| extendedReporting = FALSE, | ||
| scenario_title = scenario, model_name = "REMIND", | ||
| gdx = file.path(outputdir, "fulldata.gdx")) | ||
|
|
||
| write.mif(EDGET_output, remind_reporting_file, append = TRUE) | ||
| deletePlus(remind_reporting_file, writemif = TRUE) | ||
|
|
||
| message("end generation of EDGE-T reporting") | ||
| } | ||
|
|
||
| magpie_reporting_file <- envir$cfg$pathToMagpieReport | ||
| if (! is.null(magpie_reporting_file) && file.exists(magpie_reporting_file)) { | ||
| message("add MAgPIE reporting from ", magpie_reporting_file) | ||
| tmp_rem <- quitte::as.quitte(remind_reporting_file) | ||
| tmp_mag <- dplyr::filter(quitte::as.quitte(magpie_reporting_file), .data$period %in% unique(tmp_rem$period)) | ||
| # remove population from magpie reporting to avoid duplication (units "million" vs. "million people") | ||
| sharedvariables <- intersect(tmp_mag$variable, tmp_rem$variable) | ||
| if (length(sharedvariables) > 0) { | ||
| message("The following variables will be dropped from MAgPIE reporting because they are in REMIND reporting: ", paste(sharedvariables, collapse = ", ")) | ||
| tmp_mag <- dplyr::filter(tmp_mag, ! .data$variable %in% sharedvariables) | ||
| } | ||
| # harmonize scenario name from -mag-xx to -rem-xx | ||
| tmp_mag$scenario <- paste0(scenario) | ||
| tmp_rem_mag <- rbind(tmp_rem, tmp_mag) | ||
| quitte::write.mif(tmp_rem_mag, path = remind_reporting_file) | ||
| deletePlus(remind_reporting_file, writemif = TRUE) | ||
| } | ||
|
|
||
| message("### end generation of mif files at ", Sys.time()) | ||
|
|
||
| ## produce REMIND LCOE reporting *.csv based on gdx information | ||
| message("### start generation of LCOE reporting at ", Sys.time()) | ||
| tmp <- try(convGDX2CSV_LCOE(gdx,file=LCOE_reporting_file,scen=scenario)) # execute convGDX2MIF_LCOE | ||
| message("### end generation of LCOE reporting at ", Sys.time()) | ||
|
|
||
| ## generate DIETER reporting if it is needed | ||
| ## the reporting is appended to REMIND_generic_<scenario>.MIF in "DIETER" Sub Directory | ||
| DIETERGDX <- "report_DIETER.gdx" | ||
| if(file.exists(file.path(outputdir, DIETERGDX))){ | ||
| message("start generation of DIETER reporting at ", Sys.time()) | ||
| remind2::reportDIETER(DIETERGDX,outputdir) | ||
| message("end generation of DIETER reporting at ", Sys.time()) | ||
| } | ||
|
|
||
| message("### reporting finished.") |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Some generated files are not rendered by default. Learn more about how customized files appear on GitHub.
Oops, something went wrong.