! Started logfile: .//J1407+2827_X_20180722_dfm.log on Wed Aug 8 14:47:34 2018 float field_size field_size = 1024 float field_cell field_cell = 0.2000 float taper_size taper_size = 70.0 float freq freq = 7.39200 obs .//J1407+2827_X_20180722_uva.fits ! Reading UV FITS file: .//J1407+2827_X_20180722_uva.fits ! AN table 1: 17 integrations on 15 of 15 possible baselines. ! Apparent sampling: 0.945098 visibilities/baseline/integration-bin. ! Found source: J1407+2827 ! ! There are 8 IFs, and a total of 8 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 7.392e+09 3.2e+07 1 3.2e+07 ! 02 2 7.424e+09 3.2e+07 1 3.2e+07 ! 03 3 7.456e+09 3.2e+07 1 3.2e+07 ! 04 4 7.552e+09 3.2e+07 1 3.2e+07 ! 05 5 7.744e+09 3.2e+07 1 3.2e+07 ! 06 6 7.776e+09 3.2e+07 1 3.2e+07 ! 07 7 7.808e+09 3.2e+07 1 3.2e+07 ! 08 8 7.84e+09 3.2e+07 1 3.2e+07 ! ! Polarization(s): RR ! ! Read 85 lines of history. ! ! Reading 1928 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 8 channels. ! Selecting polarization: RR, channels: 1..8 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 ! Reading IF 5 channels: 5..5 ! Reading IF 6 channels: 6..6 ! Reading IF 7 channels: 7..7 ! Reading IF 8 channels: 8..8 print "field_size =", field_size ! field_size = 1024 print "field_cell =", field_cell ! field_cell = 0.2 print "taper_size =", taper_size ! taper_size = 70 ![@/opt64/bin/pima_mupet_01.dfm J1407+2827_X] 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.200x0.200 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 20 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 20 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 20 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 20 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 5. ! A total of 20 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 6. ! A total of 3 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 7. ! A total of 20 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 8. ! A total of 3 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=1.253715Jy sigma=28.654416 ! Fit after self-cal, rms=0.262206Jy sigma=6.167582 ! 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.489 mas, bmaj=3.179 mas, bpa=56.94 degrees ! Estimated noise=0.988489 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=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.601411 Jy ! Component: 100 - total flux cleaned = 0.706833 Jy ! Total flux subtracted in 100 components = 0.706833 Jy ! Clean residual min=-0.070444 max=0.070633 Jy/beam ! Clean residual mean=-0.000265 rms=0.025884 Jy/beam ! Combined flux in latest and established models = 0.706833 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 18 components and 0.706833 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.150269Jy sigma=3.459151 ! Fit after self-cal, rms=0.150057Jy sigma=3.455003 ! 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.642 mas, bmaj=3.05 mas, bpa=50.11 degrees ! Estimated noise=0.826629 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.00319841 Jy ! Component: 100 - total flux cleaned = -0.00968028 Jy ! Component: 150 - total flux cleaned = 0.00168975 Jy ! Component: 200 - total flux cleaned = 0.0113451 Jy ! Total flux subtracted in 200 components = 0.0113451 Jy ! Clean residual min=-0.019397 max=0.019517 Jy/beam ! Clean residual mean=-0.000110 rms=0.008047 Jy/beam ! Combined flux in latest and established models = 0.718178 Jy ! Performing phase self-cal ! Adding 66 model components to the UV plane model. ! The established model now contains 80 components and 0.718178 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.057659Jy sigma=1.322047 ! Fit after self-cal, rms=0.057010Jy sigma=1.312382 ! 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.642 x 3.05 at 50.11 degrees (North through East) ! Clean map min=-0.04393 max=0.72202 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.96 FD 1.00* KP 0.98 NL 1.04 ! OV 0.96 PT 0.95 ! ! ! Correcting IF 2. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.96 FD 1.00* KP 0.99 NL 1.03 ! OV 0.95 PT 0.96 ! ! ! Correcting IF 3. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.98 FD 1.00* KP 0.99 NL 1.01 ! OV 0.95 PT 0.97 ! ! ! Correcting IF 4. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.01 FD 1.00* KP 1.00 NL 1.04 ! OV 0.95 PT 0.97 ! ! ! Correcting IF 5. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 FD 1.00* KP 1.02 NL 1.00 ! OV 0.97 PT 0.98 ! ! ! Correcting IF 6. ! Telescope amplitude corrections in sub-array 1: ! BR 0.97 FD 1.02 KP 1.00 NL 1.02 ! OV 0.96 PT 0.97 ! ! ! Correcting IF 7. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.97 FD 1.00* KP 1.01 NL 1.03 ! OV 0.95 PT 1.00 ! ! ! Correcting IF 8. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00 FD 1.00 KP 1.01 NL 1.03 ! OV 0.95 PT 0.99 ! ! ! Fit before self-cal, rms=0.057010Jy sigma=1.312382 ! Fit after self-cal, rms=0.043910Jy sigma=1.011162 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.644 mas, bmaj=3.085 mas, bpa=49.28 degrees ! Estimated noise=0.80721 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.0203526 Jy ! Component: 100 - total flux cleaned = 0.02309 Jy ! Component: 150 - total flux cleaned = 0.0242333 Jy ! Component: 200 - total flux cleaned = 0.024523 Jy ! Total flux subtracted in 200 components = 0.024523 Jy ! Clean residual min=-0.004389 max=0.004344 Jy/beam ! Clean residual mean=-0.000010 rms=0.001597 Jy/beam ! Combined flux in latest and established models = 0.742701 Jy ! Performing phase self-cal ! Adding 49 model components to the UV plane model. ! The established model now contains 128 components and 0.742701 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.029066Jy sigma=0.668664 ! Fit after self-cal, rms=0.028729Jy sigma=0.661772 ! 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. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.028729Jy sigma=0.661772 ! Fit after self-cal, rms=0.028163Jy sigma=0.653461 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.645 mas, bmaj=3.085 mas, bpa=49.38 degrees ! Estimated noise=0.804558 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00357706 Jy ! Component: 100 - total flux cleaned = 0.00482216 Jy ! Component: 150 - total flux cleaned = 0.00482391 Jy ! Component: 200 - total flux cleaned = 0.00442121 Jy ! Total flux subtracted in 200 components = 0.00442121 Jy ! Clean residual min=-0.002078 max=0.002077 Jy/beam ! Clean residual mean=-0.000001 rms=0.000759 Jy/beam ! Combined flux in latest and established models = 0.747123 Jy selfcal ! Performing phase self-cal ! Adding 52 model components to the UV plane model. ! The established model now contains 175 components and 0.747123 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.025918Jy sigma=0.599063 ! Fit after self-cal, rms=0.025859Jy sigma=0.597914 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.000465736 Jy ! Component: 100 - total flux cleaned = -0.000849127 Jy ! Component: 150 - total flux cleaned = -0.00110505 Jy ! Component: 200 - total flux cleaned = -0.00133526 Jy ! Total flux subtracted in 200 components = -0.00133526 Jy ! Clean residual min=-0.001235 max=0.001229 Jy/beam ! Clean residual mean=-0.000000 rms=0.000454 Jy/beam ! Combined flux in latest and established models = 0.745787 Jy ! Performing phase self-cal ! Adding 61 model components to the UV plane model. ! The established model now contains 223 components and 0.745787 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.025278Jy sigma=0.583660 ! Fit after self-cal, rms=0.025257Jy sigma=0.583247 ! 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. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.025257Jy sigma=0.583247 ! Fit after self-cal, rms=0.025095Jy sigma=0.580602 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.644 mas, bmaj=3.086 mas, bpa=49.35 degrees ! Estimated noise=0.804123 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -7.77415e-05 Jy ! Component: 100 - total flux cleaned = 0.000122855 Jy ! Component: 150 - total flux cleaned = -5.06889e-05 Jy ! Component: 200 - total flux cleaned = -0.000317103 Jy ! Total flux subtracted in 200 components = -0.000317103 Jy ! Clean residual min=-0.000842 max=0.000839 Jy/beam ! Clean residual mean=0.000000 rms=0.000302 Jy/beam ! Combined flux in latest and established models = 0.74547 Jy selfcal ! Performing phase self-cal ! Adding 53 model components to the UV plane model. ! The established model now contains 267 components and 0.74547 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.024817Jy sigma=0.573553 ! Fit after self-cal, rms=0.024808Jy sigma=0.573376 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.000591727 Jy ! Component: 100 - total flux cleaned = -0.000801132 Jy ! Component: 150 - total flux cleaned = -0.00106761 Jy ! Component: 200 - total flux cleaned = -0.00131522 Jy ! Total flux subtracted in 200 components = -0.00131522 Jy ! Clean residual min=-0.000574 max=0.000570 Jy/beam ! Clean residual mean=0.000000 rms=0.000205 Jy/beam ! Combined flux in latest and established models = 0.744155 Jy ! Performing phase self-cal ! Adding 69 model components to the UV plane model. ! The established model now contains 325 components and 0.744154 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.024698Jy sigma=0.570657 ! Fit after self-cal, rms=0.024695Jy sigma=0.570591 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 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. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.024695Jy sigma=0.570591 ! Fit after self-cal, rms=0.024660Jy sigma=0.570094 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.852 mas, bmaj=3.265 mas, bpa=52.04 degrees ! Estimated noise=0.819595 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000327003 Jy ! Component: 100 - total flux cleaned = -0.000297414 Jy ! Component: 150 - total flux cleaned = -0.000326948 Jy ! Component: 200 - total flux cleaned = -0.000461661 Jy ! Total flux subtracted in 200 components = -0.000461661 Jy ! Clean residual min=-0.000447 max=0.000448 Jy/beam ! Clean residual mean=0.000000 rms=0.000162 Jy/beam ! Combined flux in latest and established models = 0.743693 Jy selfcal ! Performing phase self-cal ! Adding 58 model components to the UV plane model. ! The established model now contains 366 components and 0.743693 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.024588Jy sigma=0.568305 ! Fit after self-cal, rms=0.024587Jy sigma=0.568262 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.000243523 Jy ! Component: 100 - total flux cleaned = -0.000314492 Jy ! Component: 150 - total flux cleaned = -0.000337067 Jy ! Component: 200 - total flux cleaned = -0.000422415 Jy ! Total flux subtracted in 200 components = -0.000422415 Jy ! Clean residual min=-0.000347 max=0.000348 Jy/beam ! Clean residual mean=0.000000 rms=0.000123 Jy/beam ! Combined flux in latest and established models = 0.74327 Jy ! Performing phase self-cal ! Adding 52 model components to the UV plane model. ! The established model now contains 401 components and 0.74327 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.024550Jy sigma=0.567373 ! Fit after self-cal, rms=0.024550Jy sigma=0.567355 ! 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.644 mas, bmaj=3.086 mas, bpa=49.34 degrees ! Estimated noise=0.803998 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.644 x 3.086 at 49.34 degrees (North through East) ! Clean map min=-0.024344 max=0.70463 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. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.024550Jy sigma=0.567355 ! Fit after self-cal, rms=0.020477Jy sigma=0.466622 ! Inverting map and beam ! Estimated beam: bmin=2.644 mas, bmaj=3.086 mas, bpa=49.37 degrees ! Estimated noise=0.803845 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 6.94771e-05 Jy ! Component: 100 - total flux cleaned = 0.000155501 Jy ! Component: 150 - total flux cleaned = 0.000277844 Jy ! Component: 200 - total flux cleaned = 0.000374453 Jy ! Total flux subtracted in 200 components = 0.000374453 Jy ! Clean residual min=-0.000318 max=0.000315 Jy/beam ! Clean residual mean=0.000001 rms=0.000101 Jy/beam ! Combined flux in latest and established models = 0.743645 Jy ! Performing phase self-cal ! Adding 39 model components to the UV plane model. ! The established model now contains 431 components and 0.743645 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.020437Jy sigma=0.465687 ! Fit after self-cal, rms=0.020437Jy sigma=0.465667 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000109513 Jy ! Component: 100 - total flux cleaned = 0.000194303 Jy ! Component: 150 - total flux cleaned = 0.000194553 Jy ! Component: 200 - total flux cleaned = 0.000194135 Jy ! Total flux subtracted in 200 components = 0.000194135 Jy ! Clean residual min=-0.000250 max=0.000253 Jy/beam ! Clean residual mean=0.000000 rms=0.000083 Jy/beam ! Combined flux in latest and established models = 0.743839 Jy ! Performing phase self-cal ! Adding 50 model components to the UV plane model. ! The established model now contains 468 components and 0.743839 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.020415Jy sigma=0.465091 ! Fit after self-cal, rms=0.020415Jy sigma=0.465078 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.020415Jy sigma=0.465078 ! Fit after self-cal, rms=0.020409Jy sigma=0.465003 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.492 mas, bmaj=3.2 mas, bpa=55.51 degrees ! Estimated noise=0.964076 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.587028 Jy ! Component: 100 - total flux cleaned = 0.688622 Jy ! Component: 150 - total flux cleaned = 0.675925 Jy ! Component: 200 - total flux cleaned = 0.655953 Jy ! Total flux subtracted in 200 components = 0.655953 Jy ! Clean residual min=-0.031850 max=0.031905 Jy/beam ! Clean residual mean=-0.000132 rms=0.012531 Jy/beam ! Combined flux in latest and established models = 0.655952 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 40 model components to the UV plane model. ! The established model now contains 40 components and 0.655952 Jy ! Inverting map and beam ! Estimated beam: bmin=2.644 mas, bmaj=3.086 mas, bpa=49.38 degrees ! Estimated noise=0.803754 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.0145819 Jy ! Component: 100 - total flux cleaned = -0.0209765 Jy ! Component: 150 - total flux cleaned = -0.0199538 Jy ! Component: 200 - total flux cleaned = -0.0199857 Jy ! Total flux subtracted in 200 components = -0.0199857 Jy ! Clean residual min=-0.010126 max=0.010186 Jy/beam ! Clean residual mean=-0.000051 rms=0.004169 Jy/beam ! Combined flux in latest and established models = 0.635967 Jy ! Adding 67 model components to the UV plane model. ! The established model now contains 102 components and 0.635967 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.852 mas, bmaj=3.265 mas, bpa=52.09 degrees ! Estimated noise=0.819359 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00118528 Jy ! Component: 100 - total flux cleaned = 0.00656627 Jy ! Component: 150 - total flux cleaned = 0.00967658 Jy ! Component: 200 - total flux cleaned = 0.00997858 Jy ! Total flux subtracted in 200 components = 0.00997858 Jy ! Clean residual min=-0.004069 max=0.004047 Jy/beam ! Clean residual mean=-0.000020 rms=0.001668 Jy/beam ! Combined flux in latest and established models = 0.645945 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.000705061 Jy ! Component: 100 - total flux cleaned = -0.00090921 Jy ! Component: 150 - total flux cleaned = -0.00110012 Jy ! Component: 200 - total flux cleaned = -0.00160208 Jy ! Total flux subtracted in 200 components = -0.00160208 Jy ! Clean residual min=-0.002642 max=0.002675 Jy/beam ! Clean residual mean=-0.000011 rms=0.000961 Jy/beam ! Combined flux in latest and established models = 0.644343 Jy ! Adding 120 model components to the UV plane model. ! The established model now contains 209 components and 0.644343 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.644 mas, bmaj=3.086 mas, bpa=49.38 degrees ! Estimated noise=0.803754 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000292183 Jy ! Component: 100 - total flux cleaned = -0.000419942 Jy ! Component: 150 - total flux cleaned = -0.0001331 Jy ! Component: 200 - total flux cleaned = -5.81179e-05 Jy ! Total flux subtracted in 200 components = -5.81179e-05 Jy ! Clean residual min=-0.001265 max=0.001259 Jy/beam ! Clean residual mean=-0.000001 rms=0.000456 Jy/beam ! Combined flux in latest and established models = 0.644285 Jy selfcal ! Performing phase self-cal ! Adding 53 model components to the UV plane model. ! The established model now contains 261 components and 0.644285 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.021279Jy sigma=0.486514 ! Fit after self-cal, rms=0.021135Jy sigma=0.482651 wmodel J1407+2827_X_map.mod ! Writing 261 model components to file: J1407+2827_X_map.mod wobs J1407+2827_X_uvs.fits ! Writing UV FITS file: J1407+2827_X_uvs.fits wwins J1407+2827_X_map.win ! wwins: No CLEAN windows to write. ! Error occured in command: wwins ![Exited script file: /opt64/bin/pima_mupet_01.dfm] ! wmap .//J1407+2827_X_20180722_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.644 x 3.086 at 49.38 degrees (North through East) ! Clean map min=-0.025343 max=0.6888 Jy/beam ! Writing clean map to FITS file: .//J1407+2827_X_20180722_map.fits wmodel .//J1407+2827_X_20180722_map.mod ! Writing 261 model components to file: .//J1407+2827_X_20180722_map.mod wobs .//J1407+2827_X_20180722_uvs.fits ! Writing UV FITS file: .//J1407+2827_X_20180722_uvs.fits wwins .//J1407+2827_X_20180722_map.win ! wwins: No CLEAN windows to write. ! Error occured in command: wwins ! float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.688796 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000426685 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.18584 0.18584 0.37168 0.743359 1.48672 2.97344 5.94687 11.8937 23.7875 47.575 95.15 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0012 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.688 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 1614.29 quit ! Quitting program ! Log file .//J1407+2827_X_20180722_dfm.log closed on Wed Aug 8 14:47:37 2018