Merge pull request remindmodel#1576 from 0UmfHxcvx5J7JoaOhFSs5mncnisT…

…JJ6q/fix/neg_fossil_ind_emi fix negative fossil emissions in industry
tabeado · Mar 1, 2024 · a41963e · a41963e
2 parents 0e2dddb + a00ef4d
commit a41963e
Show file tree

Hide file tree

Showing 2 changed files with 16 additions and 0 deletions.
diff --git a/modules/37_industry/subsectors/declarations.gms b/modules/37_industry/subsectors/declarations.gms
@@ -107,6 +107,7 @@ $endif.no_calibration
   q37_nonIncineratedPlastics(ttot,all_regi,all_enty,all_enty,all_emiMkt)    "calculate carbon contained in plastics that are not incinerated [GtC]"
   q37_feedstockEmiUnknownFate(ttot,all_regi,all_enty,all_enty,all_emiMkt)   "calculate carbon contained in chemical feedstock with unknown fate [GtC]"
   q37_feedstocksLimit(ttot,all_regi,all_enty,all_enty,all_emiMkt)           "restrict feedstocks flow to total energy flows into industry"
+  q37_feedstocksShares(ttot,all_regi,all_enty,all_enty,all_emiMkt)   "identical fossil/biomass/synfuel shares for FE and feedstocks"
 
   !! process-based implementation
   q37_demMatPrc(tall,all_regi,mat)                        "Material demand of processes"

diff --git a/modules/37_industry/subsectors/equations.gms b/modules/37_industry/subsectors/equations.gms
@@ -239,6 +239,21 @@ q37_feedstocksLimit(t,regi,entySe,entyFe,emiMkt)$(
   vm_demFENonEnergySector(t,regi,entySe,entyFe,"indst",emiMkt)
 ;
 
+*' Feedstocks have identical fossil/biomass/synfuel shares as industry FE
+q37_feedstocksShares(t,regi,entySE,entyFE,emiMkt)$(
+                                          sum(te, se2fe(entySE,entyFE,te)) ) ..
+    vm_demFEsector_afterTax(t,regi,entySE,entyFE,"indst",emiMkt)
+  * sum(se2fe(entySE2,entyFE,te),
+      vm_demFENonEnergySector(t,regi,entySE2,entyFE,"indst",emiMkt)
+    )
+  =e=
+    vm_demFENonEnergySector(t,regi,entySE,entyFE,"indst",emiMkt)
+  * sum(se2fe2(entySE2,entyFE,te),
+      vm_demFEsector_afterTax(t,regi,entySE2,entyFE,"indst",emiMkt)
+    )
+;
+
+
 *' Calculate mass of carbon contained in chemical feedstocks
 q37_FeedstocksCarbon(t,regi,sefe(entySe,entyFe),emiMkt)$(
                          entyFE2sector2emiMkt_NonEn(entyFe,"indst",emiMkt) ) ..