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	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);