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

figure
 subplot(2,2,1)
  plot(t,S,'k-',t,S2,'k--','LineWidth',LW)
  title('Susceptable (-) Remote (--)','FontSize',FS,'FontWeight',FW)
  ylabel('number/m^2','FontSize',FS,'FontWeight',FW)
 subplot(2,2,2)
  plot(t,I,'k-','LineWidth',LW)
  title(['Infected  Case ' PAR.case],...
        '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'])

infP=PAR.IPinfect .* P .* S;
infI=PAR.Iinfect  .* I .* S;
infD=PAR.Dinfect  .* D .* S;
figure
 plot(t,infP,'k-',t,infI,'k--',t,infD,'k-.','LineWidth',LW)
 title(['Infection rate (num/day): P(-) I(--) D(-.) ' PAR.case],...
       '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'])

K=(1 - (S+I)/PAR.Carry);i=find(K<0);K(i)=0;
repS=PAR.ReproS .* K .* S;
repI=PAR.ReproI .* K .* I;
K2=(1 - (S2)/PAR.Carry2);i=find(K2<0);K2(i)=0;
rep2=PAR.Repro2 .* K .* S2;
figure
 plot(t,repS,'k-',t,repI,'k--',t,rep2,'k-.','LineWidth',LW)
 title(['Reproduction rate (num/day): S(-) I(--), S2(-.) ' PAR.case] ,...
       '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 'RepRate.png'])

D=S2-S;if D < 0; D=0; end
inpI=PAR.Imm * S2; 
inpD=PAR.Diff * D;
figure
 plot(t,inpI,'k-',t,inpD,'k--','LineWidth',LW)
 title(['Immigration rate (num/day): Im(-) Dif(--) ' PAR.case],...
       '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 'ImmRate.png'])