function plotPP(src,event,figPP,figCO2,CO2,ha)
gx=get(gca,'XLim');
gy=get(gca,'YLim');
f=find(CO2(:,3)>=gx(1) & CO2(:,3)<=gx(2) & CO2(:,2)>=gy(1) & CO2(:,2)<=gy(2));
figure(figPP); scatter(CO2(f,3),CO2(f,2),15,'k','filled');
figure(figCO2); axes(ha(5)); plot(CO2(f,1),CO2(f,3),'k.');
axes(ha(7)); plot(CO2(f,1),CO2(f,2),'k.');