Implement sclopf

etrago/appl.py

@@ -28,6 +28,7 @@
 import datetime
 import os
 import os.path
+import pandas as pd
 
 __copyright__ = (
     "Flensburg University of Applied Sciences, "
@@ -51,7 +52,7 @@
     "db": "egon-data",  # database session
     "gridversion": None,  # None for model_draft or Version number
     "method": {  # Choose method and settings for optimization
-        "type": "market_grid",  # type of optimization, 'lopf' or 'market_grid'
+        "type": "market_grid",  # type of optimization, 'lopf', 'sclopf or 'market_grid'
         "n_iter": 1,  # abort criterion of iterative optimization, 'n_iter' or 'threshold'
         "pyomo": True,  # set if pyomo is used for model building
         "formulation": "pyomo",
@@ -736,6 +737,9 @@ def run_etrago(args, json_path):
     # skip snapshots
     etrago.skip_snapshots()
 
+    # Temporary drop DLR as it is currently not working with sclopf
+    etrago.network.lines_t.s_max_pu = pd.DataFrame(index=etrago.network.snapshots, columns=etrago.network.lines.index, data=1.)
+
     # start linear optimal powerflow calculations
 
     etrago.network.storage_units.cyclic_state_of_charge = True

etrago/execute/__init__.py

@@ -418,6 +418,15 @@ def optimize(self):
 
         self.grid_optimization()
 
+    elif self.args["method"]["type"] == "sclopf":
+        self.sclopf(
+            post_lopf=False,
+            extra_functionality=None,
+            n_process=4,
+            delta=0.01,
+            n_overload=0,
+            div_ext_lines=False,
+        )
     else:
         print("Method not defined")
 
@@ -579,7 +588,8 @@ def import_gen_from_links(network, drop_small_capacities=True):
     df["control"] = "PV"
     df.reset_index(inplace=True)
 
-    df.index = df.bus + " " + df.carrier
+    if not df.empty:
+        df.index = df.bus + " " + df.carrier
 
     # Aggregate disptach time series for new generators
     gas_to_add["bus1_carrier"] = gas_to_add.bus + " " + gas_to_add.carrier

etrago/execute/grid_optimization.py

@@ -51,7 +51,13 @@ def grid_optimization(self):
     #         self.network.links.index.str.contains('ramp')
     #         ].index, inplace=True)
     logger.info("Start solving grid optimization model")
-    self.lopf()
+    self.sclopf(
+        post_lopf=False,
+        n_process=4,
+        delta=0.01,
+        n_overload=0,
+        div_ext_lines=False,
+    )
 
 
 def fix_chp_generation(self):