19 |
vecCaxis=[[3 3 1];[0.5 0.5 0.2]]; |
vecCaxis=[[3 3 1];[0.5 0.5 0.2]]; |
20 |
%if nargin>6; vecCaxis=varargin{1}; else; vecCaxis=[]; end; |
%if nargin>6; vecCaxis=varargin{1}; else; vecCaxis=[]; end; |
21 |
|
|
22 |
%load paths and grid: |
gcmfaces_global; |
23 |
%-------------------- |
global myparms; |
24 |
global mygrid; if isempty(mygrid); grid_load('GRID/',5,'compact'); end; |
|
|
|
|
|
%set more optional paramaters to default values |
|
|
global myparms; year0=myparms.yearInAve(1); year1=myparms.yearInAve(2); |
|
|
listData={'argo_*_1992_to_2007*','argo_*_2008_to_2010*',... |
|
|
'argo_*_2011_to_2012*','climode_*',... |
|
|
'ctd_*','itp_*','seals_*','xbt_*'}; |
|
25 |
listVar={'T','S'}; |
listVar={'T','S'}; |
26 |
%shorter test case: listData={'argo_in*'}; listVar={'T'}; |
%shorter test case: listData={'argo_in*'}; listVar={'T'}; |
27 |
dirData='mat/profiles/output/'; |
dirData='mat/profiles/output/'; |
44 |
end; |
end; |
45 |
end; |
end; |
46 |
|
|
47 |
|
if isempty(myparms)|isempty(mygrid); |
48 |
|
load([dirMat 'diags_grid_parms.mat']); |
49 |
|
end; |
50 |
|
|
51 |
|
if isempty(whos('year0')); |
52 |
|
year0=myparms.yearInAve(1); year1=myparms.yearInAve(2); |
53 |
|
end; |
54 |
|
|
55 |
|
if isempty(whos('listData')); |
56 |
|
listData=dir([dirMat 'profiles/output/*.nc']); |
57 |
|
listData={listData(:).name}; |
58 |
|
for ff=1:length(listData); listData{ff}=[listData{ff}(1:end-3) '*']; end; |
59 |
|
end; |
60 |
|
|
61 |
if doComp; |
if doComp; |
62 |
|
|
63 |
%initialize delaunay triangulation: |
%initialize delaunay triangulation: |