anchor.diff:           enddo
         endif 
      else
  -      m%UR_bar        = 0.0_ReKi
  -      m%UL            = 0.0_ReKi
         ! --- CB modes contribution to motion (L-DOF only), NO STATIC IMPROVEMENT
         if (bElastic .and. p%nDOFM > 0) then
            m%UL = matmul( p%PhiM,  x%qm    )
  +      else
  +         m%UL = 0.0_ReKi
         end if
         ! --- Adding Guyan contribution to R and L DOFs
         if (bGuyan .and. .not.p%Floating) then
  -         m%UR_bar        =                       matmul( p%TI      , m%u_TP       )
  -         m%UL            =   m%UL            +   matmul( p%PhiRb_TI, m%u_TP       ) 
  +         m%UR_bar =         matmul( p%TI      , m%u_TP       )
  +         m%UL     = m%UL +  matmul( p%PhiRb_TI, m%u_TP       ) 
         else
            ! Guyan modes are rigid body modes, we will add them in the "Full system" later
  +         m%UR_bar = 0.0_ReKi
         endif
         ! --- Build original DOF vectors (DOF before the CB reduction)
  -      m%U_red       (p%IDI__) = m%UR_bar
  -      m%U_red       (p%ID__L) = m%UL     
  -      m%U_red       (p%IDC_Rb)= 0    ! NOTE: for now we don't have leader DOF at "C" (bottom)
  -      m%U_red       (p%ID__F) = 0
  +      m%U_red(p%IDI__) = m%UR_bar
  +      m%U_red(p%ID__L) = m%UL     
  +      m%U_red(p%IDC_Rb)= 0    ! NOTE: for now we don't have leader DOF at "C" (bottom)
  +      m%U_red(p%ID__F) = 0
         if (p%reduced) then
  -         DU_full        = matmul(p%T_red, m%U_red)
  +         DU_full = matmul(p%T_red, m%U_red)
         else
  -         DU_full        = m%U_red
  +         DU_full = m%U_red
         endif
         ! --- Adding Guyan contribution for rigid body
         if (bGuyan .and. p%Floating) then
  
positive.code:       endif
  
        ! Compute external force on internal (F_L) and interface nodes (F_I)
        call GetExtForceOnInternalDOF(u, p, x, m, m%F_L, ErrStat2, ErrMsg2, ExtraMoment=(p%ExtraMoment), RotateLoads=.False., U_full=m%U_full); if(Failed()) return
        call GetExtForceOnInterfaceDOF(p, m%Fext, F_I)
  
        ! Compute reaction/coupling force at TP
  
negative.code:             rIP0(1:3)   = p%DP0(1:3, iSDNode)
              rIP(1:3)    = matmul(Rb2g, rIP0)
              duP(1:3)    = rIP - rIP0 ! NOTE: without m%u_TP(1:3)
              ! Full diplacements Guyan + rotated CB (if asked)
              DU_full(DOFList(1:3)) = matmul(Rb2g, DU_full(DOFList(1:3))) + duP(1:3)       
              DU_full(DOFList(4:6)) = matmul(Rb2g, DU_full(DOFList(4:6))) + rotations(1:3)
           enddo
  