course-idma2015/computing/mitgcm_plot_plume.m

function []=mitgcm_plot_plume(fil,cc);

gcmfaces_global;

if isempty(whos('fil'));
    %cd MITgcm/verification/tutorial_plume_on_slope/run;
    grid_load('./',1,'straight'); 
    cc=[-0.1 0.02];
    figureL;
    for tt=300:300:8640; 
        fil=sprintf('T.%010d',tt);
        mitgcm_plot_plume(fil,cc);
        pause(0.1);
    end;
    return;
end;

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
3			gcmfaces_global;
4
5			if isempty(whos('fil'));
6			%cd MITgcm/verification/tutorial_plume_on_slope/run;
7			grid_load('./',1,'straight');
8			cc=[-0.1 0.02];
9			figureL;
10			for tt=300:300:8640;
11			fil=sprintf('T.%010d',tt);
12			mitgcm_plot_plume(fil,cc);
13			pause(0.1);
14			end;
15			return;
16			end;
17
18			nx=320; nr=60;
19			x=squeeze(mygrid.XC{1})*ones(1,nr);
20			y=ones(nx,1)*mygrid.RC';
21			m=squeeze(mygrid.mskC{1});
22
23			c=rdmds2gcmfaces(fil);
24			c=squeeze(c{1});
25
26			% figureL;
27			depthStretchPlot('pcolor',{x,y,m.*c}); shading flat;
28			if isempty(whos('cc')); caxis(cc); end;
29			colorbar; title(fil);