Mosaicing

mospsf needs to estimate the ``average'' beam appropriate for restoring.

  % invert ... beam=xxx.bm options=systemp,double,mosaic imsize=129,129 cell=1,1 
  % imfit in=xxx.bm object=beam

  % mossdi ...
or
  % mosmem ...

  % restor ... fwhm=8,6 pa=40


TODO: needs more explanation

Even for a single pointing observation, your beam (dataset xxx.bm in the example) will currently contain 3 maps (i.e. an image cube). The first plane is mostl likely the OVRO-OVRO beam, followed by the OVRO-BIMA beam, and finally the BIMA=BIMA beam.

It is also important to set the area to be cleaned carefully. Use a mosaic sensitivity map (see the task mossen) and use something like a 1.5$\sigma$ cutoff. The mask thus generated can be copied into the cube to be cleaned.