course-idma2016/matlab/mitgcm_plot_plume.m

function []=mitgcm_plot_plume(fil,cc);
%Example:
%    dirExp='MITgcm/verification/tutorial_plume_on_slope/run/';
%    grid_load(dirExp,1,'straight'); 
%    cc=[-0.1 0.02];
%    figure;
%    for tt=300:300:8640; 
%        fil=sprintf('%sT.%010d',dirExp,tt);
%        mitgcm_plot_plume(fil,cc);
%        pause(0.1);
%    end;
%
%note: this requires that the gcmfaces software has been 
% installed and added to the matlab path as explained in
% MITgcm_contrib/gael/matlab_class/gcmfaces.pdf

gcmfaces_global;

nx=320; nr=60;
x=squeeze(mygrid.XC{1})*ones(1,nr);
y=ones(nx,1)*mygrid.RC';
m=squeeze(mygrid.mskC{1});

c=rdmds2gcmfaces(fil);
c=squeeze(c{1});

% figureL;
depthStretchPlot('pcolor',{x,y,m.*c}); shading flat; 
if isempty(whos('cc')); caxis(cc); end;
colorbar; title(fil);

1	gforget	1.1	function []=mitgcm_plot_plume(fil,cc);
2	gforget	1.2	%Example:
3	gforget	1.3	% dirExp='MITgcm/verification/tutorial_plume_on_slope/run/';
4			% grid_load(dirExp,1,'straight');
5	gforget	1.1	% cc=[-0.1 0.02];
6	gforget	1.3	% figure;
7	gforget	1.1	% for tt=300:300:8640;
8	gforget	1.3	% fil=sprintf('%sT.%010d',dirExp,tt);
9	gforget	1.1	% mitgcm_plot_plume(fil,cc);
10			% pause(0.1);
11			% end;
12	gforget	1.2	%
13			%note: this requires that the gcmfaces software has been
14			% installed and added to the matlab path as explained in
15	gforget	1.3	% MITgcm_contrib/gael/matlab_class/gcmfaces.pdf
16	gforget	1.1
17			gcmfaces_global;
18
19			nx=320; nr=60;
20			x=squeeze(mygrid.XC{1})*ones(1,nr);
21			y=ones(nx,1)*mygrid.RC';
22			m=squeeze(mygrid.mskC{1});
23
24			c=rdmds2gcmfaces(fil);
25			c=squeeze(c{1});
26
27			% figureL;
28			depthStretchPlot('pcolor',{x,y,m.*c}); shading flat;
29			if isempty(whos('cc')); caxis(cc); end;
30			colorbar; title(fil);
31