! Started logfile: /vlbi/uf001/uva//J0625+1440_S_dfm.log on Sun Jan 13 22:59:37 2019 float field_size field_size = 1024 float field_cell field_cell = 0.5000 float taper_size taper_size = 20.0 float freq freq = 2.22000 obs /vlbi/uf001/uva//J0625+1440_S_uva.fits ! Reading UV FITS file: /vlbi/uf001/uva//J0625+1440_S_uva.fits ! AN table 1: 72 integrations on 45 of 45 possible baselines. ! Apparent sampling: 0.93179 visibilities/baseline/integration-bin. ! Found source: J0625+1440 ! ! 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 9057 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 J0625+1440_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 48 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 48 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 48 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=1.048783Jy sigma=20.565366 ! Fit after self-cal, rms=0.838768Jy sigma=16.852044 ! 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=3.068 mas, bmaj=6.247 mas, bpa=-3.371 degrees ! Estimated noise=0.534556 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.127691 Jy ! Component: 100 - total flux cleaned = 0.195701 Jy ! Total flux subtracted in 100 components = 0.195701 Jy ! Clean residual min=-0.021505 max=0.045391 Jy/beam ! Clean residual mean=-0.000229 rms=0.007159 Jy/beam ! Combined flux in latest and established models = 0.195701 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 27 components and 0.195701 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.234161Jy sigma=3.824002 ! Fit after self-cal, rms=0.233013Jy sigma=3.799823 ! Inverting map ! Added new window around map position (4.5, -1.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0485162 Jy ! Component: 100 - total flux cleaned = 0.0835908 Jy ! Total flux subtracted in 100 components = 0.0835908 Jy ! Clean residual min=-0.019596 max=0.034333 Jy/beam ! Clean residual mean=-0.000210 rms=0.006331 Jy/beam ! Combined flux in latest and established models = 0.279291 Jy ! Performing phase self-cal ! Adding 48 model components to the UV plane model. ! The established model now contains 71 components and 0.279291 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.207925Jy sigma=3.372954 ! Fit after self-cal, rms=0.207567Jy sigma=3.364544 ! 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=3.516 mas, bmaj=6.991 mas, bpa=-2.222 degrees ! Estimated noise=0.397746 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.0578481 Jy ! Component: 100 - total flux cleaned = 0.100161 Jy ! Component: 150 - total flux cleaned = 0.132608 Jy ! Component: 200 - total flux cleaned = 0.157816 Jy ! Total flux subtracted in 200 components = 0.157816 Jy ! Clean residual min=-0.024664 max=0.040646 Jy/beam ! Clean residual mean=-0.000228 rms=0.009235 Jy/beam ! Combined flux in latest and established models = 0.437107 Jy ! Performing phase self-cal ! Adding 84 model components to the UV plane model. ! The established model now contains 140 components and 0.437107 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.176738Jy sigma=2.926778 ! Fit after self-cal, rms=0.175994Jy sigma=2.907205 ! Inverting map print "============== Finished natural weighting clean ==============" ! ============== Finished natural weighting clean ============== ! restore ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.516 x 6.991 at -2.222 degrees (North through East) ! Clean map min=-0.023506 max=0.15055 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.17 FD 1.09 HN 1.34 KP 1.26 ! LA 0.66 MK 1.32 NL 1.45 OV 1.28 ! PT 0.64 SC 1.38 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 1.16 FD 1.00 HN 1.17 KP 1.07 ! LA 0.74 MK 1.37 NL 1.33 OV 1.13 ! PT 0.73 SC 1.28 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.27 FD 1.01 HN 1.07 KP 1.03 ! LA 0.76 MK 1.12 NL 1.26 OV 1.19 ! PT 0.81 SC 1.21 ! ! ! Fit before self-cal, rms=0.175994Jy sigma=2.907205 ! Fit after self-cal, rms=0.159821Jy sigma=2.592112 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.768 mas, bmaj=7.526 mas, bpa=-2.479 degrees ! Estimated noise=0.444031 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 (-20, -13). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0405033 Jy ! Component: 100 - total flux cleaned = 0.0614443 Jy ! Component: 150 - total flux cleaned = 0.0767127 Jy ! Component: 200 - total flux cleaned = 0.0886117 Jy ! Total flux subtracted in 200 components = 0.0886117 Jy ! Clean residual min=-0.021864 max=0.034520 Jy/beam ! Clean residual mean=-0.000054 rms=0.007150 Jy/beam ! Combined flux in latest and established models = 0.525719 Jy ! Performing phase self-cal ! Adding 58 model components to the UV plane model. ! The established model now contains 193 components and 0.525719 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.148924Jy sigma=2.351739 ! Fit after self-cal, rms=0.148416Jy sigma=2.322891 ! 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.148416Jy sigma=2.322891 ! Fit after self-cal, rms=0.147843Jy sigma=2.229549 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.829 mas, bmaj=7.564 mas, bpa=-1.106 degrees ! Estimated noise=0.437458 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0158887 Jy ! Component: 100 - total flux cleaned = 0.0264741 Jy ! Component: 150 - total flux cleaned = 0.0350519 Jy ! Component: 200 - total flux cleaned = 0.0422081 Jy ! Total flux subtracted in 200 components = 0.0422081 Jy ! Clean residual min=-0.021453 max=0.031280 Jy/beam ! Clean residual mean=-0.000066 rms=0.006326 Jy/beam ! Combined flux in latest and established models = 0.567927 Jy selfcal ! Performing phase self-cal ! Adding 69 model components to the UV plane model. ! The established model now contains 238 components and 0.567927 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.145326Jy sigma=2.179431 ! Fit after self-cal, rms=0.145221Jy sigma=2.173670 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 (3, -28). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0266454 Jy ! Component: 100 - total flux cleaned = 0.0388188 Jy ! Component: 150 - total flux cleaned = 0.0446794 Jy ! Component: 200 - total flux cleaned = 0.0490461 Jy ! Total flux subtracted in 200 components = 0.0490461 Jy ! Clean residual min=-0.018858 max=0.026211 Jy/beam ! Clean residual mean=-0.000055 rms=0.005754 Jy/beam ! Combined flux in latest and established models = 0.616973 Jy ! Performing phase self-cal ! Adding 59 model components to the UV plane model. ! The established model now contains 294 components and 0.616973 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.139643Jy sigma=2.067871 ! Fit after self-cal, rms=0.139147Jy sigma=2.056481 ! Inverting map ! Added new window around map position (79.5, 55.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0229749 Jy ! Component: 100 - total flux cleaned = 0.0322614 Jy ! Component: 150 - total flux cleaned = 0.0376845 Jy ! Component: 200 - total flux cleaned = 0.040423 Jy ! Total flux subtracted in 200 components = 0.040423 Jy ! Clean residual min=-0.017928 max=0.022989 Jy/beam ! Clean residual mean=-0.000073 rms=0.005181 Jy/beam ! Combined flux in latest and established models = 0.657396 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 329 components and 0.657396 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.134257Jy sigma=1.970541 ! Fit after self-cal, rms=0.133853Jy sigma=1.958670 ! Inverting map ! Added new window around map position (-1.5, 27). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0212141 Jy ! Component: 100 - total flux cleaned = 0.0326656 Jy ! Component: 150 - total flux cleaned = 0.0415237 Jy ! Component: 200 - total flux cleaned = 0.0484218 Jy ! Total flux subtracted in 200 components = 0.0484218 Jy ! Clean residual min=-0.015782 max=0.020062 Jy/beam ! Clean residual mean=-0.000059 rms=0.004616 Jy/beam ! Combined flux in latest and established models = 0.705818 Jy ! Performing phase self-cal ! Adding 48 model components to the UV plane model. ! The established model now contains 372 components and 0.705818 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.129871Jy sigma=1.871934 ! Fit after self-cal, rms=0.129163Jy sigma=1.864654 ! Inverting map ! Added new window around map position (120.5, 25). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0179611 Jy ! Component: 100 - total flux cleaned = 0.0283878 Jy ! Component: 150 - total flux cleaned = 0.0364899 Jy ! Component: 200 - total flux cleaned = 0.0415189 Jy ! Total flux subtracted in 200 components = 0.0415189 Jy ! Clean residual min=-0.016094 max=0.016318 Jy/beam ! Clean residual mean=-0.000104 rms=0.004156 Jy/beam ! Combined flux in latest and established models = 0.747337 Jy ! Performing phase self-cal ! Adding 68 model components to the UV plane model. ! The established model now contains 428 components and 0.747337 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.126075Jy sigma=1.795764 ! Fit after self-cal, rms=0.125642Jy sigma=1.788931 ! Inverting map ! Added new window around map position (101, 72). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0095647 Jy ! Component: 100 - total flux cleaned = 0.0111708 Jy ! Component: 150 - total flux cleaned = 0.0121182 Jy ! Component: 200 - total flux cleaned = 0.0127037 Jy ! Total flux subtracted in 200 components = 0.0127037 Jy ! Clean residual min=-0.014132 max=0.014986 Jy/beam ! Clean residual mean=-0.000090 rms=0.003854 Jy/beam ! Combined flux in latest and established models = 0.76004 Jy ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 458 components and 0.76004 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.122303Jy sigma=1.743383 ! Fit after self-cal, rms=0.122116Jy sigma=1.737888 ! 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.122116Jy sigma=1.737888 ! Fit after self-cal, rms=0.120753Jy sigma=1.653015 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.959 mas, bmaj=7.657 mas, bpa=0.843 degrees ! Estimated noise=0.441655 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00969135 Jy ! Component: 100 - total flux cleaned = 0.0145405 Jy ! Component: 150 - total flux cleaned = 0.0193734 Jy ! Component: 200 - total flux cleaned = 0.0224725 Jy ! Total flux subtracted in 200 components = 0.0224725 Jy ! Clean residual min=-0.012844 max=0.014301 Jy/beam ! Clean residual mean=-0.000075 rms=0.003524 Jy/beam ! Combined flux in latest and established models = 0.782513 Jy selfcal ! Performing phase self-cal ! Adding 49 model components to the UV plane model. ! The established model now contains 495 components and 0.782513 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.118987Jy sigma=1.621764 ! Fit after self-cal, rms=0.118982Jy sigma=1.619190 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.00317354 Jy ! Component: 100 - total flux cleaned = 0.00630106 Jy ! Component: 150 - total flux cleaned = 0.00855438 Jy ! Component: 200 - total flux cleaned = 0.0102615 Jy ! Total flux subtracted in 200 components = 0.0102615 Jy ! Clean residual min=-0.011937 max=0.013921 Jy/beam ! Clean residual mean=-0.000074 rms=0.003402 Jy/beam ! Combined flux in latest and established models = 0.792774 Jy ! Performing phase self-cal ! Adding 44 model components to the UV plane model. ! The established model now contains 507 components and 0.792774 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.117928Jy sigma=1.603534 ! Fit after self-cal, rms=0.117963Jy sigma=1.602605 ! 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.117963Jy sigma=1.602605 ! Fit after self-cal, rms=0.118346Jy sigma=1.597909 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.851 mas, bmaj=11.12 mas, bpa=-1.816 degrees ! Estimated noise=0.499134 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00872151 Jy ! Component: 100 - total flux cleaned = 0.0145534 Jy ! Component: 150 - total flux cleaned = 0.0179006 Jy ! Component: 200 - total flux cleaned = 0.0199034 Jy ! Total flux subtracted in 200 components = 0.0199034 Jy ! Clean residual min=-0.014880 max=0.016461 Jy/beam ! Clean residual mean=-0.000129 rms=0.004386 Jy/beam ! Combined flux in latest and established models = 0.812678 Jy selfcal ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 522 components and 0.812678 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.117534Jy sigma=1.581743 ! Fit after self-cal, rms=0.117459Jy sigma=1.580369 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.00311633 Jy ! Component: 100 - total flux cleaned = 0.00379321 Jy ! Component: 150 - total flux cleaned = 0.00443693 Jy ! Component: 200 - total flux cleaned = 0.00545654 Jy ! Total flux subtracted in 200 components = 0.00545654 Jy ! Clean residual min=-0.014363 max=0.015993 Jy/beam ! Clean residual mean=-0.000120 rms=0.004252 Jy/beam ! Combined flux in latest and established models = 0.818134 Jy ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 534 components and 0.818134 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.116993Jy sigma=1.572701 ! Fit after self-cal, rms=0.116924Jy sigma=1.572147 ! 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=3.964 mas, bmaj=7.656 mas, bpa=1.297 degrees ! Estimated noise=0.440564 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.964 x 7.656 at 1.297 degrees (North through East) ! Clean map min=-0.012396 max=0.17128 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 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=3.412 mas, bmaj=6.702 mas, bpa=0.09006 degrees ! Estimated noise=0.63853 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.145073 Jy ! Component: 100 - total flux cleaned = 0.220618 Jy ! Component: 150 - total flux cleaned = 0.271193 Jy ! Component: 200 - total flux cleaned = 0.312297 Jy ! Total flux subtracted in 200 components = 0.312297 Jy ! Clean residual min=-0.020589 max=0.028228 Jy/beam ! Clean residual mean=-0.000189 rms=0.006943 Jy/beam ! Combined flux in latest and established models = 0.312297 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 63 model components to the UV plane model. ! The established model now contains 63 components and 0.312297 Jy ! Inverting map and beam ! Estimated beam: bmin=3.964 mas, bmaj=7.656 mas, bpa=1.297 degrees ! Estimated noise=0.440564 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.0550846 Jy ! Component: 100 - total flux cleaned = 0.0970432 Jy ! Component: 150 - total flux cleaned = 0.132753 Jy ! Component: 200 - total flux cleaned = 0.164393 Jy ! Total flux subtracted in 200 components = 0.164393 Jy ! Clean residual min=-0.021461 max=0.026414 Jy/beam ! Clean residual mean=-0.000171 rms=0.007169 Jy/beam ! Combined flux in latest and established models = 0.476689 Jy ! Adding 87 model components to the UV plane model. ! The established model now contains 139 components and 0.476689 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.851 mas, bmaj=11.12 mas, bpa=-1.816 degrees ! Estimated noise=0.499134 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0466407 Jy ! Component: 100 - total flux cleaned = 0.0841908 Jy ! Component: 150 - total flux cleaned = 0.116262 Jy ! Component: 200 - total flux cleaned = 0.144367 Jy ! Total flux subtracted in 200 components = 0.144367 Jy ! Clean residual min=-0.023939 max=0.025972 Jy/beam ! Clean residual mean=-0.000263 rms=0.008047 Jy/beam ! Combined flux in latest and established models = 0.621056 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.0249147 Jy ! Component: 100 - total flux cleaned = 0.0462327 Jy ! Component: 150 - total flux cleaned = 0.0637195 Jy ! Component: 200 - total flux cleaned = 0.0787659 Jy ! Total flux subtracted in 200 components = 0.0787659 Jy ! Clean residual min=-0.020328 max=0.020660 Jy/beam ! Clean residual mean=-0.000235 rms=0.006210 Jy/beam ! Combined flux in latest and established models = 0.699822 Jy ! Adding 122 model components to the UV plane model. ! The established model now contains 238 components and 0.699822 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.964 mas, bmaj=7.656 mas, bpa=1.297 degrees ! Estimated noise=0.440564 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0130189 Jy ! Component: 100 - total flux cleaned = 0.0213137 Jy ! Component: 150 - total flux cleaned = 0.0292869 Jy ! Component: 200 - total flux cleaned = 0.0350999 Jy ! Total flux subtracted in 200 components = 0.0350999 Jy ! Clean residual min=-0.015422 max=0.015242 Jy/beam ! Clean residual mean=-0.000113 rms=0.003872 Jy/beam ! Combined flux in latest and established models = 0.734922 Jy selfcal ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 274 components and 0.734922 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.122576Jy sigma=1.687181 ! Fit after self-cal, rms=0.122055Jy sigma=1.676990 wmodel J0625+1440_S_map.mod ! Writing 274 model components to file: J0625+1440_S_map.mod wobs J0625+1440_S_uvs.fits ! Writing UV FITS file: J0625+1440_S_uvs.fits wwins J0625+1440_S_map.win ! wwins: Wrote 8 windows to J0625+1440_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.00031978 Jy ! Component: 100 - total flux cleaned = -6.053e-05 Jy ! Component: 150 - total flux cleaned = -0.000968435 Jy ! Component: 200 - total flux cleaned = -0.00140346 Jy ! Component: 250 - total flux cleaned = -0.00161241 Jy ! Component: 300 - total flux cleaned = -0.00181147 Jy ! Component: 350 - total flux cleaned = -0.00220357 Jy ! Component: 400 - total flux cleaned = -0.00220564 Jy ! Component: 450 - total flux cleaned = -0.00239245 Jy ! Component: 500 - total flux cleaned = -0.00221561 Jy ! Component: 550 - total flux cleaned = -0.00221954 Jy ! Component: 600 - total flux cleaned = -0.00205047 Jy ! Component: 650 - total flux cleaned = -0.00188449 Jy ! Component: 700 - total flux cleaned = -0.00171724 Jy ! Component: 750 - total flux cleaned = -0.00107294 Jy ! Component: 800 - total flux cleaned = -0.000915993 Jy ! Component: 850 - total flux cleaned = -0.000451072 Jy ! Component: 900 - total flux cleaned = 5.38171e-06 Jy ! Component: 950 - total flux cleaned = 0.000749845 Jy ! Component: 1000 - total flux cleaned = 0.00103935 Jy ! Component: 1050 - total flux cleaned = 0.00175245 Jy ! Component: 1100 - total flux cleaned = 0.00245653 Jy ! Component: 1150 - total flux cleaned = 0.00328472 Jy ! Component: 1200 - total flux cleaned = 0.00369731 Jy ! Component: 1250 - total flux cleaned = 0.0046342 Jy ! Component: 1300 - total flux cleaned = 0.00489887 Jy ! Component: 1350 - total flux cleaned = 0.00554934 Jy ! Component: 1400 - total flux cleaned = 0.00580259 Jy ! Component: 1450 - total flux cleaned = 0.00580323 Jy ! Component: 1500 - total flux cleaned = 0.0061757 Jy ! Component: 1550 - total flux cleaned = 0.00642276 Jy ! Component: 1600 - total flux cleaned = 0.0067866 Jy ! Component: 1650 - total flux cleaned = 0.00726424 Jy ! Component: 1700 - total flux cleaned = 0.00785084 Jy ! Component: 1750 - total flux cleaned = 0.00819893 Jy ! Component: 1800 - total flux cleaned = 0.00865532 Jy ! Component: 1850 - total flux cleaned = 0.00933369 Jy ! Component: 1900 - total flux cleaned = 0.0102247 Jy ! Component: 1950 - total flux cleaned = 0.0104422 Jy ! Component: 2000 - total flux cleaned = 0.0109858 Jy ! Component: 2050 - total flux cleaned = 0.0118451 Jy ! Component: 2100 - total flux cleaned = 0.012586 Jy ! Component: 2150 - total flux cleaned = 0.0133187 Jy ! Component: 2200 - total flux cleaned = 0.0138353 Jy ! Component: 2250 - total flux cleaned = 0.0143459 Jy ! Component: 2300 - total flux cleaned = 0.0147496 Jy ! Component: 2350 - total flux cleaned = 0.0153494 Jy ! Component: 2400 - total flux cleaned = 0.0158424 Jy ! Component: 2450 - total flux cleaned = 0.0162328 Jy ! Component: 2500 - total flux cleaned = 0.0165229 Jy ! Component: 2550 - total flux cleaned = 0.0171001 Jy ! Component: 2600 - total flux cleaned = 0.0174788 Jy ! Component: 2650 - total flux cleaned = 0.0179471 Jy ! Component: 2700 - total flux cleaned = 0.0181327 Jy ! Component: 2750 - total flux cleaned = 0.0187765 Jy ! Component: 2800 - total flux cleaned = 0.0193232 Jy ! Component: 2850 - total flux cleaned = 0.0197722 Jy ! Component: 2900 - total flux cleaned = 0.0200397 Jy ! Component: 2950 - total flux cleaned = 0.0204839 Jy ! Component: 3000 - total flux cleaned = 0.0210086 Jy ! Component: 3050 - total flux cleaned = 0.0212693 Jy ! Component: 3100 - total flux cleaned = 0.02204 Jy ! Component: 3150 - total flux cleaned = 0.0222945 Jy ! Component: 3200 - total flux cleaned = 0.0223793 Jy ! Component: 3250 - total flux cleaned = 0.0232114 Jy ! Component: 3300 - total flux cleaned = 0.023623 Jy ! Component: 3350 - total flux cleaned = 0.0238685 Jy ! Component: 3400 - total flux cleaned = 0.024515 Jy ! Component: 3450 - total flux cleaned = 0.0249153 Jy ! Component: 3500 - total flux cleaned = 0.0256301 Jy ! Component: 3550 - total flux cleaned = 0.0261803 Jy ! Component: 3600 - total flux cleaned = 0.0266464 Jy ! Component: 3650 - total flux cleaned = 0.0269554 Jy ! Component: 3700 - total flux cleaned = 0.0276449 Jy ! Component: 3750 - total flux cleaned = 0.0281001 Jy ! Component: 3800 - total flux cleaned = 0.028552 Jy ! Component: 3850 - total flux cleaned = 0.0287759 Jy ! Component: 3900 - total flux cleaned = 0.0293687 Jy ! Component: 3950 - total flux cleaned = 0.0298098 Jy ! Component: 4000 - total flux cleaned = 0.0302475 Jy ! Component: 4050 - total flux cleaned = 0.0308249 Jy ! Total flux subtracted in 4096 components = 0.0309671 Jy ! Clean residual min=-0.004480 max=0.004177 Jy/beam ! Clean residual mean=-0.000025 rms=0.001487 Jy/beam ! Combined flux in latest and established models = 0.765889 Jy keep ! Adding 388 model components to the UV plane model. ! The established model now contains 661 components and 0.765889 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /vlbi/uf001/uva//J0625+1440_S_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.964 x 7.656 at 1.297 degrees (North through East) ! Clean map min=-0.011313 max=0.1738 Jy/beam ! Writing clean map to FITS file: /vlbi/uf001/uva//J0625+1440_S_map.fits wmodel /vlbi/uf001/uva//J0625+1440_S_map.mod ! Writing 661 model components to file: /vlbi/uf001/uva//J0625+1440_S_map.mod wobs /vlbi/uf001/uva//J0625+1440_S_uvs.fits ! Writing UV FITS file: /vlbi/uf001/uva//J0625+1440_S_uvs.fits wwins /vlbi/uf001/uva//J0625+1440_S_map.win ! wwins: Wrote 9 windows to /vlbi/uf001/uva//J0625+1440_S_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.1738 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00148136 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -2.55702 2.55702 5.11403 10.2281 20.4561 40.9123 81.8245 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0044 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.173 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 117.324 quit ! Quitting program ! Log file /vlbi/uf001/uva//J0625+1440_S_dfm.log closed on Sun Jan 13 22:59:42 2019