day_Amp_flows = np.zeros((24, numLines))
day_Temps = np.zeros((24, 1))
day_GHI = np.zeros((24, 1))
day_P_bus = np.zeros((24, 9))
day_Home_kW = np.zeros((24, numBuses))
day_EV_kW = np.zeros((24, numBuses))
day_PV_kW = np.zeros((24, numBuses))
day_Slack_kW_kVAR = np.zeros((24, 2))
for hr in range(24):
#dfSys['Bus'].Pd = dfSys['Bus'].Pd * (1+(hr/100))
# Apply KDE to Home Loads
from funcKDE import funcKDE
[loadHome_kW, loadHome_kVAR, hrTempC,
hrGHI] = funcKDE(dfHomeDay, dfTempDay, dfGHIDay, hr, numHomes,
numBuses)
day_Home_kW[hr, :] = loadHome_kW
day_Temps[hr, :] = hrTempC
day_GHI[hr, :] = hrGHI
# Apply Poisson to EV Loads
from funcPoiss import funcPoiss
[loadEV_kW, loadEV_kVAR,
EVatHome] = funcPoiss(dfEV, dfSys, maxEV, chgrRate, hr, EVstoHomes,
numHomes, numBuses)
day_EV_kW[hr, :] = loadEV_kW
# PV Solar Gen (no distribution)
from funcSolar import funcSolar
[genPV_kW, genPV_kVAR] = funcSolar(maxPV, hrGHI, hr, PVtoHomes,
numHomes, numBuses)
day_Amp_flows = np.zeros((24,numLines))
day_Temps = np.zeros((24,1))
day_GHI = np.zeros((24,1))
day_P_bus = np.zeros((24,9))
day_Home_kW = np.zeros((24,numBuses))
day_EV_kW = np.zeros((24,numBuses))
day_PV_kW = np.zeros((24,numBuses))
day_Slack_kW_kVAR = np.zeros((24,2))
for hr in range(24):
#dfSys['Bus'].Pd = dfSys['Bus'].Pd * (1+(hr/100))
# Apply KDE to Home Loads
from funcKDE import funcKDE
[loadHome_kW, loadHome_kVAR, hrTempC, hrGHI] = funcKDE(dfHomeDay, dfTempDay, dfGHIDay, hr, numHomes, numBuses)
day_Home_kW [hr,:] = loadHome_kW
day_Temps[hr,:] = hrTempC
day_GHI[hr,:] = hrGHI
# Apply Poisson to EV Loads
from funcPoiss import funcPoiss
[loadEV_kW, loadEV_kVAR, EVatHome] = funcPoiss(dfEV, dfSys, maxEV, chgrRate, hr, EVstoHomes, numHomes, numBuses)
day_EV_kW[hr,:] = loadEV_kW
# PV Solar Gen (no distribution)
from funcSolar import funcSolar
[genPV_kW, genPV_kVAR] = funcSolar(maxPV, hrGHI, hr, PVtoHomes, numHomes, numBuses)
day_PV_kW[hr,:] = genPV_kW
# Sum Home and EV Loads