%N gaupar1 %F $MIR/src/subs/gaupar.for
%D Determine effective beam of the convolution of two images.
%P rjs
%: image analysis
%B
        subroutine GauPar1(lIn,bmaj2,bmin2,bpa2,
     *                          bunit,bmaj,bmin,bpa,fac)

        integer lIn
        real bmaj2,bmin2,bpa2,bmaj,bmin,bpa,fac
        character bunit*(*)

  Determine the units and effective beam of the convolution of an
  image dataset with a gaussian.

  Input:
    lIn        Handle of the image dataset.
    bmaj2,bmin2 Gaussian beam major and minor FWHM, in same units as
               the input dataset.
    bpa2       Position angle of the gaussian beam, in degrees.
  Output:
    bunit      Units of the result (usually JY/BEAM).
    bmaj,bmin  Effective beam major and minor FWHM.
    bpa        Effective beam position angle, in degrees.
    fac        Scale factor to multiply by to correct units.
