fcyhump.m

波浪数值模拟

dx = 1;
xeta = [0:39]'+dx/2;
xu = [0:40]';
xv = [0:41]' - dx/2;
y = [0:99];

nn = 80;
dt = 1;

for i=0:nn,
   ss = sprintf('load etayhump_%04d.dat;',i);
   eval(ss);
   ss = sprintf('load Vyhump_%04d.dat;  ',i);
   eval(ss);
   ss = sprintf('load Uyhump_%04d.dat;  ',i);
   eval(ss);
end;   

figure(1),clf;

for i =40:nn;
   t = i*dt;
   ss = sprintf('E = etayhump_%04d;',i,i);
   eval(ss);
   ss = sprintf('V = Vyhump_%04d;',i,i);
   eval(ss);
   ss = sprintf('U = Uyhump_%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('fCyhump:  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)')
   plot(y,E(20,:)),grid
   xlabel('y (m)')
   ylabel('\eta (m)')
   pause(0.1);
%   pause
end;