53 |
if DiagDebug, disp([' Debug -- Loading GRADS field.']); end |
if DiagDebug, disp([' Debug -- Loading GRADS field.']); end |
54 |
[data,xax,yax,zax,months,time,dim] = ... |
[data,xax,yax,zax,months,time,dim] = ... |
55 |
DiagLoadGradsData(Grads,dad,fln,DiagDebug); |
DiagLoadGradsData(Grads,dad,fln,DiagDebug); |
56 |
if DiagDebug, disp([' Debug -- ''data'' size after load: ',mat2str(size(data))]); end |
if DiagDebug, disp([' Debug -- ''data'' size after load: ',... |
57 |
|
mat2str(size(data))]); end |
58 |
data = DiagAverage(data,fln,avg,months,ddf,dim); |
data = DiagAverage(data,fln,avg,months,ddf,dim); |
59 |
if DiagDebug, disp([' Debug -- ''data'' size after averaging: ',mat2str(size(data))]); end |
if DiagDebug, disp([' Debug -- ''data'' size after averaging: ',... |
60 |
if isequal(slc,'Sur') |
mat2str(size(data))]); end |
61 |
data = data'; |
[data,xax,yax,pltslc] = ... |
62 |
pltslc='lonlat'; |
DiagSliceGradsData(fln,flu,slc,data,xax,yax,zax); |
63 |
else |
if DiagDebug, disp([' Debug -- ''data'' size after slicing: ',... |
64 |
error('Loading grads data only set to handle surface plots!'); |
mat2str(size(data))]); end |
|
end |
|
65 |
|
|
66 |
% Load monitor data. |
% Load monitor data. |
67 |
elseif isequal(dat(1:2),'Mn') |
elseif isequal(dat(1:2),'Mn') |
417 |
data = rdmds([dad,'/',fln,filesuffix],itr); |
data = rdmds([dad,'/',fln,filesuffix],itr); |
418 |
elseif isequal(dfm,'MNC') |
elseif isequal(dfm,'MNC') |
419 |
data = rdmnc_mod2([dad,mnchandle],[fln,filesuffix],'iter','T',itr); |
data = rdmnc_mod2([dad,mnchandle],[fln,filesuffix],'iter','T',itr); |
420 |
if ~isequal(itr,data.iter'), error('Missing iterations in data!'); end |
%if ~isequal(itr,data.iter'), error('Missing iterations in data!'); end |
421 |
eval(['data = data.',fln,filesuffix,';']); |
eval(['data = data.',fln,filesuffix,';']); |
422 |
else |
else |
423 |
error(['Unrecognized data type: ',dfm]); |
error(['Unrecognized data type: ',dfm]); |