! Started logfile: /scr/pima/rdv62_uvs/J1130+3815_X_uvs.log on Fri Nov 22 22:37:21 2013 obs /scr/pima/rdv62_uvs/J1130+3815_X_uva.fits ! Reading UV FITS file: /scr/pima/rdv62_uvs/J1130+3815_X_uva.fits ! AN table 1: 23 integrations on 6 of 6 possible baselines. ! AN table 2: 15 integrations on 66 of 66 possible baselines. ! AN table 3: 3 integrations on 3 of 3 possible baselines. ! AN table 4: 6 integrations on 55 of 55 possible baselines. ! AN table 5: 7 integrations on 28 of 28 possible baselines. ! AN table 6: 9 integrations on 21 of 21 possible baselines. ! Apparent sampling: 0.458423 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1130+3815 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 8.40599e+09 8e+06 1 8e+06 ! 02 2 8.47599e+09 8e+06 1 8e+06 ! 03 3 8.79099e+09 8e+06 1 8e+06 ! 04 4 8.89599e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 3396 visibilities. select RR ! Selecting polarization: RR, channels: 1..4 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 print "field_size =",field_size ! field_size = 1024 print "field_cell =",field_cell ! field_cell = 0.2 print "taper_size =",taper_size ! taper_size = 70 ![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv62_uvs/J1130+3815_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:HH self-cal status=correctable weight=1 ! 1:NY self-cal status=correctable weight=1 ! 1:WZ self-cal status=correctable weight=1 ! 1:ZC self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:HH self-cal status=correctable weight=1 ! 2:HN self-cal status=correctable weight=1 ! 2:KP self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:NY self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 2:ZC self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:NY self-cal status=correctable weight=1 ! 3:OV self-cal status=correctable weight=1 ! 4:BR self-cal status=correctable weight=1 ! 4:FD self-cal status=correctable weight=1 ! 4:HN self-cal status=correctable weight=1 ! 4:KK self-cal status=correctable weight=1 ! 4:KP self-cal status=correctable weight=1 ! 4:LA self-cal status=correctable weight=1 ! 4:MK 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:HN self-cal status=correctable weight=1 ! 5:KK self-cal status=correctable weight=1 ! 5:KP 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:PT self-cal status=correctable weight=1 ! 6:BR self-cal status=correctable weight=1 ! 6:KK self-cal status=correctable weight=1 ! 6:KP self-cal status=correctable weight=1 ! 6:LA self-cal status=correctable weight=1 ! 6:MK self-cal status=correctable weight=1 ! 6:OV self-cal status=correctable weight=1 ! 6:PT self-cal status=correctable weight=1 system "rm /tmp/ff.ff" mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.200x0.200 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 67 telescope corrections were flagged in sub-array 1. ! A total of 88 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 16 telescope corrections were flagged in sub-array 4. ! A total of 30 telescope corrections were flagged in sub-array 5. ! A total of 23 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 67 telescope corrections were flagged in sub-array 1. ! A total of 88 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 16 telescope corrections were flagged in sub-array 4. ! A total of 30 telescope corrections were flagged in sub-array 5. ! A total of 23 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 67 telescope corrections were flagged in sub-array 1. ! A total of 88 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 16 telescope corrections were flagged in sub-array 4. ! A total of 30 telescope corrections were flagged in sub-array 5. ! A total of 23 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 71 telescope corrections were flagged in sub-array 1. ! A total of 94 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 16 telescope corrections were flagged in sub-array 4. ! A total of 30 telescope corrections were flagged in sub-array 5. ! A total of 23 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=1.455448Jy sigma=17.246909 ! Fit after self-cal, rms=0.249863Jy sigma=2.475698 ! 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:HH ! Flagging baselines of antenna 1:HH in the currently selected channels. flag 1:NY ! Flagging baselines of antenna 1:NY in the currently selected channels. flag 1:WZ ! Flagging baselines of antenna 1:WZ in the currently selected channels. flag 1:ZC ! Flagging baselines of antenna 1:ZC in the currently selected channels. flag 2:HH ! Flagging baselines of antenna 2:HH in the currently selected channels. flag 2:NY ! Flagging baselines of antenna 2:NY in the currently selected channels. flag 2:WZ ! Flagging baselines of antenna 2:WZ in the currently selected channels. flag 2:ZC ! Flagging baselines of antenna 2:ZC in the currently selected channels. flag 3:NY ! Flagging baselines of antenna 3:NY in the currently selected channels. flag 4:KK ! Flagging baselines of antenna 4:KK in the currently selected channels. flag 5:KK ! Flagging baselines of antenna 5:KK in the currently selected channels. flag 6:KK ! Flagging baselines of antenna 6: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=1.012 mas, bmaj=1.337 mas, bpa=38.81 degrees ! Estimated noise=1.77396 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 = 0.854473 Jy ! Component: 100 - total flux cleaned = 1.05603 Jy ! Total flux subtracted in 100 components = 1.05603 Jy ! Clean residual min=-0.021815 max=0.074038 Jy/beam ! Clean residual mean=-0.000352 rms=0.006723 Jy/beam ! Combined flux in latest and established models = 1.05603 Jy ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 7 components and 1.05603 Jy ! ! Correcting IF 1. ! A total of 25 telescope corrections were flagged in sub-array 1. ! A total of 34 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 3 telescope corrections were flagged in sub-array 5. ! A total of 5 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 25 telescope corrections were flagged in sub-array 1. ! A total of 34 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 3 telescope corrections were flagged in sub-array 5. ! A total of 5 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 25 telescope corrections were flagged in sub-array 1. ! A total of 34 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 3 telescope corrections were flagged in sub-array 5. ! A total of 5 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 21 telescope corrections were flagged in sub-array 1. ! A total of 28 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 3 telescope corrections were flagged in sub-array 5. ! A total of 5 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.168456Jy sigma=2.307384 ! Fit after self-cal, rms=0.168533Jy sigma=2.306778 ! Inverting map and beam ! Estimated beam: bmin=1.006 mas, bmaj=1.332 mas, bpa=38.78 degrees ! Estimated noise=1.78411 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0776928 Jy ! Component: 100 - total flux cleaned = 0.118168 Jy ! Total flux subtracted in 100 components = 0.118168 Jy ! Clean residual min=-0.019744 max=0.019613 Jy/beam ! Clean residual mean=-0.000165 rms=0.004128 Jy/beam ! Combined flux in latest and established models = 1.17419 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 17 components and 1.17419 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.115457Jy sigma=1.509927 ! Fit after self-cal, rms=0.115197Jy sigma=1.507312 ! 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=1.27 mas, bmaj=1.837 mas, bpa=41.08 degrees ! Estimated noise=1.21968 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.0334877 Jy ! Component: 100 - total flux cleaned = 0.047751 Jy ! Component: 150 - total flux cleaned = 0.0546164 Jy ! Component: 200 - total flux cleaned = 0.0576002 Jy ! Total flux subtracted in 200 components = 0.0576002 Jy ! Clean residual min=-0.014990 max=0.017813 Jy/beam ! Clean residual mean=-0.000046 rms=0.003958 Jy/beam ! Combined flux in latest and established models = 1.23179 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 28 components and 1.23179 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.107547Jy sigma=1.375540 ! Fit after self-cal, rms=0.107462Jy sigma=1.374708 ! 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: ! HH 0.00* NY 0.00* WZ 0.00* ZC 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.98 FD 0.98 HH 1.00* HN 0.99 ! KP 0.99 LA 0.94 NL 0.99 NY 1.00* ! PT 0.98 SC 1.05 WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* NY 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 0.97 FD 0.99 HN 1.05 KK 1.00* ! KP 1.00 LA 0.94 MK 0.98 NL 0.99 ! OV 0.98 PT 0.96 SC 1.01 ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! FD 1.00 HN 0.95 KK 1.00* KP 1.00 ! LA 1.01 MK 0.99 NL 0.99 PT 0.97 ! ! A total of 12 telescope corrections were flagged in sub-array 6. ! Telescope amplitude corrections in sub-array 6: ! BR 1.01 KK 1.00* KP 1.01 LA 1.04 ! MK 0.97 OV 1.00 PT 1.02 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* NY 0.00* WZ 0.00* ZC 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.99 FD 0.98 HH 1.00* HN 1.00 ! KP 0.94 LA 0.94 NL 0.96 NY 1.00* ! PT 1.04 SC 1.00 WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* NY 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 0.97 FD 0.98 HN 1.07 KK 1.00* ! KP 0.97 LA 0.93 MK 1.01 NL 0.97 ! OV 0.97 PT 1.01 SC 0.98 ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! FD 0.98 HN 0.98 KK 1.00* KP 0.96 ! LA 1.00 MK 1.02 NL 0.98 PT 1.03 ! ! A total of 12 telescope corrections were flagged in sub-array 6. ! Telescope amplitude corrections in sub-array 6: ! BR 0.99 KK 1.00* KP 0.96 LA 1.06 ! MK 0.98 OV 0.98 PT 1.11 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* NY 0.00* WZ 0.00* ZC 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.03 FD 1.03 HH 1.00* HN 1.01 ! KP 0.91 LA 1.00 NL 1.00 NY 1.00* ! PT 1.08 SC 0.98 WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* NY 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 1.02 FD 1.04 HN 1.05 KK 1.00* ! KP 0.90 LA 1.01 MK 1.00 NL 1.00 ! OV 1.03 PT 1.04 SC 0.99 ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! FD 1.04 HN 0.95 KK 1.00* KP 0.88 ! LA 1.08 MK 1.04 NL 1.01 PT 1.08 ! ! A total of 12 telescope corrections were flagged in sub-array 6. ! Telescope amplitude corrections in sub-array 6: ! BR 1.07 KK 1.00* KP 0.88 LA 1.13 ! MK 1.02 OV 1.07 PT 1.13 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* NY 0.00* WZ 0.00* ZC 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.02 FD 1.06 HH 1.00* HN 1.06 ! KP 0.96 LA 0.98 NL 1.05 NY 1.00* ! PT 1.07 SC 1.05 WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* NY 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 1.01 FD 1.05 HN 1.12 KK 1.00* ! KP 0.95 LA 0.98 MK 1.06 NL 1.04 ! OV 1.00 PT 1.04 SC 1.04 ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! FD 1.07 HN 1.03 KK 1.00* KP 0.95 ! LA 1.07 MK 1.04 NL 1.04 PT 1.06 ! ! A total of 12 telescope corrections were flagged in sub-array 6. ! Telescope amplitude corrections in sub-array 6: ! BR 1.03 KK 1.00* KP 0.95 LA 1.12 ! MK 1.01 OV 1.04 PT 1.13 ! ! ! Fit before self-cal, rms=0.107462Jy sigma=1.374708 ! Fit after self-cal, rms=0.076041Jy sigma=0.953041 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.259 mas, bmaj=1.831 mas, bpa=40.82 degrees ! Estimated noise=1.23483 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.000684083 Jy ! Component: 100 - total flux cleaned = -0.000851137 Jy ! Component: 150 - total flux cleaned = -0.000989867 Jy ! Component: 200 - total flux cleaned = -0.00105158 Jy ! Total flux subtracted in 200 components = -0.00105158 Jy ! Clean residual min=-0.004712 max=0.005592 Jy/beam ! Clean residual mean=0.000003 rms=0.001226 Jy/beam ! Combined flux in latest and established models = 1.23074 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 36 components and 1.23074 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.075818Jy sigma=0.950990 ! Fit after self-cal, rms=0.075790Jy sigma=0.950383 ! 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.075790Jy sigma=0.950383 ! Fit after self-cal, rms=0.075806Jy sigma=0.950240 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.259 mas, bmaj=1.832 mas, bpa=40.84 degrees ! Estimated noise=1.23522 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000437941 Jy ! Component: 100 - total flux cleaned = -0.000579301 Jy ! Component: 150 - total flux cleaned = -0.000640527 Jy ! Component: 200 - total flux cleaned = -0.000748845 Jy ! Total flux subtracted in 200 components = -0.000748845 Jy ! Clean residual min=-0.004582 max=0.005469 Jy/beam ! Clean residual mean=0.000003 rms=0.001214 Jy/beam ! Combined flux in latest and established models = 1.22999 Jy selfcal ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 36 components and 1.22999 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.075658Jy sigma=0.948899 ! Fit after self-cal, rms=0.075613Jy sigma=0.948627 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 (-32.6, 6.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00329508 Jy ! Component: 100 - total flux cleaned = 0.00385848 Jy ! Component: 150 - total flux cleaned = 0.00428751 Jy ! Component: 200 - total flux cleaned = 0.00421778 Jy ! Total flux subtracted in 200 components = 0.00421778 Jy ! Clean residual min=-0.004311 max=0.004274 Jy/beam ! Clean residual mean=-0.000010 rms=0.001148 Jy/beam ! Combined flux in latest and established models = 1.23421 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 54 components and 1.23421 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.075213Jy sigma=0.944202 ! Fit after self-cal, rms=0.075206Jy sigma=0.944073 ! 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.075206Jy sigma=0.944073 ! Fit after self-cal, rms=0.075197Jy sigma=0.943790 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.259 mas, bmaj=1.832 mas, bpa=40.86 degrees ! Estimated noise=1.23537 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000214472 Jy ! Component: 100 - total flux cleaned = 0.000151674 Jy ! Component: 150 - total flux cleaned = 0.000151753 Jy ! Component: 200 - total flux cleaned = 9.7656e-05 Jy ! Total flux subtracted in 200 components = 9.7656e-05 Jy ! Clean residual min=-0.004302 max=0.004268 Jy/beam ! Clean residual mean=-0.000009 rms=0.001133 Jy/beam ! Combined flux in latest and established models = 1.23431 Jy selfcal ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 60 components and 1.23431 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.075078Jy sigma=0.942689 ! Fit after self-cal, rms=0.075066Jy sigma=0.942618 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 = -5.6356e-05 Jy ! Component: 100 - total flux cleaned = -5.8674e-05 Jy ! Component: 150 - total flux cleaned = -0.000107096 Jy ! Component: 200 - total flux cleaned = -5.95763e-05 Jy ! Total flux subtracted in 200 components = -5.95763e-05 Jy ! Clean residual min=-0.004299 max=0.004273 Jy/beam ! Clean residual mean=-0.000010 rms=0.001127 Jy/beam ! Combined flux in latest and established models = 1.23425 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 66 components and 1.23425 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.074990Jy sigma=0.941928 ! Fit after self-cal, rms=0.074979Jy sigma=0.941873 ! 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.074979Jy sigma=0.941873 ! Fit after self-cal, rms=0.074978Jy sigma=0.941813 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.651 mas, bmaj=3.075 mas, bpa=52.96 degrees ! Estimated noise=1.40453 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000427641 Jy ! Component: 100 - total flux cleaned = -0.000390882 Jy ! Component: 150 - total flux cleaned = -0.00050228 Jy ! Component: 200 - total flux cleaned = -0.000525183 Jy ! Total flux subtracted in 200 components = -0.000525183 Jy ! Clean residual min=-0.005478 max=0.004388 Jy/beam ! Clean residual mean=-0.000017 rms=0.001391 Jy/beam ! Combined flux in latest and established models = 1.23372 Jy selfcal ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 76 components and 1.23373 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.074968Jy sigma=0.941635 ! Fit after self-cal, rms=0.074931Jy sigma=0.941374 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 = -9.213e-08 Jy ! Component: 100 - total flux cleaned = -2.5658e-05 Jy ! Component: 150 - total flux cleaned = -5.00272e-05 Jy ! Component: 200 - total flux cleaned = -5.98776e-06 Jy ! Total flux subtracted in 200 components = -5.98776e-06 Jy ! Clean residual min=-0.005520 max=0.004352 Jy/beam ! Clean residual mean=-0.000017 rms=0.001387 Jy/beam ! Combined flux in latest and established models = 1.23372 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 79 components and 1.23372 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.074933Jy sigma=0.941338 ! Fit after self-cal, rms=0.074908Jy sigma=0.941163 ! 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 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:NL,true ! 2:NL 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 selfant 3:BR,true ! 3:BR self-cal status=fixed weight=1 selfant 3:OV,true ! 3:OV self-cal status=fixed weight=1 selfant 4:BR,true ! 4:BR self-cal status=fixed weight=1 selfant 4:FD,true ! 4:FD self-cal status=fixed weight=1 selfant 4:HN,true ! 4:HN self-cal status=fixed weight=1 selfant 4:KP,true ! 4:KP self-cal status=fixed weight=1 selfant 4:LA,true ! 4:LA self-cal status=fixed weight=1 selfant 4:MK,true ! 4:MK 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:HN,true ! 5:HN self-cal status=fixed weight=1 selfant 5:KP,true ! 5:KP 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:PT,true ! 5:PT self-cal status=fixed weight=1 selfant 6:BR,true ! 6:BR self-cal status=fixed weight=1 selfant 6:KP,true ! 6:KP self-cal status=fixed weight=1 selfant 6:LA,true ! 6:LA self-cal status=fixed weight=1 selfant 6:MK,true ! 6:MK self-cal status=fixed weight=1 selfant 6:OV,true ! 6:OV self-cal status=fixed weight=1 selfant 6:PT,true ! 6:PT self-cal status=fixed weight=1 ![Exited script file: /tmp/ff.ff] selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 76 telescope corrections were flagged in sub-array 1. ! A total of 35 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 76 telescope corrections were flagged in sub-array 1. ! A total of 35 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 76 telescope corrections were flagged in sub-array 1. ! A total of 35 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 92 telescope corrections were flagged in sub-array 1. ! A total of 41 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.283859Jy sigma=2.332488 ! Fit after self-cal, rms=0.205900Jy sigma=2.204711 selfant "",false ! 1:HH self-cal status=correctable weight=1 ! 1:NY self-cal status=correctable weight=1 ! 1:WZ self-cal status=correctable weight=1 ! 1:ZC self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:HH self-cal status=correctable weight=1 ! 2:HN self-cal status=correctable weight=1 ! 2:KP self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:NY self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 2:ZC self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:NY self-cal status=correctable weight=1 ! 3:OV self-cal status=correctable weight=1 ! 4:BR self-cal status=correctable weight=1 ! 4:FD self-cal status=correctable weight=1 ! 4:HN self-cal status=correctable weight=1 ! 4:KK self-cal status=correctable weight=1 ! 4:KP self-cal status=correctable weight=1 ! 4:LA self-cal status=correctable weight=1 ! 4:MK 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:HN self-cal status=correctable weight=1 ! 5:KK self-cal status=correctable weight=1 ! 5:KP 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:PT self-cal status=correctable weight=1 ! 6:BR self-cal status=correctable weight=1 ! 6:KK self-cal status=correctable weight=1 ! 6:KP self-cal status=correctable weight=1 ! 6:LA self-cal status=correctable weight=1 ! 6:MK self-cal status=correctable weight=1 ! 6:OV self-cal status=correctable weight=1 ! 6:PT self-cal status=correctable weight=1 selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 62 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 13 telescope corrections were flagged in sub-array 4. ! A total of 29 telescope corrections were flagged in sub-array 5. ! A total of 28 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 62 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 13 telescope corrections were flagged in sub-array 4. ! A total of 29 telescope corrections were flagged in sub-array 5. ! A total of 28 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 62 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 13 telescope corrections were flagged in sub-array 4. ! A total of 29 telescope corrections were flagged in sub-array 5. ! A total of 28 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 62 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 13 telescope corrections were flagged in sub-array 4. ! A total of 29 telescope corrections were flagged in sub-array 5. ! A total of 28 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.205900Jy sigma=2.204711 ! Fit after self-cal, rms=0.143989Jy sigma=0.952315 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 ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Estimated beam: bmin=0.7245 mas, bmaj=1.074 mas, bpa=15.34 degrees ! Estimated noise=1.84385 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (-3, -4.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00466603 Jy ! Component: 100 - total flux cleaned = 0.00437679 Jy ! Component: 150 - total flux cleaned = 0.00356262 Jy ! Component: 200 - total flux cleaned = 0.00280395 Jy ! Total flux subtracted in 200 components = 0.00280395 Jy ! Clean residual min=-0.007735 max=0.008086 Jy/beam ! Clean residual mean=-0.000003 rms=0.001790 Jy/beam ! Combined flux in latest and established models = 1.23652 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 91 components and 1.23652 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.142324Jy sigma=0.947016 ! Fit after self-cal, rms=0.142017Jy sigma=0.945458 ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Added new window around map position (2.8, 36). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00223928 Jy ! Component: 100 - total flux cleaned = 0.00227094 Jy ! Component: 150 - total flux cleaned = 0.00205294 Jy ! Component: 200 - total flux cleaned = 0.00162079 Jy ! Total flux subtracted in 200 components = 0.00162079 Jy ! Clean residual min=-0.007322 max=0.007608 Jy/beam ! Clean residual mean=0.000001 rms=0.001717 Jy/beam ! Combined flux in latest and established models = 1.23814 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 106 components and 1.23814 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.141234Jy sigma=0.942497 ! Fit after self-cal, rms=0.141119Jy sigma=0.942133 ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Added new window around map position (-14.2, -8.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00480999 Jy ! Component: 100 - total flux cleaned = 0.00525196 Jy ! Component: 150 - total flux cleaned = 0.00546244 Jy ! Component: 200 - total flux cleaned = 0.00507356 Jy ! Total flux subtracted in 200 components = 0.00507356 Jy ! Clean residual min=-0.007101 max=0.007438 Jy/beam ! Clean residual mean=0.000001 rms=0.001658 Jy/beam ! Combined flux in latest and established models = 1.24322 Jy ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 120 components and 1.24322 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.140468Jy sigma=0.939318 ! Fit after self-cal, rms=0.140349Jy sigma=0.938983 ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Added new window around map position (46.4, 35.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00475785 Jy ! Component: 100 - total flux cleaned = 0.00593703 Jy ! Component: 150 - total flux cleaned = 0.00631376 Jy ! Component: 200 - total flux cleaned = 0.00649018 Jy ! Total flux subtracted in 200 components = 0.00649018 Jy ! Clean residual min=-0.006896 max=0.006818 Jy/beam ! Clean residual mean=-0.000001 rms=0.001602 Jy/beam ! Combined flux in latest and established models = 1.24971 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 131 components and 1.24971 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.139775Jy sigma=0.937209 ! Fit after self-cal, rms=0.139681Jy sigma=0.936885 ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Added new window around map position (-28.6, -27.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00373891 Jy ! Component: 100 - total flux cleaned = 0.00447849 Jy ! Component: 150 - total flux cleaned = 0.00465071 Jy ! Component: 200 - total flux cleaned = 0.0048117 Jy ! Total flux subtracted in 200 components = 0.0048117 Jy ! Clean residual min=-0.006743 max=0.006768 Jy/beam ! Clean residual mean=-0.000005 rms=0.001556 Jy/beam ! Combined flux in latest and established models = 1.25452 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 141 components and 1.25452 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.139232Jy sigma=0.935031 ! Fit after self-cal, rms=0.139176Jy sigma=0.934761 ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Added new window around map position (29.6, -0.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00232036 Jy ! Component: 100 - total flux cleaned = 0.00300703 Jy ! Component: 150 - total flux cleaned = 0.00300785 Jy ! Component: 200 - total flux cleaned = 0.00332245 Jy ! Total flux subtracted in 200 components = 0.00332245 Jy ! Clean residual min=-0.006485 max=0.006630 Jy/beam ! Clean residual mean=-0.000006 rms=0.001508 Jy/beam ! Combined flux in latest and established models = 1.25784 Jy ! Performing phase self-cal ! Adding 41 model components to the UV plane model. ! The established model now contains 149 components and 1.25784 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.138716Jy sigma=0.933117 ! Fit after self-cal, rms=0.138639Jy sigma=0.932872 ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Added new window around map position (18.6, 4.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00318448 Jy ! Component: 100 - total flux cleaned = 0.00368574 Jy ! Component: 150 - total flux cleaned = 0.00446431 Jy ! Component: 200 - total flux cleaned = 0.00416078 Jy ! Total flux subtracted in 200 components = 0.00416078 Jy ! Clean residual min=-0.006488 max=0.006186 Jy/beam ! Clean residual mean=-0.000007 rms=0.001465 Jy/beam ! Combined flux in latest and established models = 1.262 Jy ! Performing phase self-cal ! Adding 39 model components to the UV plane model. ! The established model now contains 162 components and 1.262 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.138206Jy sigma=0.930936 ! Fit after self-cal, rms=0.138202Jy sigma=0.930747 ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec 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 ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Estimated beam: bmin=0.9032 mas, bmaj=1.358 mas, bpa=18.35 degrees ! Estimated noise=1.24096 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.6, 2.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00105504 Jy ! Component: 100 - total flux cleaned = 0.000906182 Jy ! Component: 150 - total flux cleaned = 0.00089941 Jy ! Component: 200 - total flux cleaned = 0.00101522 Jy ! Component: 250 - total flux cleaned = 0.00101734 Jy ! Component: 300 - total flux cleaned = 0.00101599 Jy ! Component: 350 - total flux cleaned = 0.00101643 Jy ! Component: 400 - total flux cleaned = 0.000815546 Jy ! Total flux subtracted in 400 components = 0.000815546 Jy ! Clean residual min=-0.004531 max=0.004155 Jy/beam ! Clean residual mean=-0.000015 rms=0.001017 Jy/beam ! Combined flux in latest and established models = 1.26282 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 199 components and 1.26282 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.137496Jy sigma=0.926531 ! Fit after self-cal, rms=0.137502Jy sigma=0.926198 ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec 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.000555479 Jy ! Component: 100 - total flux cleaned = 0.000652402 Jy ! Component: 150 - total flux cleaned = 0.00046078 Jy ! Component: 200 - total flux cleaned = 0.000457868 Jy ! Component: 250 - total flux cleaned = 0.000185595 Jy ! Component: 300 - total flux cleaned = 0.000275963 Jy ! Component: 350 - total flux cleaned = 0.000188152 Jy ! Component: 400 - total flux cleaned = 0.000189177 Jy ! Total flux subtracted in 400 components = 0.000189177 Jy ! Clean residual min=-0.004418 max=0.004002 Jy/beam ! Clean residual mean=-0.000012 rms=0.000989 Jy/beam ! Combined flux in latest and established models = 1.26301 Jy ! Performing phase self-cal ! Adding 64 model components to the UV plane model. ! The established model now contains 216 components and 1.26301 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.137092Jy sigma=0.924187 ! Fit after self-cal, rms=0.137049Jy sigma=0.924045 ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec 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. 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 76 telescope corrections were flagged in sub-array 1. ! A total of 97 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 16 telescope corrections were flagged in sub-array 4. ! A total of 33 telescope corrections were flagged in sub-array 5. ! A total of 32 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 76 telescope corrections were flagged in sub-array 1. ! A total of 97 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 16 telescope corrections were flagged in sub-array 4. ! A total of 33 telescope corrections were flagged in sub-array 5. ! A total of 32 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 76 telescope corrections were flagged in sub-array 1. ! A total of 97 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 16 telescope corrections were flagged in sub-array 4. ! A total of 33 telescope corrections were flagged in sub-array 5. ! A total of 32 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 92 telescope corrections were flagged in sub-array 1. ! A total of 103 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 16 telescope corrections were flagged in sub-array 4. ! A total of 33 telescope corrections were flagged in sub-array 5. ! A total of 32 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.274164Jy sigma=2.173521 ! Fit after self-cal, rms=0.135021Jy sigma=0.921963 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Estimated beam: bmin=0.9002 mas, bmaj=1.357 mas, bpa=18.36 degrees ! Estimated noise=1.24052 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00113876 Jy ! Component: 100 - total flux cleaned = 0.00113853 Jy ! Component: 150 - total flux cleaned = 0.00122312 Jy ! Component: 200 - total flux cleaned = 0.00131161 Jy ! Component: 250 - total flux cleaned = 0.00139583 Jy ! Component: 300 - total flux cleaned = 0.00156015 Jy ! Component: 350 - total flux cleaned = 0.00156057 Jy ! Component: 400 - total flux cleaned = 0.0018007 Jy ! Total flux subtracted in 400 components = 0.0018007 Jy ! Clean residual min=-0.004088 max=0.003910 Jy/beam ! Clean residual mean=-0.000007 rms=0.000948 Jy/beam ! Combined flux in latest and established models = 1.26481 Jy selfcal ! Performing phase self-cal ! Adding 56 model components to the UV plane model. ! The established model now contains 229 components and 1.26481 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.134552Jy sigma=0.920137 ! Fit after self-cal, rms=0.134432Jy sigma=0.919981 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec 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 = 1.92893e-05 Jy ! Component: 100 - total flux cleaned = 9.98555e-05 Jy ! Component: 150 - total flux cleaned = 0.000100753 Jy ! Component: 200 - total flux cleaned = 0.000252794 Jy ! Component: 250 - total flux cleaned = 0.000252542 Jy ! Component: 300 - total flux cleaned = 0.000253865 Jy ! Component: 350 - total flux cleaned = 0.000398974 Jy ! Component: 400 - total flux cleaned = 0.000469394 Jy ! Total flux subtracted in 400 components = 0.000469394 Jy ! Clean residual min=-0.004052 max=0.003933 Jy/beam ! Clean residual mean=-0.000007 rms=0.000933 Jy/beam ! Combined flux in latest and established models = 1.26528 Jy ! Performing phase self-cal ! Adding 65 model components to the UV plane model. ! The established model now contains 239 components and 1.26528 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.134113Jy sigma=0.918797 ! Fit after self-cal, rms=0.134033Jy sigma=0.918702 ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec 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.134033Jy sigma=0.918702 ! Fit after self-cal, rms=0.133622Jy sigma=0.918459 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Estimated beam: bmin=0.8993 mas, bmaj=1.357 mas, bpa=18.33 degrees ! Estimated noise=1.24084 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00087481 Jy ! Component: 100 - total flux cleaned = 0.000799798 Jy ! Component: 150 - total flux cleaned = 0.000800418 Jy ! Component: 200 - total flux cleaned = 0.000869159 Jy ! Component: 250 - total flux cleaned = 0.000939776 Jy ! Component: 300 - total flux cleaned = 0.00093795 Jy ! Component: 350 - total flux cleaned = 0.000803594 Jy ! Component: 400 - total flux cleaned = 0.000869803 Jy ! Total flux subtracted in 400 components = 0.000869803 Jy ! Clean residual min=-0.003972 max=0.003959 Jy/beam ! Clean residual mean=-0.000006 rms=0.000920 Jy/beam ! Combined flux in latest and established models = 1.26615 Jy selfcal ! Performing phase self-cal ! Adding 69 model components to the UV plane model. ! The established model now contains 248 components and 1.26615 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.133306Jy sigma=0.917315 ! Fit after self-cal, rms=0.133232Jy sigma=0.917227 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec 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 (46.6, -42.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00286763 Jy ! Component: 100 - total flux cleaned = 0.00305454 Jy ! Component: 150 - total flux cleaned = 0.00297264 Jy ! Component: 200 - total flux cleaned = 0.00297077 Jy ! Component: 250 - total flux cleaned = 0.00297287 Jy ! Component: 300 - total flux cleaned = 0.00310965 Jy ! Component: 350 - total flux cleaned = 0.00324223 Jy ! Component: 400 - total flux cleaned = 0.00337279 Jy ! Total flux subtracted in 400 components = 0.00337279 Jy ! Clean residual min=-0.003776 max=0.003790 Jy/beam ! Clean residual mean=-0.000013 rms=0.000900 Jy/beam ! Combined flux in latest and established models = 1.26952 Jy ! Performing phase self-cal ! Adding 70 model components to the UV plane model. ! The established model now contains 277 components and 1.26952 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.132915Jy sigma=0.915120 ! Fit after self-cal, rms=0.132896Jy sigma=0.914993 ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Added new window around map position (-15.2, -27.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00181797 Jy ! Component: 100 - total flux cleaned = 0.00201951 Jy ! Component: 150 - total flux cleaned = 0.00244408 Jy ! Component: 200 - total flux cleaned = 0.00276622 Jy ! Component: 250 - total flux cleaned = 0.00299604 Jy ! Component: 300 - total flux cleaned = 0.00336015 Jy ! Component: 350 - total flux cleaned = 0.00363887 Jy ! Component: 400 - total flux cleaned = 0.00370836 Jy ! Total flux subtracted in 400 components = 0.00370836 Jy ! Clean residual min=-0.003610 max=0.003826 Jy/beam ! Clean residual mean=-0.000015 rms=0.000857 Jy/beam ! Combined flux in latest and established models = 1.27323 Jy ! Performing phase self-cal ! Adding 57 model components to the UV plane model. ! The established model now contains 303 components and 1.27323 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.132538Jy sigma=0.912308 ! Fit after self-cal, rms=0.132517Jy sigma=0.912244 ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Added new window around map position (17, 28). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00205676 Jy ! Component: 100 - total flux cleaned = 0.00240466 Jy ! Component: 150 - total flux cleaned = 0.002796 Jy ! Component: 200 - total flux cleaned = 0.00360822 Jy ! Component: 250 - total flux cleaned = 0.00417359 Jy ! Component: 300 - total flux cleaned = 0.00479088 Jy ! Component: 350 - total flux cleaned = 0.00532668 Jy ! Component: 400 - total flux cleaned = 0.00565462 Jy ! Total flux subtracted in 400 components = 0.00565462 Jy ! Clean residual min=-0.003463 max=0.003588 Jy/beam ! Clean residual mean=-0.000017 rms=0.000829 Jy/beam ! Combined flux in latest and established models = 1.27888 Jy ! Performing phase self-cal ! Adding 62 model components to the UV plane model. ! The established model now contains 329 components and 1.27888 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.132226Jy sigma=0.910114 ! Fit after self-cal, rms=0.132203Jy sigma=0.910037 ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Added new window around map position (44.4, -8.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00169783 Jy ! Component: 100 - total flux cleaned = 0.00126488 Jy ! Component: 150 - total flux cleaned = 0.00126128 Jy ! Component: 200 - total flux cleaned = 0.00140107 Jy ! Component: 250 - total flux cleaned = 0.00146949 Jy ! Component: 300 - total flux cleaned = 0.00147073 Jy ! Component: 350 - total flux cleaned = 0.00173549 Jy ! Component: 400 - total flux cleaned = 0.00179982 Jy ! Total flux subtracted in 400 components = 0.00179982 Jy ! Clean residual min=-0.003430 max=0.003349 Jy/beam ! Clean residual mean=-0.000021 rms=0.000802 Jy/beam ! Combined flux in latest and established models = 1.28068 Jy ! Performing phase self-cal ! Adding 68 model components to the UV plane model. ! The established model now contains 353 components and 1.28068 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.131931Jy sigma=0.907937 ! Fit after self-cal, rms=0.131922Jy sigma=0.907844 ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec 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.131922Jy sigma=0.907844 ! Fit after self-cal, rms=0.129476Jy sigma=0.887653 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Estimated beam: bmin=2.595 mas, bmaj=3.015 mas, bpa=50.41 degrees ! Estimated noise=1.39051 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000413553 Jy ! Component: 100 - total flux cleaned = 0.000938317 Jy ! Component: 150 - total flux cleaned = 0.00140775 Jy ! Component: 200 - total flux cleaned = 0.00184094 Jy ! Component: 250 - total flux cleaned = 0.00220481 Jy ! Component: 300 - total flux cleaned = 0.00250859 Jy ! Component: 350 - total flux cleaned = 0.00284339 Jy ! Component: 400 - total flux cleaned = 0.00308566 Jy ! Total flux subtracted in 400 components = 0.00308566 Jy ! Clean residual min=-0.003830 max=0.002915 Jy/beam ! Clean residual mean=-0.000040 rms=0.000879 Jy/beam ! Combined flux in latest and established models = 1.28377 Jy selfcal ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 385 components and 1.28377 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.129472Jy sigma=0.886883 ! Fit after self-cal, rms=0.129483Jy sigma=0.886823 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec 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.000303651 Jy ! Component: 100 - total flux cleaned = 0.000710453 Jy ! Component: 150 - total flux cleaned = 0.000826851 Jy ! Component: 200 - total flux cleaned = 0.00108855 Jy ! Component: 250 - total flux cleaned = 0.00116201 Jy ! Component: 300 - total flux cleaned = 0.00126985 Jy ! Component: 350 - total flux cleaned = 0.00144305 Jy ! Component: 400 - total flux cleaned = 0.00151123 Jy ! Total flux subtracted in 400 components = 0.00151123 Jy ! Clean residual min=-0.003765 max=0.002936 Jy/beam ! Clean residual mean=-0.000040 rms=0.000850 Jy/beam ! Combined flux in latest and established models = 1.28528 Jy ! Performing phase self-cal ! Adding 52 model components to the UV plane model. ! The established model now contains 406 components and 1.28528 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.129494Jy sigma=0.886450 ! Fit after self-cal, rms=0.129498Jy sigma=0.886430 ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec 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 ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Estimated beam: bmin=0.899 mas, bmaj=1.358 mas, bpa=18.33 degrees ! Estimated noise=1.24198 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.899 x 1.358 at 18.33 degrees (North through East) ! Clean map min=-0.0046689 max=1.0929 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.129498Jy sigma=0.886430 ! Fit after self-cal, rms=0.113214Jy sigma=0.770501 ! Inverting map and beam ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Estimated beam: bmin=0.8989 mas, bmaj=1.358 mas, bpa=18.31 degrees ! Estimated noise=1.24216 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000154208 Jy ! Component: 100 - total flux cleaned = 0.000637886 Jy ! Component: 150 - total flux cleaned = 0.00101715 Jy ! Component: 200 - total flux cleaned = 0.00145271 Jy ! Component: 250 - total flux cleaned = 0.00180116 Jy ! Component: 300 - total flux cleaned = 0.00186703 Jy ! Component: 350 - total flux cleaned = 0.00193157 Jy ! Component: 400 - total flux cleaned = 0.00180286 Jy ! Total flux subtracted in 400 components = 0.00180286 Jy ! Clean residual min=-0.003263 max=0.002974 Jy/beam ! Clean residual mean=-0.000025 rms=0.000738 Jy/beam ! Combined flux in latest and established models = 1.28708 Jy ! Performing phase self-cal ! Adding 60 model components to the UV plane model. ! The established model now contains 430 components and 1.28708 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.112847Jy sigma=0.768896 ! Fit after self-cal, rms=0.112812Jy sigma=0.768842 ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000265393 Jy ! Component: 100 - total flux cleaned = 0.000457035 Jy ! Component: 150 - total flux cleaned = 0.000456829 Jy ! Component: 200 - total flux cleaned = 0.000335278 Jy ! Component: 250 - total flux cleaned = 0.000517457 Jy ! Component: 300 - total flux cleaned = 0.000398121 Jy ! Component: 350 - total flux cleaned = 0.000573084 Jy ! Component: 400 - total flux cleaned = 0.000455567 Jy ! Total flux subtracted in 400 components = 0.000455567 Jy ! Clean residual min=-0.003263 max=0.002965 Jy/beam ! Clean residual mean=-0.000026 rms=0.000729 Jy/beam ! Combined flux in latest and established models = 1.28754 Jy ! Performing phase self-cal ! Adding 95 model components to the UV plane model. ! The established model now contains 455 components and 1.28754 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.112559Jy sigma=0.767790 ! Fit after self-cal, rms=0.112538Jy sigma=0.767767 ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.112538Jy sigma=0.767767 ! Fit after self-cal, rms=0.112473Jy sigma=0.767655 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 ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Estimated beam: bmin=0.7226 mas, bmaj=1.076 mas, bpa=15.32 degrees ! Estimated noise=1.84817 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.836842 Jy ! Component: 100 - total flux cleaned = 1.04003 Jy ! Component: 150 - total flux cleaned = 1.12445 Jy ! Component: 200 - total flux cleaned = 1.16982 Jy ! Component: 250 - total flux cleaned = 1.19614 Jy ! Component: 300 - total flux cleaned = 1.21216 Jy ! Component: 350 - total flux cleaned = 1.22282 Jy ! Component: 400 - total flux cleaned = 1.23081 Jy ! Total flux subtracted in 400 components = 1.23081 Jy ! Clean residual min=-0.007354 max=0.007135 Jy/beam ! Clean residual mean=-0.000011 rms=0.001629 Jy/beam ! Combined flux in latest and established models = 1.23081 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 39 model components to the UV plane model. ! The established model now contains 39 components and 1.23081 Jy ! Inverting map and beam ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Estimated beam: bmin=0.8986 mas, bmaj=1.359 mas, bpa=18.32 degrees ! Estimated noise=1.2425 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.6, -1.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00873439 Jy ! Component: 100 - total flux cleaned = 0.0148424 Jy ! Component: 150 - total flux cleaned = 0.0193633 Jy ! Component: 200 - total flux cleaned = 0.0231558 Jy ! Component: 250 - total flux cleaned = 0.0256004 Jy ! Component: 300 - total flux cleaned = 0.0280927 Jy ! Component: 350 - total flux cleaned = 0.0299621 Jy ! Component: 400 - total flux cleaned = 0.0322388 Jy ! Total flux subtracted in 400 components = 0.0322388 Jy ! Clean residual min=-0.004242 max=0.004372 Jy/beam ! Clean residual mean=0.000001 rms=0.001046 Jy/beam ! Combined flux in latest and established models = 1.26304 Jy ! Adding 62 model components to the UV plane model. ! The established model now contains 89 components and 1.26304 Jy ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Added new window around map position (-14.6, -33.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00263315 Jy ! Component: 100 - total flux cleaned = 0.00446567 Jy ! Component: 150 - total flux cleaned = 0.00607687 Jy ! Component: 200 - total flux cleaned = 0.00721403 Jy ! Component: 250 - total flux cleaned = 0.00844967 Jy ! Component: 300 - total flux cleaned = 0.00937758 Jy ! Component: 350 - total flux cleaned = 0.0102818 Jy ! Component: 400 - total flux cleaned = 0.011034 Jy ! Total flux subtracted in 400 components = 0.011034 Jy ! Clean residual min=-0.003769 max=0.003963 Jy/beam ! Clean residual mean=-0.000011 rms=0.000945 Jy/beam ! Combined flux in latest and established models = 1.27408 Jy ! Adding 81 model components to the UV plane model. ! The established model now contains 125 components and 1.27408 Jy ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Added new window around map position (-1.8, -16.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00176467 Jy ! Component: 100 - total flux cleaned = 0.00340328 Jy ! Component: 150 - total flux cleaned = 0.00412992 Jy ! Component: 200 - total flux cleaned = 0.00495131 Jy ! Component: 250 - total flux cleaned = 0.00540692 Jy ! Component: 300 - total flux cleaned = 0.00585344 Jy ! Component: 350 - total flux cleaned = 0.00607124 Jy ! Component: 400 - total flux cleaned = 0.00671895 Jy ! Total flux subtracted in 400 components = 0.00671895 Jy ! Clean residual min=-0.003566 max=0.003658 Jy/beam ! Clean residual mean=-0.000007 rms=0.000874 Jy/beam ! Combined flux in latest and established models = 1.2808 Jy ! Adding 85 model components to the UV plane model. ! The established model now contains 155 components and 1.2808 Jy ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Added new window around map position (-18.4, 4.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00243736 Jy ! Component: 100 - total flux cleaned = 0.00397541 Jy ! Component: 150 - total flux cleaned = 0.00513365 Jy ! Component: 200 - total flux cleaned = 0.00595252 Jy ! Component: 250 - total flux cleaned = 0.00645453 Jy ! Component: 300 - total flux cleaned = 0.00704242 Jy ! Component: 350 - total flux cleaned = 0.00733132 Jy ! Component: 400 - total flux cleaned = 0.00752401 Jy ! Total flux subtracted in 400 components = 0.00752401 Jy ! Clean residual min=-0.003520 max=0.003220 Jy/beam ! Clean residual mean=-0.000008 rms=0.000817 Jy/beam ! Combined flux in latest and established models = 1.28832 Jy ! Adding 99 model components to the UV plane model. ! The established model now contains 191 components and 1.28832 Jy ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec 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 ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Estimated beam: bmin=2.595 mas, bmaj=3.016 mas, bpa=50.49 degrees ! Estimated noise=1.39074 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000956103 Jy ! Component: 100 - total flux cleaned = 0.00232535 Jy ! Component: 150 - total flux cleaned = 0.00331812 Jy ! Component: 200 - total flux cleaned = 0.00402608 Jy ! Component: 250 - total flux cleaned = 0.0046224 Jy ! Component: 300 - total flux cleaned = 0.00518779 Jy ! Component: 350 - total flux cleaned = 0.00578654 Jy ! Component: 400 - total flux cleaned = 0.00641699 Jy ! Total flux subtracted in 400 components = 0.00641699 Jy ! Clean residual min=-0.003437 max=0.002899 Jy/beam ! Clean residual mean=-0.000013 rms=0.000827 Jy/beam ! Combined flux in latest and established models = 1.29474 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.000608618 Jy ! Component: 100 - total flux cleaned = 0.00103532 Jy ! Component: 150 - total flux cleaned = 0.00150286 Jy ! Component: 200 - total flux cleaned = 0.00200584 Jy ! Component: 250 - total flux cleaned = 0.00239665 Jy ! Component: 300 - total flux cleaned = 0.00258603 Jy ! Component: 350 - total flux cleaned = 0.00295879 Jy ! Component: 400 - total flux cleaned = 0.00323038 Jy ! Total flux subtracted in 400 components = 0.00323038 Jy ! Clean residual min=-0.003396 max=0.002790 Jy/beam ! Clean residual mean=-0.000012 rms=0.000786 Jy/beam ! Combined flux in latest and established models = 1.29797 Jy ! Adding 79 model components to the UV plane model. ! The established model now contains 260 components and 1.29797 Jy ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! Estimated beam: bmin=0.8986 mas, bmaj=1.359 mas, bpa=18.32 degrees ! Estimated noise=1.2425 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000525807 Jy ! Component: 100 - total flux cleaned = -0.000333729 Jy ! Component: 150 - total flux cleaned = -0.000333038 Jy ! Component: 200 - total flux cleaned = -0.000334901 Jy ! Component: 250 - total flux cleaned = -0.000248777 Jy ! Component: 300 - total flux cleaned = 6.09847e-06 Jy ! Component: 350 - total flux cleaned = -7.93039e-05 Jy ! Component: 400 - total flux cleaned = 8.33048e-05 Jy ! Total flux subtracted in 400 components = 8.33048e-05 Jy ! Clean residual min=-0.003359 max=0.002945 Jy/beam ! Clean residual mean=-0.000008 rms=0.000735 Jy/beam ! Combined flux in latest and established models = 1.29805 Jy selfcal ! Performing phase self-cal ! Adding 75 model components to the UV plane model. ! The established model now contains 288 components and 1.29805 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.112275Jy sigma=0.770063 ! Fit after self-cal, rms=0.111761Jy sigma=0.768765 wmodel /scr/pima/rdv62_uvs/J1130+3815_X_map.mod ! Writing 288 model components to file: /scr/pima/rdv62_uvs/J1130+3815_X_map.mod wobs /scr/pima/rdv62_uvs/J1130+3815_X_uvs.fits ! Applying 1321 buffered edits. ! Writing UV FITS file: /scr/pima/rdv62_uvs/J1130+3815_X_uvs.fits wwins /scr/pima/rdv62_uvs/J1130+3815_X_map.win ! wwins: Wrote 18 windows to /scr/pima/rdv62_uvs/J1130+3815_X_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.00149953 Jy ! Component: 100 - total flux cleaned = -0.00227047 Jy ! Component: 150 - total flux cleaned = -0.0026921 Jy ! Component: 200 - total flux cleaned = -0.0030513 Jy ! Component: 250 - total flux cleaned = -0.0034031 Jy ! Component: 300 - total flux cleaned = -0.00379682 Jy ! Component: 350 - total flux cleaned = -0.00413565 Jy ! Component: 400 - total flux cleaned = -0.00456438 Jy ! Component: 450 - total flux cleaned = -0.00493965 Jy ! Component: 500 - total flux cleaned = -0.00531038 Jy ! Component: 550 - total flux cleaned = -0.00567593 Jy ! Component: 600 - total flux cleaned = -0.00590138 Jy ! Component: 650 - total flux cleaned = -0.00616885 Jy ! Component: 700 - total flux cleaned = -0.00638979 Jy ! Component: 750 - total flux cleaned = -0.00660868 Jy ! Component: 800 - total flux cleaned = -0.00686833 Jy ! Component: 850 - total flux cleaned = -0.0069115 Jy ! Component: 900 - total flux cleaned = -0.00716638 Jy ! Component: 950 - total flux cleaned = -0.00733475 Jy ! Component: 1000 - total flux cleaned = -0.00737637 Jy ! Component: 1050 - total flux cleaned = -0.00741741 Jy ! Component: 1100 - total flux cleaned = -0.00754104 Jy ! Component: 1150 - total flux cleaned = -0.00758186 Jy ! Component: 1200 - total flux cleaned = -0.00790569 Jy ! Component: 1250 - total flux cleaned = -0.00798573 Jy ! Component: 1300 - total flux cleaned = -0.00806531 Jy ! Component: 1350 - total flux cleaned = -0.00834235 Jy ! Component: 1400 - total flux cleaned = -0.00834234 Jy ! Component: 1450 - total flux cleaned = -0.00842064 Jy ! Component: 1500 - total flux cleaned = -0.00845935 Jy ! Component: 1550 - total flux cleaned = -0.00865187 Jy ! Component: 1600 - total flux cleaned = -0.00865207 Jy ! Component: 1650 - total flux cleaned = -0.00884191 Jy ! Component: 1700 - total flux cleaned = -0.0089172 Jy ! Component: 1750 - total flux cleaned = -0.00895438 Jy ! Component: 1800 - total flux cleaned = -0.00914081 Jy ! Component: 1850 - total flux cleaned = -0.00925218 Jy ! Component: 1900 - total flux cleaned = -0.00932603 Jy ! Component: 1950 - total flux cleaned = -0.00936269 Jy ! Component: 2000 - total flux cleaned = -0.00936238 Jy ! Component: 2050 - total flux cleaned = -0.00936252 Jy ! Component: 2100 - total flux cleaned = -0.0094706 Jy ! Component: 2150 - total flux cleaned = -0.00947052 Jy ! Component: 2200 - total flux cleaned = -0.00954183 Jy ! Component: 2250 - total flux cleaned = -0.0095772 Jy ! Component: 2300 - total flux cleaned = -0.0097182 Jy ! Component: 2350 - total flux cleaned = -0.00982365 Jy ! Component: 2400 - total flux cleaned = -0.00989374 Jy ! Component: 2450 - total flux cleaned = -0.00982398 Jy ! Component: 2500 - total flux cleaned = -0.00992772 Jy ! Component: 2550 - total flux cleaned = -0.0100997 Jy ! Component: 2600 - total flux cleaned = -0.0102025 Jy ! Component: 2650 - total flux cleaned = -0.0103049 Jy ! Component: 2700 - total flux cleaned = -0.0102027 Jy ! Component: 2750 - total flux cleaned = -0.0102704 Jy ! Component: 2800 - total flux cleaned = -0.0102368 Jy ! Component: 2850 - total flux cleaned = -0.0103038 Jy ! Component: 2900 - total flux cleaned = -0.0103369 Jy ! Component: 2950 - total flux cleaned = -0.0103703 Jy ! Component: 3000 - total flux cleaned = -0.0103372 Jy ! Component: 3050 - total flux cleaned = -0.0102055 Jy ! Component: 3100 - total flux cleaned = -0.0104024 Jy ! Component: 3150 - total flux cleaned = -0.0102717 Jy ! Component: 3200 - total flux cleaned = -0.0104347 Jy ! Component: 3250 - total flux cleaned = -0.0104348 Jy ! Component: 3300 - total flux cleaned = -0.0105316 Jy ! Component: 3350 - total flux cleaned = -0.0105316 Jy ! Component: 3400 - total flux cleaned = -0.0105637 Jy ! Component: 3450 - total flux cleaned = -0.010468 Jy ! Component: 3500 - total flux cleaned = -0.0106586 Jy ! Component: 3550 - total flux cleaned = -0.0106269 Jy ! Component: 3600 - total flux cleaned = -0.0106898 Jy ! Component: 3650 - total flux cleaned = -0.0105327 Jy ! Component: 3700 - total flux cleaned = -0.0105954 Jy ! Component: 3750 - total flux cleaned = -0.0107514 Jy ! Component: 3800 - total flux cleaned = -0.0108136 Jy ! Component: 3850 - total flux cleaned = -0.0106589 Jy ! Component: 3900 - total flux cleaned = -0.0107824 Jy ! Component: 3950 - total flux cleaned = -0.0107825 Jy ! Component: 4000 - total flux cleaned = -0.0106294 Jy ! Component: 4050 - total flux cleaned = -0.0106295 Jy ! Total flux subtracted in 4096 components = -0.0106601 Jy ! Clean residual min=-0.001848 max=0.001605 Jy/beam ! Clean residual mean=-0.000002 rms=0.000497 Jy/beam ! Combined flux in latest and established models = 1.28739 Jy keep ! Adding 313 model components to the UV plane model. ! The established model now contains 600 components and 1.28739 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv62_uvs/J1130+3815_X\_map.fits" ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.8986 x 1.359 at 18.32 degrees (North through East) ! Clean map min=-0.0034391 max=1.0931 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv62_uvs/J1130+3815_X_map.fits float peak_flux peak_flux = peak(flux,max) invert ! Inverting map ! Your choice of large map pixels excluded 7.13% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.191 milli-arcsec float image_peak image_peak = int(10000*peak_flux)/10000 float image_rms image_rms = int(100000*imstat(rms))/100000 print "=========================================" ! ========================================= print "The image peak= ", image_peak, " rms= ", image_rms, " Jy" ! The image peak= 1.093 rms= 0.00049 Jy quit ! Quitting program ! Log file /scr/pima/rdv62_uvs/J1130+3815_X_uvs.log closed on Fri Nov 22 22:37:28 2013