course-idma2016/matlab/mitgcm_plot_adju_cs32.m

function []=mitgcm_plot_adju_cs32(fil,cc);
%Example:
%    dirExp='MITgcm/verification/adjustment.cs-32x32x1/run/';
%    grid_load(dirExp,6,'cube');%awnser: 'tile*.mitgrid'
%    cc=[-1000:100:1000]/500; 
%    fil0=[dirExp 'dynDiag']; rec=3;
%    figure; set(gcf,'Renderer','zbuffer');
%    for tt=22:24:190;
%        fil={sprintf('%s.%010d',fil0,tt),rec};
%        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
% MITgcm_contrib/gael/matlab_class/gcmfaces.pdf

typePlot=0;
if ~isempty(which('gcmfaces_sphere')); typePlot=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;
if typePlot==1;
  gcmfaces_sphere(c.*m,cc,[],{100,20},1);
else;
  qwckplot(c.*m); gcmfaces_cmap_cbar(cc);
end;
title(fil);

1	gforget	1.1	function []=mitgcm_plot_adju_cs32(fil,cc);
2			%Example:
3	gforget	1.3	% dirExp='MITgcm/verification/adjustment.cs-32x32x1/run/';
4			% grid_load(dirExp,6,'cube');%awnser: 'tile*.mitgrid'
5			% cc=[-1000:100:1000]/500;
6			% fil0=[dirExp 'dynDiag']; rec=3;
7			% figure; set(gcf,'Renderer','zbuffer');
8	gforget	1.1	% for tt=22:24:190;
9	gforget	1.3	% fil={sprintf('%s.%010d',fil0,tt),rec};
10	gforget	1.1	% 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	gforget	1.3	% MITgcm_contrib/gael/matlab_class/gcmfaces.pdf
17
18			typePlot=0;
19			if ~isempty(which('gcmfaces_sphere')); typePlot=1; end;
20	gforget	1.1
21			gcmfaces_global;
22
23			m=squeeze(mygrid.mskC(:,:,1));
24
25			if ~iscell(fil);
26			c=rdmds2gcmfaces(fil);
27			else;
28			c=rdmds2gcmfaces(fil{1},'rec',fil{2});
29			end;
30			c=c(:,:,1);
31
32			clf;
33			if isempty(whos('cc')); cc=[]; end;
34	gforget	1.3	if typePlot==1;
35			gcmfaces_sphere(c.*m,cc,[],{100,20},1);
36			else;
37			qwckplot(c.*m); gcmfaces_cmap_cbar(cc);
38			end;
39	gforget	1.1	title(fil);
40