clear all
close all
dz=nc_varget('H04S5_3d_2007.09.nc','dz');
h04s5t1=nc_varget('H04S5T1_trc_2007.09.nc','trc');
h04s5=0;
for i=1:29
	h04s5=h04s5+nc_varget('H04S5_trc_2007.09.nc',['trc',int2str(i)]);
end

figure(1)
pcolor(squeeze(h04s5(1,:,:)));shading flat;colorbar;
% caxis([0 15000]);
title('29 tracer accumlated concentration');
figure(2)
pcolor(squeeze(h04s5t1(1,:,:)));shading flat;colorbar;
% caxis([0 15000]);
title('one single tracer concentration');
figure(3)
diff=h04s5-h04s5t1;
pcolor(squeeze(diff(1,:,:)));shading flat;colorbar;
title('differences between accumulated and single one');
% caxis([-100 100]);

%concentration percent
sum1=sum(sum(nansum(h04s5t1)));
sum2=sum(sum(nansum(diff)));
perc=abs(sum2/sum1)*100
% perc =%     0.0567

%mass percent
sum1=sum(sum(nansum(h04s5t1.*dz)));
sum2=sum(sum(nansum(diff.*dz)));
perc=abs(sum2/sum1)*100
% perc =%    0.0166