anchor.diff:         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
  +   !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
  +   c__=0; 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 
  +          c__=c__+1
  +          p%IDC__(c__) = p%NodesDOFtilde(iNode)%List(J) ! DOF number 
             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 
  
positive.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          ! counters
     integer(IntKi)              :: iNode, iiNode
     integer(IntKi)              :: nNodes_R
     integer(IntKi), allocatable :: IDAll(:)
  
negative.code:    ! --------------------------------------------------------------------------------
  
     ! Distribute the interface DOFs into R,F
     c__=0; c_B=0;  c_F=0 ! Counters over R and F dofs
     do iiNode= 1,p%nNodes_I !Loop on interface nodes
        iNode = p%Nodes_I(iiNode,1)
        do J = 1, 6 ! DOFs: ItfTDXss    ItfTDYss    ItfTDZss    ItfRDXss    ItfRDYss    ItfRDZss
            c__=c__+1
            p%IDI__(c__) = p%NodesDOFtilde(iNode)%List(J) ! DOF number 
            if (p%Nodes_I(iiNode, J+1)==idBC_Leader) then
               c_B=c_B+1
               p%IDI_Rb(c_B) = p%NodesDOFtilde(iNode)%List(J) ! DOF number 
  