anchor.diff:        end do
               
         do i = 1,p%numBladeNodes ! Loop through the blade nodes / elements
  -         
  -         NodeTxt = '(node '//trim(num2lstr(i))//', blade '//trim(num2lstr(j))//')'
            
               ! Set the active blade element for UnsteadyAero
            m%UA%iBladeNode = i
  
positive.code:          if (p%DBEMT_Mod /= DBEMT_none) then
              do j = 1,p%numBlades ! Loop through all blades
                 do i = 1,p%numBladeNodes ! Loop through the blade nodes / elements
                    NodeTxt = '(node '//trim(num2lstr(i))//', blade '//trim(num2lstr(j))//')'
  
                    call calculate_Inductions_from_BEMT(i, j, p, y%phi(i,j), u, OtherState, AFInfo, y%axInduction(i,j), y%tanInduction(i,j), ErrStat2, ErrMsg2)
                       call SetErrStat(ErrStat2,ErrMsg2,ErrStat,ErrMsg,RoutineName//trim(NodeTxt))
                       if (errStat >= AbortErrLev) return 
  
                 end do
              end do
  
negative.code:                   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)
                          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
                    
  