/[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.10 by gforget, Fri Mar 20 23:20:04 2015 UTC
# Line 1  Line 1 
1    
2  %environment variables:  %to load the ECCO v4 / LLC90 grid:
3  global myenv;  %
4  myenv.gcmfaces_verbose=0;  %nF=5;
5  myenv.gcmfaces_dir=[pwd '/gcmfaces/'];  %fileFormat='compact';
6  myenv.MITprof_verbose=0;  %dirGrid=fullfile([myenv.gcmfaces_dir '../']);
7  myenv.MITprof_dir=[pwd '/gcmfaces/'];  %dirGrid=fullfile(dirGrid,'GRID/');
8  myenv.useNativeMatlabNetcdf = ~isempty(which('netcdf.open'));  %
9    %to instead load e.g. your MITgcm lat-lon grid,
10  %re-init mygrid:  %you would edit the example routines to look like:
11  global mygrid; mygrid=[];  %
12    % nF=1;
13  choiceGrid='v4';  % fileFormat='straight';
14    % dirGrid='./';
15  %transform those to functions with v3/v4 arguments:  %
16  basic_diags_compute_v3_or_v4(choiceGrid);  %to activate example_transports:
17  basic_diags_display_v3_or_v4(choiceGrid);  %
18    %mkdir release1
19  example_bin_average(choiceGrid,1);%incl. call to example_smooth.m;  %wget --recursive ftp://mit.ecco-group.org/ecco_for_las/version_4/release1/nctiles_climatology
20  example_griddata(choiceGrid);  %mv mit.ecco-group.org/gforget/nctiles_climatology release1/.
21  example_interp(choiceGrid);  %rm -rf mit.ecco-group.org
22    %
23  plot_one_field(choiceGrid,0);%%incl. call to m_map_gcmfaces;  %to activate example_budget:
24    %
25    %wget --recursive ftp://mit.ecco-group.org/gforget/nctiles_budget
26    %mv mit.ecco-group.org/gforget/nctiles_budget gcmfaces/sample_input/.
27    %rm -rf mit.ecco-group.org
28    %
29    %to activate example_remap:
30    %
31    %wget ftp://mit.ecco-group.org/gforget/testcase_remap.mat
32    %mv testcase_remap.mat gcmfaces/sample_input/.
33    
34    %choose verbose level
35    fprintf('\n Please set the amount of explanatory text display :\n');
36    fprintf('    0: none.\n');
37    fprintf('    1: comments.\n');
38    fprintf('    2: comments preceeded with calling sequence.\n');
39    fprintf('    3: same as 2, but preceeded with pause.\n');
40    verbose=input(' and/or type return. 0 is the default. \n');
41    if isempty(verbose); verbose=0; end;
42    
43    %initialize environment variables and mygrid
44    addpath gcmfaces;
45    gcmfaces_global;
46    myenv.issueWarnings=0;
47    myenv.verbose=verbose;
48    if myenv.verbose>0;
49        gcmfaces_msg('* set path and environment variables (myenv) by calling gcmfaces_global');
50    end;
51    
52    addpath MITprof;
53    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    fprintf('\n\n');
61    gcmfaces_msg('/////////////////////////////////////');
62    gcmfaces_msg('demo of plotting routines','// PART 1 :');
63    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    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.1  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22