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(1,:),'.',x,so4s(1,:),'-',x,so4sh(1,:),'-.')
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,'49\circN 16\circE')
title('SO^=_4','FontSize',15)
ylabel('Svratouch','FontSize',13)
h=subplot(4,3,4),plot(x,so4o(2,:),'.',x,so4s(2,:),'-',x,so4sh(2,:),'-.')
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 12\circE')
ylabel('Montelibretti','FontSize',13)
h=subplot(4,3,7),plot(x,so4o(3,:),'.',x,so4s(3,:),'-',x,so4sh(3,:),'-.')
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,'53\circN 6\circE')
ylabel('Kollumerwaard','FontSize',13)
h=subplot(4,3,10),plot(x,so4o(4,:),'.',x,so4s(4,:),'-',x,so4sh(4,:),'-.')
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,'69\circN 29\circE')
ylabel('Janiskoski','FontSize',13)

h=subplot(4,3,2),plot(x,no3o(1,:),'.',x,no3s(1,:),'-')
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(2,:),'.',x,no3s(2,:),'-')
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(3,:),'.',x,no3s(3,:),'-')
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(4,:),'.',x,no3s(4,:),'-')
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(1,:),'.',x,nh4s(1,:),'-')
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(2,:),'.',x,nh4s(2,:),'-')
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(3,:),'.',x,nh4s(3,:),'-')
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(4,:),'.',x,nh4s(4,:),'-')
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')
