s=dir('sit*');

sitabs=[];

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

fs=[datenum([2024 06 19 09 41 00;...
2024 06 19 10 05 00;...
2024 06 19 10 26 00;...
2024 06 19 10 48 00;...
2024 06 19 11 13 00;
2024 06 19 11 33 00;...
2024 06 19 11 55 00;...
2024 06 19 12 17 00;...
2024 06 19 12 41 00;...
2024 06 19 13 03 00;...
2024 06 19 13 25 00;...
2024 06 19 13 45 00;...
2024 06 19 14 09 00;...
2024 06 19 14 30 00;...
2024 06 19 14 51 00]),...
[77;119;137;199;217;223;196;216;225;238;220;150;123;139;97]];

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