anchor.diff:     call BEMT_UpdateStates(t, n, m%BEMT_u(1), m%BEMT_u(2),  p%BEMT, x%BEMT, xd%BEMT, z%BEMT, OtherState%BEMT, p%AFI%AFInfo, m%BEMT, errStat2, errMsg2)
         call SetErrStat(ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName)
            
  +      ! Call AeroAcoustics updates states
  +   if ( p%CompAA ) then
  +      ! We need the outputs from BEMT as inputs to AeroAcoustics module
  +      ! Also,  SetInputs() [called above] calls SetInputsForBEMT() which in turn establishes current versions of the Global to local transformations we need as inputs to AA
  +      call SetInputsForAA(p, u(1), m, errStat2, errMsg2)  
  +         call SetErrStat(ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName)
  +      call AA_UpdateStates(t,  n,m%AA, m%AA_u, p%AA, xd%AA,  errStat2, errMsg2)
  +         call SetErrStat(ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName)
  +   end if       
  +   
              
      call Cleanup()
      
  
positive.code:       call SetErrStat(ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName)
                    
     call SetOutputsFromBEMT(p, m, y )
                            
     if ( p%TwrAero ) then
        call ADTwr_CalcOutput(p, u, m, y, ErrStat2, ErrMsg2 )
           call SetErrStat(ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName)      
  
negative.code:    use BEMT
     use AirfoilInfo
     use NWTC_LAPACK
     use AeroAcoustics
     
     
     implicit none
  