anchor.diff:     Eref = 0.0D0
   
      DO i=1,niter
  -!WRITE(*,*) "N-R Iteration #", i
          StifK = 0.0D0
          RHS = 0.0D0
          MassM = 0.0D0
  
positive.code:    INTEGER(IntKi),     INTENT(IN   ):: node_total
     INTEGER(IntKi),     INTENT(IN   ):: ngp
     INTEGER(IntKi),     INTENT(IN   ):: niter
     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
  
  
negative.code:        ENDDO
     ENDDO
  
     CALL BD_InputGlobalLocal(p,u_tmp,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 )
  