enderton/Diagnostics/DiagLoadMonitor.m

function [data,time] = ...
    DiagLoadMonitor(fln,mnchandle,dad,itr,tst,SecPerYear,DiagDebug);

% Read in files names.
%files = dir([dad,'/',mnchandle]);
filesin=ls([dad,'/',mnchandle]);
index=1;
while ~isempty(filesin)
    [token,filesin] = strtok(filesin);
    if ~isempty(filesin)
        files(index).name=token; index=index+1;
    end
end


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                             Read in data                                %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Initialize result arrays.
if isempty(itr),
    time=[];
    data=[];
else
    time=NaN.*zeros(size(itr));
    data=NaN.*zeros(size(itr));
end

for ifile = 1:length(files)
    
    % Open monitor file.
    if DiagDebug, disp(['Entering file:  ',files(ifile).name]); end
    nc=netcdf(files(ifile).name,'read');
    
    % Read time and data information.
    nciter=nc{'T'}; if isempty(nciter), nciter=nc{'iter'}; end
    ncdata=nc{fln};
    if isempty(nciter), error('Monitor time axis not found!'); end
    if isempty(ncdata), error(['Monitor field not found:  ',fln]); end
    
    % Load desired indecies.
    if isempty(itr)
        if isempty(time),
            time=nciter(:);
            data=ncdata(:);
        else
            index=find(~ismember(nciter(:),time));
            nciter=nciter(index); time=[time;nciter];
            ncdata=ncdata(index); data=[data;ncdata];
        end
    else
        [test,loc]=ismember(itr,nciter(:));
        index=find(test);
        if ~isempty(index)
            time(index)=nciter(loc(index));
            data(index)=ncdata(loc(index));
        end
    end
    
    close(nc);
end

% Check for missing indecies.
if ~isempty(find(isnan(time))) || ~isempty(find(isnan(data)))
    error('Missing monitor data for specified indecies!');
end

time=time*tst./SecPerYear;
1	enderton	1.1	function [data,time] = ...
2			DiagLoadMonitor(fln,mnchandle,dad,itr,tst,SecPerYear,DiagDebug);
3
4			% Read in files names.
5			%files = dir([dad,'/',mnchandle]);
6			filesin=ls([dad,'/',mnchandle]);
7			index=1;
8			while ~isempty(filesin)
9			[token,filesin] = strtok(filesin);
10			if ~isempty(filesin)
11			files(index).name=token; index=index+1;
12			end
13			end
14
15
16			%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
17			% Read in data %
18			%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
19
20			% Initialize result arrays.
21			if isempty(itr),
22			time=[];
23			data=[];
24			else
25			time=NaN.*zeros(size(itr));
26			data=NaN.*zeros(size(itr));
27			end
28
29			for ifile = 1:length(files)
30
31			% Open monitor file.
32			if DiagDebug, disp(['Entering file: ',files(ifile).name]); end
33			nc=netcdf(files(ifile).name,'read');
34
35			% Read time and data information.
36			nciter=nc{'T'}; if isempty(nciter), nciter=nc{'iter'}; end
37			ncdata=nc{fln};
38			if isempty(nciter), error('Monitor time axis not found!'); end
39			if isempty(ncdata), error(['Monitor field not found: ',fln]); end
40
41			% Load desired indecies.
42			if isempty(itr)
43			if isempty(time),
44			time=nciter(:);
45			data=ncdata(:);
46			else
47			index=find(~ismember(nciter(:),time));
48			nciter=nciter(index); time=[time;nciter];
49			ncdata=ncdata(index); data=[data;ncdata];
50			end
51			else
52			[test,loc]=ismember(itr,nciter(:));
53			index=find(test);
54			if ~isempty(index)
55			time(index)=nciter(loc(index));
56			data(index)=ncdata(loc(index));
57			end
58			end
59
60			close(nc);
61			end
62
63			% Check for missing indecies.
64			if ~isempty(find(isnan(time))) \|\| ~isempty(find(isnan(data)))
65			error('Missing monitor data for specified indecies!');
66			end
67
68			time=time*tst./SecPerYear;