! Started logfile: /scr/pima/rdv93_uvs/J1427-4206_X_uvs.log on Tue Oct 8 15:48:10 2013 obs /scr/pima/rdv93_uvs/J1427-4206_X_uva.fits ! Reading UV FITS file: /scr/pima/rdv93_uvs/J1427-4206_X_uva.fits ! AN table 1: 5 integrations on 1 of 1 possible baselines. ! AN table 2: 3 integrations on 3 of 3 possible baselines. ! AN table 3: 3 integrations on 6 of 6 possible baselines. ! AN table 4: 6 integrations on 21 of 21 possible baselines. ! AN table 5: 8 integrations on 28 of 28 possible baselines. ! Apparent sampling: 0.685864 visibilities/baseline/integration-bin. ! Found source: J1427-4206 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 8.40599e+09 1.6e+07 1 1.6e+07 ! 02 2 8.47599e+09 1.6e+07 1 1.6e+07 ! 03 3 8.79099e+09 1.6e+07 1 1.6e+07 ! 04 4 8.89599e+09 1.6e+07 1 1.6e+07 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 1048 visibilities. select RR ! Selecting polarization: RR, channels: 1..4 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 print "field_size =",field_size ! field_size = 1024 print "field_cell =",field_cell ! field_cell = 0.2 print "taper_size =",taper_size ! taper_size = 70 ![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv93_uvs/J1427-4206_X] integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 800000 float soltime2 soltime2 = 15 float soltime3 soltime3 = 1 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float old_peak float new_peak float flux_cutoff selfant "",false ! 1:FT self-cal status=correctable weight=1 ! 1:SC self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:FT self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 ! 3:FT self-cal status=correctable weight=1 ! 3:KP self-cal status=correctable weight=1 ! 3:LA self-cal status=correctable weight=1 ! 3:PT self-cal status=correctable weight=1 ! 4:FD self-cal status=correctable weight=1 ! 4:FT self-cal status=correctable weight=1 ! 4:LA self-cal status=correctable weight=1 ! 4:NL self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 4:PT self-cal status=correctable weight=1 ! 4:SC self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:FT self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:MK self-cal status=correctable weight=1 ! 5:NL self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 5:SC self-cal status=correctable weight=1 system "rm /tmp/ff.ff" mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.200x0.200 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 36 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 36 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 36 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 36 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=2.460492Jy sigma=95.173777 ! Fit after self-cal, rms=1.268829Jy sigma=52.903723 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. print "=== Flag out IVS stations" ! === Flag out IVS stations integer lun string cfg_file lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_flag_process.csh /tmp/ff.ff" ![@/tmp/ff.ff \] flag 1:FT ! Flagging baselines of antenna 1:FT in the currently selected channels. flag 2:FT ! Flagging baselines of antenna 2:FT in the currently selected channels. flag 3:FT ! Flagging baselines of antenna 3:FT in the currently selected channels. flag 4:FT ! Flagging baselines of antenna 4:FT in the currently selected channels. flag 5:FT ! Flagging baselines of antenna 5:FT in the currently selected channels. ![Exited script file: /tmp/ff.ff] uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.6575 mas, bmaj=2.229 mas, bpa=2.004 degrees ! Estimated noise=1.04502 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.75175 Jy ! Component: 100 - total flux cleaned = 2.13647 Jy ! Total flux subtracted in 100 components = 2.13647 Jy ! Clean residual min=-0.102226 max=0.212290 Jy/beam ! Clean residual mean=0.001612 rms=0.029937 Jy/beam ! Combined flux in latest and established models = 2.13647 Jy ! Performing phase self-cal ! Adding 3 model components to the UV plane model. ! The established model now contains 3 components and 2.13647 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.492099Jy sigma=16.702590 ! Fit after self-cal, rms=0.492183Jy sigma=16.696883 ! Inverting map ! Added new window around map position (-1.6, -0.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.185715 Jy ! Component: 100 - total flux cleaned = 0.276009 Jy ! Total flux subtracted in 100 components = 0.276009 Jy ! Clean residual min=-0.072717 max=0.165887 Jy/beam ! Clean residual mean=0.000921 rms=0.022081 Jy/beam ! Combined flux in latest and established models = 2.41248 Jy ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 10 components and 2.41247 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.413530Jy sigma=12.560130 ! Fit after self-cal, rms=0.409781Jy sigma=11.940547 ! Inverting map print "=== Finished uniform weighting clean ===" ! === Finished uniform weighting clean === uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.8108 mas, bmaj=2.595 mas, bpa=1.227 degrees ! Estimated noise=0.813828 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.127183 Jy ! Component: 100 - total flux cleaned = 0.197135 Jy ! Component: 150 - total flux cleaned = 0.237273 Jy ! Component: 200 - total flux cleaned = 0.24422 Jy ! Total flux subtracted in 200 components = 0.24422 Jy ! Clean residual min=-0.080703 max=0.131487 Jy/beam ! Clean residual mean=0.000452 rms=0.021213 Jy/beam ! Combined flux in latest and established models = 2.65669 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 19 components and 2.65669 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.389132Jy sigma=9.299801 ! Fit after self-cal, rms=0.391573Jy sigma=9.104247 ! Inverting map print "=== Finished natural weighting clean for VLBI sites ===" ! === Finished natural weighting clean for VLBI sites === gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! FT 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.00* FT 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FT 0.00* KP 0.00* LA 0.00* PT 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! FD 0.95 FT 1.00* LA 0.99 NL 1.25 ! OV 0.96 PT 0.99 SC 1.04 ! ! Telescope amplitude corrections in sub-array 5: ! FD 0.95 FT 1.00* LA 1.00 MK 1.01 ! NL 1.52 OV 0.93 PT 1.00 SC 1.02 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! FT 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.00* FT 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FT 0.00* KP 0.00* LA 0.00* PT 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! FD 0.95 FT 1.00* LA 0.96 NL 1.26 ! OV 0.94 PT 1.00 SC 1.07 ! ! Telescope amplitude corrections in sub-array 5: ! FD 0.95 FT 1.00* LA 0.97 MK 1.02 ! NL 1.55 OV 0.90 PT 1.00 SC 1.05 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! FT 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.00* FT 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FT 0.00* KP 0.00* LA 0.00* PT 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! FD 0.93 FT 1.00* LA 1.01 NL 1.31 ! OV 0.98 PT 0.98 SC 1.03 ! ! Telescope amplitude corrections in sub-array 5: ! FD 0.93 FT 1.00* LA 1.03 MK 0.99 ! NL 1.59 OV 0.94 PT 0.99 SC 1.01 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! FT 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.00* FT 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FT 0.00* KP 0.00* LA 0.00* PT 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! FD 0.98 FT 1.00* LA 1.00 NL 1.34 ! OV 0.94 PT 1.02 SC 1.04 ! ! Telescope amplitude corrections in sub-array 5: ! FD 0.98 FT 1.00* LA 1.01 MK 1.03 ! NL 1.62 OV 0.92 PT 1.01 SC 1.02 ! ! ! Fit before self-cal, rms=0.391573Jy sigma=9.104247 ! Fit after self-cal, rms=0.245182Jy sigma=4.723061 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.8164 mas, bmaj=2.614 mas, bpa=1.316 degrees ! Estimated noise=0.827357 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 (1.4, 2.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0122528 Jy ! Component: 100 - total flux cleaned = 0.00427032 Jy ! Component: 150 - total flux cleaned = -0.000215866 Jy ! Component: 200 - total flux cleaned = -0.00119433 Jy ! Total flux subtracted in 200 components = -0.00119433 Jy ! Clean residual min=-0.030116 max=0.027121 Jy/beam ! Clean residual mean=-0.000060 rms=0.006774 Jy/beam ! Combined flux in latest and established models = 2.6555 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 38 components and 2.6555 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.208874Jy sigma=4.124135 ! Fit after self-cal, rms=0.202908Jy sigma=3.995718 ! 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.202908Jy sigma=3.995718 ! Fit after self-cal, rms=0.200917Jy sigma=3.949825 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.8153 mas, bmaj=2.612 mas, bpa=1.329 degrees ! Estimated noise=0.822984 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0012145 Jy ! Component: 100 - total flux cleaned = -0.00358391 Jy ! Component: 150 - total flux cleaned = -0.00428254 Jy ! Component: 200 - total flux cleaned = -0.00429044 Jy ! Total flux subtracted in 200 components = -0.00429044 Jy ! Clean residual min=-0.027947 max=0.023577 Jy/beam ! Clean residual mean=-0.000052 rms=0.005893 Jy/beam ! Combined flux in latest and established models = 2.65121 Jy selfcal ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 50 components and 2.65121 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.183433Jy sigma=3.793376 ! Fit after self-cal, rms=0.180275Jy sigma=3.732943 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.000766461 Jy ! Component: 100 - total flux cleaned = -8.40745e-05 Jy ! Component: 150 - total flux cleaned = -0.000614234 Jy ! Component: 200 - total flux cleaned = -0.000101681 Jy ! Total flux subtracted in 200 components = -0.000101681 Jy ! Clean residual min=-0.027418 max=0.021129 Jy/beam ! Clean residual mean=-0.000049 rms=0.005499 Jy/beam ! Combined flux in latest and established models = 2.65111 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 58 components and 2.65111 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.170589Jy sigma=3.645498 ! Fit after self-cal, rms=0.168550Jy sigma=3.609850 ! 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.168550Jy sigma=3.609850 ! Fit after self-cal, rms=0.166265Jy sigma=3.395575 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.8148 mas, bmaj=2.611 mas, bpa=1.355 degrees ! Estimated noise=0.821014 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00639912 Jy ! Component: 100 - total flux cleaned = -0.00592048 Jy ! Component: 150 - total flux cleaned = -0.00546549 Jy ! Component: 200 - total flux cleaned = -0.00545958 Jy ! Total flux subtracted in 200 components = -0.00545958 Jy ! Clean residual min=-0.020423 max=0.019313 Jy/beam ! Clean residual mean=-0.000028 rms=0.004261 Jy/beam ! Combined flux in latest and established models = 2.64565 Jy selfcal ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 67 components and 2.64565 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.158770Jy sigma=3.324314 ! Fit after self-cal, rms=0.157213Jy sigma=3.298604 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.000551539 Jy ! Component: 100 - total flux cleaned = 0.000223992 Jy ! Component: 150 - total flux cleaned = 0.000217117 Jy ! Component: 200 - total flux cleaned = 0.000207702 Jy ! Total flux subtracted in 200 components = 0.000207702 Jy ! Clean residual min=-0.019432 max=0.018646 Jy/beam ! Clean residual mean=-0.000027 rms=0.004019 Jy/beam ! Combined flux in latest and established models = 2.64586 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 74 components and 2.64586 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.152647Jy sigma=3.257135 ! Fit after self-cal, rms=0.151589Jy sigma=3.241188 ! 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.151589Jy sigma=3.241188 ! Fit after self-cal, rms=0.152174Jy sigma=3.230595 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.573 mas, bmaj=10.27 mas, bpa=-1.417 degrees ! Estimated noise=0.98131 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00793412 Jy ! Component: 100 - total flux cleaned = -0.00821602 Jy ! Component: 150 - total flux cleaned = -0.00800986 Jy ! Component: 200 - total flux cleaned = -0.00799841 Jy ! Total flux subtracted in 200 components = -0.00799841 Jy ! Clean residual min=-0.016483 max=0.015626 Jy/beam ! Clean residual mean=-0.000028 rms=0.004140 Jy/beam ! Combined flux in latest and established models = 2.63786 Jy selfcal ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 77 components and 2.63786 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.150593Jy sigma=3.244677 ! Fit after self-cal, rms=0.150182Jy sigma=3.234694 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.000521938 Jy ! Component: 100 - total flux cleaned = -0.000529307 Jy ! Component: 150 - total flux cleaned = -0.000363295 Jy ! Component: 200 - total flux cleaned = -0.000518783 Jy ! Total flux subtracted in 200 components = -0.000518783 Jy ! Clean residual min=-0.016177 max=0.015523 Jy/beam ! Clean residual mean=-0.000026 rms=0.004100 Jy/beam ! Combined flux in latest and established models = 2.63734 Jy ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 78 components and 2.63734 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.149450Jy sigma=3.254748 ! Fit after self-cal, rms=0.149344Jy sigma=3.251020 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "=== Finished amplitude self-cal for VLBA stations ===" ! === Finished amplitude self-cal for VLBA stations === unflag * ! Unflagging all baselines in the currently selected channels. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. integer lun lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_selfant_process.csh /tmp/ff.ff" ![@/tmp/ff.ff] selfant 1:SC,true ! 1:SC self-cal status=fixed weight=1 selfant 2:FD,true ! 2:FD self-cal status=fixed weight=1 selfant 2:SC,true ! 2:SC self-cal status=fixed weight=1 selfant 3:KP,true ! 3:KP self-cal status=fixed weight=1 selfant 3:LA,true ! 3:LA self-cal status=fixed weight=1 selfant 3:PT,true ! 3:PT self-cal status=fixed weight=1 selfant 4:FD,true ! 4:FD self-cal status=fixed weight=1 selfant 4:LA,true ! 4:LA self-cal status=fixed weight=1 selfant 4:NL,true ! 4:NL self-cal status=fixed weight=1 selfant 4:OV,true ! 4:OV self-cal status=fixed weight=1 selfant 4:PT,true ! 4:PT self-cal status=fixed weight=1 selfant 4:SC,true ! 4:SC self-cal status=fixed weight=1 selfant 5:FD,true ! 5:FD self-cal status=fixed weight=1 selfant 5:LA,true ! 5:LA self-cal status=fixed weight=1 selfant 5:MK,true ! 5:MK self-cal status=fixed weight=1 selfant 5:NL,true ! 5:NL self-cal status=fixed weight=1 selfant 5:OV,true ! 5:OV self-cal status=fixed weight=1 selfant 5:PT,true ! 5:PT self-cal status=fixed weight=1 selfant 5:SC,true ! 5:SC self-cal status=fixed weight=1 ![Exited script file: /tmp/ff.ff] selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 5 telescope corrections were flagged in sub-array 1. ! A total of 3 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 5 telescope corrections were flagged in sub-array 1. ! A total of 3 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 5 telescope corrections were flagged in sub-array 1. ! A total of 3 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 5 telescope corrections were flagged in sub-array 1. ! A total of 3 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.685527Jy sigma=22.511017 ! Fit after self-cal, rms=0.685527Jy sigma=22.511017 selfant "",false ! 1:FT self-cal status=correctable weight=1 ! 1:SC self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:FT self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 ! 3:FT self-cal status=correctable weight=1 ! 3:KP self-cal status=correctable weight=1 ! 3:LA self-cal status=correctable weight=1 ! 3:PT self-cal status=correctable weight=1 ! 4:FD self-cal status=correctable weight=1 ! 4:FT self-cal status=correctable weight=1 ! 4:LA self-cal status=correctable weight=1 ! 4:NL self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 4:PT self-cal status=correctable weight=1 ! 4:SC self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:FT self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:MK self-cal status=correctable weight=1 ! 5:NL self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 5:SC self-cal status=correctable weight=1 selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 5 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 28 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 5 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 28 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 5 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 28 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 5 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 28 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.685527Jy sigma=22.511017 ! Fit after self-cal, rms=0.149506Jy sigma=3.243311 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.6558 mas, bmaj=2.226 mas, bpa=2.134 degrees ! Estimated noise=1.0679 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.00507396 Jy ! Component: 100 - total flux cleaned = 0.00519511 Jy ! Component: 150 - total flux cleaned = 0.00567644 Jy ! Component: 200 - total flux cleaned = 0.00603453 Jy ! Total flux subtracted in 200 components = 0.00603453 Jy ! Clean residual min=-0.019839 max=0.020194 Jy/beam ! Clean residual mean=-0.000026 rms=0.004632 Jy/beam ! Combined flux in latest and established models = 2.64337 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 91 components and 2.64337 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.145898Jy sigma=3.188662 ! Fit after self-cal, rms=0.144722Jy sigma=3.168026 ! Inverting map uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.815 mas, bmaj=2.611 mas, bpa=1.343 degrees ! Estimated noise=0.820331 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 (11.6, 7.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.008911 Jy ! Component: 100 - total flux cleaned = 0.00805521 Jy ! Component: 150 - total flux cleaned = 0.00770312 Jy ! Component: 200 - total flux cleaned = 0.00704061 Jy ! Component: 250 - total flux cleaned = 0.00613373 Jy ! Component: 300 - total flux cleaned = 0.00529955 Jy ! Component: 350 - total flux cleaned = 0.00451965 Jy ! Component: 400 - total flux cleaned = 0.0040296 Jy ! Total flux subtracted in 400 components = 0.0040296 Jy ! Clean residual min=-0.014263 max=0.011909 Jy/beam ! Clean residual mean=-0.000028 rms=0.003027 Jy/beam ! Combined flux in latest and established models = 2.6474 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 118 components and 2.6474 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.138704Jy sigma=3.020378 ! Fit after self-cal, rms=0.138427Jy sigma=3.011451 ! 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.00118065 Jy ! Component: 100 - total flux cleaned = 0.00119427 Jy ! Component: 150 - total flux cleaned = 0.00162926 Jy ! Component: 200 - total flux cleaned = 0.00162917 Jy ! Component: 250 - total flux cleaned = 0.00162186 Jy ! Component: 300 - total flux cleaned = 0.00143743 Jy ! Component: 350 - total flux cleaned = 0.00144367 Jy ! Component: 400 - total flux cleaned = 0.00144443 Jy ! Total flux subtracted in 400 components = 0.00144443 Jy ! Clean residual min=-0.013691 max=0.011084 Jy/beam ! Clean residual mean=-0.000027 rms=0.002898 Jy/beam ! Combined flux in latest and established models = 2.64885 Jy ! Performing phase self-cal ! Adding 41 model components to the UV plane model. ! The established model now contains 138 components and 2.64885 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.135853Jy sigma=2.981144 ! Fit after self-cal, rms=0.135618Jy sigma=2.977104 ! Inverting map print "=== Finished natural weighting clean for al sites ===" ! === Finished natural weighting clean for al sites === unflag * ! Unflagging all baselines in the currently selected channels. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 36 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 36 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 36 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 36 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.683585Jy sigma=22.503313 ! Fit after self-cal, rms=0.135825Jy sigma=2.961648 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.8148 mas, bmaj=2.61 mas, bpa=1.348 degrees ! Estimated noise=0.819684 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00140822 Jy ! Component: 100 - total flux cleaned = 0.00396811 Jy ! Component: 150 - total flux cleaned = 0.00448966 Jy ! Component: 200 - total flux cleaned = 0.00482261 Jy ! Component: 250 - total flux cleaned = 0.00498588 Jy ! Component: 300 - total flux cleaned = 0.00530275 Jy ! Component: 350 - total flux cleaned = 0.00561081 Jy ! Component: 400 - total flux cleaned = 0.00576563 Jy ! Total flux subtracted in 400 components = 0.00576563 Jy ! Clean residual min=-0.013269 max=0.009403 Jy/beam ! Clean residual mean=-0.000026 rms=0.002697 Jy/beam ! Combined flux in latest and established models = 2.65461 Jy selfcal ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 153 components and 2.65461 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.133760Jy sigma=2.939471 ! Fit after self-cal, rms=0.133524Jy sigma=2.936279 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.000271106 Jy ! Component: 100 - total flux cleaned = 0.00102825 Jy ! Component: 150 - total flux cleaned = 0.00131243 Jy ! Component: 200 - total flux cleaned = 0.00145111 Jy ! Component: 250 - total flux cleaned = 0.00172105 Jy ! Component: 300 - total flux cleaned = 0.00171823 Jy ! Component: 350 - total flux cleaned = 0.00172023 Jy ! Component: 400 - total flux cleaned = 0.00184724 Jy ! Total flux subtracted in 400 components = 0.00184724 Jy ! Clean residual min=-0.013199 max=0.009127 Jy/beam ! Clean residual mean=-0.000027 rms=0.002622 Jy/beam ! Combined flux in latest and established models = 2.65646 Jy ! Performing phase self-cal ! Adding 41 model components to the UV plane model. ! The established model now contains 159 components and 2.65646 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.132071Jy sigma=2.923871 ! Fit after self-cal, rms=0.131871Jy sigma=2.921285 ! 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.131871Jy sigma=2.921285 ! Fit after self-cal, rms=0.131901Jy sigma=2.916027 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.8147 mas, bmaj=2.61 mas, bpa=1.347 degrees ! Estimated noise=0.819254 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000450618 Jy ! Component: 100 - total flux cleaned = 0.00188784 Jy ! Component: 150 - total flux cleaned = 0.0022803 Jy ! Component: 200 - total flux cleaned = 0.0025329 Jy ! Component: 250 - total flux cleaned = 0.00277579 Jy ! Component: 300 - total flux cleaned = 0.00289777 Jy ! Component: 350 - total flux cleaned = 0.0030176 Jy ! Component: 400 - total flux cleaned = 0.0030194 Jy ! Total flux subtracted in 400 components = 0.0030194 Jy ! Clean residual min=-0.012937 max=0.009132 Jy/beam ! Clean residual mean=-0.000025 rms=0.002526 Jy/beam ! Combined flux in latest and established models = 2.65948 Jy selfcal ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 163 components and 2.65948 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.130636Jy sigma=2.905214 ! Fit after self-cal, rms=0.130456Jy sigma=2.903195 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.000251856 Jy ! Component: 100 - total flux cleaned = 0.000841114 Jy ! Component: 150 - total flux cleaned = 0.00106885 Jy ! Component: 200 - total flux cleaned = 0.00117715 Jy ! Component: 250 - total flux cleaned = 0.00139438 Jy ! Component: 300 - total flux cleaned = 0.00150102 Jy ! Component: 350 - total flux cleaned = 0.00150129 Jy ! Component: 400 - total flux cleaned = 0.00170687 Jy ! Total flux subtracted in 400 components = 0.00170687 Jy ! Clean residual min=-0.012743 max=0.009081 Jy/beam ! Clean residual mean=-0.000025 rms=0.002479 Jy/beam ! Combined flux in latest and established models = 2.66119 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 169 components and 2.66118 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.129485Jy sigma=2.895679 ! Fit after self-cal, rms=0.129339Jy sigma=2.894055 ! 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.129339Jy sigma=2.894055 ! Fit after self-cal, rms=0.129027Jy sigma=2.877631 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.574 mas, bmaj=10.27 mas, bpa=-1.428 degrees ! Estimated noise=0.979439 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 1.66915e-05 Jy ! Component: 100 - total flux cleaned = 0.000828231 Jy ! Component: 150 - total flux cleaned = 0.0012205 Jy ! Component: 200 - total flux cleaned = 0.001645 Jy ! Component: 250 - total flux cleaned = 0.00177148 Jy ! Component: 300 - total flux cleaned = 0.00189254 Jy ! Component: 350 - total flux cleaned = 0.00194801 Jy ! Component: 400 - total flux cleaned = 0.00194731 Jy ! Total flux subtracted in 400 components = 0.00194731 Jy ! Clean residual min=-0.013236 max=0.010195 Jy/beam ! Clean residual mean=-0.000039 rms=0.002941 Jy/beam ! Combined flux in latest and established models = 2.66313 Jy selfcal ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 176 components and 2.66313 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.128856Jy sigma=2.884250 ! Fit after self-cal, rms=0.128903Jy sigma=2.882477 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.000215026 Jy ! Component: 100 - total flux cleaned = 0.000389906 Jy ! Component: 150 - total flux cleaned = 0.000505396 Jy ! Component: 200 - total flux cleaned = 0.000612924 Jy ! Component: 250 - total flux cleaned = 0.000662142 Jy ! Component: 300 - total flux cleaned = 0.000709684 Jy ! Component: 350 - total flux cleaned = 0.000758451 Jy ! Component: 400 - total flux cleaned = 0.000756751 Jy ! Total flux subtracted in 400 components = 0.000756751 Jy ! Clean residual min=-0.013383 max=0.009936 Jy/beam ! Clean residual mean=-0.000038 rms=0.002919 Jy/beam ! Combined flux in latest and established models = 2.66389 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 180 components and 2.66389 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.128836Jy sigma=2.890997 ! Fit after self-cal, rms=0.128882Jy sigma=2.890398 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "=== Finished amplitude self-cal for all stations ===" ! === Finished amplitude self-cal for all stations === restore ! Inverting map and beam ! Estimated beam: bmin=0.8147 mas, bmaj=2.609 mas, bpa=1.342 degrees ! Estimated noise=0.818911 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.8147 x 2.609 at 1.342 degrees (North through East) ! Clean map min=-0.035254 max=2.1723 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.128882Jy sigma=2.890398 ! Fit after self-cal, rms=0.121493Jy sigma=2.497112 ! Inverting map and beam ! Estimated beam: bmin=0.8148 mas, bmaj=2.609 mas, bpa=1.327 degrees ! Estimated noise=0.819063 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000976094 Jy ! Component: 100 - total flux cleaned = 0.00174519 Jy ! Component: 150 - total flux cleaned = 0.00187045 Jy ! Component: 200 - total flux cleaned = 0.00186673 Jy ! Component: 250 - total flux cleaned = 0.00186773 Jy ! Component: 300 - total flux cleaned = 0.00186433 Jy ! Component: 350 - total flux cleaned = 0.00186063 Jy ! Component: 400 - total flux cleaned = 0.00196954 Jy ! Total flux subtracted in 400 components = 0.00196954 Jy ! Clean residual min=-0.012544 max=0.008983 Jy/beam ! Clean residual mean=-0.000027 rms=0.002440 Jy/beam ! Combined flux in latest and established models = 2.66586 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 185 components and 2.66586 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.120478Jy sigma=2.478190 ! Fit after self-cal, rms=0.120301Jy sigma=2.476513 ! Inverting map ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000470242 Jy ! Component: 100 - total flux cleaned = -0.000577423 Jy ! Component: 150 - total flux cleaned = -0.000471124 Jy ! Component: 200 - total flux cleaned = -0.000269485 Jy ! Component: 250 - total flux cleaned = -0.000270831 Jy ! Component: 300 - total flux cleaned = -0.000269576 Jy ! Component: 350 - total flux cleaned = -8.11977e-05 Jy ! Component: 400 - total flux cleaned = 9.81227e-06 Jy ! Total flux subtracted in 400 components = 9.81227e-06 Jy ! Clean residual min=-0.012245 max=0.008622 Jy/beam ! Clean residual mean=-0.000026 rms=0.002399 Jy/beam ! Combined flux in latest and established models = 2.66587 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 189 components and 2.66587 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.119548Jy sigma=2.468119 ! Fit after self-cal, rms=0.119435Jy sigma=2.466948 ! 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.119435Jy sigma=2.466948 ! Fit after self-cal, rms=0.119506Jy sigma=2.466496 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.6559 mas, bmaj=2.225 mas, bpa=2.12 degrees ! Estimated noise=1.06644 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 1.78588 Jy ! Component: 100 - total flux cleaned = 2.18311 Jy ! Component: 150 - total flux cleaned = 2.32261 Jy ! Component: 200 - total flux cleaned = 2.41302 Jy ! Component: 250 - total flux cleaned = 2.47784 Jy ! Component: 300 - total flux cleaned = 2.52563 Jy ! Component: 350 - total flux cleaned = 2.561 Jy ! Component: 400 - total flux cleaned = 2.58541 Jy ! Total flux subtracted in 400 components = 2.58541 Jy ! Clean residual min=-0.026856 max=0.024242 Jy/beam ! Clean residual mean=0.000116 rms=0.005309 Jy/beam ! Combined flux in latest and established models = 2.58541 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 34 model components to the UV plane model. ! The established model now contains 34 components and 2.58541 Jy ! Inverting map and beam ! Estimated beam: bmin=0.8148 mas, bmaj=2.609 mas, bpa=1.326 degrees ! Estimated noise=0.818852 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.0268054 Jy ! Component: 100 - total flux cleaned = 0.0453606 Jy ! Component: 150 - total flux cleaned = 0.0557977 Jy ! Component: 200 - total flux cleaned = 0.0591695 Jy ! Component: 250 - total flux cleaned = 0.0613013 Jy ! Component: 300 - total flux cleaned = 0.0644025 Jy ! Component: 350 - total flux cleaned = 0.0661747 Jy ! Component: 400 - total flux cleaned = 0.0674995 Jy ! Total flux subtracted in 400 components = 0.0674995 Jy ! Clean residual min=-0.018399 max=0.015460 Jy/beam ! Clean residual mean=0.000002 rms=0.003562 Jy/beam ! Combined flux in latest and established models = 2.65291 Jy ! Adding 46 model components to the UV plane model. ! The established model now contains 65 components and 2.65291 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.575 mas, bmaj=10.28 mas, bpa=-1.435 degrees ! Estimated noise=0.979166 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00311758 Jy ! Component: 100 - total flux cleaned = 0.00619484 Jy ! Component: 150 - total flux cleaned = 0.00909409 Jy ! Component: 200 - total flux cleaned = 0.0113043 Jy ! Component: 250 - total flux cleaned = 0.0130073 Jy ! Component: 300 - total flux cleaned = 0.0138014 Jy ! Component: 350 - total flux cleaned = 0.0143633 Jy ! Component: 400 - total flux cleaned = 0.0147308 Jy ! Total flux subtracted in 400 components = 0.0147308 Jy ! Clean residual min=-0.016691 max=0.011461 Jy/beam ! Clean residual mean=-0.000047 rms=0.003807 Jy/beam ! Combined flux in latest and established models = 2.66764 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.000175469 Jy ! Component: 100 - total flux cleaned = 0.000510352 Jy ! Component: 150 - total flux cleaned = 0.00084346 Jy ! Component: 200 - total flux cleaned = 0.00100278 Jy ! Component: 250 - total flux cleaned = 0.00131448 Jy ! Component: 300 - total flux cleaned = 0.00131312 Jy ! Component: 350 - total flux cleaned = 0.00131964 Jy ! Component: 400 - total flux cleaned = 0.00145894 Jy ! Total flux subtracted in 400 components = 0.00145894 Jy ! Clean residual min=-0.015928 max=0.010658 Jy/beam ! Clean residual mean=-0.000046 rms=0.003570 Jy/beam ! Combined flux in latest and established models = 2.6691 Jy ! Adding 16 model components to the UV plane model. ! The established model now contains 79 components and 2.6691 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=0.8148 mas, bmaj=2.609 mas, bpa=1.326 degrees ! Estimated noise=0.818852 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00169796 Jy ! Component: 100 - total flux cleaned = -0.00346822 Jy ! Component: 150 - total flux cleaned = -0.00406805 Jy ! Component: 200 - total flux cleaned = -0.00434217 Jy ! Component: 250 - total flux cleaned = -0.00410673 Jy ! Component: 300 - total flux cleaned = -0.00365268 Jy ! Component: 350 - total flux cleaned = -0.00321551 Jy ! Component: 400 - total flux cleaned = -0.00260058 Jy ! Total flux subtracted in 400 components = -0.00260058 Jy ! Clean residual min=-0.014089 max=0.011145 Jy/beam ! Clean residual mean=-0.000025 rms=0.002888 Jy/beam ! Combined flux in latest and established models = 2.6665 Jy selfcal ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 96 components and 2.6665 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.123658Jy sigma=2.580024 ! Fit after self-cal, rms=0.121855Jy sigma=2.533079 wmodel /scr/pima/rdv93_uvs/J1427-4206_X_map.mod ! Writing 96 model components to file: /scr/pima/rdv93_uvs/J1427-4206_X_map.mod wobs /scr/pima/rdv93_uvs/J1427-4206_X_uvs.fits ! Applying 789 buffered edits. ! Writing UV FITS file: /scr/pima/rdv93_uvs/J1427-4206_X_uvs.fits wwins /scr/pima/rdv93_uvs/J1427-4206_X_map.win ! wwins: Wrote 4 windows to /scr/pima/rdv93_uvs/J1427-4206_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.00569252 Jy ! Component: 100 - total flux cleaned = -0.00996077 Jy ! Component: 150 - total flux cleaned = -0.0119148 Jy ! Component: 200 - total flux cleaned = -0.0140977 Jy ! Component: 250 - total flux cleaned = -0.0158557 Jy ! Component: 300 - total flux cleaned = -0.0176913 Jy ! Component: 350 - total flux cleaned = -0.0197491 Jy ! Component: 400 - total flux cleaned = -0.0214536 Jy ! Component: 450 - total flux cleaned = -0.0233816 Jy ! Component: 500 - total flux cleaned = -0.0249829 Jy ! Component: 550 - total flux cleaned = -0.0265384 Jy ! Component: 600 - total flux cleaned = -0.027925 Jy ! Component: 650 - total flux cleaned = -0.0295161 Jy ! Component: 700 - total flux cleaned = -0.0304673 Jy ! Component: 750 - total flux cleaned = -0.0311618 Jy ! Component: 800 - total flux cleaned = -0.0317307 Jy ! Component: 850 - total flux cleaned = -0.0320659 Jy ! Component: 900 - total flux cleaned = -0.0327172 Jy ! Component: 950 - total flux cleaned = -0.0334604 Jy ! Component: 1000 - total flux cleaned = -0.0338763 Jy ! Component: 1050 - total flux cleaned = -0.0342856 Jy ! Component: 1100 - total flux cleaned = -0.0343865 Jy ! Component: 1150 - total flux cleaned = -0.0348831 Jy ! Component: 1200 - total flux cleaned = -0.0352713 Jy ! Component: 1250 - total flux cleaned = -0.0354618 Jy ! Component: 1300 - total flux cleaned = -0.0357451 Jy ! Component: 1350 - total flux cleaned = -0.0360243 Jy ! Component: 1400 - total flux cleaned = -0.0362987 Jy ! Component: 1450 - total flux cleaned = -0.0363891 Jy ! Component: 1500 - total flux cleaned = -0.0366577 Jy ! Component: 1550 - total flux cleaned = -0.0370095 Jy ! Component: 1600 - total flux cleaned = -0.0372689 Jy ! Component: 1650 - total flux cleaned = -0.0377808 Jy ! Component: 1700 - total flux cleaned = -0.0379487 Jy ! Component: 1750 - total flux cleaned = -0.0383652 Jy ! Component: 1800 - total flux cleaned = -0.0388563 Jy ! Component: 1850 - total flux cleaned = -0.0392617 Jy ! Component: 1900 - total flux cleaned = -0.0394219 Jy ! Component: 1950 - total flux cleaned = -0.0398975 Jy ! Component: 2000 - total flux cleaned = -0.0401325 Jy ! Component: 2050 - total flux cleaned = -0.0405973 Jy ! Component: 2100 - total flux cleaned = -0.0404453 Jy ! Component: 2150 - total flux cleaned = -0.0406729 Jy ! Component: 2200 - total flux cleaned = -0.0407481 Jy ! Component: 2250 - total flux cleaned = -0.0407489 Jy ! Component: 2300 - total flux cleaned = -0.0409691 Jy ! Component: 2350 - total flux cleaned = -0.0408964 Jy ! Component: 2400 - total flux cleaned = -0.0407516 Jy ! Component: 2450 - total flux cleaned = -0.0407518 Jy ! Component: 2500 - total flux cleaned = -0.0407518 Jy ! Component: 2550 - total flux cleaned = -0.0406818 Jy ! Component: 2600 - total flux cleaned = -0.0408209 Jy ! Component: 2650 - total flux cleaned = -0.0406835 Jy ! Component: 2700 - total flux cleaned = -0.0409581 Jy ! Component: 2750 - total flux cleaned = -0.0408905 Jy ! Component: 2800 - total flux cleaned = -0.0411604 Jy ! Component: 2850 - total flux cleaned = -0.0410933 Jy ! Component: 2900 - total flux cleaned = -0.0412928 Jy ! Component: 2950 - total flux cleaned = -0.041227 Jy ! Component: 3000 - total flux cleaned = -0.0413587 Jy ! Component: 3050 - total flux cleaned = -0.0415539 Jy ! Component: 3100 - total flux cleaned = -0.041683 Jy ! Component: 3150 - total flux cleaned = -0.0416833 Jy ! Component: 3200 - total flux cleaned = -0.0416838 Jy ! Component: 3250 - total flux cleaned = -0.0420014 Jy ! Component: 3300 - total flux cleaned = -0.0421908 Jy ! Component: 3350 - total flux cleaned = -0.042316 Jy ! Component: 3400 - total flux cleaned = -0.0425028 Jy ! Component: 3450 - total flux cleaned = -0.0426885 Jy ! Component: 3500 - total flux cleaned = -0.0427506 Jy ! Component: 3550 - total flux cleaned = -0.0429348 Jy ! Component: 3600 - total flux cleaned = -0.042996 Jy ! Component: 3650 - total flux cleaned = -0.0431781 Jy ! Component: 3700 - total flux cleaned = -0.0431178 Jy ! Component: 3750 - total flux cleaned = -0.0432986 Jy ! Component: 3800 - total flux cleaned = -0.0435387 Jy ! Component: 3850 - total flux cleaned = -0.0435384 Jy ! Component: 3900 - total flux cleaned = -0.0437166 Jy ! Component: 3950 - total flux cleaned = -0.0437165 Jy ! Component: 4000 - total flux cleaned = -0.0436572 Jy ! Component: 4050 - total flux cleaned = -0.0437164 Jy ! Total flux subtracted in 4096 components = -0.0435992 Jy ! Clean residual min=-0.003499 max=0.003246 Jy/beam ! Clean residual mean=-0.000006 rms=0.001019 Jy/beam ! Combined flux in latest and established models = 2.6229 Jy keep ! Adding 243 model components to the UV plane model. ! The established model now contains 339 components and 2.6229 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv93_uvs/J1427-4206_X\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.8148 x 2.609 at 1.326 degrees (North through East) ! Clean map min=-0.019854 max=2.2091 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv93_uvs/J1427-4206_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= 2.209 rms= 0.00098 Jy quit ! Quitting program ! Log file /scr/pima/rdv93_uvs/J1427-4206_X_uvs.log closed on Tue Oct 8 15:48:16 2013