function [T,Ks]=kp_read_fun(file);
% чтение индексов геомагнитной активности Ks из месячного файла в TAB-формате
%
%    [Tp,Ks]=kp_read_fun(file);
%
% file - имя файла в виде path_data\kp0801.tab
%
% T     - сериальная дата (часы = 1,4,7,10,13,16,19,22)
% Ks    - индекс

% 27.08.2010    - С.Ю.Хомутов, Геофизическая обсерватория "Ключи"
%                   Новосибирск

NL=length(file);
year0=2000+str2num(file(NL-7:NL-6));
mon0=str2num(file(NL-5:NL-4));
nday=datenum(year0,mon0+1,1)-datenum(year0,mon0,1);

Ks_array=ones(nday,8)*NaN;

fid=fopen(file,'r');
for kd=1:nday,
    str=fgetl(fid);
    year=2000+str2num(str(1:2));
    mon=str2num(str(3:4));
    day=str2num(str(5:6));
    R=str(8:32);
    for ki=1:8,
        [ss,R]=strtok(R);
        switch ss(2),
            case '-',
                fr=-0.33;
            case 'o',
                fr=0;
            case '+',
                fr=+0.33;
        end
        Ks_array(kd,ki)=str2num(ss(1))+fr;
    end
end
Ks=Ks_array'; Ks=Ks(:);
T=datenum(year0,mon0,1,(1:3:24*nday)',0,0);
