% This is a matlab script that generates the input data
% taken from cfpaola

% Dimensions of grid
NX=1024;
NY=1024;

% Nominal depth of model (meters)
H=2000;

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

% Topog
% - channel
hbot=-H*ones(NX,NY);
% Walls
% hbot(:,end)=0;
fid=fopen('cg2d_bathy.bin','w',ieee); 
fwrite(fid,hbot,accuracy); 
fclose(fid);