150 Computational Modelling in Hydraulic and Coastal Engineering
% Input data;
g=9.81;
ho=4;
fbo=0.001;
zbmax=2;
nd=4;
Dx=2;
Dt=0.1;
nx=200;
nt=240;
r=g*(Dt/Dx)^2/2;
nx2=nx/2;
% Option for variable bed friction;
for i=1:nx
fb(i)=fbo;
end
for i=1:nx
z(i)=0;
zo(i)=0;
zm(i)=0;
zn(i)=0;
h(i)=ho;
bed(i)=0;
end
0.45
0.35
0.25
0.15
0.05
–0.05
0
0.4
0.3
0.2
0.1
0
20
40
60
80
100
120
140
160
180
200
Number of spatial steps × 2 (m)
Location of seafloor movement
Wave height (m)
Time = nt*Dt
Time = nt*Dt/2
Figure 6.5 Propagation of a tsunami wave.