208 Computational Modelling in Hydraulic and Coastal Engineering
csource=1;
dc=0.00001;
imm=12;
jmm=12;
Dd=2000;
Dt=600;
nx=22;
ny=23;
nt=288;
% 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 ThermD.txt
load ThermU.txt
load ThermV.txt
load ThermE.txt
fDep=ThermD;
fCoU=ThermU;
fCoV=ThermV;
fEdd=ThermE;
% 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;