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

DUSTSF_temp=ncread('~/Aerosol_Variaiton/IRUN/read_nc/ncfile/read_nc_DUSTSF.nc','DUSTSF');

DUSTSF=zeros(144,96);
DUSTSF(:,:)=mean(mean(DUSTSF_temp,4),3);

DUSTSF_Global=sum(sum(DUSTSF.*area,1),2)*1e3*60*60*24*365/1e12;

DUSTSF_Asia=sum(sum(DUSTSF(29:55,54:80).*area(29:55,54:80),1),2)*1e3*60*60*24*365/1e12;