matlab_class/sample_analysis/plot_std_field.m

function []=plot_std_field(choiceV3orV4);
%object:    compute and display a standard deviation field
%inputs:    choiceV3orV4 ('v3' or 'v4') selects the sample GRID

input_list_check('plot_std_field',nargin);


%%%%%%%%%%%%%%%%%
%load parameters:
%%%%%%%%%%%%%%%%%
gcmfaces_global;
myenv.nctiles=1;
if strcmp(choiceV3orV4,'v4');
    nF=5;
    fileFormat='compact';
    dir0=myenv.gcmfaces_dir;
    dirGrid=fullfile(dir0,'../GRID/');
    myenv.nctilesdir=fullfile(dir0,'sample_input/nctiles_climatology/');
else;
    error('files are missing');
    nF=1;
    dir0=fullfile(myenv.gcmfaces_dir,'sample_input/');
    fileFormat='straight';
    dirGrid=fullfile(dir0,'grid_occa/');
    myenv.nctilesdir=fullfile(dir0,'nctiles_occa/');
end;

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.9	function []=plot_std_field(choiceV3orV4);
2			%object: compute and display a standard deviation field
3			%inputs: choiceV3orV4 ('v3' or 'v4') selects the sample GRID
4
5			input_list_check('plot_std_field',nargin);
6	gforget	1.1
7	gforget	1.11
8	gforget	1.1	%%%%%%%%%%%%%%%%%
9			%load parameters:
10			%%%%%%%%%%%%%%%%%
11	gforget	1.10	gcmfaces_global;
12	gforget	1.12	myenv.nctiles=1;
13	gforget	1.11	if strcmp(choiceV3orV4,'v4');
14			nF=5;
15			fileFormat='compact';
16	gforget	1.12	dir0=myenv.gcmfaces_dir;
17			dirGrid=fullfile(dir0,'../GRID/');
18			myenv.nctilesdir=fullfile(dir0,'sample_input/nctiles_climatology/');
19	gforget	1.11	else;
20			error('files are missing');
21			nF=1;
22	gforget	1.12	dir0=fullfile(myenv.gcmfaces_dir,'sample_input/');
23	gforget	1.11	fileFormat='straight';
24	gforget	1.12	dirGrid=fullfile(dir0,'grid_occa/');
25			myenv.nctilesdir=fullfile(dir0,'nctiles_occa/');
26			end;
27
28			if ~isdir(myenv.nctilesdir);
29			warning(['skipping plot_std_field (missing ' myenv.nctilesdir ')']);
30			return;
31	gforget	1.11	end;
32
33			mygrid=[]; grid_load(dirGrid,nF,fileFormat,1);
34	gforget	1.1
35			%%%%%%%%%%%
36			%get field:
37			%%%%%%%%%%%
38	gforget	1.12	fld=read_nctiles([myenv.nctilesdir 'ETAN/ETAN'],'ETAN');
39	gforget	1.11	fld=std(fld,[],3); fld(find(fld==0))=NaN;
40			cc=[0:0.1:1]*0.10;
41
42	gforget	1.1	%%%%%%%%%%%%
43			%plot field:
44			%%%%%%%%%%%%
45	gforget	1.10	if ~myenv.lessplot;
46	gforget	1.11	gcmfaces_sphere(fld,cc,[],[],1);
47	gforget	1.10	end;
48	gforget	1.1
49