! Started logfile: /scr/pima/rv125_uvs/J0755+2542_S_dfm.log on Sun Dec 23 06:57:02 2018 float field_size field_size = 1024 float field_cell field_cell = 0.5000 float taper_size taper_size = 20.0 float freq freq = 2.23275 obs /scr/pima/rv125_uvs/J0755+2542_S_uva.fits ! Reading UV FITS file: /scr/pima/rv125_uvs/J0755+2542_S_uva.fits ! AN table 1: 53 integrations on 28 of 28 possible baselines. ! Apparent sampling: 0.708895 visibilities/baseline/integration-bin. ! Found source: J0755+2542 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 2.23275e+09 8e+06 1 8e+06 ! 02 2 2.26275e+09 8e+06 1 8e+06 ! 03 3 2.35275e+09 8e+06 1 8e+06 ! 04 4 2.37275e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 69 lines of history. ! ! Reading 4208 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 J0755+2542_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 72 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 72 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 424 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 72 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=1.002769Jy sigma=20.477668 ! Fit after self-cal, rms=0.971895Jy sigma=19.866005 ! 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.333 mas, bmaj=7.046 mas, bpa=27.12 degrees ! Estimated noise=1.00692 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.0237445 Jy ! Component: 100 - total flux cleaned = 0.0283156 Jy ! Total flux subtracted in 100 components = 0.0283156 Jy ! Clean residual min=-0.003746 max=0.004429 Jy/beam ! Clean residual mean=0.000002 rms=0.001091 Jy/beam ! Combined flux in latest and established models = 0.0283156 Jy ! Performing phase self-cal ! Adding 4 model components to the UV plane model. ! The established model now contains 4 components and 0.0283156 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.047973Jy sigma=0.943046 ! Fit after self-cal, rms=0.046128Jy sigma=0.905211 ! 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=4.045 mas, bmaj=7.731 mas, bpa=24.32 degrees ! Estimated noise=0.628057 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.00266597 Jy ! Component: 100 - total flux cleaned = 0.00309376 Jy ! Component: 150 - total flux cleaned = 0.00312809 Jy ! Component: 200 - total flux cleaned = 0.00312922 Jy ! Total flux subtracted in 200 components = 0.00312922 Jy ! Clean residual min=-0.002556 max=0.002735 Jy/beam ! Clean residual mean=-0.000003 rms=0.000700 Jy/beam ! Combined flux in latest and established models = 0.0314449 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 18 components and 0.0314449 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.045940Jy sigma=0.902565 ! Fit after self-cal, rms=0.045928Jy sigma=0.902312 ! Inverting map print "============== Finished natural weighting clean ==============" ! ============== Finished natural weighting clean ============== ! restore ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 4.045 x 7.731 at 24.32 degrees (North through East) ! Clean map min=-0.0024678 max=0.03146 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.43 HN 1.01 KP 0.86 LA 0.84 ! NL 1.30 OV 0.88 PT 0.98 SC 1.13 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 1.18 HN 0.88 KP 0.99 LA 0.83 ! NL 1.14 OV 1.03 PT 1.20 SC 0.91 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* HN 0.00* KP 0.00* LA 0.00* ! NL 0.00* OV 0.00* PT 0.00* SC 0.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 HN 0.86 KP 0.98 LA 0.96 ! NL 1.26 OV 0.83 PT 1.04 SC 0.90 ! ! ! Fit before self-cal, rms=0.045928Jy sigma=0.902312 ! Fit after self-cal, rms=0.045862Jy sigma=0.893727 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.991 mas, bmaj=7.669 mas, bpa=23.06 degrees ! Estimated noise=0.623473 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 = -6.36043e-05 Jy ! Component: 100 - total flux cleaned = -0.000144202 Jy ! Component: 150 - total flux cleaned = -0.000142148 Jy ! Component: 200 - total flux cleaned = -0.000107757 Jy ! Total flux subtracted in 200 components = -0.000107757 Jy ! Clean residual min=-0.002039 max=0.002017 Jy/beam ! Clean residual mean=-0.000002 rms=0.000511 Jy/beam ! Combined flux in latest and established models = 0.0313371 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 30 components and 0.0313371 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.045831Jy sigma=0.893311 ! Fit after self-cal, rms=0.046047Jy sigma=0.891314 ! 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.046047Jy sigma=0.891314 ! Fit after self-cal, rms=0.047483Jy sigma=0.889525 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=4.01 mas, bmaj=7.596 mas, bpa=22.83 degrees ! Estimated noise=0.628969 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000104134 Jy ! Component: 100 - total flux cleaned = 0.000102212 Jy ! Component: 150 - total flux cleaned = 0.000121176 Jy ! Component: 200 - total flux cleaned = 0.000171868 Jy ! Total flux subtracted in 200 components = 0.000171868 Jy ! Clean residual min=-0.002040 max=0.001897 Jy/beam ! Clean residual mean=-0.000003 rms=0.000494 Jy/beam ! Combined flux in latest and established models = 0.031509 Jy selfcal ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 39 components and 0.031509 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.047462Jy sigma=0.889252 ! Fit after self-cal, rms=0.047568Jy sigma=0.888484 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.000159369 Jy ! Component: 100 - total flux cleaned = -0.000192556 Jy ! Component: 150 - total flux cleaned = -0.000207587 Jy ! Component: 200 - total flux cleaned = -0.000220693 Jy ! Total flux subtracted in 200 components = -0.000220693 Jy ! Clean residual min=-0.001935 max=0.001810 Jy/beam ! Clean residual mean=-0.000003 rms=0.000490 Jy/beam ! Combined flux in latest and established models = 0.0312883 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 44 components and 0.0312883 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.047553Jy sigma=0.888343 ! Fit after self-cal, rms=0.047562Jy sigma=0.888312 ! 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.047562Jy sigma=0.888312 ! Fit after self-cal, rms=0.048166Jy sigma=0.888150 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=4.009 mas, bmaj=7.589 mas, bpa=22.97 degrees ! Estimated noise=0.632388 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 3.65404e-06 Jy ! Component: 100 - total flux cleaned = 2.65855e-06 Jy ! Component: 150 - total flux cleaned = -1.16184e-05 Jy ! Component: 200 - total flux cleaned = -1.16824e-05 Jy ! Total flux subtracted in 200 components = -1.16824e-05 Jy ! Clean residual min=-0.001960 max=0.001737 Jy/beam ! Clean residual mean=-0.000003 rms=0.000488 Jy/beam ! Combined flux in latest and established models = 0.0312766 Jy selfcal ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 45 components and 0.0312766 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.048153Jy sigma=0.888011 ! Fit after self-cal, rms=0.048183Jy sigma=0.887754 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.00012224 Jy ! Component: 100 - total flux cleaned = -0.00012249 Jy ! Component: 150 - total flux cleaned = -0.000122427 Jy ! Component: 200 - total flux cleaned = -0.000111148 Jy ! Total flux subtracted in 200 components = -0.000111148 Jy ! Clean residual min=-0.001951 max=0.001694 Jy/beam ! Clean residual mean=-0.000003 rms=0.000491 Jy/beam ! Combined flux in latest and established models = 0.0311654 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 47 components and 0.0311654 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.048174Jy sigma=0.887663 ! Fit after self-cal, rms=0.048174Jy sigma=0.887645 ! 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.048174Jy sigma=0.887645 ! Fit after self-cal, rms=0.048581Jy sigma=0.887592 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.003 mas, bmaj=10.28 mas, bpa=-16.08 degrees ! Estimated noise=0.767006 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 4.25209e-05 Jy ! Component: 100 - total flux cleaned = 4.30739e-05 Jy ! Component: 150 - total flux cleaned = 2.93882e-05 Jy ! Component: 200 - total flux cleaned = 2.93732e-05 Jy ! Total flux subtracted in 200 components = 2.93732e-05 Jy ! Clean residual min=-0.001900 max=0.002005 Jy/beam ! Clean residual mean=-0.000004 rms=0.000558 Jy/beam ! Combined flux in latest and established models = 0.0311948 Jy selfcal ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 55 components and 0.0311948 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.048579Jy sigma=0.887597 ! Fit after self-cal, rms=0.048595Jy sigma=0.887565 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 = 2.74983e-05 Jy ! Component: 100 - total flux cleaned = 1.57062e-05 Jy ! Component: 150 - total flux cleaned = 1.01215e-05 Jy ! Component: 200 - total flux cleaned = 5.41143e-07 Jy ! Total flux subtracted in 200 components = 5.41143e-07 Jy ! Clean residual min=-0.001888 max=0.002013 Jy/beam ! Clean residual mean=-0.000004 rms=0.000557 Jy/beam ! Combined flux in latest and established models = 0.0311954 Jy ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 57 components and 0.0311954 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.048593Jy sigma=0.887587 ! Fit after self-cal, rms=0.048592Jy sigma=0.887580 ! 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=4.007 mas, bmaj=7.595 mas, bpa=23.15 degrees ! Estimated noise=0.635871 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 4.007 x 7.595 at 23.15 degrees (North through East) ! Clean map min=-0.0019006 max=0.03254 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 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.326 mas, bmaj=6.723 mas, bpa=25.62 degrees ! Estimated noise=0.998273 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0257312 Jy ! Component: 100 - total flux cleaned = 0.0316523 Jy ! Component: 150 - total flux cleaned = 0.0324911 Jy ! Component: 200 - total flux cleaned = 0.032422 Jy ! Total flux subtracted in 200 components = 0.032422 Jy ! Clean residual min=-0.003179 max=0.003356 Jy/beam ! Clean residual mean=-0.000000 rms=0.000859 Jy/beam ! Combined flux in latest and established models = 0.032422 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.032422 Jy ! Inverting map and beam ! Estimated beam: bmin=4.007 mas, bmaj=7.595 mas, bpa=23.15 degrees ! Estimated noise=0.635871 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.00107908 Jy ! Component: 100 - total flux cleaned = -0.00121754 Jy ! Component: 150 - total flux cleaned = -0.00140615 Jy ! Component: 200 - total flux cleaned = -0.001452 Jy ! Total flux subtracted in 200 components = -0.001452 Jy ! Clean residual min=-0.002061 max=0.001723 Jy/beam ! Clean residual mean=-0.000003 rms=0.000502 Jy/beam ! Combined flux in latest and established models = 0.03097 Jy ! Adding 13 model components to the UV plane model. ! The established model now contains 34 components and 0.03097 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.003 mas, bmaj=10.28 mas, bpa=-16.08 degrees ! Estimated noise=0.767006 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000117419 Jy ! Component: 100 - total flux cleaned = 0.000127537 Jy ! Component: 150 - total flux cleaned = 0.00012739 Jy ! Component: 200 - total flux cleaned = 0.000127109 Jy ! Total flux subtracted in 200 components = 0.000127109 Jy ! Clean residual min=-0.001935 max=0.002037 Jy/beam ! Clean residual mean=-0.000004 rms=0.000560 Jy/beam ! Combined flux in latest and established models = 0.0310971 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 = 5.4266e-06 Jy ! Component: 100 - total flux cleaned = 5.44971e-06 Jy ! Component: 150 - total flux cleaned = 1.05808e-05 Jy ! Component: 200 - total flux cleaned = 1.05478e-05 Jy ! Total flux subtracted in 200 components = 1.05478e-05 Jy ! Clean residual min=-0.001931 max=0.002037 Jy/beam ! Clean residual mean=-0.000004 rms=0.000559 Jy/beam ! Combined flux in latest and established models = 0.0311076 Jy ! Adding 16 model components to the UV plane model. ! The established model now contains 45 components and 0.0311076 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=4.007 mas, bmaj=7.595 mas, bpa=23.15 degrees ! Estimated noise=0.635871 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000173246 Jy ! Component: 100 - total flux cleaned = -0.000190942 Jy ! Component: 150 - total flux cleaned = -0.000191563 Jy ! Component: 200 - total flux cleaned = -0.000206477 Jy ! Total flux subtracted in 200 components = -0.000206477 Jy ! Clean residual min=-0.001975 max=0.001714 Jy/beam ! Clean residual mean=-0.000003 rms=0.000497 Jy/beam ! Combined flux in latest and established models = 0.0309012 Jy selfcal ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 52 components and 0.0309012 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.048605Jy sigma=0.887687 ! Fit after self-cal, rms=0.048604Jy sigma=0.887654 wmodel J0755+2542_S_map.mod ! Writing 52 model components to file: J0755+2542_S_map.mod wobs J0755+2542_S_uvs.fits ! Writing UV FITS file: J0755+2542_S_uvs.fits wwins J0755+2542_S_map.win ! wwins: Wrote 1 windows to J0755+2542_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.000907632 Jy ! Component: 100 - total flux cleaned = -0.00149262 Jy ! Component: 150 - total flux cleaned = -0.00194967 Jy ! Component: 200 - total flux cleaned = -0.00226913 Jy ! Component: 250 - total flux cleaned = -0.00243641 Jy ! Component: 300 - total flux cleaned = -0.00264991 Jy ! Component: 350 - total flux cleaned = -0.00270119 Jy ! Component: 400 - total flux cleaned = -0.0028003 Jy ! Component: 450 - total flux cleaned = -0.00280046 Jy ! Component: 500 - total flux cleaned = -0.00284802 Jy ! Component: 550 - total flux cleaned = -0.00280233 Jy ! Component: 600 - total flux cleaned = -0.00282486 Jy ! Component: 650 - total flux cleaned = -0.00284663 Jy ! Component: 700 - total flux cleaned = -0.00288901 Jy ! Component: 750 - total flux cleaned = -0.00284687 Jy ! Component: 800 - total flux cleaned = -0.00288785 Jy ! Component: 850 - total flux cleaned = -0.00296785 Jy ! Component: 900 - total flux cleaned = -0.00290867 Jy ! Component: 950 - total flux cleaned = -0.00292791 Jy ! Component: 1000 - total flux cleaned = -0.00296559 Jy ! Component: 1050 - total flux cleaned = -0.00292834 Jy ! Component: 1100 - total flux cleaned = -0.00292851 Jy ! Component: 1150 - total flux cleaned = -0.00289265 Jy ! Component: 1200 - total flux cleaned = -0.00285737 Jy ! Component: 1250 - total flux cleaned = -0.00280528 Jy ! Component: 1300 - total flux cleaned = -0.00275401 Jy ! Component: 1350 - total flux cleaned = -0.00273711 Jy ! Component: 1400 - total flux cleaned = -0.00273691 Jy ! Component: 1450 - total flux cleaned = -0.0027198 Jy ! Component: 1500 - total flux cleaned = -0.00270361 Jy ! Component: 1550 - total flux cleaned = -0.00270374 Jy ! Component: 1600 - total flux cleaned = -0.0027512 Jy ! Component: 1650 - total flux cleaned = -0.0027199 Jy ! Component: 1700 - total flux cleaned = -0.00278185 Jy ! Component: 1750 - total flux cleaned = -0.0028122 Jy ! Component: 1800 - total flux cleaned = -0.00279716 Jy ! Component: 1850 - total flux cleaned = -0.00284213 Jy ! Component: 1900 - total flux cleaned = -0.00288652 Jy ! Component: 1950 - total flux cleaned = -0.00293045 Jy ! Component: 2000 - total flux cleaned = -0.00294483 Jy ! Component: 2050 - total flux cleaned = -0.00297352 Jy ! Component: 2100 - total flux cleaned = -0.00303013 Jy ! Component: 2150 - total flux cleaned = -0.00301605 Jy ! Component: 2200 - total flux cleaned = -0.00309963 Jy ! Component: 2250 - total flux cleaned = -0.00311329 Jy ! Component: 2300 - total flux cleaned = -0.00314067 Jy ! Component: 2350 - total flux cleaned = -0.00320844 Jy ! Component: 2400 - total flux cleaned = -0.0032353 Jy ! Component: 2450 - total flux cleaned = -0.00330187 Jy ! Component: 2500 - total flux cleaned = -0.00332839 Jy ! Component: 2550 - total flux cleaned = -0.00338063 Jy ! Component: 2600 - total flux cleaned = -0.00340653 Jy ! Component: 2650 - total flux cleaned = -0.00339361 Jy ! Component: 2700 - total flux cleaned = -0.00343219 Jy ! Component: 2750 - total flux cleaned = -0.00349576 Jy ! Component: 2800 - total flux cleaned = -0.00355885 Jy ! Component: 2850 - total flux cleaned = -0.0036463 Jy ! Component: 2900 - total flux cleaned = -0.00368338 Jy ! Component: 2950 - total flux cleaned = -0.003708 Jy ! Component: 3000 - total flux cleaned = -0.00376925 Jy ! Component: 3050 - total flux cleaned = -0.0038056 Jy ! Component: 3100 - total flux cleaned = -0.00384178 Jy ! Component: 3150 - total flux cleaned = -0.00385375 Jy ! Component: 3200 - total flux cleaned = -0.00387761 Jy ! Component: 3250 - total flux cleaned = -0.00388946 Jy ! Component: 3300 - total flux cleaned = -0.00392474 Jy ! Component: 3350 - total flux cleaned = -0.00394808 Jy ! Component: 3400 - total flux cleaned = -0.00399465 Jy ! Component: 3450 - total flux cleaned = -0.00402947 Jy ! Component: 3500 - total flux cleaned = -0.0040063 Jy ! Component: 3550 - total flux cleaned = -0.00405233 Jy ! Component: 3600 - total flux cleaned = -0.00408673 Jy ! Component: 3650 - total flux cleaned = -0.00414379 Jy ! Component: 3700 - total flux cleaned = -0.00416651 Jy ! Component: 3750 - total flux cleaned = -0.00420048 Jy ! Component: 3800 - total flux cleaned = -0.00422303 Jy ! Component: 3850 - total flux cleaned = -0.00424544 Jy ! Component: 3900 - total flux cleaned = -0.00429018 Jy ! Component: 3950 - total flux cleaned = -0.00426794 Jy ! Component: 4000 - total flux cleaned = -0.00432343 Jy ! Component: 4050 - total flux cleaned = -0.00434552 Jy ! Total flux subtracted in 4096 components = -0.00435653 Jy ! Clean residual min=-0.000692 max=0.000752 Jy/beam ! Clean residual mean=-0.000004 rms=0.000194 Jy/beam ! Combined flux in latest and established models = 0.0265446 Jy keep ! Adding 267 model components to the UV plane model. ! The established model now contains 319 components and 0.0265446 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /scr/pima/rv125_uvs/J0755+2542_S_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 4.007 x 7.595 at 23.15 degrees (North through East) ! Clean map min=-0.0012967 max=0.032357 Jy/beam ! Writing clean map to FITS file: /scr/pima/rv125_uvs/J0755+2542_S_map.fits wmodel /scr/pima/rv125_uvs/J0755+2542_S_map.mod ! Writing 319 model components to file: /scr/pima/rv125_uvs/J0755+2542_S_map.mod wobs /scr/pima/rv125_uvs/J0755+2542_S_uvs.fits ! Writing UV FITS file: /scr/pima/rv125_uvs/J0755+2542_S_uvs.fits wwins /scr/pima/rv125_uvs/J0755+2542_S_map.win ! wwins: Wrote 2 windows to /scr/pima/rv125_uvs/J0755+2542_S_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.0323567 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000190784 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -1.76888 1.76888 3.53776 7.07552 14.151 28.3021 56.6042 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0005 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.032 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 169.598 quit ! Quitting program ! Log file /scr/pima/rv125_uvs/J0755+2542_S_dfm.log closed on Sun Dec 23 06:57:06 2018