--- MITgcm_contrib/bbl/lookat_output.m	2010/11/22 22:29:55	1.2
+++ MITgcm_contrib/bbl/lookat_output.m	2011/03/06 02:11:04	1.3
@@ -14,7 +14,7 @@
 % compare baseline, downslope, and bbl
 colormap(cmap), cx=[29.99 30.1];
 s1=zeros(25,15,25); s2=s1; s3=s1;
-for ts=216:216:21600
+for ts=216:216:7200
  s1(:,:,1:25)=35+readbin(['run/SALTanom.' myint2str(ts,10) '.data'],[25 15 25]);
  s2(:,:,1:25)=35+readbin(['run.down_slope/SALTanom.' myint2str(ts,10) '.data'],[25 15 25]);
  s3(:,:,1:25)=35+readbin(['run.bbl/SALTanom.' myint2str(ts,10) '.data'],[25 15 25]);
@@ -30,12 +30,12 @@
  mypcolor(lat,-dpt/1e3,squeeze(s3(10,:,:))'); caxis(cx), thincolorbar
  title(['bbl salinity section on day ' int2str(ts*1200/60/60/24)])
  hold on, plot(lat,-Depth(10,:)/1e3,'r.-'), ylabel('depth (km)')
- xlabel('latitude (deg)'), pause(.1)
+ xlabel('latitude (deg)'), pause(.5)
 end
 
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 % look at some diagnostics from bbl integration
-for ts=216:216:21600
+for ts=216:216:7200
  s1=35+readbin(['run.bbl/SALTanom.' myint2str(ts,10) '.data'],[25 15 25]);
  clf, subplot(321), mypcolor(lat,-dpt/1e3,squeeze(s1(10,:,:))'); caxis(cx), thincolorbar
  title(['salinity section on day ' int2str(ts*1200/60/60/24)])
@@ -61,7 +61,7 @@
   RAC(i,j)=RAC(i,j)*hFacC(i,j,kbot(i,j));
  end
 end
-for ts=216:216:21600
+for ts=216:216:7200
  eta=readbin(['run.bbl/BBLeta.' myint2str(ts,10) '.data'],[25 15]);
  tendS=readbin(['run.bbl/BBLtendS.' myint2str(ts,10) '.data'],[25 15]);
  tmp=RAC.*tendS;