glu=[datenum(string(readcell('Clarity_Export_GRATA_PAMELA_2024-01-21_144750',...
    'Range','B3503:B3611'))),readmatrix('Clarity_Export_GRATA_PAMELA_2024-01-21_144750','Range','H3503:H3611')];
fs=[datenum([2024 01 20 12 57 00;...
    2024 01 20 13 32 00;...
    2024 01 20 13 53 00;...
    2024 01 20 14 14 00;...
    2024 01 20 14 37 00;...
    2024 01 20 15 00 00;...
    2024 01 20 15 22 00;...
    2024 01 20 15 45 00;...
    2024 01 20 16 08 00;...
    2024 01 20 16 29 00;...
    2024 01 20 16 52 00;...
    2024 01 20 17 14 00;...
    2024 01 20 17 38 00;...
    2024 01 20 18 02 00;...
    2024 01 20 18 23 00;...
    2024 01 20 18 45 00;...
    2024 01 20 19 09 00;...
    2024 01 20 19 31 00;...
    2024 01 20 19 52 00]),...
    [180;183;202;189;204;218;256;287;289;281;241;232;220;174;136;118;104;85;95]];

tnom=-30+5.12/2:.0025:-5.12/2;
p=dir('pam*');

for i=1:length(p)
    filets(i)=datenum(str2num(p(i).name(5:8)),...
        str2num(p(i).name(10:11)),...
        str2num(p(i).name(13:14)),...
        str2num(p(i).name(16:17)),...
        str2num(p(i).name(19:20)),...
        str2num(p(i).name(22:23)));
    spects(i,:)=filets(i)+tnom/24/3600;
end

t=spects';
glufit=t(:);
glufit(:,2)=splinapp([glu(:,1)-30/24/60;fs(:,1)],[glu(:,2);fs(:,2)],ones(size(glu,1)+size(fs,1),1)*8e6,glufit(:,1));
f3=splinapp(fs(:,1),fs(:,2),ones(size(fs,1),1)*5e6,fs(:,1));
glufit(:,3)=splinapp(fs(:,1),fs(:,2),ones(size(fs,1),1)*5e6,glufit(:,1));
p2p=max(glufit(:,2))-min(glufit(:,2));
glufit(:,4)=(glufit(:,2)-p2p)*1.1+p2p-11;
plot(glu(:,1)-30/24/60,glu(:,2),'o:',fs(:,1),fs(:,2),'o:',glufit(:,1),glufit(:,4),glufit(:,1),glufit(:,3));grid;legend('Dex','FS','DexFit','FSFit');shg
save glu glu fs filets spects glufit f3
