--- MITgcm_contrib/enderton/Diagnostics/DiagLoad.m 2005/08/18 15:50:24 1.12 +++ MITgcm_contrib/enderton/Diagnostics/DiagLoad.m 2005/09/05 18:48:27 1.13 @@ -53,17 +53,15 @@ if DiagDebug, disp([' Debug -- Loading GRADS field.']); end [data,xax,yax,zax,months,time,dim] = ... DiagLoadGradsData(Grads,dad,fln,DiagDebug); - if DiagDebug, disp([' Debug -- ''data'' size after load: ',mat2str(size(data))]); end + if DiagDebug, disp([' Debug -- ''data'' size after load: ',... + mat2str(size(data))]); end data = DiagAverage(data,fln,avg,months,ddf,dim); - if DiagDebug, disp([' Debug -- ''data'' size after averaging: ',mat2str(size(data))]); end -% if isequal(slc,'Sur') -% data = data'; -% pltslc='lonlat'; -% else -% error('Loading grads data only set to handle surface plots!'); -% end + if DiagDebug, disp([' Debug -- ''data'' size after averaging: ',... + mat2str(size(data))]); end [data,xax,yax,pltslc] = ... DiagSliceGradsData(fln,flu,slc,data,xax,yax,zax); + if DiagDebug, disp([' Debug -- ''data'' size after slicing: ',... + mat2str(size(data))]); end % Load monitor data. elseif isequal(dat(1:2),'Mn') @@ -419,7 +417,7 @@ data = rdmds([dad,'/',fln,filesuffix],itr); elseif isequal(dfm,'MNC') data = rdmnc_mod2([dad,mnchandle],[fln,filesuffix],'iter','T',itr); - if ~isequal(itr,data.iter'), error('Missing iterations in data!'); end + %if ~isequal(itr,data.iter'), error('Missing iterations in data!'); end eval(['data = data.',fln,filesuffix,';']); else error(['Unrecognized data type: ',dfm]);