load ~/NC/grid
load ~/NC/china.mat
load ~/NC/area_1.9_2.5.mat
year=1:150;

ADE_temp=ncread('~/Aerosol_Variaiton/IRUN/read_nc/ncfile/read_nc_ADE.nc','ADE');

ADE=zeros(144,96);
ADE(:,:)=mean(mean(ADE_temp,4),3);
%ADE(:,:)=mean(ADE_temp(:,:,4,:),4);

ADE_Global=sum(sum(ADE.*area,1),2)/sum(area(:));
ADE_China=sum(sum(ADE.*area.*china,1),2)/sum(sum(area.*china,1),2);

