function Instmeta = readAQDHeader(aqdFile)
%function to get instrument metadata about AWAC deployment from header file

hdrFile = strcat(aqdFile,'.hdr');
hdr = fopen(hdrFile);
str='User setup';
while (~strncmp(str,'Hardware configuration',22));
     str=fgetl(hdr);
         if (strncmp(str,'Measurement interval',20))
             is=findstr(str,'sec');
             Instmeta.AQDMeasurementInterval = str2num(str(39:is-2));
         elseif (strncmp(str,'Average interval',15))
             is=findstr(str,'sec');
             Instmeta.AQDAverageInterval = str2num(str(39:is-2));
          elseif (strncmp(str,'Measurement Load',20))
             is=findstr(str,'%');
             Instmeta.AQDMeasurementLoad = str2num(str(39:is-2));
         elseif (strncmp(str,'Transmit pulse length',17))
             is=findstr(str,'m');
             Instmeta.AQDTransmitPulseLength = str2num(str(39:42));
         elseif (strncmp(str,'Blanking distance',17))
             is=findstr(str,'m');
             Instmeta.AQDBlankingDistance = str2num(str(39:is-2));
         elseif (strncmp(str,'Diagnostics measurements',17))
             Instmeta.AQDDiagMeas = (str(39:end));
         elseif (strncmp(str,'Diagnostics - Interval',20))
             is=findstr(str,'sec');
             Instmeta.AQDDiagInterval = str2num(str(39:is-2));
         elseif (strncmp(str,'Diagnostics - Number of samples',26))
             Instmeta.AQDDiagNumberOfSamples = str2num(str(39:end));
         elseif (strncmp(str,'Diagnostics - Number of pings',26))
             Instmeta.AQDDiagNumberOfPings = str2num(str(39:end));
         elseif (strncmp(str,'Coordinate system',15))
             Instmeta.AQDCoordinateSystem = str(39:end);
         elseif (strncmp(str,'Software version',15))
             Instmeta.AQDSoftwareVersion = num2str(str(39:42));
         elseif (strncmp(str,'Salinity',8))
             is=findstr(str,'ppt');
             Instmeta.AQDAssumedSalinity = str2num(str(39:is-2));
         elseif (strncmp(str,'Deployment name',15))
             Instmeta.AQDDeploymentName = str(39:44);
         elseif (strncmp(str,'Deployment time',15))
             Instmeta.AQDDeploymentTime = str(39:59);
         elseif (strncmp(str,'Comments',8))
             Instmeta.AQDComment = str(39:end);
         end
end
while (~strncmp(str,'Head configuration',15));
     str=fgetl(hdr);
         if (strncmp(str,'Serial number',13))
             Instmeta.AQDInstrumentSerialNumber = str(39:46);
         end
end

while (~strncmp(str,'Data file format',16));
     str=fgetl(hdr);
         if (strncmp(str,'Head frequency',12))
             is=findstr(str,'kHz');
             Instmeta.AQDFrequency = str2num(str(39:is-2));
         elseif (strncmp(str,'Serial number',13))
             Instmeta.AQDHeadSerialNumber = str(39:46);
         elseif (strncmp(str,'Transformation Matrix',15))
             Instmeta.AQDTransformationMatrix = zeros(3,3);
             Instmeta.AQDTransformationMatrix(1,:) = strread(str(39:end));
             str=fgetl(hdr);
             Instmeta.AQDTransformationMatrix(2,:) = strread(str(39:end));
             str=fgetl(hdr);
             Instmeta.AQDTransformationMatrix(3,:) = strread(str(39:end));
         end
end

fclose(hdr)