function outp = rd_ac9(thisfile)

    disp(['Reading ' thisfile ' using rd_ac9...']);
    
    slashes = sort([strfind(thisfile, '\') strfind(thisfile, '/')]);
    thisfilepath = thisfile(1:slashes(end));
    fname = thisfile((slashes(end)+1):end);

    tmppdir = ['tmp\'];
    if ~exist(tmppdir), mkdir(tmppdir); end

    err = dos(['WETPROCE ac90266.dev ' thisfile ' ' tmppdir fname '.out']);

    % Read in Wetview format (based on routine by E Boss, rd_wetview_new.m)
    fid=fopen([tmppdir fname '.out'], 'r');
    dat = [];
    for i = 1:31   %skip first 31 lines
       junk = fgetl(fid);
    end
    
    buf = sscanf(fgetl(fid), '%f');
    if (length(buf)>=19)
       dat = [dat; buf([1 8:10 14:16 2:4 17:19 5:7 11:13])'];
    end
    while(~feof(fid))
        buf = sscanf(fgetl(fid), '%f');
        if (length(buf)>=19)
            dat = [dat; buf([1 8:10 14:16 2:4 17:19 5:7 11:13])'];
        end
    end
    fclose(fid);

    outp = dat;
    
return
    
% function outp = rd_ac9(thisfile)
% 
%     disp('Using parse_ac9');
% 
%     time0 = timestamp_mininode(thisfile);
%     
%     slashes = sort([strfind(thisfile, '\') strfind(thisfile, '/')]);
%     thisfilepath = thisfile(1:slashes(end));
%     fname = thisfile((slashes(end)+1):end);
% 
%     tmppdir = ['tmp\'];
%     if ~exist(tmppdir), mkdir(tmppdir); end
% 
%     err = dos(['WETPROCE ac90266.dev ' thisfile ' ' tmppdir fname '.out']);
% 
%     % Read in Wetview format (based on routine by E Boss, rd_wetview_new.m)
%     fid=fopen([tmppdir fname '.out'], 'r');
%     dat = [];
%     for i = 1:31   %skip first 31 lines
%        junk = fgetl(fid);
%     end
%     
%     buf = sscanf(fgetl(fid), '%f');
%     if (length(buf)>=19)
%        dat = [dat; buf([1 8:10 14:16 2:4 17:19 5:7 11:13])'];
%     end
%     while(~feof(fid))
%         buf = sscanf(fgetl(fid), '%f');
%         if (length(buf)>=19)
%             dat = [dat; buf([1 8:10 14:16 2:4 17:19 5:7 11:13])'];
%         end
%     end
%     fclose(fid);
% 
%     % use time column in wetview file (millisec)
%     time = time0 + (dat(:,1)/1000)/(60*60*24);
%     
%     disp(['** ' num2str(size(dat,1)) ' samples; time ' num2str(time0) ' to ' num2str(time(end))])
% 
%     outp = [time dat];
%     
% return
%     


