Support for new switches to represent trade scenarios with limited cooperation between regions

CHANGELOG.md

@@ -10,6 +10,8 @@ The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/).
 ### added
 - **config** regex tests for many parameters
 - **testthat** test and compile all config files
+- **22_subsidizeLearning** option to define independent learning cumulative capacity for alternative region groups (`cm_altLearnRegiSet`).
+- **24_trade** option to limit specific regions exports (`cm_XportRegiLim`).
 
 ### fixed
 - fixed weights of energy carriers in `pm_IndstCO2Captured`

main.gms

@@ -1240,6 +1240,18 @@ $setGlobal cm_gas_scen  medGas         !! def = medGas  !! regexp = low|medium|h
 ***  (medCoal): medium
 ***  (highCoal): high
 $setGlobal cm_coal_scen  medCoal        !! def = medCoal  !! regexp = 0|lowCoal|medCoal|highCoal
+*** cm_XportRegiLim
+*** limit export levels on specific regions in relation to reference scenario levels
+*** example: "(2025, 2030, 2035, 2040, 2045, 2050).REF_regi.peoil 0.5" limits total exports of oil in REF to 50% of the values observed in the reference gdx
+$setGLobal cm_XportRegiLim off !! def off
+*** cm_XportRegiLimRampUp
+*** export limited levels defined at cm_XportRegiLim are gradually implemented with half of its level on the previous year
+*** example: "cm_XportRegiLimRampUp = on" and "cm_XportRegiLim = (2030, 2035, 2040, 2045, 2050).REF_regi.peoil 0.5", exports will be limited to 75% of the reference scenario by 2025 (cm_XportRegiLimRampUp), and 50% onwards
+$setGLobal cm_XportRegiLimRampUp off
+*** cm_altLearnRegiSet
+*** define an alternative group of regions that have independent learning cumulative capacity
+*** example: "(REF, CHA).(wind, spv)", REF and CHA will have their cumulative capacity for wind and solar pv learning calculations accounted idenpendently of the rest of the World.
+$setGLobal cm_altLearnRegiSet off !! def off
 *** c_ccsinjecrateRegi  "regional upper bound of the CCS injection rate, overwrites for specified regions the settings set with c_ccsinjecratescen"
 ***  ("off") no regional differentiation
 ***  ("GLO 0.005") reproduces c_ccsinjecratescen = 1

modules/22_subsidizeLearning/globallyOptimal/presolve.gms

@@ -9,9 +9,19 @@
 
 p22_deltacap0(ttot,regi,teLearn,rlf)$( (ttot.val ge 2005) and (pm_SolNonInfes(regi) eq 1)) = vm_deltaCap.l(ttot,regi,teLearn,rlf);
 
-loop((ttot,regi,teLearn),
-pm_capCumForeign(ttot,regi,teLearn)$(ttot.val ge 2005) = sum(regi2$(not sameas(regi2,regi)), pm_capCum0(ttot,regi2,teLearn) );  
-);
+$ifthen.altLearnRegiSet not "%cm_altLearnRegiSet%" == "off"
+  loop((ttot,regi,teLearn),
+    if(altLearnRegi22(regi,teLearn),
+      pm_capCumForeign(ttot,regi,teLearn)$(ttot.val ge 2005) = sum(regi2$((not sameas(regi2,regi)) and (altLearnRegi22(regi2,teLearn))), pm_capCum0(ttot,regi2,teLearn) );  
+    else
+      pm_capCumForeign(ttot,regi,teLearn)$(ttot.val ge 2005) = sum(regi2$((not sameas(regi2,regi)) and (not altLearnRegi22(regi2,teLearn))), pm_capCum0(ttot,regi2,teLearn) );  
+    );
+  );
+$else.altLearnRegiSet
+  loop((ttot,regi,teLearn),
+    pm_capCumForeign(ttot,regi,teLearn)$(ttot.val ge 2005) = sum(regi2$(not sameas(regi2,regi)), pm_capCum0(ttot,regi2,teLearn) );  
+  );
+$endif.altLearnRegiSet
 
 display pm_capCumForeign;
 

modules/22_subsidizeLearning/globallyOptimal/sets.gms

@@ -5,12 +5,16 @@
 *** |  REMIND License Exception, version 1.0 (see LICENSE file).
 *** |  Contact: [email protected]
 *** SOF ./modules/22_subsidizeLearning/globallyOptimal/sets.gms
-sets
-
-m22_learnteUpdateSub(all_te)  "learning technologies"
-
 
+sets
+  m22_learnteUpdateSub(all_te)  "learning technologies"
 ;
 
+$IFTHEN.altLearnRegiSet not "%cm_altLearnRegiSet%" == "off"
+Set
+  altLearnRegi22(all_regi,all_te) "Define regions to be accounted under an alternative learning cumulative capacity for each learning technology." 
+    / %cm_altLearnRegiSet% /
+;
+$ENDIF.altLearnRegiSet
 
 *** EOF ./modules/22_subsidizeLearning/globallyOptimal/sets.gms

modules/24_trade/standard/bounds.gms

@@ -128,4 +128,26 @@ if (cm_phaseoutBiolc eq 1,
 vm_Mport.fx(t,regi,entySe) = 0;
 vm_Xport.fx(t,regi,entySe) = 0;
 
+*** limit exports based on reference scenario results 
+$IFTHEN.XportRegiLim not "%cm_XportRegiLim%" == "off"
+Execute_Loadpoint 'input_ref' p24_Xport_refgdx = vm_Xport.l;
+  loop((ttot,ext_regi,enty)$(p24_XportRegiLim(ttot,ext_regi,enty) and (ttot.val ge cm_startyear)),
+    loop(regi$regi_groupExt(ext_regi,regi),
+      vm_Xport.up(ttot,regi,enty) = min(vm_Xport.up(ttot,regi,enty),p24_XportRegiLim(ttot,ext_regi,enty)*p24_Xport_refgdx(ttot,regi,enty));
+    );
+  );
+$IFTHEN.cm_XportRegiLimRampUp not "%cm_XportRegiLimRampUp%" == "off"
+loop((ext_regi,enty),
+  loop(ttot$(p24_XportRegiLim(ttot,ext_regi,enty) and (ttot.val ge cm_startyear)),
+    loop(ttot2,
+        break$((ttot2.val ge ttot.val) and (ttot2.val ge cm_startyear) and (ttot2.val ge 2020));
+        s24_prevYear = ttot2.val;
+    );
+    vm_Xport.up(ttot3,regi,enty)$(ttot3.val eq s24_prevYear) = min(vm_Xport.up(ttot3,regi,enty),(1-((1-p24_XportRegiLim(ttot,ext_regi,enty))/2))*p24_Xport_refgdx(ttot3,regi,enty));
+    break;
+  );
+);
+$ENDIF.cm_XportRegiLimRampUp
+$ENDIF.XportRegiLim
+
 *** EOF ./modules/24_trade/standard/bounds.gms

modules/24_trade/standard/declarations.gms

@@ -16,7 +16,16 @@ p24_Mport2005correct(all_regi,all_enty)     "Correction factor to match fossil s
 
 pm_MPortsPrice(tall,all_regi,tradeSe)              "Secondary energy import price for region. Calculated in the postsolve and assuming that trade is distributed uniformetly according existent capacities defined at p24_seTradeCapacity [T$/TWa]"
 pm_XPortsPrice(tall,all_regi,tradeSe)              "Secondary energy export price for region. Calculated in the postsolve and corresponding to the region secondary energy price [T$/TWa]"
+
+$ifthen.XportRegiLim not "%cm_XportRegiLim%" == "off"
+  p24_Xport_refgdx(tall,all_regi,all_enty)   "Export of traded commodity in the reference gdx"
+  p24_XportRegiLim(ttot,ext_regi,all_enty)   "limit of total exports per region or region group in relation to the reference scenario [%]" / %cm_XportRegiLim% /
+$IFTHEN.cm_XportRegiLimRampUp not "%cm_XportRegiLimRampUp%" == "off"
+  s24_prevYear                               "value of the previous year to the first exports limit"
+$ENDIF.cm_XportRegiLimRampUp
+$endIf.XportRegiLim
 ;
+
 ***-------------------------------------------------------------------------------
 ***                                   VARIABLES
 ***-------------------------------------------------------------------------------