
fil='../tmp/heat.nc';
egM=ncread(fil,'NHEATM_EG');   egM=egM(1:69);
egV=ncread(fil,'NHEATSTD_EG'); egV=egV(1:69);
deM=ncread(fil,'NHEATM_DEFAULT');   deM=deM(1:69);
deV=ncread(fil,'NHEATSTD_DEFAULT'); deV=deV(1:69);
coM=ncread(fil,'NHEATM_KISOKGM600');   coM=coM(1:69);
coV=ncread(fil,'NHEATSTD_KISOKGM600'); coV=coV(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','nextplot','add','fontsize',fs,'linewidth',ln,'ygrid','on','xgrid','on','ylim',[-10 -3], ...
         'tickdir','out')
 plot(deM,'k','linewidth',lnp)
 plot(egM,'b','linewidth',lnp)
 plot(coM,'r','linewidth',lnp)
mean(deM(1:20))
mean(egM(1:20))
mean(coM(1:20))


 ylabel('heat flux [W/m^2]')
% xlabel('year')

print -depsc ../ps/heat.eps
