k=0;
for ii=is(2):24*29:is(8)-24*29;
  kk=ii:ii+24*29;
  u1=spdllp(kk);
  u2=abs(wllp(kk));
  k=k+1;
  jdk(k)=jd(ii+15*24);
  if(length(find([u1+u2]~=NaN)>1)),
  [mean1, std1, mean2, std2, cor(k), transfn(k)]=stats(u1,u2,1);
  else
     cor(k)=NaN;
     transfn(k)=NaN;
  end
end
