function [data] = read_eos_fields_airs(filename);
%
% Read EOS fields for AIRS-MODIS collocation study
%
% Input
%    filename         : eos_hdf file name; string
%
% Output
%    data : cell list containing field values in double precision
%       radiances: [2378x90x135 double]
%   spectral_freq: [2378x1 double]
%        Latitude: [90x135 double]
%            Time: [90x135 double]
%       Longitude: [90x135 double]
%         scanang: [90x135 double]
%        landFrac: [90x135 double]
%             NeN: [2378x1 double]
%       satheight: [135x1 double]
%          satzen: [90x135 double]
%          solzen: [90x135 double]
%           topog: [90x135 double]
%

if nargin == 0
  [fn,dn] = uigetfile('AIRS*');
  filename = fullfile(dn,fn);
end

swathname  = 'L1B_AIRS_Science';

start  = '[]';
stride = '[]';
edge   = '[]';
		
clear flist

% Field name list
flist = {'radiances','spectral_freq','Latitude','Time',...
         'Longitude','scanang','landFrac','NeN',...
	 'satheight','satzen','solzen','topog'};

% Loop through fieldnames
for i=1:length(flist)
  s = ['[data.' flist{i} ']=read_eos_field(''' filename ''',''' swathname ''',''' flist{i} ''', start,stride,edge);'];
  eval(s)
end
