s=dir('sit*');

sitabs=[];

for i=1:length(s)
    load(s(i).name,'pds01');
    sitabs=[sitabs;pds01];
    clear pds01;
end

fs=[datenum([2024 05 18 13 22 00;...
2024 05 18 13 44 00;...
2024 05 18 14 00 00;...
2024 05 18 14 17 00;...
2024 05 18 14 32 00;...
2024 05 18 14 47 00;...
2024 05 18 15 03 00;...
2024 05 18 15 18 00;...
2024 05 18 15 34 00;...
2024 05 18 15 50 00;...
2024 05 18 16 05 00;...
2024 05 18 16 21 00]),...
[99;134;172;157;140;151;149;151;107;98;107;94]];


load(s(i).name,'mtim');
st=ones(56,1)*mtim' + [-59:1:-4]'*ones(1,length(mtim))/(24*60*60);
st=st(:);
sg=splinapp(fs(:,1),fs(:,2),ones(size(fs,1),1)*5e7,st);

sitabs=[sitabs,st,sg];

save day2 sitabs