gael/matlab_class/gcmfaces_demo.m

function []=gcmfaces_demo();
% GCMFACES_DEMO demonstrate capabilities of the gcmfaces toolbox.
%
%expected directory structure:
%
%  gcmfaces     (codes)
%  MITprof      (codes)
%  nctiles_grid (ECCO v4 grid)
%  release1     (ECCO v4 output)
%  sample_input (additional demo material)
%
%to activate example_transports:
%
%  mkdir release1
%  wget --recursive ftp://mit.ecco-group.org/ecco_for_las/version_4/release1/nctiles_climatology
%  mv mit.ecco-group.org/ecco_for_las/version_4/release1/nctiles_climatology release1/.
%  rm -rf mit.ecco-group.org
%
%to activate example_budget:
%
%  mkdir sample_input
%  wget --recursive ftp://mit.ecco-group.org/gforget/nctiles_budget_2d
%  mv mit.ecco-group.org/gforget/nctiles_budget_2d sample_input/.
%  rm -rf mit.ecco-group.org
%
% call sequence:
%  addpath gcmfaces;
%  addpath MITprof;
%  gcmfaces_demo;

%choose verbose level
fprintf('\n Please set the amount of explanatory text display :\n');
fprintf('    0: none.\n');
fprintf('    1: comments.\n');
fprintf('    2: comments preceeded with calling sequence.\n');
fprintf('    3: same as 2, but preceeded with pause.\n');
verbose=input(' and/or type return. 0 is the default. \n');
if isempty(verbose); verbose=0; end;

fprintf('\n');

%so that gcmfaces_msg will work even before calling gcmfaces_global:
tmp1=which('gcmfaces_demo');
[PATH,NAME,EXT] = fileparts(tmp1);
addpath([PATH filesep 'gcmfaces_misc' filesep]); 
global myenv; 
myenv.issueWarnings=0;%skip warnings
myenv.verbose=verbose;%apply verbose level selected by user
%
fprintf('\n\n');
gcmfaces_msg('/////////////////////////////////////');
gcmfaces_msg('demo of gcmfaces_global and MITprof_global','// PART 0 :');
gcmfaces_msg('please hit return','// >> ');
gcmfaces_msg('/////////////////////////////////////');
fprintf('\n\n');
pause;
if myenv.verbose>0;
    gcmfaces_msg(['* gcmfaces_global: adds gcmfaces directories to path' ...
                 ' and define environment variables (see myenv)']);
end;
myenv=[];
gcmfaces_global;%this will display warning
myenv.issueWarnings=0;%skip warnings
myenv.verbose=verbose;%apply verbose level selected by user
if myenv.verbose>0;
    gcmfaces_msg('* (this warning gets resolved below by calling grid_load)');
end;

if ~isempty(which('MITprof_global'));
    MITprof_global;
    if myenv.verbose>0;
        gcmfaces_msg('* MITprof_global: adds MITprof directories to path');
    end;
end;

fprintf('\n\n');
gcmfaces_msg('/////////////////////////////////////');
gcmfaces_msg('demo of grid_load','// PART 1 :');
gcmfaces_msg('please hit return','// >> ');
gcmfaces_msg('/////////////////////////////////////');
fprintf('\n\n');
pause;
grid_load;

fprintf('\n\n');
gcmfaces_msg('/////////////////////////////////////');
gcmfaces_msg('demo of plotting routines','// PART 1 :');
gcmfaces_msg('please hit return','// >> ');
gcmfaces_msg('/////////////////////////////////////');
fprintf('\n\n');
pause;
example_display;

fprintf('\n\n');
gcmfaces_msg('///////////////////////////////////////////');
gcmfaces_msg('demo of interpolation and remapping ','// PART 2 :');
gcmfaces_msg('please hit return','// >> ');
gcmfaces_msg('///////////////////////////////////////////');
fprintf('\n\n');
pause;
example_interp;

fprintf('\n\n');
gcmfaces_msg('///////////////////////////////////////////');
gcmfaces_msg('demo of bin averaging data sample to grid','// PART 2 :');
gcmfaces_msg('please hit return','// >> ');
gcmfaces_msg('///////////////////////////////////////////');
fprintf('\n\n');
pause;
fld=example_bin_average;

fprintf('\n\n');
gcmfaces_msg('///////////////////////////////////////////');
gcmfaces_msg('demo of smoothing through diffusion','// PART 2 :');
gcmfaces_msg('please hit return','// >> ');
gcmfaces_msg('///////////////////////////////////////////');
fprintf('\n\n');
pause;
example_smooth(fld);

fprintf('\n\n');
gcmfaces_msg('/////////////////////////////////////////');
gcmfaces_msg('demo of transport computations','// PART 3 :');
gcmfaces_msg('please hit return','// >> ');
gcmfaces_msg('/////////////////////////////////////////');
fprintf('\n\n');
pause;
diags=example_transports;

if ~isempty(diags);
    fprintf('\n\n');
    gcmfaces_msg('/////////////////////////////////////////');
    gcmfaces_msg('demo of transport display','// PART 3 :');
    gcmfaces_msg('please hit return','// >> ');
    gcmfaces_msg('/////////////////////////////////////////');
    fprintf('\n\n');
    example_transports_disp(diags);
end;

fprintf('\n\n');
gcmfaces_msg('/////////////////////////////////////////');
gcmfaces_msg('demo of budget computations','// PART 4 :');
gcmfaces_msg('please hit return','// >> ');
gcmfaces_msg('/////////////////////////////////////////');
fprintf('\n\n');
pause;
example_budget;

1	gforget	1.14	function []=gcmfaces_demo();
2			% GCMFACES_DEMO demonstrate capabilities of the gcmfaces toolbox.
3	gforget	1.11	%
4			%expected directory structure:
5			%
6			% gcmfaces (codes)
7			% MITprof (codes)
8			% nctiles_grid (ECCO v4 grid)
9			% release1 (ECCO v4 output)
10	gforget	1.14	% sample_input (additional demo material)
11	gforget	1.11	%
12	gforget	1.10	%to activate example_transports:
13			%
14	gforget	1.11	% mkdir release1
15			% wget --recursive ftp://mit.ecco-group.org/ecco_for_las/version_4/release1/nctiles_climatology
16	gforget	1.16	% mv mit.ecco-group.org/ecco_for_las/version_4/release1/nctiles_climatology release1/.
17	gforget	1.11	% rm -rf mit.ecco-group.org
18	gforget	1.10	%
19			%to activate example_budget:
20			%
21	gforget	1.16	% mkdir sample_input
22	gforget	1.11	% wget --recursive ftp://mit.ecco-group.org/gforget/nctiles_budget_2d
23			% mv mit.ecco-group.org/gforget/nctiles_budget_2d sample_input/.
24			% rm -rf mit.ecco-group.org
25	gforget	1.10	%
26	gforget	1.14	% call sequence:
27			% addpath gcmfaces;
28			% addpath MITprof;
29			% gcmfaces_demo;
30	gforget	1.3
31	gforget	1.10	%choose verbose level
32	gforget	1.3	fprintf('\n Please set the amount of explanatory text display :\n');
33			fprintf(' 0: none.\n');
34			fprintf(' 1: comments.\n');
35			fprintf(' 2: comments preceeded with calling sequence.\n');
36			fprintf(' 3: same as 2, but preceeded with pause.\n');
37			verbose=input(' and/or type return. 0 is the default. \n');
38			if isempty(verbose); verbose=0; end;
39
40	gforget	1.14	fprintf('\n');
41
42			%so that gcmfaces_msg will work even before calling gcmfaces_global:
43			tmp1=which('gcmfaces_demo');
44			[PATH,NAME,EXT] = fileparts(tmp1);
45			addpath([PATH filesep 'gcmfaces_misc' filesep]);
46			global myenv;
47			myenv.issueWarnings=0;%skip warnings
48			myenv.verbose=verbose;%apply verbose level selected by user
49			%
50			fprintf('\n\n');
51			gcmfaces_msg('/////////////////////////////////////');
52			gcmfaces_msg('demo of gcmfaces_global and MITprof_global','// PART 0 :');
53			gcmfaces_msg('please hit return','// >> ');
54			gcmfaces_msg('/////////////////////////////////////');
55			fprintf('\n\n');
56			pause;
57	gforget	1.3	if myenv.verbose>0;
58	gforget	1.14	gcmfaces_msg(['* gcmfaces_global: adds gcmfaces directories to path' ...
59			' and define environment variables (see myenv)']);
60			end;
61			myenv=[];
62			gcmfaces_global;%this will display warning
63			myenv.issueWarnings=0;%skip warnings
64			myenv.verbose=verbose;%apply verbose level selected by user
65			if myenv.verbose>0;
66			gcmfaces_msg('* (this warning gets resolved below by calling grid_load)');
67	gforget	1.3	end;
68	gforget	1.1
69	gforget	1.6	if ~isempty(which('MITprof_global'));
70	gforget	1.10	MITprof_global;
71			if myenv.verbose>0;
72	gforget	1.14	gcmfaces_msg('* MITprof_global: adds MITprof directories to path');
73	gforget	1.10	end;
74	gforget	1.6	end;
75
76	gforget	1.5	fprintf('\n\n');
77			gcmfaces_msg('/////////////////////////////////////');
78	gforget	1.14	gcmfaces_msg('demo of grid_load','// PART 1 :');
79			gcmfaces_msg('please hit return','// >> ');
80			gcmfaces_msg('/////////////////////////////////////');
81			fprintf('\n\n');
82			pause;
83			grid_load;
84
85			fprintf('\n\n');
86			gcmfaces_msg('/////////////////////////////////////');
87	gforget	1.8	gcmfaces_msg('demo of plotting routines','// PART 1 :');
88	gforget	1.5	gcmfaces_msg('please hit return','// >> ');
89			gcmfaces_msg('/////////////////////////////////////');
90			fprintf('\n\n');
91			pause;
92	gforget	1.12	example_display;
93	gforget	1.1
94	gforget	1.5	fprintf('\n\n');
95			gcmfaces_msg('///////////////////////////////////////////');
96	gforget	1.15	gcmfaces_msg('demo of interpolation and remapping ','// PART 2 :');
97	gforget	1.5	gcmfaces_msg('please hit return','// >> ');
98			gcmfaces_msg('///////////////////////////////////////////');
99			fprintf('\n\n');
100			pause;
101	gforget	1.15	example_interp;
102	gforget	1.10
103	gforget	1.13	fprintf('\n\n');
104			gcmfaces_msg('///////////////////////////////////////////');
105	gforget	1.15	gcmfaces_msg('demo of bin averaging data sample to grid','// PART 2 :');
106	gforget	1.13	gcmfaces_msg('please hit return','// >> ');
107			gcmfaces_msg('///////////////////////////////////////////');
108			fprintf('\n\n');
109			pause;
110			fld=example_bin_average;
111
112			fprintf('\n\n');
113			gcmfaces_msg('///////////////////////////////////////////');
114	gforget	1.15	gcmfaces_msg('demo of smoothing through diffusion','// PART 2 :');
115	gforget	1.13	gcmfaces_msg('please hit return','// >> ');
116			gcmfaces_msg('///////////////////////////////////////////');
117			fprintf('\n\n');
118			pause;
119			example_smooth(fld);
120	gforget	1.1
121	gforget	1.6	fprintf('\n\n');
122			gcmfaces_msg('/////////////////////////////////////////');
123	gforget	1.8	gcmfaces_msg('demo of transport computations','// PART 3 :');
124			gcmfaces_msg('please hit return','// >> ');
125			gcmfaces_msg('/////////////////////////////////////////');
126			fprintf('\n\n');
127			pause;
128	gforget	1.9	diags=example_transports;
129	gforget	1.13
130	gforget	1.8	if ~isempty(diags);
131	gforget	1.13	fprintf('\n\n');
132			gcmfaces_msg('/////////////////////////////////////////');
133			gcmfaces_msg('demo of transport display','// PART 3 :');
134			gcmfaces_msg('please hit return','// >> ');
135			gcmfaces_msg('/////////////////////////////////////////');
136			fprintf('\n\n');
137	gforget	1.8	example_transports_disp(diags);
138			end;
139
140			fprintf('\n\n');
141			gcmfaces_msg('/////////////////////////////////////////');
142	gforget	1.6	gcmfaces_msg('demo of budget computations','// PART 4 :');
143			gcmfaces_msg('please hit return','// >> ');
144			gcmfaces_msg('/////////////////////////////////////////');
145			fprintf('\n\n');
146			pause;
147	gforget	1.9	example_budget;
148	gforget	1.1