anchor.diff:  !> Construct force vector on internal DOF (L) from the values on the input mesh 
   !! First, the full vector of external forces is built on the non-reduced DOF
   !! Then, the vector is reduced using the Tred matrix
  -SUBROUTINE GetExtForceOnInternalDOF( u, p, m, rotateLoads, F_L, ErrStat, ErrMsg )
  +SUBROUTINE GetExtForceOnInternalDOF( u, p, m, rotateLoads, F_L, ErrStat, ErrMsg, DoExtraMoment )
      type(SD_InputType),     intent(in   )  :: u ! Inputs
      type(SD_ParameterType), intent(in   )  :: p ! Parameters
      logical               , intent(in   )  :: rotateLoads ! If true, loads are rotated to body coordinate 
  
positive.code:          Y1 = -(   matmul(p%KBB,   m%u_TP) + matmul(p%D1_12, m%udot_TP) + matmul(p%D1_13, m%udotdot_TP) &
                   + matmul(p%D1_141, m%F_L) + matmul(p%D1_142, m%F_L)  - matmul( F_I, p%TI ) ) 
        end if
        print*,'-----------------------'
        print*,'t        ',t
        print*,'UTP      ',m%u_TP
        print*,'-----------------------'
        print*,'Y1 old2  ',Y1
  
        Y1_Utp  = - (matmul(p%KBB,   m%u_TP) + matmul(p%D1_12, m%udot_TP) + matmul(p%D1_13, m%udotdot_TP) )
        Y1_CB_L = - (matmul(p%D1_141, m%F_L))
        if ( p%nDOFM > 0) then
           Y1_CB = -( matmul(p%C1_11, x%qm)   + matmul(p%C1_12,x%qmdot))
        else
           Y1_CB = 0.0_ReKi
        endif
        Y1_Guy_R =   matmul( F_I, p%TI )
        Y1_Guy_L = - matmul(p%D1_142, m%F_L) 
        Y1 = Y1_CB + Y1_Utp + Y1_CB_L+ Y1_Guy_L + Y1_Guy_R 
        print*,'-----------------------'
        print*,'Split terms (old)'
  
negative.code:       Rg2b(1:3,1:3) = u%TPMesh%Orientation(:,:,1)  ! global 2 body coordinates
     endif
  
     if (p%ExtraMoment .or. present(DoExtraMoment)) then
        if (p%Floating) then
           ! TODO deactivate for floating 
           ! For fully floating case, we prescribe the Guyan motion as a "rigid" (non-linear) motion
  