course-idma2015/computing/mitgcm_plot_held_suarez.m

function []=mitgcm_plot_held_suarez(fil,cc);

gcmfaces_global;

if isempty(whos('fil'));
    %cd MITgcm/verification/tutorial_held_suarez_cs/run;
    grid_load('./',6,'cube'); 
%     cc=[-1000:100:1000]; suf='Eta'; num=1;
    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;
    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);

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