%N pbint %F $MIR/src/subs/pb.for
%D Get the value of an integral of a primary beam.
%P mhw
%: image-data
%B
      real function pbint(pbObj,x,y)

      integer pbObj
      real x,y
  ---------------------------------------------------------------------
  Determine the value of the integral (wrt frequency) of the
  primary beam at a given location.

  Input:
    pbObj      Handle of the primary beam object.
    x,y        Grid location of the position of interest.
  Output:
    pbInt      Value of the primary beam.
