delete empty list for flexible heat pumps

Author: CarlosEpia

edisgo/io/pypsa_csv_import.py

@@ -439,7 +439,9 @@ def _build_hv_trafos(raw: pd.DataFrame, base: pd.DataFrame) -> pd.DataFrame:
 
     # Active power - loads
     load_p = _read_ts(folder, "loads-p_set.csv", ts_slice)
-    #load_p = load_p * 10 # to force use of 14a
+    #load_p.loc[:,load_p.columns.str.contains("heat_")] = 0.1 # to force usage of 14a
+    #edisgo_obj.topology.loads_df.loc[edisgo_obj.topology.loads_df.index.str.contains("heat_"),"p_set"] = 0.1 # to force usage of 14a
+
     if not load_p.empty:
         ts.loads_active_power = _assign_timeseries(load_p, timeindex)
 

examples/loma-14a.py

@@ -434,7 +434,7 @@ def run_optimization_14a(edisgo):
     start_time = datetime.now()
 
     # Run optimization
-    edisgo.pm_optimize(opf_version=5, curtailment_14a=True, flexible_hps=[])
+    edisgo.pm_optimize(opf_version=5, curtailment_14a=True)
 
     duration = (datetime.now() - start_time).total_seconds()
 
@@ -446,7 +446,7 @@ def run_optimization_14a(edisgo):
 
 grid_path = "/home/carlos/LoMa/exec_folder/results/MGB_model_pypsa"
 
-edisgo = EDisGo(pypsa_csv_dir=grid_path, snapshot_range=(0, 3))
+edisgo = EDisGo(pypsa_csv_dir=grid_path, snapshot_range=(0, 2))
 mv_grid_geom = gpd.read_file(
     "/home/carlos/LoMa/exec_folder/data/Input_files/MV_grid_district/husum_district.shp"
 )
@@ -468,10 +468,12 @@ def run_optimization_14a(edisgo):
     edisgo.topology.buses_df.index != "HV_dummy_bus"
 ]
 
+######### This will be removed once Paul include the CP timeseries ############
 # Remove remporarily the charging points
 edisgo.topology.loads_df = edisgo.topology.loads_df[
     edisgo.topology.loads_df.type != "charging_point"
 ]
+######### This will be removed once Paul include the CP timeseries ############
 
 # Set Zero active power for batteries (let the OPF optimize dispatch freely)
 storage_names = edisgo.topology.storage_units_df.index
@@ -519,3 +521,13 @@ def run_optimization_14a(edisgo):
 for ts in edisgo.timeseries.timeindex:
     plot_network(edisgo, show=False, snapshot=str(ts))
 create_network_gif(duration=500)
+
+# edisgo.results.v_res
+# edisgo.results.s_res
+
+# load_t = edisgo.timeseries.loads_active_power
+# g = edisgo.topology.generators_df
+# gt = edisgo.timeseries.generators_active_power
+# lines = edisgo.topology.lines_df
+
+# demand_q = edisgo.timeseries.loads_reactive_power