matlab_class/sample_analysis/plot_std_field.m

function []=plot_std_field();
%object:    compute and display a standard deviation field

input_list_check('plot_std_field',nargin);

%%%%%%%%%%%%%%%%%
%load parameters:
%%%%%%%%%%%%%%%%%
gcmfaces_global;
myenv.nctiles=1;
nF=5;
fileFormat='compact';
dirGrid=fullfile(myenv.gcmfaces_dir,'../GRID/');
myenv.nctilesdir=fullfile(myenv.gcmfaces_dir,'sample_input/nctiles_climatology/');

if ~isdir(myenv.nctilesdir);
    warning(['skipping plot_std_field (missing ' myenv.nctilesdir ')']);
    return;
end;

mygrid=[]; grid_load(dirGrid,nF,fileFormat,1);

%%%%%%%%%%%
%get field:
%%%%%%%%%%%
fld=read_nctiles([myenv.nctilesdir 'ETAN/ETAN'],'ETAN');
fld=std(fld,[],3); fld(find(fld==0))=NaN;
cc=[0:0.1:1]*0.10;

%%%%%%%%%%%%
%plot field:
%%%%%%%%%%%%
if ~myenv.lessplot;
    gcmfaces_sphere(fld,cc,[],[],1);
end;


1	gforget	1.13	function []=plot_std_field();
2	gforget	1.9	%object: compute and display a standard deviation field
3
4			input_list_check('plot_std_field',nargin);
5	gforget	1.1
6			%%%%%%%%%%%%%%%%%
7			%load parameters:
8			%%%%%%%%%%%%%%%%%
9	gforget	1.10	gcmfaces_global;
10	gforget	1.12	myenv.nctiles=1;
11	gforget	1.13	nF=5;
12			fileFormat='compact';
13			dirGrid=fullfile(myenv.gcmfaces_dir,'../GRID/');
14			myenv.nctilesdir=fullfile(myenv.gcmfaces_dir,'sample_input/nctiles_climatology/');
15	gforget	1.12
16			if ~isdir(myenv.nctilesdir);
17			warning(['skipping plot_std_field (missing ' myenv.nctilesdir ')']);
18			return;
19	gforget	1.11	end;
20
21			mygrid=[]; grid_load(dirGrid,nF,fileFormat,1);
22	gforget	1.1
23			%%%%%%%%%%%
24			%get field:
25			%%%%%%%%%%%
26	gforget	1.12	fld=read_nctiles([myenv.nctilesdir 'ETAN/ETAN'],'ETAN');
27	gforget	1.11	fld=std(fld,[],3); fld(find(fld==0))=NaN;
28			cc=[0:0.1:1]*0.10;
29
30	gforget	1.1	%%%%%%%%%%%%
31			%plot field:
32			%%%%%%%%%%%%
33	gforget	1.10	if ~myenv.lessplot;
34	gforget	1.11	gcmfaces_sphere(fld,cc,[],[],1);
35	gforget	1.10	end;
36	gforget	1.1
37