matlab_class/sample_analysis/example_IO.m

function []=example_IO();
%
% EXAMPLE_IO computes and displays a standard deviation field
%
% stand-alone call:
%   addpath gcmfaces/sample_analysis/; example_IO;
%
% needed input files:
%   mkdir release2
%   wget --recursive ftp://mit.ecco-group.org/ecco_for_las/version_4/release2/nctiles_climatology/ETAN
%   mv mit.ecco-group.org/ecco_for_las/version_4/release2/nctiles_climatology release2_climatology/.

gcmfaces_global;

input_list_check('example_IO',nargin);

%expected location:
myenv.nctilesdir=fullfile(pwd,'/release2_climatology/nctiles_climatology/ETAN/');
%if ETAN is not found then try old locations:
if ~isdir(myenv.nctilesdir);
    %if not found then try old location:
    tmpdir=fullfile(pwd,'/release2_climatology/nctiles_climatology/ETAN/');
    if isdir(tmpdir); myenv.nctilesdir=tmpdir; end;
end;
%if ETAN is still not found then issue warning and skip example_IO
if ~isdir(myenv.nctilesdir);
    help example_IO;
    warning('example_IO requires release2_climatology/nctiles_climatology/ETAN that was not found ---> abort!');
    return;
end;

%%%%%%%%%%%%%%%%%
%load grid:
%%%%%%%%%%%%%%%%%

if isempty(mygrid);
   grid_load;
end;

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

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


1	gforget	1.1	function []=example_IO();
2			%
3			% EXAMPLE_IO computes and displays a standard deviation field
4			%
5	gforget	1.3	% stand-alone call:
6			% addpath gcmfaces/sample_analysis/; example_IO;
7	gforget	1.1	%
8			% needed input files:
9	gforget	1.3	% mkdir release2
10			% wget --recursive ftp://mit.ecco-group.org/ecco_for_las/version_4/release2/nctiles_climatology/ETAN
11			% mv mit.ecco-group.org/ecco_for_las/version_4/release2/nctiles_climatology release2_climatology/.
12	gforget	1.1
13			gcmfaces_global;
14
15			input_list_check('example_IO',nargin);
16
17			%expected location:
18	gforget	1.3	myenv.nctilesdir=fullfile(pwd,'/release2_climatology/nctiles_climatology/ETAN/');
19			%if ETAN is not found then try old locations:
20	gforget	1.1	if ~isdir(myenv.nctilesdir);
21			%if not found then try old location:
22	gforget	1.3	tmpdir=fullfile(pwd,'/release2_climatology/nctiles_climatology/ETAN/');
23	gforget	1.1	if isdir(tmpdir); myenv.nctilesdir=tmpdir; end;
24			end;
25			%if ETAN is still not found then issue warning and skip example_IO
26			if ~isdir(myenv.nctilesdir);
27			help example_IO;
28	gforget	1.3	warning('example_IO requires release2_climatology/nctiles_climatology/ETAN that was not found ---> abort!');
29	gforget	1.1	return;
30			end;
31
32			%%%%%%%%%%%%%%%%%
33			%load grid:
34			%%%%%%%%%%%%%%%%%
35
36			if isempty(mygrid);
37			grid_load;
38			end;
39
40			%%%%%%%%%%%
41			%get field:
42			%%%%%%%%%%%
43			fld=read_nctiles([myenv.nctilesdir 'ETAN'],'ETAN');
44			fld=std(fld,[],3); fld(find(fld==0))=NaN;
45			cc=[0:0.1:1]*0.10;
46
47			%%%%%%%%%%%%
48			%plot field:
49			%%%%%%%%%%%%
50			if ~myenv.lessplot;
51			figureL; gcmfaces_sphere(fld,cc,[],[],1);
52			end;
53
54