function [time,z0b,z0tb,taub3n,taub3e,hsb3,hlb3] = read_bulk(path,period)
load([path 'YOPP_BRW_ARM_flux_data_2018_avg_10_min.mat'])
time = datenum([2018 1 0]) + yd_blk';

Archive = [time z0b' z0tb' taub3n' taub3e' hsb3' hlb3'];
Archive(time < datenum(period(1,:)) | time > datenum(period(2,:)) + 1 - 10/24/60,:) = [];

time = Archive(:,1);
z0b = Archive(:,2);
z0tb = Archive(:,3);
taub3n = Archive(:,4);
taub3e = Archive(:,5);
hsb3 = Archive(:,6);
hlb3 = Archive(:,7);