%
% now load and plot stations
%
  load formgrid.out
  snum=formgrid(:,1);
  lat=formgrid(:,6);
  lon=formgrid(:,7);
  wdep=formgrid(:,9);
  clear formgrid;
  plot(lon,lat,'mo','markersize',20);
  if length(snum)<100;
    for j=1:length(snum);
       text(lon(j),lat(j),['  ',num2str(snum(j))],...
          'horizontalalignment','center','verticalalignment','middle');
%       text(lon(j),lat(j),'/','fontsize',30)
    end;
 end;
 wd=0;
 if wd==1;
   for j=1:length(snum);
       text(lon(j),lat(j),['  ',num2str(wdep(j))],...
       'horizontalalignment','center','verticalalignment','middle');
    end;
end;    
