function [Tk,K] = dka_read_fun(file)
% Считывание K-индексов из файлов типа *.dka (INTERMAGNET)
%
%            [Tk,K]=dka_read_fun(file);
%   
%
% Формат файла:
%                           PET
%
%               K - Index Values for 2020     ( K9-limit = 450nT )
%
%   DA-MON-YR  DAY #    1    2    3    4      5    6    7    8       SK
%
%   05-SEP-20   249     1    1    2    3      3    1    1    1       13
%   06-SEP-20   250     1    0    0    1      1    0    0    1        4
%   07-SEP-20   251     1    1    1    1      1    2    0    0        7

% 09.09.2020 - S.Y. Khomutov, Observatory "Paratunka", IKIR FEB RAS, Kamchatka

format_dka='%9c %d %d %d %d %d %d %d %d %d %d';

if ~exist(file)
    error(['!!! Не могу найти файл ... ' file])
end
fid=fopen(file,'rt');
while 1     % код станции
    str=fgetl(fid);
    if ~isempty(strtrim(str)), head.Station=strtrim(str); break, end
end
while 1     % K9
    str=fgetl(fid);
    if ~isempty(strtrim(str))
        head.K9=str2num(str(strfind(str,'=')+1:strfind(str,'nT')-1));
        break
    end
end
while 1     % ищем начало данных
    str=fgetl(fid);
    if ~isempty(strtrim(str)), break, end
end
fgetl(fid); % пропускаем пустую строку
% данные
str=char(fread(fid))';  % считываем остаток бинарно и переводим в строку
C=textscan(str,format_dka);
Ttmp=datenum(C{1},'dd-mmm-yy');   % сериальная дата
doy=C{2};   % day-of_year
tmp=[C{3} C{4} C{5} C{6} C{7} C{8} C{9} C{10}];
SK=C{11};   % суммарный K-индекс
% разворачиваем по времени
tmp=tmp'; K=tmp(:);
th=repmat(1.5:3:24,length(Ttmp),1);
% tmp=Ttmp+th/24;
tmp=repmat(Ttmp,1,8)+th/24;
tmp=tmp'; Tk=tmp(:);
end

