anchor.diff:         enddo
      enddo
      ! Indices IDC__ = [IDC_B, IDC_F, IDC_L], interface
  -   call concatenate_3lists(p%IDC_Rb, p%IDC_F, p%IDC_L, p%IDC__, ErrStat2, ErrMsg2); if(Failed()) return
  +   !call concatenate_3lists(p%IDC_Rb, p%IDC_F, p%IDC_L, p%IDC__, ErrStat2, ErrMsg2); if(Failed()) return
  +   !call sort_in_place(p%IDC__)
  +
   
      ! Indices IDR__ = [IDI__, IDC__], interface
      !call concatenate_lists(p%IDI__, p%IDC__, p%IDR__, ErrStat2, ErrMsg2); if(Failed()) return
  
positive.code:        enddo
     enddo
     ! Indices IDI__ = [IDI_B, IDI_F], interface
     call concatenate_lists(p%IDI_Rb, p%IDI_F, p%IDI__, ErrStat2, ErrMsg2); if(Failed()) return
  
     ! Distribute the reaction DOFs into R,F,L 
     c_B=0; c_F=0; c_L=0; ! Counters over R, F, L dofs
     do iiNode= 1,p%nNodes_C !Loop on interface nodes
        iNode = p%Nodes_C(iiNode,1)
        do J = 1, 6 ! DOFs 
            if (p%Nodes_C(iiNode, J+1)==idBC_Leader) then
               c_B=c_B+1
               p%IDC_Rb(c_B) = p%NodesDOFtilde(iNode)%List(J) ! DOF number 
  
negative.code:    integer(IntKi),          intent(  out) :: ErrStat     !< Error status of the operation
     character(*),            intent(  out) :: ErrMsg      !< Error message if ErrStat /= ErrID_None
     ! local variables
     integer(IntKi)              :: I, J, c_B, c_F, c_L, c__          ! counters
     integer(IntKi)              :: iNode, iiNode
     integer(IntKi)              :: nNodes_R
     integer(IntKi), allocatable :: IDAll(:)
  