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); |