anchor.diff:     CHARACTER(*),           INTENT(  OUT) :: ErrMsg         ! Error message if ErrStat /= ErrID_None
      
      !locals
  -   INTEGER                               :: I, J, Jj       ! counters into full or reduced matrix
  +   INTEGER                               :: I, J           ! counters into full or reduced matrix
      INTEGER                               :: L              ! number of DOFs to eliminate 
  +
      INTEGER, ALLOCATABLE                  :: idx(:)         ! vector to map reduced matrix to full matrix
      INTEGER                               :: NReactDOFs
      INTEGER                               :: DOF_reduced
  
positive.code:       ErrMsg = 'ReduceKMdofs:invalid matrix sizes.'
        RETURN
     END IF
     
     !........
     ! Calculate how many rows/columns need to be eliminated:
     !........
     L = 0
     DO I = 1, NReactDOFs  !Cycle on reaction DOFs
        IF (Init%BCs(I, 2) == 1) THEN
           L=L+1 !number of DOFs to eliminate
        END IF    
     END DO   
     
  
negative.code:    !       this will add additional matrices to the SubDyn summary file.
     !............................
  
     TYPE(ProgDesc), PARAMETER  :: SD_ProgDesc = ProgDesc( 'SubDyn', 'v1.01.00b-bjj', '3-Sept-2014' )
        
     ! ..... Public Subroutines ...................................................................................................
  
  