[daip] Re: your mail
Indra Bains
ib at star.herts.ac.uk
Tue May 22 09:19:36 EDT 2001
Hi again
Firstly, I have just found out that apparently the 31DEC99 version that I
used elsewhere (Jodrell actually) is in fact what you were referring to as
31DEC00 but when it was updated it was never renamed. So sorry for the
confusion there.
These are the inputs I used when I last tried to make the stokes I cube:
>inp
AIPS 1: IMAGR: Wide field imaging/Clean task
AIPS 1: Adverbs Values Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: INNAME '1827-145 ' Input UV data (name)
AIPS 1: INCLASS 'WTMOD ' Input UV data (class)
AIPS 1: INSEQ 2 Input UV data (seq. #)
AIPS 1: INDISK 2 Input UV data disk drive #
AIPS 1: SOURCES '1827-145 ' 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: 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 'I ' Stokes parameters (see HELP)
AIPS 1: BCHAN 21 Low freq. channel 0 for cont.
AIPS 1: ECHAN 235 Highest freq channel
AIPS 1: CHANNEL 0 Restart channel number
AIPS 1: NPOINTS 1 Number of chan. to average.
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 3 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: ** press RETURN for more, enter Q or next line to quit print **
#
AIPS 1: CELLSIZE 0.045 0.045 (X,Y) size of grid in asec
AIPS 1: IMSIZE 512 1024 Minimum image size
AIPS 1: NFIELD 1 Number of fields (max. 512)
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 'N ' 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 3000 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 1 Number of boxes for Clean
AIPS 1: NB: field 1 only.
AIPS 1: CLBOX 129 385 Four coordinates for each box
AIPS 1: 384 640 *rest 0
AIPS 1: BOXFILE *all ' ' Input file of field params
AIPS 1: and Clean boxes; ' ' => use
AIPS 1: FLDSIZE, RASHIFT, DECSHIFT,
AIPS 1: ** press RETURN for more, enter Q or next line to quit print **
#
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.03 Minimum Clean component (Jy)
AIPS 1: MINPATCH 255 Min. BEAM half-width in AP.
AIPS 1: BMAJ 0.3 FWHM(asec) major axis Clean
AIPS 1: restoring beam.
AIPS 1: BMIN 0.2 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: 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 0 0 Task enrichment parameters
AIPS 1: 0 0 (1) Antenna diameter (m)
AIPS 1: 1 *rest 0 (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-15) 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: MAXPIXEL 45000 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 1 2 Disks to avoid for scratch.
AIPS 1: *rest 0
and some of the messages - first they seemed OK:
zeus > IMAGR1: Create 1827-145 .IMAGR . 1 (UV) on disk 3 cno 7
zeus > IMAGR1: Beginning channel 21 through 21 with 1 IFs
zeus > IMAGR1: IMACPY: Copied 28946 visibilities to be imaged
zeus > IMAGR1: Create 1827-145 .IBM001. 3 (MA) on disk 3 cno 10
zeus > IMAGR1: Create 1827-145 .ICL001. 3 (MA) on disk 3 cno 11
zeus > IMAGR1: GRDFLT: X and Y convolution type = SPHEROIDAL
zeus > IMAGR1: GRDFLT: X and Y parms = 3.0000 1.0000
zeus > IMAGR1: GRDFLT: convolution function sampled every 1/100 of a cell
zeus > IMAGR1: GRIDUV: All 275 Rows in 1280k AP ( 314 Rows Max)
zeus > IMAGR1: GRIDUV: Max U Baseline 1185244. lambda (= 266 cells)
zeus > IMAGR1: GRIDUV: Frequency 1.667150E+09 Hz
zeus > IMAGR1: Field 1 Sum of gridding weights = 2.88582E+05
zeus > IMAGR1: Field 1 Beam min = -120.9 MilliJy, max = 1.0 Jy
zeus > IMAGR1: Field 1 fit FWHM = 489.441 x 124.027 Milliarcsec, PA=
16.6
zeus > IMAGR1: Field 1 use FWHM = 300.000 x 200.000 Milliarcsec, PA=
0.0
zeus > IMAGR1: Field 1 min = -43.0 MilliJy,max = 46.2 MilliJy
zeus > IMAGR1: BGC Clean: using 509 cell beam + residuals > 4.53
MilliJy
zeus > IMAGR1: 43105 Residual map points loaded
zeus > IMAGR1: Reached min. Clean flux density= 29.770 MilliJy iter=
33
zeus > IMAGR1: Total Cleaned flux density = 12.871 MilliJy 33
comps
zeus > IMAGR1: VISDFT: Begin DFT component subtraction
zeus > IMAGR1: VISDFT: Model components of type Point
zeus > IMAGR1: I Polarization model processed
zeus > IMAGR1: Field 1 min = -42.5 MilliJy,max = 44.7 MilliJy
zeus > IMAGR1: Peak residual in window 29.413 less than limit 30.000
MilliJy
zeus > IMAGR1: Beginning with field 1 scaling residuals by 0.81832
zeus > IMAGR1: Restoring Clean components
zeus > IMAGR1: Beginning channel 22 through 22 with 1 IFs
zeus > IMAGR1: GRIDUV: All 275 Rows in 1280k AP ( 314 Rows Max)
zeus > IMAGR1: GRIDUV: Max U Baseline 1185246. lambda (= 266 cells)
zeus > IMAGR1: GRIDUV: Frequency 1.667152E+09 Hz
zeus > IMAGR1: Field 1 Sum of gridding weights = 2.88582E+05
zeus > IMAGR1: Field 1 Beam min = -120.9 MilliJy, max = 1.0 Jy
zeus > IMAGR1: Field 1 fit FWHM = 489.440 x 124.027 Milliarcsec, PA=
16.6
zeus > IMAGR1: Field 1 use FWHM = 300.000 x 200.000 Milliarcsec, PA=
0.0
zeus > IMAGR1: Field 1 min = -51.4 MilliJy,max = 43.8 MilliJy
zeus > IMAGR1: BGC Clean: using 509 cell beam + residuals > 3.99
MilliJy
zeus > IMAGR1: 44948 Residual map points loaded
zeus > IMAGR1: Reached min. Clean flux density= -29.822 MilliJy iter=
32
zeus > IMAGR1: Total Cleaned flux density = 1.210 MilliJy 32
comps
zeus > LOGFILE FOR USER 40 GETTING LARGE: USE PRTMSG AND CLRMSG NOW!
zeus > IMAGR1: VISDFT: Begin DFT component subtraction
zeus > IMAGR1: I Polarization model processed
zeus > IMAGR1: Field 1 min = -50.8 MilliJy,max = 43.5 MilliJy
zeus > IMAGR1: Peak residual in window 29.913 less than limit 30.000
MilliJy
zeus > IMAGR1: Beginning with field 1 scaling residuals by 0.81833
zeus > IMAGR1: Restoring Clean components
zeus > IMAGR1: Beginning channel 23 through 23 with 1 IFs
zeus > IMAGR1: GRIDUV: All 275 Rows in 1280k AP ( 314 Rows Max)
zeus > IMAGR1: GRIDUV: Max U Baseline 1185247. lambda (= 266 cells)
zeus > IMAGR1: GRIDUV: Frequency 1.667154E+09 Hz
zeus > IMAGR1: Field 1 Sum of gridding weights = 2.88582E+05
zeus > IMAGR1: Field 1 Beam min = -120.9 MilliJy, max = 1.0 Jy
zeus > IMAGR1: Field 1 fit FWHM = 489.439 x 124.027 Milliarcsec, PA=
16.6
zeus > IMAGR1: Field 1 use FWHM = 300.000 x 200.000 Milliarcsec, PA=
0.0
zeus > IMAGR1: Field 1 min = -47.0 MilliJy,max = 41.6 MilliJy
zeus > IMAGR1: BGC Clean: using 509 cell beam + residuals > 3.97
MilliJy
zeus > IMAGR1: 44938 Residual map points loaded
zeus > IMAGR1: Reached min. Clean flux density= -29.917 MilliJy iter=
19
zeus > IMAGR1: Total Cleaned flux density = -23.464 MilliJy 19
comps
zeus > IMAGR1: VISDFT: Begin DFT component subtraction
zeus > IMAGR1: I Polarization model processed
zeus > IMAGR1: Field 1 min = -47.5 MilliJy,max = 40.9 MilliJy
zeus > IMAGR1: Peak residual in window 29.859 less than limit 30.000
MilliJy
zeus > IMAGR1: Beginning with field 1 scaling residuals by 0.81833
zeus > IMAGR1: Restoring Clean components
zeus > IMAGR1: Beginning channel 24 through 24 with 1 IFs
zeus > IMAGR1: GRIDUV: All 275 Rows in 1280k AP ( 314 Rows Max)
zeus > IMAGR1: GRIDUV: Max U Baseline 1185249. lambda (= 266 cells)
zeus > IMAGR1: GRIDUV: Frequency 1.667156E+09 Hz
zeus > IMAGR1: Field 1 Sum of gridding weights = 2.88582E+05
zeus > IMAGR1: Field 1 Beam min = -120.9 MilliJy, max = 1.0 Jy
zeus > IMAGR1: Field 1 fit FWHM = 489.439 x 124.026 Milliarcsec, PA=
16.6
zeus > IMAGR1: Field 1 use FWHM = 300.000 x 200.000 Milliarcsec, PA=
0.0
zeus > IMAGR1: Field 1 min = -45.1 MilliJy,max = 43.7 MilliJy
zeus > IMAGR1: BGC Clean: using 509 cell beam + residuals > 4.02
MilliJy
zeus > IMAGR1: 44148 Residual map points loaded
zeus > IMAGR1: Reached min. Clean flux density= -29.870 MilliJy iter=
21
zeus > IMAGR1: Total Cleaned flux density = -10.522 MilliJy 21
comps
zeus > IMAGR1: VISDFT: Begin DFT component subtraction
zeus > IMAGR1: I Polarization model processed
zeus > IMAGR1: Field 1 min = -44.7 MilliJy,max = 42.8 MilliJy
zeus > IMAGR1: Peak residual in window 29.386 less than limit 30.000
MilliJy
zeus > IMAGR1: Beginning with field 1 scaling residuals by 0.81833
zeus > IMAGR1: Restoring Clean components
zeus > IMAGR1: Beginning channel 25 through 25 with 1 IFs
zeus > IMAGR1: GRIDUV: All 275 Rows in 1280k AP ( 314 Rows Max)
zeus > IMAGR1: GRIDUV: Max U Baseline 1185250. lambda (= 266 cells)
zeus > IMAGR1: GRIDUV: Frequency 1.667158E+09 Hz
zeus > IMAGR1: Field 1 Sum of gridding weights = 2.88582E+05
zeus > IMAGR1: Field 1 Beam min = -120.9 MilliJy, max = 1.0 Jy
zeus > IMAGR1: Field 1 fit FWHM = 489.438 x 124.026 Milliarcsec, PA=
16.6
zeus > IMAGR1: Field 1 use FWHM = 300.000 x 200.000 Milliarcsec, PA=
0.0
zeus > IMAGR1: Field 1 min = -51.4 MilliJy,max = 44.7 MilliJy
zeus > IMAGR1: BGC Clean: using 509 cell beam + residuals > 3.90
MilliJy
zeus > IMAGR1: 43685 Residual map points loaded
zeus > IMAGR1: Reached min. Clean flux density= 29.842 MilliJy iter=
12
zeus > IMAGR1: Total Cleaned flux density = -484.727 MicroJy 12
comps
zeus > IMAGR1: VISDFT: Begin DFT component subtraction
zeus > IMAGR1: I Polarization model processed
zeus > IMAGR1: Field 1 min = -50.7 MilliJy,max = 43.9 MilliJy
zeus > IMAGR1: Peak residual in window 29.597 less than limit 30.000
MilliJy
zeus > IMAGR1: Beginning with field 1 scaling residuals by 0.81833
zeus > IMAGR1: Restoring Clean components
zeus > IMAGR1: Beginning channel 26 through 26 with 1 IFs
and then it tried to load too many again:
zeus > IMAGR1: Total Cleaned flux density = 27.023 MilliJy 12
comps
zeus > IMAGR1: VISDFT: Begin DFT component subtraction
zeus > IMAGR1: I Polarization model processed
zeus > IMAGR1: Field 1 min = -47.1 MilliJy,max = 47.5 MilliJy
zeus > IMAGR1: BGC Clean: using 511 cell beam + residuals > 3.64
MilliJy
zeus > IMAGR1: WARNING: > MAXPIXELS LOADED TO AP!!
zeus > IMAGR1: 45121 Residual map points loaded
zeus > IMAGR1: Reached min. Clean flux density= -29.686 MilliJy iter=
14
zeus > IMAGR1: Total Cleaned flux density = 21.051 MilliJy 14
comps
zeus > IMAGR1: VISDFT: Begin DFT component subtraction
zeus > IMAGR1: I Polarization model processed
zeus > IMAGR1: Field 1 min = -47.0 MilliJy,max = 47.6 MilliJy
zeus > IMAGR1: Peak residual in window 29.616 less than limit 30.000
MilliJy
zeus > IMAGR1: Beginning with field 1 scaling residuals by 0.81840
zeus > IMAGR1: Restoring Clean components
zeus > IMAGR1: Beginning channel 141 through 141 with 1 IFs
zeus > IMAGR1: GRIDUV: All 275 Rows in 1280k AP ( 314 Rows Max)
zeus > IMAGR1: GRIDUV: Max U Baseline 1185411. lambda (= 266 cells)
zeus > IMAGR1: GRIDUV: Frequency 1.667384E+09 Hz
zeus > IMAGR1: Field 1 Sum of gridding weights = 2.88582E+05
zeus > IMAGR1: Field 1 Beam min = -120.9 MilliJy, max = 1.0 Jy
zeus > IMAGR1: Field 1 fit FWHM = 489.370 x 124.009 Milliarcsec, PA=
16.6
zeus > IMAGR1: Field 1 use FWHM = 300.000 x 200.000 Milliarcsec, PA=
0.0
zeus > IMAGR1: Field 1 min = -50.5 MilliJy,max = 54.2 MilliJy
zeus > IMAGR1: BGC Clean: using 511 cell beam + residuals > 4.12
MilliJy
zeus > IMAGR1: 42528 Residual map points loaded
zeus > IMAGR1: Reached min. Clean flux density= -29.905 MilliJy iter=
13
zeus > IMAGR1: Total Cleaned flux density = -1.340 MilliJy 13
comps
zeus > IMAGR1: VISDFT: Begin DFT component subtraction
Please ask if you need any more information - would it help to see the
error messages I got when I used the default maxpix?
Thanks
Indra
--------------------------------------------------------------------------
I. Bains Phone (main) +44 (0)1707 284602
Astronomy Group Phone (desk) +44 (0)1707 286071
Department of Physical Sciences Fax +44 (0)1707 284644
University of Hertfordshire email ib at star.herts.ac.uk
College Lane
Hatfield
Hertfordshire
AL10 9AB
United Kingdom
--------------------------------------------------------------------------
On Mon, 21 May 2001, Eric Greisen wrote:
> Indra Bains writes:
>
> > We have the 31DEC01 TST version of AIPS installed here and I have been
> > trying to use it to CLEAN a MERLIN OH datacube with IMAGR. I left
> > maxpixels as the default of 20050 and the MSGSERVER complained:
> >
> > IMAGR1: WARNING: > MAXPIXELS LOADED TO AP!!
> >
>
> This is an indication of something really going wrong in a bad
> way. Your site is not on the midnight job but has taken a snapshot of
> 31DEC01 AIPS on some date. 31DEC99 does not exist - it was renamed
> 31DEC00 and so again you are talking about some possibly uncorrected
> snapshot. (This is why we have gone back to frozen releases).
>
> Can you send me your inputs to imagr and some of the messages coming
> out? Please send too much rather than too little andI will look at it
> - I vaguely remember seeing this message recently but do not find
> anything in the change.doc to remind what was done.
>
> Thanks,
>
> Eric Greisen
>
More information about the Daip
mailing list