% cd ~/mitgcm/bbl/MITgcm_contrib/bbl

% generate bathymetry file - a periodic channel with sloping bottom
bathy=372*ones(25,15);
bathy(:,1)=0;
bathy(:,2)=3000;
bathy(:,3:5)=3600;
for j=6:12
    bathy(:,j)=3600-(3600-432)*(j-5)/7;
end
bathy(:,13)=412;
bathy(:,14)=392;
clf, mypcolor(-bathy'); thincolorbar
writebin('input/bathyFile.bin',-bathy);

% generate initial salinity file
s=30*ones(25,15,25);
s(:,15,3)=30.1;
writebin('input/saltClim.bin',s);