load casnetdat
load emefsdat
load emepdat

x=0.001:100;
y=x;
y12=0.5*x;
y21=2*x;

so4o=casnetdat(:,5);
no3o=casnetdat(:,7);
nh4o=casnetdat(:,9);
so4s=casnetdat(:,6);
no3s=casnetdat(:,8);
nh4s=casnetdat(:,10);

f1=subplot(3,3,1),h=loglog(so4o,so4s,'.')
axis([0.001 10 0.001 10])
hold
%subplot(3,3,1),loglog(x,y,'LineWidth',2,x,y12,'.k',x,y21,'.k')
subplot(3,3,1),loglog(x,y,'LineWidth',2)
subplot(3,3,1),loglog(x,y12,'k',x,y21,'k')
set(f1,'XTick',[0.001 0.01 0.1 1 10])
set(f1,'XTickLabel','0.001|0.01|0.1|1|10')
set(f1,'YTick',[0.001 0.01 0.1 1 10])
set(f1,'YTickLabel','0.001|0.01|0.1|1|10')
set(h,'MarkerSize',6)
title('CASNET','FontSize',15)
ly=ylabel({'SO^=_4','model'},'FontSize',13)
set(ly,'Position',[0.000035 0.1 17])

subplot(3,3,4),h=loglog(no3o,no3s,'.')
axis([0.001 10 0.001 10])
hold
subplot(3,3,4),loglog(x,y,'LineWidth',2)
subplot(3,3,4),loglog(x,y12,'k',x,y21,'k')
set(gca,'XTick',[0.001 0.01 0.1 1 10])
set(gca,'XTickLabel','0.001|0.01|0.1|1|10')
set(gca,'YTick',[0.001 0.01 0.1 1 10])
set(gca,'YTickLabel','0.001|0.01|0.1|1|10')
set(h,'MarkerSize',6)
ly=ylabel({'NO^-_3','model'},'FontSize',13)
set(ly,'Position',[0.000035 0.1 17])

subplot(3,3,7),h=loglog(nh4o,nh4s,'.')
axis([0.001 10 0.001 10])
hold
subplot(3,3,7),loglog(x,y,'LineWidth',2)
subplot(3,3,7),loglog(x,y12,'k',x,y21,'k')
set(gca,'XTick',[0.001 0.01 0.1 1 10])
set(gca,'XTickLabel','0.001|0.01|0.1|1|10')
set(gca,'YTick',[0.001 0.01 0.1 1 10])
set(gca,'YTickLabel','0.001|0.01|0.1|1|10')
set(h,'MarkerSize',6)
ly=ylabel({'NH^+_4','model'},'FontSize',13)
set(ly,'Position',[0.000035 0.1 17])
xlabel('measure','FontSize',13)

!!!!!!!!!!!!!! EMEFS
so4o=emefsdat(:,5);
no3o=emefsdat(:,7);
nh4o=emefsdat(:,9);
so4s=emefsdat(:,6);
no3s=emefsdat(:,8);
nh4s=emefsdat(:,10);

subplot(3,3,2),h=loglog(so4o,so4s,'.')
axis([0.001 10 0.001 10])
hold
subplot(3,3,2),loglog(x,y,'LineWidth',2)
subplot(3,3,2),loglog(x,y12,'k',x,y21,'k')
set(gca,'XTick',[0.001 0.01 0.1 1 10])
set(gca,'XTickLabel','0.001|0.01|0.1|1|10')
set(gca,'YTick',[0.001 0.01 0.1 1 10])
set(gca,'YTickLabel','0.001|0.01|0.1|1|10')
set(h,'MarkerSize',6)
title('EMEFS','FontSize',15)

subplot(3,3,5),h=loglog(no3o,no3s,'.')
axis([0.001 10 0.001 10])
hold
subplot(3,3,5),loglog(x,y,'LineWidth',2)
subplot(3,3,5),loglog(x,y12,'k',x,y21,'k')
set(gca,'XTick',[0.001 0.01 0.1 1 10])
set(gca,'XTickLabel','0.001|0.01|0.1|1|10')
set(gca,'YTick',[0.001 0.01 0.1 1 10])
set(gca,'YTickLabel','0.001|0.01|0.1|1|10')
set(h,'MarkerSize',6)

subplot(3,3,8),h=loglog(nh4o,nh4s,'.')
axis([0.001 10 0.001 10])
hold
subplot(3,3,8),loglog(x,y,'LineWidth',2)
subplot(3,3,8),loglog(x,y12,'k',x,y21,'k')
set(gca,'XTick',[0.001 0.01 0.1 1 10])
set(gca,'XTickLabel','0.001|0.01|0.1|1|10')
set(gca,'YTick',[0.001 0.01 0.1 1 10])
set(gca,'YTickLabel','0.001|0.01|0.1|1|10')
set(h,'MarkerSize',6)
xlabel('measure','FontSize',13)

!!!!!!!!!!!!!! EMEP
so4o=emepdat(:,5);
no3o=emepdat(:,7);
nh4o=emepdat(:,9);
so4s=emepdat(:,6);
no3s=emepdat(:,8);
nh4s=emepdat(:,10);

subplot(3,3,3),h=loglog(so4o,so4s,'.')
axis([0.001 10 0.001 10])
hold
subplot(3,3,3),loglog(x,y,'LineWidth',2)
subplot(3,3,3),loglog(x,y12,'k',x,y21,'k')
set(gca,'XTick',[0.001 0.01 0.1 1 10])
set(gca,'XTickLabel','0.001|0.01|0.1|1|10')
set(gca,'YTick',[0.001 0.01 0.1 1 10])
set(gca,'YTickLabel','0.001|0.01|0.1|1|10')
set(h,'MarkerSize',6)
title('EMEP','FontSize',15)

subplot(3,3,6),h=loglog(no3o,no3s,'.')
axis([0.001 10 0.001 10])
hold
subplot(3,3,6),loglog(x,y,'LineWidth',2)
subplot(3,3,6),loglog(x,y12,'k',x,y21,'k')
set(gca,'XTick',[0.001 0.01 0.1 1 10])
set(gca,'XTickLabel','0.001|0.01|0.1|1|10')
set(gca,'YTick',[0.001 0.01 0.1 1 10])
set(gca,'YTickLabel','0.001|0.01|0.1|1|10')
set(h,'MarkerSize',6)

subplot(3,3,9),h=loglog(nh4o,nh4s,'.')
axis([0.001 10 0.001 10])
hold
subplot(3,3,9),loglog(x,y,'LineWidth',2)
subplot(3,3,9),loglog(x,y12,'k',x,y21,'k')
set(gca,'XTick',[0.001 0.01 0.1 1 10])
set(gca,'XTickLabel','0.001|0.01|0.1|1|10')
set(gca,'YTick',[0.001 0.01 0.1 1 10])
set(gca,'YTickLabel','0.001|0.01|0.1|1|10')
set(h,'MarkerSize',6)
xlabel('measure','FontSize',13)