189 |
V = DiagAverage(V,fln,avg,months,ddf,Dim); |
V = DiagAverage(V,fln,avg,months,ddf,Dim); |
190 |
[XC,XG,YC,YG,Ylat,ZC,ZG,RAC,drC,drF,HFacC,HFacW,HFacS,dxG,dyG,dxC,dyC] = ... |
[XC,XG,YC,YG,Ylat,ZC,ZG,RAC,drC,drF,HFacC,HFacW,HFacS,dxG,dyG,dxC,dyC] = ... |
191 |
DiagLoadGridData(LoadGridData,grd,gdf,flu,GridSuffix,ZcordFile); |
DiagLoadGridData(LoadGridData,grd,gdf,flu,GridSuffix,ZcordFile); |
192 |
%[U,V]=uvcube2latlon(XC,YC,U,V,XL,YL); |
[U,V]=uvcube2latlon(XC,YC,U,V,XL,YL); |
193 |
if ismember(fln,{'U','uVel','fizhi_U'}) |
if ismember(fln,{'U','uVel','fizhi_U'}) |
194 |
data = U; |
data = U; |
195 |
elseif ismember(fln,{'V','vVel','fizhi_V'}) |
elseif ismember(fln,{'V','vVel','fizhi_V'}) |
370 |
if DiagDebug, disp([' Debug -- ''data'' size after vector manipulation: ',mat2str(size(data))]); end |
if DiagDebug, disp([' Debug -- ''data'' size after vector manipulation: ',mat2str(size(data))]); end |
371 |
else |
else |
372 |
% Sequence here for Vector = 0 |
% Sequence here for Vector = 0 |
373 |
if ~isequal(Index,0) |
if ~isequal(Index,0) & isequal(size(Index),[1 1]) & isequal(sign(Index),1) |
374 |
if isequal(Dim,2), data = squeeze(data(:,:,Index,:)); |
if isequal(Dim,2), data = squeeze(data(:,:,Index,:)); |
375 |
elseif isequal(Dim,3), data = squeeze(data(:,:,:,Index,:)); end |
elseif isequal(Dim,3), data = squeeze(data(:,:,:,Index,:)); end |
376 |
|
elseif isequal(size(Index),[1 1]) & ~isequal(sign(Index),1) |
377 |
|
if isequal(Dim,2), data = sign(Index(1)).*squeeze(data(:,:,abs(Index(1)),:)); |
378 |
|
elseif isequal(Dim,3), data = sign(Index(1)).*squeeze(data(:,:,:,abs(Index(1)),:)); end |
379 |
|
elseif isequal(size(Index),[1 2]) |
380 |
|
if isequal(Dim,2), data = sign(Index(1)).*squeeze(data(:,:,abs(Index(1)),:)) + sign(Index(2)).*squeeze(data(:,:,abs(Index(2)),:)); |
381 |
|
elseif isequal(Dim,3), data = sign(Index(1)).*squeeze(data(:,:,:,abs(Index(1)),:)) + sign(Index(2)).*squeeze(data(:,:,:,abs(Index(2)),:)); end |
382 |
end |
end |
383 |
if DiagDebug, disp([' Debug -- ''data'' size after indexing: ',mat2str(size(data))]); end |
if DiagDebug, disp([' Debug -- ''data'' size after indexing: ',mat2str(size(data))]); end |
384 |
data = DiagAverage(data,fln,avg,months,ddf,Dim); |
data = DiagAverage(data,fln,avg,months,ddf,Dim); |