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:    real(ReKi)          ,   intent(out)    :: F_L(p%nDOF__L)  !< External force on internal nodes "L"
     integer(IntKi),         intent(  out)  :: ErrStat     !< Error status of the operation
     character(*),           intent(  out)  :: ErrMsg      !< Error message if ErrStat /= ErrID_None
     integer :: iNode ! indices of u-mesh nodes and SD nodes
     integer :: nMembers
     integer :: startDOF, I
  
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
  