bjoern/pypsa-eur
1
0
Fork 0
Code Issues Pull Requests Actions 1 Packages Projects Releases Wiki Activity

build_carbon_budget: parse input_co2

Browse Source
This commit is contained in:
Fabian Neumann 2023-08-22 14:49:44 +02:00
parent 401fbb7be7
commit 3ccc63aae4
3 changed files with 10 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
scripts/make_summary.py
View File

@ -711,5 +711,5 @@ if __name__ == "__main__":
if snakemake.params.foresight == "myopic":
cumulative_cost = calculate_cumulative_cost()
cumulative_cost.to_csv(
"results/" + snakemake.params.RDIR + "/csvs/cumulative_cost.csv"
"results/" + snakemake.params.RDIR + "csvs/cumulative_cost.csv"
)

2
scripts/plot_summary.py
View File

@ -549,7 +549,7 @@ def plot_carbon_budget_distribution(input_eurostat):
fancybox=True, fontsize=18, loc=(0.01, 0.01), facecolor="white", frameon=True
)
path_cb_plot = "results/" + snakemake.params.RDIR + "/graphs/"
path_cb_plot = "results/" + snakemake.params.RDIR + "graphs/"
plt.savefig(path_cb_plot + "carbon_budget_plot.pdf", dpi=300)

11
scripts/prepare_sector_network.py
View File

@ -220,7 +220,7 @@ def co2_emissions_year(
# TODO: move to own rule with sector-opts wildcard?
def build_carbon_budget(o, input_eurostat, fn, emissions_scope, report_year):
def build_carbon_budget(o, input_eurostat, fn, emissions_scope, report_year, input_co2):
"""
Distribute carbon budget following beta or exponential transition path.
"""
@ -3402,13 +3402,18 @@ if __name__ == "__main__":
if "cb" not in o:
continue
limit_type = "carbon budget"
fn = "results/" + snakemake.params.RDIR + "/csvs/carbon_budget_distribution.csv"
fn = "results/" + snakemake.params.RDIR + "csvs/carbon_budget_distribution.csv"
if not os.path.exists(fn):
emissions_scope = snakemake.params.emissions_scope
report_year = snakemake.params.eurostat_report_year
input_co2 = snakemake.input.co2
build_carbon_budget(
o, snakemake.input.eurostat, fn, emissions_scope, report_year
o,
snakemake.input.eurostat,
fn,
emissions_scope,
report_year,
input_co2,
)
co2_cap = pd.read_csv(fn, index_col=0).squeeze()
limit = co2_cap.loc[investment_year]