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

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

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

revision 1.7 by gforget, Thu Feb 12 21:40:26 2015 UTC revision 1.17 by gforget, Mon Mar 14 01:24:17 2016 UTC
# Line 1  Line 1 
1    function []=gcmfaces_demo();
2  choiceGrid='v4';  % GCMFACES_DEMO demonstrate capabilities of the gcmfaces toolbox.
3    %
4  %need to replace the missing v3 files...  %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  %choose grid for test  %  mkdir release1
21  % fprintf('\n Please set the demo grid : \n');  %  wget --recursive ftp://mit.ecco-group.org/ecco_for_las/version_4/release1/nctiles_climatology
22  % fprintf('    ''llc'': lat-lon-cap grid (5 faces). \n');  %  mv mit.ecco-group.org/ecco_for_las/version_4/release1/nctiles_climatology release1/.
23  % fprintf('    ''ll'' : simple lat-lon grid (1 face). \n');  %  rm -rf mit.ecco-group.org
24  % choiceGrid=input(' and type return. ''llc'' is the default.\n');  %
25  % if isempty(choiceGrid); choiceGrid='llc'; end;  %to activate example_budget:
26  % if strcmp(choiceGrid,'llc'); choiceGrid='v4';  %
27  % elseif strcmp(choiceGrid,'ll'); choiceGrid='v3';  %  mkdir sample_input
28  % else; error('wrong grid choice');  %  wget --recursive ftp://mit.ecco-group.org/gforget/nctiles_budget_2d
29  % end;  %  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 verbise level  %choose verbose level
38  fprintf('\n Please set the amount of explanatory text display :\n');  fprintf('\n Please set the amount of explanatory text display :\n');
39  fprintf('    0: none.\n');  fprintf('    0: none.\n');
40  fprintf('    1: comments.\n');  fprintf('    1: comments.\n');
# Line 23  fprintf('    3: same as 2, but preceeded Line 43  fprintf('    3: same as 2, but preceeded
43  verbose=input(' and/or type return. 0 is the default. \n');  verbose=input(' and/or type return. 0 is the default. \n');
44  if isempty(verbose); verbose=0; end;  if isempty(verbose); verbose=0; end;
45    
46  %initialize environment variables and mygrid  fprintf('\n');
 addpath gcmfaces;  
 gcmfaces_global;  
 myenv.verbose=verbose;  
 if myenv.verbose>0;  
     gcmfaces_msg('* set path and environment variables (myenv) by calling gcmfaces_global');  
 end;  
47    
48  addpath MITprof;  %so that gcmfaces_msg will work even before calling gcmfaces_global:
49  if ~isempty(which('MITprof_global'));  tmp1=which('gcmfaces_demo');
50  MITprof_global;  [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;  if myenv.verbose>0;
64      gcmfaces_msg('* set MITprof path by calling MITprof_global');      gcmfaces_msg(['* gcmfaces_global: adds gcmfaces directories to path' ...
65                     ' and define environment variables (see myenv)']);
66  end;  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;  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');  fprintf('\n\n');
83  gcmfaces_msg('/////////////////////////////////////////');  gcmfaces_msg('/////////////////////////////////////');
84  gcmfaces_msg('demo of transport computations','// PART 1 :');  gcmfaces_msg('demo of grid_load','// PART 1 :');
85  gcmfaces_msg('please hit return','// >> ');  gcmfaces_msg('please hit return','// >> ');
86  gcmfaces_msg('/////////////////////////////////////////');  gcmfaces_msg('/////////////////////////////////////');
87  fprintf('\n\n');  fprintf('\n\n');
88  pause;  pause;
89  diags=example_transports(choiceGrid);  grid_load;
 if ~isempty(diags);  
     example_transports_disp(diags);  
 end;  
90    
91  fprintf('\n\n');  fprintf('\n\n');
92  gcmfaces_msg('/////////////////////////////////////');  gcmfaces_msg('/////////////////////////////////////');
93  gcmfaces_msg('demo of plotting routines','// PART 2 :');  gcmfaces_msg('demo of plotting routines','// PART 1 :');
94  gcmfaces_msg('please hit return','// >> ');  gcmfaces_msg('please hit return','// >> ');
95  gcmfaces_msg('/////////////////////////////////////');  gcmfaces_msg('/////////////////////////////////////');
96  fprintf('\n\n');  fprintf('\n\n');
97  pause;  pause;
98  plot_one_field(choiceGrid,0);%%incl. call to m_map_gcmfaces;  example_display;
99    
100  fprintf('\n\n');  fprintf('\n\n');
101  gcmfaces_msg('///////////////////////////////////////////');  gcmfaces_msg('///////////////////////////////////////////');
102  gcmfaces_msg('demo of gridding v interpolating','// PART 3 :');  gcmfaces_msg('demo of interpolation and remapping ','// PART 2 :');
103  gcmfaces_msg('please hit return','// >> ');  gcmfaces_msg('please hit return','// >> ');
104  gcmfaces_msg('///////////////////////////////////////////');  gcmfaces_msg('///////////////////////////////////////////');
105  fprintf('\n\n');  fprintf('\n\n');
106  pause;  pause;
107  example_bin_average(choiceGrid,1);%incl. call to example_smooth.m;  example_interp;
108  example_griddata(choiceGrid);  
109  example_interp(choiceGrid);  fprintf('\n\n');
110  example_faces2latlon2faces(choiceGrid);  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');  fprintf('\n\n');
147  gcmfaces_msg('/////////////////////////////////////////');  gcmfaces_msg('/////////////////////////////////////////');
# Line 79  gcmfaces_msg('please hit return','// >> Line 150  gcmfaces_msg('please hit return','// >>
150  gcmfaces_msg('/////////////////////////////////////////');  gcmfaces_msg('/////////////////////////////////////////');
151  fprintf('\n\n');  fprintf('\n\n');
152  pause;  pause;
153  example_budget(choiceGrid);  example_budget;
154    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22