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 06 22 10 40 00;...
2024 06 22 10 56 00;...
2024 06 22 11 12 00;...
2024 06 22 11 34 00;...
2024 06 22 11 52 00;...
2024 06 22 12 08 00;...
2024 06 22 12 24 00;...
2024 06 22 12 40 00;...
2024 06 22 12 57 00;...
2024 06 22 13 16 00;...
2024 06 22 13 33 00;...
2024 06 22 13 50 00;...
2024 06 22 14 06 00;...
2024 06 22 14 22 00;...
2024 06 22 14 39 00]),...
[088;091;102;156;225;244;223;221;232;183;165;157;133;099;106],...
[101;095;114;175;210;225;243;222;228;199;178;157;136;117;103]];


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(:,3),ones(size(fs,1),1)*5e7,st);

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

save day1 sitabs sitabsorg