files=dir('st1_AA-HED187h.dat');
N=numel(files);
% group_id, azimuth, zenith, ED, Li, Lt
% id time Ed Li Lt
st1_AA=NaN;
st1_AB=NaN;
% AB=11  AB=

for k=1:N
HED=files(k).name; %  Ed dark
base=HPE(1:7);  %
HSE=[base,'HSE187h.dat'];  % progfile Ed light
HLD251c=[base,'HLD251c.dat'];  % Dark 251 C
Hsl251c=[base,'Hsl251c.dat'];    % Light 251 C
HLD250c=[base,'HLD250c.dat'];  % Dark 250c
Hsl250c=[base,'Hsl250c.dat'];
k=1
HED=files(k).name; %  Ed dark
base=HPE(1:7);  %
HSE=[base,'HSE187h.dat'];  % progfile Ed light
HLD251c=[base,'HLD251c.dat'];  % Dark 251 C
Hsl251c=[base,'Hsl251c.dat'];    % Light 251 C
HLD250c=[base,'HLD250c.dat'];  % Dark 250c
Hsl250c=[base,'Hsl250c.dat'];
HED=files(k).name; %  Ed dark
base=HED(1:7);  %
HSE=[base,'HSE187h.dat'];  % progfile Ed light
HLD251c=[base,'HLD251c.dat'];  % Dark 251 C
Hsl251c=[base,'Hsl251c.dat'];    % Light 251 C
HLD250c=[base,'HLD250c.dat'];  % Dark 250c
Hsl250c=[base,'Hsl250c.dat'];
[data,p187_time,SN,INTTIME] = readOCR_2(HED,HSE);
Ed_187 = data;
[data,p250_time,SN,INTTIME] = readOCR_2(HLD250c,Hsl250c);
Lu_250 = data;
[data,p251_time,SN,INTTIME] = readOCR_2(HLD251c,Hsl251c);
Lu_251 = data;

end

save AA_hypro.mat Ed_187 Lu_250 Lu_251 p*