clear
close all
load MOD20090404.mat  
%load MOD20130422.mat
MOD=data;
load PM_stacje.mat
STAT=data;
%load PM25_2013.mat  
load PM_2013.mat
for k=1:length(data.name)
 for i=1:length(STAT.name)
  name=cell2mat(STAT.name(i));
  name=name(2:end-2);
  if strcmp(cell2mat(data.name(k)),name);  
   data.lon(k)=STAT.lon(i);
   data.lat(k)=STAT.lat(i);
   break   
  end
 end
end

AOD=interp2(MOD.lon,MOD.lat,MOD.AOT,data.lon,data.lat)
ind=find(data.month==4&data.day==9&data.hour>11&data.hour<13);
PM25=data.PM10(ind,:);
plot(AOD,PM25,'o')

figure
pcolor(MOD.lon,MOD.lat,MOD.AOT)
plotcoast('k')
shading flat
colorbar

figure
ind=find(isnan(PM25)==0);
plotcolor(data.lon(ind),data.lat(ind),PM25(ind))
axis([14 25 48 56])
colorbar
plotcoast('k')
