1 |
% cd ~/mitgcm/bbl/MITgcm |
2 |
|
3 |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
4 |
% compare baseline, downslope, and bbl |
5 |
for ts=216:216:21600 |
6 |
s1=35+readbin(['run/SALTanom.' myint2str(ts,10) '.data'],[25 15 25]); |
7 |
s2=35+readbin(['run.down_slope/SALTanom.' myint2str(ts,10) '.data'],[25 15 25]); |
8 |
s3=35+readbin(['run.bbl/SALTanom.' myint2str(ts,10) '.data'],[25 15 25]); |
9 |
clf, subplot(311) |
10 |
mypcolor(1:15,25:-1:1,squeeze(s1(10,:,:))'); |
11 |
caxis([30 30.1]), thincolorbar |
12 |
title(['baseline salinity section on day ' int2str(ts*1200/60/60/24)]) |
13 |
subplot(312) |
14 |
mypcolor(1:15,25:-1:1,squeeze(s2(10,:,:))'); |
15 |
caxis([30 30.1]), thincolorbar |
16 |
title(['downslope salinity section on day ' int2str(ts*1200/60/60/24)]) |
17 |
subplot(313) |
18 |
mypcolor(1:15,25:-1:1,squeeze(s3(10,:,:))'); |
19 |
caxis([30 30.1]), thincolorbar |
20 |
title(['bbl salinity section on day ' int2str(ts*1200/60/60/24)]) |
21 |
pause(.1) |
22 |
end |
23 |
|
24 |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
25 |
% look at some diagnostics from bbl integration |
26 |
kbot=ones(25,15); |
27 |
s1=35+readbin(['run.bbl/SALTanom.' myint2str(216,10) '.data'],[25 15 25]); |
28 |
for k=2:25 |
29 |
kbot(find(s1(:,:,k)>0))=k; |
30 |
end |
31 |
for ts=216:216:21600 |
32 |
s1=35+readbin(['run.bbl/SALTanom.' myint2str(ts,10) '.data'],[25 15 25]); |
33 |
clf, subplot(321), caxis([30 30.1]), mypcolor(1:15,25:-1:1,squeeze(s1(10,:,:))'); |
34 |
caxis([30 30.1]), thincolorbar |
35 |
title(['salinity section on day ' int2str(ts*1200/60/60/24)]) |
36 |
s2=kbot; for i=1:25, for j=1:15, s2(i,j)=s1(i,j,kbot(i,j)); end, end |
37 |
subplot(322), mypcolor(s2'); caxis([30 30.1]), thincolorbar |
38 |
title('salinity at bottom grid') |
39 |
s3=readbin(['run.bbl/BBLsalt.' myint2str(ts,10) '.data'],[25 15]); |
40 |
subplot(323), mypcolor(s3'); caxis([30 30.1]), thincolorbar |
41 |
title('bbl salinity') |
42 |
subplot(324), mypcolor(s3'-s2'); thincolorbar |
43 |
title('bbl salinity minus bottom salinity') |
44 |
ts=readbin(['run.bbl/BBLtendS.' myint2str(ts,10) '.data'],[25 15]); |
45 |
subplot(325), mypcolor(ts'); thincolorbar |
46 |
title('salinity tendency'), pause(.1) |
47 |
end |