cheapAML_box/input/gendata.m

clear;

ieee='b';
accuracy='real*8';

H=1000;
nx=100;
ny=100;
nz=1;

% =======================  TOPOGRAPHY ========================================= 

% Flat bottom at z=-Ho

h=-H*ones(nx,ny);


% Walls
h(1,:)   = 0.0;
h(end,:) = 0.0;
h(:,1)   = 0.0;
h(:,end) = 0.0;

% land
for n1=1:nx
for n2=1:ny
  if (n2>60 && n1<n2)
    h(n1,n2) = 0.0;
  end
end
end

fid=fopen('topog.box','w',ieee); fwrite(fid,h,accuracy); fclose(fid);

1	clear;
2
3	ieee='b';
4	accuracy='real*8';
5
6	H=1000;
7	nx=100;
8	ny=100;
9	nz=1;
10
11	% ======================= TOPOGRAPHY =========================================
12
13	% Flat bottom at z=-Ho
14
15	h=-H*ones(nx,ny);
16
17
18
19	% Walls
20	h(1,:) = 0.0;
21	h(end,:) = 0.0;
22	h(:,1) = 0.0;
23	h(:,end) = 0.0;
24
25	% land
26	for n1=1:nx
27	for n2=1:ny
28	if (n2>60 && n1<n2)
29	h(n1,n2) = 0.0;
30	end
31	end
32	end
33
34	fid=fopen('topog.box','w',ieee); fwrite(fid,h,accuracy); fclose(fid);
35