load /Users/yang146/Aerosol_Variaiton/OBS/ncfile/sites.mat
n=0;
for i=1:length(siteID)
    if lon_site(i)>10000 && lon_site(i)<11000 && lat_site(i)>3000 && lat_site(i)<4500
        n=n+1;
        site_Gobi(n)=siteID(i);
        lon_Gobi(n)=lon_site(i);
        lat_Gobi(n)=lat_site(i);
    end
end

year=2008:2015;
temp=importdata('/Users/yangyang/Aerosol_Variaiton/OBS/dat/2007-2015DJF_DustStorms/stationID.txt');
stationID=temp(:,1);

Gobi_temp=zeros(length(stationID),1);
for i=1:length(stationID)
    for j=1:n
        if stationID(i)==site_Gobi(j)
            Gobi_temp(i)=1;
        end
    end
end
   
        
Gobi=(Gobi_temp==1);


dust_2008_2015=zeros(length(year),90);
load /Users/yangyang/Aerosol_Variaiton/OBS/dat/2007-2015DJF_DustStorms/DustStorm_2008DJF.mat
dust_2008_2015(1,:)=sum(sum(dust_storm_events(Gobi,1:90,5),3),1);
load /Users/yangyang/Aerosol_Variaiton/OBS/dat/2007-2015DJF_DustStorms/DustStorm_2009DJF.mat
dust_2008_2015(2,:)=sum(sum(dust_storm_events(Gobi,1:90,5),3),1);
load /Users/yangyang/Aerosol_Variaiton/OBS/dat/2007-2015DJF_DustStorms/DustStorm_2010DJF.mat
dust_2008_2015(3,:)=sum(sum(dust_storm_events(Gobi,1:90,5),3),1);
load /Users/yangyang/Aerosol_Variaiton/OBS/dat/2007-2015DJF_DustStorms/DustStorm_2011DJF.mat
dust_2008_2015(4,:)=sum(sum(dust_storm_events(Gobi,1:90,5),3),1);
load /Users/yangyang/Aerosol_Variaiton/OBS/dat/2007-2015DJF_DustStorms/DustStorm_2012DJF.mat
dust_2008_2015(5,:)=sum(sum(dust_storm_events(Gobi,1:90,5),3),1);
load /Users/yangyang/Aerosol_Variaiton/OBS/dat/2007-2015DJF_DustStorms/DustStorm_2013DJF.mat
dust_2008_2015(6,:)=sum(sum(dust_storm_events(Gobi,1:90,5),3),1);
load /Users/yangyang/Aerosol_Variaiton/OBS/dat/2007-2015DJF_DustStorms/DustStorm_2014DJF.mat
dust_2008_2015(7,:)=sum(sum(dust_storm_events(Gobi,1:90,5),3),1);
load /Users/yangyang/Aerosol_Variaiton/OBS/dat/2007-2015DJF_DustStorms/DustStorm_2015DJF.mat
dust_2008_2015(8,:)=sum(sum(dust_storm_events(Gobi,1:90,5),3),1);

lon_site=Lon;
lat_site=Lat;

save('Dust_Gobi_2008_2015_daily','dust_2008_2015','lon_site','lat_site')










                
    