! Started logfile: /scr/pima/rdv76_uvs/J1354-0206_S_dfm.log on Fri Jun 15 05:55:03 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/rdv76_uvs/J1354-0206_S_uva.fits ! Reading UV FITS file: /scr/pima/rdv76_uvs/J1354-0206_S_uva.fits ! AN table 1: 11 integrations on 6 of 6 possible baselines. ! AN table 2: 27 integrations on 36 of 36 possible baselines. ! AN table 3: 77 integrations on 45 of 45 possible baselines. ! AN table 4: 45 integrations on 3 of 3 possible baselines. ! Apparent sampling: 0.428849 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1354-0206 ! ! 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 79 lines of history. ! ! Reading 7956 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 J1354-0206_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 12 telescope corrections were flagged in sub-array 1. ! A total of 78 telescope corrections were flagged in sub-array 2. ! A total of 476 telescope corrections were flagged in sub-array 3. ! A total of 84 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 2. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 78 telescope corrections were flagged in sub-array 2. ! A total of 476 telescope corrections were flagged in sub-array 3. ! A total of 84 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 3. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 78 telescope corrections were flagged in sub-array 2. ! A total of 476 telescope corrections were flagged in sub-array 3. ! A total of 84 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 4. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 78 telescope corrections were flagged in sub-array 2. ! A total of 476 telescope corrections were flagged in sub-array 3. ! A total of 84 telescope corrections were flagged in sub-array 4. ! ! Fit before self-cal, rms=1.348935Jy sigma=14.445921 ! Fit after self-cal, rms=0.269391Jy sigma=2.671046 ! 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.523 mas, bmaj=7.089 mas, bpa=-7.871 degrees ! Estimated noise=1.14127 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.648738 Jy ! Component: 100 - total flux cleaned = 0.800101 Jy ! Total flux subtracted in 100 components = 0.800101 Jy ! Clean residual min=-0.022900 max=0.057976 Jy/beam ! Clean residual mean=0.000081 rms=0.005820 Jy/beam ! Combined flux in latest and established models = 0.800101 Jy ! Performing phase self-cal ! Adding 5 model components to the UV plane model. ! The established model now contains 5 components and 0.800101 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.227998Jy sigma=2.359606 ! Fit after self-cal, rms=0.228014Jy sigma=2.359482 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.059408 Jy ! Component: 100 - total flux cleaned = 0.0854434 Jy ! Total flux subtracted in 100 components = 0.0854434 Jy ! Clean residual min=-0.016337 max=0.025232 Jy/beam ! Clean residual mean=0.000033 rms=0.004491 Jy/beam ! Combined flux in latest and established models = 0.885544 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 11 components and 0.885544 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.207502Jy sigma=2.108470 ! Fit after self-cal, rms=0.207482Jy sigma=2.108176 ! 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.896 mas, bmaj=8.083 mas, bpa=-7.995 degrees ! Estimated noise=0.895644 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 (-3, 10). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.031708 Jy ! Component: 100 - total flux cleaned = 0.0419782 Jy ! Component: 150 - total flux cleaned = 0.0470101 Jy ! Component: 200 - total flux cleaned = 0.0512529 Jy ! Total flux subtracted in 200 components = 0.0512529 Jy ! Clean residual min=-0.016041 max=0.017881 Jy/beam ! Clean residual mean=0.000025 rms=0.004271 Jy/beam ! Combined flux in latest and established models = 0.936797 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 37 components and 0.936797 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.202205Jy sigma=2.042740 ! Fit after self-cal, rms=0.201948Jy sigma=2.040326 ! 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.896 x 8.083 at -7.995 degrees (North through East) ! Clean map min=-0.016466 max=0.84606 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! HN 0.89 NL 1.05 OV 0.99 SC 0.93 ! ! A total of 9 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.03 FD 0.87 KP 1.00 LA 1.05 ! NL 0.98 OV 0.99 PT 0.99 SC 0.96 ! WZ 1.02 ! ! A total of 27 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.05 FD 0.88 KK 1.04 KP 1.37 ! LA 1.02 MK 0.95 NL 1.01 OV 0.99 ! PT 0.97 SC 0.97 ! ! A total of 51 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! NY 0.00* WZ 0.00* ZC 0.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! HN 0.96 NL 0.99 OV 1.01 SC 0.89 ! ! A total of 9 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.04 FD 0.87 KP 1.01 LA 1.01 ! NL 1.04 OV 1.01 PT 0.98 SC 0.97 ! WZ 1.05 ! ! A total of 27 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.04 FD 0.89 KK 1.04 KP 0.98 ! LA 1.01 MK 0.99 NL 1.02 OV 0.99 ! PT 0.98 SC 0.99 ! ! A total of 51 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! NY 0.00* WZ 0.00* ZC 0.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! HN 0.95 NL 1.07 OV 0.99 SC 0.90 ! ! A total of 9 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.06 FD 0.89 KP 1.00 LA 1.02 ! NL 1.00 OV 1.00 PT 1.01 SC 0.97 ! WZ 1.06 ! ! A total of 27 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.03 FD 0.89 KK 1.07 KP 0.99 ! LA 1.01 MK 0.99 NL 1.02 OV 0.98 ! PT 0.98 SC 0.97 ! ! A total of 51 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! NY 0.00* WZ 0.00* ZC 0.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! HN 0.86 NL 1.11 OV 1.02 SC 0.91 ! ! A total of 9 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.05 FD 0.89 KP 1.00 LA 1.00 ! NL 1.00 OV 1.03 PT 1.00 SC 0.98 ! WZ 1.03 ! ! A total of 27 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.09 FD 0.87 KK 1.08 KP 1.01 ! LA 1.01 MK 0.98 NL 1.03 OV 0.99 ! PT 0.96 SC 1.04 ! ! A total of 51 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! NY 0.00* WZ 0.00* ZC 0.00* ! ! ! Fit before self-cal, rms=0.201948Jy sigma=2.040326 ! Fit after self-cal, rms=0.192733Jy sigma=1.901834 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.937 mas, bmaj=8.114 mas, bpa=-7.267 degrees ! Estimated noise=0.899821 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.00162642 Jy ! Component: 100 - total flux cleaned = 0.0022589 Jy ! Component: 150 - total flux cleaned = 0.00261757 Jy ! Component: 200 - total flux cleaned = 0.00292841 Jy ! Total flux subtracted in 200 components = 0.00292841 Jy ! Clean residual min=-0.008573 max=0.009296 Jy/beam ! Clean residual mean=-0.000001 rms=0.002173 Jy/beam ! Combined flux in latest and established models = 0.939725 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 55 components and 0.939725 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.192230Jy sigma=1.896566 ! Fit after self-cal, rms=0.192311Jy sigma=1.895699 ! 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.192311Jy sigma=1.895699 ! Fit after self-cal, rms=0.190726Jy sigma=1.885341 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.938 mas, bmaj=8.115 mas, bpa=-7.23 degrees ! Estimated noise=0.89872 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000508683 Jy ! Component: 100 - total flux cleaned = 0.000654339 Jy ! Component: 150 - total flux cleaned = 0.000651017 Jy ! Component: 200 - total flux cleaned = 0.000768785 Jy ! Total flux subtracted in 200 components = 0.000768785 Jy ! Clean residual min=-0.008227 max=0.008224 Jy/beam ! Clean residual mean=0.000006 rms=0.002095 Jy/beam ! Combined flux in latest and established models = 0.940494 Jy selfcal ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 64 components and 0.940494 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.190465Jy sigma=1.882977 ! Fit after self-cal, rms=0.190424Jy sigma=1.882793 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.000755864 Jy ! Component: 100 - total flux cleaned = 0.00106694 Jy ! Component: 150 - total flux cleaned = 0.00115492 Jy ! Component: 200 - total flux cleaned = 0.00124066 Jy ! Total flux subtracted in 200 components = 0.00124066 Jy ! Clean residual min=-0.007992 max=0.008114 Jy/beam ! Clean residual mean=0.000005 rms=0.002076 Jy/beam ! Combined flux in latest and established models = 0.941735 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 72 components and 0.941735 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.190308Jy sigma=1.881811 ! Fit after self-cal, rms=0.190292Jy sigma=1.881754 ! 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.190292Jy sigma=1.881754 ! Fit after self-cal, rms=0.182560Jy sigma=1.780283 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.938 mas, bmaj=8.118 mas, bpa=-7.161 degrees ! Estimated noise=0.898392 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00102509 Jy ! Component: 100 - total flux cleaned = 0.00114043 Jy ! Component: 150 - total flux cleaned = 0.00133865 Jy ! Component: 200 - total flux cleaned = 0.00151432 Jy ! Total flux subtracted in 200 components = 0.00151432 Jy ! Clean residual min=-0.005936 max=0.006558 Jy/beam ! Clean residual mean=0.000012 rms=0.001574 Jy/beam ! Combined flux in latest and established models = 0.943249 Jy selfcal ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 77 components and 0.943249 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.182410Jy sigma=1.778963 ! Fit after self-cal, rms=0.182521Jy sigma=1.778691 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.000616854 Jy ! Component: 100 - total flux cleaned = 0.000773674 Jy ! Component: 150 - total flux cleaned = 0.00084425 Jy ! Component: 200 - total flux cleaned = 0.000974314 Jy ! Total flux subtracted in 200 components = 0.000974314 Jy ! Clean residual min=-0.006060 max=0.006570 Jy/beam ! Clean residual mean=0.000011 rms=0.001559 Jy/beam ! Combined flux in latest and established models = 0.944223 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 85 components and 0.944223 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.182452Jy sigma=1.778134 ! Fit after self-cal, rms=0.182447Jy sigma=1.778109 ! 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.182447Jy sigma=1.778109 ! Fit after self-cal, rms=0.182396Jy sigma=1.777738 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=9.338 mas, bmaj=12.84 mas, bpa=7.653 degrees ! Estimated noise=1.10291 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000951042 Jy ! Component: 100 - total flux cleaned = 0.00104352 Jy ! Component: 150 - total flux cleaned = 0.00104251 Jy ! Component: 200 - total flux cleaned = 0.00104215 Jy ! Total flux subtracted in 200 components = 0.00104215 Jy ! Clean residual min=-0.006192 max=0.006632 Jy/beam ! Clean residual mean=0.000024 rms=0.001934 Jy/beam ! Combined flux in latest and established models = 0.945265 Jy selfcal ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 94 components and 0.945266 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.182390Jy sigma=1.777424 ! Fit after self-cal, rms=0.182381Jy sigma=1.777373 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.000211941 Jy ! Component: 100 - total flux cleaned = 0.000310243 Jy ! Component: 150 - total flux cleaned = 0.000366919 Jy ! Component: 200 - total flux cleaned = 0.000392625 Jy ! Total flux subtracted in 200 components = 0.000392625 Jy ! Clean residual min=-0.006111 max=0.006498 Jy/beam ! Clean residual mean=0.000024 rms=0.001926 Jy/beam ! Combined flux in latest and established models = 0.945658 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 101 components and 0.945658 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.182376Jy sigma=1.777251 ! Fit after self-cal, rms=0.182371Jy sigma=1.777237 ! 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.939 mas, bmaj=8.116 mas, bpa=-7.103 degrees ! Estimated noise=0.898612 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.939 x 8.116 at -7.103 degrees (North through East) ! Clean map min=-0.0055995 max=0.84401 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.182371Jy sigma=1.777237 ! Fit after self-cal, rms=0.169952Jy sigma=1.575658 ! Inverting map and beam ! Estimated beam: bmin=2.938 mas, bmaj=8.116 mas, bpa=-7.084 degrees ! Estimated noise=0.898415 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00105579 Jy ! Component: 100 - total flux cleaned = -0.000903049 Jy ! Component: 150 - total flux cleaned = -0.000903816 Jy ! Component: 200 - total flux cleaned = -0.000844584 Jy ! Total flux subtracted in 200 components = -0.000844584 Jy ! Clean residual min=-0.005716 max=0.006296 Jy/beam ! Clean residual mean=0.000012 rms=0.001470 Jy/beam ! Combined flux in latest and established models = 0.944814 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 104 components and 0.944814 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.169864Jy sigma=1.575010 ! Fit after self-cal, rms=0.169866Jy sigma=1.574967 ! Inverting map ! Added new window around map position (-23.5, 68). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00122417 Jy ! Component: 100 - total flux cleaned = 0.00170232 Jy ! Component: 150 - total flux cleaned = 0.00137 Jy ! Component: 200 - total flux cleaned = 0.00116458 Jy ! Total flux subtracted in 200 components = 0.00116458 Jy ! Clean residual min=-0.005241 max=0.006060 Jy/beam ! Clean residual mean=0.000007 rms=0.001405 Jy/beam ! Combined flux in latest and established models = 0.945978 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 119 components and 0.945978 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.169584Jy sigma=1.571729 ! Fit after self-cal, rms=0.169559Jy sigma=1.571470 ! Inverting map ! Added new window around map position (103, 80). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00414868 Jy ! Component: 100 - total flux cleaned = 0.00494015 Jy ! Component: 150 - total flux cleaned = 0.00533171 Jy ! Component: 200 - total flux cleaned = 0.00542692 Jy ! Total flux subtracted in 200 components = 0.00542692 Jy ! Clean residual min=-0.004890 max=0.005770 Jy/beam ! Clean residual mean=0.000004 rms=0.001356 Jy/beam ! Combined flux in latest and established models = 0.951405 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 134 components and 0.951405 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.169311Jy sigma=1.569490 ! Fit after self-cal, rms=0.169294Jy sigma=1.569307 ! Inverting map ! Added new window around map position (-26, -83). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00474371 Jy ! Component: 100 - total flux cleaned = 0.00639516 Jy ! Component: 150 - total flux cleaned = 0.00808532 Jy ! Component: 200 - total flux cleaned = 0.009117 Jy ! Total flux subtracted in 200 components = 0.009117 Jy ! Clean residual min=-0.005047 max=0.005840 Jy/beam ! Clean residual mean=0.000009 rms=0.001301 Jy/beam ! Combined flux in latest and established models = 0.960522 Jy ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 156 components and 0.960522 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.169056Jy sigma=1.566889 ! Fit after self-cal, rms=0.169042Jy sigma=1.566772 ! Inverting map ! Added new window around map position (48.5, 125.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00385354 Jy ! Component: 100 - total flux cleaned = 0.00577282 Jy ! Component: 150 - total flux cleaned = 0.00728107 Jy ! Component: 200 - total flux cleaned = 0.00831162 Jy ! Total flux subtracted in 200 components = 0.00831162 Jy ! Clean residual min=-0.005031 max=0.005043 Jy/beam ! Clean residual mean=0.000008 rms=0.001248 Jy/beam ! Combined flux in latest and established models = 0.968834 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 172 components and 0.968834 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.168852Jy sigma=1.564784 ! Fit after self-cal, rms=0.168836Jy sigma=1.564692 ! 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.168836Jy sigma=1.564692 ! Fit after self-cal, rms=0.168664Jy sigma=1.563805 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.493 mas, bmaj=7.096 mas, bpa=-7.471 degrees ! Estimated noise=1.14312 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.65214 Jy ! Component: 100 - total flux cleaned = 0.801969 Jy ! Component: 150 - total flux cleaned = 0.859817 Jy ! Component: 200 - total flux cleaned = 0.887767 Jy ! Total flux subtracted in 200 components = 0.887767 Jy ! Clean residual min=-0.010826 max=0.014188 Jy/beam ! Clean residual mean=0.000025 rms=0.002395 Jy/beam ! Combined flux in latest and established models = 0.887767 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 9 model components to the UV plane model. ! The established model now contains 9 components and 0.887767 Jy ! Inverting map and beam ! Estimated beam: bmin=2.94 mas, bmaj=8.116 mas, bpa=-7.048 degrees ! Estimated noise=0.897687 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.0236821 Jy ! Component: 100 - total flux cleaned = 0.0364008 Jy ! Component: 150 - total flux cleaned = 0.0433998 Jy ! Component: 200 - total flux cleaned = 0.0496541 Jy ! Total flux subtracted in 200 components = 0.0496541 Jy ! Clean residual min=-0.006148 max=0.005931 Jy/beam ! Clean residual mean=0.000024 rms=0.001581 Jy/beam ! Combined flux in latest and established models = 0.937421 Jy ! Adding 39 model components to the UV plane model. ! The established model now contains 44 components and 0.937421 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=9.342 mas, bmaj=12.86 mas, bpa=7.77 degrees ! Estimated noise=1.10108 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00644025 Jy ! Component: 100 - total flux cleaned = 0.010505 Jy ! Component: 150 - total flux cleaned = 0.0135835 Jy ! Component: 200 - total flux cleaned = 0.0164248 Jy ! Total flux subtracted in 200 components = 0.0164248 Jy ! Clean residual min=-0.005011 max=0.006362 Jy/beam ! Clean residual mean=0.000030 rms=0.001604 Jy/beam ! Combined flux in latest and established models = 0.953846 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.00247632 Jy ! Component: 100 - total flux cleaned = 0.00465322 Jy ! Component: 150 - total flux cleaned = 0.0064748 Jy ! Component: 200 - total flux cleaned = 0.00771186 Jy ! Total flux subtracted in 200 components = 0.00771186 Jy ! Clean residual min=-0.004719 max=0.005838 Jy/beam ! Clean residual mean=0.000025 rms=0.001486 Jy/beam ! Combined flux in latest and established models = 0.961558 Jy ! Adding 63 model components to the UV plane model. ! The established model now contains 104 components and 0.961558 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.94 mas, bmaj=8.116 mas, bpa=-7.048 degrees ! Estimated noise=0.897687 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 6.98634e-06 Jy ! Component: 100 - total flux cleaned = 0.00174459 Jy ! Component: 150 - total flux cleaned = 0.0025164 Jy ! Component: 200 - total flux cleaned = 0.00323031 Jy ! Total flux subtracted in 200 components = 0.00323031 Jy ! Clean residual min=-0.005080 max=0.005210 Jy/beam ! Clean residual mean=0.000010 rms=0.001312 Jy/beam ! Combined flux in latest and established models = 0.964788 Jy selfcal ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 130 components and 0.964788 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.169246Jy sigma=1.568186 ! Fit after self-cal, rms=0.169207Jy sigma=1.567660 wmodel J1354-0206_S_map.mod ! Writing 130 model components to file: J1354-0206_S_map.mod wobs J1354-0206_S_uvs.fits ! Writing UV FITS file: J1354-0206_S_uvs.fits wwins J1354-0206_S_map.win ! wwins: Wrote 6 windows to J1354-0206_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.00192571 Jy ! Component: 100 - total flux cleaned = 0.00281929 Jy ! Component: 150 - total flux cleaned = 0.00333139 Jy ! Component: 200 - total flux cleaned = 0.00365857 Jy ! Component: 250 - total flux cleaned = 0.00381728 Jy ! Component: 300 - total flux cleaned = 0.00397394 Jy ! Component: 350 - total flux cleaned = 0.0038973 Jy ! Component: 400 - total flux cleaned = 0.00397298 Jy ! Component: 450 - total flux cleaned = 0.00375309 Jy ! Component: 500 - total flux cleaned = 0.00353658 Jy ! Component: 550 - total flux cleaned = 0.00339472 Jy ! Component: 600 - total flux cleaned = 0.0033243 Jy ! Component: 650 - total flux cleaned = 0.00346272 Jy ! Component: 700 - total flux cleaned = 0.0032571 Jy ! Component: 750 - total flux cleaned = 0.00339221 Jy ! Component: 800 - total flux cleaned = 0.00332513 Jy ! Component: 850 - total flux cleaned = 0.00339087 Jy ! Component: 900 - total flux cleaned = 0.00306303 Jy ! Component: 950 - total flux cleaned = 0.00293402 Jy ! Component: 1000 - total flux cleaned = 0.00299751 Jy ! Component: 1050 - total flux cleaned = 0.00280663 Jy ! Component: 1100 - total flux cleaned = 0.00261793 Jy ! Component: 1150 - total flux cleaned = 0.00261755 Jy ! Component: 1200 - total flux cleaned = 0.00249501 Jy ! Component: 1250 - total flux cleaned = 0.00243466 Jy ! Component: 1300 - total flux cleaned = 0.0026164 Jy ! Component: 1350 - total flux cleaned = 0.00273721 Jy ! Component: 1400 - total flux cleaned = 0.00279752 Jy ! Component: 1450 - total flux cleaned = 0.00279739 Jy ! Component: 1500 - total flux cleaned = 0.0028551 Jy ! Component: 1550 - total flux cleaned = 0.00273836 Jy ! Component: 1600 - total flux cleaned = 0.00291118 Jy ! Component: 1650 - total flux cleaned = 0.00273889 Jy ! Component: 1700 - total flux cleaned = 0.00251163 Jy ! Component: 1750 - total flux cleaned = 0.00273742 Jy ! Component: 1800 - total flux cleaned = 0.00240138 Jy ! Component: 1850 - total flux cleaned = 0.00245713 Jy ! Component: 1900 - total flux cleaned = 0.00240193 Jy ! Component: 1950 - total flux cleaned = 0.00229282 Jy ! Component: 2000 - total flux cleaned = 0.00229245 Jy ! Component: 2050 - total flux cleaned = 0.00240087 Jy ! Component: 2100 - total flux cleaned = 0.00229339 Jy ! Component: 2150 - total flux cleaned = 0.00218763 Jy ! Component: 2200 - total flux cleaned = 0.00208164 Jy ! Component: 2250 - total flux cleaned = 0.00197729 Jy ! Component: 2300 - total flux cleaned = 0.00187323 Jy ! Component: 2350 - total flux cleaned = 0.00208034 Jy ! Component: 2400 - total flux cleaned = 0.0018222 Jy ! Component: 2450 - total flux cleaned = 0.00171997 Jy ! Component: 2500 - total flux cleaned = 0.00197394 Jy ! Component: 2550 - total flux cleaned = 0.00192333 Jy ! Component: 2600 - total flux cleaned = 0.00192352 Jy ! Component: 2650 - total flux cleaned = 0.00197385 Jy ! Component: 2700 - total flux cleaned = 0.00212299 Jy ! Component: 2750 - total flux cleaned = 0.00212294 Jy ! Component: 2800 - total flux cleaned = 0.00207363 Jy ! Component: 2850 - total flux cleaned = 0.00212278 Jy ! Component: 2900 - total flux cleaned = 0.00207402 Jy ! Component: 2950 - total flux cleaned = 0.00212266 Jy ! Component: 3000 - total flux cleaned = 0.00217047 Jy ! Component: 3050 - total flux cleaned = 0.0021233 Jy ! Component: 3100 - total flux cleaned = 0.00202861 Jy ! Component: 3150 - total flux cleaned = 0.00207578 Jy ! Component: 3200 - total flux cleaned = 0.00207601 Jy ! Component: 3250 - total flux cleaned = 0.00184289 Jy ! Component: 3300 - total flux cleaned = 0.00193507 Jy ! Component: 3350 - total flux cleaned = 0.00193491 Jy ! Component: 3400 - total flux cleaned = 0.00193465 Jy ! Component: 3450 - total flux cleaned = 0.00175203 Jy ! Component: 3500 - total flux cleaned = 0.00179773 Jy ! Component: 3550 - total flux cleaned = 0.00170759 Jy ! Component: 3600 - total flux cleaned = 0.00184207 Jy ! Component: 3650 - total flux cleaned = 0.0017974 Jy ! Component: 3700 - total flux cleaned = 0.00170863 Jy ! Component: 3750 - total flux cleaned = 0.00175293 Jy ! Component: 3800 - total flux cleaned = 0.00170908 Jy ! Component: 3850 - total flux cleaned = 0.00166545 Jy ! Component: 3900 - total flux cleaned = 0.00179614 Jy ! Component: 3950 - total flux cleaned = 0.00157996 Jy ! Component: 4000 - total flux cleaned = 0.00157973 Jy ! Component: 4050 - total flux cleaned = 0.00157994 Jy ! Total flux subtracted in 4096 components = 0.00166548 Jy ! Clean residual min=-0.003413 max=0.002567 Jy/beam ! Clean residual mean=0.000007 rms=0.000773 Jy/beam ! Combined flux in latest and established models = 0.966454 Jy keep ! Adding 286 model components to the UV plane model. ! The established model now contains 416 components and 0.966454 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /scr/pima/rdv76_uvs/J1354-0206_S_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.94 x 8.116 at -7.048 degrees (North through East) ! Clean map min=-0.0059945 max=0.84509 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv76_uvs/J1354-0206_S_map.fits wmodel /scr/pima/rdv76_uvs/J1354-0206_S_map.mod ! Writing 416 model components to file: /scr/pima/rdv76_uvs/J1354-0206_S_map.mod wobs /scr/pima/rdv76_uvs/J1354-0206_S_uvs.fits ! Writing UV FITS file: /scr/pima/rdv76_uvs/J1354-0206_S_uvs.fits wwins /scr/pima/rdv76_uvs/J1354-0206_S_map.win ! wwins: Wrote 7 windows to /scr/pima/rdv76_uvs/J1354-0206_S_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.845094 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000771517 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.273881 0.273881 0.547762 1.09552 2.19105 4.38209 8.76418 17.5284 35.0567 70.1135 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0023 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.845 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 1095.37 quit ! Quitting program ! Log file /scr/pima/rdv76_uvs/J1354-0206_S_dfm.log closed on Fri Jun 15 05:55:07 2018