    load gc80m.bty
    id=find(gc80m(:,1)<0);
    plot(gc80m(id,1),gc80m(id,2),'w.','markersize',4);
