anchor.diff:           
      CALL BD_TiSchmPredictorStep( x_tmp%q,x_tmp%dqdt,OS_tmp%acc,OS_tmp%xcc,             &
                                   p%coef,p%dt,x%q,x%dqdt,OtherState%acc,OtherState%xcc, &
  -                                p%node_total,p%dof_node )
  +                                p%node_total,p%dof_node, ErrStat2,ErrMsg2)
  +      call SetErrStat(ErrStat2,ErrMsg2,ErrStat,ErrMsg,RoutineName)
      ! find x at t+dt
  -   CALL BD_InputGlobalLocal(p,u_interp,0,ErrStat2,ErrMsg2)
  +   CALL BD_InputGlobalLocal(p,u_interp,ErrStat2,ErrMsg2)
         call SetErrStat(ErrStat2,ErrMsg2,ErrStat,ErrMsg,RoutineName)
      
      CALL BD_BoundaryGA2(x,p,u_interp,t+p%dt,OtherState,ErrStat2,ErrMsg2)
  
positive.code:    REAL(ReKi)                       :: feqv(dof_total-6)
     REAL(ReKi)                       :: Eref
     REAL(ReKi)                       :: Enorm
     REAL(ReKi),PARAMETER             :: TOLF = 1.0D-05
     REAL(ReKi)                       :: d
     INTEGER(IntKi)                   :: indx(dof_total-6)
     INTEGER(IntKi)                   :: i
  
negative.code:    INTEGER(IntKi),     INTENT(IN   ):: node_total
     INTEGER(IntKi),     INTENT(IN   ):: ngp
     INTEGER(IntKi),     INTENT(IN   ):: niter
     REAL(ReKi),         INTENT(IN   ):: tol
     REAL(ReKi),         INTENT(INOUT):: uuNf(:)
     REAL(ReKi),         INTENT(INOUT):: vvNf(:)
     REAL(ReKi),         INTENT(INOUT):: aaNf(:)
     REAL(ReKi),         INTENT(INOUT):: xxNf(:)
     INTEGER(IntKi),     INTENT(  OUT):: ErrStat     ! Error status of the operation
     CHARACTER(*),       INTENT(  OUT):: ErrMsg      ! Error message if ErrStat /= ErrID_None
  
  