course-idma2016/matlab/mitgcm_plot_adju_cs32.m

function []=mitgcm_plot_adju_cs32(fil,cc);
%Example:
%    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;
%
%note: this requires that the gcmfaces software has been
% installed and added to the matlab path as explained in
% http://mitgcm.org/viewvc/*checkout*/MITgcm/MITgcm_contrib/gael/matlab_class/gcmfaces.pdf

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,[],{100,20},1)
qwckplot(c.*m); gcmfaces_cmap_cbar(cc);
title(fil);

1	gforget	1.1	function []=mitgcm_plot_adju_cs32(fil,cc);
2			%Example:
3			% cd MITgcm/verification/adjustment.cs-32x32x1/run;
4			% fprintf('\n\n grid files answer is: ''tile*.mitgrid''\n\n');
5			% grid_load('./',6,'cube');
6			% cc=[-1000:100:1000]/500; suf='dynDiag'; num=3;
7			% figureL; set(gcf,'Renderer','zbuffer');
8			% for tt=22:24:190;
9			% fil={sprintf('%s.%010d',suf,tt),num};
10			% mitgcm_plot_adju_cs32(fil,cc);
11			% pause(1);
12			% end;
13	gforget	1.2	%
14			%note: this requires that the gcmfaces software has been
15			% installed and added to the matlab path as explained in
16			% http://mitgcm.org/viewvc/checkout/MITgcm/MITgcm_contrib/gael/matlab_class/gcmfaces.pdf
17	gforget	1.1
18			gcmfaces_global;
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,[],{100,20},1)
32			qwckplot(c.*m); gcmfaces_cmap_cbar(cc);
33			title(fil);
34