anchor.diff:  logical :: getPointsOutsideBoundary
   logical :: getPointsOnBoundary
   double precision :: epsilon
  +logical :: printRecord
   
   end module clip
   
  
positive.code:         i = i + 1
          call get_command_argument(i,tag)
          read(tag,*) epsilon
  
      else
          call usage('clip: no option '//trim(tag))
  
negative.code:                 getPointsInsideBoundary, &
                  getPointsOutsideBoundary, &
                  getPointsOnBoundary, &
                  epsilon, &
                  printRecord
  
  implicit none
  
  