%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 
dis=dis*4;

%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_dis=dis(:,1);
oliv10_con_dis=dis(:,2);
oliv100_con_dis=dis(:,3);
oliv1000_con_dis=dis(:,4);
olivInf_con_t_dis=dis(:,5);
oliv10_con_t_dis=dis(:,6);
oliv100_con_t_dis=dis(:,7);
oliv1000_con_t_dis=dis(:,8);
olivInf_dis30x2E4_dis=dis(:,9);
oliv10_dis30x2E4_dis=dis(:,10);
oliv100_dis30x2E4_dis=dis(:,11);
oliv1000_dis30x2E4_dis=dis(:,12);
olivInf_dis9_dis=dis(:,13);
oliv10_dis9_dis=dis(:,14);
oliv100_dis9_dis=dis(:,15);
oliv1000_dis9_dis=dis(:,16);
control_dis=dis(:,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_dis,'b--','linewidth',2);
h2=plot(t(1:4:end),oliv1000_con_dis(1:4:end),'b-x','linewidth',2,'markersize',10);
h3=plot(t,oliv100_con_t_dis,'g--','linewidth',2);
h4=plot(t(1:4:end),oliv100_con_dis(1:4:end),'g-x','linewidth',2,'markersize',10);
h5=plot(t,oliv10_con_t_dis,'r--','linewidth',2);
h6=plot(t(1:4:end),oliv10_con_dis(1:4:end),'r-x','linewidth',2,'markersize',10);
h7=plot(t,olivInf_con_t_dis,'--','color',[0.5 0.5 0.5],'linewidth',2);
h8=plot(t(1:4:end),olivInf_con_dis(1:4:end),'-x','color',[0.5 0.5 0.5],'linewidth',2,'markersize',10);
h9=plot(t(1:4:end),control_dis(1:4:end),'k-x','linewidth',2,'markersize',10);
set(gca,'Yticklabel',{' '})

set(gca,'fontsize',24)
grid on
set(figure(1),'position',[100,100,450,600]);
xlim=[2020 2100];


set(gcf,'PaperPositionMode','auto')
set(gcf,'color','w')


print('group1_dis','-depsc')

clf

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

h5=plot(t,olivInf_dis9_dis,'-','color',[0.5 0.5 0.5],'linewidth',1.5);
h5=plot(t(3:4:end),olivInf_dis9_dis(3:4:end),'o','color',[0.5 0.5 0.5],'linewidth',1.5,'markersize',8);
h6=plot(t,oliv10_dis9_dis,'-','color','r','linewidth',1.5);
h6=plot(t(1:4:end),oliv10_dis9_dis(1:4:end),'o','color','r','linewidth',1.5,'markersize',8);
h7=plot(t,oliv100_dis9_dis,'-','color','g','linewidth',1.5);
h7=plot(t(3:4:end),oliv100_dis9_dis(3:4:end),'o','color','g','linewidth',1.5,'markersize',8);
h8=plot(t,oliv1000_dis9_dis,'-','color','b','linewidth',1.5);
h8=plot(t(1:4:end),oliv1000_dis9_dis(1:4:end),'o','color','b','linewidth',1.5,'markersize',8);
h9=plot(t,control_dis,'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',[0 2.5])


set(gcf,'PaperPositionMode','auto')
set(gcf,'color','w')
print('group2_dis','-depsc')
clf






