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

ADESF_temp=ncread('~/Aerosol_Variaiton/IRUN/read_nc/ncfile/read_nc_ADESF.nc','ADESF');

ADESF=zeros(144,96);
%ADESF(:,:)=mean(mean(ADESF_temp,4),3);
ADESF(:,:)=mean(ADESF_temp(:,:,4,:),4);

ADESF_Global=sum(sum(ADESF.*area,1),2)/sum(area(:));
ADESF_China=sum(sum(ADESF.*area.*china,1),2)/sum(sum(area.*china,1),2);

