fcsolitaryy.m

波浪数值模拟

dx = 0.5;

dt = 1;

load fCsolitaryyEal.dat;   EL = fCsolitaryyEal;  
load fCsolitaryyVal.dat;   VL = fCsolitaryyVal;  


[n,m] = size(EL);

y = [0:m-1]*dx;

h2 = 0.45*ones(1,m);

ELmax = max(EL,[],2);
time = [1:n]*dt;

figure(2),clf, orient tall


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

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


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