course-idma2015/computing/mitgcm_plot_adju_128.m

function []=mitgcm_plot_adju_128(fil,cc);

gcmfaces_global;

if isempty(whos('fil'));
    %cd MITgcm/verification/adjustment.128x64x1/run;
    grid_load('./',1,'straight');
    cc=[-1000:100:1000]/20; suf='Eta'; num=1;
    figureL; set(gcf,'Renderer','zbuffer');
    for tt=0:8:24;
        fil={sprintf('%s.%010d',suf,tt),num};
        mitgcm_plot_adju_cs32(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);

clf;
if isempty(whos('cc')); cc=[]; end;
%gcmfaces_sphere(c.*m,cc,[],{100,20},1)
qwckplot(c.*m); gcmfaces_cmap_cbar(cc);
title(fil);
1	function []=mitgcm_plot_adju_128(fil,cc);
2
3	gcmfaces_global;
4
5	if isempty(whos('fil'));
6	%cd MITgcm/verification/adjustment.128x64x1/run;
7	grid_load('./',1,'straight');
8	cc=[-1000:100:1000]/20; suf='Eta'; num=1;
9	figureL; set(gcf,'Renderer','zbuffer');
10	for tt=0:8:24;
11	fil={sprintf('%s.%010d',suf,tt),num};
12	mitgcm_plot_adju_cs32(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	clf;
28	if isempty(whos('cc')); cc=[]; end;
29	%gcmfaces_sphere(c.*m,cc,[],{100,20},1)
30	qwckplot(c.*m); gcmfaces_cmap_cbar(cc);
31	title(fil);