[daip] IMAGR: multi-resolution with ATCA data

[Olaf Wucknitz]

Hello,

I am currently trying to map some ATCA observations of a source with 
structure on many scales. The multi-resolution option of IMAGR seems to be 
very appropriate in this case. Unfortunately it does not work properly at 
the moment.

The problem is that it only produces one dirty beam (for the point-source 
channel) and does not scale the fluxes of the dirty maps accordingly. All 
parameters in IMAGRPRM that depend on the beam ratio do not have an 
effect, because the beam sizes are assumed to be equal.
The reason for this behaviour seems to be related to the following piece 
of code in IMAGR.FOR (line 667 and following):


C                                       East-West not 3D
      IF ((UVTYPE.NE.'UV') .AND. (DO3DIM)) THEN
         IF (MFIELD.GT.MAXAFL) THEN
            MSGTXT = 'NFIELD > 64 NOT ALLOWED FOR -NCP DATA'
            CALL MSGWRT (7)
            IRET = 8
            GO TO 999
            END IF
         DO3DIM = .FALSE.               <-------- this line !!!!
         DIM(1) = 1
         DIM(2) = 1
         CALL OPUT ('Input', 'DO3DIMAG', OOALOG, DIM, DO3DIM, CDUMMY,
     *      IRET)
         IF (IRET.NE.0) GO TO 999
         CALL OPUT (CLEAN, 'DO3DIMAG', OOALOG, DIM, DO3DIM, CDUMMY,
     *      IRET)
         IF (IRET.NE.0) GO TO 999
         END IF

I understand why DO3DIM is not really needed for east-west arrays, but in 
this case setting it to false seriously confuses the multi-resolution 
clean which requires DO3DIM = true.
In the current case, commenting the marked line out solved the problem for 
me.


Cheers,
Olaf