anchor.diff:  
        ! Bessel Y function of second kind:
        pure subroutine zbesy(re, im, order, scaling, length, zOut_r, zOut_i, &
  -        underflows, workr, worki, ierr)
  +          underflows, workr, worki, ierr)
          implicit none
          double precision, intent(in) :: re, im, order
          integer, intent(in) :: scaling, length
  
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:        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
  