Merge pull request #69 from LaviniaBaumstark/develop

new realization of the module 45_carbonprice

core/postsolve.gms

@@ -332,15 +332,6 @@ if(cm_iterative_target_adj eq 9,
 *RP* Update tax levels/ multigasbudget values to reach the peak CO2 budget, with a linear increase afterwards given by cm_taxCO2inc_after_peakBudgYr
 *** The PeakBudgYr is found automatically by the algorithm (within the time window 204-2100)
  
-*** Save the original functional form of the CO2 price trajectory so values for all times can be accessed even if the peakBudgYr is shifted. 
-*** Then change to linear increasing CO2 price after peaking time 
-  if( iteration.val eq 1, 
-    p_taxCO2eq_until2150(t,regi) = pm_taxCO2eq(t,regi);
-	loop(t2$(t2.val eq cm_peakBudgYr),
-	  pm_taxCO2eq(t,regi)$(t.val gt cm_peakBudgYr) = p_taxCO2eq_until2150(t2,regi) + (t.val - t2.val) * cm_taxCO2inc_after_peakBudgYr * sm_DptCO2_2_TDpGtC;  !! increase by cm_taxCO2inc_after_peakBudgYr per year
-	);
-  );
-
   p_actualbudgetco2(t) =  sum(ttot$(ttot.val < t.val AND ttot.val > 2010), (sum(regi, (vm_emiTe.l(ttot,regi,"co2") + vm_emiCdr.l(ttot,regi,"co2") + vm_emiMac.l(ttot,regi,"co2"))) * sm_c_2_co2 * pm_ts(ttot)))
                           + sum(regi, (vm_emiTe.l(t,regi,"co2") + vm_emiCdr.l(t,regi,"co2") + vm_emiMac.l(t,regi,"co2")))*sm_c_2_co2 * (pm_ts(t) * 0.5 + 0.5)
                           + sum(regi, (vm_emiTe.l("2010",regi,"co2") + vm_emiCdr.l("2010",regi,"co2") + vm_emiMac.l("2010",regi,"co2")))*sm_c_2_co2 * 2;

core/preloop.gms

@@ -70,6 +70,19 @@ pm_vintage_in(regi,"1",te) = pm_vintage_in(regi,"1",te) * max((pm_histfegrowth(r
 pm_vintage_in(regi,"6",te) = pm_vintage_in(regi,"6",te) * max(((pm_histfegrowth(regi,entyFe)- 0.005 + 1/fm_dataglob("lifetime",te))/(1/fm_dataglob("lifetime",te)) + 1) * 0.75,0.2);
 );
 
+*RP
+*** First adjustment of CO2 price path for peakBudget runs (set by cm_iterative_target_adj eq 9)
+if(cm_iterative_target_adj eq 9,
+*** Save the original functional form of the CO2 price trajectory so values for all times can be accessed even if the peakBudgYr is shifted. 
+*** Then change to linear increasing CO2 price after peaking time 
+  p_taxCO2eq_until2150(t,regi) = pm_taxCO2eq(t,regi);
+  loop(t2$(t2.val eq cm_peakBudgYr),
+    pm_taxCO2eq(t,regi)$(t.val gt cm_peakBudgYr) = p_taxCO2eq_until2150(t2,regi) + (t.val - t2.val) * cm_taxCO2inc_after_peakBudgYr * sm_DptCO2_2_TDpGtC;  !! increase by cm_taxCO2inc_after_peakBudgYr per year
+  );
+);
+
+display p_taxCO2eq_until2150, pm_taxCO2eq;
+
 $ifthen setGlobal c_scaleEmiHistorical
 *re-scale MAgPie reference emissions to be inline with eurostat data (MagPie overestimates non-CO2 GHG emissions by a factor of 50% more)
 display p_macBaseMagpie;

core/sets.gms

@@ -782,7 +782,6 @@ $ENDIF.RegScenCapt
 ***######################## R SECTION START (MODULES) ###############################
 *** THIS CODE IS CREATED AUTOMATICALLY, DO NOT MODIFY THESE LINES DIRECTLY
 *** ANY DIRECT MODIFICATION WILL BE LOST AFTER NEXT MODEL START
-*** CHANGES CAN BE DONE USING THE RESPECTIVE LINES IN scripts/start_functions.R
 
 sets
 
@@ -823,7 +822,7 @@ sets
        codePerformance
        /
 
-module2realisation(modules,*) "mapping of modules and active realisations" /
+      module2realisation(modules,*) "mapping of modules and active realisations" /
        macro . %macro%
        welfare . %welfare%
        PE_FE_parameters . %PE_FE_parameters%

main.gms

@@ -313,7 +313,7 @@ cm_bioenergymaxscen = 0;         !! def = 0
 cm_tradecost_bio     = 2;         !! def = 2
 $setglobal cm_LU_emi_scen  SSP2   !! def = SSP2
 cm_1stgen_phaseout  = 0;         !! def = 0
-cm_cprice_red_factor  = 0.5;         !! def = 0.5
+cm_cprice_red_factor  = 1;         !! def = 1
 
 $setglobal cm_POPscen  pop_SSP2  !! def = pop_SSP2
 $setglobal cm_GDPscen  gdp_SSP2  !! def = gdp_SSP2

modules/45_carbonprice/diffCurvPhaseIn2Lin/datainput.gms

+*** |  (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
+*** SOF ./modules/45_carbonprice/diffCurvPhaseIn2Lin/datainput.gms
+***------------------------------------------------------------------------------------------------------------------------
+*** *BS* 20190930 linear convergence with starting points differentiated by GDP/capita, global price from 2040
+***-----------------------------------------------------------------------------------------------------------------------
+
+
+*** convergence to global CO2 price depends on GDP per capita (in 1e3 $ PPP 2005).
+p45_gdppcap2015_PPP(regi) = pm_gdp("2015",regi)/pm_shPPPMER(regi) / pm_pop("2015",regi);
+display p45_gdppcap2015_PPP;
+
+p45_phasein_2025ratio(regi)$(p45_gdppcap2015_PPP(regi) le 3) = 0.1;
+p45_phasein_2025ratio(regi)$(p45_gdppcap2015_PPP(regi) gt 3 and p45_gdppcap2015_PPP(regi) le 5) = 0.2;
+p45_phasein_2025ratio(regi)$(p45_gdppcap2015_PPP(regi) gt 5 and p45_gdppcap2015_PPP(regi) le 8) = 0.3;
+p45_phasein_2025ratio(regi)$(p45_gdppcap2015_PPP(regi) gt 8 and p45_gdppcap2015_PPP(regi) le 11) = 0.5;
+p45_phasein_2025ratio(regi)$(p45_gdppcap2015_PPP(regi) gt 11 and p45_gdppcap2015_PPP(regi) le 14) = 0.65;
+p45_phasein_2025ratio(regi)$(p45_gdppcap2015_PPP(regi) gt 14 and p45_gdppcap2015_PPP(regi) le 19) = 0.8;
+p45_phasein_2025ratio(regi)$(p45_gdppcap2015_PPP(regi) gt 19 and p45_gdppcap2015_PPP(regi) le 24) = 0.9;
+p45_phasein_2025ratio(regi)$(p45_gdppcap2015_PPP(regi) gt 24) = 1;
+display p45_phasein_2025ratio;
+
+
+*** for the current implementation, use the following trajectory for rich countries:
+*** global price is linear from 2010 until the pkBudgYr, then increases with cm_taxCO2inc_after_peakBudgYr
+if(cm_co2_tax_2020 lt 0,
+  abort "please choose a valid cm_co2_tax_2020"
+elseif cm_co2_tax_2020 ge 0,
+*** convert tax value from $/t CO2eq to T$/GtC
+  p45_CO2priceTrajDeveloped("2040")= 3 * cm_co2_tax_2020 * sm_DptCO2_2_TDpGtC;  !! shifted to 2040 to make sure that even in delay scenarios the fixpoint of the linear price path is inside the "t" range, otherwise the CO2 prices from reference run may be overwritten
+*** The factor 3 comes from shifting the 2020 value 20 years into the future at linear increase of 10% of 2020 value per year.
+);
+
+
+
+p45_CO2priceTrajDeveloped(ttot)$((ttot.val gt 2005) AND (ttot.val ge cm_startyear)) = p45_CO2priceTrajDeveloped("2040")*( 1 + 0.1/3 * (ttot.val-2040)); !! no CO2 price in 2005 and only change CO2 prices after ; 
+*** annual increase by (10/3)% of the 2040 value is the same as a 10% increase of the 2020 value is the same as a linear increase from 0 in 2010 to the 2020/2040 value
+
+
+*** Then create regional phase-in:
+loop(ttot$((ttot.val ge cm_startyear) AND (ttot.val le cm_CO2priceRegConvEndYr) ),
+  p45_regCO2priceFactor(ttot,regi) = 
+   min(1,
+       max(0, 
+	        p45_phasein_2025ratio(regi) + (1 - p45_phasein_2025ratio(regi)) 
+			                               * ( ( (ttot.val - 2025) + (cm_CO2priceRegConvEndYr - 2025) * 0.1 ) 
+                                               / ( (cm_CO2priceRegConvEndYr - 2025) * 1.1 ) 
+											 ) ** 2
+       )											 
+   );
+);
+p45_regCO2priceFactor(ttot,regi)$(ttot.val ge cm_CO2priceRegConvEndYr) = 1;
+
+
+*** transition to global price - starting point depends on GDP/cap
+pm_taxCO2eq(t,regi) = p45_regCO2priceFactor(t,regi) * p45_CO2priceTrajDeveloped(t);
+
+
+display p45_regCO2priceFactor, p45_CO2priceTrajDeveloped, pm_taxCO2eq;
+
+*** EOF ./modules/45_carbonprice/diffCurvPhaseIn2Lin/datainput.gms

modules/45_carbonprice/diffCurvPhaseIn2Lin/declarations.gms

+*** |  (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
+*** SOF ./modules/45_carbonprice/diffCurvPhaseIn2Lin/declarations.gms
+***------------------------------------------------------------------------------------------------------------------------
+*** *BS* 20190930 linear convergence with starting points differentiated by GDP/capita, global price from 2040
+***-----------------------------------------------------------------------------------------------------------------------
+
+parameters
+p45_tauCO2_ref(ttot, all_regi)              "CO2 tax path of reference policy (NDC)"
+p45_gdppcap2015_PPP(all_regi)               "2015 GDP per capita (k $ PPP 2005)"
+p45_phasein_2025ratio(all_regi)             "ratio of CO2 price to that of developed region in 2025"
+
+p45_regCO2priceFactor(ttot,all_regi)                    "regional multiplicative factor to the CO2 price of the developed countries"
+p45_CO2priceTrajDeveloped(ttot)                         "CO2 price trajectory for developed/rich countries"
+;
+
+
+*** EOF ./modules/45_carbonprice/diffCurvPhaseIn2Lin/declarations.gms

modules/45_carbonprice/diffCurvPhaseIn2Lin/not_used.txt

+name,type,reason
+sm_c_2_co2,input,questionnaire
+cm_peakBudgYr,input,questionnaire
+vm_co2eq,input,questionnaire
+vm_emiFgas,input,questionnaire
+vm_prodSe,input,questionnaire
+vm_prodFe,input,questionnaire
+pm_pvp,input,questionnaire
+pm_globalMeanTemperature,input,questionnaire
+pm_temperatureImpulseResponseCO2,input,questionnaire
+pm_GDPGross,input,questionnaire
+pm_ttot_val,input,questionnaire
+pm_ts,input,questionnaire
+pm_ttot_2_tall,input,questionnaire
+cm_emiscen,input,questionnaire
+cm_co2_tax_growth,input,questionnaire
+cm_iterative_target_adj,input,questionnaire
+cm_expoLinear_yearStart,input,questionnaire
+cm_carbonprice_temperatureLimit,input,questionnaire
+cm_taxCO2inc_after_peakBudgYr,input,questionnaire
+vm_cesIO,input,questionnaire
+pm_consPC,input,questionnaire
+pm_prtp,input,questionnaire
+pm_prodCouple,input,questionnaire
+pm_bunker_share_in_nonldv_fe,input,questionnaire

modules/45_carbonprice/diffCurvPhaseIn2Lin/postsolve.gms

+*** |  (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
+*** SOF ./modules/45_carbonprice/diffCurvPhaseIn2Lin/postsolve.gms
+***------------------------------------------------------------------------------------------------------------------------
+*** *BS* 20190930 linear convergence with starting points differentiated by GDP/capita, global price from 2040
+***-----------------------------------------------------------------------------------------------------------------------
+
+***   *** updated constant global price as scalar (regional prices are the same anyway)
+***   s45_constantCO2price = sum((ttot,regi)$(ttot.val eq s45_stageend), pm_taxCO2eq(ttot,regi))/card(regi) ;
+***   *** entire path has been shifted in update, so have to set these again
+***   pm_taxCO2eq(ttot,regi)$(ttot.val le s45_stagestart) = p45_tauCO2_ref(ttot, regi);
+***   pm_taxCO2eq(ttot,regi)$(ttot.val gt s45_stagestart and ttot.val lt s45_stageend)
+***     = s45_constantCO2price * (ttot.val - p45_phasein_zeroyear(regi))/(s45_stageend - p45_phasein_zeroyear(regi));
+***   *** price trajectory should be constant anyway but let's be explicit here
+***   pm_taxCO2eq(ttot,regi)$(ttot.val ge s45_stageend) = s45_constantCO2price;
+
+
+*** re-create the regional differentation, use path from developed countries as the basis.
+*** This doesn't need to be a loop, but it will be correct for any cycle of the loop, so also for the last cycle.
+loop(regi$(p45_gdppcap2015_PPP(regi) gt 30),
+  p45_CO2priceTrajDeveloped(t) = pm_taxCO2eq(t,regi);
+);
+
+*** linear transition to global price - starting point depends on GDP/cap
+pm_taxCO2eq(t,regi) = p45_regCO2priceFactor(t,regi) * p45_CO2priceTrajDeveloped(t);
+
+display pm_taxCO2eq;
+*** EOF ./modules/45_carbonprice/diffCurvPhaseIn2Lin/postsolve.gms

modules/45_carbonprice/diffCurvPhaseIn2Lin/realization.gms

+*** |  (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
+
+*####################### R SECTION START (PHASES) ##############################
+$Ifi "%phase%" == "declarations" $include "./modules/45_carbonprice/diffCurvPhaseIn2Lin/declarations.gms"
+$Ifi "%phase%" == "datainput" $include "./modules/45_carbonprice/diffCurvPhaseIn2Lin/datainput.gms"
+$Ifi "%phase%" == "postsolve" $include "./modules/45_carbonprice/diffCurvPhaseIn2Lin/postsolve.gms"
+*######################## R SECTION END (PHASES) ###############################

modules/45_carbonprice/module.gms

@@ -20,6 +20,7 @@ $Ifi "%carbonprice%" == "ExogSameAsPrevious" $include "./modules/45_carbonprice/
 $Ifi "%carbonprice%" == "NDC2018" $include "./modules/45_carbonprice/NDC2018/realization.gms"
 $Ifi "%carbonprice%" == "NDC2constant" $include "./modules/45_carbonprice/NDC2constant/realization.gms"
 $Ifi "%carbonprice%" == "NPi2018" $include "./modules/45_carbonprice/NPi2018/realization.gms"
+$Ifi "%carbonprice%" == "diffCurvPhaseIn2Lin" $include "./modules/45_carbonprice/diffCurvPhaseIn2Lin/realization.gms"
 $Ifi "%carbonprice%" == "diffPhaseIn2Constant" $include "./modules/45_carbonprice/diffPhaseIn2Constant/realization.gms"
 $Ifi "%carbonprice%" == "diffPhaseIn2Lin" $include "./modules/45_carbonprice/diffPhaseIn2Lin/realization.gms"
 $Ifi "%carbonprice%" == "diffPhaseIn2LinFlex" $include "./modules/45_carbonprice/diffPhaseIn2LinFlex/realization.gms"