function []=mitgcm_plot_adv_xy(fil,cc);

gcmfaces_global;

if isempty(whos('fil'));
    %cd MITgcm/verification/advect_xy/run;
    grid_load('./',1,'straight');
    cc=[-0.25:0.05:1]/10; suf='T'; num=1;
    figureL; set(gcf,'Renderer','zbuffer');
    for tt=0:16:80;
        fil={sprintf('%s.%010d',suf,tt),num};
        mitgcm_plot_adv_xy(fil,cc);
        %fprintf('paused: hit enter to proceed\n'); pause;
        pause(1);
    end;
    return;
end;

m=squeeze(mygrid.mskC(:,:,1));

if ~iscell(fil);
    c=rdmds2gcmfaces(fil);
else;
    c=rdmds2gcmfaces(fil{1},'rec',fil{2});
end;
c=c(:,:,1);

% figureL;
clf;
if isempty(whos('cc')); cc=[]; end;
qwckplot(c.*m); gcmfaces_cmap_cbar(cc);
title(fil);