anchor.diff:                                 p%Stif0_GL,p%Mass0_GL,p%gravity,u_interp,       &
                                  p%damp_flag,p%beta,                             &
                                  p%node_elem,p%dof_node,p%elem_total,p%dof_total,&
  -                               p%node_total,p%niter,p%ngp,p%coef, ErrStat2, ErrMsg2)
  +                               p%node_total,p%niter,p%tol,p%ngp,p%coef, ErrStat2, ErrMsg2)
         call SetErrStat(ErrStat2,ErrMsg2,ErrStat,ErrMsg,RoutineName)
   
      call cleanup()
  
positive.code: 
     END SUBROUTINE BD_GenerateDynamicElementGA2
  
     SUBROUTINE BD_InputGlobalLocal( p, u, flag, ErrStat, ErrMsg)
  
     TYPE(BD_ParameterType), INTENT(IN   ):: p
     TYPE(BD_InputType),     INTENT(INOUT):: u
     INTEGER(IntKi),         INTENT(IN   ):: flag            ! 0: Global to Blade;
     INTEGER(IntKi),         INTENT(  OUT):: ErrStat       ! Error status of the operation
     CHARACTER(*),           INTENT(  OUT):: ErrMsg        ! Error message if ErrStat /= ErrID_None
                                                             ! 1: Blade to Global
  
negative.code:          
     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, ErrStat2,ErrMsg2)
        call SetErrStat(ErrStat2,ErrMsg2,ErrStat,ErrMsg,RoutineName)
     ! find x at t+dt
     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)
  