anchor.diff:     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
  +   logical, intent(in   ) ,optional :: DoExtraMoment
      integer :: iNode ! indices of u-mesh nodes and SD nodes
      integer :: nMembers
      integer :: startDOF, I
  
positive.code: !> 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 )
     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 
  
negative.code:          ! NOTE: this is only needed for the few DOFs below the interface
           !       so this could be greatly optimized with the mapping computed at init only..
           !       TODO...
  !          m%U_red_SIM(:) = 0.0_ReKi
  !          m%U_red_SIM(p%ID__L) = ULS   ! Only the static contribution
  !          !m%U_red_SIM(p%ID__L) = m%UL   ! Only the static contribution
  !          if (p%reduced) then
  !             m%U_full_SIM = matmul(p%T_red, m%U_red_SIM)
  !          else
  !             m%U_full_SIM = m%U_red_SIM
  !          endif
  !          iSDNode = p%Nodes_I2(1) ! Second node of interface element
  !          DOFList => p%NodesDOF(iSDNode)%List
  !          ULS_I2 = m%U_full_SIM(DOFList(1:6))
  !          !print*,'NODE',iSDNode
  !          !print*,'DOFLIST',DOFLIST
  !          !print*,'DOFLIST', p%NodesDOF(p%Elems(1,2))%List(1:6) 
  !          !print*,'DOFLIST', p%NodesDOF(p%Elems(1,3))%List(1:6) 
  !          Y1_Guy_L =  - matmul( p%K_I2 , ULS_I2)
  !          print*,'ULs    ',ULS_I2
  !          print*,'ULe    ',m%U_full_elast(DOFList)
  !          if (any(ULS_I2>2)) then
  !             STOP
  !          endif
  !          print*,'-----------------------'
  !          print*,'Split terms (new)'
  !          print*,'Y1 UtP   ',Y1_Utp
  !          print*,'Y1 GR    ',Y1_Guy_R
  !          print*,'Y1 CB    ',Y1_CB
  !          print*,'Y1 CBL   ',Y1_CB_L
  !          print*,'Y1 GL new',Y1_Guy_L
  !          Y1_Guy_L(1:3) = matmul(Rb2g,Y1_Guy_L(1:3))
  !          Y1_Guy_L(4:6) = matmul(Rb2g,Y1_Guy_L(4:6))
  !          Y1_CB(1:3)    = matmul(Rb2g,Y1_CB(1:3))
  !          Y1_CB(4:6)    = matmul(Rb2g,Y1_CB(4:6))
  !          Y1_CB_L(1:3)  = matmul(Rb2g,Y1_CB_L(1:3))
  !          Y1_CB_L(4:6)  = matmul(Rb2g,Y1_CB_L(4:6))
  !          print*,'Y1 GL rot',Y1_Guy_L
  !          print*,'Y1 CB rot',Y1_CB
  !          print*,'Y1 CBLrot',Y1_CB_L
  
           Y1 = Y1_CB + Y1_Utp + Y1_CB_L+ Y1_Guy_L + Y1_Guy_R 
           print*,'Y1    new',Y1
  