anchor.diff:  	cov = cov / float(nt)
   	if(.not.present(lpc)) deallocate(zll,zls)
   	if(.not.present(rpc)) deallocate(zrr,zrs)
  -	print*, 'SVD cov(:2,:2):',cov(1:2,1:2)
   
   	! SVD
   	! ---
   	allocate(zleof(nsl,ns), zev(ns))
  -	if(zbLargeMatrix)then
  +	if(zbigmat)then
   		! FIXME: problem with gesdd
   		call gesdd(cov, zev, u=zleof, job='V', info=la_info)
   	else
  
positive.code: 
  	! Use divide/conquer algo?
  	! ------------------------
  	if(.not.present(bLargeMatrix))then
  		zbLargeMatrix = .false.
  	else
  		zbLargeMatrix = bLargeMatrix
  	end if
  
  	! Weights
  
negative.code: 	!	the first phase conincides with the maximmum.
  	!
  	!
  	! :Necessary arguments:
  	!	- varrec: Space-time array
  	!	- np:    Number of requested phases over the 360 degrees cycle [default:8]
  	!
  	! :Optional arguments:
  	!	- weights:    Space array of weights
  	!	- offset:     Minimal normalized amplitude of the index [default:0.]
  	!	- firstphase: Value in degrees of the first phase [default:0.]
  	!
  	! :Dependencies:
  	!	sl_pca
  
  
  