%  PLOTabalone1.m
%
%   analyze and plot the model results

LW=2;FS=14;FW='demi';
Nx=PAR.Nx;Ny=PAR.Ny;Ns=PAR.Nspecies;
it=round(.4*Nt);

figure
subplot(2,2,1)
plot(t,S(:,1,1,1),'-','LineWidth',LW)
hold on
for is=1:Ns
 for ix=1:Nx
  for iy=1:Ny
   plot(t,S(:,is,iy,ix),'-','LineWidth',LW)
   text(t(it),S(it,is,iy,ix),num2str(is))
  end
 end
end
hold off
title('Susceptable','FontSize',FS,'FontWeight',FW)
ylabel('number/m^2','FontSize',FS,'FontWeight',FW)

subplot(2,2,2)
plot(t,I(:,1,1,1),'-','LineWidth',LW)
hold on
for is=1:Ns
 for ix=1:Nx
  for iy=1:Ny
   plot(t,I(:,is,iy,ix),'-','LineWidth',LW)
   text(t(it),I(it,is,iy,ix),num2str(is))
  end
 end
end
hold off
title('Infected','FontSize',FS,'FontWeight',FW)
ylabel('number/m^2','FontSize',FS,'FontWeight',FW)

subplot(2,2,3)
plot(t,D(:,1,1,1),'-','LineWidth',LW)
hold on
for is=1:Ns
 for ix=1:Nx
  for iy=1:Ny
   plot(t,D(:,is,iy,ix),'-','LineWidth',LW)
   text(t(it),D(it,is,iy,ix),num2str(is))
  end
 end
end
hold off
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(:,1,1,1),'-','LineWidth',LW)
hold on
for is=1:Ns
 for ix=1:Nx
  for iy=1:Ny
   plot(t,P(:,is,iy,ix),'-','LineWidth',LW)
   text(t(it),P(it,is,iy,ix),num2str(is))
  end
 end
end
hold off
title('Infect Particles','FontSize',FS,'FontWeight',FW)
xlabel('day','FontSize',FS,'FontWeight',FW)
ylabel('number/m^3','FontSize',FS,'FontWeight',FW) 
print('-dpng',['Figures/Ab' PAR.case 'Num.png'])


%                 for Nx(column) by Ny(row) domain
figure
 for ir=1:Ny
  for ic=1:Nx
   subplot(Ny,Nx,(Ny-ir)*Nx+ic)
     plot(t,P(:,1,ir,ic),'LineWidth',LW,'LineWidth',LW)
     text(t(it),P(it,1,ir,ic),num2str(1))
     hold on 
     for is=2:Ns
       plot(t,P(:,is,ir,ic),'LineWidth',LW,'LineWidth',LW)
       text(t(it),P(it,is,ir,ic),num2str(is))
     end
     hold off
     datetick('x','mmmdd')
  end
 end
%   put labels
 for ic=1:Nx
  subplot(Ny,Nx,ic)
  title('Inf Particles','FontSize',FS,'FontWeight',FW)
 end
 for ic=1:Nx
  subplot(Ny,Nx,(Ny-1)*Nx+ic)
  xlabel('Days','FontSize',FS,'FontWeight',FW)
 end
% set(gca,'LineWidth',LW,'FontSize',FS,'FontWeight',FW)
 print('-dpng',['Figures/Ab' PAR.case 'Pgrid.png'])

%                 for Nx(column) by Ny(row) domain
figure
 for ir=1:Ny
  for ic=1:Nx
   subplot(Ny,Nx,(Ny-ir)*Nx+ic)
     plot(t,I(:,1,ir,ic),'LineWidth',LW)
     text(t(it),I(it,1,ir,ic),num2str(1))
     hold on
     for is=2:Ns
       plot(t,I(:,is,ir,ic),'LineWidth',LW)
       text(t(it),I(it,is,ir,ic),num2str(is))
     end
     hold off
     datetick('x','mmmdd')
  end
 end
%   put labels
 for ic=1:Nx
  subplot(Ny,Nx,ic)
  title('Infected','FontSize',FS,'FontWeight',FW)
 end
 for ic=1:Nx
  subplot(Ny,Nx,(Ny-1)*Nx+ic)
  xlabel('Days','FontSize',FS,'FontWeight',FW)
 end
% set(gca,'LineWidth',LW,'FontSize',FS,'FontWeight',FW)
 print('-dpng',['Figures/Ab' PAR.case 'Igrid.png'])

%                 for Nx(column) by Ny(row) domain
figure
 for ir=1:Ny
  for ic=1:Nx
   subplot(Ny,Nx,(Ny-ir)*Nx+ic)
     plot(t,S(:,1,ir,ic),'LineWidth',LW)
     text(t(it),S(it,1,ir,ic),num2str(1))
     hold on
     for is=2:Ns
       plot(t,S(:,is,ir,ic),'LineWidth',LW)
       text(t(it),S(it,is,ir,ic),num2str(is))
     end
     hold off
     datetick('x','mmmdd')
  end
 end
%   put labels
 for ic=1:Nx
  subplot(Ny,Nx,ic)
  title('Susceptables','FontSize',FS,'FontWeight',FW)
 end
 for ic=1:Nx
  subplot(Ny,Nx,(Ny-1)*Nx+ic)
  xlabel('Days','FontSize',FS,'FontWeight',FW)
 end
% set(gca,'LineWidth',LW,'FontSize',FS,'FontWeight',FW)
 print('-dpng',['Figures/Ab' PAR.case 'Sgrid.png'])