1 |
|
|
2 |
%to load the ECCO v4 / LLC90 grid: |
%to load the ECCO v4 / LLC90 grid: |
|
% nF=5; |
|
|
% fileFormat='compact'; |
|
|
% dirGrid=fullfile([myenv.gcmfaces_dir '../']); |
|
|
% dirGrid=fullfile(dirGrid,'GRID/'); |
|
3 |
% |
% |
4 |
%to instead load e.g. your MITgcm lat-lon grid, |
%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: |
%you would edit the example routines to look like: |
11 |
|
% |
12 |
% nF=1; |
% nF=1; |
13 |
% fileFormat='straight'; |
% fileFormat='straight'; |
14 |
% dirGrid='./'; |
% 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 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'); |
42 |
|
|
43 |
%initialize environment variables and mygrid |
%initialize environment variables and mygrid |
44 |
addpath gcmfaces; |
addpath gcmfaces; |
45 |
gcmfaces_global; |
gcmfaces_global; |
46 |
myenv.issueWarnings=0; |
myenv.issueWarnings=0; |
47 |
myenv.verbose=verbose; |
myenv.verbose=verbose; |
48 |
if myenv.verbose>0; |
if myenv.verbose>0; |
51 |
|
|
52 |
addpath MITprof; |
addpath MITprof; |
53 |
if ~isempty(which('MITprof_global')); |
if ~isempty(which('MITprof_global')); |
54 |
MITprof_global; |
MITprof_global; |
55 |
if myenv.verbose>0; |
if myenv.verbose>0; |
56 |
gcmfaces_msg('* set MITprof path by calling MITprof_global'); |
gcmfaces_msg('* set MITprof path by calling MITprof_global'); |
57 |
end; |
end; |
58 |
end; |
end; |
59 |
|
|
60 |
fprintf('\n\n'); |
fprintf('\n\n'); |
68 |
|
|
69 |
fprintf('\n\n'); |
fprintf('\n\n'); |
70 |
gcmfaces_msg('///////////////////////////////////////////'); |
gcmfaces_msg('///////////////////////////////////////////'); |
71 |
gcmfaces_msg('demo of gridding v interpolating','// PART 2 :'); |
gcmfaces_msg('demo of remapping and smoothing','// 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 |
|
|
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; |
example_bin_average(1);%incl. call to example_smooth.m; |
|
example_griddata; |
|
|
example_interp; |
|
|
example_faces2latlon2faces; |
|
85 |
|
|
86 |
fprintf('\n\n'); |
fprintf('\n\n'); |
87 |
gcmfaces_msg('/////////////////////////////////////////'); |
gcmfaces_msg('/////////////////////////////////////////'); |