function TOPlot(src,event,figTO,figCTD,CTD,ha)
figure(figTO)
gx=get(gca,'XLim');
gy=get(gca,'YLim');
f=find(CTD(:,2)>=gx(1) & CTD(:,2)<=gx(2) & CTD(:,4)>=gy(1) & CTD(:,4)<=gy(2));
figure(figTO); plot(CTD(f,2),CTD(f,4),'k.');
figure(figCTD); axes(ha(3)); plot(CTD(f,1),CTD(f,4),'k.');
figure(figTO);
end