! Started logfile: /scr/pima/rdv62_uvs/J2258-2758_X_uvs.log on Fri Nov 22 22:43:27 2013 obs /scr/pima/rdv62_uvs/J2258-2758_X_uva.fits ! Reading UV FITS file: /scr/pima/rdv62_uvs/J2258-2758_X_uva.fits ! AN table 1: 7 integrations on 28 of 28 possible baselines. ! AN table 2: 12 integrations on 3 of 3 possible baselines. ! AN table 3: 3 integrations on 1 of 1 possible baselines. ! AN table 4: 40 integrations on 3 of 3 possible baselines. ! AN table 5: 3 integrations on 1 of 1 possible baselines. ! Apparent sampling: 0.662011 visibilities/baseline/integration-bin. ! Found source: J2258-2758 ! ! 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 8e+06 1 8e+06 ! 02 2 8.47599e+09 8e+06 1 8e+06 ! 03 3 8.79099e+09 8e+06 1 8e+06 ! 04 4 8.89599e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 948 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/rdv62_uvs/J2258-2758_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:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:KP 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:LA self-cal status=correctable weight=1 ! 2:MK self-cal status=correctable weight=1 ! 2:OV self-cal status=correctable weight=1 ! 3:KK self-cal status=correctable weight=1 ! 3:MK self-cal status=correctable weight=1 ! 4:HH self-cal status=correctable weight=1 ! 4:WZ self-cal status=correctable weight=1 ! 4:ZC self-cal status=correctable weight=1 ! 5:HH self-cal status=correctable weight=1 ! 5:ZC 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 12 telescope corrections were flagged in sub-array 1. ! A total of 30 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 51 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 30 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 51 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 30 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 51 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 30 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 120 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=1.888226Jy sigma=28.218005 ! Fit after self-cal, rms=0.681792Jy sigma=10.404932 ! 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 3:KK ! Flagging baselines of antenna 3:KK in the currently selected channels. flag 4:HH ! Flagging baselines of antenna 4:HH in the currently selected channels. flag 4:WZ ! Flagging baselines of antenna 4:WZ in the currently selected channels. flag 4:ZC ! Flagging baselines of antenna 4:ZC in the currently selected channels. flag 5:HH ! Flagging baselines of antenna 5:HH in the currently selected channels. flag 5:ZC ! Flagging baselines of antenna 5:ZC 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.911 mas, bmaj=4.322 mas, bpa=-6.617 degrees ! Estimated noise=2.701 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 = 1.23646 Jy ! Component: 100 - total flux cleaned = 1.5088 Jy ! Total flux subtracted in 100 components = 1.5088 Jy ! Clean residual min=-0.052204 max=0.086903 Jy/beam ! Clean residual mean=-0.000993 rms=0.016455 Jy/beam ! Combined flux in latest and established models = 1.5088 Jy ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 8 components and 1.5088 Jy ! ! Correcting IF 1. ! A total of 69 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 2. ! A total of 69 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 3. ! A total of 69 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.224630Jy sigma=4.159823 ! Fit after self-cal, rms=0.224682Jy sigma=4.159396 ! Inverting map and beam ! Estimated beam: bmin=0.911 mas, bmaj=4.322 mas, bpa=-6.617 degrees ! Estimated noise=2.701 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=1.301 mas, bmaj=4.862 mas, bpa=-4.334 degrees ! Estimated noise=1.77652 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.131156 Jy ! Component: 100 - total flux cleaned = 0.180166 Jy ! Component: 150 - total flux cleaned = 0.203069 Jy ! Component: 200 - total flux cleaned = 0.212406 Jy ! Total flux subtracted in 200 components = 0.212406 Jy ! Clean residual min=-0.061509 max=0.063002 Jy/beam ! Clean residual mean=-0.000178 rms=0.014441 Jy/beam ! Combined flux in latest and established models = 1.72121 Jy ! Performing phase self-cal ! Adding 44 model components to the UV plane model. ! The established model now contains 49 components and 1.72121 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.158062Jy sigma=2.721777 ! Fit after self-cal, rms=0.157869Jy sigma=2.719620 ! 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. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.96 FD 0.98 KP 0.98 LA 0.96 ! MK 0.98 NL 1.10 OV 1.00 PT 0.97 ! ! A total of 6 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! LA 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! KK 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! HH 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! HH 0.00* ZC 0.00* ! ! ! Correcting IF 2. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.98 FD 0.96 KP 0.93 LA 0.95 ! MK 0.98 NL 1.08 OV 0.97 PT 1.03 ! ! A total of 6 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! LA 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! KK 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! HH 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! HH 0.00* ZC 0.00* ! ! ! Correcting IF 3. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.02 FD 1.02 KP 0.87 LA 1.02 ! MK 0.98 NL 1.11 OV 1.04 PT 1.06 ! ! A total of 6 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! LA 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! KK 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! HH 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! HH 0.00* ZC 0.00* ! ! ! Correcting IF 4. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.03 FD 1.02 KP 0.92 LA 1.00 ! MK 1.01 NL 1.14 OV 1.04 PT 1.05 ! ! A total of 6 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! LA 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! KK 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! HH 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! HH 0.00* ZC 0.00* ! ! ! Fit before self-cal, rms=0.157869Jy sigma=2.719620 ! Fit after self-cal, rms=0.103069Jy sigma=1.693372 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.311 mas, bmaj=4.768 mas, bpa=-4.686 degrees ! Estimated noise=1.83384 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 (-3.6, -3.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0291456 Jy ! Component: 100 - total flux cleaned = 0.0295959 Jy ! Component: 150 - total flux cleaned = 0.0330282 Jy ! Component: 200 - total flux cleaned = 0.0330225 Jy ! Total flux subtracted in 200 components = 0.0330225 Jy ! Clean residual min=-0.014738 max=0.017936 Jy/beam ! Clean residual mean=-0.000026 rms=0.003916 Jy/beam ! Combined flux in latest and established models = 1.75423 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 66 components and 1.75423 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.090133Jy sigma=1.500066 ! Fit after self-cal, rms=0.089541Jy sigma=1.490376 ! 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.089541Jy sigma=1.490376 ! Fit after self-cal, rms=0.087741Jy sigma=1.475145 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.31 mas, bmaj=4.772 mas, bpa=-4.685 degrees ! Estimated noise=1.8297 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00722441 Jy ! Component: 100 - total flux cleaned = 0.00805648 Jy ! Component: 150 - total flux cleaned = 0.0076806 Jy ! Component: 200 - total flux cleaned = 0.00738626 Jy ! Total flux subtracted in 200 components = 0.00738626 Jy ! Clean residual min=-0.013840 max=0.013562 Jy/beam ! Clean residual mean=-0.000037 rms=0.003090 Jy/beam ! Combined flux in latest and established models = 1.76162 Jy selfcal ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 75 components and 1.76162 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.085932Jy sigma=1.450634 ! Fit after self-cal, rms=0.085689Jy sigma=1.447608 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.0022362 Jy ! Component: 100 - total flux cleaned = 0.00211935 Jy ! Component: 150 - total flux cleaned = 0.00211611 Jy ! Component: 200 - total flux cleaned = 0.00211706 Jy ! Total flux subtracted in 200 components = 0.00211706 Jy ! Clean residual min=-0.013322 max=0.012802 Jy/beam ! Clean residual mean=-0.000028 rms=0.002959 Jy/beam ! Combined flux in latest and established models = 1.76374 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 87 components and 1.76374 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.085327Jy sigma=1.442651 ! Fit after self-cal, rms=0.085284Jy sigma=1.442316 ! 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.085284Jy sigma=1.442316 ! Fit after self-cal, rms=0.084666Jy sigma=1.437050 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.31 mas, bmaj=4.777 mas, bpa=-4.694 degrees ! Estimated noise=1.82886 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00397165 Jy ! Component: 100 - total flux cleaned = 0.00470329 Jy ! Component: 150 - total flux cleaned = 0.00481173 Jy ! Component: 200 - total flux cleaned = 0.00480943 Jy ! Total flux subtracted in 200 components = 0.00480943 Jy ! Clean residual min=-0.012353 max=0.011689 Jy/beam ! Clean residual mean=-0.000020 rms=0.002692 Jy/beam ! Combined flux in latest and established models = 1.76854 Jy selfcal ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 99 components and 1.76854 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.084263Jy sigma=1.430785 ! Fit after self-cal, rms=0.084213Jy sigma=1.430124 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.00132905 Jy ! Component: 100 - total flux cleaned = 0.00140717 Jy ! Component: 150 - total flux cleaned = 0.00134094 Jy ! Component: 200 - total flux cleaned = 0.00114396 Jy ! Total flux subtracted in 200 components = 0.00114396 Jy ! Clean residual min=-0.011916 max=0.011172 Jy/beam ! Clean residual mean=-0.000015 rms=0.002621 Jy/beam ! Combined flux in latest and established models = 1.76969 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 109 components and 1.76969 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.084077Jy sigma=1.427894 ! Fit after self-cal, rms=0.084069Jy sigma=1.427806 ! 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.084069Jy sigma=1.427806 ! Fit after self-cal, rms=0.083871Jy sigma=1.426024 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.683 mas, bmaj=8.807 mas, bpa=17.81 degrees ! Estimated noise=1.95567 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00153195 Jy ! Component: 100 - total flux cleaned = 0.00203897 Jy ! Component: 150 - total flux cleaned = 0.00188607 Jy ! Component: 200 - total flux cleaned = 0.00175294 Jy ! Total flux subtracted in 200 components = 0.00175294 Jy ! Clean residual min=-0.010609 max=0.010592 Jy/beam ! Clean residual mean=-0.000036 rms=0.002790 Jy/beam ! Combined flux in latest and established models = 1.77144 Jy selfcal ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 118 components and 1.77144 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.083750Jy sigma=1.423723 ! Fit after self-cal, rms=0.083716Jy sigma=1.423324 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.000782929 Jy ! Component: 100 - total flux cleaned = 0.000843671 Jy ! Component: 150 - total flux cleaned = 0.000951853 Jy ! Component: 200 - total flux cleaned = 0.00100003 Jy ! Total flux subtracted in 200 components = 0.00100003 Jy ! Clean residual min=-0.010306 max=0.010330 Jy/beam ! Clean residual mean=-0.000030 rms=0.002732 Jy/beam ! Combined flux in latest and established models = 1.77244 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 125 components and 1.77244 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.083712Jy sigma=1.422989 ! Fit after self-cal, rms=0.083703Jy sigma=1.422896 ! 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:FD,true ! 1:FD self-cal status=fixed weight=1 selfant 1:KP,true ! 1:KP 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:LA,true ! 2:LA self-cal status=fixed weight=1 selfant 2:MK,true ! 2:MK self-cal status=fixed weight=1 selfant 2:OV,true ! 2:OV self-cal status=fixed weight=1 selfant 3:MK,true ! 3:MK 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 3 telescope corrections were flagged in sub-array 3. ! A total of 120 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 120 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 120 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 120 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.916739Jy sigma=8.602071 ! Fit after self-cal, rms=0.580540Jy sigma=8.447277 selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:KP 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:LA self-cal status=correctable weight=1 ! 2:MK self-cal status=correctable weight=1 ! 2:OV self-cal status=correctable weight=1 ! 3:KK self-cal status=correctable weight=1 ! 3:MK self-cal status=correctable weight=1 ! 4:HH self-cal status=correctable weight=1 ! 4:WZ self-cal status=correctable weight=1 ! 4:ZC self-cal status=correctable weight=1 ! 5:HH self-cal status=correctable weight=1 ! 5:ZC 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 15 telescope corrections were flagged in sub-array 1. ! A total of 36 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 2. ! A total of 15 telescope corrections were flagged in sub-array 1. ! A total of 36 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 3. ! A total of 15 telescope corrections were flagged in sub-array 1. ! A total of 36 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 4. ! A total of 15 telescope corrections were flagged in sub-array 1. ! A total of 36 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! ! Fit before self-cal, rms=0.580540Jy sigma=8.447277 ! Fit after self-cal, rms=0.083577Jy sigma=1.421820 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.9177 mas, bmaj=4.207 mas, bpa=-7.479 degrees ! Estimated noise=2.5991 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.0015115 Jy ! Component: 100 - total flux cleaned = 0.000360513 Jy ! Component: 150 - total flux cleaned = -0.00017782 Jy ! Component: 200 - total flux cleaned = -0.000304838 Jy ! Total flux subtracted in 200 components = -0.000304838 Jy ! Clean residual min=-0.010550 max=0.010926 Jy/beam ! Clean residual mean=-0.000006 rms=0.002753 Jy/beam ! Combined flux in latest and established models = 1.77214 Jy ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 151 components and 1.77214 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.083556Jy sigma=1.421266 ! Fit after self-cal, rms=0.083530Jy sigma=1.420962 ! 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=1.308 mas, bmaj=4.783 mas, bpa=-4.717 degrees ! Estimated noise=1.82854 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.000274205 Jy ! Component: 100 - total flux cleaned = -0.000432132 Jy ! Component: 150 - total flux cleaned = -0.00063091 Jy ! Component: 200 - total flux cleaned = -0.000691094 Jy ! Component: 250 - total flux cleaned = -0.000752237 Jy ! Component: 300 - total flux cleaned = -0.000910328 Jy ! Component: 350 - total flux cleaned = -0.000959619 Jy ! Component: 400 - total flux cleaned = -0.00100709 Jy ! Total flux subtracted in 400 components = -0.00100709 Jy ! Clean residual min=-0.010106 max=0.009636 Jy/beam ! Clean residual mean=-0.000029 rms=0.002348 Jy/beam ! Combined flux in latest and established models = 1.77113 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 178 components and 1.77113 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.083293Jy sigma=1.417488 ! Fit after self-cal, rms=0.083299Jy sigma=1.417444 ! 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.000234762 Jy ! Component: 100 - total flux cleaned = 0.000275492 Jy ! Component: 150 - total flux cleaned = 0.000236398 Jy ! Component: 200 - total flux cleaned = 0.000274897 Jy ! Component: 250 - total flux cleaned = 0.000236897 Jy ! Component: 300 - total flux cleaned = 0.000236542 Jy ! Component: 350 - total flux cleaned = 0.00027125 Jy ! Component: 400 - total flux cleaned = 0.000271348 Jy ! Total flux subtracted in 400 components = 0.000271348 Jy ! Clean residual min=-0.009892 max=0.009564 Jy/beam ! Clean residual mean=-0.000028 rms=0.002320 Jy/beam ! Combined flux in latest and established models = 1.7714 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 194 components and 1.7714 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.083241Jy sigma=1.416417 ! Fit after self-cal, rms=0.083243Jy sigma=1.416409 ! 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 15 telescope corrections were flagged in sub-array 1. ! A total of 36 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 120 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 15 telescope corrections were flagged in sub-array 1. ! A total of 36 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 120 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 15 telescope corrections were flagged in sub-array 1. ! A total of 36 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 120 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 15 telescope corrections were flagged in sub-array 1. ! A total of 36 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 120 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.916259Jy sigma=8.608265 ! Fit after self-cal, rms=0.083137Jy sigma=1.415473 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.307 mas, bmaj=4.786 mas, bpa=-4.722 degrees ! Estimated noise=1.82839 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000852046 Jy ! Component: 100 - total flux cleaned = 0.000901235 Jy ! Component: 150 - total flux cleaned = 0.000775816 Jy ! Component: 200 - total flux cleaned = 0.00081467 Jy ! Component: 250 - total flux cleaned = 0.000776492 Jy ! Component: 300 - total flux cleaned = 0.000776081 Jy ! Component: 350 - total flux cleaned = 0.000741996 Jy ! Component: 400 - total flux cleaned = 0.000674391 Jy ! Total flux subtracted in 400 components = 0.000674391 Jy ! Clean residual min=-0.009778 max=0.009185 Jy/beam ! Clean residual mean=-0.000034 rms=0.002249 Jy/beam ! Combined flux in latest and established models = 1.77208 Jy selfcal ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 214 components and 1.77208 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.083044Jy sigma=1.413876 ! Fit after self-cal, rms=0.083041Jy sigma=1.413789 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.000420216 Jy ! Component: 100 - total flux cleaned = 0.000387992 Jy ! Component: 150 - total flux cleaned = 0.000388713 Jy ! Component: 200 - total flux cleaned = 0.000326423 Jy ! Component: 250 - total flux cleaned = 0.000296348 Jy ! Component: 300 - total flux cleaned = 0.000266508 Jy ! Component: 350 - total flux cleaned = 0.000236112 Jy ! Component: 400 - total flux cleaned = 0.000207417 Jy ! Total flux subtracted in 400 components = 0.000207417 Jy ! Clean residual min=-0.009641 max=0.009079 Jy/beam ! Clean residual mean=-0.000034 rms=0.002223 Jy/beam ! Combined flux in latest and established models = 1.77228 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 222 components and 1.77228 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.082999Jy sigma=1.413022 ! Fit after self-cal, rms=0.082999Jy sigma=1.413002 ! 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.082999Jy sigma=1.413002 ! Fit after self-cal, rms=0.082934Jy sigma=1.412486 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.307 mas, bmaj=4.788 mas, bpa=-4.725 degrees ! Estimated noise=1.82837 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000428301 Jy ! Component: 100 - total flux cleaned = 0.000479474 Jy ! Component: 150 - total flux cleaned = 0.000479978 Jy ! Component: 200 - total flux cleaned = 0.000480967 Jy ! Component: 250 - total flux cleaned = 0.000481613 Jy ! Component: 300 - total flux cleaned = 0.000450771 Jy ! Component: 350 - total flux cleaned = 0.000419526 Jy ! Component: 400 - total flux cleaned = 0.000389936 Jy ! Total flux subtracted in 400 components = 0.000389936 Jy ! Clean residual min=-0.009493 max=0.008800 Jy/beam ! Clean residual mean=-0.000036 rms=0.002173 Jy/beam ! Combined flux in latest and established models = 1.77267 Jy selfcal ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 229 components and 1.77267 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.082866Jy sigma=1.411288 ! Fit after self-cal, rms=0.082862Jy sigma=1.411212 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.000342518 Jy ! Component: 100 - total flux cleaned = 0.000342884 Jy ! Component: 150 - total flux cleaned = 0.000314174 Jy ! Component: 200 - total flux cleaned = 0.000314481 Jy ! Component: 250 - total flux cleaned = 0.000313602 Jy ! Component: 300 - total flux cleaned = 0.000314492 Jy ! Component: 350 - total flux cleaned = 0.000340778 Jy ! Component: 400 - total flux cleaned = 0.000341021 Jy ! Total flux subtracted in 400 components = 0.000341021 Jy ! Clean residual min=-0.009386 max=0.008726 Jy/beam ! Clean residual mean=-0.000035 rms=0.002152 Jy/beam ! Combined flux in latest and established models = 1.77301 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 232 components and 1.77301 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.082829Jy sigma=1.410601 ! Fit after self-cal, rms=0.082828Jy sigma=1.410578 ! 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.082828Jy sigma=1.410578 ! Fit after self-cal, rms=0.081836Jy sigma=1.389203 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.684 mas, bmaj=8.808 mas, bpa=17.72 degrees ! Estimated noise=1.95615 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000672712 Jy ! Component: 100 - total flux cleaned = 0.000972524 Jy ! Component: 150 - total flux cleaned = 0.00113277 Jy ! Component: 200 - total flux cleaned = 0.00127121 Jy ! Component: 250 - total flux cleaned = 0.00120956 Jy ! Component: 300 - total flux cleaned = 0.00115129 Jy ! Component: 350 - total flux cleaned = 0.00107023 Jy ! Component: 400 - total flux cleaned = 0.00104527 Jy ! Total flux subtracted in 400 components = 0.00104527 Jy ! Clean residual min=-0.008694 max=0.008147 Jy/beam ! Clean residual mean=-0.000047 rms=0.002316 Jy/beam ! Combined flux in latest and established models = 1.77406 Jy selfcal ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 246 components and 1.77406 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.081825Jy sigma=1.388745 ! Fit after self-cal, rms=0.081821Jy sigma=1.388664 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.000287806 Jy ! Component: 100 - total flux cleaned = 0.000289465 Jy ! Component: 150 - total flux cleaned = 0.000242291 Jy ! Component: 200 - total flux cleaned = 0.000198595 Jy ! Component: 250 - total flux cleaned = 0.000134115 Jy ! Component: 300 - total flux cleaned = 0.000134328 Jy ! Component: 350 - total flux cleaned = 0.000134193 Jy ! Component: 400 - total flux cleaned = 0.00015388 Jy ! Total flux subtracted in 400 components = 0.00015388 Jy ! Clean residual min=-0.008541 max=0.007960 Jy/beam ! Clean residual mean=-0.000045 rms=0.002293 Jy/beam ! Combined flux in latest and established models = 1.77421 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 254 components and 1.77421 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.081818Jy sigma=1.388551 ! Fit after self-cal, rms=0.081818Jy sigma=1.388531 ! 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=1.307 mas, bmaj=4.791 mas, bpa=-4.725 degrees ! Estimated noise=1.82881 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.307 x 4.791 at -4.725 degrees (North through East) ! Clean map min=-0.0088946 max=1.6273 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.081818Jy sigma=1.388531 ! Fit after self-cal, rms=0.074678Jy sigma=1.261770 ! Inverting map and beam ! Estimated beam: bmin=1.307 mas, bmaj=4.792 mas, bpa=-4.727 degrees ! Estimated noise=1.82848 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -7.44365e-05 Jy ! Component: 100 - total flux cleaned = -1.36427e-05 Jy ! Component: 150 - total flux cleaned = 4.3515e-05 Jy ! Component: 200 - total flux cleaned = -4.69521e-06 Jy ! Component: 250 - total flux cleaned = -4.48004e-06 Jy ! Component: 300 - total flux cleaned = -4.60603e-05 Jy ! Component: 350 - total flux cleaned = -0.000123985 Jy ! Component: 400 - total flux cleaned = -0.000160415 Jy ! Total flux subtracted in 400 components = -0.000160415 Jy ! Clean residual min=-0.009234 max=0.008020 Jy/beam ! Clean residual mean=-0.000027 rms=0.002057 Jy/beam ! Combined flux in latest and established models = 1.77405 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 278 components and 1.77405 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.074531Jy sigma=1.259613 ! Fit after self-cal, rms=0.074530Jy sigma=1.259563 ! Inverting map ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000112433 Jy ! Component: 100 - total flux cleaned = 0.000112739 Jy ! Component: 150 - total flux cleaned = 0.000143205 Jy ! Component: 200 - total flux cleaned = 0.000143662 Jy ! Component: 250 - total flux cleaned = 0.000144377 Jy ! Component: 300 - total flux cleaned = 0.000200609 Jy ! Component: 350 - total flux cleaned = 0.000200951 Jy ! Component: 400 - total flux cleaned = 0.000202301 Jy ! Total flux subtracted in 400 components = 0.000202301 Jy ! Clean residual min=-0.009209 max=0.007942 Jy/beam ! Clean residual mean=-0.000027 rms=0.002045 Jy/beam ! Combined flux in latest and established models = 1.77426 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 289 components and 1.77426 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.074485Jy sigma=1.258860 ! Fit after self-cal, rms=0.074484Jy sigma=1.258850 ! 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.074484Jy sigma=1.258850 ! Fit after self-cal, rms=0.074442Jy sigma=1.258544 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.9175 mas, bmaj=4.218 mas, bpa=-7.471 degrees ! Estimated noise=2.59944 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 1.24535 Jy ! Component: 100 - total flux cleaned = 1.52306 Jy ! Component: 150 - total flux cleaned = 1.61736 Jy ! Component: 200 - total flux cleaned = 1.66715 Jy ! Component: 250 - total flux cleaned = 1.69826 Jy ! Component: 300 - total flux cleaned = 1.71969 Jy ! Component: 350 - total flux cleaned = 1.73484 Jy ! Component: 400 - total flux cleaned = 1.74601 Jy ! Total flux subtracted in 400 components = 1.74601 Jy ! Clean residual min=-0.012972 max=0.011745 Jy/beam ! Clean residual mean=-0.000124 rms=0.003243 Jy/beam ! Combined flux in latest and established models = 1.74601 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 80 model components to the UV plane model. ! The established model now contains 80 components and 1.74601 Jy ! Inverting map and beam ! Estimated beam: bmin=1.307 mas, bmaj=4.794 mas, bpa=-4.727 degrees ! Estimated noise=1.82853 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam keep new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0121184 Jy ! Component: 100 - total flux cleaned = 0.0191507 Jy ! Component: 150 - total flux cleaned = 0.0221537 Jy ! Component: 200 - total flux cleaned = 0.0242715 Jy ! Component: 250 - total flux cleaned = 0.0253188 Jy ! Component: 300 - total flux cleaned = 0.0262325 Jy ! Component: 350 - total flux cleaned = 0.0268142 Jy ! Component: 400 - total flux cleaned = 0.0272349 Jy ! Total flux subtracted in 400 components = 0.0272349 Jy ! Clean residual min=-0.010136 max=0.008775 Jy/beam ! Clean residual mean=-0.000033 rms=0.002241 Jy/beam ! Combined flux in latest and established models = 1.77324 Jy ! Adding 50 model components to the UV plane model. ! The established model now contains 125 components and 1.77324 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.684 mas, bmaj=8.806 mas, bpa=17.69 degrees ! Estimated noise=1.95585 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00220622 Jy ! Component: 100 - total flux cleaned = 0.00323008 Jy ! Component: 150 - total flux cleaned = 0.0033861 Jy ! Component: 200 - total flux cleaned = 0.00331534 Jy ! Component: 250 - total flux cleaned = 0.00337683 Jy ! Component: 300 - total flux cleaned = 0.00337523 Jy ! Component: 350 - total flux cleaned = 0.00332572 Jy ! Component: 400 - total flux cleaned = 0.00332545 Jy ! Total flux subtracted in 400 components = 0.00332545 Jy ! Clean residual min=-0.008921 max=0.007972 Jy/beam ! Clean residual mean=-0.000029 rms=0.002258 Jy/beam ! Combined flux in latest and established models = 1.77657 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 = 7.47605e-09 Jy ! Component: 100 - total flux cleaned = 3.41761e-07 Jy ! Component: 150 - total flux cleaned = 1.20853e-06 Jy ! Component: 200 - total flux cleaned = 1.04762e-06 Jy ! Component: 250 - total flux cleaned = 1.64937e-06 Jy ! Component: 300 - total flux cleaned = 3.71001e-05 Jy ! Component: 350 - total flux cleaned = 3.72109e-05 Jy ! Component: 400 - total flux cleaned = 3.75594e-05 Jy ! Total flux subtracted in 400 components = 3.75594e-05 Jy ! Clean residual min=-0.008628 max=0.007841 Jy/beam ! Clean residual mean=-0.000028 rms=0.002229 Jy/beam ! Combined flux in latest and established models = 1.77661 Jy ! Adding 35 model components to the UV plane model. ! The established model now contains 154 components and 1.77661 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=1.307 mas, bmaj=4.794 mas, bpa=-4.727 degrees ! Estimated noise=1.82853 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00108005 Jy ! Component: 100 - total flux cleaned = -0.00107242 Jy ! Component: 150 - total flux cleaned = -0.000980406 Jy ! Component: 200 - total flux cleaned = -0.000981522 Jy ! Component: 250 - total flux cleaned = -0.000981427 Jy ! Component: 300 - total flux cleaned = -0.00098221 Jy ! Component: 350 - total flux cleaned = -0.00105175 Jy ! Component: 400 - total flux cleaned = -0.00105091 Jy ! Total flux subtracted in 400 components = -0.00105091 Jy ! Clean residual min=-0.009221 max=0.008204 Jy/beam ! Clean residual mean=-0.000023 rms=0.002061 Jy/beam ! Combined flux in latest and established models = 1.77556 Jy selfcal ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 175 components and 1.77556 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.074599Jy sigma=1.260808 ! Fit after self-cal, rms=0.074546Jy sigma=1.260187 wmodel /scr/pima/rdv62_uvs/J2258-2758_X_map.mod ! Writing 175 model components to file: /scr/pima/rdv62_uvs/J2258-2758_X_map.mod wobs /scr/pima/rdv62_uvs/J2258-2758_X_uvs.fits ! Applying 845 buffered edits. ! Writing UV FITS file: /scr/pima/rdv62_uvs/J2258-2758_X_uvs.fits wwins /scr/pima/rdv62_uvs/J2258-2758_X_map.win ! wwins: Wrote 2 windows to /scr/pima/rdv62_uvs/J2258-2758_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 = -0.00120754 Jy ! Component: 100 - total flux cleaned = -0.00185902 Jy ! Component: 150 - total flux cleaned = -0.00187035 Jy ! Component: 200 - total flux cleaned = -0.00144484 Jy ! Component: 250 - total flux cleaned = -0.00134258 Jy ! Component: 300 - total flux cleaned = -0.00115069 Jy ! Component: 350 - total flux cleaned = -0.00106319 Jy ! Component: 400 - total flux cleaned = -0.000980577 Jy ! Component: 450 - total flux cleaned = -0.00106144 Jy ! Component: 500 - total flux cleaned = -0.00113894 Jy ! Component: 550 - total flux cleaned = -0.00121414 Jy ! Component: 600 - total flux cleaned = -0.00121504 Jy ! Component: 650 - total flux cleaned = -0.00121375 Jy ! Component: 700 - total flux cleaned = -0.00100804 Jy ! Component: 750 - total flux cleaned = -0.000678047 Jy ! Component: 800 - total flux cleaned = -0.000421623 Jy ! Component: 850 - total flux cleaned = -0.000236409 Jy ! Component: 900 - total flux cleaned = 5.51311e-06 Jy ! Component: 950 - total flux cleaned = 0.000300938 Jy ! Component: 1000 - total flux cleaned = 0.000702868 Jy ! Component: 1050 - total flux cleaned = 0.000760982 Jy ! Component: 1100 - total flux cleaned = 0.00109198 Jy ! Component: 1150 - total flux cleaned = 0.001361 Jy ! Component: 1200 - total flux cleaned = 0.00172922 Jy ! Component: 1250 - total flux cleaned = 0.00183219 Jy ! Component: 1300 - total flux cleaned = 0.00208787 Jy ! Component: 1350 - total flux cleaned = 0.00244098 Jy ! Component: 1400 - total flux cleaned = 0.00283889 Jy ! Component: 1450 - total flux cleaned = 0.00303548 Jy ! Component: 1500 - total flux cleaned = 0.00318162 Jy ! Component: 1550 - total flux cleaned = 0.00337318 Jy ! Component: 1600 - total flux cleaned = 0.00375273 Jy ! Component: 1650 - total flux cleaned = 0.00398739 Jy ! Component: 1700 - total flux cleaned = 0.0042195 Jy ! Component: 1750 - total flux cleaned = 0.00435795 Jy ! Component: 1800 - total flux cleaned = 0.00449514 Jy ! Component: 1850 - total flux cleaned = 0.00463085 Jy ! Component: 1900 - total flux cleaned = 0.00472085 Jy ! Component: 1950 - total flux cleaned = 0.00494422 Jy ! Component: 2000 - total flux cleaned = 0.00507768 Jy ! Component: 2050 - total flux cleaned = 0.00521065 Jy ! Component: 2100 - total flux cleaned = 0.00521081 Jy ! Component: 2150 - total flux cleaned = 0.0053423 Jy ! Component: 2200 - total flux cleaned = 0.00560438 Jy ! Component: 2250 - total flux cleaned = 0.00569162 Jy ! Component: 2300 - total flux cleaned = 0.00577841 Jy ! Component: 2350 - total flux cleaned = 0.00586515 Jy ! Component: 2400 - total flux cleaned = 0.00599486 Jy ! Component: 2450 - total flux cleaned = 0.00629626 Jy ! Component: 2500 - total flux cleaned = 0.00646815 Jy ! Component: 2550 - total flux cleaned = 0.0067257 Jy ! Component: 2600 - total flux cleaned = 0.00694006 Jy ! Component: 2650 - total flux cleaned = 0.00702587 Jy ! Component: 2700 - total flux cleaned = 0.00728282 Jy ! Component: 2750 - total flux cleaned = 0.0074546 Jy ! Component: 2800 - total flux cleaned = 0.00766922 Jy ! Component: 2850 - total flux cleaned = 0.00779802 Jy ! Component: 2900 - total flux cleaned = 0.00809903 Jy ! Component: 2950 - total flux cleaned = 0.00840052 Jy ! Component: 3000 - total flux cleaned = 0.00852997 Jy ! Component: 3050 - total flux cleaned = 0.00870283 Jy ! Component: 3100 - total flux cleaned = 0.00883277 Jy ! Component: 3150 - total flux cleaned = 0.00909329 Jy ! Component: 3200 - total flux cleaned = 0.00926741 Jy ! Component: 3250 - total flux cleaned = 0.00952959 Jy ! Component: 3300 - total flux cleaned = 0.00970505 Jy ! Component: 3350 - total flux cleaned = 0.009925 Jy ! Component: 3400 - total flux cleaned = 0.0100577 Jy ! Component: 3450 - total flux cleaned = 0.010191 Jy ! Component: 3500 - total flux cleaned = 0.0103253 Jy ! Component: 3550 - total flux cleaned = 0.0104604 Jy ! Component: 3600 - total flux cleaned = 0.0105054 Jy ! Component: 3650 - total flux cleaned = 0.0106417 Jy ! Component: 3700 - total flux cleaned = 0.0107331 Jy ! Component: 3750 - total flux cleaned = 0.0107791 Jy ! Component: 3800 - total flux cleaned = 0.0109176 Jy ! Component: 3850 - total flux cleaned = 0.0110106 Jy ! Component: 3900 - total flux cleaned = 0.0111978 Jy ! Component: 3950 - total flux cleaned = 0.0112447 Jy ! Component: 4000 - total flux cleaned = 0.0113397 Jy ! Component: 4050 - total flux cleaned = 0.0114355 Jy ! Total flux subtracted in 4096 components = 0.0115802 Jy ! Clean residual min=-0.003174 max=0.002855 Jy/beam ! Clean residual mean=0.000040 rms=0.000796 Jy/beam ! Combined flux in latest and established models = 1.78714 Jy keep ! Adding 227 model components to the UV plane model. ! The established model now contains 402 components and 1.78714 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv62_uvs/J2258-2758_X\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.307 x 4.794 at -4.727 degrees (North through East) ! Clean map min=-0.0078944 max=1.6159 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv62_uvs/J2258-2758_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= 1.6158 rms= 0.00077 Jy quit ! Quitting program ! Log file /scr/pima/rdv62_uvs/J2258-2758_X_uvs.log closed on Fri Nov 22 22:43:32 2013