anchor.diff:         integer, intent(out) :: underflow, ierr
        end subroutine zairy
   
  +     ! Airy function Bi (or its derivative)
  +     pure subroutine zbiry(re, im, deriv, scaling, zOut_r, zOut_i, underflow, ierr)
  +       implicit none
  +       double precision, intent(in) :: re, im
  +       integer, intent(in) :: deriv, scaling
  +       double precision, intent(out) :: zOut_r, zOut_i
  +       integer, intent(out) :: underflow, ierr
  +     end subroutine zbiry
  +
     end interface
   
   end module amos
  
positive.code: 
       ! Bessel Y function of second kind:
       pure subroutine zbesy(re, im, order, scaling, length, zOut_r, zOut_i, &
          underflows, workr, worki, ierr)
         implicit none
         double precision, intent(in) :: re, im, order
         integer, intent(in) :: scaling, length
  
negative.code:   ! Technically, there is one exception: the ?1mach.f90 routines (used in
    ! finding machine constants) each contain a write statement in case the query
    ! fails. We ignore this possibility (as we trust the queries made by amos).
    ! The author of the wrapper uses this without any problems in production code.
  
    interface
  
  