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

ns=PAR.Nspecies;
it=round(.4*Nt);
figure
subplot(2,2,1)
plot(t,S,'-','LineWidth',LW)
title('Susceptable','FontSize',FS,'FontWeight',FW)
ylabel('number/m^2','FontSize',FS,'FontWeight',FW)
for s=1:ns
  text(t(it),S(it,s),num2str(s))
end

subplot(2,2,2)
plot(t,I,'-','LineWidth',LW)
title('Infected','FontSize',FS,'FontWeight',FW)
ylabel('number/m^2','FontSize',FS,'FontWeight',FW)
for s=1:ns
  text(t(it),I(it,s),num2str(s),'FontSize',FS,'FontWeight',FW)
end

subplot(2,2,3)
plot(t,D,'-','LineWidth',LW)
title('Dead Infected','FontSize',FS,'FontWeight',FW)
xlabel('day','FontSize',FS,'FontWeight',FW)
ylabel('number/m^2','FontSize',FS,'FontWeight',FW)
for s=1:ns
  text(t(it),D(it,s),num2str(s),'FontSize',FS,'FontWeight',FW)
end

subplot(2,2,4)
plot(t,P,'-','LineWidth',LW)
title('Infect Particles','FontSize',FS,'FontWeight',FW)
xlabel('day','FontSize',FS,'FontWeight',FW)
ylabel('number/m^3','FontSize',FS,'FontWeight',FW) 
for s=1:ns
  text(t(it),P(it,s),num2str(s),'FontSize',FS,'FontWeight',FW)
end
%print('-dpng',['Figures/Ab' PAR.case 'Num.png'])