%  PLOTabalone3a.m
%
%   analyze and plot the model results
 LW=2;FS=14;FW='demi';
Nx=PAR.Nx;Ny=PAR.Ny;

figure
 subplot(2,2,1)
  plot(t,S(:,1,1),'LineWidth',LW)
  title('Susceptable','FontSize',FS,'FontWeight',FW)
  ylabel('number/m^2','FontSize',FS,'FontWeight',FW)
  hold on
  for ix=1:Nx
   for iy=1:Ny
    plot(t,S(:,iy,ix),'LineWidth',LW)
   end
  end 
  hold off
%  set(gca,'LineWidth',LW,'FontSize',FS,'FontWeight',FW)
 subplot(2,2,2)
  plot(t,I(:,1,1),'LineWidth',LW)
  title('Infected','FontSize',FS,'FontWeight',FW)
  ylabel('number/m^2','FontSize',FS,'FontWeight',FW)
  hold on 
  for ix=1:Nx
   for iy=1:Ny
    plot(t,I(:,iy,ix),'LineWidth',LW)
   end
  end 
 hold off
%  set(gca,'LineWidth',LW,'FontSize',FS,'FontWeight',FW)
 subplot(2,2,3)
  plot(t,D(:,1,1),'LineWidth',LW)
  title('Dead Infected','FontSize',FS,'FontWeight',FW)
  xlabel('day','FontSize',FS,'FontWeight',FW)
  ylabel('number/m^2','FontSize',FS,'FontWeight',FW)
  hold on 
  for ix=1:Nx
   for iy=1:Ny
    plot(t,D(:,iy,ix),'LineWidth',LW)
   end
  end 
  hold off
%  set(gca,'LineWidth',LW,'FontSize',FS,'FontWeight',FW)
  subplot(2,2,4)
   plot(t,P(:,1,1),'LineWidth',LW)
   title('Infect Particles','FontSize',FS,'FontWeight',FW)
   xlabel('day','FontSize',FS,'FontWeight',FW)
   ylabel('number/m^3','FontSize',FS,'FontWeight',FW)
   hold on 
   for ix=1:Nx
    for iy=1:Ny
     plot(t,P(:,iy,ix),'LineWidth',LW)
    end
  end 
  hold off
%  set(gca,'LineWidth',LW,'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)
      %         (3,1)=(r,c)
     plot(t,P(:,ir,ic),'LineWidth',LW)
     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(:,ir,ic),'LineWidth',LW)
     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)
      %         (3,1)=(r,c)
     plot(t,S(:,ir,ic),'LineWidth',LW)
     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'])