From d58a7f86a45c4825c843d1a35f64715f7cd7c376 Mon Sep 17 00:00:00 2001
From: Fabian Neumann <fabian.neumann@outlook.de>
Date: Fri, 9 Jul 2021 13:22:00 +0200
Subject: [PATCH] allow copperplated carbon capture to be distributed freely to
 co2 stores

---
 scripts/prepare_sector_network.py | 21 ++++++++++++---------
 1 file changed, 12 insertions(+), 9 deletions(-)

diff --git a/scripts/prepare_sector_network.py b/scripts/prepare_sector_network.py
index aced4719..16eb12f6 100644
--- a/scripts/prepare_sector_network.py
+++ b/scripts/prepare_sector_network.py
@@ -1748,12 +1748,13 @@ def add_industry(n, costs):
         efficiency=1.
     )
 
-    n.add("Link",
-        "solid biomass for industry CC",
+    n.madd("Link",
+        spatial.co2.locations,
+        suffix=" solid biomass for industry CC",
         bus0="EU solid biomass",
         bus1="solid biomass for industry",
         bus2="co2 atmosphere",
-        bus3="co2 stored", # TODO co2: where to allocate if co2 is spatially resolved?
+        bus3=spatial.co2.nodes,
         carrier="solid biomass for industry CC",
         p_nom_extendable=True,
         capital_cost=costs.at["cement capture", "fixed"] * costs.at['solid biomass', 'CO2 intensity'],
@@ -1786,12 +1787,13 @@ def add_industry(n, costs):
         efficiency2=costs.at['gas', 'CO2 intensity']
     )
 
-    n.add("Link",
-        "gas for industry CC",
+    n.madd("Link",
+        spatial.co2.locations,
+        suffix=" gas for industry CC",
         bus0="EU gas",
         bus1="gas for industry",
         bus2="co2 atmosphere",
-        bus3="co2 stored",  # TODO co2: where to allocate if co2 is spatially resolved?
+        bus3=spatial.co2.nodes,
         carrier="gas for industry CC",
         p_nom_extendable=True,
         capital_cost=costs.at["cement capture", "fixed"] * costs.at['gas', 'CO2 intensity'],
@@ -1963,11 +1965,12 @@ def add_industry(n, costs):
     )
 
     #assume enough local waste heat for CC
-    n.add("Link",
-        "process emissions CC",
+    n.madd("Link",
+        spatial.co2.locations,
+        suffix=" process emissions CC",
         bus0="process emissions",
         bus1="co2 atmosphere",
-        bus2="co2 stored",  # TODO co2: where to allocate if co2 is spatially resolved?
+        bus2=spatial.co2.nodes,
         carrier="process emissions CC",
         p_nom_extendable=True,
         capital_cost=costs.at["cement capture", "fixed"],