anchor.diff:      ErrStat = ErrID_None
       ErrMsgn = ""
       ! Body of TNO 
  -    epsabs = 1e-10     !absolute accuracy
  -    epsrel = 1e-10     !relative accuracy
       band_ratio = 2.**(1./3.)
       ! Module AirfoilParams
       Mach = real(U  / p%SpdSound)
  
positive.code:     character(ErrMsgLen)    :: ErrMsg2                                                                                   ! temporary Error message
      character(*), parameter :: RoutineName                              = 'TBLTE_TNO'
      REAL(kind=4) :: a,b
      REAL(kind=4) :: epsabs,epsrel
      REAL(kind=4) :: answer
      REAL(kind=4) :: abserr,resabs,resasc
      REAL(kind=4) :: alist (limit),blist (limit),rlist (limit)
      REAL(kind=4) :: elist (limit)
      REAL(kind=4) :: freq(size(p%FreqList))
      REAL(kind=4) :: SPL_press,SPL_suction
      REAL(kind=4) :: band_width,band_ratio
      REAL(kind=4) :: Spectrum
      REAL(ReKi)   :: DBARH
      REAL(kind=4) :: P1,P2,P4
      INTEGER (4)  :: neval
      INTEGER (4)  :: ier
      INTEGER (4)  :: iord (limit)
      INTEGER (4)  :: last
      INTEGER (4)  :: n_freq,i_low,i_hi
      INTEGER (4)  :: i_omega
      real(TNOKi)  :: omega
  
  
negative.code:         band_width = freq(i_omega)*(sqrt(band_ratio)-1./sqrt(band_ratio)) * 4. * pi
          ISSUCTION = .TRUE.
          IF (Cf(1) .GT. 0.) THEN
              CALL solve_qk61(omega,f_int2,a,b,answer,abserr,resabs,resasc)   ! wrapper for qk61 routine from slatec library
              Spectrum = D/(4.*pi*R**2.)*answer
              SPL_suction = 10.*log10(Spectrum*DBARH/2.e-5/2.e-5)
              SPLS(i_omega) = SPL_suction + 10.*log10(band_width)
  