%Better do it on laptop because the dashed curve will look better on that display
% This goes to generate the time series 
% First I need to copy the time_series.exl into matlab

% load 

olivInf_con_ph=ph(:,1);
oliv10_con_ph=ph(:,2);
oliv100_con_ph=ph(:,3);
oliv1000_con_ph=ph(:,4);
olivInf_con_t_ph=ph(:,5);
oliv10_con_t_ph=ph(:,6);
oliv100_con_t_ph=ph(:,7);
oliv1000_con_t_ph=ph(:,8);
olivInf_ph30x2E4_ph=ph(:,9);
oliv10_ph30x2E4_ph=ph(:,10);
oliv100_ph30x2E4_ph=ph(:,11);
oliv1000_ph30x2E4_ph=ph(:,12);
olivInf_ph9_ph=ph(:,13);
oliv10_ph9_ph=ph(:,14);
oliv100_ph9_ph=ph(:,15);
oliv1000_ph9_ph=ph(:,16);
control_ph=ph(:,17);

% colors: 
% black
% '0.5 0.5 0.5','c','m','r','g','b' and its dashed
% [0.6 0.65 0.1],[0.5 0.75 0.5][0.7 0.2 0.2][0.5 0.75 1]

% figure 1

t=2020:2099

figure(1)
hold on
h1=plot(t,oliv1000_con_t_ph,'b--','linewidth',2);
h2=plot(t,oliv1000_con_ph,'b','linewidth',5);
h3=plot(t,oliv100_con_t_ph,'g--','linewidth',2);
h4=plot(t,oliv100_con_ph,'g','linewidth',5);
h5=plot(t,oliv10_con_t_ph,'r--','linewidth',2);
h6=plot(t,oliv10_con_ph,'r','linewidth',5);
h7=plot(t,olivInf_con_t_ph,'--','color',[0.5 0.5 0.5],'linewidth',2);
h8=plot(t,olivInf_con_ph,'color',[0.5 0.5 0.5],'linewidth',5);
h9=plot(t,control_ph,'k','linewidth',2);
set(gca,'Yticklabel',{' '})

set(gca,'fontsize',24)
grid on
set(figure(1),'position',[100,100,450,600]);
xlim=[2020 2100];
set(gca,'ylim',[7.7 8.8])
set(gcf,'PaperPositionMode','auto')
set(gcf,'color','w')

print('group1_glob_ph','-depsc')

clf

hold on
h1=plot(t,olivInf_ph30x2E4_ph,'-','color',[0.5 0.5 0.5],'linewidth',1.5);
h1=plot(t(1:4:end),olivInf_ph30x2E4_ph(1:4:end),'^','color',[0.5 0.5 0.5],'linewidth',1.5,'markersize',8);
h2=plot(t,oliv10_ph30x2E4_ph,'-','color','r','linewidth',1.5);
h2=plot(t(3:4:end),oliv10_ph30x2E4_ph(3:4:end),'^','color','r','linewidth',1.5,'markersize',8);
h3=plot(t,oliv100_ph30x2E4_ph,'-','color','g','linewidth',1.5);
h3=plot(t(1:4:end),oliv100_ph30x2E4_ph(1:4:end),'^','color','g','linewidth',1.5,'markersize',8);
h4=plot(t,oliv1000_ph30x2E4_ph,'-','color','b','linewidth',1.5);
h4=plot(t(3:4:end),oliv1000_ph30x2E4_ph(3:4:end),'^','color','b','linewidth',1.5,'markersize',8);

h5=plot(t,olivInf_ph9_ph,'-','color',[0.5 0.5 0.5],'linewidth',1.5);
h5=plot(t(3:4:end),olivInf_ph9_ph(3:4:end),'o','color',[0.5 0.5 0.5],'linewidth',1.5,'markersize',8);
h6=plot(t,oliv10_ph9_ph,'-','color','r','linewidth',1.5);
h6=plot(t(1:4:end),oliv10_ph9_ph(1:4:end),'o','color','r','linewidth',1.5,'markersize',8);
h7=plot(t,oliv100_ph9_ph,'-','color','g','linewidth',1.5);
h7=plot(t(3:4:end),oliv100_ph9_ph(3:4:end),'o','color','g','linewidth',1.5,'markersize',8);
h8=plot(t,oliv1000_ph9_ph,'-','color','b','linewidth',1.5);
h8=plot(t(1:4:end),oliv1000_ph9_ph(1:4:end),'o','color','b','linewidth',1.5,'markersize',8);
h9=plot(t,control_ph,'k','linewidth',2);
set(gca,'Yticklabel',{' '})




set(gca,'fontsize',24)
grid on
set(figure(1),'position',[100,100,450,600]);
xlim=[2020 2100];
set(gca,'ylim',[7.7 8.8])
set(gcf,'PaperPositionMode','auto')
set(gcf,'color','w')
print('group2_glob_ph','-depsc')
clf


