PICARD: take current units into account for calculating total flux

For data in mJy/beam, the total flux is in mJy/beam*arcsec**2 - this
quantity should be divided by the beam area in arcsec**2 to get a total in mJy.

Also make autophotom the default method as is it more reliable than my
home-grown attempt.