! Started logfile: /scr/pima/rdv91_uvs/J1753+2848_S_dfm.log on Fri Jun 15 11:28:57 2018 float field_size field_size = 1024 float field_cell field_cell = 0.5000 float taper_size taper_size = 20.0 float freq freq = 2.23299 obs /scr/pima/rdv91_uvs/J1753+2848_S_uva.fits ! Reading UV FITS file: /scr/pima/rdv91_uvs/J1753+2848_S_uva.fits ! AN table 1: 329 integrations on 78 of 78 possible baselines. ! Apparent sampling: 0.30477 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1753+2848 ! ! 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.23299e+09 8e+06 1 8e+06 ! 02 2 2.26299e+09 8e+06 1 8e+06 ! 03 3 2.35299e+09 8e+06 1 8e+06 ! 04 4 2.37299e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 74 lines of history. ! ! Reading 31284 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 J1753+2848_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 2736 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 2736 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 2736 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 2736 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=1.280184Jy sigma=17.526184 ! Fit after self-cal, rms=0.291587Jy sigma=3.630428 ! 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.377 mas, bmaj=3.847 mas, bpa=-21.05 degrees ! Estimated noise=0.655503 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.591382 Jy ! Component: 100 - total flux cleaned = 0.720384 Jy ! Total flux subtracted in 100 components = 0.720384 Jy ! Clean residual min=-0.010916 max=0.037171 Jy/beam ! Clean residual mean=0.000105 rms=0.003134 Jy/beam ! Combined flux in latest and established models = 0.720384 Jy ! Performing phase self-cal ! Adding 3 model components to the UV plane model. ! The established model now contains 3 components and 0.720384 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.169479Jy sigma=2.023767 ! Fit after self-cal, rms=0.169476Jy sigma=2.023743 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0378167 Jy ! Component: 100 - total flux cleaned = 0.0579098 Jy ! Total flux subtracted in 100 components = 0.0579098 Jy ! Clean residual min=-0.008267 max=0.010647 Jy/beam ! Clean residual mean=0.000003 rms=0.002184 Jy/beam ! Combined flux in latest and established models = 0.778294 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 14 components and 0.778294 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.161295Jy sigma=1.840185 ! Fit after self-cal, rms=0.161308Jy sigma=1.840109 ! Inverting map print "============== Finished uniform weighting clean ===============" ! ============== Finished uniform weighting clean =============== ! uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.716 mas, bmaj=4.624 mas, bpa=-19.06 degrees ! Estimated noise=0.358772 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.0134138 Jy ! Component: 100 - total flux cleaned = 0.0197527 Jy ! Component: 150 - total flux cleaned = 0.0227447 Jy ! Component: 200 - total flux cleaned = 0.0230226 Jy ! Total flux subtracted in 200 components = 0.0230226 Jy ! Clean residual min=-0.007568 max=0.008566 Jy/beam ! Clean residual mean=-0.000056 rms=0.001802 Jy/beam ! Combined flux in latest and established models = 0.801317 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 29 components and 0.801317 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.160517Jy sigma=1.818595 ! Fit after self-cal, rms=0.160517Jy sigma=1.818528 ! Inverting map print "============== Finished natural weighting clean ==============" ! ============== Finished natural weighting clean ============== ! restore ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.716 x 4.624 at -19.06 degrees (North through East) ! Clean map min=-0.0071271 max=0.7637 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 21 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00 FD 1.00 FT 1.00* HN 0.97 ! KP 0.98 LA 1.00 MK 0.93 NL 1.01 ! NY 1.00 OV 0.98 PT 0.99 SC 0.98 ! WZ 1.00 ! ! ! Correcting IF 2. ! A total of 21 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 FD 1.00 FT 1.00* HN 0.99 ! KP 0.98 LA 1.01 MK 0.95 NL 1.00 ! NY 1.00 OV 0.99 PT 0.99 SC 0.98 ! WZ 1.01 ! ! ! Correcting IF 3. ! A total of 21 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00 FD 1.00 FT 1.00* HN 1.01 ! KP 1.00 LA 1.00 MK 0.96 NL 1.01 ! NY 1.01 OV 0.98 PT 1.01 SC 0.97 ! WZ 1.03 ! ! ! Correcting IF 4. ! A total of 21 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.11 FD 1.02 FT 1.00* HN 1.00 ! KP 0.99 LA 1.01 MK 0.97 NL 1.02 ! NY 1.01 OV 1.01 PT 1.01 SC 0.98 ! WZ 1.03 ! ! ! Fit before self-cal, rms=0.160517Jy sigma=1.818528 ! Fit after self-cal, rms=0.158576Jy sigma=1.780524 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.711 mas, bmaj=4.625 mas, bpa=-18.92 degrees ! Estimated noise=0.358263 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.00137304 Jy ! Component: 100 - total flux cleaned = 0.00138187 Jy ! Component: 150 - total flux cleaned = 0.0013075 Jy ! Component: 200 - total flux cleaned = 0.00117805 Jy ! Total flux subtracted in 200 components = 0.00117805 Jy ! Clean residual min=-0.005890 max=0.006272 Jy/beam ! Clean residual mean=-0.000059 rms=0.001414 Jy/beam ! Combined flux in latest and established models = 0.802495 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 38 components and 0.802495 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.158466Jy sigma=1.779307 ! Fit after self-cal, rms=0.158478Jy sigma=1.779256 ! 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.158478Jy sigma=1.779256 ! Fit after self-cal, rms=0.156369Jy sigma=1.756902 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.707 mas, bmaj=4.63 mas, bpa=-18.93 degrees ! Estimated noise=0.358061 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00116581 Jy ! Component: 100 - total flux cleaned = 0.00144446 Jy ! Component: 150 - total flux cleaned = 0.00138844 Jy ! Component: 200 - total flux cleaned = 0.0013407 Jy ! Total flux subtracted in 200 components = 0.0013407 Jy ! Clean residual min=-0.003827 max=0.004581 Jy/beam ! Clean residual mean=0.000097 rms=0.000960 Jy/beam ! Combined flux in latest and established models = 0.803835 Jy selfcal ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 42 components and 0.803835 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.156301Jy sigma=1.756265 ! Fit after self-cal, rms=0.156309Jy sigma=1.756223 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 (-20.5, -92). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00414245 Jy ! Component: 100 - total flux cleaned = 0.00503801 Jy ! Component: 150 - total flux cleaned = 0.00539701 Jy ! Component: 200 - total flux cleaned = 0.00569768 Jy ! Total flux subtracted in 200 components = 0.00569768 Jy ! Clean residual min=-0.003751 max=0.003912 Jy/beam ! Clean residual mean=0.000075 rms=0.000913 Jy/beam ! Combined flux in latest and established models = 0.809533 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 57 components and 0.809533 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.156174Jy sigma=1.754535 ! Fit after self-cal, rms=0.156162Jy sigma=1.754398 ! Inverting map ! Added new window around map position (-1, 19.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0019708 Jy ! Component: 100 - total flux cleaned = 0.00234225 Jy ! Component: 150 - total flux cleaned = 0.00270111 Jy ! Component: 200 - total flux cleaned = 0.00295872 Jy ! Total flux subtracted in 200 components = 0.00295872 Jy ! Clean residual min=-0.003499 max=0.003793 Jy/beam ! Clean residual mean=0.000059 rms=0.000888 Jy/beam ! Combined flux in latest and established models = 0.812492 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 71 components and 0.812492 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.156027Jy sigma=1.753147 ! Fit after self-cal, rms=0.156021Jy sigma=1.753107 ! Inverting map ! Added new window around map position (21, -53). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00117951 Jy ! Component: 100 - total flux cleaned = 0.000305162 Jy ! Component: 150 - total flux cleaned = -0.000262405 Jy ! Component: 200 - total flux cleaned = -0.000385607 Jy ! Total flux subtracted in 200 components = -0.000385607 Jy ! Clean residual min=-0.003575 max=0.003592 Jy/beam ! Clean residual mean=0.000062 rms=0.000865 Jy/beam ! Combined flux in latest and established models = 0.812106 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 91 components and 0.812106 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.155899Jy sigma=1.751920 ! Fit after self-cal, rms=0.155897Jy sigma=1.751887 ! 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.155897Jy sigma=1.751887 ! Fit after self-cal, rms=0.155273Jy sigma=1.746297 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.705 mas, bmaj=4.63 mas, bpa=-18.94 degrees ! Estimated noise=0.358143 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 8.8502e-05 Jy ! Component: 100 - total flux cleaned = 0.000146368 Jy ! Component: 150 - total flux cleaned = 9.19309e-05 Jy ! Component: 200 - total flux cleaned = 9.05426e-05 Jy ! Total flux subtracted in 200 components = 9.05426e-05 Jy ! Clean residual min=-0.002484 max=0.003852 Jy/beam ! Clean residual mean=0.000095 rms=0.000693 Jy/beam ! Combined flux in latest and established models = 0.812197 Jy selfcal ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 105 components and 0.812197 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.155230Jy sigma=1.745685 ! Fit after self-cal, rms=0.155247Jy sigma=1.745631 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 (5.5, -4.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00311171 Jy ! Component: 100 - total flux cleaned = 0.00377649 Jy ! Component: 150 - total flux cleaned = 0.00366013 Jy ! Component: 200 - total flux cleaned = 0.00366092 Jy ! Total flux subtracted in 200 components = 0.00366092 Jy ! Clean residual min=-0.002481 max=0.003209 Jy/beam ! Clean residual mean=0.000071 rms=0.000650 Jy/beam ! Combined flux in latest and established models = 0.815858 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 127 components and 0.815857 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.155162Jy sigma=1.744267 ! Fit after self-cal, rms=0.155155Jy sigma=1.744150 ! Inverting map ! Added new window around map position (-11, 71). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00134829 Jy ! Component: 100 - total flux cleaned = 0.001465 Jy ! Component: 150 - total flux cleaned = 0.00180162 Jy ! Component: 200 - total flux cleaned = 0.00201024 Jy ! Total flux subtracted in 200 components = 0.00201024 Jy ! Clean residual min=-0.002398 max=0.002724 Jy/beam ! Clean residual mean=0.000059 rms=0.000625 Jy/beam ! Combined flux in latest and established models = 0.817868 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 138 components and 0.817868 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.155085Jy sigma=1.743282 ! Fit after self-cal, rms=0.155080Jy sigma=1.743257 ! Inverting map ! Added new window around map position (63.5, 81). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00251787 Jy ! Component: 100 - total flux cleaned = 0.00346907 Jy ! Component: 150 - total flux cleaned = 0.00361173 Jy ! Component: 200 - total flux cleaned = 0.00378849 Jy ! Total flux subtracted in 200 components = 0.00378849 Jy ! Clean residual min=-0.002470 max=0.002745 Jy/beam ! Clean residual mean=0.000052 rms=0.000603 Jy/beam ! Combined flux in latest and established models = 0.821656 Jy ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 156 components and 0.821656 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.155027Jy sigma=1.742581 ! Fit after self-cal, rms=0.155024Jy sigma=1.742533 ! Inverting map ! Added new window around map position (50, -0.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00242786 Jy ! Component: 100 - total flux cleaned = 0.00345406 Jy ! Component: 150 - total flux cleaned = 0.00396885 Jy ! Component: 200 - total flux cleaned = 0.00431619 Jy ! Total flux subtracted in 200 components = 0.00431619 Jy ! Clean residual min=-0.002352 max=0.002548 Jy/beam ! Clean residual mean=0.000034 rms=0.000580 Jy/beam ! Combined flux in latest and established models = 0.825972 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 167 components and 0.825972 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.154969Jy sigma=1.741922 ! Fit after self-cal, rms=0.154965Jy sigma=1.741908 ! Inverting map ! Added new window around map position (105.5, 21). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00153728 Jy ! Component: 100 - total flux cleaned = 0.00177788 Jy ! Component: 150 - total flux cleaned = 0.00197898 Jy ! Component: 200 - total flux cleaned = 0.00234521 Jy ! Total flux subtracted in 200 components = 0.00234521 Jy ! Clean residual min=-0.002229 max=0.002469 Jy/beam ! Clean residual mean=0.000033 rms=0.000562 Jy/beam ! Combined flux in latest and established models = 0.828317 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 182 components and 0.828317 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.154918Jy sigma=1.741227 ! Fit after self-cal, rms=0.154917Jy sigma=1.741205 ! Inverting map ! Added new window around map position (-97, -11). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00200756 Jy ! Component: 100 - total flux cleaned = 0.00238152 Jy ! Component: 150 - total flux cleaned = 0.00266521 Jy ! Component: 200 - total flux cleaned = 0.00292936 Jy ! Total flux subtracted in 200 components = 0.00292936 Jy ! Clean residual min=-0.002145 max=0.002352 Jy/beam ! Clean residual mean=0.000028 rms=0.000544 Jy/beam ! Combined flux in latest and established models = 0.831247 Jy ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 202 components and 0.831247 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.154867Jy sigma=1.740585 ! Fit after self-cal, rms=0.154866Jy sigma=1.740554 ! Inverting map ! Added new window around map position (57, 39). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0021912 Jy ! Component: 100 - total flux cleaned = 0.00265399 Jy ! Component: 150 - total flux cleaned = 0.00289788 Jy ! Component: 200 - total flux cleaned = 0.00307734 Jy ! Total flux subtracted in 200 components = 0.00307734 Jy ! Clean residual min=-0.002123 max=0.002380 Jy/beam ! Clean residual mean=0.000020 rms=0.000522 Jy/beam ! Combined flux in latest and established models = 0.834324 Jy ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 222 components and 0.834324 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.154812Jy sigma=1.739935 ! Fit after self-cal, rms=0.154809Jy sigma=1.739919 ! Inverting map ! Added new window around map position (-116, -8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000866534 Jy ! Component: 100 - total flux cleaned = 0.00120282 Jy ! Component: 150 - total flux cleaned = 0.00172789 Jy ! Component: 200 - total flux cleaned = 0.00217217 Jy ! Total flux subtracted in 200 components = 0.00217217 Jy ! Clean residual min=-0.002094 max=0.002172 Jy/beam ! Clean residual mean=0.000019 rms=0.000507 Jy/beam ! Combined flux in latest and established models = 0.836496 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 237 components and 0.836496 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.154760Jy sigma=1.739317 ! Fit after self-cal, rms=0.154757Jy sigma=1.739295 ! 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.154757Jy sigma=1.739295 ! Fit after self-cal, rms=0.154602Jy sigma=1.738997 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.728 mas, bmaj=9.224 mas, bpa=-22.25 degrees ! Estimated noise=0.492247 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00134236 Jy ! Component: 100 - total flux cleaned = 0.00208182 Jy ! Component: 150 - total flux cleaned = 0.00278046 Jy ! Component: 200 - total flux cleaned = 0.00342056 Jy ! Total flux subtracted in 200 components = 0.00342056 Jy ! Clean residual min=-0.001870 max=0.002398 Jy/beam ! Clean residual mean=0.000045 rms=0.000606 Jy/beam ! Combined flux in latest and established models = 0.839917 Jy selfcal ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 265 components and 0.839917 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.154595Jy sigma=1.738737 ! Fit after self-cal, rms=0.154595Jy sigma=1.738722 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.000806346 Jy ! Component: 100 - total flux cleaned = 0.00140897 Jy ! Component: 150 - total flux cleaned = 0.00186912 Jy ! Component: 200 - total flux cleaned = 0.00223989 Jy ! Total flux subtracted in 200 components = 0.00223989 Jy ! Clean residual min=-0.001871 max=0.002308 Jy/beam ! Clean residual mean=0.000029 rms=0.000587 Jy/beam ! Combined flux in latest and established models = 0.842157 Jy ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 280 components and 0.842157 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.154593Jy sigma=1.738632 ! Fit after self-cal, rms=0.154591Jy sigma=1.738631 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "============== Finished amplitude self-cal ==============" ! ============== Finished amplitude self-cal ============== restore ! Inverting map and beam ! Estimated beam: bmin=2.704 mas, bmaj=4.631 mas, bpa=-18.95 degrees ! Estimated noise=0.358064 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.704 x 4.631 at -18.95 degrees (North through East) ! Clean map min=-0.0021592 max=0.76211 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.154591Jy sigma=1.738631 ! Fit after self-cal, rms=0.146290Jy sigma=1.575272 ! Inverting map and beam ! Estimated beam: bmin=2.703 mas, bmaj=4.63 mas, bpa=-18.98 degrees ! Estimated noise=0.358243 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -8.35475e-05 Jy ! Component: 100 - total flux cleaned = -3.69983e-05 Jy ! Component: 150 - total flux cleaned = 0.000142563 Jy ! Component: 200 - total flux cleaned = 0.000438028 Jy ! Total flux subtracted in 200 components = 0.000438028 Jy ! Clean residual min=-0.001946 max=0.002037 Jy/beam ! Clean residual mean=0.000008 rms=0.000460 Jy/beam ! Combined flux in latest and established models = 0.842595 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 296 components and 0.842595 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.146251Jy sigma=1.574912 ! Fit after self-cal, rms=0.146250Jy sigma=1.574898 ! Inverting map ! Added new window around map position (-18, -113). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000667534 Jy ! Component: 100 - total flux cleaned = 0.000843618 Jy ! Component: 150 - total flux cleaned = 0.00113825 Jy ! Component: 200 - total flux cleaned = 0.00153891 Jy ! Total flux subtracted in 200 components = 0.00153891 Jy ! Clean residual min=-0.001928 max=0.001756 Jy/beam ! Clean residual mean=0.000004 rms=0.000449 Jy/beam ! Combined flux in latest and established models = 0.844134 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 318 components and 0.844134 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.146205Jy sigma=1.574479 ! Fit after self-cal, rms=0.146202Jy sigma=1.574463 ! 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.146202Jy sigma=1.574463 ! Fit after self-cal, rms=0.146182Jy sigma=1.574408 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.361 mas, bmaj=3.873 mas, bpa=-21.26 degrees ! Estimated noise=0.648309 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.58799 Jy ! Component: 100 - total flux cleaned = 0.716524 Jy ! Component: 150 - total flux cleaned = 0.756764 Jy ! Component: 200 - total flux cleaned = 0.776699 Jy ! Total flux subtracted in 200 components = 0.776699 Jy ! Clean residual min=-0.005648 max=0.009788 Jy/beam ! Clean residual mean=0.000084 rms=0.001427 Jy/beam ! Combined flux in latest and established models = 0.776699 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 12 model components to the UV plane model. ! The established model now contains 12 components and 0.776699 Jy ! Inverting map and beam ! Estimated beam: bmin=2.702 mas, bmaj=4.631 mas, bpa=-18.98 degrees ! Estimated noise=0.358309 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.0148531 Jy ! Component: 100 - total flux cleaned = 0.0232783 Jy ! Component: 150 - total flux cleaned = 0.0287818 Jy ! Component: 200 - total flux cleaned = 0.0330693 Jy ! Total flux subtracted in 200 components = 0.0330693 Jy ! Clean residual min=-0.002948 max=0.003355 Jy/beam ! Clean residual mean=0.000120 rms=0.000760 Jy/beam ! Combined flux in latest and established models = 0.809768 Jy ! Adding 39 model components to the UV plane model. ! The established model now contains 47 components and 0.809768 Jy ! Inverting map ! Added new window around map position (13, 5.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00391747 Jy ! Component: 100 - total flux cleaned = 0.00708249 Jy ! Component: 150 - total flux cleaned = 0.00948656 Jy ! Component: 200 - total flux cleaned = 0.0113774 Jy ! Total flux subtracted in 200 components = 0.0113774 Jy ! Clean residual min=-0.002637 max=0.002723 Jy/beam ! Clean residual mean=0.000076 rms=0.000664 Jy/beam ! Combined flux in latest and established models = 0.821146 Jy ! Adding 45 model components to the UV plane model. ! The established model now contains 70 components and 0.821146 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.722 mas, bmaj=9.219 mas, bpa=-22.25 degrees ! Estimated noise=0.493087 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00401767 Jy ! Component: 100 - total flux cleaned = 0.00584556 Jy ! Component: 150 - total flux cleaned = 0.00705367 Jy ! Component: 200 - total flux cleaned = 0.00809997 Jy ! Total flux subtracted in 200 components = 0.00809997 Jy ! Clean residual min=-0.002191 max=0.002475 Jy/beam ! Clean residual mean=0.000112 rms=0.000706 Jy/beam ! Combined flux in latest and established models = 0.829246 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.00130406 Jy ! Component: 100 - total flux cleaned = 0.00250689 Jy ! Component: 150 - total flux cleaned = 0.0036251 Jy ! Component: 200 - total flux cleaned = 0.00432413 Jy ! Total flux subtracted in 200 components = 0.00432413 Jy ! Clean residual min=-0.001998 max=0.002259 Jy/beam ! Clean residual mean=0.000083 rms=0.000631 Jy/beam ! Combined flux in latest and established models = 0.83357 Jy ! Adding 68 model components to the UV plane model. ! The established model now contains 129 components and 0.83357 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.702 mas, bmaj=4.631 mas, bpa=-18.98 degrees ! Estimated noise=0.358309 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00201208 Jy ! Component: 100 - total flux cleaned = 0.00220297 Jy ! Component: 150 - total flux cleaned = 0.00265835 Jy ! Component: 200 - total flux cleaned = 0.00308802 Jy ! Total flux subtracted in 200 components = 0.00308802 Jy ! Clean residual min=-0.002201 max=0.002400 Jy/beam ! Clean residual mean=0.000018 rms=0.000546 Jy/beam ! Combined flux in latest and established models = 0.836658 Jy selfcal ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 154 components and 0.836658 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.146546Jy sigma=1.578693 ! Fit after self-cal, rms=0.146543Jy sigma=1.578483 wmodel J1753+2848_S_map.mod ! Writing 154 model components to file: J1753+2848_S_map.mod wobs J1753+2848_S_uvs.fits ! Writing UV FITS file: J1753+2848_S_uvs.fits wwins J1753+2848_S_map.win ! wwins: Wrote 14 windows to J1753+2848_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.000676265 Jy ! Component: 100 - total flux cleaned = 0.000965556 Jy ! Component: 150 - total flux cleaned = 0.00100589 Jy ! Component: 200 - total flux cleaned = 0.00100596 Jy ! Component: 250 - total flux cleaned = 0.000858445 Jy ! Component: 300 - total flux cleaned = 0.000678886 Jy ! Component: 350 - total flux cleaned = 0.00050312 Jy ! Component: 400 - total flux cleaned = 0.000503254 Jy ! Component: 450 - total flux cleaned = 0.000434755 Jy ! Component: 500 - total flux cleaned = 0.00036774 Jy ! Component: 550 - total flux cleaned = 0.000302266 Jy ! Component: 600 - total flux cleaned = 0.000334529 Jy ! Component: 650 - total flux cleaned = 0.000237999 Jy ! Component: 700 - total flux cleaned = 7.96826e-05 Jy ! Component: 750 - total flux cleaned = 4.91269e-05 Jy ! Component: 800 - total flux cleaned = 4.91083e-05 Jy ! Component: 850 - total flux cleaned = -1.2337e-05 Jy ! Component: 900 - total flux cleaned = -4.30497e-05 Jy ! Component: 950 - total flux cleaned = -0.000103548 Jy ! Component: 1000 - total flux cleaned = -0.000133315 Jy ! Component: 1050 - total flux cleaned = -0.000250717 Jy ! Component: 1100 - total flux cleaned = -0.000279653 Jy ! Component: 1150 - total flux cleaned = -0.000336902 Jy ! Component: 1200 - total flux cleaned = -0.000422095 Jy ! Component: 1250 - total flux cleaned = -0.0003649 Jy ! Component: 1300 - total flux cleaned = -0.000337353 Jy ! Component: 1350 - total flux cleaned = -0.000393707 Jy ! Component: 1400 - total flux cleaned = -0.00050418 Jy ! Component: 1450 - total flux cleaned = -0.000531417 Jy ! Component: 1500 - total flux cleaned = -0.000531292 Jy ! Component: 1550 - total flux cleaned = -0.000639425 Jy ! Component: 1600 - total flux cleaned = -0.000719586 Jy ! Component: 1650 - total flux cleaned = -0.000693006 Jy ! Component: 1700 - total flux cleaned = -0.000798402 Jy ! Component: 1750 - total flux cleaned = -0.0008246 Jy ! Component: 1800 - total flux cleaned = -0.000902645 Jy ! Component: 1850 - total flux cleaned = -0.000954044 Jy ! Component: 1900 - total flux cleaned = -0.00100536 Jy ! Component: 1950 - total flux cleaned = -0.000980053 Jy ! Component: 2000 - total flux cleaned = -0.00118279 Jy ! Component: 2050 - total flux cleaned = -0.0011573 Jy ! Component: 2100 - total flux cleaned = -0.00110756 Jy ! Component: 2150 - total flux cleaned = -0.00113258 Jy ! Component: 2200 - total flux cleaned = -0.00133007 Jy ! Component: 2250 - total flux cleaned = -0.00125621 Jy ! Component: 2300 - total flux cleaned = -0.00137811 Jy ! Component: 2350 - total flux cleaned = -0.00137797 Jy ! Component: 2400 - total flux cleaned = -0.00140205 Jy ! Component: 2450 - total flux cleaned = -0.00140187 Jy ! Component: 2500 - total flux cleaned = -0.00159227 Jy ! Component: 2550 - total flux cleaned = -0.00156836 Jy ! Component: 2600 - total flux cleaned = -0.00154499 Jy ! Component: 2650 - total flux cleaned = -0.00156848 Jy ! Component: 2700 - total flux cleaned = -0.00170839 Jy ! Component: 2750 - total flux cleaned = -0.00177769 Jy ! Component: 2800 - total flux cleaned = -0.00182369 Jy ! Component: 2850 - total flux cleaned = -0.00186951 Jy ! Component: 2900 - total flux cleaned = -0.0019379 Jy ! Component: 2950 - total flux cleaned = -0.00198334 Jy ! Component: 3000 - total flux cleaned = -0.00207355 Jy ! Component: 3050 - total flux cleaned = -0.00225302 Jy ! Component: 3100 - total flux cleaned = -0.00225303 Jy ! Component: 3150 - total flux cleaned = -0.0023417 Jy ! Component: 3200 - total flux cleaned = -0.00249628 Jy ! Component: 3250 - total flux cleaned = -0.00249629 Jy ! Component: 3300 - total flux cleaned = -0.00256178 Jy ! Component: 3350 - total flux cleaned = -0.00249651 Jy ! Component: 3400 - total flux cleaned = -0.0025615 Jy ! Component: 3450 - total flux cleaned = -0.00264775 Jy ! Component: 3500 - total flux cleaned = -0.00264792 Jy ! Component: 3550 - total flux cleaned = -0.00273347 Jy ! Component: 3600 - total flux cleaned = -0.00281844 Jy ! Component: 3650 - total flux cleaned = -0.00283958 Jy ! Component: 3700 - total flux cleaned = -0.0028817 Jy ! Component: 3750 - total flux cleaned = -0.00296551 Jy ! Component: 3800 - total flux cleaned = -0.00298623 Jy ! Component: 3850 - total flux cleaned = -0.00306921 Jy ! Component: 3900 - total flux cleaned = -0.0031104 Jy ! Component: 3950 - total flux cleaned = -0.00317208 Jy ! Component: 4000 - total flux cleaned = -0.00315168 Jy ! Component: 4050 - total flux cleaned = -0.00329454 Jy ! Total flux subtracted in 4096 components = -0.00333522 Jy ! Clean residual min=-0.001190 max=0.001500 Jy/beam ! Clean residual mean=0.000010 rms=0.000351 Jy/beam ! Combined flux in latest and established models = 0.833322 Jy keep ! Adding 309 model components to the UV plane model. ! The established model now contains 457 components and 0.833323 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /scr/pima/rdv91_uvs/J1753+2848_S_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.702 x 4.631 at -18.98 degrees (North through East) ! Clean map min=-0.0021312 max=0.76161 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv91_uvs/J1753+2848_S_map.fits wmodel /scr/pima/rdv91_uvs/J1753+2848_S_map.mod ! Writing 457 model components to file: /scr/pima/rdv91_uvs/J1753+2848_S_map.mod wobs /scr/pima/rdv91_uvs/J1753+2848_S_uvs.fits ! Writing UV FITS file: /scr/pima/rdv91_uvs/J1753+2848_S_uvs.fits wwins /scr/pima/rdv91_uvs/J1753+2848_S_map.win ! wwins: Wrote 15 windows to /scr/pima/rdv91_uvs/J1753+2848_S_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.761609 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00034964 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.137724 0.137724 0.275449 0.550897 1.10179 2.20359 4.40718 8.81436 17.6287 35.2574 70.5149 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.001 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.761 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 2178.26 quit ! Quitting program ! Log file /scr/pima/rdv91_uvs/J1753+2848_S_dfm.log closed on Fri Jun 15 11:29:04 2018