course-idma2015/computing/mitgcm_plot_adv_cs.m

function []=mitgcm_plot_adv_cs(fil,cc);

gcmfaces_global;

if isempty(whos('fil'));
    %cd MITgcm/verification/advect_cs/run;
    grid_load('./',6,'cube');
    cc=[-0.25:0.05:1]*2; suf='ts_Diag'; num=1;
    figureL; set(gcf,'Renderer','zbuffer');
    for tt=0:32:160;
        fil={sprintf('%s.%010d',suf,tt),num};
        mitgcm_plot_adv_cs(fil,cc);
        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);
gcmfaces_sphere(c.*m,cc,[],{-100,20},1);
title(fil);
1	gforget	1.1	function []=mitgcm_plot_adv_cs(fil,cc);
2
3			gcmfaces_global;
4
5			if isempty(whos('fil'));
6			%cd MITgcm/verification/advect_cs/run;
7			grid_load('./',6,'cube');
8			cc=[-0.25:0.05:1]*2; suf='ts_Diag'; num=1;
9			figureL; set(gcf,'Renderer','zbuffer');
10			for tt=0:32:160;
11			fil={sprintf('%s.%010d',suf,tt),num};
12			mitgcm_plot_adv_cs(fil,cc);
13			pause(1);
14			end;
15			return;
16			end;
17
18			m=squeeze(mygrid.mskC(:,:,1));
19
20			if ~iscell(fil);
21			c=rdmds2gcmfaces(fil);
22			else;
23			c=rdmds2gcmfaces(fil{1},'rec',fil{2});
24			end;
25			c=c(:,:,1);
26
27			% figureL;
28			clf;
29			if isempty(whos('cc')); cc=[]; end;
30			% qwckplot(c.*m); gcmfaces_cmap_cbar(cc);
31			gcmfaces_sphere(c.*m,cc,[],{-100,20},1);
32			title(fil);