% % ==================================
% % Adjusted for region alpha --> beta
% % ==================================
% 
% load untitled.mat
% load untitled1.mat
% load untitled2.mat
% load untitled3.mat
% 
% u = find(time == 500);
% time = time(u:end)-time(u);
% 
% subplot(5,2,1)
% plot(time,region_results(u:end),'r')
% title('Region')
% axis([0 time(end) 1 2])
% set(gca,'XTick',[0 750 1500])
% set(gca,'YTick',[1 2])
% 
% subplot(5,2,3)
% plot(time,Tci_results(u:end),'b')
% title('\itd\rm_1 [°C]')
% axis([0 time(end) 65 75])
% set(gca,'XTick',[0 750 1500])
% set(gca,'YTick',[65 70 75])
% 
% subplot(5,2,4)
% plot(time,840.8*ones(1,length(time)),'b')
% title('\itd\rm_2 [t/h]')
% axis([0 time(end) 400 1000])
% set(gca,'XTick',[0 750 1500])
% set(gca,'YTick',[400 700 1000])
% 
% subplot(5,2,5)
% 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 750 1500])
% set(gca,'YTick',[125.9 126 126.1])
% 
% subplot(5,2,6)
% plot(time,Tco_results(u:end),'b')
% title('\ity\rm_5 [°C]')
% axis([0 time(end) 90 110])
% set(gca,'XTick',[0 750 1500])
% set(gca,'YTick',[90 100 110])
% 
% subplot(5,2,7)
% plot(time,uCI_results(2,u:end).*100,'g',time,uHB_results(2,u:end).*100,'k:',time,uCB_results(2,u:end).*100,'m:')
% title('\itu\rm_1 [%] (long); \itu\rm_4 [%] (solid); \itu\rm_5 [%] (short)')
% axis([0 time(end) 30 100])
% set(gca,'XTick',[0 750 1500])
% set(gca,'YTick',[30 65 100])
% 
% subplot(5,2,8)
% plot(time,wh_results(u:end),'r')
% title('\ity\rm_6 [t/h]')
% axis([0 time(end) 249 251])
% set(gca,'XTick',[0 750 1500])
% set(gca,'YTick',[249 250 251])
% 
% subplot(5,2,9)
% plot(time,P_results(2,u:end),'k')
% title('\itu\rm_6 [W]')
% axis([0 time(end) 5000 5200])
% set(gca,'XTick',[0 750 1500])
% set(gca,'YTick',[5000 5100 5200])

% ==================================
% Adjusted for region beta --> alpha
% ==================================
% 
% load untitled.mat
% load untitled1.mat
% load untitled2.mat
% load untitled3.mat
% 
% u = find(time == 500);
% time = time(u:end)-time(u);
% 
% subplot(5,2,1)
% plot(time,region_results(u:end),'r')
% title('Region')
% axis([0 time(end) 1 2])
% set(gca,'XTick',[0 750 1500])
% set(gca,'YTick',[1 2])
% 
% subplot(5,2,3)
% plot(time,Tci_results(u:end),'b')
% title('\itd\rm_1 [°C]')
% axis([0 time(end) 65 75])
% set(gca,'XTick',[0 750 1500])
% set(gca,'YTick',[65 70 75])
% 
% subplot(5,2,4)
% plot(time,840.8*ones(1,length(time)),'b')
% title('\itd\rm_2 [t/h]')
% axis([0 time(end) 400 1000])
% set(gca,'XTick',[0 750 1500])
% set(gca,'YTick',[400 700 1000])
% 
% subplot(5,2,5)
% 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 750 1500])
% set(gca,'YTick',[125.9 126 126.1])
% 
% subplot(5,2,6)
% plot(time,Tco_results(u:end),'b')
% title('\ity\rm_5 [°C]')
% axis([0 time(end) 90 110])
% set(gca,'XTick',[0 750 1500])
% set(gca,'YTick',[90 100 110])
% 
% subplot(5,2,7)
% plot(time,uCI_results(2,u:end).*100,'g',time,uHB_results(2,u:end).*100,'k:',time,uCB_results(2,u:end).*100,'m:')
% title('\itu\rm_1 [%] (long); \itu\rm_4 [%] (solid); \itu\rm_5 [%] (short)')
% axis([0 time(end) 30 100])
% set(gca,'XTick',[0 750 1500])
% set(gca,'YTick',[30 65 100])
% 
% subplot(5,2,8)
% plot(time,wh_results(u:end),'r')
% title('\ity\rm_6 [t/h]')
% axis([0 time(end) 249 251])
% set(gca,'XTick',[0 750 1500])
% set(gca,'YTick',[249 250 251])
% 
% subplot(5,2,9)
% plot(time,P_results(2,u:end),'k')
% title('\itu\rm_6 [W]')
% axis([0 time(end) 5000 5200])
% set(gca,'XTick',[0 750 1500])
% set(gca,'YTick',[5000 5100 5200])

% ==================================
% Adjusted for region beta --> gamma
% ==================================

load untitled.mat
load untitled2.mat
load untitled3.mat

% u = find(time == 500);
% time = time(u:end)-time(u);

subplot(5,2,1)
plot(time,region_results(u:end),'r')
title('Region')
axis([0 time(end) 2 3])
set(gca,'XTick',[0 750 1500])
set(gca,'YTick',[2 3])

subplot(5,2,3)
plot(time,Tci_results(u:end),'b')
title('\itd\rm_1 [°C]')
axis([0 time(end) 65 85])
set(gca,'XTick',[0 750 1500])
set(gca,'YTick',[65 75 85])

subplot(5,2,4)
plot(time,840.8*ones(1,length(time)),'b')
title('\itd\rm_2 [t/h]')
axis([0 time(end) 400 1000])
set(gca,'XTick',[0 750 1500])
set(gca,'YTick',[400 700 1000])

subplot(5,2,5)
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 750 1500])
set(gca,'YTick',[125.9 126 126.1])

subplot(5,2,6)
plot(time,Tco_results(u:end),'b')
title('\ity\rm_5 [°C]')
axis([0 time(end) 100 130])
set(gca,'XTick',[0 750 1500])
set(gca,'YTick',[100 115 130])

subplot(5,2,7)
plot(time,uHB_results(2,u:end).*100,'k:',time,uCB_results(2,u:end).*100,'m:')
title('\itu\rm_1 [%] (long); \itu\rm_5 [%] (short)')
axis([0 time(end) 0 100])
set(gca,'XTick',[0 750 1500])
set(gca,'YTick',[0 50 100])

subplot(5,2,8)
plot(time,wh_results(u:end),'r')
title('\ity\rm_6 [t/h]')
axis([0 time(end) 248 266])
set(gca,'XTick',[0 750 1500])
set(gca,'YTick',[248 257 266])

subplot(5,2,9)
plot(time,P_results(2,u:end),'k')
title('\itu\rm_6 [W]')
axis([0 time(end) 4000 8000])
set(gca,'XTick',[0 750 1500])
set(gca,'YTick',[4000 6000 8000])