anchor.diff:  	real,    intent(in)  :: offset, firstphase
   	real,    intent(out) :: phases(ns, np)
   
  -	! Internal
  -	! --------
  -	real, allocatable :: myphases(:,:)
  -
   	! Call to original subroutine
   	! ---------------------------
  -	call sl_phasecomp(ffrec, np, myphases, weights=weights, offset=offset, firstphase=firstphase)
  -
  -	! Recover results
  -	! ---------------
  -	phases = myphases
  +	call sl_phasecomp(ffrec, np, phases, weights=weights, offset=offset, firstphase=firstphase)
   
   end subroutine phasecomp
   
  
positive.code: 	real,    intent(in)  :: weights(ns)
  	integer, intent(in)  :: useteof
  
  	! Internal
  	! --------
  	real, allocatable :: mypc(:,:), myeof(:,:), myev(:)
  
  	! Call to original subroutine
  	! ---------------------------
  	call sl_pca(ff, nkeep=nkeep, xeof=myeof, pc=mypc, ev=myev, weights=weights, useteof=useteof)
  
  	! Recover results
  	! ---------------
  	xeof = myeof
  	pc   = mypc
  	ev   = myev
  
  end subroutine pca
  
  
negative.code: 	real,   intent(in)  :: ff(nchan,nt)
  	real,	  intent(out) :: steof(nchan*nwindow,nkeep), stpc(nt-nwindow+1,nkeep), ev(nkeep)
  
  	! Call to original subroutine
  	! ---------------------------
  	call sl_mssa(ff, nwindow, nkeep, steof=steof, stpc=stpc, ev=ev)
  		
  end subroutine mssa
  
  
  