function plot_cbc_cont


g=loadgrid('MESH/stw2');
x=read_cbc('BCS_CBC0.cbc');
d(1)=0;
for k=2:length(x.node),
    d(k)=d(k-1)+abs(  g.x(x.node(k))-g.x(x.node(k-1)) +i*(g.y(x.node(k))-g.y(x.node(k-1)) ));
end

t=[1:x.NTS]'*x.DT/(24*3600);
clf;
subplot(2,1,1);
%pcolor(d/1000,t,x.f);
%shading interp;caxis([-.02,.02]);colorbar 'v';
[cs,h]=contour(d/1000,t,x.f,[-.025:.0025:.025],'k');
clabel(cs,h);
%caxis([-.025,.025]);
%colormap('jet(20)');colorbar 'v';
%xlabel('distance along boundary (km)');
ylabel('time (days)');
title('1');

x1=read_cbc('BCS_CBC1.cbc');
subplot(2,1,2);
%pcolor(d/1000,t,x.f+x1.f);
%shading interp;caxis([-.02,.02]);colorbar 'v';
[cs,h]=contour(d/1000,t,x.f+x1.f,[-.025:.0025:.025],'k');
clabel(cs,h);
%caxis([-.025,.025]);
%colormap('jet(20)');colorbar 'v'
xlabel('distance along boundary (km)');
ylabel('time (days)');
title('2');