! Started logfile: /scr/pima/rv123_uvs/J1939-1525_S_dfm.log on Sun Dec 23 05:05:53 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/rv123_uvs/J1939-1525_S_uva.fits ! Reading UV FITS file: /scr/pima/rv123_uvs/J1939-1525_S_uva.fits ! AN table 1: 7 integrations on 15 of 15 possible baselines. ! AN table 2: 56 integrations on 21 of 21 possible baselines. ! AN table 3: 7 integrations on 21 of 21 possible baselines. ! AN table 4: 3 integrations on 15 of 15 possible baselines. ! AN table 5: 10 integrations on 36 of 36 possible baselines. ! Apparent sampling: 0.424441 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1939-1525 ! ! 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 89 lines of history. ! ! Reading 3112 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 J1939-1525_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 28 telescope corrections were flagged in sub-array 1. ! A total of 258 telescope corrections were flagged in sub-array 2. ! A total of 19 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 46 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 28 telescope corrections were flagged in sub-array 1. ! A total of 258 telescope corrections were flagged in sub-array 2. ! A total of 19 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 46 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 28 telescope corrections were flagged in sub-array 1. ! A total of 314 telescope corrections were flagged in sub-array 2. ! A total of 49 telescope corrections were flagged in sub-array 3. ! A total of 18 telescope corrections were flagged in sub-array 4. ! A total of 90 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 28 telescope corrections were flagged in sub-array 1. ! A total of 258 telescope corrections were flagged in sub-array 2. ! A total of 19 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 46 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=1.326255Jy sigma=13.443611 ! Fit after self-cal, rms=0.320639Jy sigma=2.931673 ! 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.242 mas, bmaj=3.011 mas, bpa=16.32 degrees ! Estimated noise=2.24785 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.596531 Jy ! Component: 100 - total flux cleaned = 0.758245 Jy ! Total flux subtracted in 100 components = 0.758245 Jy ! Clean residual min=-0.027983 max=0.059568 Jy/beam ! Clean residual mean=0.000096 rms=0.008274 Jy/beam ! Combined flux in latest and established models = 0.758245 Jy ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 7 components and 0.758245 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.287087Jy sigma=2.144060 ! Fit after self-cal, rms=0.286721Jy sigma=2.143608 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0601299 Jy ! Component: 100 - total flux cleaned = 0.0932353 Jy ! Total flux subtracted in 100 components = 0.0932353 Jy ! Clean residual min=-0.027306 max=0.040513 Jy/beam ! Clean residual mean=0.000057 rms=0.007273 Jy/beam ! Combined flux in latest and established models = 0.85148 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 18 components and 0.85148 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.256002Jy sigma=1.864019 ! Fit after self-cal, rms=0.255819Jy sigma=1.863780 ! 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.615 mas, bmaj=3.57 mas, bpa=2.93 degrees ! Estimated noise=1.95954 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 (-1, -7). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0487819 Jy ! Component: 100 - total flux cleaned = 0.0763876 Jy ! Component: 150 - total flux cleaned = 0.0933643 Jy ! Component: 200 - total flux cleaned = 0.100362 Jy ! Total flux subtracted in 200 components = 0.100362 Jy ! Clean residual min=-0.027760 max=0.037058 Jy/beam ! Clean residual mean=0.000081 rms=0.006627 Jy/beam ! Combined flux in latest and established models = 0.951843 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 33 components and 0.951843 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.236482Jy sigma=1.767955 ! Fit after self-cal, rms=0.235677Jy sigma=1.756870 ! 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.615 x 3.57 at 2.93 degrees (North through East) ! Clean map min=-0.025231 max=0.76751 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: ! HH 1.25 MC 1.12 ON 0.97 WN 1.00* ! WZ 0.92 YG 1.00* ! ! A total of 36 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! HH 1.08 HN 1.08 MC 1.14 ON 0.88 ! SC 1.01 WN 1.00* WZ 1.00 ! ! A total of 6 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! HN 0.94 LA 0.93 MC 1.16 NL 1.13 ! SC 0.90 WN 1.00* WZ 1.02 ! ! Telescope amplitude corrections in sub-array 4: ! BR 1.39 FD 1.10 KP 0.94 NL 1.08 ! OV 0.95 SC 0.89 ! ! A total of 6 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.08 FD 1.07 HN 0.95 KP 0.93 ! LA 0.95 MK 1.06 NL 1.10 OV 0.96 ! PT 1.09 ! ! ! Correcting IF 2. ! A total of 6 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! HH 1.02 MC 1.06 ON 1.05 WN 1.00* ! WZ 0.86 YG 1.00* ! ! A total of 36 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! HH 1.04 HN 1.14 MC 1.00 ON 1.07 ! SC 0.97 WN 1.00* WZ 0.90 ! ! A total of 6 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! HN 0.94 LA 0.95 MC 0.94 NL 1.12 ! SC 0.90 WN 1.00* WZ 0.96 ! ! Telescope amplitude corrections in sub-array 4: ! BR 2.73 FD 1.06 KP 0.97 NL 1.08 ! OV 0.75 SC 0.88 ! ! A total of 6 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.40 FD 1.10 HN 0.97 KP 0.93 ! LA 0.93 MK 1.08 NL 1.12 OV 0.74 ! PT 1.05 ! ! ! Correcting IF 3. ! A total of 6 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! HH 1.05 MC 1.23 ON 0.68 WN 1.00* ! WZ 0.89 YG 1.00* ! ! A total of 54 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! HH 1.04 HN 1.00* MC 1.11 ON 0.78 ! SC 1.00* WN 1.00* WZ 0.84 ! ! Telescope amplitude corrections in sub-array 3: ! HN 0.00* LA 0.00* MC 0.00* NL 0.00* ! SC 0.00* WN 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* FD 0.00* KP 0.00* NL 0.00* ! OV 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 0.00* FD 0.00* HN 0.00* KP 0.00* ! LA 0.00* MK 0.00* NL 0.00* OV 0.00* ! PT 0.00* ! ! ! Correcting IF 4. ! A total of 6 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! HH 0.88 MC 1.43 ON 0.74 WN 1.00* ! WZ 0.70 YG 1.00* ! ! A total of 36 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! HH 0.93 HN 1.14 MC 1.11 ON 0.80 ! SC 1.00 WN 1.00* WZ 0.88 ! ! A total of 6 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! HN 0.94 LA 0.94 MC 1.13 NL 1.10 ! SC 0.90 WN 1.00* WZ 0.79 ! ! Telescope amplitude corrections in sub-array 4: ! BR 1.66 FD 1.06 KP 0.96 NL 1.08 ! OV 0.92 SC 0.88 ! ! A total of 6 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.14 FD 1.04 HN 0.97 KP 0.95 ! LA 0.91 MK 1.06 NL 1.15 OV 0.89 ! PT 1.03 ! ! ! Fit before self-cal, rms=0.235677Jy sigma=1.756870 ! Fit after self-cal, rms=0.147655Jy sigma=0.976580 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.505 mas, bmaj=3.749 mas, bpa=1.633 degrees ! Estimated noise=2.02433 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.00333979 Jy ! Component: 100 - total flux cleaned = -0.00264556 Jy ! Component: 150 - total flux cleaned = -0.0028393 Jy ! Component: 200 - total flux cleaned = -0.00283337 Jy ! Total flux subtracted in 200 components = -0.00283337 Jy ! Clean residual min=-0.011493 max=0.011912 Jy/beam ! Clean residual mean=0.000042 rms=0.002651 Jy/beam ! Combined flux in latest and established models = 0.949009 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 51 components and 0.949009 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.146821Jy sigma=0.966468 ! Fit after self-cal, rms=0.146590Jy sigma=0.964411 ! 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.146590Jy sigma=0.964411 ! Fit after self-cal, rms=0.144604Jy sigma=0.947364 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.504 mas, bmaj=3.751 mas, bpa=1.483 degrees ! Estimated noise=2.02297 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00171631 Jy ! Component: 100 - total flux cleaned = -0.0017324 Jy ! Component: 150 - total flux cleaned = -0.00203962 Jy ! Component: 200 - total flux cleaned = -0.00203488 Jy ! Total flux subtracted in 200 components = -0.00203488 Jy ! Clean residual min=-0.011509 max=0.011057 Jy/beam ! Clean residual mean=0.000024 rms=0.002483 Jy/beam ! Combined flux in latest and established models = 0.946974 Jy selfcal ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 56 components and 0.946974 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.144112Jy sigma=0.942077 ! Fit after self-cal, rms=0.144116Jy sigma=0.941616 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.000986537 Jy ! Component: 100 - total flux cleaned = -0.000987419 Jy ! Component: 150 - total flux cleaned = -0.000985024 Jy ! Component: 200 - total flux cleaned = -0.000986466 Jy ! Total flux subtracted in 200 components = -0.000986466 Jy ! Clean residual min=-0.011608 max=0.010908 Jy/beam ! Clean residual mean=0.000024 rms=0.002457 Jy/beam ! Combined flux in latest and established models = 0.945988 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 59 components and 0.945988 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.143863Jy sigma=0.939463 ! Fit after self-cal, rms=0.143841Jy sigma=0.939257 ! 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.143841Jy sigma=0.939257 ! Fit after self-cal, rms=0.137681Jy sigma=0.880205 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.507 mas, bmaj=3.746 mas, bpa=1.391 degrees ! Estimated noise=2.02279 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00372806 Jy ! Component: 100 - total flux cleaned = -0.00353526 Jy ! Component: 150 - total flux cleaned = -0.00338016 Jy ! Component: 200 - total flux cleaned = -0.00344762 Jy ! Total flux subtracted in 200 components = -0.00344762 Jy ! Clean residual min=-0.007334 max=0.006544 Jy/beam ! Clean residual mean=0.000007 rms=0.001598 Jy/beam ! Combined flux in latest and established models = 0.94254 Jy selfcal ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 71 components and 0.94254 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.137565Jy sigma=0.877380 ! Fit after self-cal, rms=0.137650Jy sigma=0.876907 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.000933577 Jy ! Component: 100 - total flux cleaned = -0.000997484 Jy ! Component: 150 - total flux cleaned = -0.000939458 Jy ! Component: 200 - total flux cleaned = -0.000941502 Jy ! Total flux subtracted in 200 components = -0.000941502 Jy ! Clean residual min=-0.007380 max=0.006500 Jy/beam ! Clean residual mean=0.000007 rms=0.001584 Jy/beam ! Combined flux in latest and established models = 0.941599 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 75 components and 0.941599 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.137579Jy sigma=0.876097 ! Fit after self-cal, rms=0.137559Jy sigma=0.876007 ! 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.137559Jy sigma=0.876007 ! Fit after self-cal, rms=0.137231Jy sigma=0.875199 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.803 mas, bmaj=13.73 mas, bpa=-3.749 degrees ! Estimated noise=2.75327 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000590121 Jy ! Component: 100 - total flux cleaned = -0.00050977 Jy ! Component: 150 - total flux cleaned = -0.000473474 Jy ! Component: 200 - total flux cleaned = -0.000439599 Jy ! Total flux subtracted in 200 components = -0.000439599 Jy ! Clean residual min=-0.006566 max=0.007849 Jy/beam ! Clean residual mean=0.000008 rms=0.001951 Jy/beam ! Combined flux in latest and established models = 0.941159 Jy selfcal ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 82 components and 0.941159 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.137282Jy sigma=0.874885 ! Fit after self-cal, rms=0.137291Jy sigma=0.874816 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.000128018 Jy ! Component: 100 - total flux cleaned = -6.00761e-05 Jy ! Component: 150 - total flux cleaned = -2.85802e-05 Jy ! Component: 200 - total flux cleaned = 9.3681e-07 Jy ! Total flux subtracted in 200 components = 9.3681e-07 Jy ! Clean residual min=-0.006564 max=0.007795 Jy/beam ! Clean residual mean=0.000008 rms=0.001945 Jy/beam ! Combined flux in latest and established models = 0.94116 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 84 components and 0.94116 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.137332Jy sigma=0.874697 ! Fit after self-cal, rms=0.137332Jy sigma=0.874659 ! 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.507 mas, bmaj=3.739 mas, bpa=1.38 degrees ! Estimated noise=2.02063 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.507 x 3.739 at 1.38 degrees (North through East) ! Clean map min=-0.0067646 max=0.76733 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.137332Jy sigma=0.874659 ! Fit after self-cal, rms=0.126732Jy sigma=0.765960 ! Inverting map and beam ! Estimated beam: bmin=2.507 mas, bmaj=3.736 mas, bpa=1.374 degrees ! Estimated noise=2.02048 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00122576 Jy ! Component: 100 - total flux cleaned = -0.0015229 Jy ! Component: 150 - total flux cleaned = -0.00164844 Jy ! Component: 200 - total flux cleaned = -0.00164932 Jy ! Total flux subtracted in 200 components = -0.00164932 Jy ! Clean residual min=-0.007375 max=0.006230 Jy/beam ! Clean residual mean=0.000006 rms=0.001545 Jy/beam ! Combined flux in latest and established models = 0.939511 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 89 components and 0.939511 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.126629Jy sigma=0.764829 ! Fit after self-cal, rms=0.126598Jy sigma=0.764681 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000339527 Jy ! Component: 100 - total flux cleaned = -0.000500254 Jy ! Component: 150 - total flux cleaned = -0.000501821 Jy ! Component: 200 - total flux cleaned = -0.000545458 Jy ! Total flux subtracted in 200 components = -0.000545458 Jy ! Clean residual min=-0.007404 max=0.006300 Jy/beam ! Clean residual mean=0.000007 rms=0.001538 Jy/beam ! Combined flux in latest and established models = 0.938965 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 94 components and 0.938965 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.126529Jy sigma=0.764313 ! Fit after self-cal, rms=0.126527Jy sigma=0.764235 ! 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.126527Jy sigma=0.764235 ! Fit after self-cal, rms=0.126265Jy sigma=0.763848 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.207 mas, bmaj=3.129 mas, bpa=13.12 degrees ! Estimated noise=2.34881 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.608624 Jy ! Component: 100 - total flux cleaned = 0.773449 Jy ! Component: 150 - total flux cleaned = 0.834933 Jy ! Component: 200 - total flux cleaned = 0.869131 Jy ! Total flux subtracted in 200 components = 0.869131 Jy ! Clean residual min=-0.009629 max=0.018120 Jy/beam ! Clean residual mean=0.000031 rms=0.002491 Jy/beam ! Combined flux in latest and established models = 0.869132 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 38 model components to the UV plane model. ! The established model now contains 38 components and 0.869132 Jy ! Inverting map and beam ! Estimated beam: bmin=2.507 mas, bmaj=3.731 mas, bpa=1.336 degrees ! Estimated noise=2.0184 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.0283996 Jy ! Component: 100 - total flux cleaned = 0.0443927 Jy ! Component: 150 - total flux cleaned = 0.0538492 Jy ! Component: 200 - total flux cleaned = 0.0596079 Jy ! Total flux subtracted in 200 components = 0.0596079 Jy ! Clean residual min=-0.006953 max=0.006695 Jy/beam ! Clean residual mean=0.000015 rms=0.001581 Jy/beam ! Combined flux in latest and established models = 0.928739 Jy ! Adding 30 model components to the UV plane model. ! The established model now contains 59 components and 0.928739 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.808 mas, bmaj=13.74 mas, bpa=-3.761 degrees ! Estimated noise=2.75206 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00497819 Jy ! Component: 100 - total flux cleaned = 0.00686917 Jy ! Component: 150 - total flux cleaned = 0.00705925 Jy ! Component: 200 - total flux cleaned = 0.0072258 Jy ! Total flux subtracted in 200 components = 0.0072258 Jy ! Clean residual min=-0.006706 max=0.007743 Jy/beam ! Clean residual mean=0.000012 rms=0.001901 Jy/beam ! Combined flux in latest and established models = 0.935965 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.000146971 Jy ! Component: 100 - total flux cleaned = 0.000281475 Jy ! Component: 150 - total flux cleaned = 0.000375477 Jy ! Component: 200 - total flux cleaned = 0.000492238 Jy ! Total flux subtracted in 200 components = 0.000492238 Jy ! Clean residual min=-0.006686 max=0.007728 Jy/beam ! Clean residual mean=0.000011 rms=0.001896 Jy/beam ! Combined flux in latest and established models = 0.936457 Jy ! Adding 30 model components to the UV plane model. ! The established model now contains 83 components and 0.936458 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.507 mas, bmaj=3.731 mas, bpa=1.336 degrees ! Estimated noise=2.0184 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00221363 Jy ! Component: 100 - total flux cleaned = -0.00211862 Jy ! Component: 150 - total flux cleaned = -0.00203215 Jy ! Component: 200 - total flux cleaned = -0.00195575 Jy ! Total flux subtracted in 200 components = -0.00195575 Jy ! Clean residual min=-0.007218 max=0.006326 Jy/beam ! Clean residual mean=0.000007 rms=0.001548 Jy/beam ! Combined flux in latest and established models = 0.934502 Jy selfcal ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 91 components and 0.934502 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.126481Jy sigma=0.765939 ! Fit after self-cal, rms=0.126417Jy sigma=0.764991 wmodel J1939-1525_S_map.mod ! Writing 91 model components to file: J1939-1525_S_map.mod wobs J1939-1525_S_uvs.fits ! Writing UV FITS file: J1939-1525_S_uvs.fits wwins J1939-1525_S_map.win ! wwins: Wrote 2 windows to J1939-1525_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.00260297 Jy ! Component: 100 - total flux cleaned = -0.00248985 Jy ! Component: 150 - total flux cleaned = -0.00259738 Jy ! Component: 200 - total flux cleaned = -0.00270419 Jy ! Component: 250 - total flux cleaned = -0.00301047 Jy ! Component: 300 - total flux cleaned = -0.00340384 Jy ! Component: 350 - total flux cleaned = -0.00311194 Jy ! Component: 400 - total flux cleaned = -0.00330061 Jy ! Component: 450 - total flux cleaned = -0.00330137 Jy ! Component: 500 - total flux cleaned = -0.0038437 Jy ! Component: 550 - total flux cleaned = -0.00375374 Jy ! Component: 600 - total flux cleaned = -0.00384274 Jy ! Component: 650 - total flux cleaned = -0.0040161 Jy ! Component: 700 - total flux cleaned = -0.00410103 Jy ! Component: 750 - total flux cleaned = -0.00410302 Jy ! Component: 800 - total flux cleaned = -0.00402159 Jy ! Component: 850 - total flux cleaned = -0.00426777 Jy ! Component: 900 - total flux cleaned = -0.00426821 Jy ! Component: 950 - total flux cleaned = -0.00403001 Jy ! Component: 1000 - total flux cleaned = -0.00426506 Jy ! Component: 1050 - total flux cleaned = -0.00395519 Jy ! Component: 1100 - total flux cleaned = -0.00411078 Jy ! Component: 1150 - total flux cleaned = -0.00411106 Jy ! Component: 1200 - total flux cleaned = -0.00403667 Jy ! Component: 1250 - total flux cleaned = -0.00440825 Jy ! Component: 1300 - total flux cleaned = -0.00404048 Jy ! Component: 1350 - total flux cleaned = -0.00433083 Jy ! Component: 1400 - total flux cleaned = -0.00375345 Jy ! Component: 1450 - total flux cleaned = -0.00411054 Jy ! Component: 1500 - total flux cleaned = -0.00361538 Jy ! Component: 1550 - total flux cleaned = -0.00389559 Jy ! Component: 1600 - total flux cleaned = -0.00333991 Jy ! Component: 1650 - total flux cleaned = -0.0032028 Jy ! Component: 1700 - total flux cleaned = -0.00333942 Jy ! Component: 1750 - total flux cleaned = -0.00313756 Jy ! Component: 1800 - total flux cleaned = -0.00293631 Jy ! Component: 1850 - total flux cleaned = -0.00306905 Jy ! Component: 1900 - total flux cleaned = -0.00287056 Jy ! Component: 1950 - total flux cleaned = -0.00280437 Jy ! Component: 2000 - total flux cleaned = -0.0028031 Jy ! Component: 2050 - total flux cleaned = -0.00267441 Jy ! Component: 2100 - total flux cleaned = -0.00229144 Jy ! Component: 2150 - total flux cleaned = -0.00229194 Jy ! Component: 2200 - total flux cleaned = -0.00241678 Jy ! Component: 2250 - total flux cleaned = -0.00229247 Jy ! Component: 2300 - total flux cleaned = -0.00235455 Jy ! Component: 2350 - total flux cleaned = -0.0025378 Jy ! Component: 2400 - total flux cleaned = -0.00241526 Jy ! Component: 2450 - total flux cleaned = -0.00241568 Jy ! Component: 2500 - total flux cleaned = -0.00265573 Jy ! Component: 2550 - total flux cleaned = -0.00247717 Jy ! Component: 2600 - total flux cleaned = -0.00247686 Jy ! Component: 2650 - total flux cleaned = -0.00247694 Jy ! Component: 2700 - total flux cleaned = -0.00247601 Jy ! Component: 2750 - total flux cleaned = -0.00253413 Jy ! Component: 2800 - total flux cleaned = -0.00264905 Jy ! Component: 2850 - total flux cleaned = -0.00236361 Jy ! Component: 2900 - total flux cleaned = -0.00253422 Jy ! Component: 2950 - total flux cleaned = -0.00253333 Jy ! Component: 3000 - total flux cleaned = -0.00264481 Jy ! Component: 3050 - total flux cleaned = -0.00225507 Jy ! Component: 3100 - total flux cleaned = -0.00219999 Jy ! Component: 3150 - total flux cleaned = -0.00241952 Jy ! Component: 3200 - total flux cleaned = -0.00220173 Jy ! Component: 3250 - total flux cleaned = -0.00220174 Jy ! Component: 3300 - total flux cleaned = -0.00252522 Jy ! Component: 3350 - total flux cleaned = -0.0022573 Jy ! Component: 3400 - total flux cleaned = -0.00236343 Jy ! Component: 3450 - total flux cleaned = -0.00231089 Jy ! Component: 3500 - total flux cleaned = -0.00246791 Jy ! Component: 3550 - total flux cleaned = -0.00246798 Jy ! Component: 3600 - total flux cleaned = -0.00231203 Jy ! Component: 3650 - total flux cleaned = -0.00241513 Jy ! Component: 3700 - total flux cleaned = -0.00236374 Jy ! Component: 3750 - total flux cleaned = -0.00251627 Jy ! Component: 3800 - total flux cleaned = -0.00216103 Jy ! Component: 3850 - total flux cleaned = -0.00236276 Jy ! Component: 3900 - total flux cleaned = -0.0022623 Jy ! Component: 3950 - total flux cleaned = -0.00221234 Jy ! Component: 4000 - total flux cleaned = -0.00226165 Jy ! Component: 4050 - total flux cleaned = -0.00201573 Jy ! Total flux subtracted in 4096 components = -0.00201605 Jy ! Clean residual min=-0.002609 max=0.003205 Jy/beam ! Clean residual mean=0.000003 rms=0.000837 Jy/beam ! Combined flux in latest and established models = 0.932486 Jy keep ! Adding 300 model components to the UV plane model. ! The established model now contains 391 components and 0.932485 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /scr/pima/rv123_uvs/J1939-1525_S_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.507 x 3.731 at 1.336 degrees (North through East) ! Clean map min=-0.0049594 max=0.77847 Jy/beam ! Writing clean map to FITS file: /scr/pima/rv123_uvs/J1939-1525_S_map.fits wmodel /scr/pima/rv123_uvs/J1939-1525_S_map.mod ! Writing 391 model components to file: /scr/pima/rv123_uvs/J1939-1525_S_map.mod wobs /scr/pima/rv123_uvs/J1939-1525_S_uvs.fits ! Writing UV FITS file: /scr/pima/rv123_uvs/J1939-1525_S_uvs.fits wwins /scr/pima/rv123_uvs/J1939-1525_S_map.win ! wwins: Wrote 3 windows to /scr/pima/rv123_uvs/J1939-1525_S_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.778471 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000828158 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.319148 0.319148 0.638296 1.27659 2.55318 5.10637 10.2127 20.4255 40.8509 81.7019 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0024 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.778 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 940.003 quit ! Quitting program ! Log file /scr/pima/rv123_uvs/J1939-1525_S_dfm.log closed on Sun Dec 23 05:05:58 2018