--- MITgcm_contrib/enderton/Diagnostics/DiagLoadMonitor.m 2005/02/07 16:22:30 1.2 +++ MITgcm_contrib/enderton/Diagnostics/DiagLoadMonitor.m 2005/09/05 18:50:15 1.3 @@ -33,7 +33,7 @@ nc=netcdf(files(ifile).name,'read'); % Read time and data information. - nciter=nc{'T'}; if isempty(nciter), nciter=nc{'iter'}; end + nciter=nc{'T'}; if isempty(nciter), nciter=nc{'iter'}.*tst; end ncdata=nc{fln}; if isempty(nciter), error('Monitor time axis not found!'); end if isempty(ncdata), error(['Monitor field not found: ',fln]); end @@ -45,9 +45,9 @@ data=ncdata(:); else index=find(~ismember(nciter(:),time)); - if ~isempty(index) - nciter=nciter(index); time=[time;nciter]; - ncdata=ncdata(index); data=[data;ncdata]; + if ~isempty(index) + nciter=nciter(index); time=[time;NaN;nciter]; + ncdata=ncdata(index); data=[data;NaN;ncdata]; end end else @@ -63,8 +63,8 @@ end % Check for missing indecies. -if ~isempty(find(isnan(time))) || ~isempty(find(isnan(data))) - error('Missing monitor data for specified indecies!'); -end +% if ~isempty(find(isnan(time))) || ~isempty(find(isnan(data))) +% error('Missing monitor data for specified indecies!'); +% end -time=time*tst./SecPerYear; +time=time./SecPerYear;