From 76d57731ec12042d990e46dcdb0a972c0c5ef3f2 Mon Sep 17 00:00:00 2001 From: Fabian Neumann Date: Mon, 6 Mar 2023 17:18:07 +0100 Subject: [PATCH] fix all lv references to ll, pickup right location of default lifetime --- scripts/add_existing_baseyear.py | 4 ++-- scripts/make_summary.py | 14 +++++++------- scripts/plot_network.py | 6 +++--- 3 files changed, 12 insertions(+), 12 deletions(-) diff --git a/scripts/add_existing_baseyear.py b/scripts/add_existing_baseyear.py index bd6183e5..60974362 100644 --- a/scripts/add_existing_baseyear.py +++ b/scripts/add_existing_baseyear.py @@ -160,7 +160,7 @@ def add_power_capacities_installed_before_baseyear(n, grouping_years, costs, bas mean = df_agg.loc[biomass_i, "DateIn"].mean() df_agg.loc[biomass_i, "DateIn"] = df_agg.loc[biomass_i, "DateIn"].fillna(int(mean)) # Fill missing DateOut - dateout = df_agg.loc[biomass_i, "DateIn"] + snakemake.config["costs"]["lifetime"] + dateout = df_agg.loc[biomass_i, "DateIn"] + snakemake.config["costs"]['fill_values']["lifetime"] df_agg.loc[biomass_i, "DateOut"] = df_agg.loc[biomass_i, "DateOut"].fillna(dateout) # drop assets which are already phased out / decommissioned @@ -654,7 +654,7 @@ if __name__ == "__main__": .to_pandas() .reindex(index=n.snapshots) ) - default_lifetime = snakemake.config["costs"]["lifetime"] + default_lifetime = snakemake.config["costs"]['fill_values']["lifetime"] add_heating_capacities_installed_before_baseyear( n, baseyear, diff --git a/scripts/make_summary.py b/scripts/make_summary.py index 8322dafc..858881cc 100644 --- a/scripts/make_summary.py +++ b/scripts/make_summary.py @@ -208,15 +208,15 @@ def calculate_cumulative_cost(): # integrate cost throughout the transition path for r in cumulative_cost.columns: for cluster in cumulative_cost.index.get_level_values(level=0).unique(): - for lv in cumulative_cost.index.get_level_values(level=1).unique(): + for ll in cumulative_cost.index.get_level_values(level=1).unique(): for sector_opts in cumulative_cost.index.get_level_values( level=2 ).unique(): cumulative_cost.loc[ - (cluster, lv, sector_opts, "cumulative cost"), r + (cluster, ll, sector_opts, "cumulative cost"), r ] = np.trapz( cumulative_cost.loc[ - idx[cluster, lv, sector_opts, planning_horizons], r + idx[cluster, ll, sector_opts, planning_horizons], r ].values, x=planning_horizons, ) @@ -639,7 +639,7 @@ def make_summaries(networks_dict): ] columns = pd.MultiIndex.from_tuples( - networks_dict.keys(), names=["cluster", "lv", "opt", "planning_horizon"] + networks_dict.keys(), names=["cluster", "ll", "opt", "planning_horizon"] ) df = {} @@ -676,14 +676,14 @@ if __name__ == "__main__": logging.basicConfig(level=snakemake.config["logging"]["level"]) networks_dict = { - (cluster, lv, opt + sector_opt, planning_horizon): "results/" + (cluster, ll, opt + sector_opt, planning_horizon): "results/" + snakemake.params.RDIR - + f"/postnetworks/elec_s{simpl}_{cluster}_lv{lv}_{opt}_{sector_opt}_{planning_horizon}.nc" + + f"/postnetworks/elec_s{simpl}_{cluster}_l{ll}_{opt}_{sector_opt}_{planning_horizon}.nc" for simpl in snakemake.config["scenario"]["simpl"] for cluster in snakemake.config["scenario"]["clusters"] for opt in snakemake.config["scenario"]["opts"] for sector_opt in snakemake.config["scenario"]["sector_opts"] - for lv in snakemake.config["scenario"]["lv"] + for ll in snakemake.config["scenario"]["ll"] for planning_horizon in snakemake.config["scenario"]["planning_horizons"] } diff --git a/scripts/plot_network.py b/scripts/plot_network.py index 33a92b75..0b327a56 100644 --- a/scripts/plot_network.py +++ b/scripts/plot_network.py @@ -136,7 +136,7 @@ def plot_map( ac_color = "rosybrown" dc_color = "darkseagreen" - if snakemake.wildcards["lv"] == "1.0": + if snakemake.wildcards["ll"] == "v1.0": # should be zero line_widths = n.lines.s_nom_opt - n.lines.s_nom link_widths = n.links.p_nom_opt - n.links.p_nom @@ -698,7 +698,7 @@ def plot_map_without(network): to_drop = n.links.index[(n.links.carrier != "DC") & (n.links.carrier != "B2B")] n.links.drop(to_drop, inplace=True) - if snakemake.wildcards["lv"] == "1.0": + if snakemake.wildcards["ll"] == "v1.0": line_widths = n.lines.s_nom link_widths = n.links.p_nom else: @@ -895,7 +895,7 @@ def plot_series(network, carrier="AC", name="test"): "{}/{RDIR}maps/series-{}-{}-{}-{}-{}.pdf".format( "results", snakemake.params.RDIR, - snakemake.wildcards["lv"], + snakemake.wildcards["ll"], carrier, start, stop,