anchor.diff:                 IDOFOld = m%NodesDOF(iNode)%List(1:6)
               endif
            else
  -            ErrMsg2='Wrong joint type'; ErrStat2=ErrID_Fatal
  +            allocate(IDOFOld(1:len(m%NodesDOF(iNode))))
  +            IDOFOld(:) = m%NodesDOF(iNode)%List(:)
  +            phat = Init%Nodes(iNode, iJointDir:iJointDir+2)
  +            call JointElimination(Init%NodesConnE(iNode,:), JType, phat, Tc, Init, p, ErrStat2, ErrMsg2)
            endif
            if (allocated(Tc)) then
               nc=size(Tc,2) 
               call init_list(m%NodesDOFtilde(iNode), nc, 0, ErrStat2, ErrMsg2)
               m%NodesDOFtilde(iNode)%List(1:nc) = (/ (iprev + i, i=1,nc) /)
               IDOFNew => m%NodesDOFtilde(iNode)%List(1:nc) ! alias to shorten notations
  -            !print*,'N',iNode,'I ',IDOFOld
  -            !print*,'N',iNode,'It',IDOFNew
  +            print*,'N',iNode,'I ',IDOFOld
  +            print*,'N',iNode,'It',IDOFNew
               Tred(IDOFOld, IDOFNew) = Tc
               iPrev = iPrev + nc
            else
  
positive.code:       integer(IntKi) :: iNode
        integer(IntKi) :: ia ! Index on rigid link assembly
        integer(IntKi) :: m  ! Number of elements connected to a joint
        nDOF_ConstraintReduced = 0
  
        ! Rigid assemblies contribution
  
negative.code:       ! Contribution from all the other joints
        do iNode = 1, Init%NNode
           m = Init%NodesConnE(iNode,1) ! Col1: number of elements connected to this joint
           NodeType = Init%Nodes(iNode,iJointType)
  
           if    (NodeType == idJointPin ) then
              nDOF_ConstraintReduced = nDOF_ConstraintReduced + 5 + 1*m
              print*,'Node',iNode, 'is a pin joint, number of members involved: ', m
  
           elseif(NodeType == idJointUniversal ) then
              nDOF_ConstraintReduced = nDOF_ConstraintReduced + 4 + 2*m
              print*,'Node',iNode, 'is an universal joint, number of members involved: ', m
  
           elseif(NodeType == idJointBall ) then
              nDOF_ConstraintReduced = nDOF_ConstraintReduced + 3 + 3*m
              print*,'Node',iNode, 'is a ball joint, number of members involved: ', m
  
           elseif(NodeType == idJointCantilever ) then
              if ( NodeHasRigidElem(iNode, Init, p)) then
                 ! This joint is involved in a rigid link assembly, we skip it (accounted for above)
                 print*,'Node',iNode, 'is involved in a RA'
  