! Started logfile: /scr/pima/rv125_uvs/J0136+4751_S_dfm.log on Sun Dec 23 06:55:46 2018 float field_size field_size = 1024 float field_cell field_cell = 0.5000 float taper_size taper_size = 20.0 float freq freq = 2.23275 obs /scr/pima/rv125_uvs/J0136+4751_S_uva.fits ! Reading UV FITS file: /scr/pima/rv125_uvs/J0136+4751_S_uva.fits ! AN table 1: 2 integrations on 3 of 3 possible baselines. ! AN table 2: 39 integrations on 21 of 21 possible baselines. ! AN table 3: 118 integrations on 45 of 45 possible baselines. ! Apparent sampling: 0.243847 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J0136+4751 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 2.23275e+09 8e+06 1 8e+06 ! 02 2 2.26275e+09 8e+06 1 8e+06 ! 03 3 2.35275e+09 8e+06 1 8e+06 ! 04 4 2.37275e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 69 lines of history. ! ! Reading 5984 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 4 channels. ! Selecting polarization: RR, channels: 1..4 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 print "field_size =", field_size ! field_size = 1024 print "field_cell =", field_cell ! field_cell = 0.5 print "taper_size =", taper_size ! taper_size = 20 ![@/opt64/bin/pima_mupet_01.dfm J0136+4751_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 255 telescope corrections were flagged in sub-array 2. ! A total of 818 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 2. ! A total of 255 telescope corrections were flagged in sub-array 2. ! A total of 818 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 3. ! A total of 6 telescope corrections were flagged in sub-array 1. ! A total of 273 telescope corrections were flagged in sub-array 2. ! A total of 1180 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 4. ! A total of 255 telescope corrections were flagged in sub-array 2. ! A total of 818 telescope corrections were flagged in sub-array 3. ! ! Fit before self-cal, rms=1.495441Jy sigma=18.835080 ! Fit after self-cal, rms=0.451285Jy sigma=5.642425 ! 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.848 mas, bmaj=3.071 mas, bpa=-20.08 degrees ! Estimated noise=1.25443 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.705735 Jy ! Component: 100 - total flux cleaned = 0.972221 Jy ! Total flux subtracted in 100 components = 0.972221 Jy ! Clean residual min=-0.029868 max=0.120435 Jy/beam ! Clean residual mean=0.000082 rms=0.008100 Jy/beam ! Combined flux in latest and established models = 0.972222 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 13 components and 0.972222 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.383571Jy sigma=4.580028 ! Fit after self-cal, rms=0.383035Jy sigma=4.573624 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.141565 Jy ! Component: 100 - total flux cleaned = 0.229755 Jy ! Total flux subtracted in 100 components = 0.229755 Jy ! Clean residual min=-0.023149 max=0.048470 Jy/beam ! Clean residual mean=0.000045 rms=0.005999 Jy/beam ! Combined flux in latest and established models = 1.20198 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 31 components and 1.20198 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.281452Jy sigma=3.385900 ! Fit after self-cal, rms=0.280964Jy sigma=3.379248 ! Inverting map ! Added new window around map position (-1.5, 3). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0603946 Jy ! Component: 100 - total flux cleaned = 0.101113 Jy ! Total flux subtracted in 100 components = 0.101113 Jy ! Clean residual min=-0.022172 max=0.033328 Jy/beam ! Clean residual mean=0.000028 rms=0.005520 Jy/beam ! Combined flux in latest and established models = 1.30309 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 49 components and 1.30309 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.251033Jy sigma=3.063925 ! Fit after self-cal, rms=0.250855Jy sigma=3.061348 ! 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.089 mas, bmaj=3.683 mas, bpa=-30.8 degrees ! Estimated noise=1.04162 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.0623872 Jy ! Component: 100 - total flux cleaned = 0.101724 Jy ! Component: 150 - total flux cleaned = 0.127383 Jy ! Component: 200 - total flux cleaned = 0.144201 Jy ! Total flux subtracted in 200 components = 0.144201 Jy ! Clean residual min=-0.021088 max=0.038394 Jy/beam ! Clean residual mean=0.000167 rms=0.006460 Jy/beam ! Combined flux in latest and established models = 1.44729 Jy ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 71 components and 1.44729 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.228089Jy sigma=2.870838 ! Fit after self-cal, rms=0.227476Jy sigma=2.863400 ! Inverting map ! Added new window around map position (-7, 8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0397531 Jy ! Component: 100 - total flux cleaned = 0.0622316 Jy ! Component: 150 - total flux cleaned = 0.076675 Jy ! Component: 200 - total flux cleaned = 0.0834715 Jy ! Total flux subtracted in 200 components = 0.0834715 Jy ! Clean residual min=-0.022268 max=0.029628 Jy/beam ! Clean residual mean=0.000093 rms=0.005881 Jy/beam ! Combined flux in latest and established models = 1.53076 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 93 components and 1.53076 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.219825Jy sigma=2.764629 ! Fit after self-cal, rms=0.218629Jy sigma=2.751141 ! 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.089 x 3.683 at -30.8 degrees (North through East) ! Clean map min=-0.020542 max=0.92299 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 6 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* MC 0.00* MK 0.00* ! ! A total of 6 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! HN 0.88 MC 1.14 MK 1.08 NY 1.41 ! OV 0.89 SC 0.98 WN 1.00* ! ! A total of 60 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.15 HN 0.92 KP 0.94 LA 0.92 ! MK 1.02 NL 1.13 NY 1.35 OV 0.96 ! PT 1.04 SC 0.94 ! ! ! Correcting IF 2. ! A total of 6 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* MC 0.00* MK 0.00* ! ! A total of 6 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! HN 0.89 MC 0.95 MK 1.11 NY 1.25 ! OV 0.82 SC 0.90 WN 1.00* ! ! A total of 60 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.14 HN 0.91 KP 0.90 LA 0.91 ! MK 1.01 NL 1.14 NY 1.20 OV 1.16 ! PT 1.05 SC 0.87 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* MC 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! HN 0.00* MC 0.00* MK 0.00* NY 0.00* ! OV 0.00* SC 0.00* WN 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* HN 0.00* KP 0.00* LA 0.00* ! MK 0.00* NL 0.00* NY 0.00* OV 0.00* ! PT 0.00* SC 0.00* ! ! ! Correcting IF 4. ! A total of 6 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* MC 0.00* MK 0.00* ! ! A total of 6 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! HN 0.82 MC 1.13 MK 1.08 NY 1.19 ! OV 0.86 SC 0.86 WN 1.00* ! ! A total of 60 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.08 HN 0.89 KP 0.89 LA 0.89 ! MK 0.96 NL 1.10 NY 1.19 OV 1.00 ! PT 1.01 SC 0.87 ! ! ! Fit before self-cal, rms=0.218629Jy sigma=2.751141 ! Fit after self-cal, rms=0.179147Jy sigma=2.092045 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.276 mas, bmaj=3.75 mas, bpa=-35.54 degrees ! Estimated noise=1.09855 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.00223487 Jy ! Component: 100 - total flux cleaned = 0.00180988 Jy ! Component: 150 - total flux cleaned = -2.06585e-06 Jy ! Component: 200 - total flux cleaned = -0.00124562 Jy ! Total flux subtracted in 200 components = -0.00124562 Jy ! Clean residual min=-0.017958 max=0.014694 Jy/beam ! Clean residual mean=0.000037 rms=0.004042 Jy/beam ! Combined flux in latest and established models = 1.52952 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 109 components and 1.52952 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.177492Jy sigma=2.067791 ! Fit after self-cal, rms=0.177425Jy sigma=2.065830 ! Inverting map selfcal true, true, soltime1 ! Performing amp+phase self-cal over 120 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.177425Jy sigma=2.065830 ! Fit after self-cal, rms=0.132361Jy sigma=1.388233 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.278 mas, bmaj=3.739 mas, bpa=-34.91 degrees ! Estimated noise=1.10488 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00311902 Jy ! Component: 100 - total flux cleaned = 0.00314052 Jy ! Component: 150 - total flux cleaned = 0.00250081 Jy ! Component: 200 - total flux cleaned = 0.00163849 Jy ! Total flux subtracted in 200 components = 0.00163849 Jy ! Clean residual min=-0.008247 max=0.009101 Jy/beam ! Clean residual mean=0.000093 rms=0.002237 Jy/beam ! Combined flux in latest and established models = 1.53115 Jy selfcal ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 127 components and 1.53115 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.130584Jy sigma=1.362433 ! Fit after self-cal, rms=0.130510Jy sigma=1.360831 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.000517711 Jy ! Component: 100 - total flux cleaned = -0.00128803 Jy ! Component: 150 - total flux cleaned = -0.00198247 Jy ! Component: 200 - total flux cleaned = -0.00260931 Jy ! Total flux subtracted in 200 components = -0.00260931 Jy ! Clean residual min=-0.008078 max=0.008993 Jy/beam ! Clean residual mean=0.000096 rms=0.002188 Jy/beam ! Combined flux in latest and established models = 1.52855 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 137 components and 1.52854 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.129715Jy sigma=1.348352 ! Fit after self-cal, rms=0.129651Jy sigma=1.347859 ! Inverting map selfcal true, true, soltime2 ! Performing amp+phase self-cal over 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.129651Jy sigma=1.347859 ! Fit after self-cal, rms=0.096164Jy sigma=0.951008 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.282 mas, bmaj=3.746 mas, bpa=-34.24 degrees ! Estimated noise=1.10419 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000495511 Jy ! Component: 100 - total flux cleaned = -0.00112647 Jy ! Component: 150 - total flux cleaned = -0.00181247 Jy ! Component: 200 - total flux cleaned = -0.00195782 Jy ! Total flux subtracted in 200 components = -0.00195782 Jy ! Clean residual min=-0.004710 max=0.006729 Jy/beam ! Clean residual mean=0.000035 rms=0.001154 Jy/beam ! Combined flux in latest and established models = 1.52659 Jy selfcal ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 151 components and 1.52659 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.095543Jy sigma=0.938518 ! Fit after self-cal, rms=0.095446Jy sigma=0.937481 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 (0.5, 14.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00698942 Jy ! Component: 100 - total flux cleaned = 0.00838298 Jy ! Component: 150 - total flux cleaned = 0.00892333 Jy ! Component: 200 - total flux cleaned = 0.00928039 Jy ! Total flux subtracted in 200 components = 0.00928039 Jy ! Clean residual min=-0.004240 max=0.005148 Jy/beam ! Clean residual mean=0.000026 rms=0.001067 Jy/beam ! Combined flux in latest and established models = 1.53587 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 166 components and 1.53587 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.094819Jy sigma=0.927536 ! Fit after self-cal, rms=0.094756Jy sigma=0.926392 ! Inverting map ! Added new window around map position (0, 7.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00412999 Jy ! Component: 100 - total flux cleaned = 0.00450829 Jy ! Component: 150 - total flux cleaned = 0.00461138 Jy ! Component: 200 - total flux cleaned = 0.00461087 Jy ! Total flux subtracted in 200 components = 0.00461087 Jy ! Clean residual min=-0.004276 max=0.005151 Jy/beam ! Clean residual mean=0.000021 rms=0.001030 Jy/beam ! Combined flux in latest and established models = 1.54048 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 179 components and 1.54048 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.094435Jy sigma=0.920976 ! Fit after self-cal, rms=0.094361Jy sigma=0.920376 ! Inverting map ! Added new window around map position (-4.5, 28.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00450236 Jy ! Component: 100 - total flux cleaned = 0.00624607 Jy ! Component: 150 - total flux cleaned = 0.00708529 Jy ! Component: 200 - total flux cleaned = 0.00745298 Jy ! Total flux subtracted in 200 components = 0.00745298 Jy ! Clean residual min=-0.003966 max=0.004192 Jy/beam ! Clean residual mean=0.000014 rms=0.000988 Jy/beam ! Combined flux in latest and established models = 1.54793 Jy ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 198 components and 1.54793 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.093989Jy sigma=0.915179 ! Fit after self-cal, rms=0.093959Jy sigma=0.914799 ! Inverting map ! Added new window around map position (21, 24.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00400802 Jy ! Component: 100 - total flux cleaned = 0.00538663 Jy ! Component: 150 - total flux cleaned = 0.00595211 Jy ! Component: 200 - total flux cleaned = 0.00639091 Jy ! Total flux subtracted in 200 components = 0.00639091 Jy ! Clean residual min=-0.003938 max=0.003971 Jy/beam ! Clean residual mean=0.000008 rms=0.000962 Jy/beam ! Combined flux in latest and established models = 1.55432 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 212 components and 1.55432 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.093668Jy sigma=0.911188 ! Fit after self-cal, rms=0.093641Jy sigma=0.911004 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime2 ! Performing amp+phase self-cal over 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.093641Jy sigma=0.911004 ! Fit after self-cal, rms=0.093363Jy sigma=0.905344 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.129 mas, bmaj=9.411 mas, bpa=-79.51 degrees ! Estimated noise=1.60906 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00435138 Jy ! Component: 100 - total flux cleaned = 0.004867 Jy ! Component: 150 - total flux cleaned = 0.0054966 Jy ! Component: 200 - total flux cleaned = 0.00609229 Jy ! Total flux subtracted in 200 components = 0.00609229 Jy ! Clean residual min=-0.003937 max=0.003922 Jy/beam ! Clean residual mean=0.000024 rms=0.001094 Jy/beam ! Combined flux in latest and established models = 1.56041 Jy selfcal ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 231 components and 1.56041 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.093142Jy sigma=0.902425 ! Fit after self-cal, rms=0.093112Jy sigma=0.901833 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.00132923 Jy ! Component: 100 - total flux cleaned = 0.00171559 Jy ! Component: 150 - total flux cleaned = 0.00203998 Jy ! Component: 200 - total flux cleaned = 0.00237188 Jy ! Total flux subtracted in 200 components = 0.00237188 Jy ! Clean residual min=-0.003897 max=0.003813 Jy/beam ! Clean residual mean=0.000018 rms=0.001061 Jy/beam ! Combined flux in latest and established models = 1.56279 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 241 components and 1.56279 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.093049Jy sigma=0.900906 ! Fit after self-cal, rms=0.093044Jy sigma=0.900732 ! 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.285 mas, bmaj=3.752 mas, bpa=-34.44 degrees ! Estimated noise=1.10469 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.285 x 3.752 at -34.44 degrees (North through East) ! Clean map min=-0.003504 max=0.94883 Jy/beam if(peak(flux) > thresh) selfcal true, true clean clean_niter,clean_gain selfcal flux_cutoff = imstat(rms) * dynam repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) selfcal true, true end if ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.093044Jy sigma=0.900732 ! Fit after self-cal, rms=0.078613Jy sigma=0.777897 ! Inverting map and beam ! Estimated beam: bmin=3.287 mas, bmaj=3.751 mas, bpa=-34.79 degrees ! Estimated noise=1.1037 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000917088 Jy ! Component: 100 - total flux cleaned = -0.000911342 Jy ! Component: 150 - total flux cleaned = -0.00100892 Jy ! Component: 200 - total flux cleaned = -0.00101111 Jy ! Total flux subtracted in 200 components = -0.00101111 Jy ! Clean residual min=-0.003556 max=0.003733 Jy/beam ! Clean residual mean=0.000008 rms=0.000873 Jy/beam ! Combined flux in latest and established models = 1.56178 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 251 components and 1.56178 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.078387Jy sigma=0.773428 ! Fit after self-cal, rms=0.078370Jy sigma=0.773273 ! Inverting map ! Added new window around map position (-36, 125.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00294416 Jy ! Component: 100 - total flux cleaned = 0.00320372 Jy ! Component: 150 - total flux cleaned = 0.00345434 Jy ! Component: 200 - total flux cleaned = 0.00361295 Jy ! Total flux subtracted in 200 components = 0.00361295 Jy ! Clean residual min=-0.003607 max=0.003505 Jy/beam ! Clean residual mean=0.000008 rms=0.000858 Jy/beam ! Combined flux in latest and established models = 1.56539 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 261 components and 1.56539 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.078165Jy sigma=0.769767 ! Fit after self-cal, rms=0.078148Jy sigma=0.769618 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.078148Jy sigma=0.769618 ! Fit after self-cal, rms=0.078211Jy sigma=0.768926 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.916 mas, bmaj=3.165 mas, bpa=-32.23 degrees ! Estimated noise=1.3178 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.725409 Jy ! Component: 100 - total flux cleaned = 0.998087 Jy ! Component: 150 - total flux cleaned = 1.13798 Jy ! Component: 200 - total flux cleaned = 1.22456 Jy ! Total flux subtracted in 200 components = 1.22456 Jy ! Clean residual min=-0.011151 max=0.047039 Jy/beam ! Clean residual mean=0.000070 rms=0.003405 Jy/beam ! Combined flux in latest and established models = 1.22456 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 1.22456 Jy ! Inverting map and beam ! Estimated beam: bmin=3.289 mas, bmaj=3.755 mas, bpa=-34.91 degrees ! Estimated noise=1.10391 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.0969049 Jy ! Component: 100 - total flux cleaned = 0.157404 Jy ! Component: 150 - total flux cleaned = 0.197496 Jy ! Component: 200 - total flux cleaned = 0.225469 Jy ! Total flux subtracted in 200 components = 0.225469 Jy ! Clean residual min=-0.006745 max=0.015710 Jy/beam ! Clean residual mean=0.000132 rms=0.002031 Jy/beam ! Combined flux in latest and established models = 1.45003 Jy ! Adding 87 model components to the UV plane model. ! The established model now contains 103 components and 1.45003 Jy ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0204016 Jy ! Component: 100 - total flux cleaned = 0.035938 Jy ! Component: 150 - total flux cleaned = 0.0481408 Jy ! Component: 200 - total flux cleaned = 0.057905 Jy ! Total flux subtracted in 200 components = 0.057905 Jy ! Clean residual min=-0.004756 max=0.006784 Jy/beam ! Clean residual mean=0.000064 rms=0.001181 Jy/beam ! Combined flux in latest and established models = 1.50793 Jy ! Adding 62 model components to the UV plane model. ! The established model now contains 146 components and 1.50793 Jy ! Inverting map ! Added new window around map position (1, -3.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0082347 Jy ! Component: 100 - total flux cleaned = 0.015136 Jy ! Component: 150 - total flux cleaned = 0.0210539 Jy ! Component: 200 - total flux cleaned = 0.0262222 Jy ! Total flux subtracted in 200 components = 0.0262222 Jy ! Clean residual min=-0.004076 max=0.004435 Jy/beam ! Clean residual mean=0.000038 rms=0.000966 Jy/beam ! Combined flux in latest and established models = 1.53415 Jy ! Adding 52 model components to the UV plane model. ! The established model now contains 178 components and 1.53415 Jy ! Inverting map ! Added new window around map position (7.5, 5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00430811 Jy ! Component: 100 - total flux cleaned = 0.00733242 Jy ! Component: 150 - total flux cleaned = 0.00995572 Jy ! Component: 200 - total flux cleaned = 0.0123749 Jy ! Total flux subtracted in 200 components = 0.0123749 Jy ! Clean residual min=-0.003822 max=0.003923 Jy/beam ! Clean residual mean=0.000026 rms=0.000894 Jy/beam ! Combined flux in latest and established models = 1.54653 Jy ! Adding 53 model components to the UV plane model. ! The established model now contains 195 components and 1.54653 Jy ! Inverting map ! Added new window around map position (4, 23.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00332076 Jy ! Component: 100 - total flux cleaned = 0.00507156 Jy ! Component: 150 - total flux cleaned = 0.00670881 Jy ! Component: 200 - total flux cleaned = 0.00835331 Jy ! Total flux subtracted in 200 components = 0.00835331 Jy ! Clean residual min=-0.003641 max=0.003329 Jy/beam ! Clean residual mean=0.000019 rms=0.000856 Jy/beam ! Combined flux in latest and established models = 1.55488 Jy ! Adding 52 model components to the UV plane model. ! The established model now contains 215 components and 1.55488 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.123 mas, bmaj=9.418 mas, bpa=-79.58 degrees ! Estimated noise=1.60563 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00403006 Jy ! Component: 100 - total flux cleaned = 0.00672443 Jy ! Component: 150 - total flux cleaned = 0.00882806 Jy ! Component: 200 - total flux cleaned = 0.0105145 Jy ! Total flux subtracted in 200 components = 0.0105145 Jy ! Clean residual min=-0.003674 max=0.003632 Jy/beam ! Clean residual mean=0.000023 rms=0.001003 Jy/beam ! Combined flux in latest and established models = 1.5654 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.00126108 Jy ! Component: 100 - total flux cleaned = 0.00191411 Jy ! Component: 150 - total flux cleaned = 0.00236938 Jy ! Component: 200 - total flux cleaned = 0.00274238 Jy ! Total flux subtracted in 200 components = 0.00274238 Jy ! Clean residual min=-0.003529 max=0.003620 Jy/beam ! Clean residual mean=0.000017 rms=0.000996 Jy/beam ! Combined flux in latest and established models = 1.56814 Jy ! Adding 86 model components to the UV plane model. ! The established model now contains 287 components and 1.56814 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.289 mas, bmaj=3.755 mas, bpa=-34.91 degrees ! Estimated noise=1.10391 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00121501 Jy ! Component: 100 - total flux cleaned = -0.00121493 Jy ! Component: 150 - total flux cleaned = -0.000971622 Jy ! Component: 200 - total flux cleaned = -0.000742335 Jy ! Total flux subtracted in 200 components = -0.000742335 Jy ! Clean residual min=-0.003439 max=0.003335 Jy/beam ! Clean residual mean=0.000007 rms=0.000824 Jy/beam ! Combined flux in latest and established models = 1.5674 Jy selfcal ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 299 components and 1.5674 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.077923Jy sigma=0.759285 ! Fit after self-cal, rms=0.077866Jy sigma=0.756568 wmodel J0136+4751_S_map.mod ! Writing 299 model components to file: J0136+4751_S_map.mod wobs J0136+4751_S_uvs.fits ! Writing UV FITS file: J0136+4751_S_uvs.fits wwins J0136+4751_S_map.win ! wwins: Wrote 11 windows to J0136+4751_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 = -7.44032e-05 Jy ! Component: 100 - total flux cleaned = 0.000220467 Jy ! Component: 150 - total flux cleaned = 0.000277995 Jy ! Component: 200 - total flux cleaned = 0.000164646 Jy ! Component: 250 - total flux cleaned = 0.000109799 Jy ! Component: 300 - total flux cleaned = 5.57386e-05 Jy ! Component: 350 - total flux cleaned = -0.00010357 Jy ! Component: 400 - total flux cleaned = 5.56132e-05 Jy ! Component: 450 - total flux cleaned = 5.49763e-05 Jy ! Component: 500 - total flux cleaned = -0.000100609 Jy ! Component: 550 - total flux cleaned = -0.000253718 Jy ! Component: 600 - total flux cleaned = -0.000154403 Jy ! Component: 650 - total flux cleaned = -0.000404074 Jy ! Component: 700 - total flux cleaned = -0.000502372 Jy ! Component: 750 - total flux cleaned = -0.00045373 Jy ! Component: 800 - total flux cleaned = -0.000599304 Jy ! Component: 850 - total flux cleaned = -0.000647243 Jy ! Component: 900 - total flux cleaned = -0.000742859 Jy ! Component: 950 - total flux cleaned = -0.000930894 Jy ! Component: 1000 - total flux cleaned = -0.00102486 Jy ! Component: 1050 - total flux cleaned = -0.00107125 Jy ! Component: 1100 - total flux cleaned = -0.0011635 Jy ! Component: 1150 - total flux cleaned = -0.0013467 Jy ! Component: 1200 - total flux cleaned = -0.00143738 Jy ! Component: 1250 - total flux cleaned = -0.00152794 Jy ! Component: 1300 - total flux cleaned = -0.00161755 Jy ! Component: 1350 - total flux cleaned = -0.00188444 Jy ! Component: 1400 - total flux cleaned = -0.00179579 Jy ! Component: 1450 - total flux cleaned = -0.00183986 Jy ! Component: 1500 - total flux cleaned = -0.00201482 Jy ! Component: 1550 - total flux cleaned = -0.0019276 Jy ! Component: 1600 - total flux cleaned = -0.00192762 Jy ! Component: 1650 - total flux cleaned = -0.00214258 Jy ! Component: 1700 - total flux cleaned = -0.0021854 Jy ! Component: 1750 - total flux cleaned = -0.00214238 Jy ! Component: 1800 - total flux cleaned = -0.00231141 Jy ! Component: 1850 - total flux cleaned = -0.00226951 Jy ! Component: 1900 - total flux cleaned = -0.00226979 Jy ! Component: 1950 - total flux cleaned = -0.00226969 Jy ! Component: 2000 - total flux cleaned = -0.0024355 Jy ! Component: 2050 - total flux cleaned = -0.00243593 Jy ! Component: 2100 - total flux cleaned = -0.00243606 Jy ! Component: 2150 - total flux cleaned = -0.00239503 Jy ! Component: 2200 - total flux cleaned = -0.0026393 Jy ! Component: 2250 - total flux cleaned = -0.00272 Jy ! Component: 2300 - total flux cleaned = -0.00263928 Jy ! Component: 2350 - total flux cleaned = -0.00263927 Jy ! Component: 2400 - total flux cleaned = -0.00279937 Jy ! Component: 2450 - total flux cleaned = -0.00275943 Jy ! Component: 2500 - total flux cleaned = -0.00275945 Jy ! Component: 2550 - total flux cleaned = -0.00275969 Jy ! Component: 2600 - total flux cleaned = -0.00283835 Jy ! Component: 2650 - total flux cleaned = -0.00283829 Jy ! Component: 2700 - total flux cleaned = -0.00283831 Jy ! Component: 2750 - total flux cleaned = -0.00295523 Jy ! Component: 2800 - total flux cleaned = -0.00295509 Jy ! Component: 2850 - total flux cleaned = -0.00299363 Jy ! Component: 2900 - total flux cleaned = -0.00303209 Jy ! Component: 2950 - total flux cleaned = -0.00295532 Jy ! Component: 3000 - total flux cleaned = -0.00314649 Jy ! Component: 3050 - total flux cleaned = -0.00314655 Jy ! Component: 3100 - total flux cleaned = -0.00310872 Jy ! Component: 3150 - total flux cleaned = -0.0034109 Jy ! Component: 3200 - total flux cleaned = -0.0032603 Jy ! Component: 3250 - total flux cleaned = -0.00341035 Jy ! Component: 3300 - total flux cleaned = -0.00337285 Jy ! Component: 3350 - total flux cleaned = -0.0033729 Jy ! Component: 3400 - total flux cleaned = -0.00337271 Jy ! Component: 3450 - total flux cleaned = -0.00340988 Jy ! Component: 3500 - total flux cleaned = -0.00352053 Jy ! Component: 3550 - total flux cleaned = -0.00341023 Jy ! Component: 3600 - total flux cleaned = -0.00352005 Jy ! Component: 3650 - total flux cleaned = -0.00373897 Jy ! Component: 3700 - total flux cleaned = -0.00359335 Jy ! Component: 3750 - total flux cleaned = -0.00362963 Jy ! Component: 3800 - total flux cleaned = -0.00355751 Jy ! Component: 3850 - total flux cleaned = -0.00359346 Jy ! Component: 3900 - total flux cleaned = -0.00370102 Jy ! Component: 3950 - total flux cleaned = -0.00359345 Jy ! Component: 4000 - total flux cleaned = -0.00377171 Jy ! Component: 4050 - total flux cleaned = -0.0037004 Jy ! Total flux subtracted in 4096 components = -0.00362922 Jy ! Clean residual min=-0.002218 max=0.002021 Jy/beam ! Clean residual mean=0.000001 rms=0.000638 Jy/beam ! Combined flux in latest and established models = 1.56377 Jy keep ! Adding 338 model components to the UV plane model. ! The established model now contains 637 components and 1.56377 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /scr/pima/rv125_uvs/J0136+4751_S_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.289 x 3.755 at -34.91 degrees (North through East) ! Clean map min=-0.0028006 max=0.95193 Jy/beam ! Writing clean map to FITS file: /scr/pima/rv125_uvs/J0136+4751_S_map.fits wmodel /scr/pima/rv125_uvs/J0136+4751_S_map.mod ! Writing 637 model components to file: /scr/pima/rv125_uvs/J0136+4751_S_map.mod wobs /scr/pima/rv125_uvs/J0136+4751_S_uvs.fits ! Writing UV FITS file: /scr/pima/rv125_uvs/J0136+4751_S_uvs.fits wwins /scr/pima/rv125_uvs/J0136+4751_S_map.win ! wwins: Wrote 12 windows to /scr/pima/rv125_uvs/J0136+4751_S_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.951929 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000636722 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.200663 0.200663 0.401325 0.80265 1.6053 3.2106 6.4212 12.8424 25.6848 51.3696 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0019 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.951 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 1495.05 quit ! Quitting program ! Log file /scr/pima/rv125_uvs/J0136+4751_S_dfm.log closed on Sun Dec 23 06:55:52 2018