From 5d6d01f8c02f3f49e98ad47b279ce8b33e037c47 Mon Sep 17 00:00:00 2001
From: Fabian Neumann <fabian.neumann@outlook.de>
Date: Thu, 23 Nov 2023 12:11:12 +0100
Subject: [PATCH] make CC in blast furnaces optional

---
 config/config.default.yaml        |  1 +
 doc/configtables/sector.csv       |  1 +
 scripts/prepare_sector_network.py | 40 ++++++++++++++++---------------
 3 files changed, 23 insertions(+), 19 deletions(-)

diff --git a/config/config.default.yaml b/config/config.default.yaml
index 7bae66d1..e0eb4c4c 100644
--- a/config/config.default.yaml
+++ b/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
diff --git a/doc/configtables/sector.csv b/doc/configtables/sector.csv
index d610c862..2ddf8584 100644
--- a/doc/configtables/sector.csv
+++ b/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>`_
diff --git a/scripts/prepare_sector_network.py b/scripts/prepare_sector_network.py
index eb6fb221..99daba83 100644
--- a/scripts/prepare_sector_network.py
+++ b/scripts/prepare_sector_network.py
@@ -2940,25 +2940,27 @@ def add_industry(n, costs):
             efficiency2=costs.at["coal", "CO2 intensity"],
         )
 
-        n.madd(
-            "Link",
-            spatial.coal.industry_cc,
-            bus0=spatial.coal.nodes,
-            bus1=spatial.coal.industry,
-            bus2="co2 atmosphere",
-            bus3=spatial.co2.nodes,
-            carrier="coal for industry CC",
-            p_min_pu=1.0,
-            p_nom_extendable=True,
-            capital_cost=costs.at["cement capture", "fixed"]
-            * costs.at["coal", "CO2 intensity"],
-            efficiency=0.9,
-            efficiency2=costs.at["coal", "CO2 intensity"]
-            * (1 - costs.at["cement capture", "capture_rate"]),
-            efficiency3=costs.at["coal", "CO2 intensity"]
-            * costs.at["cement capture", "capture_rate"],
-            lifetime=costs.at["cement capture", "lifetime"],
-        )
+        if options.get("coal_for_industry_cc", False):
+
+            n.madd(
+                "Link",
+                spatial.coal.industry_cc,
+                bus0=spatial.coal.nodes,
+                bus1=spatial.coal.industry,
+                bus2="co2 atmosphere",
+                bus3=spatial.co2.nodes,
+                carrier="coal for industry CC",
+                p_min_pu=1.0,
+                p_nom_extendable=True,
+                capital_cost=costs.at["cement capture", "fixed"]
+                * costs.at["coal", "CO2 intensity"],
+                efficiency=0.9,
+                efficiency2=costs.at["coal", "CO2 intensity"]
+                * (1 - costs.at["cement capture", "capture_rate"]),
+                efficiency3=costs.at["coal", "CO2 intensity"]
+                * costs.at["cement capture", "capture_rate"],
+                lifetime=costs.at["cement capture", "lifetime"],
+            )
 
 
 def add_waste_heat(n):