fcslope1.m

来自「波浪数值模拟」· M 代码 · 共 37 行

dx = 0.5;
dt = 0.1;

load fCslope1Ecr.dat;   EC = fCslope1Ecr;  
load fCslope1eddycr.dat;   NUC = fCslope1eddycr;  

[n,m] = size(EC);

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

load fCslope2.depth;   h2 = fCslope2;

close all

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


for i=2:n,
   subplot(3,1,1)
   plot(xeta,EC(i,:),[18 18],[-1 1],'r--',[230 230],[-1 1],'r--'),grid
   ylabel('\eta (m)')
   axis([0 max(xeta) -0.75 1])
   s = sprintf('fCslope1.init:  Time = %4.1f (sec)',dt*(i-1));
   title(s);
   subplot(3,1,2)
   plot(xeta,NUC(i,:),[18 18],[-1 1],'r--',[230 230],[-1 1],'r--'),grid
   ylabel('eddy viscosity (m^2/s)')
   axis([0 max(xeta) -0.5 2])
 %  s = sprintf('fCslope1.init:  Time = %4.1f (sec)',dt*(i-1));
 %  title(s);
   pause(0.1);
end;