! Started logfile: .//J0136+4751_L_dfm.log on Wed Aug 15 15:24:57 2018 float field_size field_size = 1024 float field_cell field_cell = 0.5000 float taper_size taper_size = 20.0 float freq freq = 1.31600 obs .//J0136+4751_L_uva.fits ! Reading UV FITS file: .//J0136+4751_L_uva.fits ! AN table 1: 79 integrations on 45 of 45 possible baselines. ! Apparent sampling: 1 visibilities/baseline/integration-bin. ! Found source: J0136+4751 ! ! 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 1.316e+09 3.2e+07 1 3.2e+07 ! 02 2 1.348e+09 3.2e+07 1 3.2e+07 ! 03 3 1.38e+09 3.2e+07 1 3.2e+07 ! 04 4 1.412e+09 3.2e+07 1 3.2e+07 ! 05 5 1.444e+09 3.2e+07 1 3.2e+07 ! 06 6 1.476e+09 3.2e+07 1 3.2e+07 ! 07 7 1.508e+09 3.2e+07 1 3.2e+07 ! 08 8 1.54e+09 3.2e+07 1 3.2e+07 ! ! Polarization(s): I ! ! Read 3 lines of history. ! ! Reading 28440 visibilities. select rr, 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16 ! Polarization RR is unavailable. ! Restricting channel ranges to the available 8 channels. ! Selecting polarization: I, 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.5 print "taper_size =", taper_size ! taper_size = 20 ![@/opt64/bin/pima_mupet_01.dfm J0136+4751_L] 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 10 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 10 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 10 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 10 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 5. ! A total of 88 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 6. ! A total of 88 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 7. ! A total of 88 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 8. ! A total of 88 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=1.663081Jy sigma=40.754665 ! Fit after self-cal, rms=0.528315Jy sigma=12.060340 ! 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=4.359 mas, bmaj=7.223 mas, bpa=-10.63 degrees ! Estimated noise=0.266339 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.80779 Jy ! Component: 100 - total flux cleaned = 1.07392 Jy ! Total flux subtracted in 100 components = 1.07392 Jy ! Clean residual min=-0.046662 max=0.116871 Jy/beam ! Clean residual mean=-0.000147 rms=0.011501 Jy/beam ! Combined flux in latest and established models = 1.07392 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 14 components and 1.07392 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.483632Jy sigma=10.674526 ! Fit after self-cal, rms=0.483342Jy sigma=10.662473 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.127566 Jy ! Component: 100 - total flux cleaned = 0.200864 Jy ! Total flux subtracted in 100 components = 0.200864 Jy ! Clean residual min=-0.044900 max=0.054121 Jy/beam ! Clean residual mean=-0.000101 rms=0.010699 Jy/beam ! Combined flux in latest and established models = 1.27479 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 26 components and 1.27479 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.397909Jy sigma=9.233373 ! Fit after self-cal, rms=0.397473Jy sigma=9.222758 ! 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=5.929 mas, bmaj=9.188 mas, bpa=-10.6 degrees ! Estimated noise=0.177737 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 (-5, 6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.112132 Jy ! Component: 100 - total flux cleaned = 0.185472 Jy ! Component: 150 - total flux cleaned = 0.23817 Jy ! Component: 200 - total flux cleaned = 0.276919 Jy ! Total flux subtracted in 200 components = 0.276919 Jy ! Clean residual min=-0.032999 max=0.060007 Jy/beam ! Clean residual mean=-0.000367 rms=0.009153 Jy/beam ! Combined flux in latest and established models = 1.5517 Jy ! Performing phase self-cal ! Adding 79 model components to the UV plane model. ! The established model now contains 102 components and 1.5517 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.343594Jy sigma=8.250254 ! Fit after self-cal, rms=0.341629Jy sigma=8.202984 ! Inverting map print "============== Finished natural weighting clean ==============" ! ============== Finished natural weighting clean ============== ! restore ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 5.929 x 9.188 at -10.6 degrees (North through East) ! Clean map min=-0.03348 max=1.1255 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.08 FD 0.98 HN 1.29 KP 0.65 ! LA 0.72 MK 1.01 NL 1.04 OV 1.00 ! PT 0.88 SC 0.69 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 1.07 FD 1.02 HN 1.45 KP 1.00 ! LA 1.11 MK 1.09 NL 0.97 OV 0.99 ! PT 0.91 SC 0.71 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.07 FD 1.01 HN 1.41 KP 0.98 ! LA 1.05 MK 1.04 NL 1.03 OV 0.93 ! PT 0.91 SC 0.95 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 1.02 FD 1.04 HN 1.30 KP 0.96 ! LA 0.98 MK 1.04 NL 0.92 OV 0.92 ! PT 0.96 SC 0.80 ! ! ! Correcting IF 5. ! Telescope amplitude corrections in sub-array 1: ! BR 1.05 FD 1.06 HN 1.40 KP 1.07 ! LA 0.96 MK 0.75 NL 1.07 OV 1.01 ! PT 0.96 SC 1.00* ! ! ! Correcting IF 6. ! Telescope amplitude corrections in sub-array 1: ! BR 1.06 FD 1.04 HN 1.42 KP 1.04 ! LA 0.99 MK 1.02 NL 0.99 OV 1.02 ! PT 0.94 SC 1.00* ! ! ! Correcting IF 7. ! Telescope amplitude corrections in sub-array 1: ! BR 1.03 FD 0.97 HN 1.30 KP 1.02 ! LA 0.89 MK 0.96 NL 1.06 OV 0.95 ! PT 0.88 SC 1.00* ! ! ! Correcting IF 8. ! Telescope amplitude corrections in sub-array 1: ! BR 0.88 FD 0.90 HN 1.01 KP 0.96 ! LA 0.88 MK 0.85 NL 0.95 OV 0.93 ! PT 0.82 SC 1.00* ! ! ! Fit before self-cal, rms=0.341629Jy sigma=8.202984 ! Fit after self-cal, rms=0.244681Jy sigma=6.038563 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=5.943 mas, bmaj=9.205 mas, bpa=-7.757 degrees ! Estimated noise=0.17746 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 (9.5, 11.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00313132 Jy ! Component: 100 - total flux cleaned = -0.00531899 Jy ! Component: 150 - total flux cleaned = -0.00937184 Jy ! Component: 200 - total flux cleaned = -0.0120026 Jy ! Total flux subtracted in 200 components = -0.0120026 Jy ! Clean residual min=-0.028814 max=0.032425 Jy/beam ! Clean residual mean=0.000008 rms=0.006404 Jy/beam ! Combined flux in latest and established models = 1.5397 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 129 components and 1.5397 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.233603Jy sigma=5.772202 ! Fit after self-cal, rms=0.233226Jy sigma=5.760260 ! Inverting map ! Added new window around map position (24, 17.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00841403 Jy ! Component: 100 - total flux cleaned = 0.0103528 Jy ! Component: 150 - total flux cleaned = 0.00867134 Jy ! Component: 200 - total flux cleaned = 0.00633021 Jy ! Total flux subtracted in 200 components = 0.00633021 Jy ! Clean residual min=-0.027440 max=0.030807 Jy/beam ! Clean residual mean=0.000014 rms=0.006087 Jy/beam ! Combined flux in latest and established models = 1.54603 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 166 components and 1.54603 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.225737Jy sigma=5.569706 ! Fit after self-cal, rms=0.225400Jy sigma=5.561608 ! Inverting map ! Added new window around map position (4, 72.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00563887 Jy ! Component: 100 - total flux cleaned = 0.00558698 Jy ! Component: 150 - total flux cleaned = 0.00221021 Jy ! Component: 200 - total flux cleaned = -0.00011463 Jy ! Total flux subtracted in 200 components = -0.00011463 Jy ! Clean residual min=-0.027238 max=0.029089 Jy/beam ! Clean residual mean=0.000014 rms=0.005741 Jy/beam ! Combined flux in latest and established models = 1.54592 Jy ! Performing phase self-cal ! Adding 70 model components to the UV plane model. ! The established model now contains 226 components and 1.54592 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.218426Jy sigma=5.344187 ! Fit after self-cal, rms=0.218187Jy sigma=5.337222 ! 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.218187Jy sigma=5.337222 ! Fit after self-cal, rms=4.376530Jy sigma=5.733053 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=6.022 mas, bmaj=9.111 mas, bpa=-10.25 degrees ! Estimated noise=0.184464 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0137002 Jy ! Component: 100 - total flux cleaned = 0.0163692 Jy ! Component: 150 - total flux cleaned = 0.0214425 Jy ! Component: 200 - total flux cleaned = 0.0202872 Jy ! Total flux subtracted in 200 components = 0.0202872 Jy ! Clean residual min=-0.024181 max=0.026630 Jy/beam ! Clean residual mean=-0.000012 rms=0.005997 Jy/beam ! Combined flux in latest and established models = 1.5662 Jy selfcal ! Performing phase self-cal ! Adding 62 model components to the UV plane model. ! The established model now contains 279 components and 1.5662 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=4.361284Jy sigma=5.714559 ! Fit after self-cal, rms=4.361336Jy sigma=5.709895 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 (-10, 70). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00174456 Jy ! Component: 100 - total flux cleaned = -0.00279809 Jy ! Component: 150 - total flux cleaned = -0.00380623 Jy ! Component: 200 - total flux cleaned = -0.0057312 Jy ! Total flux subtracted in 200 components = -0.0057312 Jy ! Clean residual min=-0.023251 max=0.025637 Jy/beam ! Clean residual mean=-0.000015 rms=0.005854 Jy/beam ! Combined flux in latest and established models = 1.56047 Jy ! Performing phase self-cal ! Adding 69 model components to the UV plane model. ! The established model now contains 330 components and 1.56047 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=4.351528Jy sigma=5.733922 ! Fit after self-cal, rms=4.351481Jy sigma=5.730051 ! Inverting map ! Added new window around map position (-24, -27). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00184828 Jy ! Component: 100 - total flux cleaned = -0.00377509 Jy ! Component: 150 - total flux cleaned = -0.00559325 Jy ! Component: 200 - total flux cleaned = -0.00820305 Jy ! Total flux subtracted in 200 components = -0.00820305 Jy ! Clean residual min=-0.022957 max=0.024527 Jy/beam ! Clean residual mean=-0.000024 rms=0.005743 Jy/beam ! Combined flux in latest and established models = 1.55227 Jy ! Performing phase self-cal ! Adding 74 model components to the UV plane model. ! The established model now contains 370 components and 1.55227 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=4.343041Jy sigma=5.755343 ! Fit after self-cal, rms=4.343032Jy sigma=5.753242 ! 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=4.343032Jy sigma=5.753242 ! Fit after self-cal, rms=0.527747Jy sigma=1.838703 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=5.994 mas, bmaj=9.143 mas, bpa=-10.61 degrees ! Estimated noise=0.183714 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00230536 Jy ! Component: 100 - total flux cleaned = 0.00286893 Jy ! Component: 150 - total flux cleaned = 0.00431896 Jy ! Component: 200 - total flux cleaned = 0.00520911 Jy ! Total flux subtracted in 200 components = 0.00520911 Jy ! Clean residual min=-0.007897 max=0.006899 Jy/beam ! Clean residual mean=-0.000034 rms=0.001801 Jy/beam ! Combined flux in latest and established models = 1.55748 Jy selfcal ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 390 components and 1.55748 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.526637Jy sigma=1.685583 ! Fit after self-cal, rms=0.526546Jy sigma=1.679633 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.00193785 Jy ! Component: 100 - total flux cleaned = -0.00115636 Jy ! Component: 150 - total flux cleaned = -0.00115384 Jy ! Component: 200 - total flux cleaned = -0.00115132 Jy ! Total flux subtracted in 200 components = -0.00115132 Jy ! Clean residual min=-0.006870 max=0.006092 Jy/beam ! Clean residual mean=-0.000040 rms=0.001656 Jy/beam ! Combined flux in latest and established models = 1.55633 Jy ! Performing phase self-cal ! Adding 48 model components to the UV plane model. ! The established model now contains 414 components and 1.55633 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.526012Jy sigma=1.595177 ! Fit after self-cal, rms=0.526004Jy sigma=1.591881 ! 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. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.526004Jy sigma=1.591881 ! Fit after self-cal, rms=0.541507Jy sigma=1.564273 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=9.075 mas, bmaj=11.62 mas, bpa=-7.857 degrees ! Estimated noise=0.200772 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00560216 Jy ! Component: 100 - total flux cleaned = 0.00560131 Jy ! Component: 150 - total flux cleaned = 0.00589256 Jy ! Component: 200 - total flux cleaned = 0.00535279 Jy ! Total flux subtracted in 200 components = 0.00535279 Jy ! Clean residual min=-0.005912 max=0.005921 Jy/beam ! Clean residual mean=-0.000046 rms=0.001646 Jy/beam ! Combined flux in latest and established models = 1.56168 Jy selfcal ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 449 components and 1.56168 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.541124Jy sigma=1.494983 ! Fit after self-cal, rms=0.541137Jy sigma=1.493468 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.000539587 Jy ! Component: 100 - total flux cleaned = -0.00126359 Jy ! Component: 150 - total flux cleaned = -0.00170776 Jy ! Component: 200 - total flux cleaned = -0.00191981 Jy ! Total flux subtracted in 200 components = -0.00191981 Jy ! Clean residual min=-0.005320 max=0.005868 Jy/beam ! Clean residual mean=-0.000051 rms=0.001564 Jy/beam ! Combined flux in latest and established models = 1.55976 Jy ! Performing phase self-cal ! Adding 58 model components to the UV plane model. ! The established model now contains 470 components and 1.55976 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.540910Jy sigma=1.450992 ! Fit after self-cal, rms=0.540907Jy sigma=1.450025 ! 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=5.986 mas, bmaj=9.157 mas, bpa=-10.89 degrees ! Estimated noise=0.183556 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 5.986 x 9.157 at -10.89 degrees (North through East) ! Clean map min=-0.015102 max=1.1545 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.540907Jy sigma=1.450025 ! Fit after self-cal, rms=0.998517Jy sigma=1.407449 ! Inverting map and beam ! Estimated beam: bmin=5.977 mas, bmaj=9.166 mas, bpa=-11.05 degrees ! Estimated noise=0.183625 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000109484 Jy ! Component: 100 - total flux cleaned = -0.000512511 Jy ! Component: 150 - total flux cleaned = -0.000798196 Jy ! Component: 200 - total flux cleaned = 8.56449e-06 Jy ! Total flux subtracted in 200 components = 8.56449e-06 Jy ! Clean residual min=-0.005150 max=0.005764 Jy/beam ! Clean residual mean=-0.000040 rms=0.001370 Jy/beam ! Combined flux in latest and established models = 1.55977 Jy ! Performing phase self-cal ! Adding 44 model components to the UV plane model. ! The established model now contains 497 components and 1.55977 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.998311Jy sigma=1.341761 ! Fit after self-cal, rms=0.998315Jy sigma=1.340227 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00115197 Jy ! Component: 100 - total flux cleaned = -0.000412671 Jy ! Component: 150 - total flux cleaned = -0.000650696 Jy ! Component: 200 - total flux cleaned = -0.000647111 Jy ! Total flux subtracted in 200 components = -0.000647111 Jy ! Clean residual min=-0.004946 max=0.005521 Jy/beam ! Clean residual mean=-0.000042 rms=0.001298 Jy/beam ! Combined flux in latest and established models = 1.55912 Jy ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 508 components and 1.55912 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.998214Jy sigma=1.296359 ! Fit after self-cal, rms=0.998221Jy sigma=1.295375 ! 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.998221Jy sigma=1.295375 ! Fit after self-cal, rms=1.014018Jy sigma=1.286474 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=4.335 mas, bmaj=7.09 mas, bpa=-10.32 degrees ! Estimated noise=0.276131 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.840324 Jy ! Component: 100 - total flux cleaned = 1.12431 Jy ! Component: 150 - total flux cleaned = 1.25959 Jy ! Component: 200 - total flux cleaned = 1.33784 Jy ! Total flux subtracted in 200 components = 1.33784 Jy ! Clean residual min=-0.013709 max=0.041240 Jy/beam ! Clean residual mean=-0.000063 rms=0.003171 Jy/beam ! Combined flux in latest and established models = 1.33783 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 37 model components to the UV plane model. ! The established model now contains 37 components and 1.33783 Jy ! Inverting map and beam ! Estimated beam: bmin=5.971 mas, bmaj=9.174 mas, bpa=-11.21 degrees ! Estimated noise=0.183619 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.0827845 Jy ! Component: 100 - total flux cleaned = 0.130097 Jy ! Component: 150 - total flux cleaned = 0.161002 Jy ! Component: 200 - total flux cleaned = 0.182306 Jy ! Total flux subtracted in 200 components = 0.182306 Jy ! Clean residual min=-0.007384 max=0.012244 Jy/beam ! Clean residual mean=-0.000108 rms=0.001899 Jy/beam ! Combined flux in latest and established models = 1.52014 Jy ! Adding 91 model components to the UV plane model. ! The established model now contains 115 components and 1.52014 Jy ! Inverting map ! Added new window around map position (-0.5, 21). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0155883 Jy ! Component: 100 - total flux cleaned = 0.0233855 Jy ! Component: 150 - total flux cleaned = 0.028303 Jy ! Component: 200 - total flux cleaned = 0.0322923 Jy ! Total flux subtracted in 200 components = 0.0322923 Jy ! Clean residual min=-0.005615 max=0.006669 Jy/beam ! Clean residual mean=-0.000057 rms=0.001383 Jy/beam ! Combined flux in latest and established models = 1.55243 Jy ! Adding 53 model components to the UV plane model. ! The established model now contains 156 components and 1.55243 Jy ! Inverting map ! Added new window around map position (53.5, -69.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00400989 Jy ! Component: 100 - total flux cleaned = 0.00696595 Jy ! Component: 150 - total flux cleaned = 0.00914158 Jy ! Component: 200 - total flux cleaned = 0.0109163 Jy ! Total flux subtracted in 200 components = 0.0109163 Jy ! Clean residual min=-0.004549 max=0.006300 Jy/beam ! Clean residual mean=-0.000044 rms=0.001232 Jy/beam ! Combined flux in latest and established models = 1.56335 Jy ! Adding 47 model components to the UV plane model. ! The established model now contains 185 components and 1.56335 Jy ! Inverting map ! Added new window around map position (-8, 39.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00339678 Jy ! Component: 100 - total flux cleaned = 0.00453167 Jy ! Component: 150 - total flux cleaned = 0.00581583 Jy ! Component: 200 - total flux cleaned = 0.00642746 Jy ! Total flux subtracted in 200 components = 0.00642746 Jy ! Clean residual min=-0.004133 max=0.004326 Jy/beam ! Clean residual mean=-0.000039 rms=0.001136 Jy/beam ! Combined flux in latest and established models = 1.56978 Jy ! Adding 53 model components to the UV plane model. ! The established model now contains 214 components and 1.56978 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=9.055 mas, bmaj=11.62 mas, bpa=-8.357 degrees ! Estimated noise=0.201035 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00307787 Jy ! Component: 100 - total flux cleaned = 0.00431221 Jy ! Component: 150 - total flux cleaned = 0.00495871 Jy ! Component: 200 - total flux cleaned = 0.00555528 Jy ! Total flux subtracted in 200 components = 0.00555528 Jy ! Clean residual min=-0.004126 max=0.004471 Jy/beam ! Clean residual mean=-0.000049 rms=0.001203 Jy/beam ! Combined flux in latest and established models = 1.57533 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.00013544 Jy ! Component: 100 - total flux cleaned = 3.6948e-06 Jy ! Component: 150 - total flux cleaned = -0.000121122 Jy ! Component: 200 - total flux cleaned = -0.0001167 Jy ! Total flux subtracted in 200 components = -0.0001167 Jy ! Clean residual min=-0.003944 max=0.004278 Jy/beam ! Clean residual mean=-0.000048 rms=0.001164 Jy/beam ! Combined flux in latest and established models = 1.57522 Jy ! Adding 107 model components to the UV plane model. ! The established model now contains 316 components and 1.57522 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=5.971 mas, bmaj=9.174 mas, bpa=-11.21 degrees ! Estimated noise=0.183619 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000266678 Jy ! Component: 100 - total flux cleaned = 0.000256732 Jy ! Component: 150 - total flux cleaned = 0.000104555 Jy ! Component: 200 - total flux cleaned = -3.90611e-05 Jy ! Total flux subtracted in 200 components = -3.90611e-05 Jy ! Clean residual min=-0.003759 max=0.003951 Jy/beam ! Clean residual mean=-0.000033 rms=0.001016 Jy/beam ! Combined flux in latest and established models = 1.57518 Jy selfcal ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 334 components and 1.57518 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=1.013850Jy sigma=1.159004 ! Fit after self-cal, rms=1.013930Jy sigma=1.118989 wmodel J0136+4751_L_map.mod ! Writing 334 model components to file: J0136+4751_L_map.mod wobs J0136+4751_L_uvs.fits ! Writing UV FITS file: J0136+4751_L_uvs.fits wwins J0136+4751_L_map.win ! wwins: Wrote 10 windows to J0136+4751_L_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.000750485 Jy ! Component: 100 - total flux cleaned = 0.00174121 Jy ! Component: 150 - total flux cleaned = 0.00254495 Jy ! Component: 200 - total flux cleaned = 0.00312603 Jy ! Component: 250 - total flux cleaned = 0.00362709 Jy ! Component: 300 - total flux cleaned = 0.00405315 Jy ! Component: 350 - total flux cleaned = 0.00452969 Jy ! Component: 400 - total flux cleaned = 0.00499634 Jy ! Component: 450 - total flux cleaned = 0.00533973 Jy ! Component: 500 - total flux cleaned = 0.00578877 Jy ! Component: 550 - total flux cleaned = 0.00622993 Jy ! Component: 600 - total flux cleaned = 0.00650152 Jy ! Component: 650 - total flux cleaned = 0.00692911 Jy ! Component: 700 - total flux cleaned = 0.00719225 Jy ! Component: 750 - total flux cleaned = 0.00745174 Jy ! Component: 800 - total flux cleaned = 0.00765674 Jy ! Component: 850 - total flux cleaned = 0.00790862 Jy ! Component: 900 - total flux cleaned = 0.00800797 Jy ! Component: 950 - total flux cleaned = 0.00825387 Jy ! Component: 1000 - total flux cleaned = 0.00844782 Jy ! Component: 1050 - total flux cleaned = 0.00864005 Jy ! Component: 1100 - total flux cleaned = 0.00901886 Jy ! Component: 1150 - total flux cleaned = 0.00920668 Jy ! Component: 1200 - total flux cleaned = 0.00943942 Jy ! Component: 1250 - total flux cleaned = 0.00971553 Jy ! Component: 1300 - total flux cleaned = 0.00989805 Jy ! Component: 1350 - total flux cleaned = 0.0101242 Jy ! Component: 1400 - total flux cleaned = 0.0103937 Jy ! Component: 1450 - total flux cleaned = 0.0106157 Jy ! Component: 1500 - total flux cleaned = 0.0108364 Jy ! Component: 1550 - total flux cleaned = 0.0111421 Jy ! Component: 1600 - total flux cleaned = 0.011272 Jy ! Component: 1650 - total flux cleaned = 0.0115302 Jy ! Component: 1700 - total flux cleaned = 0.011701 Jy ! Component: 1750 - total flux cleaned = 0.0119129 Jy ! Component: 1800 - total flux cleaned = 0.0120394 Jy ! Component: 1850 - total flux cleaned = 0.0121234 Jy ! Component: 1900 - total flux cleaned = 0.0124974 Jy ! Component: 1950 - total flux cleaned = 0.0125387 Jy ! Component: 2000 - total flux cleaned = 0.0126621 Jy ! Component: 2050 - total flux cleaned = 0.0128662 Jy ! Component: 2100 - total flux cleaned = 0.0131088 Jy ! Component: 2150 - total flux cleaned = 0.0131493 Jy ! Component: 2200 - total flux cleaned = 0.0132293 Jy ! Component: 2250 - total flux cleaned = 0.0133482 Jy ! Component: 2300 - total flux cleaned = 0.0135852 Jy ! Component: 2350 - total flux cleaned = 0.0135459 Jy ! Component: 2400 - total flux cleaned = 0.0138192 Jy ! Component: 2450 - total flux cleaned = 0.0138972 Jy ! Component: 2500 - total flux cleaned = 0.01409 Jy ! Component: 2550 - total flux cleaned = 0.01432 Jy ! Component: 2600 - total flux cleaned = 0.0143202 Jy ! Component: 2650 - total flux cleaned = 0.0143962 Jy ! Component: 2700 - total flux cleaned = 0.0145847 Jy ! Component: 2750 - total flux cleaned = 0.0146598 Jy ! Component: 2800 - total flux cleaned = 0.0146971 Jy ! Component: 2850 - total flux cleaned = 0.0148826 Jy ! Component: 2900 - total flux cleaned = 0.0149562 Jy ! Component: 2950 - total flux cleaned = 0.0149192 Jy ! Component: 3000 - total flux cleaned = 0.0149925 Jy ! Component: 3050 - total flux cleaned = 0.015065 Jy ! Component: 3100 - total flux cleaned = 0.0152096 Jy ! Component: 3150 - total flux cleaned = 0.0151377 Jy ! Component: 3200 - total flux cleaned = 0.0153526 Jy ! Component: 3250 - total flux cleaned = 0.0153885 Jy ! Component: 3300 - total flux cleaned = 0.0155657 Jy ! Component: 3350 - total flux cleaned = 0.0154952 Jy ! Component: 3400 - total flux cleaned = 0.0157062 Jy ! Component: 3450 - total flux cleaned = 0.0157759 Jy ! Component: 3500 - total flux cleaned = 0.015915 Jy ! Component: 3550 - total flux cleaned = 0.0160187 Jy ! Component: 3600 - total flux cleaned = 0.0160874 Jy ! Component: 3650 - total flux cleaned = 0.0161216 Jy ! Component: 3700 - total flux cleaned = 0.0161217 Jy ! Component: 3750 - total flux cleaned = 0.0163257 Jy ! Component: 3800 - total flux cleaned = 0.0163596 Jy ! Component: 3850 - total flux cleaned = 0.0163599 Jy ! Component: 3900 - total flux cleaned = 0.0165612 Jy ! Component: 3950 - total flux cleaned = 0.0165615 Jy ! Component: 4000 - total flux cleaned = 0.0167606 Jy ! Component: 4050 - total flux cleaned = 0.0167939 Jy ! Total flux subtracted in 4096 components = 0.016794 Jy ! Clean residual min=-0.001784 max=0.001793 Jy/beam ! Clean residual mean=-0.000024 rms=0.000696 Jy/beam ! Combined flux in latest and established models = 1.59197 Jy keep ! Adding 307 model components to the UV plane model. ! The established model now contains 635 components and 1.59197 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap .//J0136+4751_L_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 5.971 x 9.174 at -11.21 degrees (North through East) ! Clean map min=-0.0065854 max=1.1643 Jy/beam ! Writing clean map to FITS file: .//J0136+4751_L_map.fits wmodel .//J0136+4751_L_map.mod ! Writing 635 model components to file: .//J0136+4751_L_map.mod wobs .//J0136+4751_L_uvs.fits ! Writing UV FITS file: .//J0136+4751_L_uvs.fits wwins .//J0136+4751_L_map.win ! wwins: Wrote 11 windows to .//J0136+4751_L_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 1.16426 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000695931 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.179324 0.179324 0.358648 0.717296 1.43459 2.86918 5.73837 11.4767 22.9535 45.9069 91.8139 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.002 print "peak= ", int(1000*peak_flux)/1000 ! peak= 1.164 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 1672.95 quit ! Quitting program ! Log file .//J0136+4751_L_dfm.log closed on Wed Aug 15 15:25:04 2018