fchump.m

波浪数值模拟

dx = 0.5;
xeta = [0:399]'+dx;
xu = [0:400]';
xv = [0:401]' - dx;
y = [0:99];

nn = 30;
dt = 1;

%for i=0:nn,
%   ss = sprintf('load etahump_%04d.dat;',i);
%   eval(ss);
%   ss = sprintf('load Vhump_%04d.dat;  ',i);
%   eval(ss);
%   ss = sprintf('load Uhump_%04d.dat;  ',i);
%   eval(ss);
%end;   

figure(1),clf;

for i =0:nn;
   t = i*dt;
   ss = sprintf('E = etahump_%04d;',i,i);
   eval(ss);
   ss = sprintf('V = Vhump_%04d;',i,i);
   eval(ss);
   ss = sprintf('U = Uhump_%04d;',i,i);
   eval(ss);
   subplot(3,1,1)
   pcolor(xeta,y,E'), shading flat, colorbar('vert');
   xlabel('x (m)')
   ylabel('y (m)')
   ss = sprintf('fChump:  t = %4.2f sec',t);
   title(ss)
   subplot(3,1,2)
   pcolor(xv,y,V'), shading flat, colorbar('vert');
   xlabel('x (m)')
   ylabel('y (m)')
   title('V (m/s)')
   subplot(3,1,3)
   pcolor(xu,y,U'), shading flat, colorbar('vert');
   xlabel('x (m)')
   ylabel('y (m)')
   title(' U (m/s)')
   pause(0.1);
%   pause
end;