% cd ~/mitgcm/bbl/MITgcm

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% compare baseline, downslope, and bbl
for ts=216:216:21600
 s1=35+readbin(['run/SALTanom.' myint2str(ts,10) '.data'],[25 15 25]);
 s2=35+readbin(['run.down_slope/SALTanom.' myint2str(ts,10) '.data'],[25 15 25]);
 s3=35+readbin(['run.bbl/SALTanom.' myint2str(ts,10) '.data'],[25 15 25]);
 clf, subplot(311)
 mypcolor(1:15,25:-1:1,squeeze(s1(10,:,:))');
 caxis([30 30.1]), thincolorbar
 title(['baseline salinity section on day ' int2str(ts*1200/60/60/24)])
 subplot(312)
 mypcolor(1:15,25:-1:1,squeeze(s2(10,:,:))');
 caxis([30 30.1]), thincolorbar
 title(['downslope salinity section on day ' int2str(ts*1200/60/60/24)])
 subplot(313)
 mypcolor(1:15,25:-1:1,squeeze(s3(10,:,:))');
 caxis([30 30.1]), thincolorbar
 title(['bbl salinity section on day ' int2str(ts*1200/60/60/24)])
 pause(.1)
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% look at some diagnostics from bbl integration
kbot=ones(25,15);
s1=35+readbin(['run.bbl/SALTanom.' myint2str(216,10) '.data'],[25 15 25]);
for k=2:25
 kbot(find(s1(:,:,k)>0))=k;
end
for ts=216:216:21600
 s1=35+readbin(['run.bbl/SALTanom.' myint2str(ts,10) '.data'],[25 15 25]);
 clf, subplot(321), caxis([30 30.1]), mypcolor(1:15,25:-1:1,squeeze(s1(10,:,:))');
 caxis([30 30.1]), thincolorbar
 title(['salinity section on day ' int2str(ts*1200/60/60/24)])
 s2=kbot; for i=1:25, for j=1:15, s2(i,j)=s1(i,j,kbot(i,j)); end, end
 subplot(322), mypcolor(s2'); caxis([30 30.1]), thincolorbar
 title('salinity at bottom grid')
 s3=readbin(['run.bbl/BBLsalt.' myint2str(ts,10) '.data'],[25 15]);
 subplot(323), mypcolor(s3'); caxis([30 30.1]), thincolorbar
 title('bbl salinity')
 subplot(324), mypcolor(s3'-s2'); thincolorbar
 title('bbl salinity minus bottom salinity')
 ts=readbin(['run.bbl/BBLtendS.' myint2str(ts,10) '.data'],[25 15]);
 subplot(325), mypcolor(ts'); thincolorbar
 title('salinity tendency'), pause(.1)
end