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