! Started logfile: /scr/pima/rv126_uvs/J1500+4751_S_dfm.log on Sun Dec 23 07:40:52 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.24075 obs /scr/pima/rv126_uvs/J1500+4751_S_uva.fits ! Reading UV FITS file: /scr/pima/rv126_uvs/J1500+4751_S_uva.fits ! AN table 1: 10 integrations on 6 of 6 possible baselines. ! AN table 2: 46 integrations on 3 of 3 possible baselines. ! AN table 3: 4 integrations on 21 of 21 possible baselines. ! AN table 4: 22 integrations on 36 of 36 possible baselines. ! AN table 5: 38 integrations on 28 of 28 possible baselines. ! Apparent sampling: 0.434518 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1500+4751 ! ! 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.24075e+09 8e+06 1 8e+06 ! 02 2 2.27075e+09 8e+06 1 8e+06 ! 03 3 2.36075e+09 8e+06 1 8e+06 ! 04 4 2.38075e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 69 lines of history. ! ! Reading 3716 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 J1500+4751_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 22 telescope corrections were flagged in sub-array 1. ! A total of 138 telescope corrections were flagged in sub-array 2. ! A total of 2 telescope corrections were flagged in sub-array 3. ! A total of 81 telescope corrections were flagged in sub-array 4. ! A total of 198 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 22 telescope corrections were flagged in sub-array 1. ! A total of 138 telescope corrections were flagged in sub-array 2. ! A total of 2 telescope corrections were flagged in sub-array 3. ! A total of 81 telescope corrections were flagged in sub-array 4. ! A total of 198 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 40 telescope corrections were flagged in sub-array 1. ! A total of 138 telescope corrections were flagged in sub-array 2. ! A total of 28 telescope corrections were flagged in sub-array 3. ! A total of 198 telescope corrections were flagged in sub-array 4. ! A total of 304 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 22 telescope corrections were flagged in sub-array 1. ! A total of 138 telescope corrections were flagged in sub-array 2. ! A total of 2 telescope corrections were flagged in sub-array 3. ! A total of 81 telescope corrections were flagged in sub-array 4. ! A total of 198 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=1.143155Jy sigma=16.894898 ! Fit after self-cal, rms=0.522896Jy sigma=7.931705 ! 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.351 mas, bmaj=4.479 mas, bpa=-41.91 degrees ! Estimated noise=1.4282 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.359571 Jy ! Component: 100 - total flux cleaned = 0.45428 Jy ! Total flux subtracted in 100 components = 0.45428 Jy ! Clean residual min=-0.020143 max=0.040873 Jy/beam ! Clean residual mean=0.000119 rms=0.005330 Jy/beam ! Combined flux in latest and established models = 0.45428 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 20 components and 0.45428 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.131156Jy sigma=1.954326 ! Fit after self-cal, rms=0.131155Jy sigma=1.954189 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.044903 Jy ! Component: 100 - total flux cleaned = 0.0687154 Jy ! Total flux subtracted in 100 components = 0.0687154 Jy ! Clean residual min=-0.017323 max=0.029171 Jy/beam ! Clean residual mean=0.000043 rms=0.004421 Jy/beam ! Combined flux in latest and established models = 0.522995 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 33 components and 0.522995 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.109157Jy sigma=1.618232 ! Fit after self-cal, rms=0.109143Jy sigma=1.617760 ! Inverting map ! Added new window around map position (7, 2.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0263677 Jy ! Component: 100 - total flux cleaned = 0.0390977 Jy ! Total flux subtracted in 100 components = 0.0390977 Jy ! Clean residual min=-0.015699 max=0.019092 Jy/beam ! Clean residual mean=0.000004 rms=0.004068 Jy/beam ! Combined flux in latest and established models = 0.562093 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 41 components and 0.562093 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.102490Jy sigma=1.511248 ! Fit after self-cal, rms=0.101528Jy sigma=1.496169 ! 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.212 mas, bmaj=5.719 mas, bpa=-43.53 degrees ! Estimated noise=1.10666 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.0131018 Jy ! Component: 100 - total flux cleaned = 0.0137559 Jy ! Component: 150 - total flux cleaned = 0.0132534 Jy ! Component: 200 - total flux cleaned = 0.0127861 Jy ! Total flux subtracted in 200 components = 0.0127861 Jy ! Clean residual min=-0.013059 max=0.015484 Jy/beam ! Clean residual mean=-0.000039 rms=0.003729 Jy/beam ! Combined flux in latest and established models = 0.574879 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 58 components and 0.574879 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.099564Jy sigma=1.461575 ! Fit after self-cal, rms=0.099341Jy sigma=1.457677 ! 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.212 x 5.719 at -43.53 degrees (North through East) ! Clean map min=-0.013389 max=0.4757 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 18 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* HN 0.00* NY 0.00* WN 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! NY 0.00* WN 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 1.01 KP 0.95 LA 0.92 MK 0.97 ! NY 1.19 OV 0.93 PT 1.05 ! ! A total of 9 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.03 FD 1.03 HN 0.87 KP 0.94 ! LA 0.93 MK 0.98 NL 1.11 OV 0.88 ! PT 1.07 ! ! A total of 54 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.01 FD 1.10 HN 0.87 KP 1.01 ! LA 0.96 NL 1.14 NY 1.26 PT 1.07 ! ! ! Correcting IF 2. ! A total of 18 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* HN 0.00* NY 0.00* WN 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! NY 0.00* WN 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 1.03 KP 0.92 LA 0.90 MK 0.92 ! NY 1.03 OV 0.85 PT 1.07 ! ! A total of 9 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.01 FD 1.02 HN 0.92 KP 0.94 ! LA 0.93 MK 1.02 NL 1.23 OV 0.85 ! PT 1.08 ! ! A total of 54 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.02 FD 1.11 HN 0.88 KP 1.03 ! LA 0.96 NL 1.13 NY 1.07 PT 1.08 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* HN 0.00* NY 0.00* WN 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! NY 0.00* WN 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 0.00* KP 0.00* LA 0.00* MK 0.00* ! NY 0.00* OV 0.00* PT 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* FD 0.00* HN 0.00* KP 0.00* ! LA 0.00* MK 0.00* NL 0.00* OV 0.00* ! PT 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 0.00* FD 0.00* HN 0.00* KP 0.00* ! LA 0.00* NL 0.00* NY 0.00* PT 0.00* ! ! ! Correcting IF 4. ! A total of 18 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* HN 0.00* NY 0.00* WN 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! NY 0.00* WN 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 1.02 KP 0.88 LA 0.93 MK 0.99 ! NY 1.18 OV 0.92 PT 1.03 ! ! A total of 9 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.02 FD 1.04 HN 0.91 KP 0.90 ! LA 0.93 MK 0.97 NL 1.09 OV 0.92 ! PT 1.05 ! ! A total of 54 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.07 FD 1.08 HN 0.89 KP 1.04 ! LA 0.93 NL 1.17 NY 1.28 PT 1.09 ! ! ! Fit before self-cal, rms=0.099341Jy sigma=1.457677 ! Fit after self-cal, rms=0.072369Jy sigma=1.040743 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=4.33 mas, bmaj=6.1 mas, bpa=-47.8 degrees ! Estimated noise=1.14895 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.00329967 Jy ! Component: 100 - total flux cleaned = -0.00347832 Jy ! Component: 150 - total flux cleaned = -0.00347858 Jy ! Component: 200 - total flux cleaned = -0.00360671 Jy ! Total flux subtracted in 200 components = -0.00360671 Jy ! Clean residual min=-0.009242 max=0.008241 Jy/beam ! Clean residual mean=0.000017 rms=0.001810 Jy/beam ! Combined flux in latest and established models = 0.571272 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 69 components and 0.571272 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.071678Jy sigma=1.027527 ! Fit after self-cal, rms=0.071747Jy sigma=1.026434 ! 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.071747Jy sigma=1.026434 ! Fit after self-cal, rms=0.066822Jy sigma=0.917340 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=4.349 mas, bmaj=6.102 mas, bpa=-47.6 degrees ! Estimated noise=1.1482 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000145608 Jy ! Component: 100 - total flux cleaned = -0.000144273 Jy ! Component: 150 - total flux cleaned = -5.10372e-05 Jy ! Component: 200 - total flux cleaned = -5.53597e-05 Jy ! Total flux subtracted in 200 components = -5.53597e-05 Jy ! Clean residual min=-0.005013 max=0.005443 Jy/beam ! Clean residual mean=0.000004 rms=0.001186 Jy/beam ! Combined flux in latest and established models = 0.571217 Jy selfcal ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 81 components and 0.571217 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.066470Jy sigma=0.912210 ! Fit after self-cal, rms=0.066456Jy sigma=0.911716 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 (16.5, 17). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00501275 Jy ! Component: 100 - total flux cleaned = 0.00692737 Jy ! Component: 150 - total flux cleaned = 0.00720559 Jy ! Component: 200 - total flux cleaned = 0.00745731 Jy ! Total flux subtracted in 200 components = 0.00745731 Jy ! Clean residual min=-0.004601 max=0.004043 Jy/beam ! Clean residual mean=-0.000000 rms=0.001077 Jy/beam ! Combined flux in latest and established models = 0.578674 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 108 components and 0.578674 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.066039Jy sigma=0.903569 ! Fit after self-cal, rms=0.066019Jy sigma=0.903291 ! 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.066019Jy sigma=0.903291 ! Fit after self-cal, rms=0.065868Jy sigma=0.899617 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=4.354 mas, bmaj=6.117 mas, bpa=-47.62 degrees ! Estimated noise=1.14727 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000406483 Jy ! Component: 100 - total flux cleaned = 0.000944138 Jy ! Component: 150 - total flux cleaned = 0.00125868 Jy ! Component: 200 - total flux cleaned = 0.00147325 Jy ! Total flux subtracted in 200 components = 0.00147325 Jy ! Clean residual min=-0.004111 max=0.003589 Jy/beam ! Clean residual mean=0.000016 rms=0.001015 Jy/beam ! Combined flux in latest and established models = 0.580148 Jy selfcal ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 113 components and 0.580148 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.065659Jy sigma=0.896010 ! Fit after self-cal, rms=0.065644Jy sigma=0.895763 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.000152678 Jy ! Component: 100 - total flux cleaned = 0.000284023 Jy ! Component: 150 - total flux cleaned = 0.000406316 Jy ! Component: 200 - total flux cleaned = 0.00051793 Jy ! Total flux subtracted in 200 components = 0.00051793 Jy ! Clean residual min=-0.003962 max=0.003451 Jy/beam ! Clean residual mean=0.000016 rms=0.000997 Jy/beam ! Combined flux in latest and established models = 0.580666 Jy ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 119 components and 0.580666 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.065530Jy sigma=0.893890 ! Fit after self-cal, rms=0.065523Jy sigma=0.893772 ! 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.065523Jy sigma=0.893772 ! Fit after self-cal, rms=0.065592Jy sigma=0.893112 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=10.15 mas, bmaj=12.03 mas, bpa=-69.05 degrees ! Estimated noise=1.34391 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 3.83813e-05 Jy ! Component: 100 - total flux cleaned = 0.000429506 Jy ! Component: 150 - total flux cleaned = 0.000629126 Jy ! Component: 200 - total flux cleaned = 0.000696665 Jy ! Total flux subtracted in 200 components = 0.000696665 Jy ! Clean residual min=-0.003613 max=0.003389 Jy/beam ! Clean residual mean=0.000042 rms=0.001068 Jy/beam ! Combined flux in latest and established models = 0.581362 Jy selfcal ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 131 components and 0.581362 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.065532Jy sigma=0.892284 ! Fit after self-cal, rms=0.065529Jy sigma=0.892163 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 = -6.05629e-06 Jy ! Component: 100 - total flux cleaned = 0.000126854 Jy ! Component: 150 - total flux cleaned = 0.000243311 Jy ! Component: 200 - total flux cleaned = 0.000321869 Jy ! Total flux subtracted in 200 components = 0.000321869 Jy ! Clean residual min=-0.003581 max=0.003407 Jy/beam ! Clean residual mean=0.000041 rms=0.001060 Jy/beam ! Combined flux in latest and established models = 0.581684 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 142 components and 0.581684 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.065500Jy sigma=0.891666 ! Fit after self-cal, rms=0.065498Jy sigma=0.891591 ! 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.357 mas, bmaj=6.131 mas, bpa=-47.62 degrees ! Estimated noise=1.14683 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 4.357 x 6.131 at -47.62 degrees (North through East) ! Clean map min=-0.0036953 max=0.48037 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.352 mas, bmaj=4.716 mas, bpa=-43.1 degrees ! Estimated noise=1.50366 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.366633 Jy ! Component: 100 - total flux cleaned = 0.45889 Jy ! Component: 150 - total flux cleaned = 0.501856 Jy ! Component: 200 - total flux cleaned = 0.527411 Jy ! Total flux subtracted in 200 components = 0.527411 Jy ! Clean residual min=-0.007111 max=0.014016 Jy/beam ! Clean residual mean=0.000068 rms=0.001913 Jy/beam ! Combined flux in latest and established models = 0.527411 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 31 model components to the UV plane model. ! The established model now contains 31 components and 0.527411 Jy ! Inverting map and beam ! Estimated beam: bmin=4.357 mas, bmaj=6.131 mas, bpa=-47.62 degrees ! Estimated noise=1.14683 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.021929 Jy ! Component: 100 - total flux cleaned = 0.0331114 Jy ! Component: 150 - total flux cleaned = 0.0401493 Jy ! Component: 200 - total flux cleaned = 0.0450071 Jy ! Total flux subtracted in 200 components = 0.0450071 Jy ! Clean residual min=-0.004212 max=0.004131 Jy/beam ! Clean residual mean=0.000032 rms=0.001075 Jy/beam ! Combined flux in latest and established models = 0.572418 Jy ! Adding 42 model components to the UV plane model. ! The established model now contains 60 components and 0.572418 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=10.15 mas, bmaj=12.03 mas, bpa=-69.05 degrees ! Estimated noise=1.34391 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0043663 Jy ! Component: 100 - total flux cleaned = 0.00710482 Jy ! Component: 150 - total flux cleaned = 0.00893836 Jy ! Component: 200 - total flux cleaned = 0.00996709 Jy ! Total flux subtracted in 200 components = 0.00996709 Jy ! Clean residual min=-0.003536 max=0.003538 Jy/beam ! Clean residual mean=0.000039 rms=0.001086 Jy/beam ! Combined flux in latest and established models = 0.582385 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.000125211 Jy ! Component: 100 - total flux cleaned = 0.00012542 Jy ! Component: 150 - total flux cleaned = 0.000125347 Jy ! Component: 200 - total flux cleaned = 0.000125257 Jy ! Total flux subtracted in 200 components = 0.000125257 Jy ! Clean residual min=-0.003500 max=0.003481 Jy/beam ! Clean residual mean=0.000039 rms=0.001075 Jy/beam ! Combined flux in latest and established models = 0.58251 Jy ! Adding 40 model components to the UV plane model. ! The established model now contains 94 components and 0.582511 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=4.357 mas, bmaj=6.131 mas, bpa=-47.62 degrees ! Estimated noise=1.14683 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000620687 Jy ! Component: 100 - total flux cleaned = -0.000275409 Jy ! Component: 150 - total flux cleaned = -9.52481e-05 Jy ! Component: 200 - total flux cleaned = 7.54747e-05 Jy ! Total flux subtracted in 200 components = 7.54747e-05 Jy ! Clean residual min=-0.003828 max=0.003763 Jy/beam ! Clean residual mean=0.000023 rms=0.000981 Jy/beam ! Combined flux in latest and established models = 0.582586 Jy selfcal ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 108 components and 0.582586 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.065502Jy sigma=0.892527 ! Fit after self-cal, rms=0.065395Jy sigma=0.891047 wmodel J1500+4751_S_map.mod ! Writing 108 model components to file: J1500+4751_S_map.mod wobs J1500+4751_S_uvs.fits ! Writing UV FITS file: J1500+4751_S_uvs.fits wwins J1500+4751_S_map.win ! wwins: Wrote 3 windows to J1500+4751_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.000149061 Jy ! Component: 100 - total flux cleaned = -0.00060369 Jy ! Component: 150 - total flux cleaned = -0.000915366 Jy ! Component: 200 - total flux cleaned = -0.00109637 Jy ! Component: 250 - total flux cleaned = -0.00121423 Jy ! Component: 300 - total flux cleaned = -0.00127124 Jy ! Component: 350 - total flux cleaned = -0.00132602 Jy ! Component: 400 - total flux cleaned = -0.00127236 Jy ! Component: 450 - total flux cleaned = -0.001326 Jy ! Component: 500 - total flux cleaned = -0.00117065 Jy ! Component: 550 - total flux cleaned = -0.00101919 Jy ! Component: 600 - total flux cleaned = -0.00106957 Jy ! Component: 650 - total flux cleaned = -0.00106925 Jy ! Component: 700 - total flux cleaned = -0.00106776 Jy ! Component: 750 - total flux cleaned = -0.00087607 Jy ! Component: 800 - total flux cleaned = -0.000687621 Jy ! Component: 850 - total flux cleaned = -0.000686978 Jy ! Component: 900 - total flux cleaned = -0.00045678 Jy ! Component: 950 - total flux cleaned = -0.000366319 Jy ! Component: 1000 - total flux cleaned = -5.27205e-05 Jy ! Component: 1050 - total flux cleaned = 0.000123759 Jy ! Component: 1100 - total flux cleaned = 0.000212084 Jy ! Component: 1150 - total flux cleaned = 0.000515099 Jy ! Component: 1200 - total flux cleaned = 0.000643578 Jy ! Component: 1250 - total flux cleaned = 0.000770716 Jy ! Component: 1300 - total flux cleaned = 0.000980861 Jy ! Component: 1350 - total flux cleaned = 0.00106364 Jy ! Component: 1400 - total flux cleaned = 0.00126915 Jy ! Component: 1450 - total flux cleaned = 0.00143195 Jy ! Component: 1500 - total flux cleaned = 0.00155328 Jy ! Component: 1550 - total flux cleaned = 0.00179331 Jy ! Component: 1600 - total flux cleaned = 0.00199138 Jy ! Component: 1650 - total flux cleaned = 0.00206945 Jy ! Component: 1700 - total flux cleaned = 0.00222504 Jy ! Component: 1750 - total flux cleaned = 0.00245612 Jy ! Component: 1800 - total flux cleaned = 0.00264694 Jy ! Component: 1850 - total flux cleaned = 0.00276011 Jy ! Component: 1900 - total flux cleaned = 0.00291027 Jy ! Component: 1950 - total flux cleaned = 0.00313381 Jy ! Component: 2000 - total flux cleaned = 0.0032451 Jy ! Component: 2050 - total flux cleaned = 0.00346528 Jy ! Component: 2100 - total flux cleaned = 0.00346539 Jy ! Component: 2150 - total flux cleaned = 0.00368236 Jy ! Component: 2200 - total flux cleaned = 0.00382586 Jy ! Component: 2250 - total flux cleaned = 0.00396835 Jy ! Component: 2300 - total flux cleaned = 0.00400386 Jy ! Component: 2350 - total flux cleaned = 0.00425006 Jy ! Component: 2400 - total flux cleaned = 0.00431976 Jy ! Component: 2450 - total flux cleaned = 0.00442362 Jy ! Component: 2500 - total flux cleaned = 0.00459511 Jy ! Component: 2550 - total flux cleaned = 0.00459464 Jy ! Component: 2600 - total flux cleaned = 0.00462846 Jy ! Component: 2650 - total flux cleaned = 0.00469617 Jy ! Component: 2700 - total flux cleaned = 0.00479646 Jy ! Component: 2750 - total flux cleaned = 0.00486251 Jy ! Component: 2800 - total flux cleaned = 0.00486262 Jy ! Component: 2850 - total flux cleaned = 0.00496087 Jy ! Component: 2900 - total flux cleaned = 0.00509091 Jy ! Component: 2950 - total flux cleaned = 0.00512336 Jy ! Component: 3000 - total flux cleaned = 0.00531578 Jy ! Component: 3050 - total flux cleaned = 0.00525199 Jy ! Component: 3100 - total flux cleaned = 0.00534696 Jy ! Component: 3150 - total flux cleaned = 0.00540989 Jy ! Component: 3200 - total flux cleaned = 0.00540999 Jy ! Component: 3250 - total flux cleaned = 0.00550297 Jy ! Component: 3300 - total flux cleaned = 0.00559545 Jy ! Component: 3350 - total flux cleaned = 0.00562611 Jy ! Component: 3400 - total flux cleaned = 0.00574789 Jy ! Component: 3450 - total flux cleaned = 0.00574775 Jy ! Component: 3500 - total flux cleaned = 0.00583791 Jy ! Component: 3550 - total flux cleaned = 0.00586733 Jy ! Component: 3600 - total flux cleaned = 0.00592679 Jy ! Component: 3650 - total flux cleaned = 0.0059562 Jy ! Component: 3700 - total flux cleaned = 0.00601505 Jy ! Component: 3750 - total flux cleaned = 0.0059855 Jy ! Component: 3800 - total flux cleaned = 0.00604357 Jy ! Component: 3850 - total flux cleaned = 0.00595708 Jy ! Component: 3900 - total flux cleaned = 0.00607204 Jy ! Component: 3950 - total flux cleaned = 0.00615768 Jy ! Component: 4000 - total flux cleaned = 0.00618613 Jy ! Component: 4050 - total flux cleaned = 0.00629898 Jy ! Total flux subtracted in 4096 components = 0.00624286 Jy ! Clean residual min=-0.001399 max=0.001741 Jy/beam ! Clean residual mean=0.000007 rms=0.000520 Jy/beam ! Combined flux in latest and established models = 0.588829 Jy keep ! Adding 282 model components to the UV plane model. ! The established model now contains 390 components and 0.588829 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /scr/pima/rv126_uvs/J1500+4751_S_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 4.357 x 6.131 at -47.62 degrees (North through East) ! Clean map min=-0.0030872 max=0.48036 Jy/beam ! Writing clean map to FITS file: /scr/pima/rv126_uvs/J1500+4751_S_map.fits wmodel /scr/pima/rv126_uvs/J1500+4751_S_map.mod ! Writing 390 model components to file: /scr/pima/rv126_uvs/J1500+4751_S_map.mod wobs /scr/pima/rv126_uvs/J1500+4751_S_uvs.fits ! Writing UV FITS file: /scr/pima/rv126_uvs/J1500+4751_S_uvs.fits wwins /scr/pima/rv126_uvs/J1500+4751_S_map.win ! wwins: Wrote 4 windows to /scr/pima/rv126_uvs/J1500+4751_S_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.480357 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00051765 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.323291 0.323291 0.646582 1.29316 2.58633 5.17265 10.3453 20.6906 41.3812 82.7625 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0015 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.48 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 927.956 quit ! Quitting program ! Log file /scr/pima/rv126_uvs/J1500+4751_S_dfm.log closed on Sun Dec 23 07:40:57 2018