function deglich(r,t)
for h=1:6
for lr=1:4
lrmean_r=mean(r((h-1)*48+(lr-1)*12+1:(h-1)*48+(lr-1)*12+12));
lrmean_t=mean(t((h-1)*48+(lr-1)*12+1:(h-1)*48+(lr-1)*12+12));
siglr_r=std(r((h-1)*48+(lr-1)*12+1:(h-1)*48+(lr-1)*12+12));
siglr_t=std(t((h-1)*48+(lr-1)*12+1:(h-1)*48+(lr-1)*12+12));
for int=1:12
Chauvenet(int)=max(abs(r((h-1)*48+(lr-1)*12+int)-lrmean_r)/siglr_r,
abs(t((h-1)*48+(lr-1)*12+int)-lrmean_t)/siglr_t);
end
z=find(Chauvenet(:)<2);
nint=max(size(z));
R(h,lr,1:nint)=r((h-1)*48+(lr-1)*12+z);
T(h,lr,1:nint)=t((h-1)*48+(lr-1)*12+z);
R(h,lr,13)=nint;
T(h,lr,13)=0;
end
end
save
R
save T