[daip] bug in IMFIT

Lorant Sjouwerman, NRAO lsjouwer at aoc.nrao.edu
Fri Jul 11 11:46:14 EDT 2003


Eric,


There is a bug in the CC-type assignment in IMFIT, probably in all versions 
as I noticed this earlier years ago. But now I know how to report it ;-)


I ran IMFIT with DOMODEL 1:

  DOMODEL.....If true (> 0), put the deconvolved solutions in a new CC
              file attached to the input image.

EURO  > IMFIT2: Writing CC file version   3 attached to input image
EURO  > IMFIT2:   Integral intensity= 4.8692E-01 +/-  2.58E-03 JANSKYS
EURO  > IMFIT2:   Major axis        =  0.00103 +/-0.00000 asec
EURO  > IMFIT2:   Minor axis        =  0.00063 +/-0.00000 asec
EURO  > IMFIT2:   Position angle    =  20.503 +/-  0.276 degrees
EURO  > IMFIT2:                   Deconvolution of component in asec
EURO  > IMFIT2:                        Nominal     Minimum     Maximum
EURO  > IMFIT2:     Major ax          0.00057     0.00056     0.00057
EURO  > IMFIT2:     Minor ax          0.00022     0.00021     0.00023
EURO  > IMFIT2:     Pos ang          55.25546    54.41662    56.10579

>go prtcc
 euro      PRTCC(31DEC03)    146     10-JUL-2003  17:31:33    Page    1
SGRA        .LCL001.   4  Disk= 2    CC table version   3
Title: AIPS CC CLEAN COMPONENTS TABLE
Created by      IMFIT on 10-JUL-2003 17:30:41
Last written by IMFIT on 10-JUL-2003 17:30:41
Ncomps       1    sort on: FLUX     (DESCAB)

   Comp   Deltax   Deltay  Majax  Minax    Pa  Type  Flux       Total
      #   MicroS   MicroS  Micro  Micro   Degr        Jy         Jy
      1     2.69     5.37**************   55.3  CGau 0.486918   0.48692
AIPS 2: Resumes

  So it thinks the type is 'CGau', ie convolved (TYPE OBJ CODE = 2)


>ine'cc'
>go prtab
 euro      PRTAB(31DEC03)    146     11-JUL-2003  08:53:07    Page    1
SGRA        .LCL001.   4  Disk= 2    CC Table version   3
Title: AIPS CC CLEAN COMPONENTS TABLE
Created by      IMFIT on 10-JUL-2003 17:30:41
Last written by IMFIT on 10-JUL-2003 17:30:41
Ncol   7  Nrow       1    Sort cols:   1 (DESCAB)
   Table can be written as a FITS ASCII table

COL. NO.      1          2          3          4          5          6
     ROW  FLUX       DELTAX     DELTAY     MAJOR AX   MINOR AX   POSANGLE
  NUMBER  JY         DEGREES    DEGREES    DEGREES    DEGREES    DEGREES
       1  4.869E-01  7.473E-10  1.493E-09  5.671E-04  2.187E-04  5.526E+01
COL. NO.          7
     ROW      TYPE OBJ
  NUMBER      CODE
       1      2.000E+00


but the values are for the DEconvolved type 'Gaus'. Previously I just
edited the column 7 (2->1), but now I thought I should report it


Regards,

Loránt Sjouwerman - VLBA Scientific Services - lsjouwerman at aoc.nrao.edu
-----------------------------------------------------------------------
c/o NRAO Array Operations Center       Phone:  +1-505-835-7332
    P.O. Box 0 (1003 Lopezville Rd)    Switch: +1-505-835-7000
    Socorro NM 87801                   Fax:    +1-505-835-7027




More information about the Daip mailing list