| 1 | % cd ~/mitgcm/bbl/MITgcm_contrib/bbl | 
| 2 |  | 
| 3 | % generate bathymetry file - a periodic channel with sloping bottom | 
| 4 | bathy=372*ones(25,15); | 
| 5 | bathy(:,1)=0; | 
| 6 | bathy(:,2)=3000; | 
| 7 | bathy(:,3:5)=3600; | 
| 8 | for j=6:12 | 
| 9 | bathy(:,j)=3600-(3600-432)*(j-5)/7; | 
| 10 | end | 
| 11 | bathy(:,13)=412; | 
| 12 | bathy(:,14)=392; | 
| 13 | clf, mypcolor(-bathy'); thincolorbar | 
| 14 | writebin('input/bathyFile.bin',-bathy); | 
| 15 |  | 
| 16 | % generate sea surface salinity file | 
| 17 | sss=30*ones(25,15,2); | 
| 18 | sss(:,15,:)=30.1; | 
| 19 | clf, mypcolor(sss(:,:,2)'); thincolorbar | 
| 20 | writebin('input/saltClim.bin',sss); |