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