Add plotting script for FE KPIs

[virio-andreyana]

Closes #34

New Features

• Curtailment map
• Line loading map
• Energy balance time series
• Emission bar (detail and overview)
• System cost bar (detail and overview)
• Capacity bar (detail and overview), also added the option to chose install, optimal, expand
• Annual Generation bar (detail and overview)

The time series and the bar are created based on the config file, for example

  country_system_cost:
    extract: system cost
    exclude: ["EU"]
    carrier_filter:
    group_carrier: sector
    plot: overview
    plot_kw:
      title: Country system cost
      ylim: [0,150]
      ylabel: Billion EUR/a

automaticly creates this plot and save the file based on the name country_system_cost. The documentation is under construction.

Other example:

  DE_electricity_plus_winter_energy_balance:
    extract: time series
    start_date: '2013-02-01'
    end_date: '2013-03-01'
    include: ["DE"]
    carrier_filter: electricity+
    group_carrier: pretty
    plot_kw:
      title: Germany winter high and low voltage energy balance
      ylim: [-50,200]
      ylabel: MW

Checklist

• I tested my contribution locally and it works as intended.
• Code and workflow changes are sufficiently documented.
• Changes in configuration options are added in config/config.default.yaml.
• Changes in configuration options are documented in doc/configtables/*.csv.
• Sources of newly added data are documented in doc/data_sources.rst.
• A release note doc/release_notes.rst is added.

[virio-andreyana]

Currently I have a problem in creating a snakemake output so that the output can expand based on the KPI list. The exception is "Not all output, log and benchmark files of rule plot_KPIs contain the same wildcards. This is crucial though, in order to avoid that two or more jobs write to the same file." This is referring to unpack(create_kpi_path). Do you have any suggestion @daniel-rdt @tgilon ?

def create_kpi_path(w):
    return {
        kpi_fn: RESULTS
        + "maps/base_s_{clusters}_l{ll}_{opts}_{sector_opts}-"
        + kpi_fn
        + "_{planning_horizons}.pdf"
        for kpi_fn in list(config_provider("kpi").keys())
    }

rule plot_KPIs:
    params:
        plotting=config_provider("plotting"),
        kpi=config_provider("plotting"),
    input:
        network=RESULTS
        + "postnetworks/base_s_{clusters}_l{ll}_{opts}_{sector_opts}_{planning_horizons}.nc",
        regions=resources("regions_onshore_base_s_{clusters}.geojson"),
        nodal_costs=RESULTS + "csvs/nodal_costs.csv",
        nodal_capacity=RESULTS + "csvs/nodal_capacities.csv",
    output:
        unpack(create_kpi_path),
...

[virio-andreyana]

@daniel-rdt there is a simple workaround: The plots is created without going through the snakefile. I can extract the path from the still functioning curtailment_map or line_loading and put it there. The disadvantage is that you cannot request for it through snakemake, only infered request from curtailment_map or line_loading works. Do you want me to do that?

[daniel-rdt]

@daniel-rdt there is a simple workaround: The plots is created without going through the snakefile. I can extract the path from the still functioning curtailment_map or line_loading and put it there. The disadvantage is that you cannot request for it through snakemake, only infered request from curtailment_map or line_loading works. Do you want me to do that?

Hm, for the workflow to function that could be a valid plan B in my opinion. Let's have a quick chat about the snakemake issue and try to solve it.

[virio-andreyana]

TODO: This needs to be recalculated because it disregards the impact of the time resolution. Also don't forget to change the unit to GWh