anchor.diff:           ! 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
  -         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
  +!          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
  
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:       endif
  
        ! Extra moment dm = Delta u x (fe + fg)
        if (p%ExtraMoment .or. present(DoExtraMoment)) then
           du = m%DU_full(p%NodesDOF(iNode)%List(1:3)) ! Lever arm
           Moment(1) = Moment(1) + du(2) * Force(3) - du(3) * Force(2)
           Moment(2) = Moment(2) + du(3) * Force(1) - du(1) * Force(3)
  