clear

tic
listing = dir;

delimiters = [',' char(9)];

ODVfidBuoy = fopen('WF-CANONCruiseSept2010-Buoy.txt', 'w');
ODVfidESP = fopen('WF-CANONCruiseSept2010-ESP.txt', 'w');
ODVfidAPR = fopen('WF-CANONCruiseSept2010_PAR.txt', 'w');
fprintf(ODVfidBuoy,'//This data has been parsed from the original ESP drifter files generated by Labview .VI\r\n');
fprintf(ODVfidESP,'//This data has been parsed from the original ESP drifter files generated by Labview .VI\r\n');
fprintf(ODVfidPAR,'//This data has been parsed from the original ESP drifter files generated by Labview .VI\r\n');
fprintf(ODVfidBuoy,'Cruise\tStation\tType\tyyyy-mm-ddThh:mm:ss.sss\tLongitude[degrees_east]\tLatitude[degrees_north]\tBot. Depth[m]\tGPS Time[hh:mm:ss.ss]\tGPS Long.[degrees E/W]\tGPS Lat.[degrees N/S]');
fprintf(ODVfid,'\tCable Tension[lbf]\tBus Voltage[VDC]');
fprintf(ODVfid,'\tBuoy Roll Angle[deg]\tBuoy Pitch Angle[deg]\tBuoy Yaw Angle[deg True N]\tBuoy X Rate[deg/sec]\tBuoy Y Rate[deg/sec]\tBuoy Z Rate[deg/sec]');
fprintf(ODVfid,'\tBuoy X Accel.[g]\tBuoy Y Accel.[g]\tBuoy Z Accel.[g]\tBuoy North Vel.[m/s]\tBuoy East Vel.[m/s]\tBuoy West Vel.[m/s]\tBuoy xBow Long.[deg]\tBuoy xBow Lat.[deg]\tBuoy xBow Alt.[m]');
fprintf(ODVfid,'\tESP X Accel.[g]\tESP Y Accel.[t]\tESP Z Accel.[g]\tESP X Rate[rad/sec]\tESP Y Rate[rad/sec]\tESP Z Rate[rad/sec]');
fprintf(ODVfid,'\tPAR CH1-412nm\tPAR CH1-443nm\tPAR CH1-490nm\tPAR CH1-510nm\tPAR CH1-555nm\tPAR CH1-665nm\tPAR CH1-PAR\tPAR CH1-Gnd.\tPAR Gain Err.\tPAR Stale');
fprintf(ODVfid,'\r\n');


cruise = 'WF-CANON-Sept2010';
station = 1;
depth = 'na';

numLCTabs = 3;
numXBowTabs = 5;
numILinkTabs = 20;
numPRRTabs = 26;

lat = '';
long = '';

for i=3:length(listing)
%for i=48:49
    LVfid = fopen(listing(i).name, 'r');
    fprintf(1,'Processing file:  %d %s\r\n', i, listing(i).name);
    date = datestr(listing(i).datenum, 29);
    
    lineNum = 0;
    while ~feof(LVfid)
        line = fgetl(LVfid);
        lineNum = lineNum + 1;
        [token, r] = strtok(line, delimiters);
        
        switch token
            case 'Tens/Bus'
                [timeStamp, r] = strtok(r, delimiters);    %Labview Time stamp
                [t, r] = strtok(r, delimiters);
                tension = str2num(t);
                [t, r] = strtok(r, delimiters);
                busVoltage = str2num(t);
                
                type = 'Load Cell';
                if ~strcmp(lat,'')
                    fprintf(ODVfid,'%s\t%d\t%s\t%s\t%13.9f\t%13.9f\t%s\t',cruise, station, type, strcat(date,'T',timeStamp), long, lat, depth);
                    for i=1:numLCTabs;
                        fprintf(ODVfid,'\t');
                    end
                    fprintf(ODVfid,'%5.0f\t%5.1f\t\r\n', tension, busVoltage);
                    station = station + 1;
                end
            case 'xBow GPS'
                if length(line) > 82
                    [timeStamp, r] = strtok(r, delimiters); % Labview time stamp
                    [GGAHeader, r] = strtok(r, delimiters);
                    [temp, r] = strtok(r, delimiters);
                    GGATime = strcat(temp(1:2),':',temp(3:4),':',temp(5:length(temp)));
                    [latTemp, r] = strtok(r, delimiters);
                    if length(latTemp) > 5
                        lat = str2num(latTemp(1:2)) + str2num(latTemp(3:length(latTemp)))/60.0;
                        [NS, r] = strtok(r, delimiters);
                        [longTemp, r] = strtok(r, delimiters);
                        if length(longTemp) > 5
                            long = str2num(longTemp(1:3)) + str2num(longTemp(4:length(longTemp)))/60.0;
                            [EW, r] = strtok(r, delimiters);
                            type = 'xBow GPS';
                            fprintf(ODVfid,'%s\t%d\t%s\t%s\t%13.9f\t%13.9f\t%s\t',cruise, station, type, strcat(date,'T',timeStamp), long, lat, depth);
                            fprintf(ODVfid,'%s\t%13.9f %1s\t%13.9f %1s\t\r\n', GGATime, long, EW, lat, NS);
                            station = station + 1;
                        else
                            fprintf(1,'Skipping line %d in file %s: %s\r\n', lineNum, listing(i).name, line);
                        end
                    else
                            fprintf(1,'Skipping line %d in file %s: %s\r\n', lineNum, listing(i).name, line);
                    end
                else
                    fprintf(1,'Skipping line %d in file %s: %s\r\n', lineNum, listing(i).name, line);
                end
            case 'xBow'
                [timeStamp, r] = strtok(r, delimiters); % Labview time stamp
                [temp, r] = strtok(r, delimiters);
                roll = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                pitch = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                yaw = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                xRate = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                yRate = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                zRate = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                xAcc = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                yAcc = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                zAcc = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                xVeloc = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                yVeloc = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                zVeloc = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                long = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                lat = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                alt = str2num(temp);
                
                type = 'Buoy xBow';
                if ~strcmp(lat,'')
                    fprintf(ODVfid,'%s\t%d\t%s\t%s\t%13.9f\t%13.9f\t%s\t',cruise, station, type, strcat(date,'T',timeStamp), long, lat, depth);
                    for i=1:numXBowTabs;
                        fprintf(ODVfid,'\t');
                    end
                    fprintf(ODVfid,'%9.5f\t%9.5f\t%9.5f\t%9.5f\t%9.5f\t%9.5f\t%9.5f\t%9.5f\t%9.5f\t%9.5f\t%9.5f\t%9.5f\t%13.9f\t%13.9f\t%9.5f\t\r\n', ...
                        roll, pitch, yaw, xRate,yRate, zRate, xAcc, yAcc, zAcc, xVeloc, yVeloc, zVeloc, long, lat, alt);
                    station = station + 1;
                end
            case 'iLink'
                [timeStamp, r] = strtok(r, delimiters); % Labview time stamp
                [temp, r] = strtok(r, delimiters);
                xAcc = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                yAcc = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                zAcc = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                xRate = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                yRate = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                zRate = str2num(temp);
                
                type = 'ESP iLink';
                if ~strcmp(lat,'')
                    fprintf(ODVfid,'%s\t%d\t%s\t%s\t%13.9f\t%13.9f\t%s\t',cruise, station, type, strcat(date,'T',timeStamp), long, lat, depth);
                    for i=1:numILinkTabs;
                        fprintf(ODVfid,'\t');
                    end
                    fprintf(ODVfid,'%9.5f\t%9.5f\t%9.5f\t%9.5f\t%9.5f\t%9.5f\r\n', ...
                        xAcc, yAcc, zAcc, xRate,yRate, zRate);
                    station = station + 1;
                end
            case 'PRR'
                [timeStamp, r] = strtok(r, delimiters); % Labview time stamp
                [temp, r] = strtok(r, delimiters);
                ch1 = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                ch2 = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                ch3 = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                ch4 = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                ch5 = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                ch6 = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                ch7 = str2num(temp);
                [temp, r] = strtok(r, delimiters);
                ch8 = str2num(temp);
                [PARGainErr, r] = strtok(r, delimiters);
                [PARStale, r] = strtok(r, delimiters);
                
                type = 'Biosph. PRR610';
                if ~strcmp(lat,'')
                    fprintf(ODVfid,'%s\t%d\t%s\t%s\t%13.9f\t%13.9f\t%s\t',cruise, station, type, strcat(date,'T',timeStamp), long, lat, depth);
                    for i=1:numPRRTabs;
                        fprintf(ODVfid,'\t');
                    end
                    fprintf(ODVfid,'%8.5f\t%8.5f\t%8.5f\t%8.5f\t%8.5f\t%8.5f\t%8.5f\t%8.5f\t%s\t%s\r\n', ...
                        ch1, ch2, ch3, ch4, ch5, ch6, ch7, ch8, PARGainErr, PARStale);
                    station = station + 1;
                end
        end
    end
    fclose(LVfid);
end
fclose(ODVfid);
toc