course-idma2016/matlab/mitgcm_plot_held_suarez.m

function []=mitgcm_plot_held_suarez(fil,cc);
%Example:
%    cd MITgcm/verification/tutorial_held_suarez_cs/run;
%    grid_load('./',6,'cube'); 
%    cc=[260:2:310]; suf='T'; num=1;
%    figureL; set(gcf,'Renderer','zbuffer');
%    for tt=282240:192*30:345600;
%        fil={sprintf('%s.%010d',suf,tt),num};
%        mitgcm_plot_held_suarez(fil,cc);
%        shading interp;
%        pause(0.1);
%    end;

gcmfaces_global;

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

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

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

1	gforget	1.1	function []=mitgcm_plot_held_suarez(fil,cc);
2			%Example:
3			% cd MITgcm/verification/tutorial_held_suarez_cs/run;
4			% grid_load('./',6,'cube');
5			% cc=[260:2:310]; suf='T'; num=1;
6			% figureL; set(gcf,'Renderer','zbuffer');
7			% for tt=282240:192*30:345600;
8			% fil={sprintf('%s.%010d',suf,tt),num};
9			% mitgcm_plot_held_suarez(fil,cc);
10			% shading interp;
11			% pause(0.1);
12			% end;
13
14			gcmfaces_global;
15
16			m=squeeze(mygrid.mskC(:,:,1));
17
18			if ~iscell(fil);
19			c=rdmds2gcmfaces(fil);
20			else;
21			c=rdmds2gcmfaces(fil{1},'rec',fil{2});
22			end;
23			c=c(:,:,1);
24
25			clf;
26			if isempty(whos('cc')); cc=[]; end;
27			%gcmfaces_sphere(c.*m,cc,[],[],1);
28			qwckplot(c.*m); gcmfaces_cmap_cbar(cc);
29			title(fil);
30