204 Computational Modelling in Hydraulic and Coastal Engineering
jmm=25;
Dd=5;
Dt=5;
nx=70;
ny=40;
nt=720;
% Initialization of variables;
for i=1:nx
for j=1:ny
c(i,j)=0;
cmax(i,j)=0;
cn(i,j)=0;
end
end
% Importing data for depths, velocities and diffusion
coefficients;
load Depths.txt
load CoastU.txt
load CoastV.txt
load EddyC.txt
fDep=Depths;
fCoU=CoastU;
fCoV=CoastV;
fEdd=EddyC;
% Subrouting for depths h=fm(i,j,fDep);
% Subroutine for velocity u=fm(i,j,fCoU);
% Subroutine for velocity v=fm(i,j,fCoV);
% Subroutine for diffusion coefficients ed=fm(i,j,fEdd);
% Main program;
for k=1:nt
for j=2:ny-2
for i=2:nx-2
h=fm(i,j,fDep);
ho=fm(i,j+1,fDep);
hu=fm(i,j-1,fDep);
hl=fm(i-1,j,fDep);
hr=fm(i+1,j,fDep);
u=fm(i,j,fCoU);
ur=fm(i+1,j,fCoU);
v=fm(i,j,fCoV);
vo=fm(i,j+1,fCoV);
ed=fm(i,j,fEdd);
if h>0.1
advx=0;
if u+ur>0
advx=u*(c(i,j)-c(i-1,j))/Dd;
else
advx=ur*(c(i+1,j)-c(i,j))/Dd;
end
advy=0;