[daip] SCIMG,SCMAP aborts

Larry Maddox maddox at nhn.ou.edu
Mon Apr 25 13:16:05 EDT 2005


when trying to image some data using SCMAP or SCIMG, the task abruptly
aborts when i choose "CONTINUE SELF CAL" option from the TV menu.

this occurs after a little cleaning, or after no cleaning.  my inputs are
as follows:

AIPS 1: SCMAP:  Imaging plus self-calibration
AIPS 1: Adverbs     Values                 Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: INNAME     'M83         '          Input UV data (name)
AIPS 1: INCLASS    'C98_SP'                Input UV data (class)
AIPS 1: INSEQ         2                    Input UV data (seq. #)
AIPS 1: INDISK        3                    Input UV data disk drive #
AIPS 1: OUTNAME    'M83_C98     '          Output image name
AIPS 1: OUTSEQ        0                    Output image seq. no.
AIPS 1: OUTDISK       3                    Output image disk drive #
AIPS 1: OUT2NAME   'M83_SC      '          Output UV file name
AIPS 1: OUT2CLAS   '      '                Output UV file class
AIPS 1: OUT2SEQ       0                    Output UV file seq. #
AIPS 1: OUT2DISK      3                    Output UV file disk drive #
AIPS 1:                                        Imaging info
AIPS 1: CELLSIZE      0.3         0.3      (X,Y) size of grid in asec
AIPS 1: IMSIZE     2048        2048        Image size (power of 2)
AIPS 1: SHIFT         0           0        (X,Y) image shift in asec
AIPS 1:                                    NOT a "3D" shift.
AIPS 1: UVTAPER       0           0        (U,V) Gaussian taper
AIPS 1:                                      units are kilo-lambda
AIPS 1: UVRANGE       0           0        Min & max baseline (klambda)
AIPS 1: UVWTFN     '  '                    UV dist. weight function
AIPS 1:                                      blank => uniform
AIPS 1: UVSIZE        0           0        Array size for doing uniform
AIPS 1:                                    weights.  0 -> actual field
AIPS 1:                                    size.  ??
AIPS 1: ROBUST        0                    Robustness power: -4 -> pure
AIPS 1:                                    uniform weights, 4 => natural
AIPS 1: UVBOX         0                    Additional rows and columns
AIPS 1:                                    used in weighting.
AIPS 1: UVBXFN        1                    Box function type when UVBOX
AIPS 1:                                    > 0.  0 -> 1 round pill box.
AIPS 1: ROTATE        0                    Rotate image CCW from N by
AIPS 1:                                    ROTATE degrees
AIPS 1: ZEROSP     *all 0                  0-spacing fluxes and weights
AIPS 1: XTYPE         5                    Conv. function type in x
AIPS 1:                                      default spheroidal
AIPS 1: YTYPE         5                    Conv. function type in y
AIPS 1:                                      default spheroidal
AIPS 1: XPARM      *all 0                  Conv. function parms for x
AIPS 1: YPARM      *all 0                  Conv. function parms for y
AIPS 1:                                         Cleaning info
AIPS 1: GUARD         0           0        x,y guard band fractional
AIPS 1:                                    radius
AIPS 1: NBOXES        0                    Number of boxes for CLEAN
AIPS 1: CLBOX      *all 0                  Four coordinates for each box
AIPS 1: OBOXFILE   'FITS:M83C98SC.BOX      Output file for final Clean
AIPS 1:
AIPS 1:               '
AIPS 1:                                    boxes
AIPS 1: GAIN          0.1                  CLEAN loop gain
AIPS 1: FLUX          0                    Minimum CLEAN component (Jy)
AIPS 1: MINPATCH     51                    Min. BEAM half-width in AP.
AIPS 1: NITER      10000                   Maximum # of CLEAN components
AIPS 1: BMAJ          0                    FWHM(asec) major axis CLEAN
AIPS 1:                                    restoring beam.
AIPS 1: BMIN          0                    FWHM(asec) minor axis CLEAN
AIPS 1:                                    restoring beam.
AIPS 1: BPA           0                    CLEAN beam position angle
AIPS 1: FACTOR        0                    Speedup factor see HELP
AIPS 1: MAXPIXEL   20000                   Maximum pixels searched in
AIPS 1:                                    each major cycle.
AIPS 1: CMETHOD    '    '                  Modeling method:
AIPS 1:                                    'DFT','GRID','    '
AIPS 1:                                         Self cal info
AIPS 1: NMAPS         4                    Number of self-cal loops.
AIPS 1: SMODEL     *all 0                  Initial model, 1=flux,2=x,3=y
AIPS 1: REFANT        4                    Reference antenna
AIPS 1: SOLINT        0                    Solution interval (min)
AIPS 1: APARM      *all 0                  Control parameters
AIPS 1:                                       1=min. no. antennas
AIPS 1:                                       2 >0 => normalize gain
AIPS 1:                                       3 > 0 => avg. RR,LL
AIPS 1:                                       5 > 0 => avg. IFs.
AIPS 1:                                       6=print level, 1=some
AIPS 1:                                       7=soln SNR cutoff (0=>5)
AIPS 1:                                       8=SIGMA (see HELP SCMAP)
AIPS 1:                                       9=Edit averaging time
AIPS 1:                                      10=Amplitude smoothing time
AIPS 1: SOLTYPE    '    '                  Soln type,'  ','L1','GCON',
AIPS 1:                                       'R', 'L1R', 'GCOR'
AIPS 1: SOLMODE    'A&P '                  Soln. mode: 'A&P','P','P!A',
AIPS 1:                                    'GCON'
AIPS 1: SOLCON        0                    Gain constraint factor.
AIPS 1: ANTWT      *all 0                  Ant. weights (0=>1.0)
AIPS 1: GAINERR    *all 0                  Std. Dev. of antenna gains.
AIPS 1: WTUV          0.1                  Weight outside full weight
AIPS 1:                                    range in self-cal
AIPS 1: WEIGHTIT      0                    Modify data weights function
AIPS 1: DOCALIB      -1                    Restart with previous SN
AIPS 1:                                    = 2 calibrate weights
AIPS 1: GAINUSE       0                    SN table to apply
AIPS 1: DOPOL        -1                    If >0 correct polarization.
AIPS 1: BLVER        -1                    BL table to apply.
AIPS 1: FLAGVER       0                    Flag table version
AIPS 1: DOBAND       -1                    If >0 apply bandpass cal.
AIPS 1:                                    Method used depends on value
AIPS 1:                                    of DOBAND (see HELP file).
AIPS 1: BPVER        -1                    Bandpass table version
AIPS 1: DOTV         -1                    Display residuals on TV ?
AIPS 1:                                    0 -> no, 1 -> after each
AIPS 1:                                    major cycle and each full
AIPS 1:                                    Clean, -1 only after each
AIPS 1:                                    full Clean.  If you can, do
AIPS 1:                                    at least -1.
AIPS 1: DOWEIGHT      1                    Rescale weights by DOWEIGHT
AIPS 1:                                    for editing display only
AIPS 1: DOTWO        -1                    True => do second observable
AIPS 1:                                    plot of main baseline for
AIPS 1:                                    display only
AIPS 1: BADDISK       1           2        Disks to avoid for scratch.
AIPS 1:               4        *rest 0
AIPS 1: BPARM      *all 0                  Subarray weight factors
AIPS 1: CPARM      *all 0                  Subarray ref. ant.

here is some of the output from the message file (this was done with no
boxes and no cleaning, just immediate selfcal):

AIPS 1 12:02:14 go
SCMAP1 12:02:14 Task SCMAP  (release of 31DEC05) begins
SCMAP1 12:02:15 Create M83_SC      .SCMAP .   2 (UV)  on disk  3  cno  428
SCMAP1 12:02:16 Create M83_C98     .ICL001.   5 (MA)  on disk  3  cno  429
SCMAP1 12:02:16 Create M83_C98     .IBM001.   5 (MA)  on disk  3  cno  430
SCMAP1 12:02:21 SCUV2S: Copied   470997 visibilities to be imaged
SCMAP1 12:02:22 Copied AN file from vol/cno/vers  3  424   1 to  3  431
1
SCMAP1 12:02:22 ******* Beginning self-calibration loop    1 *****
SCMAP1 12:02:25 UVWAIT: begin finding uniform weights
SCMAP1 12:02:27 UVWAIT: Average grid weight  1.719E+05
SCMAP1 12:02:27 UVWAIT: Adding temperance S  3.439E+04
SCMAP1 12:02:27 UVWAIT: begin applying uniform or other weights
SCMAP1 12:02:29 UVWAIT: Sum of weights in  8.411E+08 and out  1.819E+09
SCMAP1 12:02:29 UVWAIT: Noise is increased by a factor 1.198 due to
weighting
SCMAP1 12:02:29 UVWAIT: Average summed weight  1.719E+05 over     931742
vis
SCMAP1 12:02:30 GRDFLT: X and Y convolution type = SPHEROIDAL
AIPS 1: ** press RETURN for more, enter Q or next line to quit print **
#
 aragorn   AIPS (31DEC05)      2     25-APR-2005  12:04:43    Page   10
task     time   messages for user    2
SCMAP1 12:02:30 GRDFLT: X and Y parms =  3.0000  1.0000
SCMAP1 12:02:30 GRDFLT: convolution function sampled every 1/100 of a cell
SCMAP1 12:02:31 GRDMEM: Ave    2 Channels; 4.885100E+09 to 4.835100E+09 Hz
SCMAP1 12:02:34 Field   1 Sum of gridding weights =  1.93062E+10
SCMAP1 12:02:34 Field   1 Beam min =  -86.7 MilliJy, max =    1.0      Jy
SCMAP1 12:02:34 Field   1 fit FWHM =    3.544 x    3.034      arcsec, PA=
73.4
SCMAP1 12:02:34 CLBHIS: minimum component 0.500 of current peak
SCMAP1 12:02:43 Field   1 min = -965.1 MicroJy,max =   14.3 MilliJy
SCMAP1 12:02:46 Loading field   1 from -9.65E-04 to  1.43E-02 every  3
pixel
SCMAP1 12:02:50 Stop cleaning - done enough I guess!
SCMAP1 12:02:51 Checking image max/min
SCMAP1 12:02:51 Will reuse merged components brighter than    2.316
MilliJy
SCMAP1 12:02:51 RMS CLEAN residual =  231.641 MicroJy
SCMAP1 12:02:51 RMS CLEAN residual/flux =  231.641 Microratio
SCMAP1 12:02:53 Loading field  1 to TV from -9.65E-04 to  1.43E-02 every 3
pixel
SCMAP1 12:02:56 Self-cal continuing
SCMAP1 12:02:56 Determining new calibration from CLEAN model
SCMAP1 12:02:56 Copied CC file from vol/cno/vers  3  429   1 to  3  429
2
SCMAP1 12:02:56 Using       0 components with flux =   0.000      Jy
SCMAP1 12:02:56 Copied AN file from vol/cno/vers  3  424   1 to  3  428
1
SCMAP1 12:02:56 Divide data by model - first compute model by summing
SCMAP1 12:03:21 Self-cal uv range =  9.494E+04 to  9.499E+04 wavelengths
SCMAP1 12:03:21 NXTAVG: average interval set to   10.00 seconds
SCMAP1 12:03:22 ZABORS: signal 11 received
SCMAP1 12:03:22 ABORT!
AIPS 1 12:03:22 TASK QUIT WITHOUT RESUMING ME
AIPS 1 12:03:22 Begin check for any 'standard' scratch files
AIPS 1 12:03:22 Scratch files -- destroyed:    5  still active:    0
AIPS 1 12:03:22 Resumes

any ideas?

larry maddox
lam




More information about the Daip mailing list