% for modifying the mask

p=ginput;

if(0)
    for i=1:length(p(:,1))
        [I,J] = near2d(rlon,rlat,p(i,1),p(i,2));
        mask2(I,J) = 1;
    end
else
    p = floor(p);
    for i=1:length(p)
        mask2(p(i,2),p(i,1)) = 1;
    end
end


