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;
dir0=[myenv.gcmfaces_dir '/sample_input/'];
if strcmp(choiceV3orV4,'v4');
    nF=5;
    fileFormat='compact';
    dirGrid=[dir0 '../../GRID/'];
    dirIn=[dir0 '/nctiles_climatology/'];
else;
    error('files are missing');
    nF=1;
    fileFormat='straight';
    dirGrid=[dir0 '/GRID' choiceV3orV4  '/'];
    dirIn=[dir0 '/SAMPLE' choiceV3orV4  '/'];
end;

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

%%%%%%%%%%%
%get field:
%%%%%%%%%%%
fld=read_nctiles([dirIn '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			dir0=[myenv.gcmfaces_dir '/sample_input/'];
13	gforget	1.11	if strcmp(choiceV3orV4,'v4');
14			nF=5;
15			fileFormat='compact';
16			dirGrid=[dir0 '../../GRID/'];
17			dirIn=[dir0 '/nctiles_climatology/'];
18			else;
19			error('files are missing');
20			nF=1;
21			fileFormat='straight';
22			dirGrid=[dir0 '/GRID' choiceV3orV4 '/'];
23			dirIn=[dir0 '/SAMPLE' choiceV3orV4 '/'];
24			end;
25
26			mygrid=[]; grid_load(dirGrid,nF,fileFormat,1);
27	gforget	1.1
28			%%%%%%%%%%%
29			%get field:
30			%%%%%%%%%%%
31	gforget	1.11	fld=read_nctiles([dirIn 'ETAN/ETAN'],'ETAN');
32			fld=std(fld,[],3); fld(find(fld==0))=NaN;
33			cc=[0:0.1:1]*0.10;
34
35	gforget	1.1	%%%%%%%%%%%%
36			%plot field:
37			%%%%%%%%%%%%
38	gforget	1.10	if ~myenv.lessplot;
39	gforget	1.11	gcmfaces_sphere(fld,cc,[],[],1);
40	gforget	1.10	end;
41	gforget	1.1
42