! Started logfile: /scr/pima/rdv93_uvs/J1107-4449_X_uvs.log on Tue Oct 8 15:46:35 2013 obs /scr/pima/rdv93_uvs/J1107-4449_X_uva.fits ! Reading UV FITS file: /scr/pima/rdv93_uvs/J1107-4449_X_uva.fits ! AN table 1: 5 integrations on 1 of 1 possible baselines. ! AN table 2: 2 integrations on 1 of 1 possible baselines. ! AN table 3: 3 integrations on 1 of 1 possible baselines. ! AN table 4: 2 integrations on 3 of 3 possible baselines. ! AN table 5: 2 integrations on 10 of 10 possible baselines. ! Apparent sampling: 0.666667 visibilities/baseline/integration-bin. ! Found source: J1107-4449 ! ! 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 8.40599e+09 1.6e+07 1 1.6e+07 ! 02 2 8.47599e+09 1.6e+07 1 1.6e+07 ! 03 3 8.79099e+09 1.6e+07 1 1.6e+07 ! 04 4 8.89599e+09 1.6e+07 1 1.6e+07 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 96 visibilities. select RR ! 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.2 print "taper_size =",taper_size ! taper_size = 70 ![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv93_uvs/J1107-4449_X] integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 800000 float soltime2 soltime2 = 15 float soltime3 soltime3 = 1 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float old_peak float new_peak float flux_cutoff selfant "",false ! 1:FT self-cal status=correctable weight=1 ! 1:SC self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 3:FD self-cal status=correctable weight=1 ! 3:PT self-cal status=correctable weight=1 ! 4:KP self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:MK self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 system "rm /tmp/ff.ff" mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.200x0.200 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 10 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 5 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 5 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 5 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 5 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=1.577523Jy sigma=37.086873 ! Fit after self-cal, rms=0.382003Jy sigma=9.293137 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. print "=== Flag out IVS stations" ! === Flag out IVS stations integer lun string cfg_file lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_flag_process.csh /tmp/ff.ff" ![@/tmp/ff.ff \] flag 1:FT ! Flagging baselines of antenna 1:FT in the currently selected channels. ![Exited script file: /tmp/ff.ff] 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=0.6774 mas, bmaj=5.933 mas, bpa=-17.23 degrees ! Estimated noise=7.36934 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=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.875622 Jy ! Component: 100 - total flux cleaned = 0.971334 Jy ! Total flux subtracted in 100 components = 0.971334 Jy ! Clean residual min=-0.133019 max=0.132632 Jy/beam ! Clean residual mean=0.000225 rms=0.055633 Jy/beam ! Combined flux in latest and established models = 0.971334 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 16 components and 0.971334 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.262407Jy sigma=7.116211 ! Fit after self-cal, rms=0.261792Jy sigma=7.109637 ! 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=0.7602 mas, bmaj=6.342 mas, bpa=-18.52 degrees ! Estimated noise=5.87184 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.0599509 Jy ! Component: 100 - total flux cleaned = 0.12896 Jy ! Component: 150 - total flux cleaned = 0.139047 Jy ! Component: 200 - total flux cleaned = 0.122288 Jy ! Total flux subtracted in 200 components = 0.122288 Jy ! Clean residual min=-0.030207 max=0.030052 Jy/beam ! Clean residual mean=-0.000038 rms=0.018738 Jy/beam ! Combined flux in latest and established models = 1.09362 Jy ! Performing phase self-cal ! Adding 146 model components to the UV plane model. ! The established model now contains 161 components and 1.09362 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.069557Jy sigma=1.906195 ! Fit after self-cal, rms=0.069342Jy sigma=1.903759 ! Inverting map print "=== Finished natural weighting clean for VLBI sites ===" ! === Finished natural weighting clean for VLBI sites === gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! FT 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.00* LA 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 0.00* PT 0.00* ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! KP 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! FD 0.97 LA 0.94 MK 1.01 OV 0.98 ! PT 0.95 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! FT 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.00* LA 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 0.00* PT 0.00* ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! KP 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! FD 1.04 LA 0.92 MK 1.05 OV 0.96 ! PT 0.94 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! FT 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.00* LA 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 0.00* PT 0.00* ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! KP 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! FD 1.00 LA 0.95 MK 1.03 OV 0.98 ! PT 0.98 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! FT 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.00* LA 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 0.00* PT 0.00* ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! KP 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! FD 1.08 LA 0.97 MK 1.09 OV 0.92 ! PT 0.93 ! ! ! Fit before self-cal, rms=0.069342Jy sigma=1.903759 ! Fit after self-cal, rms=0.025851Jy sigma=0.225566 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.7576 mas, bmaj=8.716 mas, bpa=-13.32 degrees ! Estimated noise=8.4203 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.00119421 Jy ! Component: 100 - total flux cleaned = -0.00431806 Jy ! Component: 150 - total flux cleaned = -0.00673004 Jy ! Component: 200 - total flux cleaned = -0.00671618 Jy ! Total flux subtracted in 200 components = -0.00671618 Jy ! Clean residual min=-0.004903 max=0.004902 Jy/beam ! Clean residual mean=-0.000090 rms=0.001437 Jy/beam ! Combined flux in latest and established models = 1.08691 Jy ! Performing phase self-cal ! Adding 48 model components to the UV plane model. ! The established model now contains 209 components and 1.08691 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.018019Jy sigma=0.155181 ! Fit after self-cal, rms=0.017695Jy sigma=0.149282 ! Inverting map selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.017695Jy sigma=0.149282 ! Fit after self-cal, rms=0.016097Jy sigma=0.140939 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.7573 mas, bmaj=8.712 mas, bpa=-13.31 degrees ! Estimated noise=8.41074 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00121871 Jy ! Component: 100 - total flux cleaned = -0.00121978 Jy ! Component: 150 - total flux cleaned = -0.00117083 Jy ! Component: 200 - total flux cleaned = -0.00112251 Jy ! Total flux subtracted in 200 components = -0.00112251 Jy ! Clean residual min=-0.007212 max=0.006985 Jy/beam ! Clean residual mean=-0.000011 rms=0.001580 Jy/beam ! Combined flux in latest and established models = 1.08578 Jy selfcal ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 239 components and 1.08578 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.012584Jy sigma=0.124894 ! Fit after self-cal, rms=0.012255Jy sigma=0.106204 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.000310407 Jy ! Component: 100 - total flux cleaned = 0.00116193 Jy ! Component: 150 - total flux cleaned = 0.0024692 Jy ! Component: 200 - total flux cleaned = 0.00226138 Jy ! Total flux subtracted in 200 components = 0.00226138 Jy ! Clean residual min=-0.001794 max=0.001778 Jy/beam ! Clean residual mean=0.000011 rms=0.000626 Jy/beam ! Combined flux in latest and established models = 1.08805 Jy ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 279 components and 1.08805 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.009104Jy sigma=0.079531 ! Fit after self-cal, rms=0.009023Jy sigma=0.078686 ! Inverting map selfcal true, true, soltime2 ! Performing amp+phase self-cal over 15 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.009023Jy sigma=0.078686 ! Fit after self-cal, rms=0.007972Jy sigma=0.068599 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.7571 mas, bmaj=8.707 mas, bpa=-13.31 degrees ! Estimated noise=8.40541 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000322979 Jy ! Component: 100 - total flux cleaned = -0.000539781 Jy ! Component: 150 - total flux cleaned = -0.000470954 Jy ! Component: 200 - total flux cleaned = -0.000573923 Jy ! Total flux subtracted in 200 components = -0.000573923 Jy ! Clean residual min=-0.001895 max=0.001920 Jy/beam ! Clean residual mean=-0.000008 rms=0.000515 Jy/beam ! Combined flux in latest and established models = 1.08747 Jy selfcal ! Performing phase self-cal ! Adding 55 model components to the UV plane model. ! The established model now contains 324 components and 1.08747 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.006009Jy sigma=0.054313 ! Fit after self-cal, rms=0.005859Jy sigma=0.051146 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.00055192 Jy ! Component: 100 - total flux cleaned = -0.000849271 Jy ! Component: 150 - total flux cleaned = -0.0017255 Jy ! Component: 200 - total flux cleaned = -0.00209318 Jy ! Total flux subtracted in 200 components = -0.00209318 Jy ! Clean residual min=-0.000926 max=0.000930 Jy/beam ! Clean residual mean=-0.000011 rms=0.000317 Jy/beam ! Combined flux in latest and established models = 1.08538 Jy ! Performing phase self-cal ! Adding 51 model components to the UV plane model. ! The established model now contains 366 components and 1.08538 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.004404Jy sigma=0.038778 ! Fit after self-cal, rms=0.004365Jy sigma=0.038465 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime2 ! Performing amp+phase self-cal over 15 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.004365Jy sigma=0.038465 ! Fit after self-cal, rms=0.003649Jy sigma=0.032861 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.785 mas, bmaj=10.41 mas, bpa=-3.1 degrees ! Estimated noise=9.34722 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00032115 Jy ! Component: 100 - total flux cleaned = -0.000411826 Jy ! Component: 150 - total flux cleaned = -0.000410627 Jy ! Component: 200 - total flux cleaned = -0.000373722 Jy ! Total flux subtracted in 200 components = -0.000373722 Jy ! Clean residual min=-0.008545 max=0.008036 Jy/beam ! Clean residual mean=-0.000015 rms=0.001474 Jy/beam ! Combined flux in latest and established models = 1.085 Jy selfcal ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 391 components and 1.085 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.005146Jy sigma=0.069979 ! Fit after self-cal, rms=0.004069Jy sigma=0.039576 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 = -7.82711e-05 Jy ! Component: 100 - total flux cleaned = -9.46264e-05 Jy ! Component: 150 - total flux cleaned = -0.000105386 Jy ! Component: 200 - total flux cleaned = -0.000137741 Jy ! Total flux subtracted in 200 components = -0.000137741 Jy ! Clean residual min=-0.003764 max=0.004030 Jy/beam ! Clean residual mean=0.000012 rms=0.000696 Jy/beam ! Combined flux in latest and established models = 1.08487 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 414 components and 1.08487 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.004412Jy sigma=0.046226 ! Fit after self-cal, rms=0.004182Jy sigma=0.040217 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "=== Finished amplitude self-cal for VLBA stations ===" ! === Finished amplitude self-cal for VLBA stations === unflag * ! Unflagging all baselines in the currently selected channels. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. integer lun lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_selfant_process.csh /tmp/ff.ff" ![@/tmp/ff.ff] selfant 1:SC,true ! 1:SC self-cal status=fixed weight=1 selfant 2:FD,true ! 2:FD self-cal status=fixed weight=1 selfant 2:LA,true ! 2:LA self-cal status=fixed weight=1 selfant 3:FD,true ! 3:FD self-cal status=fixed weight=1 selfant 3:PT,true ! 3:PT self-cal status=fixed weight=1 selfant 4:KP,true ! 4:KP self-cal status=fixed weight=1 selfant 4:MK,true ! 4:MK self-cal status=fixed weight=1 selfant 4:OV,true ! 4:OV self-cal status=fixed weight=1 selfant 5:FD,true ! 5:FD self-cal status=fixed weight=1 selfant 5:LA,true ! 5:LA self-cal status=fixed weight=1 selfant 5:MK,true ! 5:MK self-cal status=fixed weight=1 selfant 5:OV,true ! 5:OV self-cal status=fixed weight=1 selfant 5:PT,true ! 5:PT self-cal status=fixed weight=1 ![Exited script file: /tmp/ff.ff] selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 5 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 5 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 5 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 5 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=1.194427Jy sigma=32.246120 ! Fit after self-cal, rms=1.194427Jy sigma=32.246120 selfant "",false ! 1:FT self-cal status=correctable weight=1 ! 1:SC self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 3:FD self-cal status=correctable weight=1 ! 3:PT self-cal status=correctable weight=1 ! 4:KP self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:MK self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 5 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 5 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 5 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 5 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 5 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 5 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 5 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 5 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=1.194427Jy sigma=32.246120 ! Fit after self-cal, rms=0.003654Jy sigma=0.034332 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=0.6583 mas, bmaj=8.044 mas, bpa=-13.51 degrees ! Estimated noise=10.7501 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 = -4.61699e-05 Jy ! Component: 100 - total flux cleaned = -0.000193549 Jy ! Component: 150 - total flux cleaned = 0.000286414 Jy ! Component: 200 - total flux cleaned = 0.00058291 Jy ! Total flux subtracted in 200 components = 0.00058291 Jy ! Clean residual min=-0.000604 max=0.000607 Jy/beam ! Clean residual mean=0.000005 rms=0.000245 Jy/beam ! Combined flux in latest and established models = 1.08545 Jy ! Performing phase self-cal ! Adding 44 model components to the UV plane model. ! The established model now contains 456 components and 1.08545 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.002231Jy sigma=0.022128 ! Fit after self-cal, rms=0.002182Jy sigma=0.021452 ! Inverting map 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=0.7569 mas, bmaj=8.706 mas, bpa=-13.3 degrees ! Estimated noise=8.40088 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=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 7.79532e-06 Jy ! Component: 100 - total flux cleaned = 7.85552e-06 Jy ! Component: 150 - total flux cleaned = -2.34915e-05 Jy ! Component: 200 - total flux cleaned = -6.76265e-05 Jy ! Component: 250 - total flux cleaned = -5.52936e-05 Jy ! Component: 300 - total flux cleaned = -4.40412e-05 Jy ! Component: 350 - total flux cleaned = -1.28808e-05 Jy ! Component: 400 - total flux cleaned = 2.79861e-05 Jy ! Total flux subtracted in 400 components = 2.79861e-05 Jy ! Clean residual min=-0.008568 max=0.008337 Jy/beam ! Clean residual mean=-0.000004 rms=0.001657 Jy/beam ! Combined flux in latest and established models = 1.08548 Jy ! Performing phase self-cal ! Adding 66 model components to the UV plane model. ! The established model now contains 509 components and 1.08548 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.005700Jy sigma=0.088989 ! Fit after self-cal, rms=0.003813Jy sigma=0.040026 ! Inverting map dynam = dynam - 0.5 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) ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000297071 Jy ! Component: 100 - total flux cleaned = 0.000343468 Jy ! Component: 150 - total flux cleaned = 0.000275181 Jy ! Component: 200 - total flux cleaned = 0.000262712 Jy ! Component: 250 - total flux cleaned = 0.000226105 Jy ! Component: 300 - total flux cleaned = 0.000159258 Jy ! Component: 350 - total flux cleaned = 6.95236e-05 Jy ! Component: 400 - total flux cleaned = -3.9759e-05 Jy ! Total flux subtracted in 400 components = -3.9759e-05 Jy ! Clean residual min=-0.022636 max=0.023242 Jy/beam ! Clean residual mean=0.000009 rms=0.004498 Jy/beam ! Combined flux in latest and established models = 1.08544 Jy ! Performing phase self-cal ! Adding 55 model components to the UV plane model. ! The established model now contains 539 components and 1.08544 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.015044Jy sigma=0.239650 ! Fit after self-cal, rms=0.008117Jy sigma=0.092940 ! Inverting map print "=== Finished natural weighting clean for al sites ===" ! === Finished natural weighting clean for al sites === unflag * ! Unflagging all baselines in the currently selected channels. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 5 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 5 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 5 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 5 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=1.196970Jy sigma=32.325427 ! Fit after self-cal, rms=0.007465Jy sigma=0.080690 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.7569 mas, bmaj=8.706 mas, bpa=-13.31 degrees ! Estimated noise=8.39741 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000128768 Jy ! Component: 100 - total flux cleaned = 0.000236039 Jy ! Component: 150 - total flux cleaned = 0.000319599 Jy ! Component: 200 - total flux cleaned = 0.000509178 Jy ! Component: 250 - total flux cleaned = 0.000530694 Jy ! Component: 300 - total flux cleaned = 0.000613508 Jy ! Component: 350 - total flux cleaned = 0.000761055 Jy ! Component: 400 - total flux cleaned = 0.000945716 Jy ! Total flux subtracted in 400 components = 0.000945716 Jy ! Clean residual min=-0.039701 max=0.038224 Jy/beam ! Clean residual mean=-0.000014 rms=0.007576 Jy/beam ! Combined flux in latest and established models = 1.08638 Jy selfcal ! Performing phase self-cal ! Adding 60 model components to the UV plane model. ! The established model now contains 567 components and 1.08638 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.025154Jy sigma=0.400344 ! Fit after self-cal, rms=0.013040Jy sigma=0.147115 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 (24.6, -45.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000769884 Jy ! Component: 100 - total flux cleaned = 0.000919316 Jy ! Component: 150 - total flux cleaned = 0.00102422 Jy ! Component: 200 - total flux cleaned = 0.000948099 Jy ! Component: 250 - total flux cleaned = 0.000947327 Jy ! Component: 300 - total flux cleaned = 0.000943479 Jy ! Component: 350 - total flux cleaned = 0.000939749 Jy ! Component: 400 - total flux cleaned = 0.000897285 Jy ! Total flux subtracted in 400 components = 0.000897285 Jy ! Clean residual min=-0.002590 max=0.002347 Jy/beam ! Clean residual mean=0.000006 rms=0.000656 Jy/beam ! Combined flux in latest and established models = 1.08728 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 609 components and 1.08728 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.006033Jy sigma=0.071565 ! Fit after self-cal, rms=0.005841Jy sigma=0.065866 ! Inverting map ! Added new window around map position (-23.8, 20.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000347471 Jy ! Component: 100 - total flux cleaned = 0.000286069 Jy ! Component: 150 - total flux cleaned = 9.78309e-05 Jy ! Component: 200 - total flux cleaned = -8.72147e-05 Jy ! Component: 250 - total flux cleaned = -0.000207438 Jy ! Component: 300 - total flux cleaned = -0.000283147 Jy ! Component: 350 - total flux cleaned = -0.000326285 Jy ! Component: 400 - total flux cleaned = -0.000325767 Jy ! Total flux subtracted in 400 components = -0.000325767 Jy ! Clean residual min=-0.000948 max=0.000975 Jy/beam ! Clean residual mean=0.000002 rms=0.000241 Jy/beam ! Combined flux in latest and established models = 1.08696 Jy ! Performing phase self-cal ! Adding 55 model components to the UV plane model. ! The established model now contains 664 components and 1.08696 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.002881Jy sigma=0.036230 ! Fit after self-cal, rms=0.002730Jy sigma=0.034756 ! Inverting map selfcal true, true, soltime2 ! Performing amp+phase self-cal over 15 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.002730Jy sigma=0.034756 ! Fit after self-cal, rms=0.001971Jy sigma=0.020947 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.7568 mas, bmaj=8.708 mas, bpa=-13.31 degrees ! Estimated noise=8.39575 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 7.06103e-05 Jy ! Component: 100 - total flux cleaned = 5.79574e-05 Jy ! Component: 150 - total flux cleaned = 4.62346e-05 Jy ! Component: 200 - total flux cleaned = 3.64773e-05 Jy ! Component: 250 - total flux cleaned = 3.66513e-05 Jy ! Component: 300 - total flux cleaned = 4.45404e-05 Jy ! Component: 350 - total flux cleaned = 4.45014e-05 Jy ! Component: 400 - total flux cleaned = 5.16554e-05 Jy ! Total flux subtracted in 400 components = 5.16554e-05 Jy ! Clean residual min=-0.000438 max=0.000448 Jy/beam ! Clean residual mean=0.000000 rms=0.000106 Jy/beam ! Combined flux in latest and established models = 1.08701 Jy selfcal ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 708 components and 1.08701 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.001658Jy sigma=0.018146 ! Fit after self-cal, rms=0.001658Jy sigma=0.017963 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 (36.8, -4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 1.96255e-05 Jy ! Component: 100 - total flux cleaned = 7.49542e-06 Jy ! Component: 150 - total flux cleaned = -3.06055e-06 Jy ! Component: 200 - total flux cleaned = -3.31887e-06 Jy ! Component: 250 - total flux cleaned = 5.05318e-06 Jy ! Component: 300 - total flux cleaned = 1.29709e-05 Jy ! Component: 350 - total flux cleaned = 5.77541e-06 Jy ! Component: 400 - total flux cleaned = -8.05591e-06 Jy ! Total flux subtracted in 400 components = -8.05591e-06 Jy ! Clean residual min=-0.000290 max=0.000303 Jy/beam ! Clean residual mean=0.000001 rms=0.000074 Jy/beam ! Combined flux in latest and established models = 1.087 Jy ! Performing phase self-cal ! Adding 63 model components to the UV plane model. ! The established model now contains 768 components and 1.087 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.001204Jy sigma=0.014430 ! Fit after self-cal, rms=0.001181Jy sigma=0.014010 ! Inverting map ! Added new window around map position (-51.2, -29.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -4.05387e-07 Jy ! Component: 100 - total flux cleaned = 1.71044e-05 Jy ! Component: 150 - total flux cleaned = 8.30065e-06 Jy ! Component: 200 - total flux cleaned = -1.24292e-05 Jy ! Component: 250 - total flux cleaned = -6.84758e-06 Jy ! Component: 300 - total flux cleaned = -1.57328e-06 Jy ! Component: 350 - total flux cleaned = 2.77862e-05 Jy ! Component: 400 - total flux cleaned = 9.9487e-05 Jy ! Total flux subtracted in 400 components = 9.9487e-05 Jy ! Clean residual min=-0.000258 max=0.000205 Jy/beam ! Clean residual mean=0.000001 rms=0.000047 Jy/beam ! Combined flux in latest and established models = 1.0871 Jy ! Performing phase self-cal ! Adding 39 model components to the UV plane model. ! The established model now contains 805 components and 1.0871 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.000827Jy sigma=0.009685 ! Fit after self-cal, rms=0.000822Jy sigma=0.009334 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime3 ! Performing amp+phase self-cal over 1 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.000822Jy sigma=0.009334 ! Fit after self-cal, rms=0.000783Jy sigma=0.009897 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.784 mas, bmaj=10.41 mas, bpa=-3.105 degrees ! Estimated noise=9.34414 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -3.65793e-06 Jy ! Component: 100 - total flux cleaned = 7.03353e-06 Jy ! Component: 150 - total flux cleaned = 2.99104e-05 Jy ! Component: 200 - total flux cleaned = 4.96632e-05 Jy ! Component: 250 - total flux cleaned = 6.70227e-05 Jy ! Component: 300 - total flux cleaned = 7.49764e-05 Jy ! Component: 350 - total flux cleaned = 8.39109e-05 Jy ! Component: 400 - total flux cleaned = 9.17767e-05 Jy ! Total flux subtracted in 400 components = 9.17767e-05 Jy ! Clean residual min=-0.000057 max=0.000035 Jy/beam ! Clean residual mean=0.000000 rms=0.000009 Jy/beam ! Combined flux in latest and established models = 1.08719 Jy selfcal ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 829 components and 1.08719 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.000745Jy sigma=0.008302 ! Fit after self-cal, rms=0.000745Jy sigma=0.008282 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=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 3.2874e-06 Jy ! Component: 100 - total flux cleaned = 7.04535e-06 Jy ! Component: 150 - total flux cleaned = 1.11138e-05 Jy ! Component: 200 - total flux cleaned = 7.24306e-06 Jy ! Component: 250 - total flux cleaned = 1.59816e-06 Jy ! Component: 300 - total flux cleaned = -4.99783e-06 Jy ! Component: 350 - total flux cleaned = -1.20471e-05 Jy ! Component: 400 - total flux cleaned = -1.9633e-05 Jy ! Total flux subtracted in 400 components = -1.9633e-05 Jy ! Clean residual min=-0.000022 max=0.000053 Jy/beam ! Clean residual mean=0.000002 rms=0.000008 Jy/beam ! Combined flux in latest and established models = 1.08717 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 848 components and 1.08717 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.000723Jy sigma=0.006387 ! Fit after self-cal, rms=0.000722Jy sigma=0.006372 ! Inverting map ! Added new window around map position (50, -48.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 3.00935e-06 Jy ! Component: 100 - total flux cleaned = 5.15206e-06 Jy ! Component: 150 - total flux cleaned = 5.28645e-06 Jy ! Component: 200 - total flux cleaned = 5.40571e-06 Jy ! Component: 250 - total flux cleaned = 5.78577e-06 Jy ! Component: 300 - total flux cleaned = 6.26649e-06 Jy ! Component: 350 - total flux cleaned = 6.98718e-06 Jy ! Component: 400 - total flux cleaned = 7.98138e-06 Jy ! Total flux subtracted in 400 components = 7.98138e-06 Jy ! Clean residual min=-0.000869 max=0.000383 Jy/beam ! Clean residual mean=-0.000002 rms=0.000131 Jy/beam ! Combined flux in latest and established models = 1.08718 Jy ! Performing phase self-cal ! Adding 90 model components to the UV plane model. ! The established model now contains 929 components and 1.08718 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.000706Jy sigma=0.005808 ! Fit after self-cal, rms=0.000686Jy sigma=0.004760 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "=== Finished amplitude self-cal for all stations ===" ! === Finished amplitude self-cal for all stations === restore ! Inverting map and beam ! Estimated beam: bmin=0.7568 mas, bmaj=8.708 mas, bpa=-13.31 degrees ! Estimated noise=8.39536 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.7568 x 8.708 at -13.31 degrees (North through East) ! Clean map min=-0.037619 max=0.94301 Jy/beam if(peak(flux) > thresh) selfcal true, true clean clean_niter,clean_gain selfcal flux_cutoff = imstat(rms) * dynam repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) selfcal true, true end if ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.000686Jy sigma=0.004760 ! Fit after self-cal, rms=0.000684Jy sigma=0.007091 ! Inverting map and beam ! Estimated beam: bmin=0.7568 mas, bmaj=8.707 mas, bpa=-13.31 degrees ! Estimated noise=8.39514 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 4.81698e-05 Jy ! Component: 100 - total flux cleaned = 9.94061e-05 Jy ! Component: 150 - total flux cleaned = 0.000120211 Jy ! Component: 200 - total flux cleaned = 0.000120293 Jy ! Component: 250 - total flux cleaned = 0.000126236 Jy ! Component: 300 - total flux cleaned = 0.000118366 Jy ! Component: 350 - total flux cleaned = 0.000118023 Jy ! Component: 400 - total flux cleaned = 0.000106744 Jy ! Total flux subtracted in 400 components = 0.000106744 Jy ! Clean residual min=-0.000316 max=0.000290 Jy/beam ! Clean residual mean=-0.000000 rms=0.000064 Jy/beam ! Combined flux in latest and established models = 1.08728 Jy ! Performing phase self-cal ! Adding 55 model components to the UV plane model. ! The established model now contains 973 components and 1.08728 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.000499Jy sigma=0.003648 ! Fit after self-cal, rms=0.000483Jy sigma=0.001394 ! Inverting map ! Added new window around map position (-41, -3.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 3.16518e-05 Jy ! Component: 100 - total flux cleaned = 0.000107276 Jy ! Component: 150 - total flux cleaned = 0.000189089 Jy ! Component: 200 - total flux cleaned = 0.000268597 Jy ! Component: 250 - total flux cleaned = 0.000344218 Jy ! Component: 300 - total flux cleaned = 0.000416505 Jy ! Component: 350 - total flux cleaned = 0.000485972 Jy ! Component: 400 - total flux cleaned = 0.0005529 Jy ! Total flux subtracted in 400 components = 0.0005529 Jy ! Clean residual min=-0.000058 max=0.000057 Jy/beam ! Clean residual mean=0.000003 rms=0.000018 Jy/beam ! Combined flux in latest and established models = 1.08784 Jy ! Performing phase self-cal ! Adding 54 model components to the UV plane model. ! The established model now contains 1021 components and 1.08784 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.000330Jy sigma=0.003771 ! Fit after self-cal, rms=0.000329Jy sigma=0.003689 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.000329Jy sigma=0.003689 ! Fit after self-cal, rms=0.000257Jy sigma=0.003986 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=0.6582 mas, bmaj=8.046 mas, bpa=-13.51 degrees ! Estimated noise=10.7426 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0885365 Jy ! Component: 100 - total flux cleaned = 0.0662195 Jy ! Component: 150 - total flux cleaned = 0.131335 Jy ! Component: 200 - total flux cleaned = 0.204835 Jy ! Component: 250 - total flux cleaned = 0.237026 Jy ! Component: 300 - total flux cleaned = 0.23683 Jy ! Component: 350 - total flux cleaned = 0.226163 Jy ! Component: 400 - total flux cleaned = 0.226247 Jy ! Total flux subtracted in 400 components = 0.226247 Jy ! Clean residual min=-0.242505 max=0.202599 Jy/beam ! Clean residual mean=0.002429 rms=0.042688 Jy/beam ! Combined flux in latest and established models = 0.226247 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 84 model components to the UV plane model. ! The established model now contains 84 components and 0.226247 Jy ! Inverting map and beam ! Estimated beam: bmin=0.7568 mas, bmaj=8.707 mas, bpa=-13.31 degrees ! Estimated noise=8.3951 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 (0.8, 3.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.120364 Jy ! Component: 100 - total flux cleaned = 0.126783 Jy ! Component: 150 - total flux cleaned = 0.126709 Jy ! Component: 200 - total flux cleaned = 0.129515 Jy ! Component: 250 - total flux cleaned = 0.13254 Jy ! Component: 300 - total flux cleaned = 0.132407 Jy ! Component: 350 - total flux cleaned = 0.132486 Jy ! Component: 400 - total flux cleaned = 0.137989 Jy ! Total flux subtracted in 400 components = 0.137989 Jy ! Clean residual min=-0.112151 max=0.181280 Jy/beam ! Clean residual mean=-0.000435 rms=0.034001 Jy/beam ! Combined flux in latest and established models = 0.364236 Jy ! Adding 72 model components to the UV plane model. ! The established model now contains 155 components and 0.364236 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.784 mas, bmaj=10.41 mas, bpa=-3.104 degrees ! Estimated noise=9.34382 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0622941 Jy ! Component: 100 - total flux cleaned = -0.113687 Jy ! Component: 150 - total flux cleaned = -0.160894 Jy ! Component: 200 - total flux cleaned = -0.206691 Jy ! Component: 250 - total flux cleaned = -0.249875 Jy ! Component: 300 - total flux cleaned = -0.291055 Jy ! Component: 350 - total flux cleaned = -0.330889 Jy ! Component: 400 - total flux cleaned = -0.371921 Jy ! Total flux subtracted in 400 components = -0.371921 Jy ! Clean residual min=-0.040519 max=0.027113 Jy/beam ! Clean residual mean=-0.007605 rms=0.015015 Jy/beam ! Combined flux in latest and established models = -0.00768527 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=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0438478 Jy ! Component: 100 - total flux cleaned = -0.0909094 Jy ! Component: 150 - total flux cleaned = -0.141549 Jy ! Component: 200 - total flux cleaned = -0.196211 Jy ! Component: 250 - total flux cleaned = -0.255277 Jy ! Component: 300 - total flux cleaned = -0.319171 Jy ! Component: 350 - total flux cleaned = -0.388338 Jy ! Component: 400 - total flux cleaned = -0.463243 Jy ! Total flux subtracted in 400 components = -0.463243 Jy ! Clean residual min=-0.081017 max=0.048480 Jy/beam ! Clean residual mean=-0.015013 rms=0.028576 Jy/beam ! Combined flux in latest and established models = -0.470928 Jy ! Adding 171 model components to the UV plane model. ! The established model now contains 325 components and -0.470928 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=0.7568 mas, bmaj=8.707 mas, bpa=-13.31 degrees ! Estimated noise=8.3951 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0408329 Jy ! Component: 100 - total flux cleaned = 0.0484495 Jy ! Component: 150 - total flux cleaned = 0.049526 Jy ! Component: 200 - total flux cleaned = 0.0477069 Jy ! Component: 250 - total flux cleaned = 0.0441649 Jy ! Component: 300 - total flux cleaned = 0.0432402 Jy ! Component: 350 - total flux cleaned = 0.0443829 Jy ! Component: 400 - total flux cleaned = 0.0444685 Jy ! Total flux subtracted in 400 components = 0.0444685 Jy ! Clean residual min=-0.030917 max=0.036804 Jy/beam ! Clean residual mean=-0.000006 rms=0.008867 Jy/beam ! Combined flux in latest and established models = -0.42646 Jy selfcal ! Performing phase self-cal ! Adding 63 model components to the UV plane model. ! The established model now contains 386 components and -0.42646 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.085989Jy sigma=1.099310 ! Fit after self-cal, rms=0.083513Jy sigma=1.056276 wmodel /scr/pima/rdv93_uvs/J1107-4449_X_map.mod ! Writing 386 model components to file: /scr/pima/rdv93_uvs/J1107-4449_X_map.mod wobs /scr/pima/rdv93_uvs/J1107-4449_X_uvs.fits ! Applying 77 buffered edits. ! Writing UV FITS file: /scr/pima/rdv93_uvs/J1107-4449_X_uvs.fits wwins /scr/pima/rdv93_uvs/J1107-4449_X_map.win ! wwins: Wrote 7 windows to /scr/pima/rdv93_uvs/J1107-4449_X_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 = 3.4716e-05 Jy ! Component: 100 - total flux cleaned = 3.45074e-05 Jy ! Component: 150 - total flux cleaned = 0.000392687 Jy ! Component: 200 - total flux cleaned = 0.000395539 Jy ! Component: 250 - total flux cleaned = 5.54554e-05 Jy ! Component: 300 - total flux cleaned = -0.000942567 Jy ! Component: 350 - total flux cleaned = -0.00464602 Jy ! Component: 400 - total flux cleaned = -0.00740235 Jy ! Component: 450 - total flux cleaned = -0.0102371 Jy ! Component: 500 - total flux cleaned = -0.0109609 Jy ! Component: 550 - total flux cleaned = -0.012114 Jy ! Component: 600 - total flux cleaned = -0.0129288 Jy ! Component: 650 - total flux cleaned = -0.0133757 Jy ! Component: 700 - total flux cleaned = -0.0133917 Jy ! Component: 750 - total flux cleaned = -0.0138708 Jy ! Component: 800 - total flux cleaned = -0.0143919 Jy ! Component: 850 - total flux cleaned = -0.0143644 Jy ! Component: 900 - total flux cleaned = -0.0143831 Jy ! Component: 950 - total flux cleaned = -0.0122586 Jy ! Component: 1000 - total flux cleaned = -0.0122772 Jy ! Component: 1050 - total flux cleaned = -0.0122974 Jy ! Component: 1100 - total flux cleaned = -0.0123327 Jy ! Component: 1150 - total flux cleaned = -0.0124262 Jy ! Component: 1200 - total flux cleaned = -0.0125619 Jy ! Component: 1250 - total flux cleaned = -0.0126712 Jy ! Component: 1300 - total flux cleaned = -0.0127543 Jy ! Component: 1350 - total flux cleaned = -0.0128646 Jy ! Component: 1400 - total flux cleaned = -0.0129837 Jy ! Component: 1450 - total flux cleaned = -0.0130661 Jy ! Component: 1500 - total flux cleaned = -0.0131904 Jy ! Component: 1550 - total flux cleaned = -0.0131788 Jy ! Component: 1600 - total flux cleaned = -0.0131178 Jy ! Component: 1650 - total flux cleaned = -0.013116 Jy ! Component: 1700 - total flux cleaned = -0.013248 Jy ! Component: 1750 - total flux cleaned = -0.0135865 Jy ! Component: 1800 - total flux cleaned = -0.0139599 Jy ! Component: 1850 - total flux cleaned = -0.0145325 Jy ! Component: 1900 - total flux cleaned = -0.015362 Jy ! Component: 1950 - total flux cleaned = -0.0162047 Jy ! Component: 2000 - total flux cleaned = -0.0173174 Jy ! Component: 2050 - total flux cleaned = -0.018547 Jy ! Component: 2100 - total flux cleaned = -0.019916 Jy ! Component: 2150 - total flux cleaned = -0.0218553 Jy ! Component: 2200 - total flux cleaned = -0.0240336 Jy ! Component: 2250 - total flux cleaned = -0.0264539 Jy ! Component: 2300 - total flux cleaned = -0.0296649 Jy ! Component: 2350 - total flux cleaned = -0.0337478 Jy ! Component: 2400 - total flux cleaned = -0.0377847 Jy ! Component: 2450 - total flux cleaned = -0.0434124 Jy ! Component: 2500 - total flux cleaned = -0.0491596 Jy ! Component: 2550 - total flux cleaned = -0.0562998 Jy ! Component: 2600 - total flux cleaned = -0.0659051 Jy ! Component: 2650 - total flux cleaned = -0.0753188 Jy ! Component: 2700 - total flux cleaned = -0.0865705 Jy ! Component: 2750 - total flux cleaned = -0.103256 Jy ! Component: 2800 - total flux cleaned = -0.11891 Jy ! Component: 2850 - total flux cleaned = -0.139758 Jy ! Component: 2900 - total flux cleaned = -0.165716 Jy ! Component: 2950 - total flux cleaned = -0.192391 Jy ! Component: 3000 - total flux cleaned = -0.227096 Jy ! Component: 3050 - total flux cleaned = -0.271621 Jy ! Component: 3100 - total flux cleaned = -0.316326 Jy ! Component: 3150 - total flux cleaned = -0.379215 Jy ! Component: 3200 - total flux cleaned = -0.457901 Jy ! Component: 3250 - total flux cleaned = -0.53179 Jy ! Component: 3300 - total flux cleaned = -0.644342 Jy ! Component: 3350 - total flux cleaned = -0.769364 Jy ! Component: 3400 - total flux cleaned = -0.914035 Jy ! Component: 3450 - total flux cleaned = -1.08345 Jy ! Component: 3500 - total flux cleaned = -1.29716 Jy ! Component: 3550 - total flux cleaned = -1.51949 Jy ! Component: 3600 - total flux cleaned = -1.83967 Jy ! Component: 3650 - total flux cleaned = -2.17702 Jy ! Component: 3700 - total flux cleaned = -2.58771 Jy ! Component: 3750 - total flux cleaned = -3.06782 Jy ! Component: 3800 - total flux cleaned = -3.65744 Jy ! Component: 3850 - total flux cleaned = -4.34552 Jy ! Component: 3900 - total flux cleaned = -5.23024 Jy ! Component: 3950 - total flux cleaned = -6.17768 Jy ! Component: 4000 - total flux cleaned = -7.33305 Jy ! Component: 4050 - total flux cleaned = -8.8851 Jy ! Total flux subtracted in 4096 components = -10.463 Jy ! Clean residual min=-2609.897949 max=2194.493652 Jy/beam ! Clean residual mean=-5.535024 rms=456.689758 Jy/beam ! Combined flux in latest and established models = -10.8894 Jy keep ! Adding 77 model components to the UV plane model. ! The established model now contains 461 components and -10.8895 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv93_uvs/J1107-4449_X\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.7568 x 8.707 at -13.31 degrees (North through East) ! Clean map min=-1276 max=1280.1 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv93_uvs/J1107-4449_X_map.fits float peak_flux peak_flux = peak(flux,max) invert ! Inverting map float image_peak image_peak = int(10000*peak_flux)/10000 float image_rms image_rms = int(100000*imstat(rms))/100000 print "=========================================" ! ========================================= print "The image peak= ", image_peak, " rms= ", image_rms, " Jy" ! The image peak= 1280.1 rms= 245.929 Jy quit ! Quitting program ! Log file /scr/pima/rdv93_uvs/J1107-4449_X_uvs.log closed on Tue Oct 8 15:46:43 2013