anchor.diff:    ! 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).
  -
  -  ! TODO: note on the *mach.f90 routines
  +  ! The author of the wrapper uses this without any problems in production code.
   
     interface
   
  
positive.code:        integer, intent(out) :: underflows, ierr
       end subroutine zbesh
  
       ! Airy function (or its derivative):
       pure subroutine zairy(re, im, deriv, scaling, zOut_r, zOut_i, &
            underflow, ierr)
         implicit none
         double precision, intent(in) :: re, im
         integer, intent(in) :: deriv, scaling
  
negative.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
  