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

fix bug carbon budget

Browse Source
This commit is contained in:
lisazeyen 2023-08-28 10:21:02 +02:00
parent b362c009f5
commit 369eaf3457
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
scripts/solve_network.py
View File

@ -247,12 +247,13 @@ def add_carbon_budget_constraint(n, snapshots):
n.stores["carrier"] = n.stores.bus.map(n.buses.carrier) n.stores["carrier"] = n.stores.bus.map(n.buses.carrier)
stores = n.stores.query("carrier in @emissions.index and not e_cyclic") stores = n.stores.query("carrier in @emissions.index and not e_cyclic")
time_valid = int(glc.loc["investment_period"]) time_valid = int(glc.loc["investment_period"])
weighting = n.investment_period_weightings.loc[time_valid, "years"]
if not stores.empty: if not stores.empty:
last = n.snapshot_weightings.reset_index().groupby("period").last() last = n.snapshot_weightings.reset_index().groupby("period").last()
last_i = last.set_index([last.index, last.timestep]).index last_i = last.set_index([last.index, last.timestep]).index
final_e = n.model["Store-e"].loc[last_i, stores.index] final_e = n.model["Store-e"].loc[last_i, stores.index]
time_i = pd.IndexSlice[time_valid, :] time_i = pd.IndexSlice[time_valid, :]
lhs = final_e.loc[time_i,:] lhs = final_e.loc[time_i,:] * weighting
n.model.add_constraints(lhs <= rhs, name=f"GlobalConstraint-{name}") n.model.add_constraints(lhs <= rhs, name=f"GlobalConstraint-{name}")
@ -738,6 +739,7 @@ def extra_functionality(n, snapshots):
add_pipe_retrofit_constraint(n) add_pipe_retrofit_constraint(n)
if n._multi_invest: if n._multi_invest:
add_carbon_constraint(n, snapshots) add_carbon_constraint(n, snapshots)
add_carbon_budget_constraint(n, snapshots)
def solve_network(n, config, solving, opts="", **kwargs): def solve_network(n, config, solving, opts="", **kwargs):