% =========================
% Adjusted for region alpha
% =========================

% load untitled.mat
% load untitled1.mat
% load untitled2.mat
% 
% u = find(time == 500);
% time = time(u:end)-time(u);
% 
% subplot(4,2,1)
% plot(time,Tci_results(u:end),'b')
% title('\itd\rm_1 [°C]')
% axis([0 time(end) 60 75])
% set(gca,'XTick',[0 2000 4000 6000])
% set(gca,'YTick',[60 67.5 75])
% 
% subplot(4,2,2)
% plot(time,wc_results(u:end),'b')
% title('\itd\rm_2 [t/h]')
% axis([0 time(end) 700 1000])
% set(gca,'XTick',[0 2000 4000 6000])
% set(gca,'YTick',[700 850 1000])
% 
% subplot(4,2,3)
% plot(time,Tci_F_results(u:end),'r')
% title('\ity\rm_1 [°C]')
% axis([0 time(end) 125.9 126.1])
% set(gca,'XTick',[0 2000 4000 6000])
% set(gca,'YTick',[125.9 126 126.1])
% 
% subplot(4,2,4)
% plot(time,Tco_results(u:end),'b')
% title('\ity\rm_5 [°C]')
% axis([0 time(end) 90 110])
% set(gca,'XTick',[0 2000 4000 6000])
% set(gca,'YTick',[90 100 110])
% 
% % subplot(3,2,5)
% % plot(time,uCI_results(2,u:end).*100,'g')
% % title('\itu\rm_4 [%]')
% % axis([0 time(end) 60 100])
% % set(gca,'XTick',[0 2000 4000 6000])
% % set(gca,'YTick',[60 80 100])
% 
% subplot(4,2,5)
% plot(time,uCI_results(2,u:end).*100,'g',time,uHB_results(2,u:end).*100,'k:')
% title('\itu\rm_4 [%] (solid); \itu\rm_1 [%] (dashed)')
% axis([0 time(end) 60 100])
% set(gca,'XTick',[0 2000 4000 6000])
% set(gca,'YTick',[60 80 100])
% 
% subplot(4,2,6)
% plot(time,wh_results(u:end),'r')
% title('\ity\rm_6 [t/h]')
% axis([0 time(end) 249 251])
% set(gca,'XTick',[0 2000 4000 6000])
% set(gca,'YTick',[249 250 251])
% 
% subplot(4,2,7)
% plot(time,P_results(2,u:end),'k')
% title('\itu\rm_6 [W]')
% axis([0 time(end) 4950 5080])
% set(gca,'XTick',[0 2000 4000 6000])
% set(gca,'YTick',[4950 5015 5080])

%========================
%Adjusted for region beta
%========================

% load untitled.mat
% load untitled1.mat
% load untitled2.mat
% 
% u = find(time == 500);
% time = time(u:end)-time(u);
% 
% subplot(4,2,1)
% plot(time,Tci_results(u:end),'b')
% title('\itd\rm_1 [°C]')
% axis([0 time(end) 60 80])
% set(gca,'XTick',[0 2000 4000 6000])
% set(gca,'YTick',[60 70 80])
% 
% subplot(4,2,2)
% plot(time,wc_results(u:end),'b')
% title('\itd\rm_2 [t/h]')
% axis([0 time(end) 400 1000])
% set(gca,'XTick',[0 2000 4000 6000])
% set(gca,'YTick',[400 700 1000])
% 
% subplot(4,2,3)
% plot(time,Tci_F_results(u:end),'r')
% title('\ity\rm_1 [°C]')
% axis([0 time(end) 125.9 126.1])
% set(gca,'XTick',[0 2000 4000 6000])
% set(gca,'YTick',[125.9 126 126.1])
% 
% subplot(4,2,4)
% plot(time,Tco_results(u:end),'b')
% title('\ity\rm_5 [°C]')
% axis([0 time(end) 100 130])
% set(gca,'XTick',[0 2000 4000 6000])
% set(gca,'YTick',[100 115 130])
% 
% % subplot(3,2,5)
% % plot(time,uCB_results(2,u:end).*100,'m')
% % title('\itu\rm_5 [%]')
% % axis([0 time(end) 0 100])
% % set(gca,'XTick',[0 2000 4000 6000])
% % set(gca,'YTick',[0 50 100])
% 
% subplot(4,2,5)
% plot(time,uCB_results(2,u:end).*100,'m',time,uHB_results(2,u:end).*100,'k:')
% title('\itu\rm_5 [%] (solid); \itu\rm_1 [%] (dashed)')
% axis([0 time(end) 0 100])
% set(gca,'XTick',[0 2000 4000 6000])
% set(gca,'YTick',[0 50 100])
% 
% subplot(4,2,6)
% plot(time,wh_results(u:end),'r')
% title('\ity\rm_6 [t/h]')
% axis([0 time(end) 249 251])
% set(gca,'XTick',[0 2000 4000 6000])
% set(gca,'YTick',[249 250 251])
% 
% subplot(4,2,7)
% plot(time,P_results(2,u:end),'k')
% title('\itu\rm_6 [W]')
% axis([0 time(end) 4860 5180])
% set(gca,'XTick',[0 2000 4000 6000])
% set(gca,'YTick',[4860 5020 5180])

% =========================
% Adjusted for region gamma
% =========================

% load untitled.mat
% load untitled1.mat
% 
% u = find(time == 500);
% time = time(u:end)-time(u);
% 
% subplot(4,2,1)
% plot(time,Tci_results(u:end),'b')
% title('\itd\rm_1 [°C]')
% axis([0 time(end) 60 90])
% set(gca,'XTick',[0 3500 7000])
% set(gca,'YTick',[60 75 90])
% 
% subplot(4,2,2)
% plot(time,wc_results(u:end),'b')
% title('\itd\rm_2 [t/h]')
% axis([0 time(end) 400 1000])
% set(gca,'XTick',[0 3500 7000])
% set(gca,'YTick',[400 700 1000])
% 
% subplot(4,2,3)
% plot(time,Tci_F_results(u:end),'r')
% title('\ity\rm_1 [°C]')
% axis([0 time(end) 125.9 126.1])
% set(gca,'XTick',[0 3500 7000])
% set(gca,'YTick',[125.9 126 126.1])
% 
% subplot(4,2,4)
% plot(time,Tco_results(u:end),'b')
% title('\ity\rm_5 [°C]')
% axis([0 time(end) 100 130])
% set(gca,'XTick',[0 3500 7000])
% set(gca,'YTick',[100 115 130])
% 
% subplot(4,2,5)
% plot(time,uHB_results(2,u:end).*100,'k:')
% title('\itu\rm_1 [%]')
% axis([0 time(end) 0 100])
% set(gca,'XTick',[0 3500 7000])
% set(gca,'YTick',[0 50 100])
% 
% subplot(4,2,6)
% plot(time,wh_results(u:end),'r')
% title('\ity\rm_6 [t/h]')
% axis([0 time(end) 245 255])
% set(gca,'XTick',[0 3500 7000])
% set(gca,'YTick',[245 250 255])
% 
% subplot(4,2,7)
% plot(time,P_results(2,u:end),'k')
% title('\itu\rm_6 [W]')
% axis([0 time(end) 5000 9000])
% set(gca,'XTick',[0 2000 4000 6000])
% set(gca,'YTick',[5000 7000 9000])

% =========================
% Adjusted for region delta
% =========================
% 
% load untitled.mat
% load untitled1.mat
% load untitled2.mat
% 
% % u = find(time == 500);
% % time = time(u:end)-time(u);
% 
% subplot(4,2,1)
% plot(time,Tci_results(u:end),'b')
% title('\itd\rm_1 [°C]')
% axis([0 time(end) 60 90])
% set(gca,'XTick',[0 3000 6000])
% set(gca,'YTick',[60 75 90])
% 
% subplot(4,2,2)
% plot(time,wc_results(u:end),'b')
% title('\itd\rm_2 [t/h]')
% axis([0 time(end) 400 1000])
% set(gca,'XTick',[0 3000 6000])
% set(gca,'YTick',[400 700 1000])
% 
% subplot(4,2,3)
% plot(time,Tci_F_results(u:end),'r')
% title('\ity\rm_1 [°C]')
% axis([0 time(end) 125.7 126.3])
% set(gca,'XTick',[0 3000 6000])
% set(gca,'YTick',[125.7 126 126.3])
% 
% subplot(4,2,4)
% plot(time,Tco_results(u:end),'b')
% title('\ity\rm_5 [°C]')
% axis([0 time(end) 115 135])
% set(gca,'XTick',[0 3000 6000])
% set(gca,'YTick',[115 125 135])
% 
% subplot(4,2,5)
% plot(time,uAC_results(2,u:end).*100,'c')
% title('\itu\rm_2 [%]')
% axis([0 time(end) 0 30])
% set(gca,'XTick',[0 3000 6000])
% set(gca,'YTick',[0 15 30])
% 
% subplot(4,2,6)
% plot(time,wh_results(u:end),'r')
% title('\ity\rm_6 [t/h]')
% axis([0 time(end) 245 255])
% set(gca,'XTick',[0 3000 6000])
% set(gca,'YTick',[245 250 255])
% 
% subplot(4,2,7)
% plot(time,P_results(2,u:end),'k')
% title('\itu\rm_6 [W]')
% axis([0 time(end) 5000 9000])
% set(gca,'XTick',[0 2000 4000 6000])
% set(gca,'YTick',[5000 7000 9000])
% 
% subplot(4,2,8)
% plot(time,Pf_results(2,u:end),'k')
% title('\itu\rm_7 [W]')
% axis([0 time(end) 0 30000])
% set(gca,'XTick',[0 3000 6000])
% set(gca,'YTick',[0 15000 30000])