Contaminant and sediment transport by advection and diffusion 251
hn(i,j)=1;
u(i,j)=0;
un(i,j)=0;
v(i,j)=0;
vn(i,j)=0;
c(i,j)=0;
co(i,j)=0;
end
end
uu=0;
vv=0;
% Number of particles entering the field per unit time;
ipp=fix(ip/nt);
vol=qg*Dt/ipp;
for i=1:ipp
x(i)=xo;
y(i)=yo;
end
for k=1:nt
ipt=ipp*k;
if ipt>ip
ipt=ip;
end
D=sqrt(6*ed/Dt);
for i=ipt-ipp+1:ipt
x(i)=xo;
y(i)=yo;
end
% Estimation of pressure head;
for j=2:ny-2
for i=2:nx-2
hn(i,j)=h(i,j)-Dt/Dd/2*(u(i+1,j)-u(i,j)+v(i,j+1)v(i,j))+(c(i,j)-co(i,j))*vol/Dd^2;
end
end
for j=1:ny-1
hn(1,j)=hn(2,j);
hn(nx-1,j)=hn(nx-2,j);
end
for i=1:nx-1
hn(i,1)=hn(i,2);
hn(i,ny-1)=hn(i,ny-2);
end
for j=1:ny-1
for i=1:nx-1
h(i,j)=hn(i,j);
co(i,j)=c(i,j);
end
end
% Computation of velocities along the x-axis;