
function yNew = qcM1_NaNcb(x,y,ha)

yNew = y;
axes(ha);
gx = get(ha,'XLim');
gy = get(ha,'YLim');
f = find(x>=gx(1) & x<=gx(2) & y>=gy(1) & y<=gy(2));
plot(x(f),y(f),'r.');
yNew(f) = NaN;
