12 |
|
|
13 |
namf='channel.bin'; |
namf='channel.bin'; |
14 |
depth=H0*ones(nx,ny); depth(:,1)=0.; |
depth=H0*ones(nx,ny); depth(:,1)=0.; |
15 |
%if kwr > 0, |
if kwr > 0, |
16 |
% fprintf('write to file: %s\n',namf); |
fprintf('write to file: %s\n',namf); |
17 |
% fid=fopen(namf,'w','b'); fwrite(fid,depth,'real*8'); fclose(fid); |
fid=fopen(namf,'w','b'); fwrite(fid,depth,'real*8'); fclose(fid); |
18 |
%end |
end |
19 |
|
|
20 |
namf='bathy_3c.bin'; |
namf='bathy_3c.bin'; |
21 |
msk=abs(xc)'*ones(1,ny)+ones(nx,1)*yc; |
msk=abs(xc)'*ones(1,ny)+ones(nx,1)*yc; |
74 |
%end |
%end |
75 |
%------------------------------------------------------ |
%------------------------------------------------------ |
76 |
|
|
77 |
|
namf='edge100.bin'; w0=1.; |
78 |
|
var=w0*ones(nx,ny); |
79 |
|
var(:,2)=0*var(:,end); |
80 |
|
var(:,3)=0.01*var(:,end-1); |
81 |
|
var(:,4)=0.1*var(:,end-2); |
82 |
|
var(:,end)=0*var(:,end); |
83 |
|
var(:,end-1)=0.01*var(:,end-1); |
84 |
|
var(:,end-2)=0.1*var(:,end-2); |
85 |
|
if kwr >0, |
86 |
|
fprintf('write to file: %s\n',namf); |
87 |
|
fid=fopen(namf,'w','b'); fwrite(fid,var,'real*8'); fclose(fid); |
88 |
|
end |
89 |
|
|
90 |
|
namf='edge+20.bin'; w0=0.2; |
91 |
|
var=w0*ones(nx,ny); |
92 |
|
var(:,2)=0*var(:,end); |
93 |
|
var(:,3)=0.01*var(:,end-1); |
94 |
|
var(:,4)=0.1*var(:,end-2); |
95 |
|
var(:,end)=0*var(:,end); |
96 |
|
var(:,end-1)=0.01*var(:,end-1); |
97 |
|
var(:,end-2)=0.1*var(:,end-2); |
98 |
|
if kwr >0, |
99 |
|
fprintf('write to file: %s\n',namf); |
100 |
|
fid=fopen(namf,'w','b'); fwrite(fid,var,'real*8'); fclose(fid); |
101 |
|
end |
102 |
|
|
103 |
|
%------------------------------------------------------ |
104 |
|
|
105 |
|
|
106 |
dsw0=100; |
dsw0=100; |
107 |
namf=['dsw_',int2str(dsw0),'.bin']; |
namf=['dsw_',int2str(dsw0),'.bin']; |
108 |
fld=dsw0*ones(nx,ny,nt); |
fld=dsw0*ones(nx,ny,nt); |