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

fix lossy bidirectional link coupling countraint for myopic

Browse Source
This commit is contained in:
Fabian Neumann 2024-01-03 12:46:42 +01:00
parent 29afffb4ca
commit 05495ce484
1 changed files with 14 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
scripts/solve_network.py
View File

@ -691,13 +691,24 @@ def add_lossy_bidirectional_link_constraints(n):
if not n.links.p_nom_extendable.any() or not "reversed" in n.links.columns: if not n.links.p_nom_extendable.any() or not "reversed" in n.links.columns:
return return
reversed_links = n.links.reversed.fillna(0).astype(bool) n.links["reversed"] = n.links.reversed.fillna(0).astype(bool)
carriers = n.links.loc[reversed_links, "carrier"].unique() carriers = n.links.loc[n.links.reversed, "carrier"].unique()
forward_i = n.links.query( forward_i = n.links.query(
"carrier in @carriers and ~reversed and p_nom_extendable" "carrier in @carriers and ~reversed and p_nom_extendable"
).index ).index
backward_i = forward_i + "-reversed"
def get_backward_i(forward_i):
return pd.Index(
[
re.sub(r"-(\d{4})$", r"-reversed-\1", s)
if re.search(r"-\d{4}$", s)
else s + "-reversed"
for s in forward_i
]
)
backward_i = get_backward_i(forward_i)
lhs = n.model["Link-p_nom"].loc[backward_i] lhs = n.model["Link-p_nom"].loc[backward_i]
rhs = n.model["Link-p_nom"].loc[forward_i] rhs = n.model["Link-p_nom"].loc[forward_i]