load tsso4sim19dat
load tsno3sim19dat
load tsnh4sim19dat
load tsso4simhalodat

for n=1:8
m1=2*n-1
m2=2*n
so4s(n,:)=tsso4sim19dat(m1,5:16)
so4o(n,:)=tsso4sim19dat(m2,5:16);
no3s(n,:)=tsno3sim19dat(m1,5:16);
no3o(n,:)=tsno3sim19dat(m2,5:16);
nh4s(n,:)=tsnh4sim19dat(m1,5:16);
nh4o(n,:)=tsnh4sim19dat(m2,5:16);
so4sh(n,:)=tsso4simhalodat(m1,5:16);
end

x=1:12;
h=subplot(4,3,1),plot(x,so4o(5,:),'.',x,so4s(5,:),'-',x,so4sh(5,:),'-.')
axis([0 13 0 5])
set(h,'XTick',[1 2 3 4 5 6 7 8 9 10 11 12])
set(h,'XTickLabel','J|F|M|A|M|J|J|A|S|O|N|D')
text(1,4,'42\circN 76\circW')
title('SO^=_4','FontSize',15)
ylabel('Ithaca','FontSize',13)
h=subplot(4,3,4),plot(x,so4o(6,:),'.',x,so4s(6,:),'-',x,so4sh(6,:),'-.')
axis([0 13 0 5])
set(h,'XTick',[1 2 3 4 5 6 7 8 9 10 11 12])
set(h,'XTickLabel','J|F|M|A|M|J|J|A|S|O|N|D')
text(1,4,'43\circN 75\circW')
ylabel('Piseco','FontSize',13)
h=subplot(4,3,7),plot(x,so4o(7,:),'.',x,so4s(7,:),'-',x,so4sh(7,:),'-.')
axis([0 13 0 5])
set(h,'XTick',[1 2 3 4 5 6 7 8 9 10 11 12])
set(h,'XTickLabel','J|F|M|A|M|J|J|A|S|O|N|D')
text(1,4,'45\circN 85\circW')
ylabel('Gaylord','FontSize',13)
h=subplot(4,3,10),plot(x,so4o(8,:),'.',x,so4s(8,:),'-',x,so4sh(8,:),'-.')
axis([0 13 0 5])
set(h,'XTick',[1 2 3 4 5 6 7 8 9 10 11 12])
set(h,'XTickLabel','J|F|M|A|M|J|J|A|S|O|N|D')
text(1,4,'45\circN 79\circW')
ylabel('Dorset','FontSize',13)

h=subplot(4,3,2),plot(x,no3o(5,:),'.',x,no3s(5,:),'-')
axis([0 13 0 5])
set(h,'XTick',[1 2 3 4 5 6 7 8 9 10 11 12])
set(h,'XTickLabel','J|F|M|A|M|J|J|A|S|O|N|D')
title('NO^-_3','FontSize',15)
h=subplot(4,3,5),plot(x,no3o(6,:),'.',x,no3s(6,:),'-')
axis([0 13 0 5])
set(h,'XTick',[1 2 3 4 5 6 7 8 9 10 11 12])
set(h,'XTickLabel','J|F|M|A|M|J|J|A|S|O|N|D')
h=subplot(4,3,8),plot(x,no3o(7,:),'.',x,no3s(7,:),'-')
axis([0 13 0 5])
set(h,'XTick',[1 2 3 4 5 6 7 8 9 10 11 12])
set(h,'XTickLabel','J|F|M|A|M|J|J|A|S|O|N|D')
h=subplot(4,3,11),plot(x,no3o(8,:),'.',x,no3s(8,:),'-')
axis([0 13 0 5])
set(h,'XTick',[1 2 3 4 5 6 7 8 9 10 11 12])
set(h,'XTickLabel','J|F|M|A|M|J|J|A|S|O|N|D')

h=subplot(4,3,3),plot(x,nh4o(5,:),'.',x,nh4s(5,:),'-')
axis([0 13 0 5])
set(h,'XTick',[1 2 3 4 5 6 7 8 9 10 11 12])
set(h,'XTickLabel','J|F|M|A|M|J|J|A|S|O|N|D')
title('NH^+_4','FontSize',15)
h=subplot(4,3,6),plot(x,nh4o(6,:),'.',x,nh4s(6,:),'-')
axis([0 13 0 5])
set(h,'XTick',[1 2 3 4 5 6 7 8 9 10 11 12])
set(h,'XTickLabel','J|F|M|A|M|J|J|A|S|O|N|D')
h=subplot(4,3,9),plot(x,nh4o(7,:),'.',x,nh4s(7,:),'-')
axis([0 13 0 5])
set(h,'XTick',[1 2 3 4 5 6 7 8 9 10 11 12])
set(h,'XTickLabel','J|F|M|A|M|J|J|A|S|O|N|D')
h=subplot(4,3,12),plot(x,nh4o(8,:),'.',x,nh4s(8,:),'-')
axis([0 13 0 5])
set(h,'XTick',[1 2 3 4 5 6 7 8 9 10 11 12])
set(h,'XTickLabel','J|F|M|A|M|J|J|A|S|O|N|D')
