%

% routine to plot broadscale stations

%



load broadsta.dat

%wsn=menu('Do you want station numbers?','y','n');

wsn=input('enter 1 for dots & labels, 2 for just green dots, 3 for just labels ');


if wsn<3;
  plot(broadsta(:,3),broadsta(:,2),'g.','markersize',20);

end;

if wsn~=2;

  for it=1:length(broadsta);

    h=text(broadsta(it,3),broadsta(it,2),num2str(broadsta(it,1)));

    set(h,'fontsize',8,'color','k');

  end;

end;