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:        ENDDO
     ENDDO
  
     CALL BD_InputGlobalLocal(p,u_tmp,0,ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
     CALL BD_BoundaryGA2(x_tmp,p,u_tmp,t,OS_tmp,ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  
negative.code:        ENDDO
         CALL BD_UpdateDynamicGA2(ai,uuNf,vvNf,aaNf,xxNf,coef,node_total,dof_node)
         IF(i==1) THEN
             Eref = SQRT(abs(DOT_PRODUCT(ai_temp,feqv)))*tol
             IF(Eref .LE. TOLF) RETURN
         ENDIF
         IF(i .GT. 1) THEN
             Enorm = SQRT(abs(DOT_PRODUCT(ai_temp,feqv)))
             IF(Enorm .LE. Eref) RETURN
         ENDIF
  
         IF(i==niter) THEN
            CALL setErrStat( ErrID_Fatal, "Solution does not converge after the maximum number of iterations", ErrStat, ErrMsg, RoutineName)
  