%  PLOTabalone1.m
%
%   analyze and plot the model results
LW=2;FS=14;FW='demi';

figure
 subplot(2,2,1)
  plot(t,S,'k-','LineWidth',LW)
  title('Susceptable','FontSize',FS,'FontWeight',FW)
  ylabel('number/m^2','FontSize',FS,'FontWeight',FW)
 subplot(2,2,2)
  plot(t,I,'k-','LineWidth',LW)
  title('Infected','FontSize',FS,'FontWeight',FW)
  ylabel('number/m^2','FontSize',FS,'FontWeight',FW)
 subplot(2,2,3)
  plot(t,D,'k-','LineWidth',LW)
  title('Dead Infected','FontSize',FS,'FontWeight',FW)
  xlabel('day','FontSize',FS,'FontWeight',FW)
  ylabel('number/m^2','FontSize',FS,'FontWeight',FW)
 subplot(2,2,4)
  plot(t,P,'k-','LineWidth',LW)
  title('Infect Particles','FontSize',FS,'FontWeight',FW)
  xlabel('day','FontSize',FS,'FontWeight',FW)
  ylabel('number/m^3','FontSize',FS,'FontWeight',FW)
 set(gca,'LineWidth',LW,'FontSize',FS,'FontWeight',FW)
% print('-dpng',['Figures/Ab' PAR.case 'Num.png'])

figure
 infP=PAR.IPinfect .* P .* S;
 infI=PAR.Iinfect  .*  I .* S;
 infD=PAR.Dinfect  .* D .* S;
 plot(t,infP,'k-',t,infI,'k--',t,infD,'k-.','LineWidth',LW)
 title('Infection rate (num/day): P(-) I(--) D(-.)'...
  ,'FontSize',FS,'FontWeight',FW)
 xlabel('day','FontSize',FS,'FontWeight',FW)
 ylabel('number/day','FontSize',FS,'FontWeight',FW)
 set(gca,'LineWidth',LW,'FontSize',FS,'FontWeight',FW)
% print('-dpng',['Figures/Ab' PAR.case 'InfRate.png'])