function [L] = draw_coast(fac)

%load('/home/jmc/matlab/MATLAB/m_coasts.mat');
load('/u/u2/jmc/MATLAB/m_coasts.mat');

  X=ncst(:,1); Y=ncst(:,2);
  for i=1:length(k)-1,
    xl=X(k(i)+1:k(i+1)-1);
    fk=finite(xl);
    if any(fk),
      yl=Y(k(i)+1:k(i+1)-1);
      nx=length(xl);     
      L(i)=line(fac*xl,fac*yl);
     end
  end
 set(L,'Color',[0 0 0],'LineWidth',2);
% set(L,'LineStyle','-'); 
return