anchor.diff:  
   subroutine normalize(vector)
   implicit none
  -double precision, intent(inout) :: vector(3)
  +double precision :: vector(3)
   double precision :: magnitude
   magnitude = sqrt(vector(1)*vector(1)+vector(2)*vector(2)+vector(3)*vector(3))
   vector = vector/magnitude
  
positive.code: implicit none
  
  ! Arguments
  integer, intent(in)  :: n, np
  integer, intent(out) :: nrot, ierr
  double precision :: a(np,np), d(np), v(np,np)
  
  ! Local variables
  
negative.code: implicit none
  
  ! Arguments
  integer          :: nvalues
  integer          :: ierr
  double precision :: values(nvalues)
  double precision :: coefvar
  
  ! Local variables
  double precision :: s, m
  