! Started logfile: /image_prod/ug002_uva//J1217+5835_S_dfm.log on Thu Mar 5 00:33:45 2020 float field_size field_size = 1024 float field_cell field_cell = 0.5000 float taper_size taper_size = 20.0 float freq freq = 2.22000 obs /image_prod/ug002_uva//J1217+5835_S_uva.fits ! Reading UV FITS file: /image_prod/ug002_uva//J1217+5835_S_uva.fits ! AN table 1: 60 integrations on 45 of 45 possible baselines. ! Apparent sampling: 0.905926 visibilities/baseline/integration-bin. ! Found source: J1217+5835 ! ! There are 3 IFs, and a total of 3 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 2.22e+09 3.2e+07 1 3.2e+07 ! 02 2 2.252e+09 3.2e+07 1 3.2e+07 ! 03 3 2.284e+09 3.2e+07 1 3.2e+07 ! ! Polarization(s): RR ! ! Read 44 lines of history. ! ! Reading 7338 visibilities. select rr, 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16 ! Restricting channel ranges to the available 3 channels. ! Selecting polarization: RR, channels: 1..3 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 print "field_size =", field_size ! field_size = 1024 print "field_cell =", field_cell ! field_cell = 0.5 print "taper_size =", taper_size ! taper_size = 20 ![@/opt64/bin/pima_mupet_01.dfm J1217+5835_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 29 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 29 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 29 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=1.333807Jy sigma=7.649583 ! Fit after self-cal, rms=0.231157Jy sigma=1.328695 ! 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=3.048 mas, bmaj=4.411 mas, bpa=-3.043 degrees ! Estimated noise=2.12299 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.638666 Jy ! Component: 100 - total flux cleaned = 0.791397 Jy ! Total flux subtracted in 100 components = 0.791397 Jy ! Clean residual min=-0.027446 max=0.058491 Jy/beam ! Clean residual mean=0.000054 rms=0.007602 Jy/beam ! Combined flux in latest and established models = 0.791397 Jy ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 8 components and 0.791397 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.184316Jy sigma=0.962135 ! Fit after self-cal, rms=0.184287Jy sigma=0.961995 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0603832 Jy ! Component: 100 - total flux cleaned = 0.0910561 Jy ! Total flux subtracted in 100 components = 0.0910561 Jy ! Clean residual min=-0.024613 max=0.031715 Jy/beam ! Clean residual mean=0.000041 rms=0.006001 Jy/beam ! Combined flux in latest and established models = 0.882453 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 18 components and 0.882453 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.161563Jy sigma=0.832330 ! Fit after self-cal, rms=0.161466Jy sigma=0.831745 ! Inverting map print "============== Finished uniform weighting clean ===============" ! ============== Finished uniform weighting clean =============== ! uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.745 mas, bmaj=5.08 mas, bpa=-2.687 degrees ! Estimated noise=1.45099 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.0170044 Jy ! Component: 100 - total flux cleaned = 0.021063 Jy ! Component: 150 - total flux cleaned = 0.0213368 Jy ! Component: 200 - total flux cleaned = 0.0218068 Jy ! Total flux subtracted in 200 components = 0.0218068 Jy ! Clean residual min=-0.024356 max=0.030130 Jy/beam ! Clean residual mean=0.000007 rms=0.005943 Jy/beam ! Combined flux in latest and established models = 0.90426 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 32 components and 0.90426 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.159184Jy sigma=0.821380 ! Fit after self-cal, rms=0.159167Jy sigma=0.821290 ! Inverting map print "============== Finished natural weighting clean ==============" ! ============== Finished natural weighting clean ============== ! restore ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.745 x 5.08 at -2.687 degrees (North through East) ! Clean map min=-0.023746 max=0.83572 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.03 FD 0.98 HN 1.05 KP 0.99 ! LA 0.99 MK 1.03 NL 0.99 OV 1.17 ! PT 0.99 SC 1.01 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.88 FD 0.97 HN 1.00 KP 0.97 ! LA 1.01 MK 1.05 NL 1.03 OV 1.00 ! PT 1.01 SC 0.91 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.04 FD 0.98 HN 1.02 KP 0.99 ! LA 1.00 MK 1.04 NL 1.01 OV 1.00 ! PT 1.01 SC 1.04 ! ! ! Fit before self-cal, rms=0.159167Jy sigma=0.821290 ! Fit after self-cal, rms=0.144685Jy sigma=0.761640 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.773 mas, bmaj=5.096 mas, bpa=-2.379 degrees ! Estimated noise=1.46157 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.00076018 Jy ! Component: 100 - total flux cleaned = 0.000133599 Jy ! Component: 150 - total flux cleaned = -4.35902e-05 Jy ! Component: 200 - total flux cleaned = -0.000194712 Jy ! Total flux subtracted in 200 components = -0.000194712 Jy ! Clean residual min=-0.022687 max=0.027214 Jy/beam ! Clean residual mean=0.000007 rms=0.005694 Jy/beam ! Combined flux in latest and established models = 0.904065 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 39 components and 0.904065 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.144038Jy sigma=0.759586 ! Fit after self-cal, rms=0.144030Jy sigma=0.759498 ! Inverting map selfcal true, true, soltime1 ! Performing amp+phase self-cal over 120 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.144030Jy sigma=0.759498 ! Fit after self-cal, rms=0.080932Jy sigma=0.405338 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.767 mas, bmaj=5.11 mas, bpa=-1.85 degrees ! Estimated noise=1.45915 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00119502 Jy ! Component: 100 - total flux cleaned = 0.00119218 Jy ! Component: 150 - total flux cleaned = 0.00108193 Jy ! Component: 200 - total flux cleaned = 0.000873676 Jy ! Total flux subtracted in 200 components = 0.000873676 Jy ! Clean residual min=-0.005836 max=0.013532 Jy/beam ! Clean residual mean=-0.000005 rms=0.001818 Jy/beam ! Combined flux in latest and established models = 0.904939 Jy selfcal ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 46 components and 0.904939 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.080516Jy sigma=0.403225 ! Fit after self-cal, rms=0.080585Jy sigma=0.402905 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 (10.5, -7.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0132154 Jy ! Component: 100 - total flux cleaned = 0.0210424 Jy ! Component: 150 - total flux cleaned = 0.0257834 Jy ! Component: 200 - total flux cleaned = 0.0275869 Jy ! Total flux subtracted in 200 components = 0.0275869 Jy ! Clean residual min=-0.004525 max=0.008107 Jy/beam ! Clean residual mean=0.000013 rms=0.001343 Jy/beam ! Combined flux in latest and established models = 0.932525 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 70 components and 0.932525 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.078245Jy sigma=0.386113 ! Fit after self-cal, rms=0.077910Jy sigma=0.384408 ! Inverting map ! Added new window around map position (23, -12). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00778861 Jy ! Component: 100 - total flux cleaned = 0.0123418 Jy ! Component: 150 - total flux cleaned = 0.0149924 Jy ! Component: 200 - total flux cleaned = 0.0158337 Jy ! Total flux subtracted in 200 components = 0.0158337 Jy ! Clean residual min=-0.004444 max=0.005189 Jy/beam ! Clean residual mean=0.000018 rms=0.001120 Jy/beam ! Combined flux in latest and established models = 0.948359 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 86 components and 0.948359 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.076959Jy sigma=0.378646 ! Fit after self-cal, rms=0.076910Jy sigma=0.378371 ! Inverting map ! Added new window around map position (-102, -73.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00448557 Jy ! Component: 100 - total flux cleaned = 0.00565402 Jy ! Component: 150 - total flux cleaned = 0.00606505 Jy ! Component: 200 - total flux cleaned = 0.0061595 Jy ! Total flux subtracted in 200 components = 0.0061595 Jy ! Clean residual min=-0.004223 max=0.004787 Jy/beam ! Clean residual mean=0.000008 rms=0.001042 Jy/beam ! Combined flux in latest and established models = 0.954519 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 99 components and 0.954519 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.076515Jy sigma=0.376101 ! Fit after self-cal, rms=0.076501Jy sigma=0.376050 ! Inverting map ! Added new window around map position (30.5, -19). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00455302 Jy ! Component: 100 - total flux cleaned = 0.00513794 Jy ! Component: 150 - total flux cleaned = 0.00495339 Jy ! Component: 200 - total flux cleaned = 0.00477823 Jy ! Total flux subtracted in 200 components = 0.00477823 Jy ! Clean residual min=-0.003737 max=0.004396 Jy/beam ! Clean residual mean=0.000008 rms=0.000967 Jy/beam ! Combined flux in latest and established models = 0.959297 Jy ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 118 components and 0.959297 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.076131Jy sigma=0.373852 ! Fit after self-cal, rms=0.076121Jy sigma=0.373651 ! Inverting map ! Added new window around map position (23, 96.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00330379 Jy ! Component: 100 - total flux cleaned = 0.00415981 Jy ! Component: 150 - total flux cleaned = 0.00434013 Jy ! Component: 200 - total flux cleaned = 0.00459939 Jy ! Total flux subtracted in 200 components = 0.00459939 Jy ! Clean residual min=-0.003869 max=0.003660 Jy/beam ! Clean residual mean=0.000009 rms=0.000896 Jy/beam ! Combined flux in latest and established models = 0.963896 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 128 components and 0.963896 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.075795Jy sigma=0.371734 ! Fit after self-cal, rms=0.075786Jy sigma=0.371638 ! Inverting map selfcal true, true, soltime2 ! Performing amp+phase self-cal over 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.075786Jy sigma=0.371638 ! Fit after self-cal, rms=0.075377Jy sigma=0.369851 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.758 mas, bmaj=5.11 mas, bpa=-1.883 degrees ! Estimated noise=1.45928 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00253215 Jy ! Component: 100 - total flux cleaned = 0.00336032 Jy ! Component: 150 - total flux cleaned = 0.00385916 Jy ! Component: 200 - total flux cleaned = 0.00432052 Jy ! Total flux subtracted in 200 components = 0.00432052 Jy ! Clean residual min=-0.003228 max=0.003495 Jy/beam ! Clean residual mean=0.000006 rms=0.000784 Jy/beam ! Combined flux in latest and established models = 0.968217 Jy selfcal ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 137 components and 0.968217 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.075161Jy sigma=0.368507 ! Fit after self-cal, rms=0.075147Jy sigma=0.368441 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 (10.5, 23). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00099499 Jy ! Component: 100 - total flux cleaned = 0.00145428 Jy ! Component: 150 - total flux cleaned = 0.00177099 Jy ! Component: 200 - total flux cleaned = 0.00176851 Jy ! Total flux subtracted in 200 components = 0.00176851 Jy ! Clean residual min=-0.002898 max=0.003116 Jy/beam ! Clean residual mean=0.000007 rms=0.000730 Jy/beam ! Combined flux in latest and established models = 0.969985 Jy ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 151 components and 0.969985 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.074898Jy sigma=0.367022 ! Fit after self-cal, rms=0.074901Jy sigma=0.366997 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime2 ! Performing amp+phase self-cal over 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.074901Jy sigma=0.366997 ! Fit after self-cal, rms=0.074850Jy sigma=0.366757 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.261 mas, bmaj=9.104 mas, bpa=13.14 degrees ! Estimated noise=1.82753 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00248467 Jy ! Component: 100 - total flux cleaned = 0.00314796 Jy ! Component: 150 - total flux cleaned = 0.00345657 Jy ! Component: 200 - total flux cleaned = 0.00373064 Jy ! Total flux subtracted in 200 components = 0.00373064 Jy ! Clean residual min=-0.002902 max=0.003055 Jy/beam ! Clean residual mean=0.000017 rms=0.000763 Jy/beam ! Combined flux in latest and established models = 0.973716 Jy selfcal ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 171 components and 0.973716 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.074797Jy sigma=0.366369 ! Fit after self-cal, rms=0.074796Jy sigma=0.366314 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.00112656 Jy ! Component: 100 - total flux cleaned = 0.00129074 Jy ! Component: 150 - total flux cleaned = 0.00144561 Jy ! Component: 200 - total flux cleaned = 0.00151768 Jy ! Total flux subtracted in 200 components = 0.00151768 Jy ! Clean residual min=-0.002691 max=0.002947 Jy/beam ! Clean residual mean=0.000018 rms=0.000745 Jy/beam ! Combined flux in latest and established models = 0.975234 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 188 components and 0.975234 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.074778Jy sigma=0.366195 ! Fit after self-cal, rms=0.074781Jy sigma=0.366188 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "============== Finished amplitude self-cal ==============" ! ============== Finished amplitude self-cal ============== restore ! Inverting map and beam ! Estimated beam: bmin=3.754 mas, bmaj=5.111 mas, bpa=-1.864 degrees ! Estimated noise=1.46018 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.754 x 5.111 at -1.864 degrees (North through East) ! Clean map min=-0.0033813 max=0.84176 Jy/beam if(peak(flux) > thresh) selfcal true, true clean clean_niter,clean_gain selfcal flux_cutoff = imstat(rms) * dynam repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) selfcal true, true end if ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.074781Jy sigma=0.366188 ! Fit after self-cal, rms=0.068485Jy sigma=0.326698 ! Inverting map and beam ! Estimated beam: bmin=3.752 mas, bmaj=5.108 mas, bpa=-1.901 degrees ! Estimated noise=1.46106 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000865518 Jy ! Component: 100 - total flux cleaned = 0.000650946 Jy ! Component: 150 - total flux cleaned = 0.000652516 Jy ! Component: 200 - total flux cleaned = 0.000652822 Jy ! Total flux subtracted in 200 components = 0.000652822 Jy ! Clean residual min=-0.002467 max=0.002824 Jy/beam ! Clean residual mean=0.000008 rms=0.000644 Jy/beam ! Combined flux in latest and established models = 0.975887 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 198 components and 0.975887 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.068354Jy sigma=0.325868 ! Fit after self-cal, rms=0.068348Jy sigma=0.325844 ! Inverting map ! Added new window around map position (11, -20). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00255762 Jy ! Component: 100 - total flux cleaned = 0.00236434 Jy ! Component: 150 - total flux cleaned = 0.00224447 Jy ! Component: 200 - total flux cleaned = 0.00207398 Jy ! Total flux subtracted in 200 components = 0.00207398 Jy ! Clean residual min=-0.002334 max=0.002421 Jy/beam ! Clean residual mean=0.000010 rms=0.000605 Jy/beam ! Combined flux in latest and established models = 0.977961 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 208 components and 0.97796 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.068188Jy sigma=0.324923 ! Fit after self-cal, rms=0.068183Jy sigma=0.324894 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.068183Jy sigma=0.324894 ! Fit after self-cal, rms=0.068228Jy sigma=0.324773 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=3.061 mas, bmaj=4.433 mas, bpa=-2.146 degrees ! Estimated noise=2.14562 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.645418 Jy ! Component: 100 - total flux cleaned = 0.800194 Jy ! Component: 150 - total flux cleaned = 0.86029 Jy ! Component: 200 - total flux cleaned = 0.88924 Jy ! Total flux subtracted in 200 components = 0.88924 Jy ! Clean residual min=-0.009939 max=0.015346 Jy/beam ! Clean residual mean=0.000034 rms=0.002532 Jy/beam ! Combined flux in latest and established models = 0.88924 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 22 model components to the UV plane model. ! The established model now contains 22 components and 0.88924 Jy ! Inverting map and beam ! Estimated beam: bmin=3.749 mas, bmaj=5.109 mas, bpa=-1.876 degrees ! Estimated noise=1.46205 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.0196461 Jy ! Component: 100 - total flux cleaned = 0.033433 Jy ! Component: 150 - total flux cleaned = 0.0441449 Jy ! Component: 200 - total flux cleaned = 0.0529747 Jy ! Total flux subtracted in 200 components = 0.0529747 Jy ! Clean residual min=-0.004715 max=0.006680 Jy/beam ! Clean residual mean=0.000002 rms=0.001266 Jy/beam ! Combined flux in latest and established models = 0.942215 Jy ! Adding 36 model components to the UV plane model. ! The established model now contains 50 components and 0.942215 Jy ! Inverting map ! Added new window around map position (17, -15.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00768225 Jy ! Component: 100 - total flux cleaned = 0.0141204 Jy ! Component: 150 - total flux cleaned = 0.0197228 Jy ! Component: 200 - total flux cleaned = 0.0244412 Jy ! Total flux subtracted in 200 components = 0.0244412 Jy ! Clean residual min=-0.003723 max=0.004077 Jy/beam ! Clean residual mean=0.000006 rms=0.000948 Jy/beam ! Combined flux in latest and established models = 0.966656 Jy ! Adding 38 model components to the UV plane model. ! The established model now contains 72 components and 0.966656 Jy ! Inverting map ! Added new window around map position (4.5, 0.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00263553 Jy ! Component: 100 - total flux cleaned = 0.00506222 Jy ! Component: 150 - total flux cleaned = 0.00697838 Jy ! Component: 200 - total flux cleaned = 0.008744 Jy ! Total flux subtracted in 200 components = 0.008744 Jy ! Clean residual min=-0.003450 max=0.003174 Jy/beam ! Clean residual mean=0.000008 rms=0.000819 Jy/beam ! Combined flux in latest and established models = 0.9754 Jy ! Adding 45 model components to the UV plane model. ! The established model now contains 93 components and 0.9754 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.254 mas, bmaj=9.1 mas, bpa=13.03 degrees ! Estimated noise=1.83221 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00284374 Jy ! Component: 100 - total flux cleaned = 0.0037586 Jy ! Component: 150 - total flux cleaned = 0.00445825 Jy ! Component: 200 - total flux cleaned = 0.00507763 Jy ! Total flux subtracted in 200 components = 0.00507763 Jy ! Clean residual min=-0.002546 max=0.002491 Jy/beam ! Clean residual mean=0.000018 rms=0.000697 Jy/beam ! Combined flux in latest and established models = 0.980478 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.000388251 Jy ! Component: 100 - total flux cleaned = 0.000540228 Jy ! Component: 150 - total flux cleaned = 0.000540842 Jy ! Component: 200 - total flux cleaned = 0.000541445 Jy ! Total flux subtracted in 200 components = 0.000541445 Jy ! Clean residual min=-0.002536 max=0.002406 Jy/beam ! Clean residual mean=0.000020 rms=0.000669 Jy/beam ! Combined flux in latest and established models = 0.981019 Jy ! Adding 59 model components to the UV plane model. ! The established model now contains 149 components and 0.981019 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.749 mas, bmaj=5.109 mas, bpa=-1.876 degrees ! Estimated noise=1.46205 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000708168 Jy ! Component: 100 - total flux cleaned = -0.000382985 Jy ! Component: 150 - total flux cleaned = -7.71368e-05 Jy ! Component: 200 - total flux cleaned = 2.1613e-05 Jy ! Total flux subtracted in 200 components = 2.1613e-05 Jy ! Clean residual min=-0.003039 max=0.002747 Jy/beam ! Clean residual mean=0.000010 rms=0.000707 Jy/beam ! Combined flux in latest and established models = 0.981041 Jy selfcal ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 166 components and 0.981041 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.068823Jy sigma=0.328734 ! Fit after self-cal, rms=0.068772Jy sigma=0.328419 wmodel J1217+5835_S_map.mod ! Writing 166 model components to file: J1217+5835_S_map.mod wobs J1217+5835_S_uvs.fits ! Writing UV FITS file: J1217+5835_S_uvs.fits wwins J1217+5835_S_map.win ! wwins: Wrote 10 windows to J1217+5835_S_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -7.10775e-05 Jy ! Component: 100 - total flux cleaned = 0.00012945 Jy ! Component: 150 - total flux cleaned = 0.000321909 Jy ! Component: 200 - total flux cleaned = 0.000461537 Jy ! Component: 250 - total flux cleaned = 0.00064268 Jy ! Component: 300 - total flux cleaned = 0.000820055 Jy ! Component: 350 - total flux cleaned = 0.00108045 Jy ! Component: 400 - total flux cleaned = 0.00112219 Jy ! Component: 450 - total flux cleaned = 0.00120579 Jy ! Component: 500 - total flux cleaned = 0.00132865 Jy ! Component: 550 - total flux cleaned = 0.00148894 Jy ! Component: 600 - total flux cleaned = 0.00156845 Jy ! Component: 650 - total flux cleaned = 0.00172322 Jy ! Component: 700 - total flux cleaned = 0.00179859 Jy ! Component: 750 - total flux cleaned = 0.00176059 Jy ! Component: 800 - total flux cleaned = 0.00179772 Jy ! Component: 850 - total flux cleaned = 0.00190761 Jy ! Component: 900 - total flux cleaned = 0.00194322 Jy ! Component: 950 - total flux cleaned = 0.00194281 Jy ! Component: 1000 - total flux cleaned = 0.0019431 Jy ! Component: 1050 - total flux cleaned = 0.00201233 Jy ! Component: 1100 - total flux cleaned = 0.00204628 Jy ! Component: 1150 - total flux cleaned = 0.00204652 Jy ! Component: 1200 - total flux cleaned = 0.00208048 Jy ! Component: 1250 - total flux cleaned = 0.00214628 Jy ! Component: 1300 - total flux cleaned = 0.00214647 Jy ! Component: 1350 - total flux cleaned = 0.00201961 Jy ! Component: 1400 - total flux cleaned = 0.00217727 Jy ! Component: 1450 - total flux cleaned = 0.00214632 Jy ! Component: 1500 - total flux cleaned = 0.00220729 Jy ! Component: 1550 - total flux cleaned = 0.00220705 Jy ! Component: 1600 - total flux cleaned = 0.00220729 Jy ! Component: 1650 - total flux cleaned = 0.00214825 Jy ! Component: 1700 - total flux cleaned = 0.00232472 Jy ! Component: 1750 - total flux cleaned = 0.00223772 Jy ! Component: 1800 - total flux cleaned = 0.0022668 Jy ! Component: 1850 - total flux cleaned = 0.00226609 Jy ! Component: 1900 - total flux cleaned = 0.00220962 Jy ! Component: 1950 - total flux cleaned = 0.00220969 Jy ! Component: 2000 - total flux cleaned = 0.00223701 Jy ! Component: 2050 - total flux cleaned = 0.00226405 Jy ! Component: 2100 - total flux cleaned = 0.00226358 Jy ! Component: 2150 - total flux cleaned = 0.00220958 Jy ! Component: 2200 - total flux cleaned = 0.00223614 Jy ! Component: 2250 - total flux cleaned = 0.00218344 Jy ! Component: 2300 - total flux cleaned = 0.00220971 Jy ! Component: 2350 - total flux cleaned = 0.00218381 Jy ! Component: 2400 - total flux cleaned = 0.00223474 Jy ! Component: 2450 - total flux cleaned = 0.00215852 Jy ! Component: 2500 - total flux cleaned = 0.00213319 Jy ! Component: 2550 - total flux cleaned = 0.00213321 Jy ! Component: 2600 - total flux cleaned = 0.00213329 Jy ! Component: 2650 - total flux cleaned = 0.00218227 Jy ! Component: 2700 - total flux cleaned = 0.00218242 Jy ! Component: 2750 - total flux cleaned = 0.0022304 Jy ! Component: 2800 - total flux cleaned = 0.00225434 Jy ! Component: 2850 - total flux cleaned = 0.0023018 Jy ! Component: 2900 - total flux cleaned = 0.00227811 Jy ! Component: 2950 - total flux cleaned = 0.00220807 Jy ! Component: 3000 - total flux cleaned = 0.00225474 Jy ! Component: 3050 - total flux cleaned = 0.00234695 Jy ! Component: 3100 - total flux cleaned = 0.0023238 Jy ! Component: 3150 - total flux cleaned = 0.00234646 Jy ! Component: 3200 - total flux cleaned = 0.00241396 Jy ! Component: 3250 - total flux cleaned = 0.00227987 Jy ! Component: 3300 - total flux cleaned = 0.00239103 Jy ! Component: 3350 - total flux cleaned = 0.0024789 Jy ! Component: 3400 - total flux cleaned = 0.00245699 Jy ! Component: 3450 - total flux cleaned = 0.00250051 Jy ! Component: 3500 - total flux cleaned = 0.00258673 Jy ! Component: 3550 - total flux cleaned = 0.00254409 Jy ! Component: 3600 - total flux cleaned = 0.00254411 Jy ! Component: 3650 - total flux cleaned = 0.00254397 Jy ! Component: 3700 - total flux cleaned = 0.00254391 Jy ! Component: 3750 - total flux cleaned = 0.00246046 Jy ! Component: 3800 - total flux cleaned = 0.00248109 Jy ! Component: 3850 - total flux cleaned = 0.00244 Jy ! Component: 3900 - total flux cleaned = 0.00244003 Jy ! Component: 3950 - total flux cleaned = 0.00241981 Jy ! Component: 4000 - total flux cleaned = 0.00233927 Jy ! Component: 4050 - total flux cleaned = 0.00239915 Jy ! Total flux subtracted in 4096 components = 0.00231942 Jy ! Clean residual min=-0.001185 max=0.001045 Jy/beam ! Clean residual mean=0.000008 rms=0.000347 Jy/beam ! Combined flux in latest and established models = 0.98336 Jy keep ! Adding 292 model components to the UV plane model. ! The established model now contains 458 components and 0.983361 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /image_prod/ug002_uva//J1217+5835_S_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.749 x 5.109 at -1.876 degrees (North through East) ! Clean map min=-0.0023345 max=0.84149 Jy/beam ! Writing clean map to FITS file: /image_prod/ug002_uva//J1217+5835_S_map.fits wmodel /image_prod/ug002_uva//J1217+5835_S_map.mod ! Writing 458 model components to file: /image_prod/ug002_uva//J1217+5835_S_map.mod wobs /image_prod/ug002_uva//J1217+5835_S_uvs.fits ! Writing UV FITS file: /image_prod/ug002_uva//J1217+5835_S_uvs.fits wwins /image_prod/ug002_uva//J1217+5835_S_map.win ! wwins: Wrote 11 windows to /image_prod/ug002_uva//J1217+5835_S_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.841491 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000345432 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.12315 0.12315 0.2463 0.4926 0.985199 1.9704 3.9408 7.8816 15.7632 31.5264 63.0528 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.001 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.841 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 2436.05 quit ! Quitting program ! Log file /image_prod/ug002_uva//J1217+5835_S_dfm.log closed on Thu Mar 5 00:33:49 2020