anchor.diff:        do n1=1,15
            
            ! perturb x^S:
  -         call getRandomVector(LinVec_1)
  -         call getRandomVector(LinVec_1_a)
  -         call getRandomVector(LinVec_1_b)
  -         
  -         LinVec_2_approx = matmul( map_mod1_mod2%dm%tv, LinVec_1 ) + matmul( map_mod1_mod2%dm%mi, LinVec_1_a ) &
  +         call random_number( LinVec_1 )
  +         LinVec_1 = (2*LinVec_1-1)/real(n1**3)  !note n1 from the loop this is called in
  +         call random_number( LinVec_1_a )
  +         LinVec_1_a = (2*LinVec_1_a-1)/real(n1**3)  !note n1 from the loop this is called in
  +         call random_number( LinVec_1_b )
  +         LinVec_1_b = (2*LinVec_1_b-1)/real(n1**3)  !note n1 from the loop this is called in
  +         
  +         LinVec_2_a = matmul( map_mod1_mod2%dm%tv, LinVec_1 ) + matmul( map_mod1_mod2%dm%mi, LinVec_1_a ) &
                           + matmul( map_mod1_mod2%dm%fx_p, LinVec_1_b )! approximate delta theta^D         
                         
            j=1
  
positive.code:       ! Linearize about the operating points:
        ! ..............................................................................................................................   
        
        call AllocAry( LinVec_1, Mesh1_O%nnodes*3, 'LinVec_1', ErrStat, ErrMsg);       IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
        call AllocAry( LinVec_1_a, Mesh1_O%nnodes*3, 'LinVec_1_a', ErrStat, ErrMsg);       IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
        call AllocAry( LinVec_1_b, Mesh1_O%nnodes*3, 'LinVec_1_b', ErrStat, ErrMsg);       IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
        call AllocAry( LinVec_1_c, Mesh2_O%nnodes*3, 'LinVec_1_c', ErrStat, ErrMsg);       IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
        call AllocAry( LinVec_2, Mesh2_O%nnodes*3, 'LinVec_2', ErrStat, ErrMsg);       IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
        call AllocAry( LinVec_2_approx, Mesh2_O%nnodes*3, 'LinVec_2_approx', ErrStat, ErrMsg);       IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
        
        
        CALL MeshCopy (        SrcMesh          = mesh1_O          &
  
negative.code:             j = j+3
           end do
                    
           print *, 'Rotational Velocity: ', n1, TwoNorm( LinVec_2 - LinVec_2_a )                              
        end do            
        
        
  