anchor.diff:     coef(8) = beta * tr2 * deltat2
      coef(9) = tr2
   
  -   END SUBROUTINE BD_TiSchmComputeCoefficients
  -
  -   SUBROUTINE BD_BoundaryGA2(x,p,u,t,OtherState,ErrStat,ErrMsg)
  +END SUBROUTINE BD_TiSchmComputeCoefficients
  +
  +SUBROUTINE BD_BoundaryGA2(x,p,u,t,OtherState,ErrStat,ErrMsg)
      !------------------------------------------------------------
      ! This subroutine applies the prescribed boundary conditions
      ! into states and otherstates at the root finite element node
  
positive.code:    ENDDO
  
  
     END SUBROUTINE BD_ComputeIniCoef
  
  SUBROUTINE BD_Static(t,n,u,utimes,p,x,xd,z,OtherState,ErrStat,ErrMsg)
  
  
negative.code: 
   END SUBROUTINE BD_diffmtc
  
  SUBROUTINE BD_ComputeMemberLength(member_total,kp_member,Coef,seg_length,member_length,total_length,&
                                    ErrStat,ErrMsg)
  !----------------------------------------------------------------------------------------
  ! This subroutine computes the segment length, member length, and total length of a beam.
  ! It also computes the ration between the segment/member and total length.
  ! Segment: defined by two adjacent key points
  !----------------------------------------------------------------------------------------
     INTEGER(IntKi),INTENT(IN   ):: member_total
     INTEGER(IntKi),INTENT(IN   ):: kp_member(:)
     REAL(ReKi),    INTENT(IN   ):: Coef(:,:,:)
  