Add plotting of hydrogen network including the Kernnetz (#71)

* created local script in Ariadne, added rule to Snakefile

* add kernnetz to the hydrogen_map plotting function

* changed logic to plot all maps at once

* revert functionality to plot network after each run

* add result to ariadne_all rule to trigger plotting

* [pre-commit.ci] auto fixes from pre-commit.com hooks

for more information, see https://pre-commit.ci

* plot only capacity which is build in the investment period not cumulative

* [pre-commit.ci] auto fixes from pre-commit.com hooks

for more information, see https://pre-commit.ci

* fix missing plots

* reset to file contents of main

* remove duplicate line

* improve numerical stability

* [pre-commit.ci] auto fixes from pre-commit.com hooks

for more information, see https://pre-commit.ci

---------

Co-authored-by: pre-commit-ci[bot] <66853113+pre-commit-ci[bot]@users.noreply.github.com>
Co-authored-by: Michael Lindner <michaellindner@posteo.de>

Author: 3 people

config/config.yaml

@@ -4,7 +4,7 @@
 
 # docs in https://pypsa-eur.readthedocs.io/en/latest/configuration.html#run
 run:
-  prefix: 20240926_meoh_load
+  prefix: 20240927plotH2Kernnetz
   name:
   # - CurrentPolicies
   - KN2045_Bal_v4

workflow/Snakefile

@@ -570,6 +570,13 @@ rule ariadne_all:
             RESULTS + "ariadne/capacity_detailed.png",
             run=config_provider("run", "name"),
         ),
+        expand(
+            RESULTS
+            + "maps/elec_s{simpl}_{clusters}_l{ll}_{opts}_{sector_opts}-h2_network_incl_kernnetz_{planning_horizons}.pdf",
+            run=config_provider("run", "name"),
+            **config["scenario"],
+            allow_missing=True,
+        ),
         exported_variables=expand(
             RESULTS + "ariadne/exported_variables_full.xlsx",
             run=config_provider("run", "name"),
@@ -602,3 +609,29 @@ rule check_sector_ratios:
         "logs/check_sector_ratios.log",
     script:
         "scripts/check_sector_ratios.py"
+
+
+rule plot_hydrogen_network_incl_kernnetz:
+    params:
+        plotting=config_provider("plotting"),
+        foresight=config_provider("foresight"),
+    input:
+        network=RESULTS
+        + "postnetworks/elec_s{simpl}_{clusters}_l{ll}_{opts}_{sector_opts}_{planning_horizons}.nc",
+        regions=resources("regions_onshore_elec_s{simpl}_{clusters}.geojson"),
+    output:
+        map=RESULTS
+        + "maps/elec_s{simpl}_{clusters}_l{ll}_{opts}_{sector_opts}-h2_network_incl_kernnetz_{planning_horizons}.pdf",
+    threads: 2
+    resources:
+        mem_mb=10000,
+    log:
+        RESULTS
+        + "logs/plot_hydrogen_network_incl_kernnetz/elec_s{simpl}_{clusters}_l{ll}_{opts}_{sector_opts}_{planning_horizons}.log",
+    benchmark:
+        (
+            RESULTS
+            + "benchmarks/plot_hydrogen_network_incl_kernnetz/elec_s{simpl}_{clusters}_l{ll}_{opts}_{sector_opts}_{planning_horizons}"
+        )
+    script:
+        "scripts/plot_hydrogen_network_incl_kernnetz.py"

workflow/scripts/export_ariadne_variables.py

@@ -150,6 +150,7 @@ def _get_fuel_fractions(n, region, fuel):
         ).sum()
         - renewable_fuel_supply.get(f"DE renewable {fuel}", pd.Series(0)).sum(),
         rtol=1e-3,
+        atol=1e-5,
     )
 
     fuel_fractions = pd.Series()
@@ -4464,7 +4465,7 @@ def get_data(
 
     if "debug" == "debug":  # For debugging
         var = pd.Series()
-        idx = 0
+        idx = -1
         n = networks[idx]
         c = costs[idx]
         _industry_demand = industry_demands[idx]