[daip] IMAGR
Eric Greisen
egreisen at nrao.edu
Thu Jun 9 15:57:53 EDT 2005
Nectaria Gizani writes:
> Hi,
>
> I was running Imagr on a very faint source, EVN data. I have boxed the
> real emission and carried on cleaning. When the new residual map came up,
> with the source "gone", I have put boxes on what I thought was also real
> emission and I turned off the tv. When imagr finished after 5000
> iterations, and I saw the map, it looked fine to me, with a low
> noise. However I had a look at the message window. It was then that I saw
> that somewhere in the 300 iteration, imagr reported "something is going
> wrong, abandon clean". However it carried on for the rest of the
> iterations ending succesfully, instead of dying!! I include the message in
> what follows. Is there something to worry about?
>
Probably not - sometimes Clean appears not to be converging so it
complains, remakes the image and procedes normally. If the message
repeats a lot then maybe something is amiss.
Eric Greisen
> Many thanks in advance,
> Nectaria
>
> ------------------
> IMAGR2 10:42:45 Task IMAGR (release of 31DEC04) begins
> IMAGR2 10:42:45 Create HERC-A .IMAGR . 1 (UV) on disk 3 cno 72
> IMAGR2 10:42:45 Beginning channel 1 through 1 with 2 IFs
> IMAGR2 10:42:45 IMACPY: Copied 9275 visibilities to be imaged
> IMAGR2 10:42:45 UVWAIT: begin finding uniform weights
> IMAGR2 10:42:45 UVWAIT: Average grid weight 1.924E-01
> IMAGR2 10:42:45 UVWAIT: Adding temperance S 3.847E-02
> IMAGR2 10:42:45 UVWAIT: begin applying uniform or other weights
> IMAGR2 10:42:45 UVWAIT: Sum of weights in 9.529E+01 and out 1.124E+02
> IMAGR2 10:42:45 UVWAIT: Noise is increased by a factor 1.516 due to
> weighting
> IMAGR2 10:42:45 UVWAIT: Average summed weight 1.924E-01 over 15444
> vis
> IMAGR2 10:42:45 Create AVG30_2 .IBM001. 2 (MA) on disk 3 cno 77
> IMAGR2 10:42:45 Create AVG30_2 .ICL001. 2 (MA) on disk 3 cno 78
> IMAGR2 10:42:46 GRDFLT: X and Y convolution type = SPHEROIDAL
> IMAGR2 10:42:46 GRDFLT: X and Y parms = 3.0000 1.0000
> IMAGR2 10:42:46 GRDFLT: convolution function sampled every 1/100 of a cell
> IMAGR2 10:42:46 GRDMEM: Ave 2 Channels; 1.654459E+09 to 1.662459E+09 Hz
> IMAGR2 10:42:46 Field 1 Sum of gridding weights = 1.19345E+03
> IMAGR2 10:42:46 Field 1 Beam min = -246.5 MilliJy, max = 1.0 Jy
> IMAGR2 10:42:46 Field 1 fit FWHM = 25.328 x 15.143 Milliarcsec, PA=
> 45.4
> IMAGR2 10:42:46 CLBHIS: minimum component 0.336 of current peak
> IMAGR2 10:42:46 Field 1 min = -1.3 MilliJy,max = 3.0 MilliJy
> IMAGR2 10:42:47 Loading field 1 to TV from -1.289E-03 to 3.015E-03
> IMAGR2 10:43:35 TVDBOX: BOX( 1) = 252.00 255.00 258.00 260.00
> IMAGR2 10:43:35 TVDBOX: BOX( 2) = 253.00 258.00 260.00 264.00
> IMAGR2 10:43:37 Clean continuing
> IMAGR2 10:43:37 BGC Clean: using 1019 cell beam + residuals > 227.91
> MicroJy
> IMAGR2 10:43:37 98 Residual map points loaded
> IMAGR2 10:43:37 Reached minimum algorithm flux = 1.115 MilliJy iter=
> 16
> IMAGR2 10:43:37 Total Cleaned flux density = 2.873 MilliJy 16
> comps
> IMAGR2 10:43:37 VISDFT: Begin DFT component subtraction
> IMAGR2 10:43:37 VISDFT: Model components of type Point
> IMAGR2 10:43:37 I Polarization model processed
> IMAGR2 10:43:37 Field 1 min = -985.2 MicroJy,max = 1.2 MilliJy
> IMAGR2 10:44:16 Clean continuing
> IMAGR2 10:44:16 BGC Clean: using 1019 cell beam + residuals > 80.23
> MicroJy
> IMAGR2 10:44:16 98 Residual map points loaded
> IMAGR2 10:44:17 Reached minimum allowed flux = 353.564 MicroJy iter=
> 53
> IMAGR2 10:44:17 Total Cleaned flux density = 5.192 MilliJy 53
> comps
> IMAGR2 10:44:17 VISDFT: Begin DFT component subtraction
> IMAGR2 10:44:17 I Polarization model processed
> IMAGR2 10:44:17 Field 1 min = -906.8 MicroJy,max = 998.0 MicroJy
> IMAGR2 10:44:17 CLNTEL: TELL operation code CHAN received
> IMAGR2 10:44:17 Loading field 1 to TV from -1.289E-03 to 3.015E-03
> IMAGR2 11:08:10 TVDBOX: BOX( 1) = 248.00 249.00 253.00 254.00
> IMAGR2 11:08:10 TVDBOX: BOX( 2) = 259.00 249.00 265.00 252.00
> IMAGR2 11:08:10 TVDBOX: BOX( 3) = 264.00 251.00 268.00 256.00
> IMAGR2 11:08:10 TVDBOX: BOX( 4) = 246.00 271.00 249.00 274.00
> IMAGR2 11:08:12 Clean continuing
> IMAGR2 11:08:12 BGC Clean: using 1019 cell beam + residuals > 73.21
> MicroJy
> IMAGR2 11:08:12 110 Residual map points loaded
> IMAGR2 11:08:12 Reached minimum allowed flux = 320.505 MicroJy iter=
> 104
> IMAGR2 11:08:12 Total Cleaned flux density = 7.892 MilliJy 104
> comps
> IMAGR2 11:08:12 VISDFT: Begin DFT component subtraction
> IMAGR2 11:08:12 I Polarization model processed
> IMAGR2 11:08:13 Field 1 min = -778.9 MicroJy,max = 773.9 MicroJy
> IMAGR2 11:08:13 Loading field 1 to TV from -1.289E-03 to 3.015E-03
> IMAGR2 11:08:27 TVDBOX: BOX( 1) = 2.00 1.00 512.00 512.00
> IMAGR2 11:08:30 Turning off DOTV option - use TELL to turn on
> IMAGR2 11:08:30 BGC Clean: using 15 cell beam + residuals > 293.17
> MicroJy
> IMAGR2 11:08:30 19984 Residual map points loaded
> IMAGR2 11:08:30 SOMETHING IS GOING WRONG - ABANDON CLEAN
> IMAGR2 11:08:30 Reached minimum algorithm flux = -577.870 MicroJy iter=
> 205
> IMAGR2 11:08:30 Total Cleaned flux density = 6.302 MilliJy 205
> comps
> IMAGR2 11:08:30 VISDFT: Begin DFT component subtraction
> IMAGR2 11:08:30 I Polarization model processed
> IMAGR2 11:08:31 Field 1 min = -595.4 MicroJy,max = 561.2 MicroJy
> IMAGR2 11:08:31 BGC Clean: using 15 cell beam + residuals > 269.43
> MicroJy
> IMAGR2 11:08:31 19992 Residual map points loaded
> IMAGR2 11:08:32 Reached minimum algorithm flux = -455.806 MicroJy iter=
> 402
> IMAGR2 11:08:32 Total Cleaned flux density = 6.506 MilliJy 402
> comps
> IMAGR2 11:08:32 ALGSTB: All 215 Rows In AP (Max 523)
> IMAGR2 11:08:32 ALGSTB: Ipol gridded model subtraction, chans 1 through
> 2
> IMAGR2 11:08:32 ALGSTB: Pass 1; 204- 0 Cells, with 9275 Pts
> IMAGR2 11:08:32 Field 1 min = -443.9 MicroJy,max = 454.6 MicroJy
> IMAGR2 11:08:32 BGC Clean: using 15 cell beam + residuals > 240.34
> MicroJy
> IMAGR2 11:08:32 19989 Residual map points loaded
> IMAGR2 11:08:33 Reached minimum algorithm flux = -377.389 MicroJy iter=
> 700
> IMAGR2 11:08:33 Total Cleaned flux density = 6.123 MilliJy 700
> comps
> IMAGR2 11:08:33 ALGSTB: All 215 Rows In AP (Max 523)
> IMAGR2 11:08:33 ALGSTB: Ipol gridded model subtraction, chans 1 through
> 2
> IMAGR2 11:08:33 ALGSTB: Pass 1; 204- 0 Cells, with 9275 Pts
> IMAGR2 11:08:33 Field 1 min = -395.6 MicroJy,max = 387.9 MicroJy
> IMAGR2 11:08:33 BGC Clean: using 15 cell beam + residuals > 210.85
> MicroJy
> IMAGR2 11:08:33 19987 Residual map points loaded
> IMAGR2 11:08:34 Reached minimum algorithm flux = -316.123 MicroJy iter=
> 1135
> IMAGR2 11:08:34 Total Cleaned flux density = 7.030 MilliJy 1135
> comps
> IMAGR2 11:08:34 ALGSTB: All 215 Rows In AP (Max 523)
> IMAGR2 11:08:34 ALGSTB: Ipol gridded model subtraction, chans 1 through
> 2
> IMAGR2 11:08:34 ALGSTB: Pass 1; 204- 0 Cells, with 9275 Pts
> IMAGR2 11:08:34 Field 1 min = -354.7 MicroJy,max = 354.5 MicroJy
> IMAGR2 11:08:35 BGC Clean: using 15 cell beam + residuals > 182.26
> MicroJy
> IMAGR2 11:08:35 19988 Residual map points loaded
> IMAGR2 11:08:35 Reached minimum algorithm flux = 263.872 MicroJy iter=
> 1745
> IMAGR2 11:08:35 Total Cleaned flux density = 7.097 MilliJy 1745
> comps
> IMAGR2 11:08:36 ALGSTB: All 215 Rows In AP (Max 523)
> IMAGR2 11:08:36 ALGSTB: Ipol gridded model subtraction, chans 1 through
> 2
> IMAGR2 11:08:36 ALGSTB: Pass 1; 204- 0 Cells, with 9275 Pts
> IMAGR2 11:08:36 Field 1 min = -325.3 MicroJy,max = 304.9 MicroJy
> IMAGR2 11:08:36 BGC Clean: using 15 cell beam + residuals > 154.13
> MicroJy
> IMAGR2 11:08:36 19992 Residual map points loaded
> IMAGR2 11:08:37 Reached minimum algorithm flux = 220.323 MicroJy iter=
> 2624
> IMAGR2 11:08:37 Total Cleaned flux density = 7.128 MilliJy 2624
> comps
> IMAGR2 11:08:37 ALGSTB: All 215 Rows In AP (Max 523)
> IMAGR2 11:08:37 ALGSTB: Ipol gridded model subtraction, chans 1 through
> 2
> IMAGR2 11:08:38 ALGSTB: Pass 1; 204- 0 Cells, with 9275 Pts
> IMAGR2 11:08:38 Field 1 min = -289.1 MicroJy,max = 336.4 MicroJy
> IMAGR2 11:08:38 CLNTEL: TELL operation code CHAN received
> IMAGR2 11:08:38 BGC Clean: using 1011 cell beam + residuals > 131.15
> MicroJy
> IMAGR2 11:08:38 19993 Residual map points loaded
> IMAGR2 11:08:44 Reached minimum algorithm flux = -172.708 MicroJy iter=
> 3501
> IMAGR2 11:08:44 Total Cleaned flux density = 7.133 MilliJy 3501
> comps
> IMAGR2 11:08:44 ALGSTB: All 215 Rows In AP (Max 523)
> IMAGR2 11:08:44 ALGSTB: Ipol gridded model subtraction, chans 1 through
> 2
> IMAGR2 11:08:44 ALGSTB: Pass 1; 204- 0 Cells, with 9275 Pts
> IMAGR2 11:08:44 Field 1 min = -223.1 MicroJy,max = 213.0 MicroJy
> IMAGR2 11:08:44 CLNTEL: TELL operation code CHAN received
> IMAGR2 11:08:44 BGC Clean: using 1011 cell beam + residuals > 111.31
> MicroJy
> IMAGR2 11:08:44 19988 Residual map points loaded
> IMAGR2 11:08:52 Reached minimum algorithm flux = -145.939 MicroJy iter=
> 4687
> IMAGR2 11:08:52 Total Cleaned flux density = 7.444 MilliJy 4687
> comps
> IMAGR2 11:08:53 ALGSTB: All 215 Rows In AP (Max 523)
> IMAGR2 11:08:53 ALGSTB: Ipol gridded model subtraction, chans 1 through
> 2
> IMAGR2 11:08:53 ALGSTB: Pass 1; 204- 0 Cells, with 9275 Pts
> IMAGR2 11:08:53 Field 1 min = -165.0 MicroJy,max = 152.9 MicroJy
> IMAGR2 11:08:53 BGC Clean: using 1011 cell beam + residuals > 97.94
> MicroJy
> IMAGR2 11:08:53 19990 Residual map points loaded
> IMAGR2 11:08:55 Reached Iter. limit, Max resid = -140.889 MicroJy iter=
> 5000
> IMAGR2 11:08:55 Total Cleaned flux density = 7.020 MilliJy 5000
> comps
> IMAGR2 11:08:55 ALGSTB: All 215 Rows In AP (Max 523)
> IMAGR2 11:08:55 ALGSTB: Ipol gridded model subtraction, chans 1 through
> 2
> IMAGR2 11:08:55 ALGSTB: Pass 1; 204- 0 Cells, with 9275 Pts
> IMAGR2 11:08:56 Field 1 min = -152.6 MicroJy,max = 157.8 MicroJy
> IMAGR2 11:08:56 Total Clean components 5000 reaches limit 5000
> IMAGR2 11:08:56 Restoring Clean components
> IMAGR2 11:08:56 Checking image max/min
> IMAGR2 11:08:56 Field 1 final Clean flux 7.020 MilliJy
> IMAGR2 11:08:56 Deleting UV work file:
> IMAGR2 11:08:56 Destroyed 1 extension files of type AN
> IMAGR2 11:08:56 Destroyed 1 extension files of type FQ
> IMAGR2 11:08:56 Destroyed UV image file: catno= 72 disk= 3
> IMAGR2 11:08:56 Appears to have ended successfully
>
>
> _______________________________________________
> Daip mailing list
> Daip at listmgr.cv.nrao.edu
> http://listmgr.cv.nrao.edu/mailman/listinfo/daip
More information about the Daip
mailing list