ecco-godae/input/gendata.m

snx=90;
sny=90;

t0=zeros(snx,sny);
tb=t0;

x = [1:snx]-.5-snx/2;
y = [1:sny]-.5-sny/2;

[xx,yy]=meshgrid(x,y);

a = 1;
sigma=10;

tb = a*exp(-.5*(xx.^2+yy.^2)/sigma);

% mit_writefield('ssh.bump',tb,'real*8');
% mit_writefield('ssh.flat',t0,'real*8');

fid = fopen('ssh.bump', 'w', 'ieee-be');
fwrite(fid,tb,'real*8',0,'ieee-be');
fclose(fid);

fid = fopen('ssh.flat', 'w', 'ieee-be');
fwrite(fid,t0,'real*8',0,'ieee-be');
fclose(fid);

% topography
delR= [50., 70., 100., 140., 190.,...
       240., 290., 340., 390., 440.,...
       490., 540., 590., 640., 690.];
H = sum(delR);

% flat topography
h0 = -H*ones(snx,sny);
% $$$ fid = fopen('bathy.flat', 'w', 'ieee-be');
% $$$ fwrite(fid,h0,'real*8',0,'ieee-be');
% $$$ fclose(fid);

% southern wall
hs=h0;
hs(:,1:35) = 0;
% $$$ fid = fopen('bathy.wall', 'w', 'ieee-be');
% $$$ fwrite(fid,hs,'real*8',0,'ieee-be');
% $$$ fclose(fid);

% write files with proper orientation
for k = 1:17
  fname = sprintf('bathy.%3.3u.001.data',k);
  disp( fname )
  fid = fopen(fname, 'w', 'ieee-be');
  if ismember(k,[1 5])
    fwrite(fid,rot90(hs,0),'real*8',0,'ieee-be');
  end
  if ismember(k,[13 17])
    fwrite(fid,rot90(hs,1),'real*8',0,'ieee-be');
  end
  if ismember(k,[9])
    fwrite(fid,h0,'real*8',0,'ieee-be');
  end
  if ismember(k,[2:4  6:8])
    fwrite(fid,h0,'real*8',0,'ieee-be');
  end
  if ismember(k,[10:12  14:16])
    fwrite(fid,rot90(h0,1),'real*8',0,'ieee-be');
  end
  fclose(fid);
end


1	edhill	1.1	snx=90;
2			sny=90;
3
4			t0=zeros(snx,sny);
5			tb=t0;
6
7			x = [1:snx]-.5-snx/2;
8			y = [1:sny]-.5-sny/2;
9
10			[xx,yy]=meshgrid(x,y);
11
12			a = 1;
13			sigma=10;
14
15			tb = aexp(-.5(xx.^2+yy.^2)/sigma);
16
17			% mit_writefield('ssh.bump',tb,'real*8');
18			% mit_writefield('ssh.flat',t0,'real*8');
19
20			fid = fopen('ssh.bump', 'w', 'ieee-be');
21			fwrite(fid,tb,'real*8',0,'ieee-be');
22			fclose(fid);
23
24			fid = fopen('ssh.flat', 'w', 'ieee-be');
25			fwrite(fid,t0,'real*8',0,'ieee-be');
26			fclose(fid);
27
28			% topography
29			delR= [50., 70., 100., 140., 190.,...
30			240., 290., 340., 390., 440.,...
31			490., 540., 590., 640., 690.];
32			H = sum(delR);
33
34			% flat topography
35			h0 = -H*ones(snx,sny);
36			% $$$ fid = fopen('bathy.flat', 'w', 'ieee-be');
37			% $$$ fwrite(fid,h0,'real*8',0,'ieee-be');
38			% $$$ fclose(fid);
39
40			% southern wall
41			hs=h0;
42			hs(:,1:35) = 0;
43			% $$$ fid = fopen('bathy.wall', 'w', 'ieee-be');
44			% $$$ fwrite(fid,hs,'real*8',0,'ieee-be');
45			% $$$ fclose(fid);
46
47			% write files with proper orientation
48			for k = 1:17
49			fname = sprintf('bathy.%3.3u.001.data',k);
50			disp( fname )
51			fid = fopen(fname, 'w', 'ieee-be');
52			if ismember(k,[1 5])
53			fwrite(fid,rot90(hs,0),'real*8',0,'ieee-be');
54			end
55			if ismember(k,[13 17])
56			fwrite(fid,rot90(hs,1),'real*8',0,'ieee-be');
57			end
58			if ismember(k,[9])
59			fwrite(fid,h0,'real*8',0,'ieee-be');
60			end
61			if ismember(k,[2:4 6:8])
62			fwrite(fid,h0,'real*8',0,'ieee-be');
63			end
64			if ismember(k,[10:12 14:16])
65			fwrite(fid,rot90(h0,1),'real*8',0,'ieee-be');
66			end
67			fclose(fid);
68			end
69
70