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

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

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


Revision 1.10 - (show annotations) (download)
Fri Mar 20 23:20:04 2015 UTC (10 years, 3 months ago) by gforget
Branch: MAIN
CVS Tags: checkpoint65r, checkpoint65p, checkpoint65q
Changes since 1.9: +40 -15 lines
- gcmfaces_demo.m : add example_remap and add download directions as comments

1
2 %to load the ECCO v4 / LLC90 grid:
3 %
4 %nF=5;
5 %fileFormat='compact';
6 %dirGrid=fullfile([myenv.gcmfaces_dir '../']);
7 %dirGrid=fullfile(dirGrid,'GRID/');
8 %
9 %to instead load e.g. your MITgcm lat-lon grid,
10 %you would edit the example routines to look like:
11 %
12 % nF=1;
13 % fileFormat='straight';
14 % dirGrid='./';
15 %
16 %to activate example_transports:
17 %
18 %mkdir release1
19 %wget --recursive ftp://mit.ecco-group.org/ecco_for_las/version_4/release1/nctiles_climatology
20 %mv mit.ecco-group.org/gforget/nctiles_climatology release1/.
21 %rm -rf mit.ecco-group.org
22 %
23 %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

  ViewVC Help
Powered by ViewVC 1.1.22