%N getbeam %F $MIR/src/subs/getbeam.for
%D Get beam from image header.
%P mchw
%: utilities
%B
        subroutine GetBeam(lIn,naxis,bunit,bmaj,bmin,omega,cbof)
        implicit none
        integer lIn,naxis
        character bunit*(*)
        real bmaj,bmin,omega,cbof

  Get beam from image header.

  Inputs:
    lIn        The handle of the Image.
    naxis      Number of image axes.
  Outputs:
    bunit      The pixel units.
    bmaj,bmin  Beam major and minor axes [radians].
    omega      Beam solid angle [steradians]
    cbof       Beam oversampling factor.
