1 |
|
|
2 |
%to load the ECCO v4 / LLC90 grid: |
%call sequence: |
3 |
% |
% |
4 |
%nF=5; |
% addpath gcmfaces; gcmfaces_demo; |
5 |
%fileFormat='compact'; |
% |
6 |
%dirGrid=fullfile([myenv.gcmfaces_dir '../']); |
%expected directory structure: |
7 |
%dirGrid=fullfile(dirGrid,'GRID/'); |
% |
8 |
% |
% gcmfaces (codes) |
9 |
%to instead load e.g. your MITgcm lat-lon grid, |
% MITprof (codes) |
10 |
%you would edit the example routines to look like: |
% nctiles_grid (ECCO v4 grid) |
11 |
% |
% release1 (ECCO v4 output) |
12 |
% nF=1; |
% sample_input (test material) |
13 |
% fileFormat='straight'; |
% |
14 |
% dirGrid='./'; |
%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: |
%to activate example_transports: |
21 |
% |
% |
22 |
%mkdir release1 |
% mkdir release1 |
23 |
%wget --recursive ftp://mit.ecco-group.org/ecco_for_las/version_4/release1/nctiles_climatology |
% 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/. |
% mv mit.ecco-group.org/gforget/nctiles_climatology release1/. |
25 |
%rm -rf mit.ecco-group.org |
% rm -rf mit.ecco-group.org |
26 |
% |
% |
27 |
%to activate example_budget: |
%to activate example_budget: |
28 |
% |
% |
29 |
%wget --recursive ftp://mit.ecco-group.org/gforget/nctiles_budget |
% wget --recursive ftp://mit.ecco-group.org/gforget/nctiles_budget_2d |
30 |
%mv mit.ecco-group.org/gforget/nctiles_budget gcmfaces/sample_input/. |
% mv mit.ecco-group.org/gforget/nctiles_budget_2d sample_input/. |
31 |
%rm -rf mit.ecco-group.org |
% rm -rf mit.ecco-group.org |
32 |
% |
% |
|
%to activate example_remap: |
|
|
% |
|
|
%wget ftp://mit.ecco-group.org/gforget/testcase_remap.mat |
|
|
%mv testcase_remap.mat gcmfaces/sample_input/. |
|
33 |
|
|
34 |
%choose verbose 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'); |
64 |
gcmfaces_msg('/////////////////////////////////////'); |
gcmfaces_msg('/////////////////////////////////////'); |
65 |
fprintf('\n\n'); |
fprintf('\n\n'); |
66 |
pause; |
pause; |
67 |
plot_one_field; |
example_display; |
68 |
|
|
69 |
fprintf('\n\n'); |
fprintf('\n\n'); |
70 |
gcmfaces_msg('///////////////////////////////////////////'); |
gcmfaces_msg('///////////////////////////////////////////'); |
71 |
gcmfaces_msg('demo of remapping and smoothing','// 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'); |
fprintf('\n\n'); |
75 |
pause; |
pause; |
76 |
|
example_remap(0); |
77 |
|
|
78 |
if 1;%one method to map lon-lat fields to gcmfaces |
fprintf('\n\n'); |
79 |
example_remap(1); |
gcmfaces_msg('///////////////////////////////////////////'); |
80 |
else;%alternatives |
gcmfaces_msg('demo of bin averaging (gridding)','// PART 2 :'); |
81 |
example_griddata; |
gcmfaces_msg('please hit return','// >> '); |
82 |
example_interp; |
gcmfaces_msg('///////////////////////////////////////////'); |
83 |
example_faces2latlon2faces; |
fprintf('\n\n'); |
84 |
end; |
pause; |
85 |
example_bin_average(1);%incl. call to example_smooth.m; |
%alternatives |
86 |
|
%example_griddata; |
87 |
|
%example_interp; |
88 |
|
%example_faces2latlon2faces; |
89 |
|
fld=example_bin_average; |
90 |
|
|
91 |
|
fprintf('\n\n'); |
92 |
|
gcmfaces_msg('///////////////////////////////////////////'); |
93 |
|
gcmfaces_msg('demo of smoothing ','// PART 2 :'); |
94 |
|
gcmfaces_msg('please hit return','// >> '); |
95 |
|
gcmfaces_msg('///////////////////////////////////////////'); |
96 |
|
fprintf('\n\n'); |
97 |
|
pause; |
98 |
|
example_smooth(fld); |
99 |
|
|
100 |
fprintf('\n\n'); |
fprintf('\n\n'); |
101 |
gcmfaces_msg('/////////////////////////////////////////'); |
gcmfaces_msg('/////////////////////////////////////////'); |
105 |
fprintf('\n\n'); |
fprintf('\n\n'); |
106 |
pause; |
pause; |
107 |
diags=example_transports; |
diags=example_transports; |
108 |
|
|
109 |
if ~isempty(diags); |
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); |
example_transports_disp(diags); |
117 |
end; |
end; |
118 |
|
|