anchor.diff:        ! --- Special case for floating with extramoment
         if (p%ExtraMoment.and.p%Floating) then
            Y1_CB_L = - (matmul(p%D1_141, m%F_L)) ! Uses rotated loads
  -         !
  -         !Y1_CB_L = matmul(RRb2g, Y1_CB_L) !>>> New
         endif
   
         ! Compute external force on internal (F_L) and interface nodes (F_I)
  
positive.code:       DD(1:nY,7:12  ) = - p%CBB
        DD(1:nY,13:18 ) = - p%MBB
        if (p%nDOFM>0) then
           DD(1:nY,13:18 ) = DD(1:nY,13:18 )+ p%MBmmB
           ! TODO TODO rotate it A MBmmB A^t
        endif
     endif
     
  
negative.code:       REAL(ReKi)                   :: Y1_Guy_R(6)
        REAL(ReKi)                   :: Y1_Guy_L(6)
        REAL(ReKi)                   :: Y1_Utp(6)
        REAL(ReKi)                   :: udotdot_TP(6)
        REAL(ReKi)                   :: Y1_ExtraMoment(3) ! Lever arm moment contributions due to interface displacement
        INTEGER(IntKi), pointer      :: DOFList(:)
        INTEGER(IntKi)               :: startDOF
  