/[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.1 by gforget, Tue Jun 21 21:07:09 2011 UTC revision 1.19 by gforget, Fri Mar 17 17:47:08 2017 UTC
# Line 1  Line 1 
1    function []=gcmfaces_demo();
2  %environment variables:  % 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_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    %to activate example_budget:
35    %
36    %  mkdir sample_input
37    %  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    %
41    % call sequence:
42    %  addpath gcmfaces;
43    %  addpath MITprof;
44    %  gcmfaces_demo;
45    
46    %choose verbose level
47    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    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;  global myenv;
62  myenv.gcmfaces_verbose=0;  myenv.issueWarnings=0;%skip warnings
63  myenv.gcmfaces_dir=[pwd '/gcmfaces/'];  myenv.verbose=verbose;%apply verbose level selected by user
64  myenv.MITprof_verbose=0;  %
65  myenv.MITprof_dir=[pwd '/gcmfaces/'];  fprintf('\n\n');
66  myenv.useNativeMatlabNetcdf = ~isempty(which('netcdf.open'));  gcmfaces_msg('/////////////////////////////////////');
67    gcmfaces_msg('demo of gcmfaces_global and MITprof_global','// PART 0 :');
68  %re-init mygrid:  gcmfaces_msg('/////////////////////////////////////');
69  global mygrid; mygrid=[];  if myenv.verbose>0; gcmfaces_msg('please hit return','// >> '); pause; end;
70    fprintf('\n\n');
71  choiceGrid='v4';  if myenv.verbose>0;
72        gcmfaces_msg(['* gcmfaces_global: adds gcmfaces directories to path' ...
73  %transform those to functions with v3/v4 arguments:                   ' and define environment variables (see myenv)']);
74  basic_diags_compute_v3_or_v4(choiceGrid);  end;
75  basic_diags_display_v3_or_v4(choiceGrid);  myenv=[];
76    gcmfaces_global;%this will display warning
77  example_bin_average(choiceGrid,1);%incl. call to example_smooth.m;  myenv.issueWarnings=0;%skip warnings
78  example_griddata(choiceGrid);  myenv.verbose=verbose;%apply verbose level selected by user
79  example_interp(choiceGrid);  if myenv.verbose>0;
80        gcmfaces_msg('* (this warning gets resolved below by calling grid_load)');
81  plot_one_field(choiceGrid,0);%%incl. call to m_map_gcmfaces;  end;
82    
83    %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    
90    fprintf('\n\n');
91    gcmfaces_msg('/////////////////////////////////////');
92    gcmfaces_msg('demo of grid_load','// PART 1 :');
93    gcmfaces_msg('/////////////////////////////////////');
94    if myenv.verbose>0; gcmfaces_msg('please hit return','// >> '); pause; end;
95    fprintf('\n\n');
96    grid_load;
97    
98    fprintf('\n\n');
99    gcmfaces_msg('/////////////////////////////////////');
100    gcmfaces_msg('demo of plotting routines','// PART 1 :');
101    gcmfaces_msg('/////////////////////////////////////');
102    if myenv.verbose>0; gcmfaces_msg('please hit return','// >> '); pause; end;
103    fprintf('\n\n');
104    example_display;
105    
106    fprintf('\n\n');
107    gcmfaces_msg('///////////////////////////////////////////');
108    gcmfaces_msg('demo of interpolation and remapping ','// PART 2 :');
109    gcmfaces_msg('///////////////////////////////////////////');
110    if myenv.verbose>0; gcmfaces_msg('please hit return','// >> '); pause; end;
111    fprintf('\n\n');
112    example_interp;
113    
114    fprintf('\n\n');
115    gcmfaces_msg('///////////////////////////////////////////');
116    gcmfaces_msg('demo of bin averaging data sample to grid','// PART 2 :');
117    gcmfaces_msg('///////////////////////////////////////////');
118    if myenv.verbose>0; gcmfaces_msg('please hit return','// >> '); pause; end;
119    fprintf('\n\n');
120    fld=example_bin_average;
121    
122    fprintf('\n\n');
123    gcmfaces_msg('///////////////////////////////////////////');
124    gcmfaces_msg('demo of smoothing through diffusion','// PART 2 :');
125    gcmfaces_msg('///////////////////////////////////////////');
126    if myenv.verbose>0; gcmfaces_msg('please hit return','// >> '); pause; end;
127    fprintf('\n\n');
128    example_smooth(fld);
129    
130    fprintf('\n\n');
131    gcmfaces_msg('/////////////////////////////////////////');
132    gcmfaces_msg('demo of transport computations','// PART 3 :');
133    gcmfaces_msg('/////////////////////////////////////////');
134    if myenv.verbose>0; gcmfaces_msg('please hit return','// >> '); pause; end;
135    fprintf('\n\n');
136    diags=example_transports;
137    
138    if ~isempty(diags);
139        fprintf('\n\n');
140        gcmfaces_msg('/////////////////////////////////////////');
141        gcmfaces_msg('demo of transport display','// PART 3 :');
142        gcmfaces_msg('/////////////////////////////////////////');
143        if myenv.verbose>0; gcmfaces_msg('please hit return','// >> '); pause; end;
144        fprintf('\n\n');
145        example_transports_disp(diags);
146    end;
147    
148    fprintf('\n\n');
149    gcmfaces_msg('/////////////////////////////////////////');
150    gcmfaces_msg('demo of budget computations','// PART 4 :');
151    gcmfaces_msg('/////////////////////////////////////////');
152    if myenv.verbose>0; gcmfaces_msg('please hit return','// >> '); pause; end;
153    fprintf('\n\n');
154    example_budget;
155    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22