[daip] IMAGR in 31DEC02: interaction goes away after SELECT FLD n
Michael Rupen
mrupen at aoc.nrao.edu
Sat May 11 16:03:12 EDT 2002
Dear Daip,
when running a multi-field IMAGR this morning in TST (31DEC02) with DOTV 1,
I found that using any of the SELECT FLD menu options gave the error
ZTTYIO: LUN = 5 NOT OPEN IN FTAB
IMAGR1: CLNTV : ERROR DISPLAYING CLEAN process object
after which IMAGR behaved as if DOTV=-1. The other menu options -- TVFIDDLE,
TVPSEUDO, SET/RESET WINDOW, TVBOX, REBOX, CONTINUE CLEAN -- all worked fine,
and this error occurred whether or not I did a "REMAKE IMAGES" first
(which also worked fine). TGET IMAGR followed by VERSION 'NEW' worked
perfectly; virtually identical inputs (with the previous day's data on the
same source) last week didn't have this problem. The maps produced seem to be
fine, so this is only an interactive bug.
I append the inputs and the results of running IMAGR in TST on verdi (one
of the new Linux boxes). The data are still on verdi so I can easily put
them somewhere for you if you'd like to play around with the same data set.
Thanks very much,
Michael
> inp imagr
AIPS 1: IMAGR: Wide field imaging/Clean task
AIPS 1: Adverbs Values Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: INNAME 'X0929-31 402' Input UV data (name)
AIPS 1: INCLASS 'CSPLIT' Input UV data (class)
AIPS 1: INSEQ 1 Input UV data (seq. #)
AIPS 1: INDISK 1 Input UV data disk drive #
AIPS 1: SOURCES *all ' ' Source name
AIPS 1: QUAL -1 Calibrator qualifier -1=>all
AIPS 1: CALCODE ' ' Calibrator code ' '=>all
AIPS 1: TIMERANG *all 0 Time range to use
AIPS 1: SELBAND -1 Bandwidth to select (kHz)
AIPS 1: SELFREQ -1 Frequency to select (MHz)
AIPS 1: FREQID -1 Freq. ID to select.
AIPS 1: SUBARRAY 0 Sub-array, 0=>all
AIPS 1: DOCALIB -1 If >0 calibrate data
AIPS 1: = 2 calibrate weights
AIPS 1: GAINUSE 4 CL (or SN) table to apply
AIPS 1: DOPOL -1 If >0.5 correct polarization.
AIPS 1: BLVER -1 BL table to apply.
AIPS 1: FLAGVER 1 Flag table version
AIPS 1: DOBAND -1 If >0.5 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: SMOOTH *all 0 Smoothing function. See
AIPS 1: HELP SMOOTH for details.
AIPS 1: STOKES ' ' Stokes parameters (see HELP)
AIPS 1: BCHAN 0 Low freq. channel 0 for cont.
AIPS 1: ECHAN 0 Highest freq channel
AIPS 1: CHANNEL 0 Restart channel number
AIPS 1: NPOINTS 7 Number of chan. to average.
AIPS 1: CHINC 7 Channel incr. between maps.
AIPS 1: BIF 0 First IF in average.
AIPS 1: EIF 0 Last IF in average.
AIPS 1: OUTNAME ' ' Output image name (name)
AIPS 1: OUTDISK 2 Output image disk drive #
AIPS 1: OUTSEQ 0 Output seq. no.
AIPS 1: OUTVER 0 CC ver. no (Continuum only)
AIPS 1: *** SET OUTVER ON RESTARTS
AIPS 1: IN2NAME ' ' UV work file name
AIPS 1: IN2CLASS ' ' UV work file class
AIPS 1: IN2SEQ 0 UV work file seq
AIPS 1: *** SET TO KEEP WORK FILE
AIPS 1: IN2DISK 0 UV work file disk
AIPS 1: CELLSIZE 0.1 0.1 (X,Y) size of grid in asec
AIPS 1: IMSIZE 512 512 Minimum image size
AIPS 1: NFIELD 4 Number of fields (max 4096)
AIPS 1: DO3DIMAG 1 > 0 => use different tangent
AIPS 1: points for each field
AIPS 1: FLDSIZE *all 0 Clean size of each field.
AIPS 1: RASHIFT *all 0 RA shift per field (asec)
AIPS 1: DECSHIFT *all 0 DEC shift per field (asec)
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: GUARD 0 0 x,y guard band fractional
AIPS 1: radius
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: SEE HELP!!
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: -5 -> pure
AIPS 1: uniform weights, 5 => 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: 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: NITER 2000 Maximum # of Clean components
AIPS 1: BCOMP *all 0 Begin at BCOMP component
AIPS 1: Specify for each field.
AIPS 1: ALLOKAY 0 For restart: > 0 => beams
AIPS 1: okay, > 1 => work file too
AIPS 1: NBOXES 0 Number of boxes for Clean
AIPS 1: NB: field 1 only.
AIPS 1: CLBOX *all 0 Four coordinates for each box
AIPS 1: BOXFILE 'MPR1:x0929-31_397c.box Input file of field params
AIPS 1: es
AIPS 1: '
AIPS 1: and Clean boxes; ' ' => use
AIPS 1: FLDSIZE, RASHIFT, DECSHIFT,
AIPS 1: NBOXES, CLBOX only.
AIPS 1: OBOXFILE 'MPR1:tmp.oboxes 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 255 Min. BEAM half-width in AP.
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: OVERLAP 2 1 => restore components to
AIPS 1: overlapped fields, >=2=>
AIPS 1: expect overlaps in Cleaning
AIPS 1: PHAT 0 Prussian hat height.
AIPS 1: FACTOR 0 Speedup factor see HELP
AIPS 1: CMETHOD ' ' Modeling method:
AIPS 1: 'DFT','GRID',' '
AIPS 1: IMAGRPRM *all 0 Task enrichment parameters
AIPS 1: (1) Antenna diameter (m)
AIPS 1: (2) Source Spectral index
AIPS 1: (3) Frequency scaling factor
AIPS 1: (4) > 0 -> SDI Clean factor
AIPS 1: (5) >0 => scale residuals
AIPS 1: (6) Half-width in x of box
AIPS 1: (7) Half-width in y of box
AIPS 1: (8) Filter components whose
AIPS 1: neighborhood is weaker than
AIPS 1: IMAGRPRM(8) Jy. 0 -> don't
AIPS 1: (9) Radius in pixels for the
AIPS 1: IMAGRPRM(8) test.
AIPS 1: (10) multiplier of image size
AIPS 1: to get beam size: 0 => 2;
AIPS 1: 2, 1, 0.5 0.25 supported
AIPS 1: (11-16) Multi-resolution
AIPS 1: added controls
AIPS 1: (20) Retry factor (see help)
AIPS 1: (20) Retry factor (see help)
AIPS 1: NGAUSS 0 Number of resolutions to use
AIPS 1: WGAUSS *all 0 Resolutions in arc sec >= 0
AIPS 1: FGAUSS *all 0 Minimum flux for each resol.
AIPS 1: MAXPIXEL 20000 Maximum pixels searched in
AIPS 1: each major cycle.
AIPS 1: DOTV 1 Display residuals on TV ?
AIPS 1: Start with field = DOTV
AIPS 1: BADDISK *all 0 Disks to avoid for scratch.
>go imagr
IMAGR1: Task IMAGR (release of 31DEC02) begins
IMAGR1: ZTXOP2: using translated file name =
IMAGR1: ZTXOP2: /home/verdi/MPR1/x0929-31_397c.boxes
IMAGR1: FIELDF: Final field parameters and # read are:
IMAGR1: Field 1 1024 1024 0.0000E+00 0.0000E+00 0
IMAGR1: Field 2 512 512 -1.0961E+01 -1.9849E+02 0
IMAGR1: Field 3 512 512 2.3550E+02 -3.6830E+01 0
IMAGR1: Field 4 512 512 -1.3989E+02 6.3950E+01 0
IMAGR1: ZTXOP2: using translated file name =
IMAGR1: ZTXOP2: /home/verdi/MPR1/x0929-31_397c.boxes
IMAGR1: WINDF: Number of clean boxes/field read as follows:
IMAGR1: Field 01: 1 Field 02: 1 Field 03: 1 Field 04: 1
IMAGR1: OBOXFILE already exists, BOXFILE not copied to it
IMAGR1: Create X0929-31 402.IMAGR . 1 (UV) on disk 2 cno 162
IMAGR1: Beginning channel 1 through 1 with 2 IFs
IMAGR1: IMACPY: Copied 217950 visibilities to be imaged
IMAGR1: UVWAIT: begin finding uniform weights
IMAGR1: UVWAIT: Average grid weight 9.669E+04
IMAGR1: UVWAIT: Adding temperance S 1.934E+04
IMAGR1: UVWAIT: begin applying uniform or other weights
IMAGR1: UVWAIT: Sum of weights in 1.700E+08 and out 3.602E+08
IMAGR1: UVWAIT: Noise is increased by a factor 1.181 due to weighting
IMAGR1: UVWAIT: Average summed weight 9.671E+04 over 434400 vis
IMAGR1: Create X0929-31 402.IBM001. 1 (MA) on disk 2 cno 163
IMAGR1: Create X0929-31 402.IBM002. 1 (MA) on disk 2 cno 164
IMAGR1: Create X0929-31 402.IBM003. 1 (MA) on disk 2 cno 165
IMAGR1: Create X0929-31 402.IBM004. 1 (MA) on disk 2 cno 166
IMAGR1: Create X0929-31 402.ICL001. 1 (MA) on disk 2 cno 167
IMAGR1: Create X0929-31 402.ICL002. 1 (MA) on disk 2 cno 168
IMAGR1: Create X0929-31 402.ICL003. 1 (MA) on disk 2 cno 169
IMAGR1: Create X0929-31 402.ICL004. 1 (MA) on disk 2 cno 170
IMAGR1: GRDFLT: X and Y convolution type = SPHEROIDAL
IMAGR1: GRDFLT: X and Y parms = 3.0000 1.0000
IMAGR1: GRDFLT: convolution function sampled every 1/100 of a cell
IMAGR1: GRDMEM: Ave 2 Channels; 4.885100E+09 to 4.835100E+09 Hz
IMAGR1: Field 1 Sum of gridding weights = 3.82380E+09
IMAGR1: Field 1 Beam min = -117.7 MilliJy, max = 1.0 Jy
IMAGR1: Field 1 fit FWHM = 1.133 x 0.332 arcsec, PA= 25.1
IMAGR1: Field 2 Sum of gridding weights = 3.82377E+09
IMAGR1: Field 2 Beam min = -117.8 MilliJy, max = 1.0 Jy
IMAGR1: Field 2 fit FWHM = 1.136 x 0.332 arcsec, PA= 25.1
IMAGR1: Field 3 Sum of gridding weights = 3.82377E+09
IMAGR1: Field 3 Beam min = -117.8 MilliJy, max = 1.0 Jy
IMAGR1: Field 3 fit FWHM = 1.132 x 0.332 arcsec, PA= 25.1
IMAGR1: Field 4 Sum of gridding weights = 3.82376E+09
IMAGR1: Field 4 Beam min = -117.7 MilliJy, max = 1.0 Jy
IMAGR1: Field 4 fit FWHM = 1.133 x 0.332 arcsec, PA= 25.2
IMAGR1: Field 1 min = -199.5 MicroJy,max = 301.3 MicroJy
IMAGR1: Field 2 min = -223.8 MicroJy,max = 1.0 MilliJy
IMAGR1: Field 3 min = -193.9 MicroJy,max = 180.5 MicroJy
IMAGR1: Field 4 min = -233.8 MicroJy,max = 341.4 MicroJy
IMAGR1: Loading field 1 from -2.00E-04 to 3.01E-04 every 2 pixel
ZTTYIO: LUN = 5 NOT OPEN IN FTAB
IMAGR1: CLNTV : ERROR DISPLAYING CLEAN process object
IMAGR1: BGC Clean: using 1021 cell beam + residuals > 42.74 MicroJy
IMAGR1: 348 Residual map points loaded
IMAGR1: Field 2 min algorithm flux = 214.558 MicroJy iter= 37
IMAGR1: Field 2 Cleaned flux density= 1.724 MilliJy 37 comps
IMAGR1: Total Cleaned flux density = 1.724 MilliJy 37 comps
IMAGR1: SETDO3: imaging done with multiple tangent planes
IMAGR1:
IMAGR1: etc.
================================================================================
More information about the Daip
mailing list