anchor.diff:  	real, allocatable :: cov(:,:)
   	real, allocatable :: wff(:,:), ww(:), zeof(:,:), zff(:,:)
   	real, allocatable :: eig(:)
  -	integer           :: zuseteof = -1, znkeep, znkeepmax=100, i,j
  -        integer c1,c2,c3
  -        logical l1,l2,l3
  +	integer           :: zuseteof = -1, znkeepmax=100, i,j
   
   	! Setups
   	! ======
  
positive.code: 	! Get ST-PCs
  	! ==========
  	if(present(stpc))then
  !!$		if(not(allocated(stpc))) allocate(stpc(nt-nwindow+1, nkeep))
  		allocate(trff(nt, nchan))
  		trff = transpose(zff)
  		do im = 1, nkeep
  
negative.code: 		! Eigenvalues
  		! -----------
  		if(present(ev))then
  			ev = eig(ns:ns-nkeep+1:-1)
  		end if
  
  