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

DUST_temp=ncread('~/Aerosol_Variaiton/IRUN/read_nc/ncfile/read_nc_DUST.nc','DUST')*1e6;
DUST_temp(isnan(DUST_temp))=0;


DUST=zeros(144,96);
DUST(:,:)=mean(mean(DUST_temp,4),3);

DUST_Global=sum(sum(DUST.*area,1),2)/1e3/1e12;

DUST_Global_mean=sum(sum(DUST.*area,1),2)/sum(area(:));

