fcshoalsol.m

波浪数值模拟

dx = 0.5;

dt = 1;

load fCshoalsolEcrwk.dat;   EC = fCshoalsolEcrwk;  
load fCshoalsolUcrwk.dat;   UC = fCshoalsolUcrwk;  
load fCshoalsolEcrnwo.dat;   ECn = fCshoalsolEcrnwo;  
load fCshoalsolUcrnwo.dat;   UCn = fCshoalsolUcrnwo;  

load fCshoalsol.depth;   h = fCshoalsol;
clear fCshoalsol;

[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,h),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:133,
   subplot(4,1,1)
   plot(xeta,EC(i,:),xeta,ECn(i,:),'--'),grid
%   plot(xeta,EC(i,:)),grid
   ylabel('\eta (m)')
  axis([0 max(xeta) -0.1 2])
   s = sprintf('fCshoalsol.init:  Time = %4.1f (sec)',dt*(i-1));
   title(s);
   subplot(4,1,2)
   plot(xu,UC(i,:),xu,UCn(i,:),'--'),grid
%   plot(xu,UC(i,:)),grid
   ylabel('U (m/s)')
   axis([0 max(xeta) -3 0.1])
 %  s = sprintf('fCshoalsol.init:  Time = %4.1f (sec)',dt*(i-1));
 %  title(s);
  if (i<5)
     pause(0.08);
  else
     pause
  end;
  
 %   pause;
end;