fix regret ST plots

Author: JulianGeis

Snakefile

@@ -1197,23 +1197,27 @@ rule regret_plots_lt:
 
 rule regret_plots:
     params:
-        scenarios=get_scenarios(run),
+        scenarios=config["run"]["name"],
+        scenarios_config=get_scenarios(run),
         planning_horizons=config_provider("scenario", "planning_horizons"),
         plotting=config_provider("plotting"),
         output_dir=directory(
-            f"results/{config['run']['prefix']}/regret_plots/ST_comparison"
+            f"results/{config['run']['prefix']}/regret_plots/{{sensitivity}}/ST_comparison"
         ),
+
     input:
-        regret_networks=expand(
-            RESULTS
-            + "regret_networks/decision_{decision}_s_{clusters}_{opts}_{sector_opts}_{planning_horizons}.nc",
-            **config["scenario"],
-            allow_missing=True,
+        regret_networks=lambda wildcards: expand(
+            f"results/{config['run']['prefix']}/{{run}}/regret_networks/{wildcards.sensitivity}/decision_{{decision}}_s_{{clusters}}_{{opts}}_{{sector_opts}}_{{planning_horizons}}.nc",
             run=config["run"]["name"],
             decision=config["run"]["name"],
+            clusters=config["scenario"]["clusters"],
+            opts=config["scenario"]["opts"],
+            sector_opts=config["scenario"]["sector_opts"],
+            planning_horizons=config["scenario"]["planning_horizons"],
+            allow_missing=True,
         ),
     output:
-        elec_price_comp_de=f"results/{config['run']['prefix']}/regret_plots/ST_comparison/elec_price_comp_de.png",
+        elec_price_comp_de=f"results/{config['run']['prefix']}/regret_plots/{{sensitivity}}/ST_comparison/elec_price_comp_de.png",
     resources:
         mem_mb=32000,
     script:

config/config.de.yaml

@@ -4,7 +4,7 @@
 
 # docs in https://pypsa-eur.readthedocs.io/en/latest/configuration.html#run
 run:
-  prefix: 20250923_regret_uc_scaleCBF
+  prefix: 20250923_split_power_limit
   name:
   # - ExPol
   # - KN2045_Mix

scripts/pypsa-de/regret_plots.py

@@ -64,14 +64,16 @@ def aggregate_by_keywords(opex_comp_agg, groups):
 
         from _helpers import mock_snakemake
 
-        snakemake = mock_snakemake(
-            "regret_plots",
-        )
+    snakemake = mock_snakemake(
+        "regret_plots",
+        sensitivity="gas_price_60", 
+    )
 
     configure_logging(snakemake)
     config = snakemake.config
     planning_horizons = snakemake.params.planning_horizons
-    scenarios = ["HighDemand", "LowDemand"]
+    scenarios = snakemake.params.scenarios
+    decisions = ["decision_" + d for d in scenarios]
     tech_colors = snakemake.params.plotting["tech_colors"]
 
     # Nested dict: networks[year][scenario][decision] = Network
@@ -81,7 +83,7 @@ def aggregate_by_keywords(opex_comp_agg, groups):
         parts = fn.split(os.sep)
 
         # scenario is the folder name 2 levels up
-        scenario = parts[-3]
+        scenario = parts[-4]
         if scenario not in scenarios:
             raise ValueError(
                 f"Unexpected scenario '{scenario}' in {fn}. Allowed: {scenarios}"
@@ -111,12 +113,10 @@ def aggregate_by_keywords(opex_comp_agg, groups):
     # Plot electricity price duration curves
 
     fig, ax = plt.subplots(
-        figsize=(10, 5 * len(planning_horizons)), nrows=len(planning_horizons), ncols=1
+        figsize=(4*len(scenarios), 5 * len(planning_horizons)), nrows=len(planning_horizons), ncols=1
     )
     ax = ax.flatten()
 
-    decisions = ["decision_HighDemand", "decision_LowDemand"]
-
     for i, year in enumerate(planning_horizons):
         for scenario, decision in itertools.product(scenarios, decisions):
             n = networks[year][scenario][decision]
@@ -164,7 +164,7 @@ def aggregate_by_keywords(opex_comp_agg, groups):
     }
 
     fig, axes = plt.subplots(
-        nrows=len(planning_horizons), ncols=1, figsize=(12, 6 * len(planning_horizons))
+        nrows=len(planning_horizons), ncols=1, figsize=(6*len(scenarios), 6 * len(planning_horizons))
     )
     axes = axes.flatten()