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