
fil='../tmp/model_ice_extend.nc';
egM=ncread(fil,'COV_EG');           egM=egM(1:69);
deM=ncread(fil,'COV_DEFAULT');      deM=deM(1:69);
coM=ncread(fil,'COV_KISOKGM600');   coM=coM(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',[1.3 1.5]*1e7)
 plot(deM,'k','linewidth',lnp)
 plot(egM,'b','linewidth',lnp)
 plot(coM,'r','linewidth',lnp)
% xlabel('year')
 ylabel('ice extend [km^2]')
 mean(deM(1:20))
 mean(egM(1:20))
 mean(coM(1:20))
 print -depsc ../ps/ice.eps
