anchor.diff:                    if (p%DBEMT_Mod == DBEMT_none ) then
                        
                        call calculate_Inductions_from_BEMT(i,j,p,z%phi(i,j),u1,OtherState,AFInfo,m%axInduction(i,j), m%tanInduction(i,j), ErrStat2,ErrMsg2)
  -                        call SetErrStat(ErrStat2,ErrMsg2,ErrStat,ErrMsg,RoutineName//trim(NodeTxt))
  -                        if (errStat >= AbortErrLev) return
  +                        if (ErrStat2 /= ErrID_None) then
  +                           call SetErrStat(ErrStat2,ErrMsg2,ErrStat,ErrMsg,RoutineName//trim(NodeText(i,j)))
  +                           if (errStat >= AbortErrLev) return
  +                        end if
                        
                     else ! use DBEMT
                     
  
positive.code:    integer(IntKi)                                 :: errStat2    ! temporary Error status of the operation
     character(*), parameter                        :: RoutineName = 'BEMT_CalcOutput'
     
     character(20)                                  :: NodeTxt
     
     
  #ifdef UA_OUTS
     integer(IntKi)                 :: k
  
negative.code:          
           do j = 1,p%numBlades ! Loop through all blades
              do i = 1,p%numBladeNodes ! Loop through the blade nodes / elements
                 
                 call BEMT_UnCoupledSolve(z%phi(i,j), p%numBlades, p%kinVisc, AFInfo(p%AFIndx(i,j)), u1%rlocal(i,j), p%chord(i,j), u1%theta(i,j),  &
                             u1%Vx(i,j), u1%Vy(i,j), p%useTanInd, p%useAIDrag, p%useTIDrag, p%useHubLoss, p%useTipLoss, p%hubLossConst(i,j), p%tipLossConst(i,j), &
                             p%maxIndIterations, p%aTol, OtherState%ValidPhi(i,j), m%FirstWarn_Phi, errStat2, errMsg2)
                    if (ErrStat2 /= ErrID_None) then
                       call SetErrStat(ErrStat2,ErrMsg2,ErrStat,ErrMsg,RoutineName//trim(NodeText(i,j)))
                       if (errStat >= AbortErrLev) return
                    end if
              end do
           end do
        else
  