function HTMon_data(fflfile, outfile, startday, endday)
%HTMon_data Creates TimeMon data file for LHO S5 run.
%
% HTMon_data(fflfile, outfile, startday, endday) reads frame list from
% fflfile for data from startday to endday (time unit: time passed in 
% days since the beginning of S5), and creates a data file for
% S5 LHO data for the time channels.
%
% The output file will have following vectors holding time-data:
% 'H1EX_mean','H1EX_max','H1EX_min', 'H1EX_n', 'H1EX_rms', 'H1EY_mean',
% 'H1EY_max', 'H1EY_min', 'H1EY_n', 'H1EY_rms', 'H2MX_mean', 'H2MX_max',
% 'H2MX_min', 'H2MX_n', 'H2MX_rms', 'H2MY_mean', 'H2MY_max', 'H2MY_min',
% 'H2MX_n', 'H2MX_rms', 'H1LV_mean', 'H1LV_max', 'H1LV_min', 'H1LV_n',
% 'H1LV_rms','H2LV_mean', 'H2LV_max', 'H2LV_min', 'H2LV_n', 'H2LV_rms'.
%
% Also, a time channel created by this script is saved as vector 'H_time',
% while the GPS start and end time are save as variables tHStart and
% tHLength.
%
% To create ffl: ls /home/rubab/TimeMon/LHO/H-M-815/*.gwf >> temp.ffl
%                ls /home/rubab/TimeMon/LHO/H-M-816/*.gwf >> temp.ffl
%                ...
%                /bin/FrDump -d 0 -i temp.ffl > HTimeMon.ffl                
%
% Example: HTMon_data('HTimeMon.ffl', 'HTimeMon_all.mat', 1, 151)
%
% See also: HTMon_plots, LTMon_data
%
% Rubab Khan, GECo, Columbia
% https://geco.phys.columbia.edu/~rubab/
% March 28, 2006.

warning off all;

if startday > 1000
    startday = (startday - 815158800)/(60*60*24);
end
if endday > 1000
    endday = (endday - 815158800)/(60*60*24);
end

stD = (startday - 1);
enD = endday;
%tA=815155213;                                % official S5 start
tHStart = 815158800  + stD*24*60*60;        % aligned with frames
tHLength = (enD-stD)*24*60*60;            
% tHLength = 878644800 - tHStart; % To get all data available till end
H_time = (tHStart:60:(tHStart+tHLength-1))';  % Creating time channel

% Original channel names
ChanData=[
    'H1:DMT-TIME_EX.   ';
    'H1:DMT-TIME_EY.   ';
    'H2:DMT-TIME_MX.   ';
    'H2:DMT-TIME_MY.   ';
    'H1:DMT-TIME_LVEA. ';
    'H2:DMT-TIME_LVEA. '];

% Channel variations
DataType=[
        'mean  ';
        'max   ';
        'min   ';
        'n     ';
        'rms   '];

CD = cellstr(ChanData);
DT = cellstr(DataType);

% Initiate the vectors
H1EX_mean = []; H1EX_max = []; H1EX_min = []; H1EX_n = []; H1EX_rms = [];
H1EY_mean = []; H1EY_max = []; H1EY_min = []; H1EY_n = []; H1EY_rms = [];
H2MX_mean = []; H2MX_max = []; H2MX_min = []; H2MX_n = []; H2MX_rms = [];
H2MY_mean = []; H2MY_max = []; H2MY_min = []; H2MY_n = []; H2MY_rms = [];
H1LV_mean = []; H1LV_max = []; H1LV_min = []; H1LV_n = []; H1LV_rms = [];
H2LV_mean = []; H2LV_max = []; H2LV_min = []; H2LV_n = []; H2LV_rms = [];

for i = 1:6    % For six channels
    for j = 1:5 % For five variations
         Chn = sprintf('%s', char(CD(i)), char(DT(j))); % Specific channel name
         t = frgetvect(fflfile, Chn, tHStart, tHLength); % Read data
         if i == 1
             if j == 1 H1EX_mean = t;
             elseif j == 2 H1EX_max = t;
             elseif j == 3 H1EX_min = t;
             elseif j == 4 H1EX_rms = t;
             elseif j == 5 H1EX_n = t;
             end
         elseif i ==2
             if j == 1 H1EY_mean = t;
             elseif j == 2 H1EY_max = t;
             elseif j == 3 H1EY_min = t;
             elseif j == 4 H1EY_rms = t;
             elseif j == 5 H1EY_n = t;
             end
         elseif i ==3
             if j == 1 H2MX_mean = t;
             elseif j == 2 H2MX_max = t;
             elseif j == 3 H2MX_min = t;
             elseif j == 4 H2MX_rms = t;
             elseif j == 5 H2MX_n = t;
             end
         elseif i == 4
             if j == 1 H2MY_mean = t;
             elseif j == 2 H2MY_max = t;
             elseif j == 3 H2MY_min = t;
             elseif j == 4 H2MY_rms = t;
             elseif j == 5 H2MY_n = t;
             end
         elseif i ==5
             if j == 1 H1LV_mean = t;
             elseif j == 2 H1LV_max = t;
             elseif j == 3 H1LV_min = t;
             elseif j == 4 H1LV_rms = t;
             elseif j == 5 H1LV_n = t;
             end
         elseif i ==6
             if j == 1 H2LV_mean = t;
             elseif j == 2 H2LV_max = t;
             elseif j == 3 H2LV_min = t;
             elseif j == 4 H2LV_rms = t;
             elseif j == 5 H2LV_n = t;
             end
         end
    end
end
   
save(outfile, 'tHStart', 'tHLength', 'H_time', 'H1EX_mean','H1EX_max','H1EX_min', 'H1EX_n', 'H1EX_rms', 'H1EY_mean', 'H1EY_max', 'H1EY_min', 'H1EY_n', 'H1EY_rms', 'H2MX_mean', 'H2MX_max', 'H2MX_min', 'H2MX_n', 'H2MX_rms', 'H2MY_mean', 'H2MY_max', 'H2MY_min', 'H2MX_n', 'H2MX_rms', 'H1LV_mean', 'H1LV_max', 'H1LV_min', 'H1LV_n', 'H1LV_rms','H2LV_mean', 'H2LV_max', 'H2LV_min', 'H2LV_n', 'H2LV_rms', '-MAT');

return