anchor.diff:  
      if (ExtraMoment) then
         ! Compute node displacements "DU_full" for lever arm
  -      call LeverArm(u, p, x, m, m%DU_full, bGuyan=.True., bElastic=.True., U_full=U_full)
  +      call LeverArm(u, p, x, m, m%DU_full, bGuyan=.True., bElastic=.False.)!, U_full=U_full)
      endif
   
      ! --- Build vector of external forces (including gravity) (Moment done below)  
  
positive.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)
              !U_full(DOFList(1:3)) = U_full(DOFList(1:3)) + duP(1:3)       
              ! >>> New
              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
  
negative.code:          enddo
        endif 
     else
        ! --- 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       ) 
        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
        if (p%reduced) then
           DU_full = matmul(p%T_red, m%U_red)
        else
           DU_full = m%U_red
        endif
        ! --- Adding Guyan contribution for rigid body
        if (bGuyan .and. p%Floating) then
  