





function out = draftConvert128_3toAOs(filename)
  
  if nargin == 0
    filename = '/data/home/indiep/DataAnalysis/a109721/DMUDMU_A/LSM10030_tmpk/LSM10030_tmpk_200910052051Z_F200909301020Z_T200909301021Z_hex.txt';
  end
  
  fid = fopen(filename);
  
  format128_3 = '%s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %s %*[^\n]';
  
  C = textscan(fid, format128_3, 'Delimiter', '\t');
  
  fclose(fid);
  
  timeFormat = java.text.SimpleDateFormat('dd-MM-yyyy HH:mm:ss.SSS');
  timeFormat.setTimeZone(java.util.TimeZone.getTimeZone('UTC'));
  timeFormat.parse('01-01-2000 00:00:00.000');
  
  toff = timeFormat.getCalendar.getTimeInMillis;
  
  timeFormat = java.text.SimpleDateFormat('dd MMM yyyy HH:mm:ss.SSS');
  timeFormat.setTimeZone(java.util.TimeZone.getTimeZone('UTC'));
  
  %%%%%%%%%%%%%%%%%%%%%%%%%%   Convert time string   %%%%%%%%%%%%%%%%%%%%%%%%%%
  
  % Column 2
  N     = 2;
  tUDMS = [];
  for ii = 3:numel(C{N})
    timeFormat.parse(C{N}{ii});
    tUDMS = [tUDMS; (timeFormat.getCalendar.getTimeInMillis - toff)/1000];
  end
  
  % Column 3
  N = 3;
  tDMU = hex2x_time(char(C{N}{3:end}));
  
  DMU_UDMS_TIME = ao(xydata(tDMU, tUDMS));
  DMU_UDMS_TIME.setName('DMU_UDMS_TIME', 'internal');
  
  t0 = time();
  t0.setTimezone('UTC');
  t0.setEpochtime(round(toff + tDMU(1)*1000));
  
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%   Convert flags   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  
  % Column 5
  N = 5;
  DMU_X1_C1_GT_5 = createFlagAO(N, @hex2dec);
  
  % Column 6
  N = 6;
  DMU_X1_C1_GT_42 = createFlagAO(N, @hex2dec);
  
  % Column 7
  N = 7;
  DMU_X1_C1_LT_38 = createFlagAO(N, @hex2dec);
  
  % Column 8
  N = 8;
  DMU_X1_C1_GT_70 = createFlagAO(N, @hex2dec);
  
  % Column 9
  N = 9;
  DMU_X1_ERR1 = createFlagAO(N, @hex2dec);
  
  % Column 10
  N = 10;
  DMU_X1_ERR2 = createFlagAO(N, @hex2dec);
  
  % Column 11
  N = 11;
  DMU_X1_ERR3 = createFlagAO(N, @hex2dec);
  
  % Column 12
  N = 12;
  DMU_X1_ETA_PHI = createFlagAO(N, @hex2dec);
  
  % Column 13
  N = 13;
  DMU_X1_VALIDITY = createFlagAO(N, @hex2dec);
  
  % Column 14
  N = 14;
  DMU_X1_ETA_PHI_VALIDITY = createFlagAO(N, @hex2dec);
  
  % Column 15
  N = 15;
  DMU_X1_C1_GT_10 = createFlagAO(N, @hex2dec);
  
  % Column 16
  N = 16;
  DMU_X1_C1_LT_6 = createFlagAO(N, @hex2dec);
  
  % Column 19
  N = 19;
  DMU_X12_C12_GT_5 = createFlagAO(N, @hex2dec);
  
  % Column 20
  N = 20;
  DMU_X12_C12_GT_42 = createFlagAO(N, @hex2dec);
  
  % Column 21
  N = 21;
  DMU_X12_C12_LT_38 = createFlagAO(N, @hex2dec);
  
  % Column 22
  N = 22;
  DMU_X12_C12_GT_70 = createFlagAO(N, @hex2dec);
  
  % Column 23
  N = 23;
  DMU_X12_ERR4 = createFlagAO(N, @hex2dec);
  
  % Column 24
  N = 24;
  DMU_X12_ERR5 = createFlagAO(N, @hex2dec);
  
  % Column 25
  N = 25;
  DMU_X12_ERR6 = createFlagAO(N, @hex2dec);
  
  % Column 26
  N = 26;
  DMU_X12_ETA_PHI = createFlagAO(N, @hex2dec);
  
  % Column 27
  N = 27;
  DMU_X12_VALIDITY = createFlagAO(N, @hex2dec);
  
  % Column 28
  N = 28;
  DMU_X12_ETA_PHI_VALIDITY = createFlagAO(N, @hex2dec);
  
  % Column 29
  N = 29;
  DMU_X12_C1_GT_10 = createFlagAO(N, @hex2dec);
  
  % Column 30
  N = 30;
  DMU_X12_C1_LT_6 = createFlagAO(N, @hex2dec);
  
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%   Convert data   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  
  % Column 17
  N = 17;
  DMU_X1_FILT = createDataAO(N, @hex2num, 'm');
  
  % Column 31
  N = 31;
  DMU_X12_FILT = createDataAO(N, @hex2num, 'm');
  
  % Column 32
  N = 32;
  DMU_DC_PHI_1_FILT = createDataAO(N, @hex2single, 'rad');
  
  % Column 33
  N = 33;
  DMU_DC_ETA_1_FILT = createDataAO(N, @hex2single, 'rad');
  
  % Column 34
  N = 34;
  DMU_DC_PHI_2_FILT = createDataAO(N, @hex2single, 'rad');
  
  % Column 35
  N = 35;
  DMU_DC_ETA_2_FILT = createDataAO(N, @hex2single, 'rad');
  
  % Column 36
  N = 36;
  DMU_DWS_PHI_1_FILT = createDataAO(N, @hex2single, 'rad');
  
  % Column 37
  N = 37;
  DMU_DWS_ETA_1_FILT = createDataAO(N, @hex2single, 'rad');
  
  % Column 38
  N = 38;
  DMU_DWS_PHI_2_FILT = createDataAO(N, @hex2single, 'rad');
  
  % Column 39
  N = 39;
  DMU_DWS_ETA_2_FILT = createDataAO(N, @hex2single, 'rad');
  
  % Column 40
  N = 40;
  DMU_PSI_F_FILT = createDataAO(N, @hex2num, 'rad');
  
  % Column 41
  N = 41;
  DMU_PSI_R_FILT = createDataAO(N, @hex2num, 'rad');
  
  % Column 42
  N = 42;
  DMU_SIGMA_F_FILT = createDataAO(N, @hex2single, '');
  
  % Column 43
  N = 43;
  DMU_SIGMA_R_FILT = createDataAO(N, @hex2single, '');
  
  out = [...
    DMU_X1_FILT, ...
    DMU_X12_FILT, ...
    DMU_DC_PHI_1_FILT, ...
    DMU_DC_ETA_1_FILT, ...
    DMU_DC_PHI_2_FILT, ...
    DMU_DC_ETA_2_FILT, ...
    DMU_DWS_PHI_1_FILT, ...
    DMU_DWS_ETA_1_FILT, ...
    DMU_DWS_PHI_2_FILT, ...
    DMU_DWS_ETA_2_FILT, ...
    DMU_PSI_F_FILT, ...
    DMU_PSI_R_FILT, ...
    DMU_SIGMA_F_FILT, ...
    DMU_SIGMA_R_FILT, ...
    DMU_X1_C1_GT_5, ...
    DMU_X1_C1_GT_42, ...
    DMU_X1_C1_LT_38, ...
    DMU_X1_C1_GT_70, ...
    DMU_X1_ERR1, ...
    DMU_X1_ERR2, ...
    DMU_X1_ERR3, ...
    DMU_X1_ETA_PHI, ...
    DMU_X1_VALIDITY, ...
    DMU_X1_ETA_PHI_VALIDITY, ...
    DMU_X1_C1_GT_10, ...
    DMU_X1_C1_LT_6, ...
    DMU_X12_C12_GT_5, ...
    DMU_X12_C12_GT_42, ...
    DMU_X12_C12_LT_38, ...
    DMU_X12_C12_GT_70, ...
    DMU_X12_ERR4, ...
    DMU_X12_ERR5, ...
    DMU_X12_ERR6, ...
    DMU_X12_ETA_PHI, ...
    DMU_X12_VALIDITY, ...
    DMU_X12_ETA_PHI_VALIDITY, ...
    DMU_X12_C1_GT_10, ...
    DMU_X12_C1_LT_6, ...
    DMU_UDMS_TIME];
  
  out = addHistoryStep(out, plist('filename', filename), '$Id: draftConvert128_3toAOs.m,v 1.3 2009/10/19 10:35:12 ingo Exp $', []);
  
  out.submit();
  
  function out = createFlagAO(N, convert_fcn)
    d = cdata(convert_fcn(char(C{N}{3:end})));
    out = ao(d);
    out.setName(C{N}{2}, 'internal');
    out.setDescription(C{N}{1}, 'internal');
  end
  
  function out = createDataAO(N, convert_fcn, unit)
    d = tsdata(convert_fcn(char(C{N}{3:end})), 10, t0);
    out = ao(d);
    out.setName(C{N}{2}, 'internal');
    out.setDescription(C{N}{1}, 'internal');
    out.setXunits('s',  'internal');
    out.setYunits(unit, 'internal');
  end
  
end


function x = hex2x_time(s)
  
  % The first 4 bytes are the seconds from the 1 January 2000
  s1 = s(:, 1:8);
  
  % The last byte divided by 256 in milliseconds
  s2 = s(:, 9:10);
  
  x = hex2dec(s1) + hex2dec(s2)/256;
  
end

%HEX2SINGLE Convert single precision IEEE hexadecimal string to number.
%   HEX2SINGLE(S), where S is a 8 character string containing
%   a hexadecimal number, returns a double type number
%   equal to the IEEE single precision
%   floating point number it represents.  Fewer than 8
%   characters are padded on the right with zeros.
%
%   If S is a character array, each row is interpreted as a single
%   precision number (and returned as a double).
%
%   NaNs, infinities and denorms are handled correctly.
%
%   Example:
%       hexsingle2num('40490fdb') returns Pi.
%       hexsingle2num('bf8') returns -1.
%
%   See also HEX2NUM.

% Based on Matlab's hex2num.
% Note: IEEE Standard 754 for floating point numbers
%
%  Floating point numbers are represented as:
%  x = +/- (1+f)*2^e
%
%  doubles: 64 bits
%           Bit 63       (1 bit)  = sign (0=positive, 1=negative)
%           Bit 62 to 52 (11 bits)= exponent biased by 1023
%           Bit 51 to 0  (52 bits)= fraction f of the number 1.f
%  singles: 32 bits
%           Bit 31       (1 bit)  = sign (0=positive, 1=negative)
%           Bit 30 to 23 (8 bits) = exponent biased by 127
%           Bit 22 to 0  (23 bits)= fraction f of the number 1.f
%
% Original file hexsingle2num from Mark Lubinski
% Changed on 19-may-05 by Matthias Noell: denormalized power set 2^-126
function x = hex2single(s)
  if iscellstr(s), s = char(s); end
  if ~ischar(s)
    error('Input to hexsingle2num must be a string.')
  end
  if isempty(s), x = []; return, end
  
  [row,col] = size(s);
  blanks = find(s==' '); % Find the blanks at the end
  if ~isempty(blanks), s(blanks) = '0'; end % Zero pad the shorter hex numbers.
  
  % Convert characters to numeric digits.
  % More than 8 characters are ignored
  % For double: d = zeros(row,16);
  d = zeros(row,8);
  d(:,1:col) = abs(lower(s)) - '0';
  d = d + ('0'+10-'a').*(d>9);
  neg = d(:,1) > 7;
  d(:,1) = d(:,1)-8*neg;
  
  if any(d > 15) | any(d < 0)
    error('Input string to hexsingle2num should have just 0-9, a-f, or A-F.')
  end
  
  % Floating point exponent.
  % For double: e = 16*(16*(d(:,1)-4) + d(:,2)) + d(:,3) + 1;
  % For double: e = 256*d(:,1) + 16*d(:,2) + d(:,3) - 1023;
  expBit = (d(:,3) > 7);
  e = 32*d(:,1) + 2*d(:,2) + expBit - 127;
  d(:,3) = d(:,3)-8*expBit;  % Remove most sig. bit of d(:,3) which belongs to exponent
  
  % Floating point fraction.
  % For double: sixteens = [16;256;4096;65536;1048576;16777216;268435456];
  % For double: sixteens2 = 268435456*sixteens(1:6);
  % For double: multiplier = 1./[sixteens;sixteens2];
  % For double: f = d(:,4:16)*multiplier;
  sixteens = [16;256;4096;65536;1048576;16777216];
  multiplier = 2./[sixteens];
  f = d(:,3:8)*multiplier;
  
  x = zeros(row,1);
  % Scale the fraction by 2 to the exponent.
  % For double: overinf = find((e>1023) & (f==0));
  overinf = find((e>127) & (f==0));
  if ~isempty(overinf), x(overinf) = inf; end
  
  % For double: overNaN = find((e>1023) & (f~=0));
  overNaN = find((e>127) & (f~=0));
  if ~isempty(overNaN), x(overNaN) = NaN; end
  
  % For double: underflow = find(e<-1022);
  underflow = find(e<-126);
  if ~isempty(underflow), x(underflow) = pow2(f(underflow),-126); end
  
  % For double: allothers = find((e<=1023) & (e>=-1022));
  allothers = find((e<=127) & (e>=-126));
  if ~isempty(allothers), x(allothers) = pow2(1+f(allothers),e(allothers)); end
  
  negatives = find(neg);
  if ~isempty(negatives), x(negatives) = -x(negatives); end
  
end
