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