[daip] IMAGR fails when re-boxing
Michael Rupen
mrupen at aoc.nrao.edu
Mon Mar 17 10:34:00 EST 2003
Dear DAIP,
running IMAGR in 31DEC03 this morning, I get the following catastrophic
failure when I REBOX interactively:
...
IMAGR1: ************* Move cursor to corner to be reset
IMAGR1: then push button A or B to do resetting
IMAGR1: or push C to go to next box *************
IMAGR1: As always, push button D to exit
IMAGR1: TVDBOX: BOX( 1) = -1.00 7.00 391.00 572.00
IMAGR1: TVDBOX: BOX( 2) = -1.00 7.00 484.00 295.00
IMAGR1: TVDBOX: BOX( 3) = -1.00 7.00 561.00 572.00
IMAGR1: TVDBOX: BOX( 4) = -1.00 7.00 932.00 702.00
IMAGR1: TVDBOX: BOX( 5) = -1.00 6.00 511.00 514.00
IMAGR1: You have 30 seconds to select a menu item by:
IMAGR1: Press buttons A, B, or C to choose an operation
IMAGR1: Press button D for on-line help
IMAGR1: Clean continuing
IMAGR1: ZTXOP2: using translated file name =
IMAGR1: ZTXOP2: /home/verdi/MPR1/tmp.oboxes
IMAGR1: OBPUT: BAD DIMENSION ARRAY
IMAGR1: KEYWORD =FIELDRES
IMAGR1: OBPUT: PROBLEM WITH OBJECT: CLEAN field number 0001
IMAGR1: CLRMAX: ERROR FINDING MAX ABS PIXEL VALUE
IMAGR1: BGC Clean: using 2001 cell beam + residuals > 49.45 MicroJy
IMAGR1: 539 Residual map points loaded
IMAGR1: Field 1 min algorithm flux = 91.912 MicroJy iter= 173
IMAGR1: Total Cleaned flux density = 5.710 MilliJy 173 comps
IMAGR1: VISDFT: Begin DFT component subtraction
IMAGR1: I Polarization model processed
IMAGR1: Field 1 min = -315.1 MicroJy,max = 405.0 MicroJy
IMAGR1: OBPUT: BAD DIMENSION ARRAY
IMAGR1: KEYWORD =FIELDRES
IMAGR1: OBPUT: PROBLEM WITH OBJECT: CLEAN field number 0001
IMAGR1: CLRMAX: ERROR FINDING MAX ABS PIXEL VALUE
IMAGR1: CLNUV : ERROR CLEANING CLEAN process object
IMAGR1: Deleting UV work file:
IMAGR1: Destroyed 1 extension files of type AN
IMAGR1: Destroyed 1 extension files of type FQ
IMAGR1: Destroyed UV image file: catno= 1308 disk= 2
IMAGR1: Purports to die of UNNATURAL causes
IMAGR1: verdi 31DEC03 TST: Cpu= 25.3 Real= 75
This has happened for several different data sets, all of which work fine
(with identical inputs) in version 'NEW'. Version 'TST' works fine if you
don't REBOX.
All runs were on verdi, a Dell Linux box.
Any thoughts?
Michael
p.s. I don't think it matters, since I happily ran with similar
inputs a few days ago, but here are the IMAGR inputs:
AIPS 1: IMAGR: Wide field imaging/Clean task
AIPS 1: Adverbs Values Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: INNAME 'X1720B 715' 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 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 0 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: NCHAV 1 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 2 2 (X,Y) size of grid in asec
AIPS 1: IMSIZE 1024 1024 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 *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 1000 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 4 Number of boxes for Clean
AIPS 1: NB: field 1 only.
AIPS 1: CLBOX -1 7 Four coordinates for each box
AIPS 1: 391 572 -1 7
AIPS 1: 484 295 -1 7
AIPS 1: 561 572 -1 7
AIPS 1: 932 702 *rest 0
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 '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 121 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.
More information about the Daip
mailing list