! Started logfile: /scr/pima/rdv62_uvs/J1330+2509_S_uvs.log on Fri Nov 22 22:37:21 2013 obs /scr/pima/rdv62_uvs/J1330+2509_S_uva.fits ! Reading UV FITS file: /scr/pima/rdv62_uvs/J1330+2509_S_uva.fits ! AN table 1: 38 integrations on 45 of 45 possible baselines. ! AN table 2: 42 integrations on 55 of 55 possible baselines. ! Apparent sampling: 0.699751 visibilities/baseline/integration-bin. ! Found source: J1330+2509 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 2.23299e+09 8e+06 1 8e+06 ! 02 2 2.26299e+09 8e+06 1 8e+06 ! 03 3 2.35299e+09 8e+06 1 8e+06 ! 04 4 2.37299e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 11252 visibilities. select RR ! Selecting polarization: RR, channels: 1..4 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 print "field_size =",field_size ! field_size = 1024 print "field_cell =",field_cell ! field_cell = 0.4 print "taper_size =",taper_size ! taper_size = 20 ![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv62_uvs/J1330+2509_S] integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 800000 float soltime2 soltime2 = 15 float soltime3 soltime3 = 1 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float old_peak float new_peak float flux_cutoff selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:HN 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 ! 1:SC self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:HN self-cal status=correctable weight=1 ! 2:KK self-cal status=correctable weight=1 ! 2:KP self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:MK self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:OV self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 system "rm /tmp/ff.ff" mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.400x0.400 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 17 telescope corrections were flagged in sub-array 1. ! A total of 100 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 2. ! A total of 17 telescope corrections were flagged in sub-array 1. ! A total of 100 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 3. ! A total of 17 telescope corrections were flagged in sub-array 1. ! A total of 100 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 4. ! A total of 17 telescope corrections were flagged in sub-array 1. ! A total of 100 telescope corrections were flagged in sub-array 2. ! ! Fit before self-cal, rms=1.344756Jy sigma=23.835597 ! Fit after self-cal, rms=0.938472Jy sigma=13.397416 ! 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 2:KK ! Flagging baselines of antenna 2:KK in the currently selected channels. ![Exited script file: /tmp/ff.ff] uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=7.749 mas, bmaj=11.26 mas, bpa=-24.54 degrees ! Estimated noise=3.26425 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 1.14726 Jy ! Component: 100 - total flux cleaned = 1.53755 Jy ! Total flux subtracted in 100 components = 1.53755 Jy ! Clean residual min=-0.295076 max=0.296262 Jy/beam ! Clean residual mean=-0.014497 rms=0.115672 Jy/beam ! Combined flux in latest and established models = 1.53755 Jy ! Performing phase self-cal ! Adding 92 model components to the UV plane model. ! The established model now contains 92 components and 1.53755 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.439316Jy sigma=9.430210 ! Fit after self-cal, rms=0.437536Jy sigma=9.424902 ! 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=9.774 mas, bmaj=12.18 mas, bpa=-21.92 degrees ! Estimated noise=1.30567 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.180675 Jy ! Component: 100 - total flux cleaned = 0.308676 Jy ! Component: 150 - total flux cleaned = 0.3707 Jy ! Component: 200 - total flux cleaned = 0.412542 Jy ! Total flux subtracted in 200 components = 0.412542 Jy ! Clean residual min=-0.102537 max=0.101733 Jy/beam ! Clean residual mean=-0.004662 rms=0.048413 Jy/beam ! Combined flux in latest and established models = 1.95009 Jy ! Performing phase self-cal ! Adding 181 model components to the UV plane model. ! The established model now contains 269 components and 1.95009 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.271048Jy sigma=4.723557 ! Fit after self-cal, rms=0.268155Jy sigma=4.693286 ! 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 1.66 FD 0.97 HN 0.72 KP 1.20 ! LA 1.04 MK 1.50 NL 1.34 OV 0.77 ! PT 0.97 SC 2.42 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.04 FD 1.03 HN 0.67 KK 1.00* ! KP 1.05 LA 1.07 MK 1.83 NL 1.11 ! OV 1.24 PT 0.92 SC 0.96 ! ! ! Correcting IF 2. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.62 FD 1.04 HN 0.99 KP 1.13 ! LA 0.97 MK 1.53 NL 1.23 OV 0.85 ! PT 0.97 SC 1.84 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.06 FD 1.09 HN 0.78 KK 1.00* ! KP 1.07 LA 1.00 MK 2.33 NL 1.18 ! OV 1.29 PT 0.92 SC 1.37 ! ! ! Correcting IF 3. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.48 FD 1.06 HN 1.38 KP 1.15 ! LA 0.84 MK 1.72 NL 1.17 OV 0.95 ! PT 0.98 SC 0.96 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.15 FD 1.09 HN 0.84 KK 1.00* ! KP 1.06 LA 0.84 MK 2.34 NL 1.30 ! OV 1.26 PT 0.96 SC 1.73 ! ! ! Correcting IF 4. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.28 FD 1.04 HN 1.36 KP 1.06 ! LA 0.71 MK 1.69 NL 1.10 OV 0.82 ! PT 0.93 SC 0.86 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.03 FD 1.05 HN 0.83 KK 1.00* ! KP 1.00 LA 0.71 MK 3.02 NL 1.19 ! OV 1.05 PT 0.91 SC 1.98 ! ! ! Fit before self-cal, rms=0.268155Jy sigma=4.693286 ! Fit after self-cal, rms=0.221548Jy sigma=1.665496 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=10.9 mas, bmaj=12.94 mas, bpa=-29.13 degrees ! Estimated noise=1.3015 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.0144356 Jy ! Component: 100 - total flux cleaned = -0.0219222 Jy ! Component: 150 - total flux cleaned = -0.034461 Jy ! Component: 200 - total flux cleaned = -0.0360406 Jy ! Total flux subtracted in 200 components = -0.0360406 Jy ! Clean residual min=-0.027955 max=0.028027 Jy/beam ! Clean residual mean=0.000003 rms=0.011709 Jy/beam ! Combined flux in latest and established models = 1.91405 Jy ! Performing phase self-cal ! Adding 95 model components to the UV plane model. ! The established model now contains 364 components and 1.91405 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.195739Jy sigma=1.105204 ! Fit after self-cal, rms=0.194992Jy sigma=1.091073 ! 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.194992Jy sigma=1.091073 ! Fit after self-cal, rms=0.198396Jy sigma=1.082003 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=11.03 mas, bmaj=12.98 mas, bpa=-30.51 degrees ! Estimated noise=1.30426 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00433996 Jy ! Component: 100 - total flux cleaned = -0.00600096 Jy ! Component: 150 - total flux cleaned = -0.00598095 Jy ! Component: 200 - total flux cleaned = -0.00938822 Jy ! Total flux subtracted in 200 components = -0.00938822 Jy ! Clean residual min=-0.018145 max=0.018196 Jy/beam ! Clean residual mean=-0.000023 rms=0.007752 Jy/beam ! Combined flux in latest and established models = 1.90466 Jy selfcal ! Performing phase self-cal ! Adding 81 model components to the UV plane model. ! The established model now contains 430 components and 1.90466 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.183578Jy sigma=0.888884 ! Fit after self-cal, rms=0.183070Jy sigma=0.879583 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.00793208 Jy ! Component: 100 - total flux cleaned = -0.0108402 Jy ! Component: 150 - total flux cleaned = -0.0126534 Jy ! Component: 200 - total flux cleaned = -0.0135529 Jy ! Total flux subtracted in 200 components = -0.0135529 Jy ! Clean residual min=-0.013526 max=0.013584 Jy/beam ! Clean residual mean=0.000002 rms=0.005767 Jy/beam ! Combined flux in latest and established models = 1.89111 Jy ! Performing phase self-cal ! Adding 93 model components to the UV plane model. ! The established model now contains 507 components and 1.89111 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.173510Jy sigma=0.793607 ! Fit after self-cal, rms=0.173367Jy sigma=0.790107 ! 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.173367Jy sigma=0.790107 ! Fit after self-cal, rms=0.176496Jy sigma=0.776648 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=11.16 mas, bmaj=13.01 mas, bpa=-32.11 degrees ! Estimated noise=1.30078 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00363278 Jy ! Component: 100 - total flux cleaned = -0.00365803 Jy ! Component: 150 - total flux cleaned = -0.00440773 Jy ! Component: 200 - total flux cleaned = -0.0031164 Jy ! Total flux subtracted in 200 components = -0.0031164 Jy ! Clean residual min=-0.010691 max=0.010731 Jy/beam ! Clean residual mean=-0.000039 rms=0.004499 Jy/beam ! Combined flux in latest and established models = 1.88799 Jy selfcal ! Performing phase self-cal ! Adding 67 model components to the UV plane model. ! The established model now contains 554 components and 1.888 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.169293Jy sigma=0.731219 ! Fit after self-cal, rms=0.169281Jy sigma=0.729840 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.00326407 Jy ! Component: 100 - total flux cleaned = -0.00621354 Jy ! Component: 150 - total flux cleaned = -0.00335667 Jy ! Component: 200 - total flux cleaned = -0.00284399 Jy ! Total flux subtracted in 200 components = -0.00284399 Jy ! Clean residual min=-0.008788 max=0.008830 Jy/beam ! Clean residual mean=-0.000026 rms=0.003724 Jy/beam ! Combined flux in latest and established models = 1.88515 Jy ! Performing phase self-cal ! Adding 70 model components to the UV plane model. ! The established model now contains 600 components and 1.88515 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.164501Jy sigma=0.703481 ! Fit after self-cal, rms=0.164521Jy sigma=0.702760 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime2 ! Performing amp+phase self-cal over 15 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.164521Jy sigma=0.702760 ! Fit after self-cal, rms=0.165169Jy sigma=0.693049 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=14.4 mas, bmaj=16.61 mas, bpa=-68.46 degrees ! Estimated noise=0.895195 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00162254 Jy ! Component: 100 - total flux cleaned = -0.00163151 Jy ! Component: 150 - total flux cleaned = -0.00341883 Jy ! Component: 200 - total flux cleaned = -0.00709612 Jy ! Total flux subtracted in 200 components = -0.00709612 Jy ! Clean residual min=-0.006668 max=0.006666 Jy/beam ! Clean residual mean=-0.000013 rms=0.002830 Jy/beam ! Combined flux in latest and established models = 1.87805 Jy selfcal ! Performing phase self-cal ! Adding 65 model components to the UV plane model. ! The established model now contains 652 components and 1.87805 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.162721Jy sigma=0.674192 ! Fit after self-cal, rms=0.162723Jy sigma=0.673581 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.00362343 Jy ! Component: 100 - total flux cleaned = -0.00512408 Jy ! Component: 150 - total flux cleaned = -0.00869975 Jy ! Component: 200 - total flux cleaned = -0.0117836 Jy ! Total flux subtracted in 200 components = -0.0117836 Jy ! Clean residual min=-0.005637 max=0.005628 Jy/beam ! Clean residual mean=-0.000015 rms=0.002416 Jy/beam ! Combined flux in latest and established models = 1.86627 Jy ! Performing phase self-cal ! Adding 80 model components to the UV plane model. ! The established model now contains 707 components and 1.86627 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.160848Jy sigma=0.662814 ! Fit after self-cal, rms=0.160839Jy sigma=0.662526 ! 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. ! Applying 2560 buffered edits. 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:HN,true ! 1:HN 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 1:SC,true ! 1:SC self-cal status=fixed weight=1 selfant 2:BR,true ! 2:BR self-cal status=fixed weight=1 selfant 2:FD,true ! 2:FD self-cal status=fixed weight=1 selfant 2:HN,true ! 2:HN self-cal status=fixed weight=1 selfant 2:KP,true ! 2:KP self-cal status=fixed weight=1 selfant 2:LA,true ! 2:LA self-cal status=fixed weight=1 selfant 2:MK,true ! 2:MK self-cal status=fixed weight=1 selfant 2:NL,true ! 2:NL self-cal status=fixed weight=1 selfant 2:OV,true ! 2:OV self-cal status=fixed weight=1 selfant 2:PT,true ! 2:PT self-cal status=fixed weight=1 selfant 2:SC,true ! 2: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 42 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 2. ! A total of 42 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 3. ! A total of 42 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 4. ! A total of 42 telescope corrections were flagged in sub-array 2. ! ! Fit before self-cal, rms=0.161339Jy sigma=0.661864 ! Fit after self-cal, rms=0.161339Jy sigma=0.661864 selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:HN 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 ! 1:SC self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:HN self-cal status=correctable weight=1 ! 2:KK self-cal status=correctable weight=1 ! 2:KP self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:MK self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:OV self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2: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 20 telescope corrections were flagged in sub-array 1. ! A total of 58 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 2. ! A total of 20 telescope corrections were flagged in sub-array 1. ! A total of 58 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 3. ! A total of 20 telescope corrections were flagged in sub-array 1. ! A total of 58 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 4. ! A total of 20 telescope corrections were flagged in sub-array 1. ! A total of 58 telescope corrections were flagged in sub-array 2. ! ! Fit before self-cal, rms=0.161339Jy sigma=0.661864 ! Fit after self-cal, rms=0.161886Jy sigma=0.661719 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=9.187 mas, bmaj=12.09 mas, bpa=-32.31 degrees ! Estimated noise=3.37118 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.00092133 Jy ! Component: 100 - total flux cleaned = 0.00118412 Jy ! Component: 150 - total flux cleaned = 0.00365699 Jy ! Component: 200 - total flux cleaned = 0.00601919 Jy ! Total flux subtracted in 200 components = 0.00601919 Jy ! Clean residual min=-0.007600 max=0.007569 Jy/beam ! Clean residual mean=0.000055 rms=0.002770 Jy/beam ! Combined flux in latest and established models = 1.87229 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 749 components and 1.87229 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.158430Jy sigma=0.654427 ! Fit after self-cal, rms=0.158231Jy sigma=0.653762 ! 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=11.19 mas, bmaj=13.02 mas, bpa=-32.7 degrees ! Estimated noise=1.29637 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.00426947 Jy ! Component: 100 - total flux cleaned = -0.00680825 Jy ! Component: 150 - total flux cleaned = -0.0088523 Jy ! Component: 200 - total flux cleaned = -0.00980725 Jy ! Component: 250 - total flux cleaned = -0.0110431 Jy ! Component: 300 - total flux cleaned = -0.0122365 Jy ! Component: 350 - total flux cleaned = -0.0119525 Jy ! Component: 400 - total flux cleaned = -0.0142373 Jy ! Total flux subtracted in 400 components = -0.0142373 Jy ! Clean residual min=-0.004685 max=0.004670 Jy/beam ! Clean residual mean=0.000038 rms=0.001943 Jy/beam ! Combined flux in latest and established models = 1.85805 Jy ! Performing phase self-cal ! Adding 83 model components to the UV plane model. ! The established model now contains 824 components and 1.85805 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.153648Jy sigma=0.639565 ! Fit after self-cal, rms=0.153581Jy sigma=0.639276 ! 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.00223127 Jy ! Component: 100 - total flux cleaned = -0.0038712 Jy ! Component: 150 - total flux cleaned = -0.00574159 Jy ! Component: 200 - total flux cleaned = -0.00757263 Jy ! Component: 250 - total flux cleaned = -0.0086045 Jy ! Component: 300 - total flux cleaned = -0.0096181 Jy ! Component: 350 - total flux cleaned = -0.00987156 Jy ! Component: 400 - total flux cleaned = -0.00938111 Jy ! Total flux subtracted in 400 components = -0.00938111 Jy ! Clean residual min=-0.004081 max=0.004070 Jy/beam ! Clean residual mean=0.000048 rms=0.001730 Jy/beam ! Combined flux in latest and established models = 1.84867 Jy ! Performing phase self-cal ! Adding 85 model components to the UV plane model. ! The established model now contains 863 components and 1.84867 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.150648Jy sigma=0.631048 ! Fit after self-cal, rms=0.150666Jy sigma=0.630816 ! 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. ! Applying 2560 buffered edits. ! Applying 2560 buffered edits. 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 20 telescope corrections were flagged in sub-array 1. ! A total of 100 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 2. ! A total of 20 telescope corrections were flagged in sub-array 1. ! A total of 100 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 3. ! A total of 20 telescope corrections were flagged in sub-array 1. ! A total of 100 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 4. ! A total of 20 telescope corrections were flagged in sub-array 1. ! A total of 100 telescope corrections were flagged in sub-array 2. ! ! Fit before self-cal, rms=0.151247Jy sigma=0.630179 ! Fit after self-cal, rms=0.149156Jy sigma=0.630058 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=11.19 mas, bmaj=13 mas, bpa=-32.32 degrees ! Estimated noise=1.29574 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00308958 Jy ! Component: 100 - total flux cleaned = -0.00477118 Jy ! Component: 150 - total flux cleaned = -0.00523464 Jy ! Component: 200 - total flux cleaned = -0.00658408 Jy ! Component: 250 - total flux cleaned = -0.00658143 Jy ! Component: 300 - total flux cleaned = -0.00701863 Jy ! Component: 350 - total flux cleaned = -0.00636696 Jy ! Component: 400 - total flux cleaned = -0.00615278 Jy ! Total flux subtracted in 400 components = -0.00615278 Jy ! Clean residual min=-0.003587 max=0.003562 Jy/beam ! Clean residual mean=0.000038 rms=0.001520 Jy/beam ! Combined flux in latest and established models = 1.84252 Jy selfcal ! Performing phase self-cal ! Adding 89 model components to the UV plane model. ! The established model now contains 915 components and 1.84252 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.146696Jy sigma=0.623434 ! Fit after self-cal, rms=0.146687Jy sigma=0.623183 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (-68.8, -21.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00127565 Jy ! Component: 100 - total flux cleaned = -0.00228268 Jy ! Component: 150 - total flux cleaned = -0.00350003 Jy ! Component: 200 - total flux cleaned = -0.00420857 Jy ! Component: 250 - total flux cleaned = -0.00474352 Jy ! Component: 300 - total flux cleaned = -0.00537892 Jy ! Component: 350 - total flux cleaned = -0.00586464 Jy ! Component: 400 - total flux cleaned = -0.00632655 Jy ! Total flux subtracted in 400 components = -0.00632655 Jy ! Clean residual min=-0.003946 max=0.004583 Jy/beam ! Clean residual mean=0.000019 rms=0.001414 Jy/beam ! Combined flux in latest and established models = 1.83619 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 959 components and 1.83619 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.145154Jy sigma=0.621644 ! Fit after self-cal, rms=0.145125Jy sigma=0.621176 ! Inverting map ! Added new window around map position (30.8, 59.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -8.8644e-05 Jy ! Component: 100 - total flux cleaned = -0.000992287 Jy ! Component: 150 - total flux cleaned = -0.00198014 Jy ! Component: 200 - total flux cleaned = -0.00274995 Jy ! Component: 250 - total flux cleaned = -0.00349761 Jy ! Component: 300 - total flux cleaned = -0.0043739 Jy ! Component: 350 - total flux cleaned = -0.00536638 Jy ! Component: 400 - total flux cleaned = -0.00605117 Jy ! Total flux subtracted in 400 components = -0.00605117 Jy ! Clean residual min=-0.004263 max=0.003960 Jy/beam ! Clean residual mean=0.000013 rms=0.001283 Jy/beam ! Combined flux in latest and established models = 1.83014 Jy ! Performing phase self-cal ! Adding 66 model components to the UV plane model. ! The established model now contains 1012 components and 1.83014 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.143718Jy sigma=0.617798 ! Fit after self-cal, rms=0.143705Jy sigma=0.617627 ! 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.143705Jy sigma=0.617627 ! Fit after self-cal, rms=0.143528Jy sigma=0.611265 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=11.21 mas, bmaj=12.99 mas, bpa=-32.27 degrees ! Estimated noise=1.29685 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000169912 Jy ! Component: 100 - total flux cleaned = -0.000527426 Jy ! Component: 150 - total flux cleaned = -0.00133949 Jy ! Component: 200 - total flux cleaned = -0.00212168 Jy ! Component: 250 - total flux cleaned = -0.00300957 Jy ! Component: 300 - total flux cleaned = -0.00399332 Jy ! Component: 350 - total flux cleaned = -0.0047124 Jy ! Component: 400 - total flux cleaned = -0.00541548 Jy ! Total flux subtracted in 400 components = -0.00541548 Jy ! Clean residual min=-0.004154 max=0.003606 Jy/beam ! Clean residual mean=0.000007 rms=0.001192 Jy/beam ! Combined flux in latest and established models = 1.82473 Jy selfcal ! Performing phase self-cal ! Adding 62 model components to the UV plane model. ! The established model now contains 1046 components and 1.82473 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.142481Jy sigma=0.609298 ! Fit after self-cal, rms=0.142465Jy sigma=0.609160 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.000378022 Jy ! Component: 100 - total flux cleaned = -0.000949377 Jy ! Component: 150 - total flux cleaned = -0.00139509 Jy ! Component: 200 - total flux cleaned = -0.00216075 Jy ! Component: 250 - total flux cleaned = -0.00248244 Jy ! Component: 300 - total flux cleaned = -0.00300836 Jy ! Component: 350 - total flux cleaned = -0.00321119 Jy ! Component: 400 - total flux cleaned = -0.00341278 Jy ! Total flux subtracted in 400 components = -0.00341278 Jy ! Clean residual min=-0.004151 max=0.003624 Jy/beam ! Clean residual mean=0.000002 rms=0.001139 Jy/beam ! Combined flux in latest and established models = 1.82131 Jy ! Performing phase self-cal ! Adding 72 model components to the UV plane model. ! The established model now contains 1081 components and 1.82131 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.141614Jy sigma=0.608005 ! Fit after self-cal, rms=0.141602Jy sigma=0.607901 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime3 ! Performing amp+phase self-cal over 1 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.141602Jy sigma=0.607901 ! Fit after self-cal, rms=0.145213Jy sigma=0.539463 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=14.44 mas, bmaj=16.61 mas, bpa=-69.22 degrees ! Estimated noise=0.895313 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000244487 Jy ! Component: 100 - total flux cleaned = -0.000142 Jy ! Component: 150 - total flux cleaned = -0.000509769 Jy ! Component: 200 - total flux cleaned = -0.000954709 Jy ! Component: 250 - total flux cleaned = -0.00129505 Jy ! Component: 300 - total flux cleaned = -0.00170541 Jy ! Component: 350 - total flux cleaned = -0.00210359 Jy ! Component: 400 - total flux cleaned = -0.00249289 Jy ! Total flux subtracted in 400 components = -0.00249289 Jy ! Clean residual min=-0.003433 max=0.002839 Jy/beam ! Clean residual mean=-0.000002 rms=0.000907 Jy/beam ! Combined flux in latest and established models = 1.81882 Jy selfcal ! Performing phase self-cal ! Adding 58 model components to the UV plane model. ! The established model now contains 1127 components and 1.81882 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.144876Jy sigma=0.537902 ! Fit after self-cal, rms=0.144796Jy sigma=0.537833 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.000731926 Jy ! Component: 100 - total flux cleaned = -0.00111121 Jy ! Component: 150 - total flux cleaned = -0.00155228 Jy ! Component: 200 - total flux cleaned = -0.00191282 Jy ! Component: 250 - total flux cleaned = -0.00247848 Jy ! Component: 300 - total flux cleaned = -0.00289307 Jy ! Component: 350 - total flux cleaned = -0.0033679 Jy ! Component: 400 - total flux cleaned = -0.00376679 Jy ! Total flux subtracted in 400 components = -0.00376679 Jy ! Clean residual min=-0.003327 max=0.003052 Jy/beam ! Clean residual mean=-0.000003 rms=0.000870 Jy/beam ! Combined flux in latest and established models = 1.81505 Jy ! Performing phase self-cal ! Adding 63 model components to the UV plane model. ! The established model now contains 1164 components and 1.81505 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.144647Jy sigma=0.536872 ! Fit after self-cal, rms=0.144636Jy sigma=0.536820 ! 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=11.23 mas, bmaj=13.01 mas, bpa=-32.22 degrees ! Estimated noise=1.29845 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 11.23 x 13.01 at -32.22 degrees (North through East) ! Clean map min=-0.20466 max=0.96651 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.144636Jy sigma=0.536820 ! Fit after self-cal, rms=0.220849Jy sigma=0.370269 ! Inverting map and beam ! Estimated beam: bmin=11.21 mas, bmaj=13 mas, bpa=-31.7 degrees ! Estimated noise=1.29839 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000539224 Jy ! Component: 100 - total flux cleaned = 0.00044263 Jy ! Component: 150 - total flux cleaned = 0.000441955 Jy ! Component: 200 - total flux cleaned = 0.000533716 Jy ! Component: 250 - total flux cleaned = 0.000535214 Jy ! Component: 300 - total flux cleaned = 0.000447535 Jy ! Component: 350 - total flux cleaned = 0.000363315 Jy ! Component: 400 - total flux cleaned = 0.00019391 Jy ! Total flux subtracted in 400 components = 0.00019391 Jy ! Clean residual min=-0.003739 max=0.003488 Jy/beam ! Clean residual mean=-0.000001 rms=0.001003 Jy/beam ! Combined flux in latest and established models = 1.81525 Jy ! Performing phase self-cal ! Adding 65 model components to the UV plane model. ! The established model now contains 1214 components and 1.81525 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.220337Jy sigma=0.368980 ! Fit after self-cal, rms=0.220179Jy sigma=0.368895 ! Inverting map ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000101735 Jy ! Component: 100 - total flux cleaned = -1.55595e-05 Jy ! Component: 150 - total flux cleaned = -0.000179724 Jy ! Component: 200 - total flux cleaned = -0.000180863 Jy ! Component: 250 - total flux cleaned = -0.00042173 Jy ! Component: 300 - total flux cleaned = -0.000500711 Jy ! Component: 350 - total flux cleaned = -0.000501287 Jy ! Component: 400 - total flux cleaned = -0.000655651 Jy ! Total flux subtracted in 400 components = -0.000655651 Jy ! Clean residual min=-0.003773 max=0.003431 Jy/beam ! Clean residual mean=-0.000007 rms=0.000974 Jy/beam ! Combined flux in latest and established models = 1.81459 Jy ! Performing phase self-cal ! Adding 56 model components to the UV plane model. ! The established model now contains 1233 components and 1.81459 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.219746Jy sigma=0.367840 ! Fit after self-cal, rms=0.219717Jy sigma=0.367747 ! 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.219717Jy sigma=0.367747 ! Fit after self-cal, rms=0.224723Jy sigma=0.367528 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=9.21 mas, bmaj=12.11 mas, bpa=-32.05 degrees ! Estimated noise=3.38704 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.578607 Jy ! Component: 100 - total flux cleaned = -0.96061 Jy ! Component: 150 - total flux cleaned = -1.26788 Jy ! Component: 200 - total flux cleaned = -1.52194 Jy ! Component: 250 - total flux cleaned = -1.74867 Jy ! Component: 300 - total flux cleaned = -1.97532 Jy ! Component: 350 - total flux cleaned = -2.17011 Jy ! Component: 400 - total flux cleaned = -2.33965 Jy ! Total flux subtracted in 400 components = -2.33965 Jy ! Clean residual min=-0.409208 max=0.770551 Jy/beam ! Clean residual mean=-0.017382 rms=0.152422 Jy/beam ! Combined flux in latest and established models = -2.33965 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 310 model components to the UV plane model. ! The established model now contains 310 components and -2.33965 Jy ! Inverting map and beam ! Estimated beam: bmin=11.22 mas, bmaj=13.01 mas, bpa=-31.78 degrees ! Estimated noise=1.2992 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam keep new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (-0.8, -0.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.722211 Jy ! Component: 100 - total flux cleaned = 0.866201 Jy ! Component: 150 - total flux cleaned = 0.948765 Jy ! Component: 200 - total flux cleaned = 1.0312 Jy ! Component: 250 - total flux cleaned = 1.1238 Jy ! Component: 300 - total flux cleaned = 1.20103 Jy ! Component: 350 - total flux cleaned = 1.26269 Jy ! Component: 400 - total flux cleaned = 1.31553 Jy ! Total flux subtracted in 400 components = 1.31553 Jy ! Clean residual min=-0.156128 max=0.114866 Jy/beam ! Clean residual mean=-0.004084 rms=0.046862 Jy/beam ! Combined flux in latest and established models = -1.02412 Jy ! Adding 219 model components to the UV plane model. ! The established model now contains 528 components and -1.02412 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=14.45 mas, bmaj=16.61 mas, bpa=-69.24 degrees ! Estimated noise=0.895845 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.045871 Jy ! Component: 100 - total flux cleaned = 0.110724 Jy ! Component: 150 - total flux cleaned = 0.161913 Jy ! Component: 200 - total flux cleaned = 0.212681 Jy ! Component: 250 - total flux cleaned = 0.257392 Jy ! Component: 300 - total flux cleaned = 0.299265 Jy ! Component: 350 - total flux cleaned = 0.33875 Jy ! Component: 400 - total flux cleaned = 0.373775 Jy ! Total flux subtracted in 400 components = 0.373775 Jy ! Clean residual min=-0.106484 max=0.080493 Jy/beam ! Clean residual mean=-0.003276 rms=0.034993 Jy/beam ! Combined flux in latest and established models = -0.650346 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.0270671 Jy ! Component: 100 - total flux cleaned = 0.0567562 Jy ! Component: 150 - total flux cleaned = 0.0887856 Jy ! Component: 200 - total flux cleaned = 0.119322 Jy ! Component: 250 - total flux cleaned = 0.146793 Jy ! Component: 300 - total flux cleaned = 0.171498 Jy ! Component: 350 - total flux cleaned = 0.190512 Jy ! Component: 400 - total flux cleaned = 0.214898 Jy ! Total flux subtracted in 400 components = 0.214898 Jy ! Clean residual min=-0.072679 max=0.060541 Jy/beam ! Clean residual mean=-0.002482 rms=0.024969 Jy/beam ! Combined flux in latest and established models = -0.435449 Jy ! Adding 293 model components to the UV plane model. ! The established model now contains 774 components and -0.435448 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=11.22 mas, bmaj=13.01 mas, bpa=-31.78 degrees ! Estimated noise=1.2992 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0163587 Jy ! Component: 100 - total flux cleaned = 0.0407072 Jy ! Component: 150 - total flux cleaned = 0.0654545 Jy ! Component: 200 - total flux cleaned = 0.0906484 Jy ! Component: 250 - total flux cleaned = 0.106008 Jy ! Component: 300 - total flux cleaned = 0.12581 Jy ! Component: 350 - total flux cleaned = 0.143715 Jy ! Component: 400 - total flux cleaned = 0.16108 Jy ! Total flux subtracted in 400 components = 0.16108 Jy ! Clean residual min=-0.056212 max=0.044462 Jy/beam ! Clean residual mean=-0.001552 rms=0.017778 Jy/beam ! Combined flux in latest and established models = -0.274368 Jy selfcal ! Performing phase self-cal ! Adding 88 model components to the UV plane model. ! The established model now contains 848 components and -0.274368 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.329410Jy sigma=1.613641 ! Fit after self-cal, rms=0.314461Jy sigma=1.340719 wmodel /scr/pima/rdv62_uvs/J1330+2509_S_map.mod ! Writing 848 model components to file: /scr/pima/rdv62_uvs/J1330+2509_S_map.mod wobs /scr/pima/rdv62_uvs/J1330+2509_S_uvs.fits ! Applying 402 buffered edits. ! Writing UV FITS file: /scr/pima/rdv62_uvs/J1330+2509_S_uvs.fits wwins /scr/pima/rdv62_uvs/J1330+2509_S_map.win ! wwins: Wrote 3 windows to /scr/pima/rdv62_uvs/J1330+2509_S_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.010826 Jy ! Component: 100 - total flux cleaned = -0.0220539 Jy ! Component: 150 - total flux cleaned = -0.0342787 Jy ! Component: 200 - total flux cleaned = -0.0460646 Jy ! Component: 250 - total flux cleaned = -0.0546871 Jy ! Component: 300 - total flux cleaned = -0.0600667 Jy ! Component: 350 - total flux cleaned = -0.0657281 Jy ! Component: 400 - total flux cleaned = -0.0712173 Jy ! Component: 450 - total flux cleaned = -0.0765359 Jy ! Component: 500 - total flux cleaned = -0.081691 Jy ! Component: 550 - total flux cleaned = -0.0867076 Jy ! Component: 600 - total flux cleaned = -0.0920175 Jy ! Component: 650 - total flux cleaned = -0.0967556 Jy ! Component: 700 - total flux cleaned = -0.102195 Jy ! Component: 750 - total flux cleaned = -0.1075 Jy ! Component: 800 - total flux cleaned = -0.111492 Jy ! Component: 850 - total flux cleaned = -0.114243 Jy ! Component: 900 - total flux cleaned = -0.117725 Jy ! Component: 950 - total flux cleaned = -0.120393 Jy ! Component: 1000 - total flux cleaned = -0.121147 Jy ! Component: 1050 - total flux cleaned = -0.121897 Jy ! Component: 1100 - total flux cleaned = -0.123396 Jy ! Component: 1150 - total flux cleaned = -0.124888 Jy ! Component: 1200 - total flux cleaned = -0.126748 Jy ! Component: 1250 - total flux cleaned = -0.128603 Jy ! Component: 1300 - total flux cleaned = -0.130082 Jy ! Component: 1350 - total flux cleaned = -0.132308 Jy ! Component: 1400 - total flux cleaned = -0.134907 Jy ! Component: 1450 - total flux cleaned = -0.137878 Jy ! Component: 1500 - total flux cleaned = -0.139734 Jy ! Component: 1550 - total flux cleaned = -0.141225 Jy ! Component: 1600 - total flux cleaned = -0.143839 Jy ! Component: 1650 - total flux cleaned = -0.144966 Jy ! Component: 1700 - total flux cleaned = -0.145341 Jy ! Component: 1750 - total flux cleaned = -0.147604 Jy ! Component: 1800 - total flux cleaned = -0.147982 Jy ! Component: 1850 - total flux cleaned = -0.149506 Jy ! Component: 1900 - total flux cleaned = -0.149889 Jy ! Component: 1950 - total flux cleaned = -0.151046 Jy ! Component: 2000 - total flux cleaned = -0.152211 Jy ! Component: 2050 - total flux cleaned = -0.154558 Jy ! Component: 2100 - total flux cleaned = -0.156131 Jy ! Component: 2150 - total flux cleaned = -0.158512 Jy ! Component: 2200 - total flux cleaned = -0.160512 Jy ! Component: 2250 - total flux cleaned = -0.162939 Jy ! Component: 2300 - total flux cleaned = -0.164567 Jy ! Component: 2350 - total flux cleaned = -0.167028 Jy ! Component: 2400 - total flux cleaned = -0.168678 Jy ! Component: 2450 - total flux cleaned = -0.17244 Jy ! Component: 2500 - total flux cleaned = -0.174979 Jy ! Component: 2550 - total flux cleaned = -0.175836 Jy ! Component: 2600 - total flux cleaned = -0.177548 Jy ! Component: 2650 - total flux cleaned = -0.180582 Jy ! Component: 2700 - total flux cleaned = -0.183199 Jy ! Component: 2750 - total flux cleaned = -0.18496 Jy ! Component: 2800 - total flux cleaned = -0.188074 Jy ! Component: 2850 - total flux cleaned = -0.190775 Jy ! Component: 2900 - total flux cleaned = -0.193044 Jy ! Component: 2950 - total flux cleaned = -0.195796 Jy ! Component: 3000 - total flux cleaned = -0.199513 Jy ! Component: 3050 - total flux cleaned = -0.200928 Jy ! Component: 3100 - total flux cleaned = -0.203776 Jy ! Component: 3150 - total flux cleaned = -0.206179 Jy ! Component: 3200 - total flux cleaned = -0.208123 Jy ! Component: 3250 - total flux cleaned = -0.21157 Jy ! Component: 3300 - total flux cleaned = -0.214059 Jy ! Component: 3350 - total flux cleaned = -0.216585 Jy ! Component: 3400 - total flux cleaned = -0.218629 Jy ! Component: 3450 - total flux cleaned = -0.222258 Jy ! Component: 3500 - total flux cleaned = -0.223834 Jy ! Component: 3550 - total flux cleaned = -0.226491 Jy ! Component: 3600 - total flux cleaned = -0.230267 Jy ! Component: 3650 - total flux cleaned = -0.231366 Jy ! Component: 3700 - total flux cleaned = -0.23358 Jy ! Component: 3750 - total flux cleaned = -0.233589 Jy ! Component: 3800 - total flux cleaned = -0.234731 Jy ! Component: 3850 - total flux cleaned = -0.235883 Jy ! Component: 3900 - total flux cleaned = -0.237059 Jy ! Component: 3950 - total flux cleaned = -0.236469 Jy ! Component: 4000 - total flux cleaned = -0.237665 Jy ! Component: 4050 - total flux cleaned = -0.238877 Jy ! Total flux subtracted in 4096 components = -0.239488 Jy ! Clean residual min=-0.036297 max=0.037429 Jy/beam ! Clean residual mean=-0.001034 rms=0.014253 Jy/beam ! Combined flux in latest and established models = -0.513856 Jy keep ! Adding 377 model components to the UV plane model. ! The established model now contains 1225 components and -0.513856 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv62_uvs/J1330+2509_S\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 11.22 x 13.01 at -31.78 degrees (North through East) ! Clean map min=-0.249 max=0.711 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv62_uvs/J1330+2509_S_map.fits float peak_flux peak_flux = peak(flux,max) invert ! Inverting map float image_peak image_peak = int(10000*peak_flux)/10000 float image_rms image_rms = int(100000*imstat(rms))/100000 print "=========================================" ! ========================================= print "The image peak= ", image_peak, " rms= ", image_rms, " Jy" ! The image peak= 0.711 rms= 0.01124 Jy quit ! Quitting program ! Log file /scr/pima/rdv62_uvs/J1330+2509_S_uvs.log closed on Fri Nov 22 22:37:30 2013