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