! Started logfile: /vlbi/uf001/uvm//J1955+5131_S_b_dfm.log on Sun Jan 13 22:23:59 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/uvm//J1955+5131_S_b_uvm.fits ! Reading UV FITS file: /vlbi/uf001/uvm//J1955+5131_S_b_uvm.fits ! AN table 1: 5 integrations on 36 of 36 possible baselines. ! AN table 2: 19 integrations on 28 of 28 possible baselines. ! AN table 3: 14 integrations on 36 of 36 possible baselines. ! Apparent sampling: 0.714638 visibilities/baseline/integration-bin. ! Found source: J1955+5131 ! ! 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 2607 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 J1955+5131_S_b] 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 12 telescope corrections were flagged in sub-array 1. ! A total of 59 telescope corrections were flagged in sub-array 2. ! A total of 75 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 2. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 59 telescope corrections were flagged in sub-array 2. ! A total of 75 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 3. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 59 telescope corrections were flagged in sub-array 2. ! A total of 75 telescope corrections were flagged in sub-array 3. ! ! Fit before self-cal, rms=1.285840Jy sigma=21.259698 ! Fit after self-cal, rms=0.301680Jy sigma=5.323133 ! 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.001 mas, bmaj=4.729 mas, bpa=33.89 degrees ! Estimated noise=1.28094 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.553971 Jy ! Component: 100 - total flux cleaned = 0.694838 Jy ! Total flux subtracted in 100 components = 0.694838 Jy ! Clean residual min=-0.029062 max=0.059586 Jy/beam ! Clean residual mean=-0.000012 rms=0.007915 Jy/beam ! Combined flux in latest and established models = 0.694838 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 11 components and 0.694838 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.197427Jy sigma=2.947904 ! Fit after self-cal, rms=0.197252Jy sigma=2.944199 ! Inverting map ! Added new window around map position (-4, 3.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0686275 Jy ! Component: 100 - total flux cleaned = 0.109617 Jy ! Total flux subtracted in 100 components = 0.109617 Jy ! Clean residual min=-0.024602 max=0.048175 Jy/beam ! Clean residual mean=-0.000009 rms=0.006315 Jy/beam ! Combined flux in latest and established models = 0.804455 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 31 components and 0.804455 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.153038Jy sigma=2.362450 ! Fit after self-cal, rms=0.152118Jy sigma=2.341342 ! Inverting map ! Added new window around map position (-12.5, 4.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0441448 Jy ! Component: 100 - total flux cleaned = 0.0730117 Jy ! Total flux subtracted in 100 components = 0.0730117 Jy ! Clean residual min=-0.021575 max=0.024624 Jy/beam ! Clean residual mean=-0.000007 rms=0.005272 Jy/beam ! Combined flux in latest and established models = 0.877466 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 54 components and 0.877466 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.130374Jy sigma=2.042266 ! Fit after self-cal, rms=0.129544Jy sigma=2.025981 ! 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=4.895 mas, bmaj=5.564 mas, bpa=30.55 degrees ! Estimated noise=1.07855 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.0366887 Jy ! Component: 100 - total flux cleaned = 0.0617941 Jy ! Component: 150 - total flux cleaned = 0.0801625 Jy ! Component: 200 - total flux cleaned = 0.0940097 Jy ! Total flux subtracted in 200 components = 0.0940097 Jy ! Clean residual min=-0.019038 max=0.020264 Jy/beam ! Clean residual mean=-0.000060 rms=0.004813 Jy/beam ! Combined flux in latest and established models = 0.971476 Jy ! Performing phase self-cal ! Adding 59 model components to the UV plane model. ! The established model now contains 92 components and 0.971476 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.112661Jy sigma=1.832925 ! Fit after self-cal, rms=0.111993Jy sigma=1.819869 ! Inverting map print "============== Finished natural weighting clean ==============" ! ============== Finished natural weighting clean ============== ! restore ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 4.895 x 5.564 at 30.55 degrees (North through East) ! Clean map min=-0.017066 max=0.72954 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.97 HN 1.34 KP 0.97 LA 0.94 ! MK 0.94 NL 0.97 OV 0.95 PT 0.96 ! SC 1.00 ! ! A total of 9 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.03 HN 0.95 KP 0.98 ! LA 0.95 NL 1.08 OV 0.97 PT 1.00 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.03 FD 1.01 HN 1.01 KP 1.03 ! LA 0.91 MK 0.98 NL 1.15 OV 1.01 ! PT 1.02 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 1.05 HN 1.37 KP 0.98 LA 0.96 ! MK 1.01 NL 0.96 OV 0.77 PT 0.94 ! SC 0.99 ! ! A total of 9 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.02 HN 0.95 KP 0.98 ! LA 0.94 NL 1.06 OV 0.90 PT 1.00 ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.78 FD 1.04 HN 1.00 KP 1.02 ! LA 0.93 MK 1.03 NL 1.18 OV 0.96 ! PT 1.04 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.01 HN 1.39 KP 0.99 LA 0.91 ! MK 1.04 NL 0.95 OV 0.94 PT 0.94 ! SC 0.97 ! ! A total of 9 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.01 HN 0.98 KP 0.97 ! LA 0.95 NL 1.04 OV 0.96 PT 1.01 ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.94 FD 1.02 HN 1.03 KP 1.04 ! LA 0.90 MK 1.04 NL 1.25 OV 1.07 ! PT 1.05 ! ! ! Fit before self-cal, rms=0.111993Jy sigma=1.819869 ! Fit after self-cal, rms=0.074470Jy sigma=1.206212 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=4.911 mas, bmaj=5.57 mas, bpa=33 degrees ! Estimated noise=1.08146 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.00676392 Jy ! Component: 100 - total flux cleaned = 0.00748977 Jy ! Component: 150 - total flux cleaned = 0.00732449 Jy ! Component: 200 - total flux cleaned = 0.00731955 Jy ! Total flux subtracted in 200 components = 0.00731955 Jy ! Clean residual min=-0.010341 max=0.010904 Jy/beam ! Clean residual mean=-0.000014 rms=0.002530 Jy/beam ! Combined flux in latest and established models = 0.978796 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 104 components and 0.978796 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.073079Jy sigma=1.184694 ! Fit after self-cal, rms=0.073010Jy sigma=1.182135 ! 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.073010Jy sigma=1.182135 ! Fit after self-cal, rms=0.064108Jy sigma=1.003836 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=4.908 mas, bmaj=5.566 mas, bpa=32.14 degrees ! Estimated noise=1.07892 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00367033 Jy ! Component: 100 - total flux cleaned = 0.0033459 Jy ! Component: 150 - total flux cleaned = 0.0033475 Jy ! Component: 200 - total flux cleaned = 0.00311851 Jy ! Total flux subtracted in 200 components = 0.00311851 Jy ! Clean residual min=-0.007959 max=0.008774 Jy/beam ! Clean residual mean=0.000003 rms=0.001672 Jy/beam ! Combined flux in latest and established models = 0.981914 Jy selfcal ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 117 components and 0.981914 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.063307Jy sigma=0.992065 ! Fit after self-cal, rms=0.063260Jy sigma=0.990772 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 (-37.5, 19). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00803601 Jy ! Component: 100 - total flux cleaned = 0.0127139 Jy ! Component: 150 - total flux cleaned = 0.0159766 Jy ! Component: 200 - total flux cleaned = 0.017924 Jy ! Total flux subtracted in 200 components = 0.017924 Jy ! Clean residual min=-0.007175 max=0.008251 Jy/beam ! Clean residual mean=0.000008 rms=0.001491 Jy/beam ! Combined flux in latest and established models = 0.999838 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 131 components and 0.999838 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.061929Jy sigma=0.971646 ! Fit after self-cal, rms=0.061837Jy sigma=0.970463 ! Inverting map ! Added new window around map position (-20, 9). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00728555 Jy ! Component: 100 - total flux cleaned = 0.0113281 Jy ! Component: 150 - total flux cleaned = 0.0144155 Jy ! Component: 200 - total flux cleaned = 0.0162923 Jy ! Total flux subtracted in 200 components = 0.0162923 Jy ! Clean residual min=-0.006788 max=0.005425 Jy/beam ! Clean residual mean=0.000019 rms=0.001317 Jy/beam ! Combined flux in latest and established models = 1.01613 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 149 components and 1.01613 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.060760Jy sigma=0.954347 ! Fit after self-cal, rms=0.060654Jy sigma=0.952934 ! 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.060654Jy sigma=0.952934 ! Fit after self-cal, rms=0.060325Jy sigma=0.949227 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=4.904 mas, bmaj=5.567 mas, bpa=32.4 degrees ! Estimated noise=1.07727 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0023761 Jy ! Component: 100 - total flux cleaned = 0.00402174 Jy ! Component: 150 - total flux cleaned = 0.00505337 Jy ! Component: 200 - total flux cleaned = 0.00608126 Jy ! Total flux subtracted in 200 components = 0.00608126 Jy ! Clean residual min=-0.005909 max=0.005091 Jy/beam ! Clean residual mean=0.000021 rms=0.001214 Jy/beam ! Combined flux in latest and established models = 1.02221 Jy selfcal ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 165 components and 1.02221 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.059897Jy sigma=0.942939 ! Fit after self-cal, rms=0.059871Jy sigma=0.942592 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.000928484 Jy ! Component: 100 - total flux cleaned = 0.00152083 Jy ! Component: 150 - total flux cleaned = 0.00207494 Jy ! Component: 200 - total flux cleaned = 0.00251953 Jy ! Total flux subtracted in 200 components = 0.00251953 Jy ! Clean residual min=-0.005796 max=0.005056 Jy/beam ! Clean residual mean=0.000022 rms=0.001188 Jy/beam ! Combined flux in latest and established models = 1.02473 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 173 components and 1.02473 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.059663Jy sigma=0.939397 ! Fit after self-cal, rms=0.059652Jy sigma=0.939240 ! 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.059652Jy sigma=0.939240 ! Fit after self-cal, rms=0.059552Jy sigma=0.938491 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.64 mas, bmaj=9.157 mas, bpa=30.89 degrees ! Estimated noise=1.29477 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00232199 Jy ! Component: 100 - total flux cleaned = 0.00327249 Jy ! Component: 150 - total flux cleaned = 0.00388809 Jy ! Component: 200 - total flux cleaned = 0.00408359 Jy ! Total flux subtracted in 200 components = 0.00408359 Jy ! Clean residual min=-0.006150 max=0.006131 Jy/beam ! Clean residual mean=0.000041 rms=0.001600 Jy/beam ! Combined flux in latest and established models = 1.02881 Jy selfcal ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 214 components and 1.02881 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.059432Jy sigma=0.937135 ! Fit after self-cal, rms=0.059422Jy sigma=0.936978 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.000469512 Jy ! Component: 100 - total flux cleaned = 0.000693292 Jy ! Component: 150 - total flux cleaned = 0.000773118 Jy ! Component: 200 - total flux cleaned = 0.000810579 Jy ! Total flux subtracted in 200 components = 0.000810579 Jy ! Clean residual min=-0.006181 max=0.006155 Jy/beam ! Clean residual mean=0.000042 rms=0.001585 Jy/beam ! Combined flux in latest and established models = 1.02963 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 229 components and 1.02963 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.059378Jy sigma=0.936428 ! Fit after self-cal, rms=0.059376Jy sigma=0.936369 ! 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=4.904 mas, bmaj=5.568 mas, bpa=32.53 degrees ! Estimated noise=1.07642 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 4.904 x 5.568 at 32.53 degrees (North through East) ! Clean map min=-0.0053363 max=0.72366 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.059376Jy sigma=0.936369 ! Fit after self-cal, rms=0.055685Jy sigma=0.867601 ! Inverting map and beam ! Estimated beam: bmin=4.902 mas, bmaj=5.567 mas, bpa=32.58 degrees ! Estimated noise=1.07633 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00125525 Jy ! Component: 100 - total flux cleaned = -0.0012534 Jy ! Component: 150 - total flux cleaned = -0.00131717 Jy ! Component: 200 - total flux cleaned = -0.00119538 Jy ! Total flux subtracted in 200 components = -0.00119538 Jy ! Clean residual min=-0.005145 max=0.004811 Jy/beam ! Clean residual mean=0.000017 rms=0.001129 Jy/beam ! Combined flux in latest and established models = 1.02843 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 241 components and 1.02843 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.055522Jy sigma=0.865043 ! Fit after self-cal, rms=0.055519Jy sigma=0.864943 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000133102 Jy ! Component: 100 - total flux cleaned = -2.00339e-05 Jy ! Component: 150 - total flux cleaned = 3.31772e-05 Jy ! Component: 200 - total flux cleaned = 8.52501e-05 Jy ! Total flux subtracted in 200 components = 8.52501e-05 Jy ! Clean residual min=-0.004982 max=0.004755 Jy/beam ! Clean residual mean=0.000017 rms=0.001119 Jy/beam ! Combined flux in latest and established models = 1.02851 Jy ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 252 components and 1.02851 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.055416Jy sigma=0.863305 ! Fit after self-cal, rms=0.055413Jy sigma=0.863235 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.055413Jy sigma=0.863235 ! Fit after self-cal, rms=0.055383Jy sigma=0.863008 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.016 mas, bmaj=4.772 mas, bpa=35.06 degrees ! Estimated noise=1.28666 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.553942 Jy ! Component: 100 - total flux cleaned = 0.700082 Jy ! Component: 150 - total flux cleaned = 0.771893 Jy ! Component: 200 - total flux cleaned = 0.819431 Jy ! Total flux subtracted in 200 components = 0.819431 Jy ! Clean residual min=-0.013414 max=0.026943 Jy/beam ! Clean residual mean=-0.000008 rms=0.003781 Jy/beam ! Combined flux in latest and established models = 0.819432 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 28 model components to the UV plane model. ! The established model now contains 28 components and 0.819432 Jy ! Inverting map and beam ! Estimated beam: bmin=4.903 mas, bmaj=5.567 mas, bpa=32.7 degrees ! Estimated noise=1.07596 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.0468284 Jy ! Component: 100 - total flux cleaned = 0.0786889 Jy ! Component: 150 - total flux cleaned = 0.10133 Jy ! Component: 200 - total flux cleaned = 0.117817 Jy ! Total flux subtracted in 200 components = 0.117817 Jy ! Clean residual min=-0.008764 max=0.009507 Jy/beam ! Clean residual mean=-0.000024 rms=0.002148 Jy/beam ! Combined flux in latest and established models = 0.937249 Jy ! Adding 75 model components to the UV plane model. ! The established model now contains 90 components and 0.937249 Jy ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0127866 Jy ! Component: 100 - total flux cleaned = 0.0231925 Jy ! Component: 150 - total flux cleaned = 0.0318001 Jy ! Component: 200 - total flux cleaned = 0.0391954 Jy ! Total flux subtracted in 200 components = 0.0391954 Jy ! Clean residual min=-0.007157 max=0.006420 Jy/beam ! Clean residual mean=-0.000006 rms=0.001541 Jy/beam ! Combined flux in latest and established models = 0.976444 Jy ! Adding 45 model components to the UV plane model. ! The established model now contains 123 components and 0.976444 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=8.645 mas, bmaj=9.158 mas, bpa=31.33 degrees ! Estimated noise=1.29404 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00977261 Jy ! Component: 100 - total flux cleaned = 0.0173934 Jy ! Component: 150 - total flux cleaned = 0.0235175 Jy ! Component: 200 - total flux cleaned = 0.0285261 Jy ! Total flux subtracted in 200 components = 0.0285261 Jy ! Clean residual min=-0.007659 max=0.005974 Jy/beam ! Clean residual mean=0.000014 rms=0.001757 Jy/beam ! Combined flux in latest and established models = 1.00497 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.00414672 Jy ! Component: 100 - total flux cleaned = 0.0076137 Jy ! Component: 150 - total flux cleaned = 0.0105332 Jy ! Component: 200 - total flux cleaned = 0.0130019 Jy ! Total flux subtracted in 200 components = 0.0130019 Jy ! Clean residual min=-0.006779 max=0.005926 Jy/beam ! Clean residual mean=0.000024 rms=0.001619 Jy/beam ! Combined flux in latest and established models = 1.01797 Jy ! Adding 147 model components to the UV plane model. ! The established model now contains 247 components and 1.01797 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=4.903 mas, bmaj=5.567 mas, bpa=32.7 degrees ! Estimated noise=1.07596 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00163299 Jy ! Component: 100 - total flux cleaned = 0.00340028 Jy ! Component: 150 - total flux cleaned = 0.00462649 Jy ! Component: 200 - total flux cleaned = 0.00558642 Jy ! Total flux subtracted in 200 components = 0.00558642 Jy ! Clean residual min=-0.005251 max=0.005125 Jy/beam ! Clean residual mean=0.000015 rms=0.001152 Jy/beam ! Combined flux in latest and established models = 1.02356 Jy selfcal ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 272 components and 1.02356 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.055621Jy sigma=0.869369 ! Fit after self-cal, rms=0.055504Jy sigma=0.866520 wmodel J1955+5131_S_b_map.mod ! Writing 272 model components to file: J1955+5131_S_b_map.mod wobs J1955+5131_S_b_uvs.fits ! Writing UV FITS file: J1955+5131_S_b_uvs.fits wwins J1955+5131_S_b_map.win ! wwins: Wrote 5 windows to J1955+5131_S_b_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 = -4.10283e-05 Jy ! Component: 100 - total flux cleaned = -0.000789933 Jy ! Component: 150 - total flux cleaned = -0.00110149 Jy ! Component: 200 - total flux cleaned = -0.00139204 Jy ! Component: 250 - total flux cleaned = -0.00139403 Jy ! Component: 300 - total flux cleaned = -0.00139414 Jy ! Component: 350 - total flux cleaned = -0.00120337 Jy ! Component: 400 - total flux cleaned = -0.000956887 Jy ! Component: 450 - total flux cleaned = -0.000537079 Jy ! Component: 500 - total flux cleaned = -0.000304223 Jy ! Component: 550 - total flux cleaned = -1.92355e-05 Jy ! Component: 600 - total flux cleaned = -2.02203e-05 Jy ! Component: 650 - total flux cleaned = 0.000142866 Jy ! Component: 700 - total flux cleaned = 0.00035561 Jy ! Component: 750 - total flux cleaned = 0.000669694 Jy ! Component: 800 - total flux cleaned = 0.000669399 Jy ! Component: 850 - total flux cleaned = 0.000923903 Jy ! Component: 900 - total flux cleaned = 0.00117309 Jy ! Component: 950 - total flux cleaned = 0.00131966 Jy ! Component: 1000 - total flux cleaned = 0.00151244 Jy ! Component: 1050 - total flux cleaned = 0.00156081 Jy ! Component: 1100 - total flux cleaned = 0.00174822 Jy ! Component: 1150 - total flux cleaned = 0.00193215 Jy ! Component: 1200 - total flux cleaned = 0.00202309 Jy ! Component: 1250 - total flux cleaned = 0.00211258 Jy ! Component: 1300 - total flux cleaned = 0.002423 Jy ! Component: 1350 - total flux cleaned = 0.0024666 Jy ! Component: 1400 - total flux cleaned = 0.00259658 Jy ! Component: 1450 - total flux cleaned = 0.00272506 Jy ! Component: 1500 - total flux cleaned = 0.00293616 Jy ! Component: 1550 - total flux cleaned = 0.00301884 Jy ! Component: 1600 - total flux cleaned = 0.00330721 Jy ! Component: 1650 - total flux cleaned = 0.00314425 Jy ! Component: 1700 - total flux cleaned = 0.0035077 Jy ! Component: 1750 - total flux cleaned = 0.00358733 Jy ! Component: 1800 - total flux cleaned = 0.00350796 Jy ! Component: 1850 - total flux cleaned = 0.00370448 Jy ! Component: 1900 - total flux cleaned = 0.00393698 Jy ! Component: 1950 - total flux cleaned = 0.00401375 Jy ! Component: 2000 - total flux cleaned = 0.00424107 Jy ! Component: 2050 - total flux cleaned = 0.00424126 Jy ! Component: 2100 - total flux cleaned = 0.00431616 Jy ! Component: 2150 - total flux cleaned = 0.00439039 Jy ! Component: 2200 - total flux cleaned = 0.00460966 Jy ! Component: 2250 - total flux cleaned = 0.00457336 Jy ! Component: 2300 - total flux cleaned = 0.00482416 Jy ! Component: 2350 - total flux cleaned = 0.00493103 Jy ! Component: 2400 - total flux cleaned = 0.00486055 Jy ! Component: 2450 - total flux cleaned = 0.00500003 Jy ! Component: 2500 - total flux cleaned = 0.00517313 Jy ! Component: 2550 - total flux cleaned = 0.00524164 Jy ! Component: 2600 - total flux cleaned = 0.00541132 Jy ! Component: 2650 - total flux cleaned = 0.00547923 Jy ! Component: 2700 - total flux cleaned = 0.00561282 Jy ! Component: 2750 - total flux cleaned = 0.00577919 Jy ! Component: 2800 - total flux cleaned = 0.0059769 Jy ! Component: 2850 - total flux cleaned = 0.00600947 Jy ! Component: 2900 - total flux cleaned = 0.00617128 Jy ! Component: 2950 - total flux cleaned = 0.00626724 Jy ! Component: 3000 - total flux cleaned = 0.00636347 Jy ! Component: 3050 - total flux cleaned = 0.00648989 Jy ! Component: 3100 - total flux cleaned = 0.00645878 Jy ! Component: 3150 - total flux cleaned = 0.00658316 Jy ! Component: 3200 - total flux cleaned = 0.00676769 Jy ! Component: 3250 - total flux cleaned = 0.00688981 Jy ! Component: 3300 - total flux cleaned = 0.0067986 Jy ! Component: 3350 - total flux cleaned = 0.0070401 Jy ! Component: 3400 - total flux cleaned = 0.00712992 Jy ! Component: 3450 - total flux cleaned = 0.00721898 Jy ! Component: 3500 - total flux cleaned = 0.00718966 Jy ! Component: 3550 - total flux cleaned = 0.00730682 Jy ! Component: 3600 - total flux cleaned = 0.00745263 Jy ! Component: 3650 - total flux cleaned = 0.00748148 Jy ! Component: 3700 - total flux cleaned = 0.00753878 Jy ! Component: 3750 - total flux cleaned = 0.00753903 Jy ! Component: 3800 - total flux cleaned = 0.00753879 Jy ! Component: 3850 - total flux cleaned = 0.00762329 Jy ! Component: 3900 - total flux cleaned = 0.00770681 Jy ! Component: 3950 - total flux cleaned = 0.00767884 Jy ! Component: 4000 - total flux cleaned = 0.00762353 Jy ! Component: 4050 - total flux cleaned = 0.00770576 Jy ! Total flux subtracted in 4096 components = 0.00781421 Jy ! Clean residual min=-0.001754 max=0.001848 Jy/beam ! Clean residual mean=0.000011 rms=0.000517 Jy/beam ! Combined flux in latest and established models = 1.03137 Jy keep ! Adding 305 model components to the UV plane model. ! The established model now contains 577 components and 1.03137 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /vlbi/uf001/uvm//J1955+5131_S_b_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 4.903 x 5.567 at 32.7 degrees (North through East) ! Clean map min=-0.0029863 max=0.72517 Jy/beam ! Writing clean map to FITS file: /vlbi/uf001/uvm//J1955+5131_S_b_map.fits wmodel /vlbi/uf001/uvm//J1955+5131_S_b_map.mod ! Writing 577 model components to file: /vlbi/uf001/uvm//J1955+5131_S_b_map.mod wobs /vlbi/uf001/uvm//J1955+5131_S_b_uvs.fits ! Writing UV FITS file: /vlbi/uf001/uvm//J1955+5131_S_b_uvs.fits wwins /vlbi/uf001/uvm//J1955+5131_S_b_map.win ! wwins: Wrote 6 windows to /vlbi/uf001/uvm//J1955+5131_S_b_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.725172 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000512814 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.212148 0.212148 0.424296 0.848593 1.69719 3.39437 6.78874 13.5775 27.155 54.3099 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0015 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.725 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 1414.1 quit ! Quitting program ! Log file /vlbi/uf001/uvm//J1955+5131_S_b_dfm.log closed on Sun Jan 13 22:24:04 2019