Circuit sharing and pausing for Panels project TEA

bin/postprocess.py

@@ -0,0 +1,65 @@
+# -*- coding: utf-8 -*-
+"""
+Created on Thu Jul 18 11:40:03 2024
+
+@author: jwang5
+"""
+
+import pandas as pd
+import os
+import numpy as np
+
+
+#%% identify unsimulated cases except timed out ones
+input_path=os.getcwd()
+
+df=pd.read_csv(os.path.join(input_path, 'simulation_log.csv'))
+df2=pd.read_csv(os.path.join(input_path, 'control_scenarios_processed.csv'))
+
+# order by upgrade and then building id
+df = df.sort_values(by=["Upgrade", "Building ID"], ascending=[True, True])
+
+# double check duplicates
+consecutive_duplicates = df['Building ID'] == df['Building ID'].shift()
+
+for i in range(len(consecutive_duplicates)):
+    if consecutive_duplicates.iloc[i]:
+        print(f'Found duplicate at index {i}')
+        # df = df.drop(i)
+
+# df.to_csv(os.path.join(input_path, 'simulation_log.csv'), index=False) 
+
+# generate a spreadsheet for failed simulations to rerun
+filtered_df = df[df["Status"].str.contains("Completed")]
+# check if ochre_complete exists in result folder
+for i in range(len(filtered_df)):
+    upgrade = filtered_df['Upgrade'].iloc[i]
+    bldg_id = filtered_df['Building ID'].iloc[i]
+    # result_path = os.path.join(os.getcwd(), '..', 'upgrade'+str(upgrade), 'results', 'simulation_output', 'up01',  'bldg'+str(bldg_id).zfill(7), 'run', 'ochre_output')
+    result_path = os.path.join(os.getcwd(), 'upgrade'+str(upgrade), str(bldg_id), 'ochre_output')
+    if not os.path.exists(os.path.join(result_path, 'OCHRE_complete')):
+        filtered_df['Status'].iloc[i] = 'Failed'
+
+filtered_df = filtered_df[filtered_df["Status"].str.contains("Failed")]
+
+filtered_df2 = df2[df2["building_id"].isin(filtered_df["Building ID"])]
+filtered_df2.to_csv(os.path.join(input_path, 'control_scenarios_cs_failed.csv'), index=False)
+
+
+#%% code for analyzing results
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+

bin/run_dwelling.py

@@ -23,8 +23,10 @@
     'time_zone': None,                            # option to specify daylight savings, in development
 
     # Input parameters - Sample building (uses HPXML file and time series schedule file)
-    'hpxml_file': os.path.join(default_input_path, 'Input Files', 'sample_resstock_properties.xml'),
-    'schedule_input_file': os.path.join(default_input_path, 'Input Files', 'sample_resstock_schedule.csv'),
+    #'hpxml_file': os.path.join(default_input_path, 'Input Files', 'sample_resstock_properties.xml'),
+    #'schedule_input_file': os.path.join(default_input_path, 'Input Files', 'sample_resstock_schedule.csv'),
+    'hpxml_file': 'C:/OCHRE-120V-hpwh/in.xml',
+    'schedule_input_file': 'C:/OCHRE-120V-hpwh/schedules.csv',
 
     # Input parameters - weather (note weather_path can be used when Weather Station is specified in HPXML file)
     # 'weather_path': weather_path,

bin/run_equipment.py

@@ -17,10 +17,10 @@
 default_args = {
     'start_time': dt.datetime(2018, 1, 1, 0, 0),  # year, month, day, hour, minute
     'time_res': dt.timedelta(minutes=15),
-    'duration': dt.timedelta(days=10),
-    'verbosity': 6,  # verbosity of results (1-9)
-    'save_results': False,  # if True, must specify output_path
-    # 'output_path': os.getcwd(),
+    'duration': dt.timedelta(days=365),
+    'verbosity': 7,  # verbosity of results (1-9)
+    'save_results': True,  # if True, must specify output_path
+    # 'output_path': os.path.join(os.getcwd(),'EVProfiles','Level2')
 }
 
 
@@ -185,25 +185,29 @@ def run_hvac():
     CreateFigures.plt.show()
 
 
-def run_ev():
+def run_ev(seed):
     equipment_args = {
         # Equipment parameters
-        'vehicle_type': 'BEV',
-        'charging_level': 'Level 1',
-        'mileage': 100,
+        "vehicle_type": "BEV",
+        "charging_level": "Level 1",
+        "equipment_event_file": "pdf_Veh4_Level2.csv",
+        "capacity": 57.5,
+        "seed": seed,
+        "output_path": os.path.join(os.getcwd(), "EVProfiles", "Level2", str(seed)),
         **default_args,
     }
 
     # Initialize equipment
     equipment = ElectricVehicle(**equipment_args)
 
     # Simulate equipment
+    equipment.main_simulator = True
     df = equipment.simulate()
 
-    print(df.head())
-    CreateFigures.plot_daily_profile(df, 'EV Electric Power (kW)', plot_max=False, plot_min=False)
-    CreateFigures.plot_time_series_detailed((df['EV SOC (-)'],))
-    CreateFigures.plt.show()
+    # print(df.head())
+    # CreateFigures.plot_daily_profile(df, 'EV Electric Power (kW)', plot_max=False, plot_min=False)
+    # CreateFigures.plot_time_series_detailed((df['EV SOC (-)'],))
+    # CreateFigures.plt.show()
 
 
 def run_equipment_from_house_model():
@@ -228,9 +232,10 @@ def run_equipment_from_house_model():
 if __name__ == '__main__':
     # Choose a scenario to run:
 
-    run_battery()
+    # run_battery()
     # run_battery_controlled()
     # run_water_heater()
     # run_hvac()
-    # run_ev()
+    for i in range(1,2):
+        run_ev(i)
     # run_equipment_from_house_model()