[daip] IMAGR and UVRANGE - no residuals?
Adam Ginsburg
keflavich at gmail.com
Wed Jan 3 18:21:37 EST 2007
Dear DAIP,
I have a problem running imagr on one of my data sets that
perplexes me because it didn't happen the last time I ran
approximately the same operation. The error I'm getting is with
CLEANing. I'm trying to map a dataset with UVRANGE 1.5 1000 (because
there are errors at UV < 1.5 kl), and it works fine with uniform
weighting (I've tried robust -5 and robust 0 successfully), but when I
try robust 5 or uvwtfn 'n', it fails. Also, it worked fine when I
made a low resolution, naturally weighted map of the center of the uv
data (cellsize 1, imsize 512). It generates a dirty beam and a dirty
image just fine, as far as I can tell, but then doesn't come up with
any residuals. Any idea what I'm doing wrong?
Thanks,
Adam Ginsburg
(IMAGR output and input below)
localh> IMAGR1: Task IMAGR (release of 31DEC07) begins
localh> IMAGR1: Create C BAND N3627.IMAGR . 1 (UV) on disk 3 cno 1
localh> IMAGR1: Beginning channel 1 through 1 with 2 IFs
localh> IMAGR1: IMACPY: Copied 2266271 visibilities to be imaged
localh> IMAGR1: Create C BAND N3627.IBM001. 4 (MA) on disk 1 cno 59
localh> IMAGR1: Create C BAND N3627.ICL001. 9 (MA) on disk 1 cno 60
localh> IMAGR1: QINIT: did a GET of 5120 Kwords, OFF -344121653
localh> IMAGR1: GRDFLT: X and Y convolution type = SPHEROIDAL
localh> IMAGR1: GRDFLT: X and Y parms = 3.0000 1.0000
localh> IMAGR1: GRDFLT: convolution function sampled every 1/100 of a cell
localh> IMAGR1: GRDMEM: Ave 2 Channels; 4.885100E+09 to 4.835100E+09 Hz
localh> IMAGR1: Field 1 Sum of gridding weights = 7.75343E+11
localh> IMAGR1: Field 1 Beam min = -42.0 MilliJy, max = 1.4 Jy
localh> IMAGR1: Field 1 fit FWHM = 920.024 x 745.552 Milliarcsec, PA= -82.0
localh> IMAGR1: CLBHIS: minimum component 0.500 of current peak
localh> IMAGR1: Field 1 min = -535.9 MicroJy,max = 2.8 MilliJy
localh> IMAGR1: Loading field 1 to TV from -5.359E-04 to 2.752E-03
localh> IMAGR1: You have 600 seconds to select a menu item by:
localh> IMAGR1: Press buttons A, B, or C to choose an operation
localh> IMAGR1: Press button D for on-line help
localh> IMAGR1: Clean continuing
localh> IMAGR1: BGC Clean: using 1019 cell beam + residuals > 2.75 MilliJy
localh> IMAGR1: 0 Residual map points loaded
localh> IMAGR1: CLMPAC: NO RESIDUAL POINTS LOADED INTO AP, I QUIT
localh> IMAGR1: CLMPAC: ERROR LOADING RESIDUALS FOR CLEAN field number 0001
localh> IMAGR1: CLNCYC: ERROR CLEANING CLEAN process object
localh> IMAGR1: CLNUV1 : ERROR CLEANING CLEAN process object
localh> IMAGR1: CLNUV : ERROR CLEANING CLEAN process object
localh> IMAGR1: Deleting UV work file:
localh> IMAGR1: Destroyed 11 extension files of type AN
localh> IMAGR1: Destroyed 1 extension files of type FQ
localh> IMAGR1: Destroyed UV image file: catno= 1 disk= 3
localh> IMAGR1: Purports to die of UNNATURAL causes
localh> IMAGR1: localhost 31DEC07 TST: Cpu= 18.2 Real= 32
AIPS 1: IMAGR: Wide field imaging/Clean task
AIPS 1: Adverbs Values Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: INNAME 'C BAND N3627' Input UV data (name)
AIPS 1: INCLASS 'DBCON' 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: ANTENNAS *all 0 Antennas to plot
AIPS 1: BASELINE *all 0 Baselines with ANTENNAS
AIPS 1: DOCALIB -1 > 0 calibrate data & weights
AIPS 1: > 99 do NOT calibrate weights
AIPS 1: GAINUSE 0 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 1 Low freq. channel 0 for cont.
AIPS 1: ECHAN 0 Highest freq channel
AIPS 1: CHANNEL 0 Restart channel number
AIPS 1: NCHAV 1 Number of chan. to average.
AIPS 1: ** press RETURN for more, enter Q or next line to quit print **
#
AIPS 1: CHINC 1 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 1 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.06 0.06 (X,Y) size of grid in asec
AIPS 1: IMSIZE 512 512 Minimum image size
AIPS 1: NFIELD 1 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 21.726 21.726 RA shift per field (asec)
AIPS 1: 21.726 *rest 0
AIPS 1: DECSHIFT -46.5 -46.5 DEC shift per field (asec)
AIPS 1: -46.5 *rest 0
AIPS 1: UVTAPER 0 0 (U,V) Gaussian taper
AIPS 1: units are kilo-lambda
AIPS 1: UVRANGE 1.5 1000 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 'N' UV dist. weight function
AIPS 1: UVSIZE 0 0 Array size for doing uniform
AIPS 1: ** press RETURN for more, enter Q or next line to quit print **
#
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 1200 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 *all ' ' Input file of field params
AIPS 1: and Clean boxes; ' ' => use
AIPS 1: FLDSIZE, RASHIFT, DECSHIFT,
AIPS 1: NBOXES, CLBOX only.
AIPS 1: OBOXFILE *all ' ' Output file for final Clean
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: BMAJ 0 FWHM(asec) major axis Clean
AIPS 1: restoring beam.
AIPS 1: BMIN 0 FWHM(asec) minor axis Clean
AIPS 1: ** press RETURN for more, enter Q or next line to quit print **
#
AIPS 1: restoring beam.
AIPS 1: BPA 0 Clean beam position angle
AIPS 1: OVERLAP 0 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: (17) spectral index radius
AIPS 1: 0 -> no correction
AIPS 1: (19) Dynamic range limit
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: ** press RETURN for more, enter Q or next line to quit print **
#
AIPS 1: MAXPIXEL 0 Maximum pixels searched in
AIPS 1: each major cycle.
AIPS 1: IN3NAME ' ' Spectral index image name
AIPS 1: IN3CLASS ' ' Spectral index image class
AIPS 1: IN3SEQ 0 Spectral index image sequence
AIPS 1: number
AIPS 1: IN3DISK 0 Spectral index image disk
AIPS 1: IN4NAME ' ' Spectral curvature name
AIPS 1: IN4CLASS ' ' Spectral curvature class
AIPS 1: IN4SEQ 0 Spectral curvature sequence
AIPS 1: number
AIPS 1: IN4DISK 0 Spectral curvature disk
AIPS 1: FQTOL -1 Frequency tolerance in kHz
AIPS 1: (primary beam & spec index)
AIPS 1: DOTV 1 Display residuals on TV ?
AIPS 1: Start with field = DOTV
AIPS 1: GRCHAN 0 Graphics channel of boundary
AIPS 1: BADDISK *all 0 Disks to avoid for scratch.
More information about the Daip
mailing list