fcshoal.m

波浪数值模拟

load fCshoal.depth;

dt = 1;

h = fCshoal;
clear fCshoal;

%load eta_cross.dat;   EC = eta_cross;
%[n,m] = size(EC);

dx = 1;

xeta = [0:399]*dx + dx/2;
xv = [0:401]*dx - dx/2;
xu = [0:400]*dx;

dx = 1;
dy = 1;
y = [0:99]*dy;
   
nn = 100;

%for i=0:nn,
%   ss = sprintf('load etashoal_%04d.dat;',i);
%   eval(ss);
%   ss = sprintf('load Vshoal_%04d.dat;  ',i);
%   eval(ss);
%   ss = sprintf('load Ushoal_%04d.dat;  ',i);
 %  eval(ss);
%   ss = sprintf('load vortshoal_%04d.dat;  ',i);
%   eval(ss);
%end;   

figure(1),clf;

for i =1:100;
   t = i*dt;
   ss = sprintf('E = etashoal_%04d;',i,i);
   eval(ss);
   ss = sprintf('V = Vshoal_%04d;',i,i);
   eval(ss);
   ss = sprintf('U = Ushoal_%04d;',i,i);
   eval(ss);
   ss = sprintf('VORT = vortshoal_%04d;',i,i);
   eval(ss);
   subplot(2,2,1)
   contour(xeta,y,h')
   hold on
   pcolor(xeta,y,E'), shading flat
   caxis([-0.04 0.04])
   colorbar('vert');
   xlabel('x (m)')
   ylabel('y (m)')
   ss = sprintf('fCshoal:  t = %4.2f sec',t);
   title(ss)
   hold off
   subplot(2,2,2)
   pcolor(xu,y,VORT'), shading flat, colorbar('vert');
   xlabel('x (m)')
   ylabel('y (m)')
   title('Vorticity (1/s)')
   subplot(2,2,3)
   pcolor(xv,y,V'), shading flat, colorbar('vert');
   xlabel('x (m)')
   ylabel('y (m)')
   title('V (m/s)')
   subplot(2,2,4)
   pcolor(xu,y,U'), shading flat, colorbar('vert');
   xlabel('x (m)')
   ylabel('y (m)')
   title(' U (m/s)')
   pause(0.1);
%   pause
end;