180 Computational Modelling in Hydraulic and Coastal Engineering
% Subroutine N=fm(i,j,index): For definition of the b.c.
% indices nn;
nx=30;
ny=30;
Dx=50;
Dy=10;
nn=1500;
% Initialization of the variable;
for i=1:nx
for j=1:ny
f(i,j)=0;
end
end
chk=0;
for k=1:nn
for i=1:nx
for j=1:ny
N=fm(i,j,index);
if N>0
chk=1000;
end
if chk~=1000
No=fm(i,j+1,index);
Nu=fm(i,j-1,index);
Nr=fm(i+1,j,index);
Nl=fm(i-1,j,index);
fo=f(i,j+1);
fu=f(i,j-1);
fr=f(i+1,j);
fl=f(i-1,j);
% Impervious boundary;
if No==1
fo=f(i,j);
end
if Nu==1
fu=f(i,j);
end
if Nr==1
fr=f(i,j);
end
if Nl==1
fl=f(i,j);
end
% Constant head boundary;
if No>1
fo=No;
end
if Nu>1
fu=Nu;
end