/[MITgcm]/MITgcm_contrib/gael/matlab_class/gcmfaces_demo.m
ViewVC logotype

Annotation of /MITgcm_contrib/gael/matlab_class/gcmfaces_demo.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph

Revision 1.19 - (hide annotations) (download)
Fri Mar 17 17:47:08 2017 UTC (8 years, 3 months ago) by gforget
Branch: MAIN
CVS Tags: checkpoint66e
Changes since 1.18: +9 -0 lines 
add directions to inputs needed for example_MITprof
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.17 %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 gforget 1.10 %to activate example_transports:
19     %
20 gforget 1.11 % mkdir release1
21     % wget --recursive ftp://mit.ecco-group.org/ecco_for_las/version_4/release1/nctiles_climatology
22 gforget 1.16 % mv mit.ecco-group.org/ecco_for_las/version_4/release1/nctiles_climatology release1/.
23 gforget 1.11 % rm -rf mit.ecco-group.org
24 gforget 1.10 %
25 gforget 1.19 %to activate example_MITprof:
26     %
27     % mkdir release2_climatology
28     % wget --recursive ftp://mit.ecco-group.org/ecco_for_las/version_4/release2/nctiles_climatology
29     % mv mit.ecco-group.org/ecco_for_las/version_4/release1/nctiles_climatology release2_climatology/.
30     % wget --recursive ftp://mit.ecco-group.org/ecco_for_las/version_4/release2/profiles
31     % mv mit.ecco-group.org/ecco_for_las/version_4/release1/profiles release2_climatology/.
32     % rm -rf mit.ecco-group.org
33     %
34 gforget 1.10 %to activate example_budget:
35     %
36 gforget 1.16 % mkdir sample_input
37 gforget 1.11 % wget --recursive ftp://mit.ecco-group.org/gforget/nctiles_budget_2d
38     % mv mit.ecco-group.org/gforget/nctiles_budget_2d sample_input/.
39     % rm -rf mit.ecco-group.org
40 gforget 1.10 %
41 gforget 1.14 % call sequence:
42     % addpath gcmfaces;
43     % addpath MITprof;
44     % gcmfaces_demo;
45 gforget 1.3
46 gforget 1.10 %choose verbose level
47 gforget 1.3 fprintf('\n Please set the amount of explanatory text display :\n');
48     fprintf(' 0: none.\n');
49     fprintf(' 1: comments.\n');
50     fprintf(' 2: comments preceeded with calling sequence.\n');
51     fprintf(' 3: same as 2, but preceeded with pause.\n');
52     verbose=input(' and/or type return. 0 is the default. \n');
53     if isempty(verbose); verbose=0; end;
54    
55 gforget 1.14 fprintf('\n');
56    
57     %so that gcmfaces_msg will work even before calling gcmfaces_global:
58     tmp1=which('gcmfaces_demo');
59     [PATH,NAME,EXT] = fileparts(tmp1);
60     addpath([PATH filesep 'gcmfaces_misc' filesep]);
61     global myenv;
62     myenv.issueWarnings=0;%skip warnings
63     myenv.verbose=verbose;%apply verbose level selected by user
64     %
65     fprintf('\n\n');
66     gcmfaces_msg('/////////////////////////////////////');
67     gcmfaces_msg('demo of gcmfaces_global and MITprof_global','// PART 0 :');
68     gcmfaces_msg('/////////////////////////////////////');
69 gforget 1.18 if myenv.verbose>0; gcmfaces_msg('please hit return','// >> '); pause; end;
70 gforget 1.14 fprintf('\n\n');
71 gforget 1.3 if myenv.verbose>0;
72 gforget 1.14 gcmfaces_msg(['* gcmfaces_global: adds gcmfaces directories to path' ...
73     ' and define environment variables (see myenv)']);
74     end;
75     myenv=[];
76     gcmfaces_global;%this will display warning
77     myenv.issueWarnings=0;%skip warnings
78     myenv.verbose=verbose;%apply verbose level selected by user
79     if myenv.verbose>0;
80     gcmfaces_msg('* (this warning gets resolved below by calling grid_load)');
81 gforget 1.3 end;
82 gforget 1.1
83 gforget 1.17 %if ~isempty(which('MITprof_global'));
84     % MITprof_global;
85     % if myenv.verbose>0;
86     % gcmfaces_msg('* MITprof_global: adds MITprof directories to path');
87     % end;
88     %end;
89 gforget 1.6
90 gforget 1.5 fprintf('\n\n');
91     gcmfaces_msg('/////////////////////////////////////');
92 gforget 1.14 gcmfaces_msg('demo of grid_load','// PART 1 :');
93     gcmfaces_msg('/////////////////////////////////////');
94 gforget 1.18 if myenv.verbose>0; gcmfaces_msg('please hit return','// >> '); pause; end;
95 gforget 1.14 fprintf('\n\n');
96     grid_load;
97    
98     fprintf('\n\n');
99     gcmfaces_msg('/////////////////////////////////////');
100 gforget 1.8 gcmfaces_msg('demo of plotting routines','// PART 1 :');
101 gforget 1.5 gcmfaces_msg('/////////////////////////////////////');
102 gforget 1.18 if myenv.verbose>0; gcmfaces_msg('please hit return','// >> '); pause; end;
103 gforget 1.5 fprintf('\n\n');
104 gforget 1.12 example_display;
105 gforget 1.1
106 gforget 1.5 fprintf('\n\n');
107     gcmfaces_msg('///////////////////////////////////////////');
108 gforget 1.15 gcmfaces_msg('demo of interpolation and remapping ','// PART 2 :');
109 gforget 1.5 gcmfaces_msg('///////////////////////////////////////////');
110 gforget 1.18 if myenv.verbose>0; gcmfaces_msg('please hit return','// >> '); pause; end;
111 gforget 1.5 fprintf('\n\n');
112 gforget 1.15 example_interp;
113 gforget 1.10
114 gforget 1.13 fprintf('\n\n');
115     gcmfaces_msg('///////////////////////////////////////////');
116 gforget 1.15 gcmfaces_msg('demo of bin averaging data sample to grid','// PART 2 :');
117 gforget 1.13 gcmfaces_msg('///////////////////////////////////////////');
118 gforget 1.18 if myenv.verbose>0; gcmfaces_msg('please hit return','// >> '); pause; end;
119 gforget 1.13 fprintf('\n\n');
120     fld=example_bin_average;
121    
122     fprintf('\n\n');
123     gcmfaces_msg('///////////////////////////////////////////');
124 gforget 1.15 gcmfaces_msg('demo of smoothing through diffusion','// PART 2 :');
125 gforget 1.13 gcmfaces_msg('///////////////////////////////////////////');
126 gforget 1.18 if myenv.verbose>0; gcmfaces_msg('please hit return','// >> '); pause; end;
127 gforget 1.13 fprintf('\n\n');
128     example_smooth(fld);
129 gforget 1.1
130 gforget 1.6 fprintf('\n\n');
131     gcmfaces_msg('/////////////////////////////////////////');
132 gforget 1.8 gcmfaces_msg('demo of transport computations','// PART 3 :');
133     gcmfaces_msg('/////////////////////////////////////////');
134 gforget 1.18 if myenv.verbose>0; gcmfaces_msg('please hit return','// >> '); pause; end;
135 gforget 1.8 fprintf('\n\n');
136 gforget 1.9 diags=example_transports;
137 gforget 1.13
138 gforget 1.8 if ~isempty(diags);
139 gforget 1.13 fprintf('\n\n');
140     gcmfaces_msg('/////////////////////////////////////////');
141     gcmfaces_msg('demo of transport display','// PART 3 :');
142     gcmfaces_msg('/////////////////////////////////////////');
143 gforget 1.18 if myenv.verbose>0; gcmfaces_msg('please hit return','// >> '); pause; end;
144 gforget 1.13 fprintf('\n\n');
145 gforget 1.8 example_transports_disp(diags);
146     end;
147    
148     fprintf('\n\n');
149     gcmfaces_msg('/////////////////////////////////////////');
150 gforget 1.6 gcmfaces_msg('demo of budget computations','// PART 4 :');
151     gcmfaces_msg('/////////////////////////////////////////');
152 gforget 1.18 if myenv.verbose>0; gcmfaces_msg('please hit return','// >> '); pause; end;
153 gforget 1.6 fprintf('\n\n');
154 gforget 1.9 example_budget;
155 gforget 1.1
  ViewVC Help
Powered by ViewVC 1.1.22