! Started logfile: /scr/pima/rv121_uvs/J0405-1308_S_dfm.log on Sat Jul 28 00:08:29 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/rv121_uvs/J0405-1308_S_uva.fits ! Reading UV FITS file: /scr/pima/rv121_uvs/J0405-1308_S_uva.fits ! AN table 1: 52 integrations on 3 of 3 possible baselines. ! AN table 2: 16 integrations on 3 of 3 possible baselines. ! AN table 3: 11 integrations on 28 of 28 possible baselines. ! AN table 4: 6 integrations on 45 of 45 possible baselines. ! AN table 5: 25 integrations on 45 of 45 possible baselines. ! Apparent sampling: 0.611956 visibilities/baseline/integration-bin. ! Found source: J0405-1308 ! ! 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 74 lines of history. ! ! Reading 4668 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 J0405-1308_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 54 telescope corrections were flagged in sub-array 1. ! A total of 18 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 16 telescope corrections were flagged in sub-array 4. ! A total of 114 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 54 telescope corrections were flagged in sub-array 1. ! A total of 18 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 16 telescope corrections were flagged in sub-array 4. ! A total of 114 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 54 telescope corrections were flagged in sub-array 1. ! A total of 48 telescope corrections were flagged in sub-array 2. ! A total of 38 telescope corrections were flagged in sub-array 3. ! A total of 34 telescope corrections were flagged in sub-array 4. ! A total of 147 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 54 telescope corrections were flagged in sub-array 1. ! A total of 18 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 16 telescope corrections were flagged in sub-array 4. ! A total of 114 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=1.057107Jy sigma=14.296324 ! Fit after self-cal, rms=0.661772Jy sigma=8.721094 ! 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.384 mas, bmaj=4.238 mas, bpa=20.74 degrees ! Estimated noise=1.25045 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.28133 Jy ! Component: 100 - total flux cleaned = 0.357337 Jy ! Total flux subtracted in 100 components = 0.357337 Jy ! Clean residual min=-0.010835 max=0.025115 Jy/beam ! Clean residual mean=0.000122 rms=0.002701 Jy/beam ! Combined flux in latest and established models = 0.357337 Jy ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 8 components and 0.357337 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.137429Jy sigma=1.722063 ! Fit after self-cal, rms=0.136220Jy sigma=1.721423 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0224186 Jy ! Component: 100 - total flux cleaned = 0.0315975 Jy ! Total flux subtracted in 100 components = 0.0315975 Jy ! Clean residual min=-0.009655 max=0.015037 Jy/beam ! Clean residual mean=0.000091 rms=0.002410 Jy/beam ! Combined flux in latest and established models = 0.388935 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 19 components and 0.388935 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.131589Jy sigma=1.642020 ! Fit after self-cal, rms=0.131423Jy sigma=1.641945 ! Inverting map ! Added new window around map position (6, -17). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.012815 Jy ! Component: 100 - total flux cleaned = 0.0185749 Jy ! Total flux subtracted in 100 components = 0.0185749 Jy ! Clean residual min=-0.009937 max=0.009548 Jy/beam ! Clean residual mean=0.000080 rms=0.002296 Jy/beam ! Combined flux in latest and established models = 0.40751 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 35 components and 0.40751 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.129614Jy sigma=1.613280 ! Fit after self-cal, rms=0.129421Jy sigma=1.610545 ! 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.042 mas, bmaj=4.578 mas, bpa=24.81 degrees ! Estimated noise=0.947866 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.0117111 Jy ! Component: 100 - total flux cleaned = 0.0186985 Jy ! Component: 150 - total flux cleaned = 0.0232736 Jy ! Component: 200 - total flux cleaned = 0.0263247 Jy ! Total flux subtracted in 200 components = 0.0263247 Jy ! Clean residual min=-0.008982 max=0.010456 Jy/beam ! Clean residual mean=0.000033 rms=0.002838 Jy/beam ! Combined flux in latest and established models = 0.433834 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 50 components and 0.433834 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.128289Jy sigma=1.593188 ! Fit after self-cal, rms=0.128129Jy sigma=1.591285 ! 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.042 x 4.578 at 24.81 degrees (North through East) ! Clean map min=-0.0089471 max=0.34594 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 102 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* ON 0.00* WZ 0.00* ! ! A total of 30 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! HH 0.00* ON 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 0.98 HN 1.00 KP 1.03 LA 1.12 ! NL 0.98 ON 1.06 PT 0.97 SC 1.02 ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.04 FD 1.03 HH 1.08 HN 1.06 ! KP 1.00 LA 1.05 NL 1.04 OV 0.98 ! PT 1.02 SC 0.98 ! ! A total of 18 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.02 FD 1.02 HN 1.11 KP 0.98 ! LA 0.98 MK 1.00 NL 1.03 OV 1.04 ! PT 0.90 SC 1.00 ! ! ! Correcting IF 2. ! A total of 102 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* ON 0.00* WZ 0.00* ! ! A total of 30 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! HH 0.00* ON 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 1.02 HN 1.00 KP 1.04 LA 1.18 ! NL 0.97 ON 0.96 PT 1.12 SC 0.93 ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.06 FD 0.93 HH 1.09 HN 0.97 ! KP 1.05 LA 1.13 NL 0.99 OV 1.05 ! PT 1.01 SC 0.98 ! ! A total of 18 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.02 FD 1.00 HN 1.12 KP 0.97 ! LA 0.98 MK 1.03 NL 1.01 OV 1.00 ! PT 0.97 SC 1.03 ! ! ! Correcting IF 3. ! A total of 102 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* ON 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! HH 0.00* ON 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 1.10 HN 1.01 KP 1.00* LA 0.93 ! NL 1.00* ON 1.05 PT 0.95 SC 1.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 1.09 FD 0.89 HH 0.91 HN 1.09 ! KP 1.00* LA 0.97 NL 1.00* OV 0.96 ! PT 0.92 SC 1.00* ! ! A total of 21 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.02 FD 0.98 HN 1.10 KP 1.00* ! LA 1.05 MK 0.91 NL 1.00* OV 1.02 ! PT 0.99 SC 1.00* ! ! ! Correcting IF 4. ! A total of 102 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* ON 0.00* WZ 0.00* ! ! A total of 30 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! HH 0.00* ON 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 1.07 HN 1.03 KP 0.94 LA 0.89 ! NL 1.08 ON 1.04 PT 0.93 SC 1.01 ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.04 FD 0.96 HH 0.94 HN 1.08 ! KP 1.02 LA 0.97 NL 1.03 OV 0.92 ! PT 0.91 SC 1.00 ! ! A total of 18 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 0.99 FD 0.97 HN 1.12 KP 0.94 ! LA 0.97 MK 1.04 NL 1.01 OV 0.99 ! PT 0.96 SC 1.04 ! ! ! Fit before self-cal, rms=0.128129Jy sigma=1.591285 ! Fit after self-cal, rms=0.124775Jy sigma=1.670468 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.133 mas, bmaj=6.239 mas, bpa=15.33 degrees ! Estimated noise=0.961135 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.00268144 Jy ! Component: 100 - total flux cleaned = 0.00348635 Jy ! Component: 150 - total flux cleaned = 0.00368483 Jy ! Component: 200 - total flux cleaned = 0.00383046 Jy ! Total flux subtracted in 200 components = 0.00383046 Jy ! Clean residual min=-0.009790 max=0.010906 Jy/beam ! Clean residual mean=0.000041 rms=0.003079 Jy/beam ! Combined flux in latest and established models = 0.437665 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 55 components and 0.437665 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.124736Jy sigma=1.669503 ! Fit after self-cal, rms=0.124870Jy sigma=1.669027 ! 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.124870Jy sigma=1.669027 ! Fit after self-cal, rms=0.133740Jy sigma=1.577709 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.133 mas, bmaj=6.251 mas, bpa=15.39 degrees ! Estimated noise=0.969795 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00292247 Jy ! Component: 100 - total flux cleaned = 0.00363414 Jy ! Component: 150 - total flux cleaned = 0.00368215 Jy ! Component: 200 - total flux cleaned = 0.00368274 Jy ! Total flux subtracted in 200 components = 0.00368274 Jy ! Clean residual min=-0.009556 max=0.012798 Jy/beam ! Clean residual mean=0.000004 rms=0.002848 Jy/beam ! Combined flux in latest and established models = 0.441348 Jy selfcal ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 61 components and 0.441348 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.133699Jy sigma=1.576384 ! Fit after self-cal, rms=0.133785Jy sigma=1.575928 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 (-3.5, -45). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.011801 Jy ! Component: 100 - total flux cleaned = 0.0176801 Jy ! Component: 150 - total flux cleaned = 0.0208926 Jy ! Component: 200 - total flux cleaned = 0.0224416 Jy ! Total flux subtracted in 200 components = 0.0224416 Jy ! Clean residual min=-0.008479 max=0.010475 Jy/beam ! Clean residual mean=0.000006 rms=0.002608 Jy/beam ! Combined flux in latest and established models = 0.463789 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 86 components and 0.463789 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.133014Jy sigma=1.554690 ! Fit after self-cal, rms=0.132905Jy sigma=1.553964 ! 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.132905Jy sigma=1.553964 ! Fit after self-cal, rms=0.148087Jy sigma=1.484038 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.143 mas, bmaj=6.272 mas, bpa=15.3 degrees ! Estimated noise=0.969642 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00217391 Jy ! Component: 100 - total flux cleaned = 0.00317898 Jy ! Component: 150 - total flux cleaned = 0.00351345 Jy ! Component: 200 - total flux cleaned = 0.00366012 Jy ! Total flux subtracted in 200 components = 0.00366012 Jy ! Clean residual min=-0.008141 max=0.010465 Jy/beam ! Clean residual mean=-0.000043 rms=0.002531 Jy/beam ! Combined flux in latest and established models = 0.467449 Jy selfcal ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 101 components and 0.467449 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.148004Jy sigma=1.480532 ! Fit after self-cal, rms=0.148014Jy sigma=1.480378 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.000113957 Jy ! Component: 100 - total flux cleaned = 0.000243203 Jy ! Component: 150 - total flux cleaned = 0.00041709 Jy ! Component: 200 - total flux cleaned = 0.000469485 Jy ! Total flux subtracted in 200 components = 0.000469485 Jy ! Clean residual min=-0.008214 max=0.010505 Jy/beam ! Clean residual mean=-0.000041 rms=0.002535 Jy/beam ! Combined flux in latest and established models = 0.467919 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 114 components and 0.467919 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.147966Jy sigma=1.479590 ! Fit after self-cal, rms=0.147969Jy sigma=1.479563 ! 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.147969Jy sigma=1.479563 ! Fit after self-cal, rms=0.148197Jy sigma=1.479047 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.733 mas, bmaj=14.53 mas, bpa=-3.543 degrees ! Estimated noise=1.10609 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00136997 Jy ! Component: 100 - total flux cleaned = 0.00143301 Jy ! Component: 150 - total flux cleaned = 0.00143467 Jy ! Component: 200 - total flux cleaned = 0.0013937 Jy ! Total flux subtracted in 200 components = 0.0013937 Jy ! Clean residual min=-0.010148 max=0.012058 Jy/beam ! Clean residual mean=-0.000080 rms=0.003634 Jy/beam ! Combined flux in latest and established models = 0.469312 Jy selfcal ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 117 components and 0.469312 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.148243Jy sigma=1.477882 ! Fit after self-cal, rms=0.148237Jy sigma=1.477820 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 = 9.24156e-05 Jy ! Component: 100 - total flux cleaned = 0.000131489 Jy ! Component: 150 - total flux cleaned = 9.44127e-05 Jy ! Component: 200 - total flux cleaned = 9.46533e-05 Jy ! Total flux subtracted in 200 components = 9.46533e-05 Jy ! Clean residual min=-0.010205 max=0.012044 Jy/beam ! Clean residual mean=-0.000078 rms=0.003636 Jy/beam ! Combined flux in latest and established models = 0.469407 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 119 components and 0.469407 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.148245Jy sigma=1.477497 ! Fit after self-cal, rms=0.148234Jy sigma=1.477455 ! 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.139 mas, bmaj=6.271 mas, bpa=15.38 degrees ! Estimated noise=0.970513 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.139 x 6.271 at 15.38 degrees (North through East) ! Clean map min=-0.008067 max=0.36874 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=2.357 mas, bmaj=4.806 mas, bpa=18.84 degrees ! Estimated noise=1.3291 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.291692 Jy ! Component: 100 - total flux cleaned = 0.368416 Jy ! Component: 150 - total flux cleaned = 0.392628 Jy ! Component: 200 - total flux cleaned = 0.407468 Jy ! Total flux subtracted in 200 components = 0.407468 Jy ! Clean residual min=-0.009036 max=0.013267 Jy/beam ! Clean residual mean=0.000017 rms=0.002867 Jy/beam ! Combined flux in latest and established models = 0.407468 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.407468 Jy ! Inverting map and beam ! Estimated beam: bmin=3.139 mas, bmaj=6.271 mas, bpa=15.38 degrees ! Estimated noise=0.970513 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) ! Added new window around map position (67.5, -108). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0174797 Jy ! Component: 100 - total flux cleaned = 0.0293317 Jy ! Component: 150 - total flux cleaned = 0.0384455 Jy ! Component: 200 - total flux cleaned = 0.0457461 Jy ! Total flux subtracted in 200 components = 0.0457461 Jy ! Clean residual min=-0.007604 max=0.010648 Jy/beam ! Clean residual mean=0.000015 rms=0.002476 Jy/beam ! Combined flux in latest and established models = 0.453214 Jy ! Adding 48 model components to the UV plane model. ! The established model now contains 68 components and 0.453214 Jy ! Inverting map ! Added new window around map position (111.5, -120). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00906476 Jy ! Component: 100 - total flux cleaned = 0.0148457 Jy ! Component: 150 - total flux cleaned = 0.0197198 Jy ! Component: 200 - total flux cleaned = 0.0238791 Jy ! Total flux subtracted in 200 components = 0.0238791 Jy ! Clean residual min=-0.007445 max=0.008840 Jy/beam ! Clean residual mean=-0.000041 rms=0.002395 Jy/beam ! Combined flux in latest and established models = 0.477093 Jy ! Adding 35 model components to the UV plane model. ! The established model now contains 87 components and 0.477093 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.733 mas, bmaj=14.53 mas, bpa=-3.543 degrees ! Estimated noise=1.10609 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00553271 Jy ! Component: 100 - total flux cleaned = 0.00923987 Jy ! Component: 150 - total flux cleaned = 0.0109617 Jy ! Component: 200 - total flux cleaned = 0.0125296 Jy ! Total flux subtracted in 200 components = 0.0125296 Jy ! Clean residual min=-0.010064 max=0.010590 Jy/beam ! Clean residual mean=-0.000098 rms=0.003430 Jy/beam ! Combined flux in latest and established models = 0.489623 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.00108764 Jy ! Component: 100 - total flux cleaned = 0.0020466 Jy ! Component: 150 - total flux cleaned = 0.00255238 Jy ! Component: 200 - total flux cleaned = 0.00302573 Jy ! Total flux subtracted in 200 components = 0.00302573 Jy ! Clean residual min=-0.010185 max=0.010618 Jy/beam ! Clean residual mean=-0.000106 rms=0.003436 Jy/beam ! Combined flux in latest and established models = 0.492649 Jy ! Adding 54 model components to the UV plane model. ! The established model now contains 133 components and 0.492649 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.139 mas, bmaj=6.271 mas, bpa=15.38 degrees ! Estimated noise=0.970513 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00145838 Jy ! Component: 100 - total flux cleaned = -0.00135562 Jy ! Component: 150 - total flux cleaned = -0.0012652 Jy ! Component: 200 - total flux cleaned = -0.00117695 Jy ! Total flux subtracted in 200 components = -0.00117695 Jy ! Clean residual min=-0.007968 max=0.009594 Jy/beam ! Clean residual mean=-0.000063 rms=0.002399 Jy/beam ! Combined flux in latest and established models = 0.491472 Jy selfcal ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 150 components and 0.491472 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.146245Jy sigma=1.457649 ! Fit after self-cal, rms=0.145847Jy sigma=1.452545 wmodel J0405-1308_S_map.mod ! Writing 150 model components to file: J0405-1308_S_map.mod wobs J0405-1308_S_uvs.fits ! Writing UV FITS file: J0405-1308_S_uvs.fits wwins J0405-1308_S_map.win ! wwins: Wrote 5 windows to J0405-1308_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.0042025 Jy ! Component: 100 - total flux cleaned = 0.00743255 Jy ! Component: 150 - total flux cleaned = 0.00975448 Jy ! Component: 200 - total flux cleaned = 0.011711 Jy ! Component: 250 - total flux cleaned = 0.0131995 Jy ! Component: 300 - total flux cleaned = 0.0145146 Jy ! Component: 350 - total flux cleaned = 0.0157964 Jy ! Component: 400 - total flux cleaned = 0.016922 Jy ! Component: 450 - total flux cleaned = 0.0180205 Jy ! Component: 500 - total flux cleaned = 0.0187367 Jy ! Component: 550 - total flux cleaned = 0.019555 Jy ! Component: 600 - total flux cleaned = 0.0204703 Jy ! Component: 650 - total flux cleaned = 0.0210304 Jy ! Component: 700 - total flux cleaned = 0.0220207 Jy ! Component: 750 - total flux cleaned = 0.0227788 Jy ! Component: 800 - total flux cleaned = 0.0234177 Jy ! Component: 850 - total flux cleaned = 0.0242543 Jy ! Component: 900 - total flux cleaned = 0.0249729 Jy ! Component: 950 - total flux cleaned = 0.0252772 Jy ! Component: 1000 - total flux cleaned = 0.0259752 Jy ! Component: 1050 - total flux cleaned = 0.0266639 Jy ! Component: 1100 - total flux cleaned = 0.0272472 Jy ! Component: 1150 - total flux cleaned = 0.028204 Jy ! Component: 1200 - total flux cleaned = 0.0286755 Jy ! Component: 1250 - total flux cleaned = 0.0291416 Jy ! Component: 1300 - total flux cleaned = 0.0295106 Jy ! Component: 1350 - total flux cleaned = 0.0297847 Jy ! Component: 1400 - total flux cleaned = 0.030506 Jy ! Component: 1450 - total flux cleaned = 0.0309506 Jy ! Component: 1500 - total flux cleaned = 0.0309504 Jy ! Component: 1550 - total flux cleaned = 0.0316486 Jy ! Component: 1600 - total flux cleaned = 0.0319933 Jy ! Component: 1650 - total flux cleaned = 0.0322502 Jy ! Component: 1700 - total flux cleaned = 0.0325891 Jy ! Component: 1750 - total flux cleaned = 0.0327559 Jy ! Component: 1800 - total flux cleaned = 0.033253 Jy ! Component: 1850 - total flux cleaned = 0.0332528 Jy ! Component: 1900 - total flux cleaned = 0.0334158 Jy ! Component: 1950 - total flux cleaned = 0.0337372 Jy ! Component: 2000 - total flux cleaned = 0.0341365 Jy ! Component: 2050 - total flux cleaned = 0.0340569 Jy ! Component: 2100 - total flux cleaned = 0.0344491 Jy ! Component: 2150 - total flux cleaned = 0.0343712 Jy ! Component: 2200 - total flux cleaned = 0.0345255 Jy ! Component: 2250 - total flux cleaned = 0.0347551 Jy ! Component: 2300 - total flux cleaned = 0.0349834 Jy ! Component: 2350 - total flux cleaned = 0.0352078 Jy ! Component: 2400 - total flux cleaned = 0.0354321 Jy ! Component: 2450 - total flux cleaned = 0.0355812 Jy ! Component: 2500 - total flux cleaned = 0.0359485 Jy ! Component: 2550 - total flux cleaned = 0.0358761 Jy ! Component: 2600 - total flux cleaned = 0.0363094 Jy ! Component: 2650 - total flux cleaned = 0.0365254 Jy ! Component: 2700 - total flux cleaned = 0.0366689 Jy ! Component: 2750 - total flux cleaned = 0.0370934 Jy ! Component: 2800 - total flux cleaned = 0.037304 Jy ! Component: 2850 - total flux cleaned = 0.0373033 Jy ! Component: 2900 - total flux cleaned = 0.0373721 Jy ! Component: 2950 - total flux cleaned = 0.0373717 Jy ! Component: 3000 - total flux cleaned = 0.0372343 Jy ! Component: 3050 - total flux cleaned = 0.0373028 Jy ! Component: 3100 - total flux cleaned = 0.0377086 Jy ! Component: 3150 - total flux cleaned = 0.0376407 Jy ! Component: 3200 - total flux cleaned = 0.0378421 Jy ! Component: 3250 - total flux cleaned = 0.0381083 Jy ! Component: 3300 - total flux cleaned = 0.0383729 Jy ! Component: 3350 - total flux cleaned = 0.0387014 Jy ! Component: 3400 - total flux cleaned = 0.0383747 Jy ! Component: 3450 - total flux cleaned = 0.0384404 Jy ! Component: 3500 - total flux cleaned = 0.0386344 Jy ! Component: 3550 - total flux cleaned = 0.0386986 Jy ! Component: 3600 - total flux cleaned = 0.038635 Jy ! Component: 3650 - total flux cleaned = 0.0390172 Jy ! Component: 3700 - total flux cleaned = 0.0390163 Jy ! Component: 3750 - total flux cleaned = 0.0389529 Jy ! Component: 3800 - total flux cleaned = 0.0391407 Jy ! Component: 3850 - total flux cleaned = 0.0387667 Jy ! Component: 3900 - total flux cleaned = 0.0390764 Jy ! Component: 3950 - total flux cleaned = 0.0388914 Jy ! Component: 4000 - total flux cleaned = 0.0388918 Jy ! Component: 4050 - total flux cleaned = 0.0389527 Jy ! Total flux subtracted in 4096 components = 0.0390136 Jy ! Clean residual min=-0.003966 max=0.004113 Jy/beam ! Clean residual mean=-0.000044 rms=0.001313 Jy/beam ! Combined flux in latest and established models = 0.530485 Jy keep ! Adding 362 model components to the UV plane model. ! The established model now contains 512 components and 0.530486 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /scr/pima/rv121_uvs/J0405-1308_S_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.139 x 6.271 at 15.38 degrees (North through East) ! Clean map min=-0.0062278 max=0.37196 Jy/beam ! Writing clean map to FITS file: /scr/pima/rv121_uvs/J0405-1308_S_map.fits wmodel /scr/pima/rv121_uvs/J0405-1308_S_map.mod ! Writing 512 model components to file: /scr/pima/rv121_uvs/J0405-1308_S_map.mod wobs /scr/pima/rv121_uvs/J0405-1308_S_uvs.fits ! Writing UV FITS file: /scr/pima/rv121_uvs/J0405-1308_S_uvs.fits wwins /scr/pima/rv121_uvs/J0405-1308_S_map.win ! wwins: Wrote 6 windows to /scr/pima/rv121_uvs/J0405-1308_S_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.371963 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.0013139 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -1.05971 1.05971 2.11941 4.23883 8.47766 16.9553 33.9106 67.8212 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0039 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.371 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 283.097 quit ! Quitting program ! Log file /scr/pima/rv121_uvs/J0405-1308_S_dfm.log closed on Sat Jul 28 00:08:33 2018