load ~/NC/grid
lat_MODIS=-89.5:1:89.5;
lon_MODIS=-179.5:1:179.5;

AOD_MODEL_temp=ncread('~/Aerosol_Variaiton/IRUN/AOD/ncfile/AOD.nc','AOD_ANN');
AOD_MODIS_temp=ncread('~/Aerosol_Variaiton/IRUN/AOD/ncfile/AOD_MODIS.nc','AOD_ANN');
AOD_MODEL_temp=mean(AOD_MODEL_temp,3);
AOD_MODIS_temp=mean(AOD_MODIS_temp,3);

ocean=ncread('/Users/yangyang/NC/OCEAN.nc','ocean');

n=0;
for i=1:144
    for j=17:80
        [~,lon_MODIS_flag]=min(abs(lon_MODIS-lon(i)));
        [~,lat_MODIS_flag]=min(abs(lat_MODIS-lat(j)));
        if (ocean(i,j)==1 && AOD_MODEL_temp(i,j)<0.3 && AOD_MODEL_temp(i,j)>0 && AOD_MODIS_temp(lon_MODIS_flag,lat_MODIS_flag)<0.3 && AOD_MODIS_temp(lon_MODIS_flag,lat_MODIS_flag)>0)
            n=n+1;
            AOD_MODEL(n)=AOD_MODEL_temp(i,j);
            AOD_MODIS(n)=AOD_MODIS_temp(lon_MODIS_flag,lat_MODIS_flag);
        end
    end
end

