6 |
|
|
7 |
% Dimensions of grid |
% Dimensions of grid |
8 |
nx=60; |
nx=60; |
9 |
ny=3; |
ny=1; |
10 |
nz=20; |
nz=20; |
11 |
% Nominal depth of model (meters) |
% Nominal depth of model (meters) |
12 |
H=200.0; |
H=200.0; |
49 |
fid=fopen('T.init','w',ieee); fwrite(fid,t,prec); fclose(fid); |
fid=fopen('T.init','w',ieee); fwrite(fid,t,prec); fclose(fid); |
50 |
|
|
51 |
% Sloping channel |
% Sloping channel |
52 |
slope=0.0 |
slope=0.03 |
53 |
offset=2.5e3; |
offset=2.5e3; |
54 |
dmax=-40.0; |
dmax=-40.0; |
55 |
d=0.0*rand([nx,ny]); |
d=0.0*rand([nx,ny]); |
65 |
end |
end |
66 |
end |
end |
67 |
end |
end |
68 |
%d(nx,:)=0.0; |
d(nx,:)=0.0; |
69 |
fid=fopen('topog.slope','w',ieee); fwrite(fid,d,prec); fclose(fid); |
fid=fopen('topog.slope','w',ieee); fwrite(fid,d,prec); fclose(fid); |
70 |
plot(x,d(:,2)) |
plot(x,d(:,1)) |
71 |
|
|
72 |
fid=fopen('delXvar','w',ieee); fwrite(fid,dx,prec); fclose(fid); |
fid=fopen('delXvar','w',ieee); fwrite(fid,dx,prec); fclose(fid); |
73 |
|
|
98 |
end |
end |
99 |
%d1(end-1:end,:)=d1(1:2,:); % debug by aja |
%d1(end-1:end,:)=d1(1:2,:); % debug by aja |
100 |
fid=fopen('topog.convex','w',ieee); fwrite(fid,d1,prec); fclose(fid); |
fid=fopen('topog.convex','w',ieee); fwrite(fid,d1,prec); fclose(fid); |
101 |
hold on |
plot(x,d1(:,1),'g') |
|
plot(x,d1(:,2),'g') |
|
102 |
|
|
103 |
%convex slope |
%convex slope |
104 |
d2=zeros(nx,ny); |
d2=zeros(nx,ny); |
125 |
%d2(end-1:end,:)=d2(1:2,:); % debug by aja |
%d2(end-1:end,:)=d2(1:2,:); % debug by aja |
126 |
fid=fopen('topog.concave','w',ieee); fwrite(fid,d2,prec); fclose(fid); |
fid=fopen('topog.concave','w',ieee); fwrite(fid,d2,prec); fclose(fid); |
127 |
hold on |
hold on |
128 |
plot(x,d2(:,2),'r') |
plot(x,d2(:,1),'r') |
129 |
|
hold off |
130 |
|
|
131 |
|
|
132 |
fid=fopen('delXvar','w',ieee); fwrite(fid,dx,prec); fclose(fid); |
fid=fopen('delXvar','w',ieee); fwrite(fid,dx,prec); fclose(fid); |