for it=1:nt
  subplot(2,1,1)
  plot(x,eta(:,35,it))
  set(gca,'xlim',[x(1) x(end)],'position',[0.13 0.58 0.67 0.35])
  title(['t=',num2str(t(it)/1e4) '/f'])
  subplot(2,1,2)
  pcolor(x,y,eta(:,:,it)');shading interp,h=gca;colorbar
  set(h,'xlim',[x(1) x(end)],'position',[0.13 0.11 0.67 0.35])
  drawnow
  pause(.1)
end

%quiver(x-nx*dx/2,y-ny*dy/2,u,v)