/[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.6 by gforget, Thu Feb 12 21:29:30 2015 UTC revision 1.13 by gforget, Sun Jan 24 15:12:33 2016 UTC
# Line 1  Line 1 
1    
2  choiceGrid='v4';  %call sequence:
3    %
4  %need to replace the missing v3 files...  %  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    %
20    %to activate example_transports:
21    %
22    %  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    %
27    %to activate example_budget:
28    %
29    %  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  %  %
 %choose grid for test  
 % fprintf('\n Please set the demo grid : \n');  
 % fprintf('    ''llc'': lat-lon-cap grid (5 faces). \n');  
 % fprintf('    ''ll'' : simple lat-lon grid (1 face). \n');  
 % choiceGrid=input(' and type return. ''llc'' is the default.\n');  
 % if isempty(choiceGrid); choiceGrid='llc'; end;  
 % if strcmp(choiceGrid,'llc'); choiceGrid='v4';  
 % elseif strcmp(choiceGrid,'ll'); choiceGrid='v3';  
 % else; error('wrong grid choice');  
 % end;  
33    
34  %choose verbise level  %choose verbose level
35  fprintf('\n Please set the amount of explanatory text display :\n');  fprintf('\n Please set the amount of explanatory text display :\n');
36  fprintf('    0: none.\n');  fprintf('    0: none.\n');
37  fprintf('    1: comments.\n');  fprintf('    1: comments.\n');
# Line 26  if isempty(verbose); verbose=0; end; Line 43  if isempty(verbose); verbose=0; end;
43  %initialize environment variables and mygrid  %initialize environment variables and mygrid
44  addpath gcmfaces;  addpath gcmfaces;
45  gcmfaces_global;  gcmfaces_global;
46    myenv.issueWarnings=0;
47  myenv.verbose=verbose;  myenv.verbose=verbose;
48  if myenv.verbose>0;  if myenv.verbose>0;
49      gcmfaces_msg('* set path and environment variables (myenv) by calling gcmfaces_global');      gcmfaces_msg('* set path and environment variables (myenv) by calling gcmfaces_global');
50  end;  end;
51    
52    addpath MITprof;
53  if ~isempty(which('MITprof_global'));  if ~isempty(which('MITprof_global'));
54  MITprof_global;      MITprof_global;
55  if myenv.verbose>0;      if myenv.verbose>0;
56      gcmfaces_msg('* set MITprof path by calling MITprof_global');          gcmfaces_msg('* set MITprof path by calling MITprof_global');
57  end;      end;
58  end;  end;
59    
60  fprintf('\n\n');  fprintf('\n\n');
61  gcmfaces_msg('/////////////////////////////////////////');  gcmfaces_msg('/////////////////////////////////////');
62  gcmfaces_msg('demo of transport computations','// PART 1 :');  gcmfaces_msg('demo of plotting routines','// PART 1 :');
63  gcmfaces_msg('please hit return','// >> ');  gcmfaces_msg('please hit return','// >> ');
64  gcmfaces_msg('/////////////////////////////////////////');  gcmfaces_msg('/////////////////////////////////////');
65  fprintf('\n\n');  fprintf('\n\n');
66  pause;  pause;
67  diags=example_transports(choiceGrid);  example_display;
 if ~isempty(diags);  
     example_transports_disp(diags);  
 end;  
68    
69  fprintf('\n\n');  fprintf('\n\n');
70  gcmfaces_msg('/////////////////////////////////////');  gcmfaces_msg('///////////////////////////////////////////');
71  gcmfaces_msg('demo of plotting routines','// PART 2 :');  gcmfaces_msg('demo of remapping ','// PART 2 :');
72  gcmfaces_msg('please hit return','// >> ');  gcmfaces_msg('please hit return','// >> ');
73  gcmfaces_msg('/////////////////////////////////////');  gcmfaces_msg('///////////////////////////////////////////');
74    fprintf('\n\n');
75    pause;
76    example_remap(0);
77    
78    fprintf('\n\n');
79    gcmfaces_msg('///////////////////////////////////////////');
80    gcmfaces_msg('demo of bin averaging (gridding)','// PART 2 :');
81    gcmfaces_msg('please hit return','// >> ');
82    gcmfaces_msg('///////////////////////////////////////////');
83  fprintf('\n\n');  fprintf('\n\n');
84  pause;  pause;
85  plot_one_field(choiceGrid,0);%%incl. call to m_map_gcmfaces;  %alternatives
86    %example_griddata;
87    %example_interp;
88    %example_faces2latlon2faces;
89    fld=example_bin_average;
90    
91  fprintf('\n\n');  fprintf('\n\n');
92  gcmfaces_msg('///////////////////////////////////////////');  gcmfaces_msg('///////////////////////////////////////////');
93  gcmfaces_msg('demo of gridding v interpolating','// PART 3 :');  gcmfaces_msg('demo of smoothing ','// PART 2 :');
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  example_bin_average(choiceGrid,1);%incl. call to example_smooth.m;  example_smooth(fld);
99  example_griddata(choiceGrid);  
100  example_interp(choiceGrid);  fprintf('\n\n');
101  example_faces2latlon2faces(choiceGrid);  gcmfaces_msg('/////////////////////////////////////////');
102    gcmfaces_msg('demo of transport computations','// PART 3 :');
103    gcmfaces_msg('please hit return','// >> ');
104    gcmfaces_msg('/////////////////////////////////////////');
105    fprintf('\n\n');
106    pause;
107    diags=example_transports;
108    
109    if ~isempty(diags);
110        fprintf('\n\n');
111        gcmfaces_msg('/////////////////////////////////////////');
112        gcmfaces_msg('demo of transport display','// PART 3 :');
113        gcmfaces_msg('please hit return','// >> ');
114        gcmfaces_msg('/////////////////////////////////////////');
115        fprintf('\n\n');
116        example_transports_disp(diags);
117    end;
118    
119  fprintf('\n\n');  fprintf('\n\n');
120  gcmfaces_msg('/////////////////////////////////////////');  gcmfaces_msg('/////////////////////////////////////////');
# Line 78  gcmfaces_msg('please hit return','// >> Line 123  gcmfaces_msg('please hit return','// >>
123  gcmfaces_msg('/////////////////////////////////////////');  gcmfaces_msg('/////////////////////////////////////////');
124  fprintf('\n\n');  fprintf('\n\n');
125  pause;  pause;
126  example_budget(choiceGrid);  example_budget;
127    
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22