From be5331c89c8aeabd64108d1c15827da47519c43a Mon Sep 17 00:00:00 2001 From: Fabian Date: Fri, 22 Sep 2023 10:51:36 +0200 Subject: [PATCH] formulate sequestration limit constraint as operational_limit constraint --- scripts/prepare_sector_network.py | 2 +- scripts/solve_network.py | 5 +---- 2 files changed, 2 insertions(+), 5 deletions(-) diff --git a/scripts/prepare_sector_network.py b/scripts/prepare_sector_network.py index feb3faef..0caf45dd 100644 --- a/scripts/prepare_sector_network.py +++ b/scripts/prepare_sector_network.py @@ -587,7 +587,7 @@ def add_co2_tracking(n, costs, options): ) n.add("Carrier", "co2 stored") - # this tracks CO2 stored, e.g. underground + # this tracks CO2 sequestered, e.g. underground sequestration_buses = spatial.co2.nodes.str.replace(" stored", " sequestered") n.madd( "Bus", diff --git a/scripts/solve_network.py b/scripts/solve_network.py index 2bbd0164..e76d4004 100644 --- a/scripts/solve_network.py +++ b/scripts/solve_network.py @@ -182,9 +182,6 @@ def add_co2_sequestration_limit(n, config, limit=200): """ Add a global constraint on the amount of Mt CO2 that can be sequestered. """ - n.carriers.loc["co2 stored", "co2_absorptions"] = -1 - n.carriers.co2_absorptions = n.carriers.co2_absorptions.fillna(0) - limit = limit * 1e6 for o in opts: if "seq" not in o: @@ -396,7 +393,7 @@ def prepare_network( if snakemake.params["sector"]["limit_max_growth"]["enable"]: n = add_max_growth(n, config) - if n.stores.carrier.eq("co2 stored").any(): + if n.stores.carrier.eq("co2 sequestered").any(): limit = co2_sequestration_potential add_co2_sequestration_limit(n, config, limit=limit)