/[MITgcm]/MITgcm_contrib/enderton/Diagnostics/DiagLoadMonitor.m
ViewVC logotype

Diff of /MITgcm_contrib/enderton/Diagnostics/DiagLoadMonitor.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by enderton, Mon Feb 7 05:31:45 2005 UTC revision 1.3 by enderton, Mon Sep 5 18:50:15 2005 UTC
# Line 33  for ifile = 1:length(files) Line 33  for ifile = 1:length(files)
33      nc=netcdf(files(ifile).name,'read');      nc=netcdf(files(ifile).name,'read');
34            
35      % Read time and data information.      % Read time and data information.
36      nciter=nc{'T'}; if isempty(nciter), nciter=nc{'iter'}; end      nciter=nc{'T'}; if isempty(nciter), nciter=nc{'iter'}.*tst; end
37      ncdata=nc{fln};      ncdata=nc{fln};
38      if isempty(nciter), error('Monitor time axis not found!'); end      if isempty(nciter), error('Monitor time axis not found!'); end
39      if isempty(ncdata), error(['Monitor field not found:  ',fln]); end      if isempty(ncdata), error(['Monitor field not found:  ',fln]); end
# Line 45  for ifile = 1:length(files) Line 45  for ifile = 1:length(files)
45              data=ncdata(:);              data=ncdata(:);
46          else          else
47              index=find(~ismember(nciter(:),time));              index=find(~ismember(nciter(:),time));
48              nciter=nciter(index); time=[time;nciter];                  if ~isempty(index)
49              ncdata=ncdata(index); data=[data;ncdata];                      nciter=nciter(index); time=[time;NaN;nciter];
50                        ncdata=ncdata(index); data=[data;NaN;ncdata];
51                end
52          end          end
53      else      else
54          [test,loc]=ismember(itr,nciter(:));          [test,loc]=ismember(itr,nciter(:));
# Line 61  for ifile = 1:length(files) Line 63  for ifile = 1:length(files)
63  end  end
64    
65  % Check for missing indecies.  % Check for missing indecies.
66  if ~isempty(find(isnan(time))) || ~isempty(find(isnan(data)))  % if ~isempty(find(isnan(time))) || ~isempty(find(isnan(data)))
67      error('Missing monitor data for specified indecies!');  %     error('Missing monitor data for specified indecies!');
68  end  % end
69    
 time=time*tst./SecPerYear;  
70    time=time./SecPerYear;
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22