Compatibility with PyPSA-Ariadne derivative

config/config.default.yaml

@@ -454,6 +454,8 @@ sector:
   hydrogen_turbine: false
   SMR: true
   SMR_cc: true
+  regional_methanol_demand: false  #set to true if regional CO2 constraints needed
+  regional_oil_demand: false  #set to true if regional CO2 constraints needed
   regional_co2_sequestration_potential:
     enable: false
     attribute: 'conservative estimate Mt'
@@ -799,6 +801,7 @@ plotting:
     Coal: '#545454'
     coal: '#545454'
     Coal marginal: '#545454'
+    coal for industry: '#343434'
     solid: '#545454'
     Lignite: '#826837'
     lignite: '#826837'

rules/solve_myopic.smk

@@ -88,11 +88,13 @@ rule solve_sector_network_myopic:
         co2_sequestration_potential=config["sector"].get(
             "co2_sequestration_potential", 200
         ),
+        countries=config["countries"],
     input:
         network=RESULTS
         + "prenetworks-brownfield/elec_s{simpl}_{clusters}_l{ll}_{opts}_{sector_opts}_{planning_horizons}.nc",
         costs="data/costs_{planning_horizons}.csv",
         config=RESULTS + "config.yaml",
+        co2_totals_name=RESOURCES + "co2_totals.csv",
     output:
         RESULTS
         + "postnetworks/elec_s{simpl}_{clusters}_l{ll}_{opts}_{sector_opts}_{planning_horizons}.nc",

scripts/add_brownfield.py

@@ -120,6 +120,39 @@ def add_brownfield(n, n_p, year):
             n.links.loc[new_pipes, "p_nom_min"] = 0.0
 
 
+def disable_grid_expansion_if_LV_limit_hit(n):
+    if not "lv_limit" in n.global_constraints.index:
+        return
+
+    # calculate minimum LV
+    attr = "nom_min"
+    dc = n.links.index[n.links.carrier == "DC"]
+    tot = (n.lines["s_" + attr] * n.lines["length"]).sum() + (
+        n.links.loc[dc, "p_" + attr] * n.links.loc[dc, "length"]
+    ).sum()
+
+    diff = n.global_constraints.at["lv_limit", "constant"] - tot
+
+    # allow small numerical differences
+    limit = 1
+
+    if diff < limit:
+        logger.info(
+            f"LV is already reached (gap {diff}), disabling expansion and LV limit"
+        )
+        expandable_acs = n.lines.index[n.lines.s_nom_extendable]
+        n.lines.loc[expandable_acs, "s_nom_extendable"] = False
+        n.lines.loc[expandable_acs, "s_nom"] = n.lines.loc[expandable_acs, "s_nom_min"]
+
+        expandable_dcs = n.links.index[
+            n.links.p_nom_extendable & (n.links.carrier == "DC")
+        ]
+        n.links.loc[expandable_dcs, "p_nom_extendable"] = False
+        n.links.loc[expandable_dcs, "p_nom"] = n.links.loc[expandable_dcs, "p_nom_min"]
+
+        n.global_constraints.drop("lv_limit", inplace=True)
+
+
 if __name__ == "__main__":
     if "snakemake" not in globals():
         from _helpers import mock_snakemake
@@ -150,5 +183,7 @@ def add_brownfield(n, n_p, year):
 
     add_brownfield(n, n_p, year)
 
+    disable_grid_expansion_if_LV_limit_hit(n)
+
     n.meta = dict(snakemake.config, **dict(wildcards=dict(snakemake.wildcards)))
     n.export_to_netcdf(snakemake.output[0])

scripts/add_existing_baseyear.py

@@ -303,7 +303,7 @@ def add_power_capacities_installed_before_baseyear(n, grouping_years, costs, bas
         else:
             bus0 = vars(spatial)[carrier[generator]].nodes
             if "EU" not in vars(spatial)[carrier[generator]].locations:
-                bus0 = bus0.intersection(capacity.index + " gas")
+                bus0 = bus0.intersection(capacity.index + " " + carrier[generator])
 
             # check for missing bus
             missing_bus = pd.Index(bus0).difference(n.buses.index)