make CC in blast furnaces optional

config/config.default.yaml

@@ -446,6 +446,7 @@ sector:
   methanation: true
   helmeth: false
   coal_cc: false
+  coal_for_industry_cc: false
   dac: true
   co2_vent: false
   allam_cycle: false

doc/configtables/sector.csv

@@ -73,6 +73,7 @@ marginal_cost_storage,currency/MWh ,float,The marginal cost of discharging batte
 methanation,--,"{true, false}",Add option for transforming hydrogen and CO2 into methane using methanation.
 helmeth,--,"{true, false}",Add option for transforming power into gas using HELMETH (Integrated High-Temperature ELectrolysis and METHanation for Effective Power to Gas Conversion)
 coal_cc,--,"{true, false}",Add option for coal CHPs with carbon capture
+coal_for_industry_cc,--,"{true, false}",Add option for carbon capture in integrated steel making (blast furnaces)
 dac,--,"{true, false}",Add option for Direct Air Capture (DAC)
 co2_vent,--,"{true, false}",Add option for vent out CO2 from storages to the atmosphere.
 allam_cycle,--,"{true, false}",Add option to include `Allam cycle gas power plants <https://en.wikipedia.org/wiki/Allam_power_cycle>`_

scripts/prepare_sector_network.py

@@ -2940,25 +2940,27 @@ def add_industry(n, costs):
             efficiency2=costs.at["coal", "CO2 intensity"],
         )
 
-        n.madd(
+        if options.get("coal_for_industry_cc", False):
-            "Link",
+
-            spatial.coal.industry_cc,
+            n.madd(
-            bus0=spatial.coal.nodes,
+                "Link",
-            bus1=spatial.coal.industry,
+                spatial.coal.industry_cc,
-            bus2="co2 atmosphere",
+                bus0=spatial.coal.nodes,
-            bus3=spatial.co2.nodes,
+                bus1=spatial.coal.industry,
-            carrier="coal for industry CC",
+                bus2="co2 atmosphere",
-            p_min_pu=1.0,
+                bus3=spatial.co2.nodes,
-            p_nom_extendable=True,
+                carrier="coal for industry CC",
-            capital_cost=costs.at["cement capture", "fixed"]
+                p_min_pu=1.0,
-            * costs.at["coal", "CO2 intensity"],
+                p_nom_extendable=True,
-            efficiency=0.9,
+                capital_cost=costs.at["cement capture", "fixed"]
-            efficiency2=costs.at["coal", "CO2 intensity"]
+                * costs.at["coal", "CO2 intensity"],
-            * (1 - costs.at["cement capture", "capture_rate"]),
+                efficiency=0.9,
-            efficiency3=costs.at["coal", "CO2 intensity"]
+                efficiency2=costs.at["coal", "CO2 intensity"]
-            * costs.at["cement capture", "capture_rate"],
+                * (1 - costs.at["cement capture", "capture_rate"]),
-            lifetime=costs.at["cement capture", "lifetime"],
+                efficiency3=costs.at["coal", "CO2 intensity"]
-        )
+                * costs.at["cement capture", "capture_rate"],
+                lifetime=costs.at["cement capture", "lifetime"],
+            )
 
 
 def add_waste_heat(n):