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	jmc	1.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