t=[0:10000:600000]';
nt=length(t);
xmax=5025; % km
ymax=3550; % km
dx=25; % km
dy=50; % km

nx=xmax/dx;
ny=ymax/dy;
x=[dx/2:dx:xmax-dx/2]';
y=[dy/2:dy:ymax-dy/2]';

eta=read_binary('eta.bin',[nx+2 ny+2 nt]);
eta=eta(2:end-1,2:end-1,:);
%imagesc(flipud(eta(:,:,5)')),colorbar
u=read_binary('u.bin',[nx+2 ny+2],'float64','ieee-be',nt,(nx+2)*(ny+2));
u=u(2:end-1,2:end-1)';
v=read_binary('v.bin',[nx+2 ny+2],'float64','ieee-be',nt,(nx+2)*(ny+2));
v=v(2:end-1,2:end-1)';
%quiver(x-nx*dx/2,y-ny*dy/2,u,v)