function [nday,X0,Y0,Z0]=ibfv200_read(file);
% Чтение базисных значений длЯ каждых суток из
% стандартного файла Internagnet в формате IBFV2.00
% 
%   [nday,X0,Y0,Z0]=ibfv200_read(file);
%
% file - файл с базисными значениЯми
% nday - порЯдковый день в году
% X0,Y0,Z0 - базисные значениЯ на каждые сутки

fid=fopen(file,'rt');
    if fid==-1, error(['!!! Не могу открыть файл ' file ' !!!']), end

% первая строка ("шапка")
str=fgetl(fid);
year=str2num(str(length(str)-3:length(str)));   % год
Nday=datenum(year+1,1,1)-datenum(year,1,1);     % кол-во дней в году

% ищем начало второй секции (Adopted)
while 1,
    str=fgetl(fid);
    if str(1)=='*', break, end
end

nday=ones(Nday,1)*NaN; X0=nday; Y0=nday; Z0=nday;
for kd=1:Nday,
    str=fgetl(fid);
    if str=='*', disp('!!! Неполный год принятых базисных !!!'); break, end
    if str==-1, error('!!! Некорректное завершение файла базисных !!!'); end
    tmp=str2num(str(1:51));
    nday(kd)=tmp(1);
    X0(kd)=tmp(2); Y0(kd)=tmp(3); Z0(kd)=tmp(4);
end
fclose(fid);
nday=nday(~isnan(nday));
X0=X0(~isnan(nday)); Y0=Y0(~isnan(nday)); Z0=Z0(~isnan(nday));
