! Started logfile: /scr/pima/rdv62_uvs/J1924-2914_S_uvs.log on Fri Nov 22 22:40:19 2013 obs /scr/pima/rdv62_uvs/J1924-2914_S_uva.fits ! Reading UV FITS file: /scr/pima/rdv62_uvs/J1924-2914_S_uva.fits ! AN table 1: 4 integrations on 21 of 21 possible baselines. ! AN table 2: 11 integrations on 36 of 36 possible baselines. ! Apparent sampling: 0.441667 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1924-2914 ! ! 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.23299e+09 8e+06 1 8e+06 ! 02 2 2.26299e+09 8e+06 1 8e+06 ! 03 3 2.35299e+09 8e+06 1 8e+06 ! 04 4 2.37299e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 848 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.4 print "taper_size =",taper_size ! taper_size = 20 ![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv62_uvs/J1924-2914_S] 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:BR self-cal status=correctable weight=1 ! 1:KK self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:MK self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:OV self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:KK self-cal status=correctable weight=1 ! 2:KP self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:MK self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:OV self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 system "rm /tmp/ff.ff" mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.400x0.400 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 1 telescope corrections were flagged in sub-array 1. ! A total of 49 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 2. ! A total of 1 telescope corrections were flagged in sub-array 1. ! A total of 49 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 3. ! A total of 1 telescope corrections were flagged in sub-array 1. ! A total of 49 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 4. ! A total of 1 telescope corrections were flagged in sub-array 1. ! A total of 49 telescope corrections were flagged in sub-array 2. ! ! Fit before self-cal, rms=4.284729Jy sigma=181.946399 ! Fit after self-cal, rms=3.375378Jy sigma=150.432021 ! 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:KK ! Flagging baselines of antenna 1:KK in the currently selected channels. flag 2:KK ! Flagging baselines of antenna 2:KK 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=3.808 mas, bmaj=15.12 mas, bpa=-4.71 degrees ! Estimated noise=1.33463 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 = 3.46987 Jy ! Component: 100 - total flux cleaned = 4.60315 Jy ! Total flux subtracted in 100 components = 4.60315 Jy ! Clean residual min=-0.238322 max=0.474369 Jy/beam ! Clean residual mean=0.005297 rms=0.114378 Jy/beam ! Combined flux in latest and established models = 4.60315 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 31 components and 4.60315 Jy ! ! Correcting IF 1. ! A total of 3 telescope corrections were flagged in sub-array 1. ! A total of 3 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 2. ! A total of 3 telescope corrections were flagged in sub-array 1. ! A total of 3 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 3. ! A total of 3 telescope corrections were flagged in sub-array 1. ! A total of 3 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 4. ! A total of 3 telescope corrections were flagged in sub-array 1. ! A total of 3 telescope corrections were flagged in sub-array 2. ! ! Fit before self-cal, rms=1.158344Jy sigma=51.513222 ! Fit after self-cal, rms=1.158759Jy sigma=51.443815 ! Inverting map and beam ! Estimated beam: bmin=3.808 mas, bmaj=15.12 mas, bpa=-4.71 degrees ! Estimated noise=1.33463 mJy/beam. 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=5.025 mas, bmaj=17.14 mas, bpa=-3.244 degrees ! Estimated noise=0.828675 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.697358 Jy ! Component: 100 - total flux cleaned = 1.03918 Jy ! Component: 150 - total flux cleaned = 1.22551 Jy ! Component: 200 - total flux cleaned = 1.32843 Jy ! Total flux subtracted in 200 components = 1.32843 Jy ! Clean residual min=-0.186810 max=0.277817 Jy/beam ! Clean residual mean=0.011242 rms=0.069827 Jy/beam ! Combined flux in latest and established models = 5.93158 Jy ! Performing phase self-cal ! Adding 79 model components to the UV plane model. ! The established model now contains 107 components and 5.93158 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.861616Jy sigma=34.582891 ! Fit after self-cal, rms=0.860788Jy sigma=34.567358 ! 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: ! BR 0.99 KK 1.00* LA 1.12 MK 0.95 ! NL 1.21 OV 1.01 PT 1.04 ! ! A total of 18 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.01 FD 0.97 KK 1.00* KP 0.99 ! LA 1.10 MK 1.12 NL 1.14 OV 0.97 ! PT 1.03 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 1.03 KK 1.00* LA 1.05 MK 0.96 ! NL 1.22 OV 1.05 PT 1.04 ! ! A total of 18 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.04 FD 1.04 KK 1.00* KP 1.02 ! LA 1.02 MK 1.08 NL 1.17 OV 1.03 ! PT 1.02 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.09 KK 1.00* LA 0.83 MK 0.91 ! NL 1.29 OV 1.12 PT 1.03 ! ! A total of 18 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.06 FD 1.05 KK 1.00* KP 1.02 ! LA 0.82 MK 0.96 NL 1.26 OV 1.07 ! PT 1.04 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 1.06 KK 1.00* LA 0.73 MK 0.90 ! NL 1.14 OV 0.92 PT 0.98 ! ! A total of 18 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.94 FD 1.02 KK 1.00* KP 0.95 ! LA 0.72 MK 0.92 NL 1.16 OV 0.88 ! PT 1.00 ! ! ! Fit before self-cal, rms=0.860788Jy sigma=34.567358 ! Fit after self-cal, rms=0.290974Jy sigma=7.086917 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=5.21 mas, bmaj=16.92 mas, bpa=-2.873 degrees ! Estimated noise=0.848668 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.0152759 Jy ! Component: 100 - total flux cleaned = 0.0189047 Jy ! Component: 150 - total flux cleaned = 0.0198991 Jy ! Component: 200 - total flux cleaned = 0.019933 Jy ! Total flux subtracted in 200 components = 0.019933 Jy ! Clean residual min=-0.059377 max=0.084988 Jy/beam ! Clean residual mean=0.000906 rms=0.017675 Jy/beam ! Combined flux in latest and established models = 5.95152 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 133 components and 5.95151 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.277881Jy sigma=6.827246 ! Fit after self-cal, rms=0.278873Jy sigma=6.791432 ! 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.278873Jy sigma=6.791432 ! Fit after self-cal, rms=0.279098Jy sigma=6.751905 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=5.204 mas, bmaj=16.88 mas, bpa=-2.928 degrees ! Estimated noise=0.846618 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00353922 Jy ! Component: 100 - total flux cleaned = 0.00443188 Jy ! Component: 150 - total flux cleaned = 0.0052099 Jy ! Component: 200 - total flux cleaned = 0.0052082 Jy ! Total flux subtracted in 200 components = 0.0052082 Jy ! Clean residual min=-0.058641 max=0.079036 Jy/beam ! Clean residual mean=0.000930 rms=0.016797 Jy/beam ! Combined flux in latest and established models = 5.95672 Jy selfcal ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 138 components and 5.95672 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.269468Jy sigma=6.604321 ! Fit after self-cal, rms=0.269256Jy sigma=6.581185 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (10.4, 6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0343708 Jy ! Component: 100 - total flux cleaned = 0.0304623 Jy ! Component: 150 - total flux cleaned = 0.0319038 Jy ! Component: 200 - total flux cleaned = 0.0344694 Jy ! Total flux subtracted in 200 components = 0.0344694 Jy ! Clean residual min=-0.055267 max=0.052310 Jy/beam ! Clean residual mean=0.000415 rms=0.013805 Jy/beam ! Combined flux in latest and established models = 5.99119 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 161 components and 5.99119 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.233246Jy sigma=5.786657 ! Fit after self-cal, rms=0.226816Jy sigma=5.667131 ! 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.226816Jy sigma=5.667131 ! Fit after self-cal, rms=0.225081Jy sigma=5.557695 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=5.197 mas, bmaj=16.85 mas, bpa=-2.994 degrees ! Estimated noise=0.845293 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0165241 Jy ! Component: 100 - total flux cleaned = 0.0140272 Jy ! Component: 150 - total flux cleaned = 0.0139531 Jy ! Component: 200 - total flux cleaned = 0.0148307 Jy ! Total flux subtracted in 200 components = 0.0148307 Jy ! Clean residual min=-0.050538 max=0.045066 Jy/beam ! Clean residual mean=0.000240 rms=0.011760 Jy/beam ! Combined flux in latest and established models = 6.00602 Jy selfcal ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 179 components and 6.00602 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.205984Jy sigma=5.188286 ! Fit after self-cal, rms=0.202731Jy sigma=5.132675 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.0054878 Jy ! Component: 100 - total flux cleaned = 0.00320277 Jy ! Component: 150 - total flux cleaned = 0.00253518 Jy ! Component: 200 - total flux cleaned = 0.00192022 Jy ! Total flux subtracted in 200 components = 0.00192022 Jy ! Clean residual min=-0.047951 max=0.042006 Jy/beam ! Clean residual mean=0.000203 rms=0.011311 Jy/beam ! Combined flux in latest and established models = 6.00794 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 201 components and 6.00794 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.193988Jy sigma=4.991889 ! Fit after self-cal, rms=0.192377Jy sigma=4.972578 ! 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 15 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.192377Jy sigma=4.972578 ! Fit after self-cal, rms=0.191251Jy sigma=4.894078 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=9.506 mas, bmaj=31.81 mas, bpa=19.27 degrees ! Estimated noise=0.800105 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00419243 Jy ! Component: 100 - total flux cleaned = 0.00419533 Jy ! Component: 150 - total flux cleaned = 0.00464929 Jy ! Component: 200 - total flux cleaned = 0.00504292 Jy ! Total flux subtracted in 200 components = 0.00504292 Jy ! Clean residual min=-0.029103 max=0.040568 Jy/beam ! Clean residual mean=0.000031 rms=0.010957 Jy/beam ! Combined flux in latest and established models = 6.01298 Jy selfcal ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 218 components and 6.01298 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.187548Jy sigma=4.749726 ! Fit after self-cal, rms=0.185374Jy sigma=4.724650 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.00614865 Jy ! Component: 100 - total flux cleaned = 0.00691938 Jy ! Component: 150 - total flux cleaned = 0.0072878 Jy ! Component: 200 - total flux cleaned = 0.0073126 Jy ! Total flux subtracted in 200 components = 0.0073126 Jy ! Clean residual min=-0.029425 max=0.039438 Jy/beam ! Clean residual mean=-0.000098 rms=0.010683 Jy/beam ! Combined flux in latest and established models = 6.0203 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 229 components and 6.0203 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.183407Jy sigma=4.651825 ! Fit after self-cal, rms=0.181940Jy sigma=4.638035 ! 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:BR,true ! 1:BR self-cal status=fixed weight=1 selfant 1:LA,true ! 1:LA self-cal status=fixed weight=1 selfant 1:MK,true ! 1:MK self-cal status=fixed weight=1 selfant 1:NL,true ! 1:NL self-cal status=fixed weight=1 selfant 1:OV,true ! 1:OV self-cal status=fixed weight=1 selfant 1:PT,true ! 1:PT self-cal status=fixed weight=1 selfant 2:BR,true ! 2:BR self-cal status=fixed weight=1 selfant 2:FD,true ! 2:FD self-cal status=fixed weight=1 selfant 2:KP,true ! 2:KP self-cal status=fixed weight=1 selfant 2:LA,true ! 2:LA self-cal status=fixed weight=1 selfant 2:MK,true ! 2:MK self-cal status=fixed weight=1 selfant 2:NL,true ! 2:NL self-cal status=fixed weight=1 selfant 2:OV,true ! 2:OV self-cal status=fixed weight=1 selfant 2:PT,true ! 2: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 1 telescope corrections were flagged in sub-array 1. ! A total of 8 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 2. ! A total of 1 telescope corrections were flagged in sub-array 1. ! A total of 8 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 3. ! A total of 1 telescope corrections were flagged in sub-array 1. ! A total of 8 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 4. ! A total of 1 telescope corrections were flagged in sub-array 1. ! A total of 8 telescope corrections were flagged in sub-array 2. ! ! Fit before self-cal, rms=0.333160Jy sigma=10.456695 ! Fit after self-cal, rms=0.333467Jy sigma=7.267241 selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:KK self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:MK self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:OV self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:KK self-cal status=correctable weight=1 ! 2:KP self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:MK self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:OV self-cal status=correctable weight=1 ! 2: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 59 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 2. ! A total of 59 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 3. ! A total of 59 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 4. ! A total of 59 telescope corrections were flagged in sub-array 2. ! ! Fit before self-cal, rms=0.333467Jy sigma=7.267241 ! Fit after self-cal, rms=0.282785Jy sigma=7.018780 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.235 mas, bmaj=15.31 mas, bpa=-6.877 degrees ! Estimated noise=1.43194 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.00379695 Jy ! Component: 100 - total flux cleaned = 0.00224801 Jy ! Component: 150 - total flux cleaned = -0.00274684 Jy ! Component: 200 - total flux cleaned = -0.00421436 Jy ! Total flux subtracted in 200 components = -0.00421436 Jy ! Clean residual min=-0.069421 max=0.060454 Jy/beam ! Clean residual mean=-0.000083 rms=0.013234 Jy/beam ! Combined flux in latest and established models = 6.01608 Jy ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 265 components and 6.01608 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.260817Jy sigma=6.848911 ! Fit after self-cal, rms=0.258567Jy sigma=6.797786 ! 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=4.463 mas, bmaj=17.12 mas, bpa=-5.144 degrees ! Estimated noise=0.843713 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 = 0.00127136 Jy ! Component: 100 - total flux cleaned = 0.00183818 Jy ! Component: 150 - total flux cleaned = 0.0013039 Jy ! Component: 200 - total flux cleaned = 0.00130288 Jy ! Component: 250 - total flux cleaned = 0.000811416 Jy ! Component: 300 - total flux cleaned = 0.000802526 Jy ! Component: 350 - total flux cleaned = 0.000341897 Jy ! Component: 400 - total flux cleaned = -8.53657e-05 Jy ! Total flux subtracted in 400 components = -8.53657e-05 Jy ! Clean residual min=-0.036669 max=0.029790 Jy/beam ! Clean residual mean=-0.000075 rms=0.008982 Jy/beam ! Combined flux in latest and established models = 6.016 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 299 components and 6.016 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.250768Jy sigma=6.731276 ! Fit after self-cal, rms=0.250155Jy sigma=6.725648 ! 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.00179821 Jy ! Component: 100 - total flux cleaned = 0.000536661 Jy ! Component: 150 - total flux cleaned = 0.000532869 Jy ! Component: 200 - total flux cleaned = 0.000136518 Jy ! Component: 250 - total flux cleaned = -0.000265522 Jy ! Component: 300 - total flux cleaned = -0.000616574 Jy ! Component: 350 - total flux cleaned = -0.000978306 Jy ! Component: 400 - total flux cleaned = -0.00134713 Jy ! Total flux subtracted in 400 components = -0.00134713 Jy ! Clean residual min=-0.034970 max=0.029087 Jy/beam ! Clean residual mean=-0.000059 rms=0.008825 Jy/beam ! Combined flux in latest and established models = 6.01465 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 312 components and 6.01465 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.245359Jy sigma=6.687746 ! Fit after self-cal, rms=0.244880Jy sigma=6.684751 ! 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 1 telescope corrections were flagged in sub-array 1. ! A total of 67 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 2. ! A total of 1 telescope corrections were flagged in sub-array 1. ! A total of 67 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 3. ! A total of 1 telescope corrections were flagged in sub-array 1. ! A total of 67 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 4. ! A total of 1 telescope corrections were flagged in sub-array 1. ! A total of 67 telescope corrections were flagged in sub-array 2. ! ! Fit before self-cal, rms=0.278381Jy sigma=6.922624 ! Fit after self-cal, rms=0.242080Jy sigma=6.649256 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=4.467 mas, bmaj=17.15 mas, bpa=-5.111 degrees ! Estimated noise=0.843115 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00333685 Jy ! Component: 100 - total flux cleaned = -0.00562414 Jy ! Component: 150 - total flux cleaned = -0.00520108 Jy ! Component: 200 - total flux cleaned = -0.00485223 Jy ! Component: 250 - total flux cleaned = -0.00485334 Jy ! Component: 300 - total flux cleaned = -0.00484552 Jy ! Component: 350 - total flux cleaned = -0.00514783 Jy ! Component: 400 - total flux cleaned = -0.00513867 Jy ! Total flux subtracted in 400 components = -0.00513867 Jy ! Clean residual min=-0.033174 max=0.027468 Jy/beam ! Clean residual mean=-0.000085 rms=0.008204 Jy/beam ! Combined flux in latest and established models = 6.00951 Jy selfcal ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 338 components and 6.00951 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.237728Jy sigma=6.604107 ! Fit after self-cal, rms=0.237149Jy sigma=6.599223 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 = 0.00180883 Jy ! Component: 100 - total flux cleaned = 0.00213546 Jy ! Component: 150 - total flux cleaned = 0.00212286 Jy ! Component: 200 - total flux cleaned = 0.0018105 Jy ! Component: 250 - total flux cleaned = 0.00121747 Jy ! Component: 300 - total flux cleaned = 0.00121912 Jy ! Component: 350 - total flux cleaned = 0.000920318 Jy ! Component: 400 - total flux cleaned = 0.000908963 Jy ! Total flux subtracted in 400 components = 0.000908963 Jy ! Clean residual min=-0.032740 max=0.027007 Jy/beam ! Clean residual mean=-0.000094 rms=0.008049 Jy/beam ! Combined flux in latest and established models = 6.01042 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 347 components and 6.01042 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.234329Jy sigma=6.570540 ! Fit after self-cal, rms=0.233792Jy sigma=6.567631 ! 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.233792Jy sigma=6.567631 ! Fit after self-cal, rms=0.233314Jy sigma=6.555673 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=4.466 mas, bmaj=17.16 mas, bpa=-5.115 degrees ! Estimated noise=0.842072 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00555571 Jy ! Component: 100 - total flux cleaned = -0.00560458 Jy ! Component: 150 - total flux cleaned = -0.00559309 Jy ! Component: 200 - total flux cleaned = -0.00589035 Jy ! Component: 250 - total flux cleaned = -0.00617362 Jy ! Component: 300 - total flux cleaned = -0.0067422 Jy ! Component: 350 - total flux cleaned = -0.00701274 Jy ! Component: 400 - total flux cleaned = -0.00702175 Jy ! Total flux subtracted in 400 components = -0.00702175 Jy ! Clean residual min=-0.030697 max=0.026582 Jy/beam ! Clean residual mean=-0.000084 rms=0.007665 Jy/beam ! Combined flux in latest and established models = 6.0034 Jy selfcal ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 357 components and 6.0034 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.230241Jy sigma=6.521222 ! Fit after self-cal, rms=0.229848Jy sigma=6.517743 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 = 0.00156949 Jy ! Component: 100 - total flux cleaned = 0.00213292 Jy ! Component: 150 - total flux cleaned = 0.00159891 Jy ! Component: 200 - total flux cleaned = 0.00135013 Jy ! Component: 250 - total flux cleaned = 0.00109584 Jy ! Component: 300 - total flux cleaned = 0.000593148 Jy ! Component: 350 - total flux cleaned = 0.000337306 Jy ! Component: 400 - total flux cleaned = 0.000336852 Jy ! Total flux subtracted in 400 components = 0.000336852 Jy ! Clean residual min=-0.030088 max=0.026320 Jy/beam ! Clean residual mean=-0.000087 rms=0.007550 Jy/beam ! Combined flux in latest and established models = 6.00374 Jy ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 365 components and 6.00374 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.227681Jy sigma=6.496496 ! Fit after self-cal, rms=0.227233Jy sigma=6.494064 ! 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, 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.227233Jy sigma=6.494064 ! Fit after self-cal, rms=0.226884Jy sigma=6.476874 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=9.596 mas, bmaj=30.37 mas, bpa=17.57 degrees ! Estimated noise=0.773768 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00401875 Jy ! Component: 100 - total flux cleaned = -0.003768 Jy ! Component: 150 - total flux cleaned = -0.00355036 Jy ! Component: 200 - total flux cleaned = -0.00370188 Jy ! Component: 250 - total flux cleaned = -0.00386733 Jy ! Component: 300 - total flux cleaned = -0.00387526 Jy ! Component: 350 - total flux cleaned = -0.00416479 Jy ! Component: 400 - total flux cleaned = -0.00444331 Jy ! Total flux subtracted in 400 components = -0.00444331 Jy ! Clean residual min=-0.023736 max=0.025785 Jy/beam ! Clean residual mean=-0.000217 rms=0.008022 Jy/beam ! Combined flux in latest and established models = 5.99929 Jy selfcal ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 387 components and 5.99929 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.226101Jy sigma=6.452675 ! Fit after self-cal, rms=0.225781Jy sigma=6.449604 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.4, 12). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0228592 Jy ! Component: 100 - total flux cleaned = 0.0307282 Jy ! Component: 150 - total flux cleaned = 0.0377085 Jy ! Component: 200 - total flux cleaned = 0.040149 Jy ! Component: 250 - total flux cleaned = 0.0430008 Jy ! Component: 300 - total flux cleaned = 0.0448921 Jy ! Component: 350 - total flux cleaned = 0.0466284 Jy ! Component: 400 - total flux cleaned = 0.0479604 Jy ! Total flux subtracted in 400 components = 0.0479604 Jy ! Clean residual min=-0.013839 max=0.012244 Jy/beam ! Clean residual mean=-0.000901 rms=0.004172 Jy/beam ! Combined flux in latest and established models = 6.04726 Jy ! Performing phase self-cal ! Adding 109 model components to the UV plane model. ! The established model now contains 496 components and 6.04726 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.221152Jy sigma=6.296877 ! Fit after self-cal, rms=0.220598Jy sigma=6.288834 ! 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=4.464 mas, bmaj=17.17 mas, bpa=-5.123 degrees ! Estimated noise=0.841124 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 4.464 x 17.17 at -5.123 degrees (North through East) ! Clean map min=-0.043341 max=4.0566 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.220598Jy sigma=6.288834 ! Fit after self-cal, rms=0.162884Jy sigma=3.305163 ! Inverting map and beam ! Estimated beam: bmin=4.462 mas, bmaj=17.17 mas, bpa=-5.138 degrees ! Estimated noise=0.840833 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00117196 Jy ! Component: 100 - total flux cleaned = 0.00347183 Jy ! Component: 150 - total flux cleaned = 0.00557245 Jy ! Component: 200 - total flux cleaned = 0.00754215 Jy ! Component: 250 - total flux cleaned = 0.00847794 Jy ! Component: 300 - total flux cleaned = 0.00936608 Jy ! Component: 350 - total flux cleaned = 0.00978927 Jy ! Component: 400 - total flux cleaned = 0.0110154 Jy ! Total flux subtracted in 400 components = 0.0110154 Jy ! Clean residual min=-0.022739 max=0.017970 Jy/beam ! Clean residual mean=-0.000702 rms=0.003913 Jy/beam ! Combined flux in latest and established models = 6.05827 Jy ! Performing phase self-cal ! Adding 53 model components to the UV plane model. ! The established model now contains 549 components and 6.05827 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.151347Jy sigma=3.154136 ! Fit after self-cal, rms=0.151428Jy sigma=3.148443 ! Inverting map ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00267775 Jy ! Component: 100 - total flux cleaned = 0.00308009 Jy ! Component: 150 - total flux cleaned = 0.00385272 Jy ! Component: 200 - total flux cleaned = 0.00421791 Jy ! Component: 250 - total flux cleaned = 0.00457166 Jy ! Component: 300 - total flux cleaned = 0.00456973 Jy ! Component: 350 - total flux cleaned = 0.00522546 Jy ! Component: 400 - total flux cleaned = 0.00522163 Jy ! Total flux subtracted in 400 components = 0.00522163 Jy ! Clean residual min=-0.021193 max=0.016828 Jy/beam ! Clean residual mean=-0.000710 rms=0.003550 Jy/beam ! Combined flux in latest and established models = 6.06349 Jy ! Performing phase self-cal ! Adding 68 model components to the UV plane model. ! The established model now contains 587 components and 6.06349 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.145083Jy sigma=3.072337 ! Fit after self-cal, rms=0.144966Jy sigma=3.069022 ! 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.144966Jy sigma=3.069022 ! Fit after self-cal, rms=0.144851Jy sigma=3.059825 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.231 mas, bmaj=15.36 mas, bpa=-6.871 degrees ! Estimated noise=1.42531 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 3.26458 Jy ! Component: 100 - total flux cleaned = 4.47266 Jy ! Component: 150 - total flux cleaned = 4.99453 Jy ! Component: 200 - total flux cleaned = 5.32076 Jy ! Component: 250 - total flux cleaned = 5.5414 Jy ! Component: 300 - total flux cleaned = 5.69159 Jy ! Component: 350 - total flux cleaned = 5.79521 Jy ! Component: 400 - total flux cleaned = 5.86727 Jy ! Total flux subtracted in 400 components = 5.86727 Jy ! Clean residual min=-0.064689 max=0.082445 Jy/beam ! Clean residual mean=0.000131 rms=0.014946 Jy/beam ! Combined flux in latest and established models = 5.86727 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 126 model components to the UV plane model. ! The established model now contains 126 components and 5.86727 Jy ! Inverting map and beam ! Estimated beam: bmin=4.462 mas, bmaj=17.17 mas, bpa=-5.125 degrees ! Estimated noise=0.840959 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 (-4.8, -4.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0779728 Jy ! Component: 100 - total flux cleaned = 0.114152 Jy ! Component: 150 - total flux cleaned = 0.122653 Jy ! Component: 200 - total flux cleaned = 0.127653 Jy ! Component: 250 - total flux cleaned = 0.132172 Jy ! Component: 300 - total flux cleaned = 0.137367 Jy ! Component: 350 - total flux cleaned = 0.142095 Jy ! Component: 400 - total flux cleaned = 0.146466 Jy ! Total flux subtracted in 400 components = 0.146466 Jy ! Clean residual min=-0.024431 max=0.019735 Jy/beam ! Clean residual mean=-0.000100 rms=0.005976 Jy/beam ! Combined flux in latest and established models = 6.01374 Jy ! Adding 95 model components to the UV plane model. ! The established model now contains 220 components and 6.01374 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=9.597 mas, bmaj=30.4 mas, bpa=17.58 degrees ! Estimated noise=0.773564 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0137619 Jy ! Component: 100 - total flux cleaned = 0.0195181 Jy ! Component: 150 - total flux cleaned = 0.0230594 Jy ! Component: 200 - total flux cleaned = 0.0262077 Jy ! Component: 250 - total flux cleaned = 0.0286092 Jy ! Component: 300 - total flux cleaned = 0.0304093 Jy ! Component: 350 - total flux cleaned = 0.0320735 Jy ! Component: 400 - total flux cleaned = 0.0330005 Jy ! Total flux subtracted in 400 components = 0.0330005 Jy ! Clean residual min=-0.012030 max=0.010470 Jy/beam ! Clean residual mean=-0.000709 rms=0.003637 Jy/beam ! Combined flux in latest and established models = 6.04674 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.000880938 Jy ! Component: 100 - total flux cleaned = 0.00170872 Jy ! Component: 150 - total flux cleaned = 0.0019768 Jy ! Component: 200 - total flux cleaned = 0.00299394 Jy ! Component: 250 - total flux cleaned = 0.00348177 Jy ! Component: 300 - total flux cleaned = 0.00395653 Jy ! Component: 350 - total flux cleaned = 0.00487286 Jy ! Component: 400 - total flux cleaned = 0.00553212 Jy ! Total flux subtracted in 400 components = 0.00553212 Jy ! Clean residual min=-0.010660 max=0.009819 Jy/beam ! Clean residual mean=-0.000752 rms=0.003135 Jy/beam ! Combined flux in latest and established models = 6.05227 Jy ! Adding 136 model components to the UV plane model. ! The established model now contains 351 components and 6.05227 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=4.462 mas, bmaj=17.17 mas, bpa=-5.125 degrees ! Estimated noise=0.840959 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00391542 Jy ! Component: 100 - total flux cleaned = -0.00576316 Jy ! Component: 150 - total flux cleaned = -0.00576871 Jy ! Component: 200 - total flux cleaned = -0.00524699 Jy ! Component: 250 - total flux cleaned = -0.00475717 Jy ! Component: 300 - total flux cleaned = -0.00428313 Jy ! Component: 350 - total flux cleaned = -0.00337871 Jy ! Component: 400 - total flux cleaned = -0.00209976 Jy ! Total flux subtracted in 400 components = -0.00209976 Jy ! Clean residual min=-0.011683 max=0.014496 Jy/beam ! Clean residual mean=-0.000461 rms=0.003238 Jy/beam ! Combined flux in latest and established models = 6.05017 Jy selfcal ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 392 components and 6.05017 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.135054Jy sigma=2.989357 ! Fit after self-cal, rms=0.130669Jy sigma=2.951775 wmodel /scr/pima/rdv62_uvs/J1924-2914_S_map.mod ! Writing 392 model components to file: /scr/pima/rdv62_uvs/J1924-2914_S_map.mod wobs /scr/pima/rdv62_uvs/J1924-2914_S_uvs.fits ! Applying 975 buffered edits. ! Writing UV FITS file: /scr/pima/rdv62_uvs/J1924-2914_S_uvs.fits wwins /scr/pima/rdv62_uvs/J1924-2914_S_map.win ! wwins: Wrote 4 windows to /scr/pima/rdv62_uvs/J1924-2914_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.0022659 Jy ! Component: 100 - total flux cleaned = 0.00168452 Jy ! Component: 150 - total flux cleaned = 0.00150259 Jy ! Component: 200 - total flux cleaned = 0.000629984 Jy ! Component: 250 - total flux cleaned = 0.000129096 Jy ! Component: 300 - total flux cleaned = 0.000285074 Jy ! Component: 350 - total flux cleaned = 0.000743024 Jy ! Component: 400 - total flux cleaned = 0.000743452 Jy ! Component: 450 - total flux cleaned = 0.000741989 Jy ! Component: 500 - total flux cleaned = 0.000316093 Jy ! Component: 550 - total flux cleaned = 0.000177212 Jy ! Component: 600 - total flux cleaned = -9.32806e-05 Jy ! Component: 650 - total flux cleaned = -9.28158e-05 Jy ! Component: 700 - total flux cleaned = -9.27603e-05 Jy ! Component: 750 - total flux cleaned = -9.34752e-05 Jy ! Component: 800 - total flux cleaned = -0.000217421 Jy ! Component: 850 - total flux cleaned = -0.000217817 Jy ! Component: 900 - total flux cleaned = -0.000581282 Jy ! Component: 950 - total flux cleaned = -0.000816975 Jy ! Component: 1000 - total flux cleaned = -0.00116412 Jy ! Component: 1050 - total flux cleaned = -0.00127871 Jy ! Component: 1100 - total flux cleaned = -0.00195914 Jy ! Component: 1150 - total flux cleaned = -0.00218226 Jy ! Component: 1200 - total flux cleaned = -0.00262431 Jy ! Component: 1250 - total flux cleaned = -0.00284186 Jy ! Component: 1300 - total flux cleaned = -0.00316478 Jy ! Component: 1350 - total flux cleaned = -0.00348532 Jy ! Component: 1400 - total flux cleaned = -0.00380145 Jy ! Component: 1450 - total flux cleaned = -0.00390569 Jy ! Component: 1500 - total flux cleaned = -0.00421732 Jy ! Component: 1550 - total flux cleaned = -0.004422 Jy ! Component: 1600 - total flux cleaned = -0.00482684 Jy ! Component: 1650 - total flux cleaned = -0.00502714 Jy ! Component: 1700 - total flux cleaned = -0.00522572 Jy ! Component: 1750 - total flux cleaned = -0.00542272 Jy ! Component: 1800 - total flux cleaned = -0.00590855 Jy ! Component: 1850 - total flux cleaned = -0.00590791 Jy ! Component: 1900 - total flux cleaned = -0.00600254 Jy ! Component: 1950 - total flux cleaned = -0.00609701 Jy ! Component: 2000 - total flux cleaned = -0.00619189 Jy ! Component: 2050 - total flux cleaned = -0.00675168 Jy ! Component: 2100 - total flux cleaned = -0.0065665 Jy ! Component: 2150 - total flux cleaned = -0.00702651 Jy ! Component: 2200 - total flux cleaned = -0.007299 Jy ! Component: 2250 - total flux cleaned = -0.00738861 Jy ! Component: 2300 - total flux cleaned = -0.0072992 Jy ! Component: 2350 - total flux cleaned = -0.00774667 Jy ! Component: 2400 - total flux cleaned = -0.00827744 Jy ! Component: 2450 - total flux cleaned = -0.00836468 Jy ! Component: 2500 - total flux cleaned = -0.00880106 Jy ! Component: 2550 - total flux cleaned = -0.00888697 Jy ! Component: 2600 - total flux cleaned = -0.00905975 Jy ! Component: 2650 - total flux cleaned = -0.00940134 Jy ! Component: 2700 - total flux cleaned = -0.00957126 Jy ! Component: 2750 - total flux cleaned = -0.0097404 Jy ! Component: 2800 - total flux cleaned = -0.0100763 Jy ! Component: 2850 - total flux cleaned = -0.00999245 Jy ! Component: 2900 - total flux cleaned = -0.0102419 Jy ! Component: 2950 - total flux cleaned = -0.0104067 Jy ! Component: 3000 - total flux cleaned = -0.0104878 Jy ! Component: 3050 - total flux cleaned = -0.0106503 Jy ! Component: 3100 - total flux cleaned = -0.0108927 Jy ! Component: 3150 - total flux cleaned = -0.0108919 Jy ! Component: 3200 - total flux cleaned = -0.0107329 Jy ! Component: 3250 - total flux cleaned = -0.0109717 Jy ! Component: 3300 - total flux cleaned = -0.0109713 Jy ! Component: 3350 - total flux cleaned = -0.010893 Jy ! Component: 3400 - total flux cleaned = -0.0109708 Jy ! Component: 3450 - total flux cleaned = -0.0105826 Jy ! Component: 3500 - total flux cleaned = -0.01066 Jy ! Component: 3550 - total flux cleaned = -0.0104298 Jy ! Component: 3600 - total flux cleaned = -0.010429 Jy ! Component: 3650 - total flux cleaned = -0.0102013 Jy ! Component: 3700 - total flux cleaned = -0.0102017 Jy ! Component: 3750 - total flux cleaned = -0.00997624 Jy ! Component: 3800 - total flux cleaned = -0.00982719 Jy ! Component: 3850 - total flux cleaned = -0.00975259 Jy ! Component: 3900 - total flux cleaned = -0.00967943 Jy ! Component: 3950 - total flux cleaned = -0.00967962 Jy ! Component: 4000 - total flux cleaned = -0.00953316 Jy ! Component: 4050 - total flux cleaned = -0.00938838 Jy ! Total flux subtracted in 4096 components = -0.00931619 Jy ! Clean residual min=-0.004449 max=0.003658 Jy/beam ! Clean residual mean=-0.000075 rms=0.001491 Jy/beam ! Combined flux in latest and established models = 6.04085 Jy keep ! Adding 228 model components to the UV plane model. ! The established model now contains 614 components and 6.04085 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv62_uvs/J1924-2914_S\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 4.462 x 17.17 at -5.125 degrees (North through East) ! Clean map min=-0.033103 max=3.8795 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv62_uvs/J1924-2914_S_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= 3.8794 rms= 0.00147 Jy quit ! Quitting program ! Log file /scr/pima/rdv62_uvs/J1924-2914_S_uvs.log closed on Fri Nov 22 22:40:26 2013