function [time,xh,Archive] = read_prsn(path,period)
Archive = [];
time = [];

for k = 1:3
    num = sprintf('%02g',k);
    z = dir([path '2018-' num '\kazr*']);
    for i = 1:length(z)
        [xt,xh,xs] = read_SR([path '2018-' num '\' z(i).name]);
        xtu = unique(xt);
        for j = 1:length(xtu)
            Archive = [Archive xs(xt == xtu(j))];
            time = [time; xtu(j)/24 + datenum(z(i).name(5:12),'yyyymmdd')];
        end
    end
end
fclose all;

% QC
Archive(Archive == -9.9) = NaN;
% convert
Archive = Archive/60/60;

xh = xh(1:111);


