anchor.diff:     INTEGER(IntKi)                   :: ErrStat2    ! Temporary Error status
      CHARACTER(ErrMsgLen)             :: ErrMsg2     ! Temporary Error message
      CHARACTER(*), PARAMETER          :: RoutineName = 'BD_DynamicSolutionGA2'         
  -   
      REAL(ReKi)                       :: errf
      REAL(ReKi)                       :: StifK(dof_total,dof_total)
      REAL(ReKi)                       :: RHS(dof_total)
  
positive.code:    END SUBROUTINE BD_TiSchmComputeCoefficients
  
     SUBROUTINE BD_BoundaryGA2(x,p,u,t,OtherState,ErrStat,ErrMsg)
  
     TYPE(BD_InputType),           INTENT(IN   )  :: u           ! Inputs at t
     REAL(DbKi),                   INTENT(IN   )  :: t           ! Inputs at t
  
negative.code:    Eref = 0.0D0
  
     DO i=1,niter
         StifK = 0.0D0
         RHS = 0.0D0
         MassM = 0.0D0
  