anchor.diff:          call usage(     'Offending line: '//trim(line))
       endif
   
  -    call pnpoly(x,y,polygon(:,1),polygon(:,2),nclip,inout,epsilon)
  +    do ipoly = 1,npoly
  +        call pnpoly(x,y,polygon(ipoly,:,1),polygon(ipoly,:,2),nclip,inout,epsilon)
  +        if (inout.eq.1) then
  +            exit
  +        endif
  +    enddo
   
       if (getPointsInsideBoundary.and.inout.eq.1) then
           if (printRecord) then
  
positive.code: implicit none
  
  ! Local variables
  integer :: i, ios, luin, luout, nclip, inout, irecord
  character(len=512) :: line
  double precision :: x, y
  double precision, allocatable :: polygon(:,:)
  
  
  ! Parse command line
  
negative.code:         endif
      elseif (labelAll) then
          if (printRecord) then
              if (multiPoly) then
                  write(luout,*) irecord,inout,ipoly
              else
                  write(luout,*) irecord,inout
              endif
          else
              if (multiPoly) then
                  write(luout,*) inout,ipoly
              else
                  write(luout,*) inout
              endif
          endif
      endif
  enddo
  