load mbs15.mat;
z = mbs15.data;

[r,c] = find(z == -9999);

for i=1:length(r)
    z(r(i),c(i)) = nan;
    if( rem(i,1000) == 0)
        fprintf(1,'%d\n', i)
    end
end
