s=dir('sit*');

sitabs=[];

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

fs=[datenum([2024 05 17 15 18 00;...
2024 05 17 15 39 00;...
2024 05 17 16 00 00;...
2024 05 17 16 23 00;...
2024 05 17 16 44 00;...
2024 05 17 17 04 00;...
2024 05 17 17 27 00;...
2024 05 17 17 48 00;...
2024 05 17 18 08 00;...
2024 05 17 18 28 00;...
2024 05 17 18 53 00;...
2024 05 17 19 14 00]),...
[73;95;112;144;108;99;98;77;85;108;110;110]];

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 day1 sitabs