anchor.diff:              j = j+3
            end do
                     
  -         print *, 'Translational Velocity: ', n1, TwoNorm( LinVec_2 - LinVec_2_approx )
  +         print *, 'Translational Velocity: ', n1, TwoNorm( LinVec_2 - LinVec_2_a )
                                 
         end do      
         
  
positive.code:    REAL(ReKi)              :: position(3)
     REAL(ReKi)              :: Angle
     REAL(ReKi), allocatable :: LinVec_1(:), LinVec_1_a(:), LinVec_1_b(:), LinVec_1_c(:)
     REAL(ReKi), allocatable :: LinVec_2(:), LinVec_2_approx(:) ! data for linearization testing
     
     CHARACTER(*), PARAMETER :: Fmt = '(ES10.3E2)'
     
  
negative.code:    
        call random_number( Vec )
        LinVec_1 = (2*Vec-1)/real(n1**3)  !note n1 from the loop this is called in
        !LinVec_1 = (Vec)/real(n1**3)  !note n1 from the loop this is called in
        
     end subroutine getRandomVector
     ! ..............................................   
  