[daip] frequency dependent correction in IMAGR
Eric Greisen
egreisen at nrao.edu
Tue Oct 14 17:04:34 EDT 2014
On 10/14/2014 11:02 AM, Nicholas Wrigley wrote:
> Hi Eric,
>
> There was an error message that I remember whilst running
it via ParselTongue - imagr error -999, but I was unsure whether
it was related to this problem.
>
I hate to tell you but I am running it now with the optimized Intel
64-bit version and this time it worked. When I ran it before I made a
mistake and set IM2PARM=2 which caused it to average channels and
average in time etc etc It blew up then when it lost track of the
maximum number of channels it could do at a time. It has worked
everywhere since including the Intel optimized. I wonder if there is
something about using ParselTongue. I include all the inputs - note
that I have build up a BOXFILE which I also attach.
I hope you are running 31DEC14 and that it is current with the MNJ -
although I have not made any serious changes to IMAGR in some time.
I would suggest trying the following
NITER=1000000
flux = 1.5e-6
But I do not know why you are doing natural weights and the use of a
UVRANGE minimum > 0 is causing deep negative bowls around the sources in
the dirty maps. They clean up but doing some baseline-length dependent
time averaging would reduce the size of the data set in the imaging and
lower the effective weight of the short spacings if you set uvrange(1)=0.
AIPS 1: IMAGR: Wide field imaging/Clean task
AIPS 1: Adverbs Values Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: INNAME 'WrigleyImagr' Input UV data (name)
AIPS 1: INCLASS 'UV3220' Input UV data (class)
AIPS 1: INSEQ 1 Input UV data (seq. #)
AIPS 1: INDISK 2 Input UV data disk drive #
AIPS 1: SOURCES 'GOODSN' Source name
AIPS 1: *rest ' '
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 include
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: PDVER 0 PD table to apply (DOPOL>0)
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 1 Low freq. channel 0 for cont.
AIPS 1: ECHAN 0 Highest freq channel
AIPS 1: CHANNEL 0 Restart channel number
AIPS 1: NCHAV 8 Number of chan. to average.
AIPS 1: CHINC 8 Channel incr. between maps.
AIPS 1: BIF 0 First IF in average.
AIPS 1: EIF 0 Last IF in average.
AIPS 1: OUTNAME 'WrigleyImagr' 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.3 0.3 (X,Y) size of grid in asec
AIPS 1: IMSIZE 4096 4096 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 7 300 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: UVSIZE 0 0 Array size for doing uniform
AIPS 1: weights. 0 -> actual field
AIPS 1: size.
AIPS 1: ROBUST 7 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 4000 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 'MYWORK:Wrigley.box' Input file of field params
AIPS 1: and Clean boxes; ' ' => use
AIPS 1: FLDSIZE, RASHIFT, DECSHIFT,
AIPS 1: NBOXES, CLBOX only.
AIPS 1: OBOXFILE 'MYWORK:Wrigley.box' 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 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 0 1 => restore components to
AIPS 1: overlapped fields, >=2=>
AIPS 1: expect overlaps in Cleaning
AIPS 1: ONEBEAM -1 > 0 use only 1 dirty beam
AIPS 1: per scale in multi-facet
AIPS 1: Cleans
AIPS 1: OVRSWTCH 0 Not 0 => switch from OVERLAP
AIPS 1: >= 2 to OVERLAP 1 - see HELP
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 25 *rest 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-scale controls
AIPS 1: (17) spectral index radius
AIPS 1: 0 -> no correction
AIPS 1: (18) Limit grids (see help)
AIPS 1: (19) Dynamic range limit
AIPS 1: (20) Retry factor (see help)
AIPS 1: IM2PARM 2 *rest 0 Yet more parameters:
AIPS 1: (1) Auto boxes: allowed #
AIPS 1: (2) : island level
AIPS 1: (3) : peak required
AIPS 1: (4) : limit wrt max
AIPS 1: (5) : extend boxes
AIPS 1: (6) : edge skip
AIPS 1: (7) reset boxes for next chan
AIPS 1: (8) TV timeout interval: init
AIPS 1: (9) timeout after 1st: in sec
AIPS 1: (11) baseline-dependent avg
AIPS 1: max time in sec
AIPS 1: (12) field size 0 -> infinite
AIPS 1: (13) Number channels averaged
AIPS 1: NGAUSS 1 Number of scales to use
AIPS 1: WGAUSS *all 0 Scales in arc sec >= 0
AIPS 1: FGAUSS *all 0 Minimum flux for each resol.
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: LTYPE 3 Label TV plots (type of label
AIPS 1: with 0, 1 => do not label)
AIPS 1: BADDISK *all 0 Disks to avoid for scratch.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Wrigley.box
Type: application/vnd.previewsystems.box
Size: 2700 bytes
Desc: not available
URL: <http://listmgr.nrao.edu/pipermail/daip/attachments/20141014/a1bf2f9c/attachment-0001.box>
More information about the Daip
mailing list