clear all; %close all
N=20; %number of points in x-axis
L_Lake=2500; %length of lake (m)
L=10000; %wavelength of atmospheric wave fronts (m)
T=1.25*60^2; %period of atmospheric wave fronts
(seconds)
rho=1000; �nsity of water (kg/m^3)
d_x=L_Lake/N;
d=2; �pth of lake (m)
a=0.01*d; %amplitude of seiche
g=9.81; % Acceleration of gravity
M_n=0; % Manning's n
c=sqrt(g*(d+abs(a)));
alpha=1;% Factor to show effect of exceeding Courant
limit;
d_t=alpha*d_x/c;
Courant=c*d_t/d_x;
M=fix(10000*(1/d_t));
iplot=2; % plots every "iplot" timestep
% Initial conditions;
t=d_t*(0:M-1);
x=-L_Lake/2+d_x*(0:N-1)+d_x/2;% This puts x at eta
points.......
for m=