anchor.diff:     call BD_Input_extrapinterp( u, utimes, u_interp, t+p%dt, ErrStat2, ErrMsg2 )
         call SetErrStat(ErrStat2,ErrMsg2,ErrStat,ErrMsg,RoutineName)
            
  -   CALL BD_InputGlobalLocal(p,u_interp,0,ErrStat2,ErrMsg2)
  +   CALL BD_InputGlobalLocal(p,u_interp,ErrStat2,ErrMsg2)
         call SetErrStat(ErrStat2,ErrMsg2,ErrStat, ErrMsg, RoutineName)
      ! Incorporate boundary conditions
      CALL BD_BoundaryGA2(x,p,u_interp,t,OtherState,ErrStat2,ErrMsg2)
  
positive.code:    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)
  
negative.code:    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
     !------------------------------------------------------------
  
     TYPE(BD_InputType),           INTENT(IN   )  :: u           ! Inputs at t
     REAL(DbKi),                   INTENT(IN   )  :: t           ! Inputs at t
  