/[MITgcm]/MITgcm_contrib/bbl/lookat_output.m
ViewVC logotype

Diff of /MITgcm_contrib/bbl/lookat_output.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.2 by dimitri, Mon Nov 22 22:29:55 2010 UTC revision 1.3 by dimitri, Sun Mar 6 02:11:04 2011 UTC
# Line 14  end Line 14  end
14  % compare baseline, downslope, and bbl  % compare baseline, downslope, and bbl
15  colormap(cmap), cx=[29.99 30.1];  colormap(cmap), cx=[29.99 30.1];
16  s1=zeros(25,15,25); s2=s1; s3=s1;  s1=zeros(25,15,25); s2=s1; s3=s1;
17  for ts=216:216:21600  for ts=216:216:7200
18   s1(:,:,1:25)=35+readbin(['run/SALTanom.' myint2str(ts,10) '.data'],[25 15 25]);   s1(:,:,1:25)=35+readbin(['run/SALTanom.' myint2str(ts,10) '.data'],[25 15 25]);
19   s2(:,:,1:25)=35+readbin(['run.down_slope/SALTanom.' myint2str(ts,10) '.data'],[25 15 25]);   s2(:,:,1:25)=35+readbin(['run.down_slope/SALTanom.' myint2str(ts,10) '.data'],[25 15 25]);
20   s3(:,:,1:25)=35+readbin(['run.bbl/SALTanom.' myint2str(ts,10) '.data'],[25 15 25]);   s3(:,:,1:25)=35+readbin(['run.bbl/SALTanom.' myint2str(ts,10) '.data'],[25 15 25]);
# Line 30  for ts=216:216:21600 Line 30  for ts=216:216:21600
30   mypcolor(lat,-dpt/1e3,squeeze(s3(10,:,:))'); caxis(cx), thincolorbar   mypcolor(lat,-dpt/1e3,squeeze(s3(10,:,:))'); caxis(cx), thincolorbar
31   title(['bbl salinity section on day ' int2str(ts*1200/60/60/24)])   title(['bbl salinity section on day ' int2str(ts*1200/60/60/24)])
32   hold on, plot(lat,-Depth(10,:)/1e3,'r.-'), ylabel('depth (km)')   hold on, plot(lat,-Depth(10,:)/1e3,'r.-'), ylabel('depth (km)')
33   xlabel('latitude (deg)'), pause(.1)   xlabel('latitude (deg)'), pause(.5)
34  end  end
35    
36  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
37  % look at some diagnostics from bbl integration  % look at some diagnostics from bbl integration
38  for ts=216:216:21600  for ts=216:216:7200
39   s1=35+readbin(['run.bbl/SALTanom.' myint2str(ts,10) '.data'],[25 15 25]);   s1=35+readbin(['run.bbl/SALTanom.' myint2str(ts,10) '.data'],[25 15 25]);
40   clf, subplot(321), mypcolor(lat,-dpt/1e3,squeeze(s1(10,:,:))'); caxis(cx), thincolorbar   clf, subplot(321), mypcolor(lat,-dpt/1e3,squeeze(s1(10,:,:))'); caxis(cx), thincolorbar
41   title(['salinity section on day ' int2str(ts*1200/60/60/24)])   title(['salinity section on day ' int2str(ts*1200/60/60/24)])
# Line 61  for i=1:25 Line 61  for i=1:25
61    RAC(i,j)=RAC(i,j)*hFacC(i,j,kbot(i,j));    RAC(i,j)=RAC(i,j)*hFacC(i,j,kbot(i,j));
62   end   end
63  end  end
64  for ts=216:216:21600  for ts=216:216:7200
65   eta=readbin(['run.bbl/BBLeta.' myint2str(ts,10) '.data'],[25 15]);   eta=readbin(['run.bbl/BBLeta.' myint2str(ts,10) '.data'],[25 15]);
66   tendS=readbin(['run.bbl/BBLtendS.' myint2str(ts,10) '.data'],[25 15]);   tendS=readbin(['run.bbl/BBLtendS.' myint2str(ts,10) '.data'],[25 15]);
67   tmp=RAC.*tendS;   tmp=RAC.*tendS;
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22