anchor.diff:        REAL(ReKi)                   :: rotations(3)
         REAL(ReKi)                   :: ULS(p%nDOFL),  UL0m(p%nDOFL),  FLt(p%nDOFL)  ! Temporary values in static improvement method
         REAL(ReKi)                   :: Y1(6)
  +      REAL(ReKi)                   :: Y1_ExtraMoment(3) ! Lever arm moment contributions due to interface displacement
         INTEGER(IntKi), pointer      :: DOFList(:)
         INTEGER(IntKi)               :: startDOF
         REAL(ReKi)                   :: DCM(3,3),junk(6,p%nNodes_L)
  
positive.code:          Y1 = -( matmul(p%KBB,   m%u_TP) + matmul(p%MBB, m%udotdot_TP) + matmul(p%D1_14, m%UFL)   &  ! -(  0*x + D1(1,1)*u(1) + 0*u(2) + D1(1,3)*u(3) + D1(1,4)*u(4)
                   - matmul( HydroForces, p%TI )  + p%FY )                                             !    + D1(1,5)*u(5) + Fy(1) )
        END IF
        
        ! values on the interface mesh are Y1 (SubDyn forces) + Hydrodynamic forces
        y%Y1Mesh%Force (:,1) = Y1(1:3) 
        y%Y1Mesh%Moment(:,1) = Y1(4:6) 
              
       !________________________________________
       ! CALCULATE OUTPUT TO BE WRITTEN TO FILE 
  
negative.code: 
     SUBROUTINE CleanUp()
        INTEGER(IntKi)             :: ErrStat3    ! The error identifier (ErrStat)
        CHARACTER(ErrMsgLen)       :: ErrMsg3     ! The error message (ErrMsg)
        CALL SD_DestroyContState( xdot,     ErrStat3, ErrMsg3 )
        CALL SD_DestroyContState( k1,       ErrStat3, ErrMsg3 )
        CALL SD_DestroyContState( k2,       ErrStat3, ErrMsg3 )
  