course-idma2015/computing/mitgcm_plot_adju_cs32.m

function []=mitgcm_plot_adju_cs32(fil,cc);

gcmfaces_global;

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