! Started logfile: /scr/pima/rv102_uvs/J1147-0724_S_dfm.log on Sat Nov 9 09:33:51 2019 float field_size field_size = 1024 float field_cell field_cell = 0.5000 float taper_size taper_size = 20.0 float freq freq = 2.23299 obs /scr/pima/rv102_uvs/J1147-0724_S_uva.fits ! Reading UV FITS file: /scr/pima/rv102_uvs/J1147-0724_S_uva.fits ! AN table 1: 18 integrations on 1 of 1 possible baselines. ! AN table 2: 29 integrations on 36 of 36 possible baselines. ! AN table 3: 11 integrations on 21 of 21 possible baselines. ! AN table 4: 26 integrations on 36 of 36 possible baselines. ! AN table 5: 10 integrations on 3 of 3 possible baselines. ! AN table 6: 4 integrations on 6 of 6 possible baselines. ! AN table 7: 5 integrations on 10 of 10 possible baselines. ! Apparent sampling: 0.782255 visibilities/baseline/integration-bin. ! Found source: J1147-0724 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 2.23299e+09 8e+06 1 8e+06 ! 02 2 2.26299e+09 8e+06 1 8e+06 ! 03 3 2.35299e+09 8e+06 1 8e+06 ! 04 4 2.37299e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 79 lines of history. ! ! Reading 7300 visibilities. select rr, 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16 ! Restricting channel ranges to the available 4 channels. ! Selecting polarization: RR, channels: 1..4 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 print "field_size =", field_size ! field_size = 1024 print "field_cell =", field_cell ! field_cell = 0.5 print "taper_size =", taper_size ! taper_size = 20 ![@/opt64/bin/pima_mupet_01.dfm J1147-0724_S] float lim_self lim_self = 120.0 float lim_val lim_va = 1.2 integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 120 float soltime2 soltime2 = 20 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float time_av time_av = 16 float old_peak float new_peak float flux_cutoff print "===============================================================" ! =============================================================== ! print "=== Difmap script: pima_mupet_01.dvm version of 2016.05.19 === " ! === Difmap script: pima_mupet_01.dvm version of 2016.05.19 === ! mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.500x0.500 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 36 telescope corrections were flagged in sub-array 1. ! A total of 43 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 58 telescope corrections were flagged in sub-array 4. ! A total of 5 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 2. ! A total of 36 telescope corrections were flagged in sub-array 1. ! A total of 72 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 58 telescope corrections were flagged in sub-array 4. ! A total of 5 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 3. ! A total of 36 telescope corrections were flagged in sub-array 1. ! A total of 43 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 58 telescope corrections were flagged in sub-array 4. ! A total of 5 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 4. ! A total of 36 telescope corrections were flagged in sub-array 1. ! A total of 72 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 58 telescope corrections were flagged in sub-array 4. ! A total of 5 telescope corrections were flagged in sub-array 7. ! ! Fit before self-cal, rms=1.326526Jy sigma=14.666583 ! Fit after self-cal, rms=0.308686Jy sigma=2.355060 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.174 mas, bmaj=7.698 mas, bpa=-7.747 degrees ! Estimated noise=1.47819 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (0, 0). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.616368 Jy ! Component: 100 - total flux cleaned = 0.78117 Jy ! Total flux subtracted in 100 components = 0.78117 Jy ! Clean residual min=-0.038750 max=0.071788 Jy/beam ! Clean residual mean=0.000363 rms=0.007861 Jy/beam ! Combined flux in latest and established models = 0.78117 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 10 components and 0.78117 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.220729Jy sigma=1.964152 ! Fit after self-cal, rms=0.220542Jy sigma=1.963865 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0805931 Jy ! Component: 100 - total flux cleaned = 0.124361 Jy ! Total flux subtracted in 100 components = 0.124361 Jy ! Clean residual min=-0.024385 max=0.031827 Jy/beam ! Clean residual mean=0.000172 rms=0.005784 Jy/beam ! Combined flux in latest and established models = 0.90553 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 17 components and 0.90553 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.174692Jy sigma=1.357746 ! Fit after self-cal, rms=0.174481Jy sigma=1.357247 ! Inverting map print "============== Finished uniform weighting clean ===============" ! ============== Finished uniform weighting clean =============== ! uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.744 mas, bmaj=9.21 mas, bpa=-7.816 degrees ! Estimated noise=0.905578 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0333096 Jy ! Component: 100 - total flux cleaned = 0.0468468 Jy ! Component: 150 - total flux cleaned = 0.0525693 Jy ! Component: 200 - total flux cleaned = 0.0528784 Jy ! Total flux subtracted in 200 components = 0.0528784 Jy ! Clean residual min=-0.011848 max=0.020107 Jy/beam ! Clean residual mean=0.001004 rms=0.003297 Jy/beam ! Combined flux in latest and established models = 0.958409 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 28 components and 0.958409 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.167153Jy sigma=1.281912 ! Fit after self-cal, rms=0.167069Jy sigma=1.281784 ! Inverting map ! Added new window around map position (-18, 14.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0174145 Jy ! Component: 100 - total flux cleaned = 0.0252918 Jy ! Component: 150 - total flux cleaned = 0.0273894 Jy ! Component: 200 - total flux cleaned = 0.0276128 Jy ! Total flux subtracted in 200 components = 0.0276128 Jy ! Clean residual min=-0.011046 max=0.014075 Jy/beam ! Clean residual mean=0.000557 rms=0.002886 Jy/beam ! Combined flux in latest and established models = 0.986022 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 51 components and 0.986022 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.165027Jy sigma=1.254395 ! Fit after self-cal, rms=0.164899Jy sigma=1.252665 ! Inverting map print "============== Finished natural weighting clean ==============" ! ============== Finished natural weighting clean ============== ! restore ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.744 x 9.21 at -7.816 degrees (North through East) ! Clean map min=-0.0095058 max=0.82805 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! MA 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 1.00 HN 1.04 KP 0.99 LA 0.99 ! MA 1.74 PT 0.99 SC 1.05 WF 1.00 ! WZ 1.02 ! ! Telescope amplitude corrections in sub-array 3: ! FD 0.98 HN 0.99 KP 0.97 LA 1.01 ! PT 1.01 SC 0.96 WF 1.00 ! ! Telescope amplitude corrections in sub-array 4: ! FD 1.02 HN 0.97 KP 0.97 LA 1.00 ! MK 1.16 NL 1.01 PT 1.02 SC 0.99 ! WF 1.00 ! ! A total of 30 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! HN 0.00* SC 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! FD 0.98 LA 1.00 NL 0.98 PT 0.99 ! ! Telescope amplitude corrections in sub-array 7: ! FD 1.06 LA 1.04 MK 0.65 NL 0.97 ! PT 0.98 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! MA 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 1.01 HN 1.01 KP 0.99 LA 1.01 ! MA 1.00* PT 1.02 SC 0.99 WF 1.03 ! WZ 1.04 ! ! Telescope amplitude corrections in sub-array 3: ! FD 1.00 HN 0.99 KP 1.01 LA 0.99 ! PT 1.03 SC 0.98 WF 1.01 ! ! Telescope amplitude corrections in sub-array 4: ! FD 0.98 HN 0.98 KP 0.94 LA 0.97 ! MK 1.04 NL 1.03 PT 1.00 SC 0.96 ! WF 1.01 ! ! A total of 30 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! HN 0.00* SC 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! FD 0.97 LA 1.00 NL 1.01 PT 1.01 ! ! Telescope amplitude corrections in sub-array 7: ! FD 1.07 LA 0.98 MK 0.73 NL 0.99 ! PT 1.01 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! MA 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 1.02 HN 1.00 KP 0.99 LA 0.98 ! MA 0.91 PT 1.03 SC 1.01 WF 1.02 ! WZ 1.03 ! ! Telescope amplitude corrections in sub-array 3: ! FD 1.02 HN 0.99 KP 1.02 LA 1.01 ! PT 1.04 SC 1.01 WF 0.98 ! ! Telescope amplitude corrections in sub-array 4: ! FD 1.01 HN 0.98 KP 0.99 LA 1.00 ! MK 1.07 NL 1.04 PT 1.02 SC 1.00 ! WF 1.01 ! ! A total of 30 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! HN 0.00* SC 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! FD 1.04 LA 1.01 NL 1.01 PT 1.02 ! ! Telescope amplitude corrections in sub-array 7: ! FD 1.49 LA 0.99 MK 0.75 NL 1.00 ! PT 1.03 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! MA 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 1.02 HN 1.05 KP 0.98 LA 1.03 ! MA 1.00* PT 1.03 SC 1.02 WF 1.01 ! WZ 1.02 ! ! Telescope amplitude corrections in sub-array 3: ! FD 1.00 HN 1.07 KP 1.01 LA 1.01 ! PT 1.01 SC 1.00 WF 1.04 ! ! Telescope amplitude corrections in sub-array 4: ! FD 1.00 HN 1.03 KP 0.97 LA 0.99 ! MK 0.75 NL 1.00 PT 1.02 SC 1.01 ! WF 1.03 ! ! A total of 30 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! HN 0.00* SC 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! FD 1.04 LA 1.02 NL 0.99 PT 0.98 ! ! Telescope amplitude corrections in sub-array 7: ! FD 2.42 LA 1.03 MK 0.51 NL 0.94 ! PT 0.98 ! ! ! Fit before self-cal, rms=0.164899Jy sigma=1.252665 ! Fit after self-cal, rms=0.169211Jy sigma=0.995949 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.725 mas, bmaj=9.416 mas, bpa=-8.028 degrees ! Estimated noise=0.914141 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (-4.5, -1). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0109307 Jy ! Component: 100 - total flux cleaned = 0.0140155 Jy ! Component: 150 - total flux cleaned = 0.014864 Jy ! Component: 200 - total flux cleaned = 0.0151511 Jy ! Total flux subtracted in 200 components = 0.0151511 Jy ! Clean residual min=-0.007526 max=0.007554 Jy/beam ! Clean residual mean=0.000622 rms=0.001751 Jy/beam ! Combined flux in latest and established models = 1.00117 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 75 components and 1.00117 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.168178Jy sigma=0.983468 ! Fit after self-cal, rms=0.168269Jy sigma=0.982396 ! Inverting map selfcal true, true, soltime1 ! Performing amp+phase self-cal over 120 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.168269Jy sigma=0.982396 ! Fit after self-cal, rms=0.167345Jy sigma=0.980857 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.721 mas, bmaj=9.395 mas, bpa=-8.016 degrees ! Estimated noise=0.913725 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00445325 Jy ! Component: 100 - total flux cleaned = 0.00556849 Jy ! Component: 150 - total flux cleaned = 0.00639326 Jy ! Component: 200 - total flux cleaned = 0.00660574 Jy ! Total flux subtracted in 200 components = 0.00660574 Jy ! Clean residual min=-0.007522 max=0.007197 Jy/beam ! Clean residual mean=0.000408 rms=0.001666 Jy/beam ! Combined flux in latest and established models = 1.00778 Jy selfcal ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 95 components and 1.00778 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.167055Jy sigma=0.977180 ! Fit after self-cal, rms=0.167060Jy sigma=0.976827 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000716504 Jy ! Component: 100 - total flux cleaned = 0.00103209 Jy ! Component: 150 - total flux cleaned = 0.00112767 Jy ! Component: 200 - total flux cleaned = 0.00121472 Jy ! Total flux subtracted in 200 components = 0.00121472 Jy ! Clean residual min=-0.007255 max=0.007010 Jy/beam ! Clean residual mean=0.000393 rms=0.001644 Jy/beam ! Combined flux in latest and established models = 1.00899 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 104 components and 1.00899 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.166913Jy sigma=0.975189 ! Fit after self-cal, rms=0.166925Jy sigma=0.975055 ! Inverting map selfcal true, true, soltime2 ! Performing amp+phase self-cal over 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.166925Jy sigma=0.975055 ! Fit after self-cal, rms=0.163633Jy sigma=0.900724 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.718 mas, bmaj=9.387 mas, bpa=-7.996 degrees ! Estimated noise=0.913673 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00177819 Jy ! Component: 100 - total flux cleaned = 0.00223687 Jy ! Component: 150 - total flux cleaned = 0.00269959 Jy ! Component: 200 - total flux cleaned = 0.00311257 Jy ! Total flux subtracted in 200 components = 0.00311257 Jy ! Clean residual min=-0.003499 max=0.004051 Jy/beam ! Clean residual mean=0.000285 rms=0.000910 Jy/beam ! Combined flux in latest and established models = 1.01211 Jy selfcal ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 119 components and 1.01211 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.163472Jy sigma=0.899384 ! Fit after self-cal, rms=0.163503Jy sigma=0.899266 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (41.5, 70). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00264358 Jy ! Component: 100 - total flux cleaned = 0.00276574 Jy ! Component: 150 - total flux cleaned = 0.00269153 Jy ! Component: 200 - total flux cleaned = 0.00275497 Jy ! Total flux subtracted in 200 components = 0.00275497 Jy ! Clean residual min=-0.003285 max=0.004087 Jy/beam ! Clean residual mean=0.000255 rms=0.000870 Jy/beam ! Combined flux in latest and established models = 1.01486 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 137 components and 1.01486 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.163367Jy sigma=0.897381 ! Fit after self-cal, rms=0.163357Jy sigma=0.897215 ! Inverting map ! Added new window around map position (-65, 122.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0022071 Jy ! Component: 100 - total flux cleaned = 0.00286186 Jy ! Component: 150 - total flux cleaned = 0.00313883 Jy ! Component: 200 - total flux cleaned = 0.00326821 Jy ! Total flux subtracted in 200 components = 0.00326821 Jy ! Clean residual min=-0.003183 max=0.003655 Jy/beam ! Clean residual mean=0.000226 rms=0.000808 Jy/beam ! Combined flux in latest and established models = 1.01813 Jy ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 159 components and 1.01813 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.163249Jy sigma=0.895394 ! Fit after self-cal, rms=0.163250Jy sigma=0.895287 ! Inverting map ! Added new window around map position (41, -7.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00114068 Jy ! Component: 100 - total flux cleaned = 0.00103275 Jy ! Component: 150 - total flux cleaned = 0.00102723 Jy ! Component: 200 - total flux cleaned = 0.00138112 Jy ! Total flux subtracted in 200 components = 0.00138112 Jy ! Clean residual min=-0.003007 max=0.003150 Jy/beam ! Clean residual mean=0.000211 rms=0.000757 Jy/beam ! Combined flux in latest and established models = 1.01951 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 171 components and 1.01951 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.163085Jy sigma=0.893318 ! Fit after self-cal, rms=0.163072Jy sigma=0.893257 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime2 ! Performing amp+phase self-cal over 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.163072Jy sigma=0.893257 ! Fit after self-cal, rms=0.162917Jy sigma=0.892854 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.291 mas, bmaj=15.6 mas, bpa=-10.61 degrees ! Estimated noise=1.16522 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00129093 Jy ! Component: 100 - total flux cleaned = 0.00166498 Jy ! Component: 150 - total flux cleaned = 0.00180433 Jy ! Component: 200 - total flux cleaned = 0.00201845 Jy ! Total flux subtracted in 200 components = 0.00201845 Jy ! Clean residual min=-0.002766 max=0.003030 Jy/beam ! Clean residual mean=0.000400 rms=0.000734 Jy/beam ! Combined flux in latest and established models = 1.02153 Jy selfcal ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 196 components and 1.02153 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.162896Jy sigma=0.892644 ! Fit after self-cal, rms=0.162894Jy sigma=0.892595 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000341215 Jy ! Component: 100 - total flux cleaned = 0.000461848 Jy ! Component: 150 - total flux cleaned = 0.000610585 Jy ! Component: 200 - total flux cleaned = 0.000714983 Jy ! Total flux subtracted in 200 components = 0.000714983 Jy ! Clean residual min=-0.002669 max=0.002945 Jy/beam ! Clean residual mean=0.000378 rms=0.000708 Jy/beam ! Combined flux in latest and established models = 1.02224 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 202 components and 1.02224 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.162887Jy sigma=0.892553 ! Fit after self-cal, rms=0.162889Jy sigma=0.892530 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "============== Finished amplitude self-cal ==============" ! ============== Finished amplitude self-cal ============== restore ! Inverting map and beam ! Estimated beam: bmin=2.717 mas, bmaj=9.382 mas, bpa=-7.992 degrees ! Estimated noise=0.913593 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.717 x 9.382 at -7.992 degrees (North through East) ! Clean map min=-0.0027591 max=0.82857 Jy/beam if(peak(flux) > thresh) selfcal true, true clean clean_niter,clean_gain selfcal flux_cutoff = imstat(rms) * dynam repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) selfcal true, true end if ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.162889Jy sigma=0.892530 ! Fit after self-cal, rms=0.180043Jy sigma=0.806813 ! Inverting map and beam ! Estimated beam: bmin=2.716 mas, bmaj=9.383 mas, bpa=-7.989 degrees ! Estimated noise=0.914095 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000788784 Jy ! Component: 100 - total flux cleaned = 0.00107304 Jy ! Component: 150 - total flux cleaned = 0.00113967 Jy ! Component: 200 - total flux cleaned = 0.00120246 Jy ! Total flux subtracted in 200 components = 0.00120246 Jy ! Clean residual min=-0.002948 max=0.003074 Jy/beam ! Clean residual mean=0.000142 rms=0.000708 Jy/beam ! Combined flux in latest and established models = 1.02345 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 211 components and 1.02345 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.179953Jy sigma=0.805995 ! Fit after self-cal, rms=0.179961Jy sigma=0.805971 ! Inverting map ! Added new window around map position (13, 93). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00175321 Jy ! Component: 100 - total flux cleaned = 0.0021868 Jy ! Component: 150 - total flux cleaned = 0.00256637 Jy ! Component: 200 - total flux cleaned = 0.00262568 Jy ! Total flux subtracted in 200 components = 0.00262568 Jy ! Clean residual min=-0.002953 max=0.002724 Jy/beam ! Clean residual mean=0.000116 rms=0.000685 Jy/beam ! Combined flux in latest and established models = 1.02607 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 235 components and 1.02607 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.179870Jy sigma=0.804722 ! Fit after self-cal, rms=0.179883Jy sigma=0.804672 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.179883Jy sigma=0.804672 ! Fit after self-cal, rms=0.180011Jy sigma=0.804567 print "============== Clearing model and starting over ==========" ! ============== Clearing model and starting over ========== clrmod true ! clrmod: Cleared the established and tentative models. uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.167 mas, bmaj=7.714 mas, bpa=-7.841 degrees ! Estimated noise=1.46535 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.609233 Jy ! Component: 100 - total flux cleaned = 0.773811 Jy ! Component: 150 - total flux cleaned = 0.8563 Jy ! Component: 200 - total flux cleaned = 0.901947 Jy ! Total flux subtracted in 200 components = 0.901947 Jy ! Clean residual min=-0.011599 max=0.022252 Jy/beam ! Clean residual mean=0.000201 rms=0.002622 Jy/beam ! Combined flux in latest and established models = 0.901947 Jy uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Adding 16 model components to the UV plane model. ! The established model now contains 16 components and 0.901947 Jy ! Inverting map and beam ! Estimated beam: bmin=2.716 mas, bmaj=9.382 mas, bpa=-7.992 degrees ! Estimated noise=0.914405 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam keep new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0409908 Jy ! Component: 100 - total flux cleaned = 0.0627193 Jy ! Component: 150 - total flux cleaned = 0.0773214 Jy ! Component: 200 - total flux cleaned = 0.0877426 Jy ! Total flux subtracted in 200 components = 0.0877426 Jy ! Clean residual min=-0.004092 max=0.007402 Jy/beam ! Clean residual mean=0.000558 rms=0.001326 Jy/beam ! Combined flux in latest and established models = 0.98969 Jy ! Adding 38 model components to the UV plane model. ! The established model now contains 48 components and 0.98969 Jy ! Inverting map ! Added new window around map position (-10, 10.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00799427 Jy ! Component: 100 - total flux cleaned = 0.0140801 Jy ! Component: 150 - total flux cleaned = 0.0191904 Jy ! Component: 200 - total flux cleaned = 0.0234374 Jy ! Total flux subtracted in 200 components = 0.0234374 Jy ! Clean residual min=-0.003671 max=0.004916 Jy/beam ! Clean residual mean=0.000280 rms=0.000920 Jy/beam ! Combined flux in latest and established models = 1.01313 Jy ! Adding 49 model components to the UV plane model. ! The established model now contains 89 components and 1.01313 Jy ! Inverting map ! Added new window around map position (22.5, -7). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00398641 Jy ! Component: 100 - total flux cleaned = 0.00527013 Jy ! Component: 150 - total flux cleaned = 0.00655643 Jy ! Component: 200 - total flux cleaned = 0.0075009 Jy ! Total flux subtracted in 200 components = 0.0075009 Jy ! Clean residual min=-0.003489 max=0.003275 Jy/beam ! Clean residual mean=0.000186 rms=0.000789 Jy/beam ! Combined flux in latest and established models = 1.02063 Jy ! Adding 42 model components to the UV plane model. ! The established model now contains 110 components and 1.02063 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.286 mas, bmaj=15.6 mas, bpa=-10.63 degrees ! Estimated noise=1.16703 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00329622 Jy ! Component: 100 - total flux cleaned = 0.00475485 Jy ! Component: 150 - total flux cleaned = 0.00519182 Jy ! Component: 200 - total flux cleaned = 0.00558192 Jy ! Total flux subtracted in 200 components = 0.00558192 Jy ! Clean residual min=-0.002791 max=0.002551 Jy/beam ! Clean residual mean=0.000255 rms=0.000714 Jy/beam ! Combined flux in latest and established models = 1.02621 Jy flux_cutoff = imstat(rms) * dynam repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam keep new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000355648 Jy ! Component: 100 - total flux cleaned = 0.000547168 Jy ! Component: 150 - total flux cleaned = 0.000681608 Jy ! Component: 200 - total flux cleaned = 0.000766774 Jy ! Total flux subtracted in 200 components = 0.000766774 Jy ! Clean residual min=-0.002633 max=0.002479 Jy/beam ! Clean residual mean=0.000232 rms=0.000675 Jy/beam ! Combined flux in latest and established models = 1.02698 Jy ! Adding 76 model components to the UV plane model. ! The established model now contains 181 components and 1.02698 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.716 mas, bmaj=9.382 mas, bpa=-7.992 degrees ! Estimated noise=0.914405 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000932165 Jy ! Component: 100 - total flux cleaned = -0.00093324 Jy ! Component: 150 - total flux cleaned = -0.00101773 Jy ! Component: 200 - total flux cleaned = -0.000940042 Jy ! Total flux subtracted in 200 components = -0.000940042 Jy ! Clean residual min=-0.003124 max=0.003005 Jy/beam ! Clean residual mean=0.000115 rms=0.000701 Jy/beam ! Combined flux in latest and established models = 1.02604 Jy selfcal ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 205 components and 1.02604 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.180170Jy sigma=0.806079 ! Fit after self-cal, rms=0.180123Jy sigma=0.805809 wmodel J1147-0724_S_map.mod ! Writing 205 model components to file: J1147-0724_S_map.mod wobs J1147-0724_S_uvs.fits ! Writing UV FITS file: J1147-0724_S_uvs.fits wwins J1147-0724_S_map.win ! wwins: Wrote 9 windows to J1147-0724_S_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.000452615 Jy ! Component: 100 - total flux cleaned = -0.000504144 Jy ! Component: 150 - total flux cleaned = -0.000361145 Jy ! Component: 200 - total flux cleaned = -3.89893e-05 Jy ! Component: 250 - total flux cleaned = 0.000229697 Jy ! Component: 300 - total flux cleaned = 0.000576515 Jy ! Component: 350 - total flux cleaned = 0.000829976 Jy ! Component: 400 - total flux cleaned = 0.00103638 Jy ! Component: 450 - total flux cleaned = 0.00132011 Jy ! Component: 500 - total flux cleaned = 0.00163774 Jy ! Component: 550 - total flux cleaned = 0.00202761 Jy ! Component: 600 - total flux cleaned = 0.00221885 Jy ! Component: 650 - total flux cleaned = 0.00233216 Jy ! Component: 700 - total flux cleaned = 0.00244444 Jy ! Component: 750 - total flux cleaned = 0.00262825 Jy ! Component: 800 - total flux cleaned = 0.00277334 Jy ! Component: 850 - total flux cleaned = 0.00305871 Jy ! Component: 900 - total flux cleaned = 0.00302319 Jy ! Component: 950 - total flux cleaned = 0.0032326 Jy ! Component: 1000 - total flux cleaned = 0.00323238 Jy ! Component: 1050 - total flux cleaned = 0.0033684 Jy ! Component: 1100 - total flux cleaned = 0.00336817 Jy ! Component: 1150 - total flux cleaned = 0.00343505 Jy ! Component: 1200 - total flux cleaned = 0.00356698 Jy ! Component: 1250 - total flux cleaned = 0.00369771 Jy ! Component: 1300 - total flux cleaned = 0.0037298 Jy ! Component: 1350 - total flux cleaned = 0.00385832 Jy ! Component: 1400 - total flux cleaned = 0.00395419 Jy ! Component: 1450 - total flux cleaned = 0.00420622 Jy ! Component: 1500 - total flux cleaned = 0.00423692 Jy ! Component: 1550 - total flux cleaned = 0.00436092 Jy ! Component: 1600 - total flux cleaned = 0.00445344 Jy ! Component: 1650 - total flux cleaned = 0.0046969 Jy ! Component: 1700 - total flux cleaned = 0.00475713 Jy ! Component: 1750 - total flux cleaned = 0.00499635 Jy ! Component: 1800 - total flux cleaned = 0.00496606 Jy ! Component: 1850 - total flux cleaned = 0.00505443 Jy ! Component: 1900 - total flux cleaned = 0.00522989 Jy ! Component: 1950 - total flux cleaned = 0.00531693 Jy ! Component: 2000 - total flux cleaned = 0.0054323 Jy ! Component: 2050 - total flux cleaned = 0.00554651 Jy ! Component: 2100 - total flux cleaned = 0.00565998 Jy ! Component: 2150 - total flux cleaned = 0.00582899 Jy ! Component: 2200 - total flux cleaned = 0.00585692 Jy ! Component: 2250 - total flux cleaned = 0.00599576 Jy ! Component: 2300 - total flux cleaned = 0.00613369 Jy ! Component: 2350 - total flux cleaned = 0.00629809 Jy ! Component: 2400 - total flux cleaned = 0.00632532 Jy ! Component: 2450 - total flux cleaned = 0.00654161 Jy ! Component: 2500 - total flux cleaned = 0.00664901 Jy ! Component: 2550 - total flux cleaned = 0.00675565 Jy ! Component: 2600 - total flux cleaned = 0.0068884 Jy ! Component: 2650 - total flux cleaned = 0.00707262 Jy ! Component: 2700 - total flux cleaned = 0.00712491 Jy ! Component: 2750 - total flux cleaned = 0.00728116 Jy ! Component: 2800 - total flux cleaned = 0.00733289 Jy ! Component: 2850 - total flux cleaned = 0.00751294 Jy ! Component: 2900 - total flux cleaned = 0.00769176 Jy ! Component: 2950 - total flux cleaned = 0.00784421 Jy ! Component: 3000 - total flux cleaned = 0.00797039 Jy ! Component: 3050 - total flux cleaned = 0.00807069 Jy ! Component: 3100 - total flux cleaned = 0.0082203 Jy ! Component: 3150 - total flux cleaned = 0.00834442 Jy ! Component: 3200 - total flux cleaned = 0.00849255 Jy ! Component: 3250 - total flux cleaned = 0.00856618 Jy ! Component: 3300 - total flux cleaned = 0.00866379 Jy ! Component: 3350 - total flux cleaned = 0.0087852 Jy ! Component: 3400 - total flux cleaned = 0.00883337 Jy ! Component: 3450 - total flux cleaned = 0.0090254 Jy ! Component: 3500 - total flux cleaned = 0.00907302 Jy ! Component: 3550 - total flux cleaned = 0.00919184 Jy ! Component: 3600 - total flux cleaned = 0.00930968 Jy ! Component: 3650 - total flux cleaned = 0.00949741 Jy ! Component: 3700 - total flux cleaned = 0.00956745 Jy ! Component: 3750 - total flux cleaned = 0.00956732 Jy ! Component: 3800 - total flux cleaned = 0.0097983 Jy ! Component: 3850 - total flux cleaned = 0.00989015 Jy ! Component: 3900 - total flux cleaned = 0.0100042 Jy ! Component: 3950 - total flux cleaned = 0.0100271 Jy ! Component: 4000 - total flux cleaned = 0.0101856 Jy ! Component: 4050 - total flux cleaned = 0.0103206 Jy ! Total flux subtracted in 4096 components = 0.0104103 Jy ! Clean residual min=-0.001423 max=0.001365 Jy/beam ! Clean residual mean=0.000027 rms=0.000403 Jy/beam ! Combined flux in latest and established models = 1.03645 Jy keep ! Adding 262 model components to the UV plane model. ! The established model now contains 466 components and 1.03645 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /scr/pima/rv102_uvs/J1147-0724_S_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.716 x 9.382 at -7.992 degrees (North through East) ! Clean map min=-0.002561 max=0.82729 Jy/beam ! Writing clean map to FITS file: /scr/pima/rv102_uvs/J1147-0724_S_map.fits wmodel /scr/pima/rv102_uvs/J1147-0724_S_map.mod ! Writing 466 model components to file: /scr/pima/rv102_uvs/J1147-0724_S_map.mod wobs /scr/pima/rv102_uvs/J1147-0724_S_uvs.fits ! Writing UV FITS file: /scr/pima/rv102_uvs/J1147-0724_S_uvs.fits wwins /scr/pima/rv102_uvs/J1147-0724_S_map.win ! wwins: Wrote 10 windows to /scr/pima/rv102_uvs/J1147-0724_S_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.827293 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000397917 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.144296 0.144296 0.288592 0.577184 1.15437 2.30874 4.61747 9.23494 18.4699 36.9398 73.8795 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0011 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.827 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 2079.06 quit ! Quitting program ! Log file /scr/pima/rv102_uvs/J1147-0724_S_dfm.log closed on Sat Nov 9 09:33:57 2019