
fil='../tmp/dic.nc';
egM=ncread(fil,'DIC_EG');                egM=egM(1:69);
deM=ncread(fil,'DIC_DEFAULT');           deM=deM(1:69);
coM=ncread(fil,'DIC_KISOKGM600');        coM=coM(1:69);
irM=ncread(fil,'DIC_DEFAULTIRONDOUBLE'); irM=irM(1:69);

%errplot_h([1:69],egM,egV.^0.5)

figure
 fs=12; ln=1.3; lnp=2;

subplot(3,1,1)
set(gca,'box','on','tickdir','out','nextplot','add','fontsize',fs,'linewidth',ln,'ygrid','on','xgrid','on') %,'ylim',[-0.1 0.2])
 plot(deM,'k','linewidth',lnp)
 plot(egM,'b','linewidth',lnp)
 plot(coM,'r','linewidth',lnp)
 plot(irM,'g','linewidth',lnp)
 
 legend('REF','E&G','CON','IRON','location','SouthWest')
 ylabel('DIC [mmol/m^3]')
% xlabel('year')
 mean(deM(1:20))
 mean(egM(1:20))
 mean(coM(1:20))

print -depsc ../ps/dic.eps
