! Started logfile: /image_prod/ug002_uva//J0735-1735_S_dfm.log on Wed Mar 4 21:53:51 2020 float field_size field_size = 1024 float field_cell field_cell = 0.5000 float taper_size taper_size = 20.0 float freq freq = 2.22000 obs /image_prod/ug002_uva//J0735-1735_S_uva.fits ! Reading UV FITS file: /image_prod/ug002_uva//J0735-1735_S_uva.fits ! AN table 1: 37 integrations on 36 of 36 possible baselines. ! AN table 2: 37 integrations on 36 of 36 possible baselines. ! Apparent sampling: 0.988739 visibilities/baseline/integration-bin. ! Found source: J0735-1735 ! ! There are 3 IFs, and a total of 3 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 2.22e+09 3.2e+07 1 3.2e+07 ! 02 2 2.252e+09 3.2e+07 1 3.2e+07 ! 03 3 2.284e+09 3.2e+07 1 3.2e+07 ! ! Polarization(s): RR ! ! Read 44 lines of history. ! ! Reading 7902 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 3 channels. ! Selecting polarization: RR, channels: 1..3 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 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 J0735-1735_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 1 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 2. ! A total of 1 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 3. ! A total of 1 telescope corrections were flagged in sub-array 2. ! ! Fit before self-cal, rms=1.444418Jy sigma=6.556775 ! Fit after self-cal, rms=0.753620Jy sigma=3.735161 ! 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.36 mas, bmaj=6.548 mas, bpa=-4.927 degrees ! Estimated noise=2.81497 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.465071 Jy ! Component: 100 - total flux cleaned = 0.714062 Jy ! Total flux subtracted in 100 components = 0.714062 Jy ! Clean residual min=-0.074490 max=0.168983 Jy/beam ! Clean residual mean=-0.000552 rms=0.025912 Jy/beam ! Combined flux in latest and established models = 0.714062 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 32 components and 0.714062 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.671797Jy sigma=2.320919 ! Fit after self-cal, rms=0.670473Jy sigma=2.315578 ! Inverting map ! Added new window around map position (-5.5, -9). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.193648 Jy ! Component: 100 - total flux cleaned = 0.345066 Jy ! Total flux subtracted in 100 components = 0.345066 Jy ! Clean residual min=-0.060000 max=0.130020 Jy/beam ! Clean residual mean=-0.000415 rms=0.019547 Jy/beam ! Combined flux in latest and established models = 1.05913 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 59 components and 1.05913 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.517261Jy sigma=1.765513 ! Fit after self-cal, rms=0.512805Jy sigma=1.750956 ! Inverting map ! Added new window around map position (3.5, 1.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.139518 Jy ! Component: 100 - total flux cleaned = 0.247183 Jy ! Total flux subtracted in 100 components = 0.247183 Jy ! Clean residual min=-0.047059 max=0.069724 Jy/beam ! Clean residual mean=-0.000314 rms=0.015307 Jy/beam ! Combined flux in latest and established models = 1.30631 Jy ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 84 components and 1.30631 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.409606Jy sigma=1.384490 ! Fit after self-cal, rms=0.407486Jy sigma=1.376726 ! 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.723 mas, bmaj=7.035 mas, bpa=-4.487 degrees ! Estimated noise=1.6417 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.140593 Jy ! Component: 100 - total flux cleaned = 0.24828 Jy ! Component: 150 - total flux cleaned = 0.3331 Jy ! Component: 200 - total flux cleaned = 0.400837 Jy ! Total flux subtracted in 200 components = 0.400837 Jy ! Clean residual min=-0.044236 max=0.066660 Jy/beam ! Clean residual mean=-0.000389 rms=0.013833 Jy/beam ! Combined flux in latest and established models = 1.70715 Jy ! Performing phase self-cal ! Adding 75 model components to the UV plane model. ! The established model now contains 139 components and 1.70715 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.276110Jy sigma=0.973061 ! Fit after self-cal, rms=0.273091Jy sigma=0.960133 ! 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.723 x 7.035 at -4.487 degrees (North through East) ! Clean map min=-0.04522 max=0.61699 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.71 FD 0.88 HN 1.07 KP 0.90 ! LA 1.11 MK 1.48 NL 0.99 OV 0.83 ! SC 0.96 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.08 FD 0.85 HN 1.21 KP 0.94 ! LA 0.95 MK 1.27 OV 0.95 PT 0.97 ! SC 1.09 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.94 FD 0.93 HN 1.07 KP 0.88 ! LA 1.25 MK 1.14 NL 1.03 OV 0.87 ! SC 1.15 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.19 FD 0.94 HN 1.19 KP 0.93 ! LA 0.99 MK 1.32 OV 0.95 PT 0.93 ! SC 1.28 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.78 FD 0.89 HN 1.13 KP 0.90 ! LA 1.07 MK 1.03 NL 1.04 OV 0.84 ! SC 0.92 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.22 FD 0.87 HN 1.23 KP 0.93 ! LA 0.89 MK 1.47 OV 0.96 PT 0.94 ! SC 1.08 ! ! ! Fit before self-cal, rms=0.273091Jy sigma=0.960133 ! Fit after self-cal, rms=0.188653Jy sigma=0.710427 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.884 mas, bmaj=7.282 mas, bpa=-3.154 degrees ! Estimated noise=1.71378 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 (12, 1.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0645165 Jy ! Component: 100 - total flux cleaned = 0.107958 Jy ! Component: 150 - total flux cleaned = 0.140445 Jy ! Component: 200 - total flux cleaned = 0.16574 Jy ! Total flux subtracted in 200 components = 0.16574 Jy ! Clean residual min=-0.021773 max=0.033094 Jy/beam ! Clean residual mean=-0.000045 rms=0.006745 Jy/beam ! Combined flux in latest and established models = 1.87289 Jy ! Performing phase self-cal ! Adding 49 model components to the UV plane model. ! The established model now contains 187 components and 1.87289 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.146458Jy sigma=0.564932 ! Fit after self-cal, rms=0.144387Jy sigma=0.554343 ! 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. ! ! Fit before self-cal, rms=0.144387Jy sigma=0.554343 ! Fit after self-cal, rms=0.132240Jy sigma=0.537660 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.891 mas, bmaj=7.297 mas, bpa=-2.746 degrees ! Estimated noise=1.71891 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0307604 Jy ! Component: 100 - total flux cleaned = 0.0527447 Jy ! Component: 150 - total flux cleaned = 0.0704514 Jy ! Component: 200 - total flux cleaned = 0.0811735 Jy ! Total flux subtracted in 200 components = 0.0811735 Jy ! Clean residual min=-0.017995 max=0.027715 Jy/beam ! Clean residual mean=-0.000021 rms=0.005096 Jy/beam ! Combined flux in latest and established models = 1.95406 Jy selfcal ! Performing phase self-cal ! Adding 54 model components to the UV plane model. ! The established model now contains 229 components and 1.95406 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.118087Jy sigma=0.492419 ! Fit after self-cal, rms=0.117419Jy sigma=0.489136 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 (21, -87). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.023829 Jy ! Component: 100 - total flux cleaned = 0.0341611 Jy ! Component: 150 - total flux cleaned = 0.0422781 Jy ! Component: 200 - total flux cleaned = 0.0481608 Jy ! Total flux subtracted in 200 components = 0.0481608 Jy ! Clean residual min=-0.015496 max=0.019100 Jy/beam ! Clean residual mean=-0.000000 rms=0.004156 Jy/beam ! Combined flux in latest and established models = 2.00222 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 255 components and 2.00222 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.109717Jy sigma=0.458681 ! Fit after self-cal, rms=0.109385Jy sigma=0.457768 ! Inverting map ! Added new window around map position (4, -53.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0114867 Jy ! Component: 100 - total flux cleaned = 0.0157717 Jy ! Component: 150 - total flux cleaned = 0.0192546 Jy ! Component: 200 - total flux cleaned = 0.0227454 Jy ! Total flux subtracted in 200 components = 0.0227454 Jy ! Clean residual min=-0.014629 max=0.013511 Jy/beam ! Clean residual mean=0.000017 rms=0.003618 Jy/beam ! Combined flux in latest and established models = 2.02497 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 273 components and 2.02497 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.105479Jy sigma=0.437681 ! Fit after self-cal, rms=0.105202Jy sigma=0.437281 ! 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. ! ! Fit before self-cal, rms=0.105202Jy sigma=0.437281 ! Fit after self-cal, rms=0.099540Jy sigma=0.426483 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.887 mas, bmaj=7.313 mas, bpa=-2.79 degrees ! Estimated noise=1.71395 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00900504 Jy ! Component: 100 - total flux cleaned = 0.0147617 Jy ! Component: 150 - total flux cleaned = 0.0189802 Jy ! Component: 200 - total flux cleaned = 0.0228064 Jy ! Total flux subtracted in 200 components = 0.0228064 Jy ! Clean residual min=-0.012164 max=0.011728 Jy/beam ! Clean residual mean=0.000008 rms=0.002963 Jy/beam ! Combined flux in latest and established models = 2.04777 Jy selfcal ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 290 components and 2.04777 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.096908Jy sigma=0.415741 ! Fit after self-cal, rms=0.096807Jy sigma=0.415510 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.00434866 Jy ! Component: 100 - total flux cleaned = 0.00694139 Jy ! Component: 150 - total flux cleaned = 0.00862039 Jy ! Component: 200 - total flux cleaned = 0.00975043 Jy ! Total flux subtracted in 200 components = 0.00975043 Jy ! Clean residual min=-0.011646 max=0.010818 Jy/beam ! Clean residual mean=0.000016 rms=0.002764 Jy/beam ! Combined flux in latest and established models = 2.05752 Jy ! Performing phase self-cal ! Adding 44 model components to the UV plane model. ! The established model now contains 309 components and 2.05753 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.095511Jy sigma=0.410038 ! Fit after self-cal, rms=0.095449Jy sigma=0.409905 ! 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. ! ! Fit before self-cal, rms=0.095449Jy sigma=0.409905 ! Fit after self-cal, rms=0.094160Jy sigma=0.408600 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=6.518 mas, bmaj=13.35 mas, bpa=-1.286 degrees ! Estimated noise=2.60762 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00902614 Jy ! Component: 100 - total flux cleaned = 0.0141843 Jy ! Component: 150 - total flux cleaned = 0.0173721 Jy ! Component: 200 - total flux cleaned = 0.0197485 Jy ! Total flux subtracted in 200 components = 0.0197485 Jy ! Clean residual min=-0.012237 max=0.012509 Jy/beam ! Clean residual mean=0.000037 rms=0.003336 Jy/beam ! Combined flux in latest and established models = 2.07727 Jy selfcal ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 329 components and 2.07727 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.093163Jy sigma=0.406846 ! Fit after self-cal, rms=0.092870Jy sigma=0.406149 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.00254997 Jy ! Component: 100 - total flux cleaned = 0.00329354 Jy ! Component: 150 - total flux cleaned = 0.00297172 Jy ! Component: 200 - total flux cleaned = 0.00283056 Jy ! Total flux subtracted in 200 components = 0.00283056 Jy ! Clean residual min=-0.011777 max=0.011632 Jy/beam ! Clean residual mean=0.000043 rms=0.003127 Jy/beam ! Combined flux in latest and established models = 2.0801 Jy ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 354 components and 2.0801 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.092715Jy sigma=0.405800 ! Fit after self-cal, rms=0.092660Jy sigma=0.405672 ! 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.877 mas, bmaj=7.307 mas, bpa=-2.95 degrees ! Estimated noise=1.70684 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.877 x 7.307 at -2.95 degrees (North through East) ! Clean map min=-0.015309 max=0.61688 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. ! ! Fit before self-cal, rms=0.092660Jy sigma=0.405672 ! Fit after self-cal, rms=0.089157Jy sigma=0.378813 ! Inverting map and beam ! Estimated beam: bmin=2.884 mas, bmaj=7.318 mas, bpa=-2.699 degrees ! Estimated noise=1.71491 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00240642 Jy ! Component: 100 - total flux cleaned = 0.00383246 Jy ! Component: 150 - total flux cleaned = 0.0051441 Jy ! Component: 200 - total flux cleaned = 0.0055347 Jy ! Total flux subtracted in 200 components = 0.0055347 Jy ! Clean residual min=-0.009941 max=0.009318 Jy/beam ! Clean residual mean=0.000006 rms=0.002300 Jy/beam ! Combined flux in latest and established models = 2.08564 Jy ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 372 components and 2.08564 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.088056Jy sigma=0.372922 ! Fit after self-cal, rms=0.088028Jy sigma=0.372792 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00167698 Jy ! Component: 100 - total flux cleaned = 0.00185375 Jy ! Component: 150 - total flux cleaned = 0.00167812 Jy ! Component: 200 - total flux cleaned = 0.00151336 Jy ! Total flux subtracted in 200 components = 0.00151336 Jy ! Clean residual min=-0.009544 max=0.008936 Jy/beam ! Clean residual mean=0.000008 rms=0.002198 Jy/beam ! Combined flux in latest and established models = 2.08715 Jy ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 381 components and 2.08715 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.087410Jy sigma=0.369525 ! Fit after self-cal, rms=0.087406Jy sigma=0.369447 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.087406Jy sigma=0.369447 ! Fit after self-cal, rms=0.086745Jy sigma=0.369070 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.49 mas, bmaj=6.709 mas, bpa=-3.506 degrees ! Estimated noise=3.01165 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.480687 Jy ! Component: 100 - total flux cleaned = 0.740285 Jy ! Component: 150 - total flux cleaned = 0.933566 Jy ! Component: 200 - total flux cleaned = 1.08031 Jy ! Total flux subtracted in 200 components = 1.08031 Jy ! Clean residual min=-0.053440 max=0.085227 Jy/beam ! Clean residual mean=-0.000434 rms=0.016393 Jy/beam ! Combined flux in latest and established models = 1.08031 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 46 model components to the UV plane model. ! The established model now contains 46 components and 1.08031 Jy ! Inverting map and beam ! Estimated beam: bmin=2.877 mas, bmaj=7.309 mas, bpa=-2.857 degrees ! Estimated noise=1.70895 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.170146 Jy ! Component: 100 - total flux cleaned = 0.293682 Jy ! Component: 150 - total flux cleaned = 0.392473 Jy ! Component: 200 - total flux cleaned = 0.474279 Jy ! Total flux subtracted in 200 components = 0.474279 Jy ! Clean residual min=-0.036444 max=0.054842 Jy/beam ! Clean residual mean=-0.000483 rms=0.011531 Jy/beam ! Combined flux in latest and established models = 1.55459 Jy ! Adding 81 model components to the UV plane model. ! The established model now contains 109 components and 1.55459 Jy ! Inverting map ! Added new window around map position (-9, -11). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0694466 Jy ! Component: 100 - total flux cleaned = 0.127654 Jy ! Component: 150 - total flux cleaned = 0.17739 Jy ! Component: 200 - total flux cleaned = 0.220197 Jy ! Total flux subtracted in 200 components = 0.220197 Jy ! Clean residual min=-0.023594 max=0.031642 Jy/beam ! Clean residual mean=-0.000267 rms=0.006957 Jy/beam ! Combined flux in latest and established models = 1.77478 Jy ! Adding 77 model components to the UV plane model. ! The established model now contains 168 components and 1.77478 Jy ! Inverting map ! Added new window around map position (-3.5, 0.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0382042 Jy ! Component: 100 - total flux cleaned = 0.0706989 Jy ! Component: 150 - total flux cleaned = 0.0977918 Jy ! Component: 200 - total flux cleaned = 0.120484 Jy ! Total flux subtracted in 200 components = 0.120484 Jy ! Clean residual min=-0.017186 max=0.016310 Jy/beam ! Clean residual mean=-0.000157 rms=0.004682 Jy/beam ! Combined flux in latest and established models = 1.89527 Jy ! Adding 69 model components to the UV plane model. ! The established model now contains 217 components and 1.89527 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=6.513 mas, bmaj=13.34 mas, bpa=-1.196 degrees ! Estimated noise=2.61201 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0481355 Jy ! Component: 100 - total flux cleaned = 0.0816082 Jy ! Component: 150 - total flux cleaned = 0.106429 Jy ! Component: 200 - total flux cleaned = 0.121775 Jy ! Total flux subtracted in 200 components = 0.121775 Jy ! Clean residual min=-0.013405 max=0.013282 Jy/beam ! Clean residual mean=-0.000140 rms=0.004497 Jy/beam ! Combined flux in latest and established models = 2.01704 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.0109676 Jy ! Component: 100 - total flux cleaned = 0.0207437 Jy ! Component: 150 - total flux cleaned = 0.0290745 Jy ! Component: 200 - total flux cleaned = 0.0348449 Jy ! Total flux subtracted in 200 components = 0.0348449 Jy ! Clean residual min=-0.011107 max=0.010415 Jy/beam ! Clean residual mean=-0.000058 rms=0.003263 Jy/beam ! Combined flux in latest and established models = 2.05189 Jy ! Adding 160 model components to the UV plane model. ! The established model now contains 356 components and 2.05189 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.877 mas, bmaj=7.309 mas, bpa=-2.857 degrees ! Estimated noise=1.70895 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00230653 Jy ! Component: 100 - total flux cleaned = 0.00528644 Jy ! Component: 150 - total flux cleaned = 0.00831489 Jy ! Component: 200 - total flux cleaned = 0.0113448 Jy ! Total flux subtracted in 200 components = 0.0113448 Jy ! Clean residual min=-0.011730 max=0.009274 Jy/beam ! Clean residual mean=-0.000015 rms=0.002434 Jy/beam ! Combined flux in latest and established models = 2.06323 Jy selfcal ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 387 components and 2.06323 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.088382Jy sigma=0.376603 ! Fit after self-cal, rms=0.087191Jy sigma=0.372819 wmodel J0735-1735_S_map.mod ! Writing 387 model components to file: J0735-1735_S_map.mod wobs J0735-1735_S_uvs.fits ! Writing UV FITS file: J0735-1735_S_uvs.fits wwins J0735-1735_S_map.win ! wwins: Wrote 8 windows to J0735-1735_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.00426328 Jy ! Component: 100 - total flux cleaned = -0.00704424 Jy ! Component: 150 - total flux cleaned = -0.0094588 Jy ! Component: 200 - total flux cleaned = -0.011578 Jy ! Component: 250 - total flux cleaned = -0.0131753 Jy ! Component: 300 - total flux cleaned = -0.0145714 Jy ! Component: 350 - total flux cleaned = -0.0150592 Jy ! Component: 400 - total flux cleaned = -0.0155258 Jy ! Component: 450 - total flux cleaned = -0.0154161 Jy ! Component: 500 - total flux cleaned = -0.0159684 Jy ! Component: 550 - total flux cleaned = -0.0162847 Jy ! Component: 600 - total flux cleaned = -0.0164872 Jy ! Component: 650 - total flux cleaned = -0.016686 Jy ! Component: 700 - total flux cleaned = -0.0168792 Jy ! Component: 750 - total flux cleaned = -0.0171617 Jy ! Component: 800 - total flux cleaned = -0.0172523 Jy ! Component: 850 - total flux cleaned = -0.0173412 Jy ! Component: 900 - total flux cleaned = -0.0175163 Jy ! Component: 950 - total flux cleaned = -0.0176855 Jy ! Component: 1000 - total flux cleaned = -0.0178498 Jy ! Component: 1050 - total flux cleaned = -0.0180938 Jy ! Component: 1100 - total flux cleaned = -0.0181731 Jy ! Component: 1150 - total flux cleaned = -0.0183288 Jy ! Component: 1200 - total flux cleaned = -0.0187813 Jy ! Component: 1250 - total flux cleaned = -0.0189287 Jy ! Component: 1300 - total flux cleaned = -0.0188565 Jy ! Component: 1350 - total flux cleaned = -0.0189995 Jy ! Component: 1400 - total flux cleaned = -0.0192765 Jy ! Component: 1450 - total flux cleaned = -0.019278 Jy ! Component: 1500 - total flux cleaned = -0.019346 Jy ! Component: 1550 - total flux cleaned = -0.0192839 Jy ! Component: 1600 - total flux cleaned = -0.019414 Jy ! Component: 1650 - total flux cleaned = -0.0192279 Jy ! Component: 1700 - total flux cleaned = -0.0193552 Jy ! Component: 1750 - total flux cleaned = -0.0194147 Jy ! Component: 1800 - total flux cleaned = -0.0193548 Jy ! Component: 1850 - total flux cleaned = -0.0193559 Jy ! Component: 1900 - total flux cleaned = -0.0193556 Jy ! Component: 1950 - total flux cleaned = -0.0194113 Jy ! Component: 2000 - total flux cleaned = -0.0195212 Jy ! Component: 2050 - total flux cleaned = -0.0193566 Jy ! Component: 2100 - total flux cleaned = -0.0193566 Jy ! Component: 2150 - total flux cleaned = -0.0193557 Jy ! Component: 2200 - total flux cleaned = -0.0191476 Jy ! Component: 2250 - total flux cleaned = -0.0189432 Jy ! Component: 2300 - total flux cleaned = -0.0187912 Jy ! Component: 2350 - total flux cleaned = -0.0187415 Jy ! Component: 2400 - total flux cleaned = -0.0183514 Jy ! Component: 2450 - total flux cleaned = -0.0184476 Jy ! Component: 2500 - total flux cleaned = -0.0183047 Jy ! Component: 2550 - total flux cleaned = -0.0178843 Jy ! Component: 2600 - total flux cleaned = -0.0179773 Jy ! Component: 2650 - total flux cleaned = -0.0178412 Jy ! Component: 2700 - total flux cleaned = -0.0177527 Jy ! Component: 2750 - total flux cleaned = -0.0177086 Jy ! Component: 2800 - total flux cleaned = -0.0177955 Jy ! Component: 2850 - total flux cleaned = -0.0178377 Jy ! Component: 2900 - total flux cleaned = -0.0177532 Jy ! Component: 2950 - total flux cleaned = -0.0175028 Jy ! Component: 3000 - total flux cleaned = -0.0176268 Jy ! Component: 3050 - total flux cleaned = -0.0175865 Jy ! Component: 3100 - total flux cleaned = -0.0174661 Jy ! Component: 3150 - total flux cleaned = -0.0174662 Jy ! Component: 3200 - total flux cleaned = -0.0175036 Jy ! Component: 3250 - total flux cleaned = -0.0174648 Jy ! Component: 3300 - total flux cleaned = -0.0175412 Jy ! Component: 3350 - total flux cleaned = -0.0177286 Jy ! Component: 3400 - total flux cleaned = -0.0176537 Jy ! Component: 3450 - total flux cleaned = -0.017654 Jy ! Component: 3500 - total flux cleaned = -0.0177266 Jy ! Component: 3550 - total flux cleaned = -0.0176186 Jy ! Component: 3600 - total flux cleaned = -0.0177604 Jy ! Component: 3650 - total flux cleaned = -0.0176901 Jy ! Component: 3700 - total flux cleaned = -0.0176213 Jy ! Component: 3750 - total flux cleaned = -0.01769 Jy ! Component: 3800 - total flux cleaned = -0.0176219 Jy ! Component: 3850 - total flux cleaned = -0.0174545 Jy ! Component: 3900 - total flux cleaned = -0.0175211 Jy ! Component: 3950 - total flux cleaned = -0.0174562 Jy ! Component: 4000 - total flux cleaned = -0.0175209 Jy ! Component: 4050 - total flux cleaned = -0.0174571 Jy ! Total flux subtracted in 4096 components = -0.0176784 Jy ! Clean residual min=-0.001667 max=0.001749 Jy/beam ! Clean residual mean=-0.000009 rms=0.000590 Jy/beam ! Combined flux in latest and established models = 2.04555 Jy keep ! Adding 336 model components to the UV plane model. ! The established model now contains 722 components and 2.04555 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /image_prod/ug002_uva//J0735-1735_S_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.877 x 7.309 at -2.857 degrees (North through East) ! Clean map min=-0.0081891 max=0.62804 Jy/beam ! Writing clean map to FITS file: /image_prod/ug002_uva//J0735-1735_S_map.fits wmodel /image_prod/ug002_uva//J0735-1735_S_map.mod ! Writing 722 model components to file: /image_prod/ug002_uva//J0735-1735_S_map.mod wobs /image_prod/ug002_uva//J0735-1735_S_uvs.fits ! Writing UV FITS file: /image_prod/ug002_uva//J0735-1735_S_uvs.fits wwins /image_prod/ug002_uva//J0735-1735_S_map.win ! wwins: Wrote 9 windows to /image_prod/ug002_uva//J0735-1735_S_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.628042 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000579928 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.277017 0.277017 0.554034 1.10807 2.21614 4.43228 8.86455 17.7291 35.4582 70.9164 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0017 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.628 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 1082.96 quit ! Quitting program ! Log file /image_prod/ug002_uva//J0735-1735_S_dfm.log closed on Wed Mar 4 21:53:56 2020