gael/matlab_class/gcmfaces_demo.m


%call sequence:
%
%  addpath gcmfaces; gcmfaces_demo;
%
%expected directory structure:
%
%  gcmfaces     (codes)
%  MITprof      (codes)
%  nctiles_grid (ECCO v4 grid)
%  release1     (ECCO v4 output)
%  sample_input (test material)
%
%to activate example_remap:
%
%  mkdir sample_input
%  wget ftp://mit.ecco-group.org/gforget/testcase_remap.mat
%  mv testcase_remap.mat sample_input/.
%
%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/gforget/nctiles_climatology release1/.
%  rm -rf mit.ecco-group.org
%
%to activate example_budget:
%
%  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
%

%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;

%initialize environment variables and mygrid
addpath gcmfaces;
gcmfaces_global;
myenv.issueWarnings=0;
myenv.verbose=verbose;
if myenv.verbose>0;
    gcmfaces_msg('* set path and environment variables (myenv) by calling gcmfaces_global');
end;

addpath MITprof;
if ~isempty(which('MITprof_global'));
    MITprof_global;
    if myenv.verbose>0;
        gcmfaces_msg('* set MITprof path by calling MITprof_global');
    end;
end;

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

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

if 1;%one method to map lon-lat fields to gcmfaces
    example_remap(1);
else;%alternatives
    example_griddata;
    example_interp;
    example_faces2latlon2faces;
end;
example_bin_average(1);%incl. call to example_smooth.m;

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);
    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.1
2	gforget	1.11	%call sequence:
3	gforget	1.9	%
4	gforget	1.11	% addpath gcmfaces; gcmfaces_demo;
5			%
6			%expected directory structure:
7			%
8			% gcmfaces (codes)
9			% MITprof (codes)
10			% nctiles_grid (ECCO v4 grid)
11			% release1 (ECCO v4 output)
12			% sample_input (test material)
13			%
14			%to activate example_remap:
15			%
16			% mkdir sample_input
17			% wget ftp://mit.ecco-group.org/gforget/testcase_remap.mat
18			% mv testcase_remap.mat sample_input/.
19	gforget	1.10	%
20			%to activate example_transports:
21			%
22	gforget	1.11	% mkdir release1
23			% wget --recursive ftp://mit.ecco-group.org/ecco_for_las/version_4/release1/nctiles_climatology
24			% mv mit.ecco-group.org/gforget/nctiles_climatology release1/.
25			% rm -rf mit.ecco-group.org
26	gforget	1.10	%
27			%to activate example_budget:
28			%
29	gforget	1.11	% wget --recursive ftp://mit.ecco-group.org/gforget/nctiles_budget_2d
30			% mv mit.ecco-group.org/gforget/nctiles_budget_2d sample_input/.
31			% rm -rf mit.ecco-group.org
32	gforget	1.10	%
33	gforget	1.3
34	gforget	1.10	%choose verbose level
35	gforget	1.3	fprintf('\n Please set the amount of explanatory text display :\n');
36			fprintf(' 0: none.\n');
37			fprintf(' 1: comments.\n');
38			fprintf(' 2: comments preceeded with calling sequence.\n');
39			fprintf(' 3: same as 2, but preceeded with pause.\n');
40			verbose=input(' and/or type return. 0 is the default. \n');
41			if isempty(verbose); verbose=0; end;
42
43	gforget	1.2	%initialize environment variables and mygrid
44	gforget	1.6	addpath gcmfaces;
45	gforget	1.10	gcmfaces_global;
46	gforget	1.8	myenv.issueWarnings=0;
47	gforget	1.3	myenv.verbose=verbose;
48			if myenv.verbose>0;
49			gcmfaces_msg('* set path and environment variables (myenv) by calling gcmfaces_global');
50			end;
51	gforget	1.1
52	gforget	1.7	addpath MITprof;
53	gforget	1.6	if ~isempty(which('MITprof_global'));
54	gforget	1.10	MITprof_global;
55			if myenv.verbose>0;
56			gcmfaces_msg('* set MITprof path by calling MITprof_global');
57			end;
58	gforget	1.6	end;
59
60	gforget	1.5	fprintf('\n\n');
61			gcmfaces_msg('/////////////////////////////////////');
62	gforget	1.8	gcmfaces_msg('demo of plotting routines','// PART 1 :');
63	gforget	1.5	gcmfaces_msg('please hit return','// >> ');
64			gcmfaces_msg('/////////////////////////////////////');
65			fprintf('\n\n');
66			pause;
67	gforget	1.9	plot_one_field;
68	gforget	1.1
69	gforget	1.5	fprintf('\n\n');
70			gcmfaces_msg('///////////////////////////////////////////');
71	gforget	1.10	gcmfaces_msg('demo of remapping and smoothing','// PART 2 :');
72	gforget	1.5	gcmfaces_msg('please hit return','// >> ');
73			gcmfaces_msg('///////////////////////////////////////////');
74			fprintf('\n\n');
75			pause;
76	gforget	1.10
77			if 1;%one method to map lon-lat fields to gcmfaces
78			example_remap(1);
79			else;%alternatives
80			example_griddata;
81			example_interp;
82			example_faces2latlon2faces;
83			end;
84	gforget	1.9	example_bin_average(1);%incl. call to example_smooth.m;
85	gforget	1.1
86	gforget	1.6	fprintf('\n\n');
87			gcmfaces_msg('/////////////////////////////////////////');
88	gforget	1.8	gcmfaces_msg('demo of transport computations','// PART 3 :');
89			gcmfaces_msg('please hit return','// >> ');
90			gcmfaces_msg('/////////////////////////////////////////');
91			fprintf('\n\n');
92			pause;
93	gforget	1.9	diags=example_transports;
94	gforget	1.8	if ~isempty(diags);
95			example_transports_disp(diags);
96			end;
97
98			fprintf('\n\n');
99			gcmfaces_msg('/////////////////////////////////////////');
100	gforget	1.6	gcmfaces_msg('demo of budget computations','// PART 4 :');
101			gcmfaces_msg('please hit return','// >> ');
102			gcmfaces_msg('/////////////////////////////////////////');
103			fprintf('\n\n');
104			pause;
105	gforget	1.9	example_budget;
106	gforget	1.1