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 |
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(:)); |
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; |