read_HOC.m

function [A, unitstr] = read_HOC(fname)
% [A, unitstr] = read_HOC(fname)
% 
% fname = name of ascii file

%HOCDIR = 'C:\Tamar\Matlab\glass data toolbox\Handbook of Optical Constants\ascii tables\';
%HOCDIR = 'C:\Users\dmarx\OneDrive\Matlab\glass data toolbox\Handbook of Optical Constants\ascii tables\';
pn = regexp(userpath,'\','split');
HOCDIR = fullfile(pn{:}, '\glass data toolbox\Handbook of Optical Constants\ascii tables\');

fid = fopen(fullfile(HOCDIR, fname),'rt');
if fid == -1, error(['cannot open file: ' HOCDIR fname]); end

lcnt = 1;
tline = fgetl(fid);
while ~strncmp(strtrim(tline),'eV',2),
   tline = fgetl(fid);
   lcnt = lcnt + 1;
end

unitstr = textscan(tline,'%s');
unitstr = unitstr{1}; % cell of cells
% Nunits = length(unitstr);

% next line is '-----------'
tline = fgetl(fid);
lcnt = lcnt + 2; 

A = textscan(fid,'%f %f %f %f %f %f','delimiter',':', 'EmptyValue', -Inf);

fclose(fid);

%A = reshape(A{1},Nunits,length(A{1})/Nunits).';


% %A = textread([HOCDIR fname],'','headerlines',lcnt,'delimiter',':');
% A = dlmread([HOCDIR fname],':',lcnt,1);

end