20230705IPC to Staff PR 4 Attachment - LOLP_Calc.pdf

% Programmed by Andres Valdepena Delgado 2022
% Email ‐ [e-mail]

% Input Description:
% Hourly net load
% Outage table
% Output Description:
% LOLPs
% LOLE
% LOLH
[LOLP_Jan_hr, states_jan] = Fun_LOLP(Outages,Capacities_Jan,Load_jan);
[LOLP_Feb_hr, states_feb] = Fun_LOLP(Outages,Capacities_Feb,Load_feb);
[LOLP_Mar_hr, states_mar] = Fun_LOLP(Outages,Capacities_Mar,Load_mar);
[LOLP_Apr_hr, states_apr] = Fun_LOLP(Outages,Capacities_Apr,Load_apr);
[LOLP_May_hr, states_may] = Fun_LOLP(Outages,Capacities_May,Load_may);
[LOLP_Jun_hr, states_jun] = Fun_LOLP(Outages,Capacities_Jun,Load_jun);
[LOLP_Jul_hr, states_jul] = Fun_LOLP(Outages,Capacities_Jul,Load_jul);
[LOLP_Aug_hr, states_aug] = Fun_LOLP(Outages,Capacities_Aug,Load_aug);
[LOLP_Sep_hr, states_sep] = Fun_LOLP(Outages,Capacities_Sep,Load_sep);
[LOLP_Oct_hr, states_oct] = Fun_LOLP(Outages,Capacities_Oct,Load_oct);
[LOLP_Nov_hr, states_nov] = Fun_LOLP(Outages,Capacities_Nov,Load_nov);
[LOLP_Dec_hr, states_dec] = Fun_LOLP(Outages,Capacities_Dec,Load_dec);
Jan_daily = reshape(LOLP_Jan_hr, [24, 31]);
Feb_daily = reshape(LOLP_Feb_hr, [24, 28]);
Mar_daily = reshape(LOLP_Mar_hr, [24, 31]);
Apr_daily = reshape(LOLP_Apr_hr, [24, 30]);
May_daily = reshape(LOLP_May_hr, [24, 31]);
Jun_daily = reshape(LOLP_Jun_hr, [24, 30]);
Jul_daily = reshape(LOLP_Jul_hr, [24, 31]);
Aug_daily = reshape(LOLP_Aug_hr, [24, 31]);
Sep_daily = reshape(LOLP_Sep_hr, [24, 30]);
Oct_daily = reshape(LOLP_Oct_hr, [24, 31]);
Nov_daily = reshape(LOLP_Nov_hr, [24, 30]);
Dec_daily = reshape(LOLP_Dec_hr, [24, 31]);
LOLP_daily = [Jan_daily Feb_daily Mar_daily Apr_daily May_daily Jun_daily Jul_daily
Aug_daily Sep_daily Oct_daily Nov_daily Dec_daily];
LOLH_Jan = sum(LOLP_Jan_hr);
LOLH_Feb = sum(LOLP_Feb_hr);
LOLH_Mar = sum(LOLP_Mar_hr);
LOLH_Apr = sum(LOLP_Apr_hr);
LOLH_May = sum(LOLP_May_hr);
LOLH_Jun = sum(LOLP_Jun_hr);
LOLH_Jul = sum(LOLP_Jul_hr);
LOLH_Aug = sum(LOLP_Aug_hr);
LOLH_Sep = sum(LOLP_Sep_hr);
LOLH_Oct = sum(LOLP_Oct_hr);
LOLH_Nov = sum(LOLP_Nov_hr);
LOLH_Dec = sum(LOLP_Dec_hr);
LOLE_Jan = sum(max(Jan_daily));
LOLE_Feb = sum(max(Feb_daily));
LOLE_Mar = sum(max(Mar_daily));
LOLE_Apr = sum(max(Apr_daily));
LOLE_May = sum(max(May_daily));
LOLE_Jun = sum(max(Jun_daily));
LOLE_Jul = sum(max(Jul_daily));
LOLE_Aug = sum(max(Aug_daily));
LOLE_Sep = sum(max(Sep_daily));
LOLE_Oct = sum(max(Oct_daily));
LOLE_Nov = sum(max(Nov_daily));
LOLE_Dec = sum(max(Dec_daily));
LOLE_Monthly = [LOLE_Jan LOLE_Feb LOLE_Mar LOLE_Apr LOLE_May LOLE_Jun LOLE_Jul 
LOLE_Aug LOLE_Sep LOLE_Oct LOLE_Nov LOLE_Dec];
LOLH_Monthly = [LOLH_Jan LOLH_Feb LOLH_Mar LOLH_Apr LOLH_May LOLH_Jun LOLH_Jul 
LOLH_Aug LOLH_Sep LOLH_Oct LOLH_Nov LOLH_Dec];
LOLE_d = sum(LOLE_Monthly);
LOLH_d = sum(LOLH_Monthly);
LOLE_h = LOLE_d
LOLH_h = LOLH_d;
LOLP_hourly = [LOLP_Jan_hr; LOLP_Feb_hr; LOLP_Mar_hr; LOLP_Apr_hr; LOLP_May_hr; 
LOLP_Jun_hr; LOLP_Jul_hr; LOLP_Aug_hr; LOLP_Sep_hr; LOLP_Oct_hr; LOLP_Nov_hr; 
LOLP_Dec_hr];