s=dir('sit*');

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

fs=[datenum([2024 07 05 11 50 00;...
2024 07 05 11 59 00;...
2024 07 05 12 08 00;...
2024 07 05 12 17 00;...
2024 07 05 12 27 00;...
2024 07 05 12 36 00;...
2024 07 05 12 44 00;...
2024 07 05 12 52 00;...
2024 07 05 13 01 00;...
2024 07 05 13 13 00;...
2024 07 05 13 20 00;...
2024 07 05 13 29 00;...
2024 07 05 13 40 00;...
2024 07 05 13 48 00;...
2024 07 05 13 56 00;...
2024 07 05 14 08 00]),...
[101;089;082;082;115;115;135;138;131;099;088;097;124;115;116;103],...%Oxiline
[106;123;106;100;136;165;178;171;161;132;123;128;149;140;143;146]];%Verio



load(s(i).name,'mtim');
st=ones(6,1)*mtim' + [-55:10:-5]'*ones(1,length(mtim))/(24*60*60);
st=st(:);
sg=splinapp(fs(:,1),fs(:,3),ones(size(fs,1),1)*5e7,st);

sitabs=[sitabs,st,sg];
sitabsorg=[sitabsorg,st,sg];

save day3cb2 sitabs sitabsorg