load /home3/ocn/jmanning/wind/stations.dat
s=stations;

for i=1:length(s);
  disp(['buoy1 ',num2str(s(i,1))]);
  for k=1:length(s);
    di=dist(s(i,3),s(i,2),s(k,3),s(k,2));
    disp(['  buoy2 ',num2str(s(k,1)),' dist = ',num2str(di)])
    pause;
  end;
end;
   