fcsolitaryx.m

波浪数值模拟

dx = 0.5;

dt = 1;

load fCsolitaryxEcr.dat;   EC = fCsolitaryxEcr;  
load fCsolitaryxUcr.dat;   UC = fCsolitaryxUcr;  


[n,m] = size(EC);

xeta = [0:m-1]*dx + dx/2;
xu = [0:m]*dx;

h2 = 0.45*ones(1,m);

ECmax = max(EC,[],2);
time = [1:n]*dt;

figure(2),clf, orient tall


subplot(4,1,4)
plot(xeta,h2),grid
set(gca,'ydir','reverse')
xlabel('Cross-shore coordinate (m)')
ylabel('Depth (m)')

subplot(4,1,3)
plot(time,ECmax),grid
xlabel('Time (s)')
ylabel('max \eta (m)')


for i=1:n,
   subplot(4,1,1)
   plot(xeta,EC(i,:)),grid
   ylabel('\eta (m)')
   axis([0 max(xeta) -0.01 0.03])
   s = sprintf('fCsolitaryx.init:  Time = %4.1f (sec)',dt*(i-1));
   title(s);
   subplot(4,1,2)
   plot(xu,UC(i,:)),grid
   ylabel('U (m/s)')
   axis([0 max(xeta) -0.1 0.1])
 %  s = sprintf('fCsolitaryx.init:  Time = %4.1f (sec)',dt*(i-1));
 %  title(s);
    pause(0.05);
 %   pause;
end;