function han = plotclbastcoast(edge)

% get plot state
nextplotstatewas = get(gca,'nextplot');

% hold whatever is already plotted
set(gca,'nextplot','add')

load('/home/wilkin/roms/cblast/grid/cblast_coast_high')

dax = [-71.325 -69.083333 40.466667 41.8];

if nargin > 0
  % set the axis limits to be the full model domain
  axis(dax+edge*[-1 1 -1 1])
end

ax = axis;
lon(find(lon<ax(1)|lon>ax(2))) = NaN;
lat(find(lat<ax(3)|lat>ax(4))) = NaN;
h = plot(lon,lat,'k-');

% add the model domain
hh = plot(dax([1 1 2 2 1]),dax([4 3 3 4 4]));
set(hh,'color',0.8*[1 1 1]);

% restore nextplotstate to what it was
set(gca,'nextplot',nextplotstatewas);

if nargout > 0
   han = h;
end
