一维水动力matlab程序,matlab 一维湖泊浅水方程，pressure driven seiche code

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=