xbdims=getdimensions
[x,H]=readvar1d('H.dat',xbdims);
[x,zs]=readvar1d('zs.dat',xbdims);
[x,zb]=readvar1d('zb.dat',xbdims);
[x,ue]=readvar1d('u.dat',xbdims);
%[x,ve]=readvar1d('ve.dat',xbdims);
figure %(1)
for i=1:xbdims.nt
subplot(411)
title(i)
plot(x,zb(:,i),x,zs(:,i));axis([0 600 -5 6])
subplot(412)
plot(x,zb(:,1),'r',x,zb(:,i));axis([0 600 -2 6])
subplot(413)
plot(x,H(:,i));axis([0 600 0 2.5])
subplot(414)
plot(x,ue(:,i));axis([0 600 -1 1])
%subplot(515)
%plot(x,ve(:,i));axis([0 1000 -1 1])
drawnow
end