remindmodel · gabriel-abrahao · May 27, 2024 · Sep 7, 2023 · Mar 19, 2024 · Mar 19, 2024
diff --git a/.gitignore b/.gitignore
@@ -19,6 +19,7 @@
 
 # Ignore everything in "input" folders
 input/
+!/scripts/input/
 
 # Ignore everything in root "output" folder
 /output/

diff --git a/config/default.cfg b/config/default.cfg
@@ -49,11 +49,12 @@ cfg$runstatistics <- "/p/projects/rd3mod/models/statistics/remind"
 # copied into remind when starting a run and required to calculate climate responses (e.g. temperature).
 cfg$magicc_template <- "/p/projects/rd3mod/magicc/"
 
-#### Folder containing files for climate-assessment that cannot be published with REMIND
-# These files currently all have permissive licenses and can be easily found online
-cfg$climate_assessment_root       <- "/p/projects/rd3mod/python/climate-assessment/src/"
-cfg$climate_assessment_files_dir  <- "/p/projects/rd3mod/climate-assessment-files/"
-cfg$climate_assessment_magicc_bin <- "/p/projects/rd3mod/climate-assessment-files/magicc-v7.5.3/bin/magicc"
+#### Folders and file paths used for climate-assessment/MAGICC7 integration & reporting
+cfg$climate_assessment_root        <- "/p/projects/rd3mod/python/climate-assessment/src/"
+cfg$climate_assessment_infiller_db <- "/p/projects/rd3mod/climate-assessment-files/1652361598937-ar6_emissions_vetted_infillerdatabase_10.5281-zenodo.6390768.csv"
+cfg$climate_assessment_magicc_bin  <- "/p/projects/rd3mod/climate-assessment-files/magicc-v7.5.3/bin/magicc"
+cfg$climate_assessment_magicc_prob_file_iteration <- "/p/projects/rd3mod/climate-assessment-files/parsets/RCP20_50.json"
+cfg$climate_assessment_magicc_prob_file_reporting <- "/p/projects/rd3mod/climate-assessment-files/parsets/0fd0f62-derived-metrics-id-f023edb-drawnset.json"
 
 #### Output folder of the modeltests generated by running config/scenario_config_AMT.csv
 # start.R searches in this folder for finished runs to be used as path_gdx...
@@ -132,6 +133,8 @@ cfg$files2export$start <- c("config/conopt3.opt",
                             "config/conopt3.op9",
                             "core/input/historical/historical.mif",
                             "scripts/input/exoGAINSAirpollutants.R",
+                            "scripts/input/climate_assessment_run.R",
+                            "scripts/input/climate_assessment_temperatureImpulseResponse.R",
                             ".Rprofile",
                             "config/mappingEDGEtoREMINDsectors.csv",
                             "modules/11_aerosols/exoGAINS/input/ef_gains.cs4r",

diff --git a/config/scenario_config_magicc.csv b/config/scenario_config_magicc.csv
@@ -0,0 +1,9 @@
+title;start;CES_parameters;optimization;slurmConfig;climate;regionmapping;extramappings_historic;cm_rcp_scen;cm_iterative_target_adj;subsidizeLearning;c_budgetCO2from2020;carbonprice;cm_co2_tax_2020;c_peakBudgYr;cm_CO2priceRegConvEndYr;cm_emiscen;c_regi_earlyreti_rate;c_tech_earlyreti_rate;cm_fetaxscen;cm_co2_tax_growth;cm_bioenergy_SustTax;cm_maxProdBiolc;c_ccsinjecratescen;c_ccscapratescen;cm_CESMkup_build;cm_CESMkup_ind;techpol;regipol;cm_implicitQttyTarget;cm_emiMktTarget;cm_NucRegiPol;cm_CoalRegiPol;cm_altFeEmiFac;cm_POPscen;cm_GDPscen;cm_demScen;cm_oil_scen;cm_gas_scen;cm_coal_scen;c_techAssumptScen;cm_nucscen;cm_so2tax_scen;cm_multigasscen;cm_LU_emi_scen;cm_tradecostBio;cm_1stgen_phaseout;c_SSP_forcing_adjust;cm_APscen;water;cm_EDGEtr_scen;cm_startyear;path_gdx;path_gdx_ref;path_gdx_bau;path_gdx_refpolicycost;description
+SSP2EU-NPi-calibrate;0;calibrate;;14;magicc;;;rcp45;;;0;NPi;1;2100;;9;;;;;;;;;;;NPi2018;regiCarbonPrice;;;on;on;EUR_regi, NEU_regi;;;;;;;;;;3;;;;;;heat;Mix1;2005;;;;;SSP2EU-NPi-calibrate: This reference policy/baseline scenario follows the Shared Socioeconomic Pathways 2 called Middle of the Road.
+SSP2EU-Base-magicc;0;;;;magicc;;;;;;0;;;2100;;;;off;;;;;;;;;;;;;;;EUR_regi, NEU_regi;;;;;;;;;;;;;;;;heat;Mix1;2005;;;;;SSP2EU-Base: This baseline calibration scenario follows the Shared Socioeconomic Pathways 2 called Middle of the Road.
+SSP2EU-NDC-magicc;0;;;;magicc;;;rcp45;3;globallyOptimal;0;NDC;1;2100;;9;;;;;;;;;;;NDC;regiCarbonPrice;2030.EUR_regi.tax.t.FE_wo_b_wo_n_e.all 1.2809;2020.2030.EUR_regi.all.year.netGHG_LULUCFGrassi 2.450;on;on;EUR_regi, NEU_regi;;;;;;;;;;3;;;;;;heat;Mix2;2025;;SSP2EU-NPi;SSP2EU-NPi;;SSP2EU-NDC: This Nationally Determined Contribution (NDC) scenario follows the Shared Socioeconomic Pathways 2 called Middle of the Road. The NDC includes all pledged policies even if not yet implemented. It assumes that the moderate and heterogeneous climate ambition reflected in the NDCs at the begining of 2021 continues over the 21st century. This scenario serves as reference for all other policy scenarios and The NPi scenario until 2020.
+SSP2EU-NPi-magicc;0;;;;magicc;;;rcp45;3;;0;NPi;1;2100;;9;;;;;;;;;;;NPi2018;regiCarbonPrice;;;on;on;EUR_regi, NEU_regi;;;;;;;;;;3;;;;;;heat;Mix1;2005;;;;;SSP2EU-NPi: This National Policies Implemented (Npi) scenario follows the Shared Socioeconomic Pathways 2 called Middle of the Road. The NPi is identical to the NDC scenario until 2020 but assumes that policies fail to achieve NDC targets in 2030. Instead, carbon prices are assumed to grow and converge more slowly, leading to emissions trajectories in line with bottom-up studies on the effect of currently implemented policies.
+SSP2EU-NPi-ar6;1;;;;ar6;;;rcp45;3;;0;NPi;1;2100;;9;;;;;;;;;;;NPi2018;regiCarbonPrice;;;on;on;EUR_regi, NEU_regi;;;;;;;;;;3;;;;;;heat;Mix1;2005;;;;;SSP2EU-NPi: This National Policies Implemented (Npi) scenario follows the Shared Socioeconomic Pathways 2 called Middle of the Road. The NPi is identical to the NDC scenario until 2020 but assumes that policies fail to achieve NDC targets in 2030. Instead, carbon prices are assumed to grow and converge more slowly, leading to emissions trajectories in line with bottom-up studies on the effect of currently implemented policies.
+SSP2EU-PkBudg500-magicc;0;;;;magicc;;;rcp20;9;globallyOptimal;500;diffCurvPhaseIn2Lin;100;2080;;9;;;;;;;;;;;NDC;regiCarbonPrice;2030.EUR_regi.tax.t.FE_wo_b_wo_n_e.all 1.2809;;on;on;EUR_regi, NEU_regi;;;;;;;;;;;;;;;;heat;Mix4;2025;;SSP2EU-NPi-magicc;SSP2EU-NPi-magicc;;SSP2EU-PkBudg500: This climate policy scenario follows the Shared Socioeconomic Pathways 2 called Middle of the Road. The stylized climate policy scenario assumes a peak budget of 500 Gt?CO2 on total CO2 emissions from 2015 to 2100. This is a 1.5??C scenario, peak warming is allowed to be at or slightly above 1.5??C, at median climate sensitivity but returns to values below 1.5??C in at least 67?% of scenarios by the end of the century.
+SSP2EU-PkBudg1150-magicc;0;;;;magicc;;;rcp26;9;globallyOptimal;1150;diffCurvPhaseIn2Lin;60;2100;;9;;;;;;;;;;;NDC;regiCarbonPrice;2030.EUR_regi.tax.t.FE_wo_b_wo_n_e.all 1.2809;;on;on;EUR_regi, NEU_regi;;;;;;;;;;;;;;;;heat;Mix3;2025;;SSP2EU-NPi-magicc;SSP2EU-NPi-magicc;;SSP2EU-PkBudg1150: This climate policy scenario follows the Shared Socioeconomic Pathways 2 called Middle of the Road. The stylized climate policy scenario assumes a peak budget of 1150 Gt?CO2 on total CO2 emissions from 2015 to 2100. This is a weel below 2??C scenario, peak warming is allowed to be at or slightly above 2??C, at median climate sensitivity but returns to values below 2??C in at least 67?% of scenarios by the end of the century.
+SSP2EU-PkBudg500-ar6;0;;;;ar6;;;rcp20;9;globallyOptimal;500;diffCurvPhaseIn2Lin;100;2080;;9;;;;;;;;;;;NDC;regiCarbonPrice;2030.EUR_regi.tax.t.FE_wo_b_wo_n_e.all 1.2809;;on;on;EUR_regi, NEU_regi;;;;;;;;;;;;;;;;heat;Mix4;2025;;SSP2EU-NPi-magicc;SSP2EU-NPi-magicc;;SSP2EU-PkBudg500: This climate policy scenario follows the Shared Socioeconomic Pathways 2 called Middle of the Road. The stylized climate policy scenario assumes a peak budget of 500 Gt?CO2 on total CO2 emissions from 2015 to 2100. This is a 1.5??C scenario, peak warming is allowed to be at or slightly above 1.5??C, at median climate sensitivity but returns to values below 1.5??C in at least 67?% of scenarios by the end of the century.
diff --git a/modules/15_climate/ar6/datainput.gms b/modules/15_climate/ar6/datainput.gms
@@ -0,0 +1,84 @@
+*** |  (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]
+*** SOF ./modules/15_climate/magicc/datainput.gms
+*** cluster rcp_scen into overshoot and not-to-exceed targets
+$if %cm_rcp_scen% == "none"    s15_rcpCluster = 1;
+$if %cm_rcp_scen% == "rcp20"   s15_rcpCluster = 1;
+$if %cm_rcp_scen% == "rcp26"   s15_rcpCluster = 1;
+$if %cm_rcp_scen% == "rcp37"   s15_rcpCluster = 1;
+$if %cm_rcp_scen% == "rcp45"   s15_rcpCluster = 0;
+$if %cm_rcp_scen% == "rcp60"   s15_rcpCluster = 0;
+$if %cm_rcp_scen% == "rcp85"   s15_rcpCluster = 0;
+
+
+s15_forcing_budgetiterationoffset = 1.5;
+$if %cm_rcp_scen% == "rcp20" s15_forcing_budgetiterationoffset = 1.2;
+s15_forcing_budgetiterationoffset_tax = 0.0;
+
+*JeS* Forcing target is now on RCP forcing instead of total forcing. The 
+*** difference is about 0.4 W/m^2, therefore the target on the RCP forcing 
+*** has to be 0.4 higher than the intended total forcing.
+$if %cm_rcp_scen% == "none"   s15_gr_forc_nte = 100;
+$if %cm_rcp_scen% == "rcp26"  s15_gr_forc_nte = 10;
+$if %cm_rcp_scen% == "rcp37"  s15_gr_forc_nte = 3.67;
+$if %cm_rcp_scen% == "rcp45"  s15_gr_forc_nte = 4.22;
+$if %cm_rcp_scen% == "rcp60"  s15_gr_forc_nte = 5.44;
+$if %cm_rcp_scen% == "rcp85"  s15_gr_forc_nte = 8.5;
+$if %cm_rcp_scen% == "rcp20"  s15_gr_forc_nte = 10;
+
+$if %cm_rcp_scen% == "none"   s15_gr_forc_os = 100;
+$if %cm_rcp_scen% == "rcp26"  s15_gr_forc_os = 2.54;
+$if %cm_rcp_scen% == "rcp37"  s15_gr_forc_os = 3.67;
+$if %cm_rcp_scen% == "rcp45"  s15_gr_forc_os = 4.23;
+$if %cm_rcp_scen% == "rcp60"  s15_gr_forc_os = 5.44;
+$if %cm_rcp_scen% == "rcp85"  s15_gr_forc_os = 8.5;
+$if %cm_rcp_scen% == "rcp20"  s15_gr_forc_os = 1.93;
+
+s15_gr_forc_kyo_nte = 0;
+s15_gr_forc_kyo     = 0;
+
+p15_forc_magicc(tall) = 0;
+
+$if  %cm_rcp_scen% == "rcp26"    $include "./modules/15_climate/magicc/input/data_oghgf_rcp3pd.inc";
+$if  %cm_rcp_scen% == "rcp37"    $include "./modules/15_climate/magicc/input/data_oghgf_rcp45.inc";
+$if  %cm_rcp_scen% == "rcp45"    $include "./modules/15_climate/magicc/input/data_oghgf_rcp45.inc";
+$if  %cm_rcp_scen% == "rcp60"    $include "./modules/15_climate/magicc/input/data_oghgf_rcp6.inc";
+$if  %cm_rcp_scen% == "rcp85"    $include "./modules/15_climate/magicc/input/data_oghgf_rcp85.inc";
+$if  %cm_rcp_scen% == "rcp20"    $include "./modules/15_climate/magicc/input/data_oghgf_rcp3pd.inc";
+$if  %cm_rcp_scen% == "none"     $include "./modules/15_climate/magicc/input/data_oghgf_rcp6.inc";
+
+pm_emicapglob(ttot) = 0;
+
+*** parameter pm_emicapglob is read in depending on cm_rcp_scen and cm_multigasscen
+if( (cm_multigasscen = 1) or (cm_multigasscen = 3),
+$offlisting
+$if %cm_rcp_scen% == "rcp20"  $include "./modules/15_climate/off/input/pm_emicapglob_450.inc";
+$if %cm_rcp_scen% == "rcp26"  $include "./modules/15_climate/off/input/pm_emicapglob_450.inc";
+$if %cm_rcp_scen% == "rcp37"  $include "./modules/15_climate/off/input/pm_emicapglob_550.inc";
+$if %cm_rcp_scen% == "rcp45"  $include "./modules/15_climate/off/input/pm_emicapglob_550.inc";
+$if %cm_rcp_scen% == "rcp60"  $include "./modules/15_climate/off/input/pm_emicapglob.inc";
+$if %cm_rcp_scen% == "rcp85"  $include "./modules/15_climate/off/input/pm_emicapglob.inc";
+$if %cm_rcp_scen% == "none"   $include "./modules/15_climate/off/input/pm_emicapglob.inc";
+$onlisting
+);
+if(cm_multigasscen = 2,
+$offlisting
+$if %cm_rcp_scen% == "rcp20"  $include "./modules/15_climate/off/input/pm_emicapglob_multigas_450.inc";
+$if %cm_rcp_scen% == "rcp26"  $include "./modules/15_climate/off/input/pm_emicapglob_multigas_450.inc";
+$if %cm_rcp_scen% == "rcp37"  $include "./modules/15_climate/off/input/pm_emicapglob_multigas_550.inc";
+$if %cm_rcp_scen% == "rcp45"  $include "./modules/15_climate/off/input/pm_emicapglob_multigas_550.inc";
+$if %cm_rcp_scen% == "rcp60"  $include "./modules/15_climate/off/input/pm_emicapglob_multigas.inc";
+$if %cm_rcp_scen% == "rcp85"  $include "./modules/15_climate/off/input/pm_emicapglob_multigas.inc";
+$if %cm_rcp_scen% == "none"   $include "./modules/15_climate/off/input/pm_emicapglob_multigas.inc";
+$onlisting
+);
+display pm_emicapglob;
+
+pm_gmt_conv=1;
+p15_gmt0(tall)=1;
+
+*** EOF ./modules/15_climate/magicc/datainput.gms
diff --git a/modules/15_climate/ar6/declarations.gms b/modules/15_climate/ar6/declarations.gms
@@ -0,0 +1,47 @@
+*** |  (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]
+*** SOF ./modules/15_climate/ar6/declarations.gms
+
+Scalar
+  s15_forcing_budgetiterationoffset "offset for the calculation of iteratively adjusted budget"
+*** to reach specified forcing level in 2100 in W/m^2 -- best guess is 2010 
+*** value of respective forcing category minus 0.7
+  s15_forcing_budgetiterationoffset_tax "offset for the calculation of iteratively adjusted budget"
+*** to reach specified forcing level in 2100 in W/m^2 -- best guess taken from analysis of ERL paper results
+s15_tempOffset2010 "mean temperature in 2010 from MAGICC"
+;
+
+parameters
+p15_oghgf_pfc    	"exogenous forcings from RCP all in W/m^2: PFCs",
+p15_oghgf_hfc	 	"exogenous forcings from RCP: HFCs",
+p15_oghgf_sf6		"exogenous forcings from RCP: SF6",
+p15_oghgf_montreal	"exogenous forcings from RCP: montreal gases",
+p15_oghgf_o3str		"exogenous forcings from RCP: stratospheric ozone",
+p15_oghgf_luc		"exogenous forcings from RCP: albedo change due to land-use change",
+p15_oghgf_crbbb		"exogenous forcings from RCP: carbonaceous aerosols from biomass burning",
+p15_oghgf_ffbc		"exogenous forcings from RCP: black carbon from fossil fuels",
+p15_oghgf_ffoc		"exogenous forcings from RCP: organic carbon from fossil fuels",
+p15_oghgf_o3trp		"exogenous forcings from RCP: tropospheric ozone",
+p15_oghgf_h2ostr	"exogenous forcings from RCP: stratospheric water vapor",
+p15_oghgf_minaer	"exogenous forcings from RCP: mineral dust",
+p15_oghgf_nitaer	"exogenous forcings from RCP: nitrates",
+p15_emicapregi(tall,all_regi)  "regional emission caps, used for calculation of global emission cap",
+p15_forc_magicc(tall) "actual radiative forcing as calculated by magicc [W/m^2]"
+pm_gmt_conv          "global mean temperature convergence"    
+p15_gmt0(tall)        "global mean temperature convergence saved for the next iteration"
+;
+
+scalars
+s15_gr_forc_kyo      "guardrail for 450 ppm Kyoto forcing, adapted between negishi iterations - dummy parameter, only needed to prevent gdx errors",
+s15_gr_forc_kyo_nte  "guardrail for 550 ppm Kyoto forcing, adapted between negishi iterations - dummy parameter, only needed to prevent gdx errors",
+
+s15_gr_forc_nte      "not to exceed radiative forcing target from 2100 on [W/m^2]",
+s15_gr_forc_os       "overshoot (os) radiative forcing target (guardrail = gr) from 2100 on [W/m^2]"
+s15_rcpCluster       " clustering of rcp_scen for the iterative traget adjustment"
+;
+
+*** EOF ./modules/15_climate/ar6/declarations.gms
diff --git a/modules/15_climate/ar6/not_used.txt b/modules/15_climate/ar6/not_used.txt
@@ -0,0 +1,7 @@
+# |  (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]
+name, type, reason