anchor.diff:  implicit none
   
   ! Arguments
  -double precision, intent(in)  :: matrix_in(3,3), angle
  -double precision, intent(out) :: matrix_out(3,3)
  -character(len=*), intent(in)  :: axis
  -integer, intent(out)          :: ierr
  +double precision :: matrix_in(3,3), angle
  +double precision :: matrix_out(3,3)
  +character(len=*) :: axis
  +integer          :: ierr
   
   ! Local variables
   double precision :: rot_matrix(3,3), rot_matrix_trans(3,3), matrix_tmp(3,3), cosa, sina
  
positive.code: implicit none
  
  ! Arguments
  integer, intent(in)          :: nvalues
  integer, intent(out)         :: ierr
  double precision, intent(in) :: values(nvalues)
  double precision             :: stdev
  
  ! Local variables
  integer :: i
  
negative.code: implicit none
  
  ! Arguments
  integer          :: nvalues
  integer          :: ierr
  double precision :: values(nvalues)
  double precision :: coefvar
  
  ! Local variables
  double precision :: s, m
  