Merge pull request #196 from mleimbach/gdp

Code changes for covid GDP correction

config/default.cfg

@@ -296,6 +296,7 @@ cfg$gms$cm_cprice_red_factor <- 1  	    # def <- 1
 
 cfg$gms$cm_POPscen            <- "pop_SSP2"     # def <- "pop_SSP2"
 cfg$gms$cm_GDPscen            <- "gdp_SSP2"     # def <- "gdp_SSP2"
+cfg$gms$cm_GDPcovid           <- 0       # def <- 0
 
 cfg$gms$cm_startyear          <- 2005    # def <- 2005 for BAU, 2015 for policy run
 cfg$gms$c_start_budget        <- 2100    # def <- 2100
@@ -366,7 +367,6 @@ cfg$gms$c_regi_nucscen <- "all" # def <- "all"
 cfg$gms$c_regi_capturescen <- "all" # def <- "all"
 
 cfg$gms$cm_EDGEtr_scen       <- "ConvCase"     # def <- "ConvCase". For calibration runs: "ConvCase", "ConvCaseWise". For standard runs:"ConvCase", "ElecEra", "HydrHype", "ConvCaseWise", "ElecEraWise", "HydrHypeWise"
-
 #-----------------------------------------------------------------------------
 #                          FLAGS
 #-----------------------------------------------------------------------------
@@ -579,6 +579,10 @@ cfg$RunsUsingTHISgdxAsBAU <- NA
 #  (gdp_SSP3):  SSP3 slowGROWTH slowCONV
 #  (gdp_SSP4):  SSP4  medGROWTH mixedCONV
 #  (gdp_SSP5):  SSP5 fastGROWTH fastCONV
+# cm_GDPcovid  "switch to turn on short-term GDP loss by covid-19"
+# *ML* emulates a schock, only feasible with start year 2020, don't use in calibration
+#  (0):  off
+#  (1):  on
 # cm_startyear "first free (or pseudo-free in case of 2005) modelling time step";
 #  (2005): standard for BAU to check if model is well calibrated
 #  (2015): standard for all policy runs (eq. to fix2010) and production baselines, especially if various baselines with varying parameters are explored
@@ -735,6 +739,7 @@ cfg$RunsUsingTHISgdxAsBAU <- NA
 # cfg$gms$c_skip_output            <- "off"    # def <- off
 # Skip output? (Useful for development)
 # cfg$gms$cm_ccsfosall             <- "off"    # def <- off
+
 ############################################################
 # here you find old and new names of the latest renaming
 #-----------------------------------------------------------

config/scenario_config_SSP2covid.csv

+title;start;cm_rcp_scen;cm_iterative_target_adj;subsidizeLearning;c_budgetCO2;carbonprice;cm_co2_tax_2020;cm_peakBudgYr;cm_emiscen;techpol;cm_OILRETIRE;cm_GDPcovid;cm_startyear;path_gdx;path_gdx_ref;path_gdx_bau
+SSP2-Base;1;none;0;off;0;none;-1;2100;1;none;off;0;2005;./config/input.gdx;;
+SSP2-NDC;1;rcp45;3;globallyOptimal;0;NDC2018;1;2100;9;NDC2018;off;0;2015;./config/input.gdx;SSP2-Base;SSP2-Base
+SSP2-NDC_covid;1;rcp45;3;globallyOptimal;0;NDC2018;1;2100;9;NDC2018;off;1;2020;./config/input.gdx;SSP2-NDC;SSP2-Base
+SSP2-Npi_covid;1;rcp45;3;off;0;NPi2018;1;2100;9;NPi2018;off;1;2025;./config/input.gdx;SSP2-NDC_covid;SSP2-Base
+SSP2-PkBudg900_covid;1;rcp20;9;globallyOptimal;900;diffCurvPhaseIn2Lin;200;2045;9;NDC2018;on;1;2025;./config/input.gdx;SSP2-NDC_covid;SSP2-Base
+SSP2-PkBudg1100_covid;1;rcp26;9;globallyOptimal;1100;diffCurvPhaseIn2Lin;100;2080;9;NDC2018;on;1;2025;./config/input.gdx;SSP2-NDC_covid;SSP2-Base
+SSP2-PkBudg1300_covid;1;rcp26;9;globallyOptimal;1300;diffCurvPhaseIn2Lin;60;2100;9;NDC2018;on;1;2025;./config/input.gdx;SSP2-NDC_covid;SSP2-Base

main.gms

@@ -290,6 +290,7 @@ cm_taxCO2inc_after_peakBudgYr "annual increase of CO2 price after the Peak Budge
 cm_CO2priceRegConvEndYr      "Year at which regional CO2 prices converge in module 45 realization diffPhaseIn2LinFlex"
 c_regi_nucscen				"regions to apply nucscen to"
 c_regi_capturescen			"region to apply ccapturescen to"
+cm_GDPcovid                  "GDP correction for covid"
 ;
 
 *** --------------------------------------------------------------------------------------------------------------------------------------------------------------------
@@ -334,6 +335,7 @@ cm_cprice_red_factor  = 1;         !! def = 1
 $setglobal cm_POPscen  pop_SSP2  !! def = pop_SSP2
 $setglobal cm_GDPscen  gdp_SSP2  !! def = gdp_SSP2
 $setglobal c_GDPpcScen  SSP2     !! def = gdp_SSP2   (automatically adjusted by start_run() based on GDPscen) 
+cm_GDPcovid      = 0;            !! def = 0
 
 *AG* and *CB* for cm_startyear greater than 2005, you have to copy the fulldata.gdx (rename it to: input_ref.gdx) from the run you want to build your new run onto.
 cm_startyear      = 2005;      !! def = 2005 for a BAU, 2015 for policy runs

modules/29_CES_parameters/calibrate/not_used.txt

+# |  (C) 2006-2019 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: remind@pik-potsdam.de
+name,type,reason
+cm_GDPcovid, switch, ???
+

modules/29_CES_parameters/load/datainput.gms

@@ -11,7 +11,9 @@
 $include "./modules/29_CES_parameters/load/input/stat_off-indu_fixed_shares-buil_simple-tran_complex-POP_pop_SSP2-GDP_gdp_SSP2-Kap_debt_limit-Reg_690d3718e1.inc"
 *###################### R SECTION END (CES INPUT) ###########################
 
-
+if (cm_GDPcovid eq 1,
+   pm_cesdata("2020",all_regi,"lab","effgr") = 0.5 * (pm_cesdata ("2015",all_regi,"lab","effgr") + pm_cesdata ("2020",all_regi,"lab","effgr"));
+);
 
 option pm_cesdata:8:3:1;
 display "loaded pm_cesdata", pm_cesdata;