/[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.4 by gforget, Tue Jul 31 19:12:44 2012 UTC revision 1.11 by gforget, Mon Jan 11 17:28:10 2016 UTC
# Line 1  Line 1 
1    
2  %choose grid for test  %call sequence:
3  fprintf('\n Please set the demo grid : \n');  %
4  fprintf('    ''llc'': lat-lon-cap grid (5 faces). \n');  %  addpath gcmfaces; gcmfaces_demo;
5  fprintf('    ''ll'' : simple lat-lon grid (1 face). \n');  %
6  choiceGrid=input(' and type return. ''llc'' is the default.\n');  %expected directory structure:
7  if isempty(choiceGrid); choiceGrid='llc'; end;  %
8  if strcmp(choiceGrid,'llc'); choiceGrid='v4';  %  gcmfaces     (codes)
9  elseif strcmp(choiceGrid,'ll'); choiceGrid='v3';  %  MITprof      (codes)
10  else; error('wrong grid choice');  %  nctiles_grid (ECCO v4 grid)
11  end;  %  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    %
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 20  verbose=input(' and/or type return. 0 is Line 41  verbose=input(' and/or type return. 0 is
41  if isempty(verbose); verbose=0; end;  if isempty(verbose); verbose=0; end;
42    
43  %initialize environment variables and mygrid  %initialize environment variables and mygrid
44  gcmfaces_global;  addpath gcmfaces;
45    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  basic_diags_compute_v3_or_v4(choiceGrid);  addpath MITprof;
53  basic_diags_display_v3_or_v4(choiceGrid);  if ~isempty(which('MITprof_global'));
54        MITprof_global;
55        if myenv.verbose>0;
56            gcmfaces_msg('* set MITprof path by calling MITprof_global');
57        end;
58    end;
59    
60  example_bin_average(choiceGrid,1);%incl. call to example_smooth.m;  fprintf('\n\n');
61  example_griddata(choiceGrid);  gcmfaces_msg('/////////////////////////////////////');
62  example_interp(choiceGrid);  gcmfaces_msg('demo of plotting routines','// PART 1 :');
63  example_faces2latlon2faces(choiceGrid);  gcmfaces_msg('please hit return','// >> ');
64    gcmfaces_msg('/////////////////////////////////////');
65    fprintf('\n\n');
66    pause;
67    plot_one_field;
68    
69    fprintf('\n\n');
70    gcmfaces_msg('///////////////////////////////////////////');
71    gcmfaces_msg('demo of remapping and smoothing','// PART 2 :');
72    gcmfaces_msg('please hit return','// >> ');
73    gcmfaces_msg('///////////////////////////////////////////');
74    fprintf('\n\n');
75    pause;
76    
77    if 1;%one method to map lon-lat fields to gcmfaces
78        example_remap(1);
79    else;%alternatives
80        example_griddata;
81        example_interp;
82        example_faces2latlon2faces;
83    end;
84    example_bin_average(1);%incl. call to example_smooth.m;
85    
86  plot_one_field(choiceGrid,0);%%incl. call to m_map_gcmfaces;  fprintf('\n\n');
87    gcmfaces_msg('/////////////////////////////////////////');
88    gcmfaces_msg('demo of transport computations','// PART 3 :');
89    gcmfaces_msg('please hit return','// >> ');
90    gcmfaces_msg('/////////////////////////////////////////');
91    fprintf('\n\n');
92    pause;
93    diags=example_transports;
94    if ~isempty(diags);
95        example_transports_disp(diags);
96    end;
97    
98    fprintf('\n\n');
99    gcmfaces_msg('/////////////////////////////////////////');
100    gcmfaces_msg('demo of budget computations','// PART 4 :');
101    gcmfaces_msg('please hit return','// >> ');
102    gcmfaces_msg('/////////////////////////////////////////');
103    fprintf('\n\n');
104    pause;
105    example_budget;
106    
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22