| 1 |
dimitri |
1.1 |
% cd ~/mitgcm/bbl/MITgcm_contrib/bbl |
| 2 |
|
|
|
| 3 |
dimitri |
1.2 |
% generate bathymetry file - a periodic channel with sloping bottom |
| 4 |
dimitri |
1.1 |
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 |
dimitri |
1.2 |
% generate sea surface salinity file |
| 17 |
dimitri |
1.1 |
sss=30*ones(25,15,2); |
| 18 |
|
|
sss(:,15,:)=30.1; |
| 19 |
|
|
clf, mypcolor(sss(:,:,2)'); thincolorbar |
| 20 |
|
|
writebin('input/saltClim.bin',sss); |