! Started logfile: /scr/pima/rdv62_uvs/J1419+5423_X_uvs.log on Fri Nov 22 22:39:12 2013 obs /scr/pima/rdv62_uvs/J1419+5423_X_uva.fits ! Reading UV FITS file: /scr/pima/rdv62_uvs/J1419+5423_X_uva.fits ! AN table 1: 6 integrations on 3 of 3 possible baselines. ! AN table 2: 14 integrations on 3 of 3 possible baselines. ! AN table 3: 4 integrations on 1 of 1 possible baselines. ! AN table 4: 8 integrations on 6 of 6 possible baselines. ! AN table 5: 4 integrations on 3 of 3 possible baselines. ! AN table 6: 5 integrations on 1 of 1 possible baselines. ! AN table 7: 6 integrations on 15 of 15 possible baselines. ! AN table 8: 6 integrations on 3 of 3 possible baselines. ! AN table 9: 10 integrations on 45 of 45 possible baselines. ! AN table 10: 26 integrations on 66 of 66 possible baselines. ! AN table 11: 13 integrations on 66 of 66 possible baselines. ! AN table 12: 9 integrations on 6 of 6 possible baselines. ! AN table 13: 11 integrations on 55 of 55 possible baselines. ! AN table 14: 8 integrations on 3 of 3 possible baselines. ! Apparent sampling: 0.485294 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1419+5423 ! ! 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 7656 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/J1419+5423_X] integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 800000 float soltime2 soltime2 = 15 float soltime3 soltime3 = 1 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float old_peak float new_peak float flux_cutoff selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:MK self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:ZC self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 4:BR self-cal status=correctable weight=1 ! 4:HN self-cal status=correctable weight=1 ! 4:NL self-cal status=correctable weight=1 ! 4:NY self-cal status=correctable weight=1 ! 5:BR self-cal status=correctable weight=1 ! 5:HN self-cal status=correctable weight=1 ! 5:NL self-cal status=correctable weight=1 ! 6:HN self-cal status=correctable weight=1 ! 6:SC self-cal status=correctable weight=1 ! 7:BR self-cal status=correctable weight=1 ! 7:HN self-cal status=correctable weight=1 ! 7:LA self-cal status=correctable weight=1 ! 7:NL self-cal status=correctable weight=1 ! 7:PT self-cal status=correctable weight=1 ! 7:SC self-cal status=correctable weight=1 ! 8:NY self-cal status=correctable weight=1 ! 8:WZ self-cal status=correctable weight=1 ! 8:ZC self-cal status=correctable weight=1 ! 9:BR self-cal status=correctable weight=1 ! 9:FD self-cal status=correctable weight=1 ! 9:HN self-cal status=correctable weight=1 ! 9:KP self-cal status=correctable weight=1 ! 9:LA self-cal status=correctable weight=1 ! 9:NL self-cal status=correctable weight=1 ! 9:OV self-cal status=correctable weight=1 ! 9:PT self-cal status=correctable weight=1 ! 9:SC self-cal status=correctable weight=1 ! 9:ZC self-cal status=correctable weight=1 ! 10:BR self-cal status=correctable weight=1 ! 10:FD self-cal status=correctable weight=1 ! 10:HN self-cal status=correctable weight=1 ! 10:KK self-cal status=correctable weight=1 ! 10:KP self-cal status=correctable weight=1 ! 10:LA self-cal status=correctable weight=1 ! 10:MK self-cal status=correctable weight=1 ! 10:NL self-cal status=correctable weight=1 ! 10:NY self-cal status=correctable weight=1 ! 10:OV self-cal status=correctable weight=1 ! 10:PT self-cal status=correctable weight=1 ! 10:ZC self-cal status=correctable weight=1 ! 11:BR self-cal status=correctable weight=1 ! 11:FD self-cal status=correctable weight=1 ! 11:HN self-cal status=correctable weight=1 ! 11:KK self-cal status=correctable weight=1 ! 11:KP self-cal status=correctable weight=1 ! 11:LA self-cal status=correctable weight=1 ! 11:MK self-cal status=correctable weight=1 ! 11:NL self-cal status=correctable weight=1 ! 11:NY self-cal status=correctable weight=1 ! 11:OV self-cal status=correctable weight=1 ! 11:PT self-cal status=correctable weight=1 ! 11:SC self-cal status=correctable weight=1 ! 12:BR self-cal status=correctable weight=1 ! 12:KK self-cal status=correctable weight=1 ! 12:MK self-cal status=correctable weight=1 ! 12:NY self-cal status=correctable weight=1 ! 13:BR self-cal status=correctable weight=1 ! 13:FD self-cal status=correctable weight=1 ! 13:HN self-cal status=correctable weight=1 ! 13:KK self-cal status=correctable weight=1 ! 13:KP self-cal status=correctable weight=1 ! 13:LA self-cal status=correctable weight=1 ! 13:MK self-cal status=correctable weight=1 ! 13:NL self-cal status=correctable weight=1 ! 13:NY self-cal status=correctable weight=1 ! 13:OV self-cal status=correctable weight=1 ! 13:PT self-cal status=correctable weight=1 ! 14:KK self-cal status=correctable weight=1 ! 14:MK self-cal status=correctable weight=1 ! 14:WZ 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 9 telescope corrections were flagged in sub-array 1. ! A total of 27 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 17 telescope corrections were flagged in sub-array 4. ! A total of 3 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 8 telescope corrections were flagged in sub-array 7. ! A total of 9 telescope corrections were flagged in sub-array 8. ! A total of 18 telescope corrections were flagged in sub-array 9. ! A total of 201 telescope corrections were flagged in sub-array 10. ! A total of 49 telescope corrections were flagged in sub-array 11. ! A total of 12 telescope corrections were flagged in sub-array 12. ! A total of 43 telescope corrections were flagged in sub-array 13. ! A total of 9 telescope corrections were flagged in sub-array 14. ! ! Correcting IF 2. ! A total of 9 telescope corrections were flagged in sub-array 1. ! A total of 27 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 17 telescope corrections were flagged in sub-array 4. ! A total of 3 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 8 telescope corrections were flagged in sub-array 7. ! A total of 9 telescope corrections were flagged in sub-array 8. ! A total of 18 telescope corrections were flagged in sub-array 9. ! A total of 201 telescope corrections were flagged in sub-array 10. ! A total of 49 telescope corrections were flagged in sub-array 11. ! A total of 12 telescope corrections were flagged in sub-array 12. ! A total of 43 telescope corrections were flagged in sub-array 13. ! A total of 9 telescope corrections were flagged in sub-array 14. ! ! Correcting IF 3. ! A total of 9 telescope corrections were flagged in sub-array 1. ! A total of 27 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 17 telescope corrections were flagged in sub-array 4. ! A total of 3 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 8 telescope corrections were flagged in sub-array 7. ! A total of 9 telescope corrections were flagged in sub-array 8. ! A total of 18 telescope corrections were flagged in sub-array 9. ! A total of 201 telescope corrections were flagged in sub-array 10. ! A total of 49 telescope corrections were flagged in sub-array 11. ! A total of 12 telescope corrections were flagged in sub-array 12. ! A total of 43 telescope corrections were flagged in sub-array 13. ! A total of 9 telescope corrections were flagged in sub-array 14. ! ! Correcting IF 4. ! A total of 9 telescope corrections were flagged in sub-array 1. ! A total of 42 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 17 telescope corrections were flagged in sub-array 4. ! A total of 3 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 8 telescope corrections were flagged in sub-array 7. ! A total of 18 telescope corrections were flagged in sub-array 8. ! A total of 27 telescope corrections were flagged in sub-array 9. ! A total of 233 telescope corrections were flagged in sub-array 10. ! A total of 49 telescope corrections were flagged in sub-array 11. ! A total of 12 telescope corrections were flagged in sub-array 12. ! A total of 43 telescope corrections were flagged in sub-array 13. ! A total of 9 telescope corrections were flagged in sub-array 14. ! ! Fit before self-cal, rms=1.168224Jy sigma=8.441891 ! Fit after self-cal, rms=0.485433Jy sigma=2.469511 ! 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:ZC ! Flagging baselines of antenna 2:ZC in the currently selected channels. flag 4:NY ! Flagging baselines of antenna 4:NY in the currently selected channels. flag 8:NY ! Flagging baselines of antenna 8:NY in the currently selected channels. flag 8:WZ ! Flagging baselines of antenna 8:WZ in the currently selected channels. flag 8:ZC ! Flagging baselines of antenna 8:ZC in the currently selected channels. flag 9:ZC ! Flagging baselines of antenna 9:ZC in the currently selected channels. flag 10:KK ! Flagging baselines of antenna 10:KK in the currently selected channels. flag 10:NY ! Flagging baselines of antenna 10:NY in the currently selected channels. flag 10:ZC ! Flagging baselines of antenna 10:ZC in the currently selected channels. flag 11:KK ! Flagging baselines of antenna 11:KK in the currently selected channels. flag 11:NY ! Flagging baselines of antenna 11:NY in the currently selected channels. flag 12:KK ! Flagging baselines of antenna 12:KK in the currently selected channels. flag 12:NY ! Flagging baselines of antenna 12:NY in the currently selected channels. flag 13:KK ! Flagging baselines of antenna 13:KK in the currently selected channels. flag 13:NY ! Flagging baselines of antenna 13:NY in the currently selected channels. flag 14:KK ! Flagging baselines of antenna 14:KK in the currently selected channels. flag 14:WZ ! Flagging baselines of antenna 14:WZ in the currently selected channels. ![Exited script file: /tmp/ff.ff] uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.9956 mas, bmaj=1.299 mas, bpa=55.81 degrees ! Estimated noise=2.14602 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.485822 Jy ! Component: 100 - total flux cleaned = 0.617116 Jy ! Total flux subtracted in 100 components = 0.617116 Jy ! Clean residual min=-0.020331 max=0.057856 Jy/beam ! Clean residual mean=0.000139 rms=0.005239 Jy/beam ! Combined flux in latest and established models = 0.617116 Jy ! Performing phase self-cal ! Adding 6 model components to the UV plane model. ! The established model now contains 6 components and 0.617116 Jy ! ! Correcting IF 1. ! A total of 15 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 9 telescope corrections were flagged in sub-array 8. ! A total of 9 telescope corrections were flagged in sub-array 9. ! A total of 52 telescope corrections were flagged in sub-array 10. ! A total of 12 telescope corrections were flagged in sub-array 11. ! A total of 24 telescope corrections were flagged in sub-array 12. ! A total of 17 telescope corrections were flagged in sub-array 13. ! A total of 15 telescope corrections were flagged in sub-array 14. ! ! Correcting IF 2. ! A total of 15 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 9 telescope corrections were flagged in sub-array 8. ! A total of 9 telescope corrections were flagged in sub-array 9. ! A total of 52 telescope corrections were flagged in sub-array 10. ! A total of 12 telescope corrections were flagged in sub-array 11. ! A total of 24 telescope corrections were flagged in sub-array 12. ! A total of 17 telescope corrections were flagged in sub-array 13. ! A total of 15 telescope corrections were flagged in sub-array 14. ! ! Correcting IF 3. ! A total of 15 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 9 telescope corrections were flagged in sub-array 8. ! A total of 9 telescope corrections were flagged in sub-array 9. ! A total of 52 telescope corrections were flagged in sub-array 10. ! A total of 12 telescope corrections were flagged in sub-array 11. ! A total of 24 telescope corrections were flagged in sub-array 12. ! A total of 17 telescope corrections were flagged in sub-array 13. ! A total of 15 telescope corrections were flagged in sub-array 14. ! ! Correcting IF 4. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 20 telescope corrections were flagged in sub-array 10. ! A total of 12 telescope corrections were flagged in sub-array 11. ! A total of 24 telescope corrections were flagged in sub-array 12. ! A total of 17 telescope corrections were flagged in sub-array 13. ! A total of 15 telescope corrections were flagged in sub-array 14. ! ! Fit before self-cal, rms=0.163409Jy sigma=1.501914 ! Fit after self-cal, rms=0.161669Jy sigma=1.501385 ! Inverting map and beam ! Estimated beam: bmin=0.9937 mas, bmaj=1.311 mas, bpa=54.68 degrees ! Estimated noise=2.16239 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0695321 Jy ! Component: 100 - total flux cleaned = 0.107137 Jy ! Total flux subtracted in 100 components = 0.107137 Jy ! Clean residual min=-0.013727 max=0.030266 Jy/beam ! Clean residual mean=0.000084 rms=0.003437 Jy/beam ! Combined flux in latest and established models = 0.724253 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 20 components and 0.724253 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.119107Jy sigma=1.034066 ! Fit after self-cal, rms=0.116838Jy sigma=1.024990 ! Inverting map ! Added new window around map position (1, -1.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0327204 Jy ! Component: 100 - total flux cleaned = 0.0504824 Jy ! Total flux subtracted in 100 components = 0.0504824 Jy ! Clean residual min=-0.011847 max=0.025190 Jy/beam ! Clean residual mean=0.000057 rms=0.002794 Jy/beam ! Combined flux in latest and established models = 0.774735 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 26 components and 0.774735 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.104388Jy sigma=0.878895 ! Fit after self-cal, rms=0.103145Jy sigma=0.873744 ! Inverting map ! Added new window around map position (4.2, -3.8). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0231852 Jy ! Component: 100 - total flux cleaned = 0.0368716 Jy ! Total flux subtracted in 100 components = 0.0368716 Jy ! Clean residual min=-0.010512 max=0.011318 Jy/beam ! Clean residual mean=0.000041 rms=0.002400 Jy/beam ! Combined flux in latest and established models = 0.811607 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 38 components and 0.811607 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.095869Jy sigma=0.785932 ! Fit after self-cal, rms=0.094572Jy sigma=0.779449 ! 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.274 mas, bmaj=1.65 mas, bpa=55.34 degrees ! Estimated noise=1.29879 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.0243989 Jy ! Component: 100 - total flux cleaned = 0.039831 Jy ! Component: 150 - total flux cleaned = 0.0505856 Jy ! Component: 200 - total flux cleaned = 0.0583145 Jy ! Total flux subtracted in 200 components = 0.0583145 Jy ! Clean residual min=-0.010614 max=0.010263 Jy/beam ! Clean residual mean=0.000017 rms=0.002490 Jy/beam ! Combined flux in latest and established models = 0.869921 Jy ! Performing phase self-cal ! Adding 41 model components to the UV plane model. ! The established model now contains 67 components and 0.869921 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.089162Jy sigma=0.702352 ! Fit after self-cal, rms=0.087798Jy sigma=0.694133 ! 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 9 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* LA 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.00* NL 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* NL 0.00* ! ! A total of 9 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* HN 0.00* NL 0.00* NY 0.00* ! ! A total of 9 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 0.00* HN 0.00* NL 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! HN 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! BR 0.97 HN 1.01 LA 0.99 NL 1.05 ! PT 1.04 SC 1.15 ! ! Telescope amplitude corrections in sub-array 8: ! NY 0.00* WZ 0.00* ZC 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 9. ! Telescope amplitude corrections in sub-array 9: ! BR 0.97 FD 0.97 HN 0.99 KP 1.00 ! LA 0.94 NL 1.00 OV 1.00 PT 0.97 ! SC 1.02 ZC 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 10. ! Telescope amplitude corrections in sub-array 10: ! BR 0.96 FD 1.00 HN 0.98 KK 1.00* ! KP 1.02 LA 0.94 MK 0.98 NL 0.99 ! NY 1.00* OV 1.01 PT 0.97 ZC 1.00* ! ! Telescope amplitude corrections in sub-array 11: ! BR 0.98 FD 1.02 HN 1.00 KK 1.00* ! KP 1.01 LA 0.93 MK 0.99 NL 0.99 ! NY 1.00* OV 0.98 PT 0.95 SC 1.02 ! ! Telescope amplitude corrections in sub-array 12: ! BR 0.00* KK 0.00* MK 0.00* NY 0.00* ! ! Telescope amplitude corrections in sub-array 13: ! BR 0.97 FD 0.99 HN 0.98 KK 1.00* ! KP 1.00 LA 0.98 MK 0.95 NL 1.03 ! NY 1.00* OV 0.98 PT 0.97 ! ! Telescope amplitude corrections in sub-array 14: ! KK 0.00* MK 0.00* WZ 0.00* ! ! ! Correcting IF 2. ! A total of 9 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* LA 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.00* NL 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* NL 0.00* ! ! A total of 9 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* HN 0.00* NL 0.00* NY 0.00* ! ! A total of 9 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 0.00* HN 0.00* NL 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! HN 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! BR 1.02 HN 0.97 LA 0.98 NL 0.98 ! PT 1.11 SC 1.10 ! ! Telescope amplitude corrections in sub-array 8: ! NY 0.00* WZ 0.00* ZC 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 9. ! Telescope amplitude corrections in sub-array 9: ! BR 0.95 FD 0.97 HN 0.99 KP 0.93 ! LA 0.94 NL 0.98 OV 0.98 PT 1.01 ! SC 0.98 ZC 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 10. ! Telescope amplitude corrections in sub-array 10: ! BR 0.96 FD 0.99 HN 1.03 KK 1.00* ! KP 0.96 LA 0.92 MK 1.01 NL 0.97 ! NY 1.00* OV 1.00 PT 1.03 ZC 1.00* ! ! Telescope amplitude corrections in sub-array 11: ! BR 0.98 FD 1.01 HN 1.01 KK 1.00* ! KP 0.97 LA 0.90 MK 0.99 NL 0.98 ! NY 1.00* OV 0.97 PT 1.01 SC 1.01 ! ! Telescope amplitude corrections in sub-array 12: ! BR 0.00* KK 0.00* MK 0.00* NY 0.00* ! ! Telescope amplitude corrections in sub-array 13: ! BR 1.01 FD 0.98 HN 1.01 KK 1.00* ! KP 0.94 LA 0.98 MK 0.96 NL 0.99 ! NY 1.00* OV 0.97 PT 1.03 ! ! Telescope amplitude corrections in sub-array 14: ! KK 0.00* MK 0.00* WZ 0.00* ! ! ! Correcting IF 3. ! A total of 9 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* LA 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.00* NL 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* NL 0.00* ! ! A total of 9 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* HN 0.00* NL 0.00* NY 0.00* ! ! A total of 9 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 0.00* HN 0.00* NL 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! HN 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! BR 1.07 HN 1.07 LA 1.03 NL 1.05 ! PT 1.14 SC 1.12 ! ! Telescope amplitude corrections in sub-array 8: ! NY 0.00* WZ 0.00* ZC 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 9. ! Telescope amplitude corrections in sub-array 9: ! BR 1.04 FD 1.02 HN 0.99 KP 0.90 ! LA 1.02 NL 1.00 OV 1.03 PT 1.04 ! SC 0.96 ZC 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 10. ! Telescope amplitude corrections in sub-array 10: ! BR 1.03 FD 1.05 HN 1.00 KK 1.00* ! KP 0.90 LA 1.02 MK 1.05 NL 1.00 ! NY 1.00* OV 1.05 PT 1.04 ZC 1.00* ! ! Telescope amplitude corrections in sub-array 11: ! BR 1.05 FD 1.05 HN 0.99 KK 1.00* ! KP 0.90 LA 0.99 MK 1.00 NL 1.00 ! NY 1.00* OV 1.03 PT 1.03 SC 1.00 ! ! Telescope amplitude corrections in sub-array 12: ! BR 0.00* KK 0.00* MK 0.00* NY 0.00* ! ! Telescope amplitude corrections in sub-array 13: ! BR 1.05 FD 1.03 HN 0.99 KK 1.00* ! KP 0.88 LA 1.04 MK 0.96 NL 1.03 ! NY 1.00* OV 1.03 PT 1.04 ! ! Telescope amplitude corrections in sub-array 14: ! KK 0.00* MK 0.00* WZ 0.00* ! ! ! Correcting IF 4. ! A total of 9 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* LA 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.00* NL 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* NL 0.00* ! ! A total of 9 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* HN 0.00* NL 0.00* NY 0.00* ! ! A total of 9 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 0.00* HN 0.00* NL 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! HN 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! BR 1.04 HN 1.03 LA 1.04 NL 1.09 ! PT 1.11 SC 1.18 ! ! Telescope amplitude corrections in sub-array 8: ! NY 0.00* WZ 0.00* ZC 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 9. ! Telescope amplitude corrections in sub-array 9: ! BR 1.03 FD 1.05 HN 1.07 KP 0.95 ! LA 1.00 NL 1.02 OV 1.01 PT 1.04 ! SC 1.01 ZC 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 10. ! Telescope amplitude corrections in sub-array 10: ! BR 0.99 FD 1.06 HN 1.01 KK 1.00* ! KP 0.96 LA 1.00 MK 1.05 NL 1.05 ! NY 1.00* OV 1.04 PT 1.02 ZC 1.00* ! ! Telescope amplitude corrections in sub-array 11: ! BR 1.04 FD 1.07 HN 1.05 KK 1.00* ! KP 0.95 LA 0.97 MK 1.04 NL 1.02 ! NY 1.00* OV 1.02 PT 1.02 SC 1.05 ! ! Telescope amplitude corrections in sub-array 12: ! BR 0.00* KK 0.00* MK 0.00* NY 0.00* ! ! Telescope amplitude corrections in sub-array 13: ! BR 1.03 FD 1.05 HN 1.03 KK 1.00* ! KP 0.95 LA 1.01 MK 0.99 NL 1.05 ! NY 1.00* OV 1.02 PT 1.03 ! ! Telescope amplitude corrections in sub-array 14: ! KK 0.00* MK 0.00* WZ 0.00* ! ! ! Fit before self-cal, rms=0.087798Jy sigma=0.694133 ! Fit after self-cal, rms=0.076825Jy sigma=0.588714 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.281 mas, bmaj=1.665 mas, bpa=56.43 degrees ! Estimated noise=1.30053 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 (5.8, -2.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.010852 Jy ! Component: 100 - total flux cleaned = 0.0172554 Jy ! Component: 150 - total flux cleaned = 0.019662 Jy ! Component: 200 - total flux cleaned = 0.020897 Jy ! Total flux subtracted in 200 components = 0.020897 Jy ! Clean residual min=-0.004615 max=0.005412 Jy/beam ! Clean residual mean=0.000015 rms=0.001230 Jy/beam ! Combined flux in latest and established models = 0.890818 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 97 components and 0.890818 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.075197Jy sigma=0.569059 ! Fit after self-cal, rms=0.074876Jy sigma=0.566835 ! 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.074876Jy sigma=0.566835 ! Fit after self-cal, rms=0.074758Jy sigma=0.565584 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.28 mas, bmaj=1.663 mas, bpa=56.49 degrees ! Estimated noise=1.3008 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00586374 Jy ! Component: 100 - total flux cleaned = 0.00774389 Jy ! Component: 150 - total flux cleaned = 0.00869874 Jy ! Component: 200 - total flux cleaned = 0.00950741 Jy ! Total flux subtracted in 200 components = 0.00950741 Jy ! Clean residual min=-0.004243 max=0.004305 Jy/beam ! Clean residual mean=0.000008 rms=0.001036 Jy/beam ! Combined flux in latest and established models = 0.900325 Jy selfcal ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 112 components and 0.900325 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.074128Jy sigma=0.558657 ! Fit after self-cal, rms=0.073939Jy sigma=0.557627 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 (2.8, -2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00295393 Jy ! Component: 100 - total flux cleaned = 0.00356917 Jy ! Component: 150 - total flux cleaned = 0.00398105 Jy ! Component: 200 - total flux cleaned = 0.0043358 Jy ! Total flux subtracted in 200 components = 0.0043358 Jy ! Clean residual min=-0.004092 max=0.004030 Jy/beam ! Clean residual mean=0.000007 rms=0.000947 Jy/beam ! Combined flux in latest and established models = 0.904661 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 128 components and 0.904661 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.073549Jy sigma=0.554344 ! Fit after self-cal, rms=0.073450Jy sigma=0.553912 ! 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.073450Jy sigma=0.553912 ! Fit after self-cal, rms=0.073347Jy sigma=0.552899 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.281 mas, bmaj=1.662 mas, bpa=56.53 degrees ! Estimated noise=1.29988 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00274512 Jy ! Component: 100 - total flux cleaned = 0.00365773 Jy ! Component: 150 - total flux cleaned = 0.00416111 Jy ! Component: 200 - total flux cleaned = 0.00438022 Jy ! Total flux subtracted in 200 components = 0.00438022 Jy ! Clean residual min=-0.003729 max=0.003645 Jy/beam ! Clean residual mean=0.000002 rms=0.000869 Jy/beam ! Combined flux in latest and established models = 0.909041 Jy selfcal ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 132 components and 0.909042 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.073093Jy sigma=0.550734 ! Fit after self-cal, rms=0.073016Jy sigma=0.550388 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.000860011 Jy ! Component: 100 - total flux cleaned = 0.00145107 Jy ! Component: 150 - total flux cleaned = 0.00176938 Jy ! Component: 200 - total flux cleaned = 0.001938 Jy ! Total flux subtracted in 200 components = 0.001938 Jy ! Clean residual min=-0.003604 max=0.003473 Jy/beam ! Clean residual mean=0.000003 rms=0.000840 Jy/beam ! Combined flux in latest and established models = 0.91098 Jy ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 135 components and 0.91098 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.072874Jy sigma=0.549375 ! Fit after self-cal, rms=0.072821Jy sigma=0.549204 ! 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.072821Jy sigma=0.549204 ! Fit after self-cal, rms=0.072759Jy sigma=0.548937 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.571 mas, bmaj=2.964 mas, bpa=61.64 degrees ! Estimated noise=1.38278 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00245018 Jy ! Component: 100 - total flux cleaned = 0.0030259 Jy ! Component: 150 - total flux cleaned = 0.0032619 Jy ! Component: 200 - total flux cleaned = 0.00334064 Jy ! Total flux subtracted in 200 components = 0.00334064 Jy ! Clean residual min=-0.004264 max=0.004057 Jy/beam ! Clean residual mean=0.000000 rms=0.001034 Jy/beam ! Combined flux in latest and established models = 0.91432 Jy selfcal ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 154 components and 0.91432 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.072737Jy sigma=0.548247 ! Fit after self-cal, rms=0.072682Jy sigma=0.548064 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.000644481 Jy ! Component: 100 - total flux cleaned = 0.000973527 Jy ! Component: 150 - total flux cleaned = 0.000975789 Jy ! Component: 200 - total flux cleaned = 0.00105654 Jy ! Total flux subtracted in 200 components = 0.00105654 Jy ! Clean residual min=-0.004230 max=0.004147 Jy/beam ! Clean residual mean=0.000001 rms=0.001020 Jy/beam ! Combined flux in latest and established models = 0.915377 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 161 components and 0.915377 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.072661Jy sigma=0.547791 ! Fit after self-cal, rms=0.072625Jy sigma=0.547712 ! 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:LA,true ! 1:LA self-cal status=fixed weight=1 selfant 1:MK,true ! 1:MK self-cal status=fixed weight=1 selfant 2:BR,true ! 2:BR self-cal status=fixed weight=1 selfant 2:NL,true ! 2:NL self-cal status=fixed weight=1 selfant 3:BR,true ! 3:BR self-cal status=fixed weight=1 selfant 3:NL,true ! 3:NL self-cal status=fixed weight=1 selfant 4:BR,true ! 4:BR self-cal status=fixed weight=1 selfant 4:HN,true ! 4:HN self-cal status=fixed weight=1 selfant 4:NL,true ! 4:NL self-cal status=fixed weight=1 selfant 5:BR,true ! 5:BR self-cal status=fixed weight=1 selfant 5:HN,true ! 5:HN self-cal status=fixed weight=1 selfant 5:NL,true ! 5:NL self-cal status=fixed weight=1 selfant 6:HN,true ! 6:HN self-cal status=fixed weight=1 selfant 6:SC,true ! 6:SC self-cal status=fixed weight=1 selfant 7:BR,true ! 7:BR self-cal status=fixed weight=1 selfant 7:HN,true ! 7:HN self-cal status=fixed weight=1 selfant 7:LA,true ! 7:LA self-cal status=fixed weight=1 selfant 7:NL,true ! 7:NL self-cal status=fixed weight=1 selfant 7:PT,true ! 7:PT self-cal status=fixed weight=1 selfant 7:SC,true ! 7:SC self-cal status=fixed weight=1 selfant 9:BR,true ! 9:BR self-cal status=fixed weight=1 selfant 9:FD,true ! 9:FD self-cal status=fixed weight=1 selfant 9:HN,true ! 9:HN self-cal status=fixed weight=1 selfant 9:KP,true ! 9:KP self-cal status=fixed weight=1 selfant 9:LA,true ! 9:LA self-cal status=fixed weight=1 selfant 9:NL,true ! 9:NL self-cal status=fixed weight=1 selfant 9:OV,true ! 9:OV self-cal status=fixed weight=1 selfant 9:PT,true ! 9:PT self-cal status=fixed weight=1 selfant 9:SC,true ! 9:SC self-cal status=fixed weight=1 selfant 10:BR,true ! 10:BR self-cal status=fixed weight=1 selfant 10:FD,true ! 10:FD self-cal status=fixed weight=1 selfant 10:HN,true ! 10:HN self-cal status=fixed weight=1 selfant 10:KP,true ! 10:KP self-cal status=fixed weight=1 selfant 10:LA,true ! 10:LA self-cal status=fixed weight=1 selfant 10:MK,true ! 10:MK self-cal status=fixed weight=1 selfant 10:NL,true ! 10:NL self-cal status=fixed weight=1 selfant 10:OV,true ! 10:OV self-cal status=fixed weight=1 selfant 10:PT,true ! 10:PT self-cal status=fixed weight=1 selfant 11:BR,true ! 11:BR self-cal status=fixed weight=1 selfant 11:FD,true ! 11:FD self-cal status=fixed weight=1 selfant 11:HN,true ! 11:HN self-cal status=fixed weight=1 selfant 11:KP,true ! 11:KP self-cal status=fixed weight=1 selfant 11:LA,true ! 11:LA self-cal status=fixed weight=1 selfant 11:MK,true ! 11:MK self-cal status=fixed weight=1 selfant 11:NL,true ! 11:NL self-cal status=fixed weight=1 selfant 11:OV,true ! 11:OV self-cal status=fixed weight=1 selfant 11:PT,true ! 11:PT self-cal status=fixed weight=1 selfant 11:SC,true ! 11:SC self-cal status=fixed weight=1 selfant 12:BR,true ! 12:BR self-cal status=fixed weight=1 selfant 12:MK,true ! 12:MK self-cal status=fixed weight=1 selfant 13:BR,true ! 13:BR self-cal status=fixed weight=1 selfant 13:FD,true ! 13:FD self-cal status=fixed weight=1 selfant 13:HN,true ! 13:HN self-cal status=fixed weight=1 selfant 13:KP,true ! 13:KP self-cal status=fixed weight=1 selfant 13:LA,true ! 13:LA self-cal status=fixed weight=1 selfant 13:MK,true ! 13:MK self-cal status=fixed weight=1 selfant 13:NL,true ! 13:NL self-cal status=fixed weight=1 selfant 13:OV,true ! 13:OV self-cal status=fixed weight=1 selfant 13:PT,true ! 13:PT self-cal status=fixed weight=1 selfant 14:MK,true ! 14:MK self-cal status=fixed weight=1 ![Exited script file: /tmp/ff.ff] selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 14 telescope corrections were flagged in sub-array 2. ! A total of 5 telescope corrections were flagged in sub-array 4. ! A total of 18 telescope corrections were flagged in sub-array 8. ! A total of 1 telescope corrections were flagged in sub-array 9. ! A total of 54 telescope corrections were flagged in sub-array 10. ! A total of 14 telescope corrections were flagged in sub-array 11. ! A total of 6 telescope corrections were flagged in sub-array 12. ! A total of 8 telescope corrections were flagged in sub-array 13. ! A total of 16 telescope corrections were flagged in sub-array 14. ! ! Correcting IF 2. ! A total of 14 telescope corrections were flagged in sub-array 2. ! A total of 5 telescope corrections were flagged in sub-array 4. ! A total of 18 telescope corrections were flagged in sub-array 8. ! A total of 1 telescope corrections were flagged in sub-array 9. ! A total of 54 telescope corrections were flagged in sub-array 10. ! A total of 14 telescope corrections were flagged in sub-array 11. ! A total of 6 telescope corrections were flagged in sub-array 12. ! A total of 8 telescope corrections were flagged in sub-array 13. ! A total of 16 telescope corrections were flagged in sub-array 14. ! ! Correcting IF 3. ! A total of 14 telescope corrections were flagged in sub-array 2. ! A total of 5 telescope corrections were flagged in sub-array 4. ! A total of 18 telescope corrections were flagged in sub-array 8. ! A total of 1 telescope corrections were flagged in sub-array 9. ! A total of 54 telescope corrections were flagged in sub-array 10. ! A total of 14 telescope corrections were flagged in sub-array 11. ! A total of 6 telescope corrections were flagged in sub-array 12. ! A total of 8 telescope corrections were flagged in sub-array 13. ! A total of 16 telescope corrections were flagged in sub-array 14. ! ! Correcting IF 4. ! A total of 14 telescope corrections were flagged in sub-array 2. ! A total of 5 telescope corrections were flagged in sub-array 4. ! A total of 18 telescope corrections were flagged in sub-array 8. ! A total of 10 telescope corrections were flagged in sub-array 9. ! A total of 62 telescope corrections were flagged in sub-array 10. ! A total of 14 telescope corrections were flagged in sub-array 11. ! A total of 6 telescope corrections were flagged in sub-array 12. ! A total of 8 telescope corrections were flagged in sub-array 13. ! A total of 16 telescope corrections were flagged in sub-array 14. ! ! Fit before self-cal, rms=0.207033Jy sigma=1.010858 ! Fit after self-cal, rms=0.336946Jy sigma=0.961171 selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:MK self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:ZC self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 4:BR self-cal status=correctable weight=1 ! 4:HN self-cal status=correctable weight=1 ! 4:NL self-cal status=correctable weight=1 ! 4:NY self-cal status=correctable weight=1 ! 5:BR self-cal status=correctable weight=1 ! 5:HN self-cal status=correctable weight=1 ! 5:NL self-cal status=correctable weight=1 ! 6:HN self-cal status=correctable weight=1 ! 6:SC self-cal status=correctable weight=1 ! 7:BR self-cal status=correctable weight=1 ! 7:HN self-cal status=correctable weight=1 ! 7:LA self-cal status=correctable weight=1 ! 7:NL self-cal status=correctable weight=1 ! 7:PT self-cal status=correctable weight=1 ! 7:SC self-cal status=correctable weight=1 ! 8:NY self-cal status=correctable weight=1 ! 8:WZ self-cal status=correctable weight=1 ! 8:ZC self-cal status=correctable weight=1 ! 9:BR self-cal status=correctable weight=1 ! 9:FD self-cal status=correctable weight=1 ! 9:HN self-cal status=correctable weight=1 ! 9:KP self-cal status=correctable weight=1 ! 9:LA self-cal status=correctable weight=1 ! 9:NL self-cal status=correctable weight=1 ! 9:OV self-cal status=correctable weight=1 ! 9:PT self-cal status=correctable weight=1 ! 9:SC self-cal status=correctable weight=1 ! 9:ZC self-cal status=correctable weight=1 ! 10:BR self-cal status=correctable weight=1 ! 10:FD self-cal status=correctable weight=1 ! 10:HN self-cal status=correctable weight=1 ! 10:KK self-cal status=correctable weight=1 ! 10:KP self-cal status=correctable weight=1 ! 10:LA self-cal status=correctable weight=1 ! 10:MK self-cal status=correctable weight=1 ! 10:NL self-cal status=correctable weight=1 ! 10:NY self-cal status=correctable weight=1 ! 10:OV self-cal status=correctable weight=1 ! 10:PT self-cal status=correctable weight=1 ! 10:ZC self-cal status=correctable weight=1 ! 11:BR self-cal status=correctable weight=1 ! 11:FD self-cal status=correctable weight=1 ! 11:HN self-cal status=correctable weight=1 ! 11:KK self-cal status=correctable weight=1 ! 11:KP self-cal status=correctable weight=1 ! 11:LA self-cal status=correctable weight=1 ! 11:MK self-cal status=correctable weight=1 ! 11:NL self-cal status=correctable weight=1 ! 11:NY self-cal status=correctable weight=1 ! 11:OV self-cal status=correctable weight=1 ! 11:PT self-cal status=correctable weight=1 ! 11:SC self-cal status=correctable weight=1 ! 12:BR self-cal status=correctable weight=1 ! 12:KK self-cal status=correctable weight=1 ! 12:MK self-cal status=correctable weight=1 ! 12:NY self-cal status=correctable weight=1 ! 13:BR self-cal status=correctable weight=1 ! 13:FD self-cal status=correctable weight=1 ! 13:HN self-cal status=correctable weight=1 ! 13:KK self-cal status=correctable weight=1 ! 13:KP self-cal status=correctable weight=1 ! 13:LA self-cal status=correctable weight=1 ! 13:MK self-cal status=correctable weight=1 ! 13:NL self-cal status=correctable weight=1 ! 13:NY self-cal status=correctable weight=1 ! 13:OV self-cal status=correctable weight=1 ! 13:PT self-cal status=correctable weight=1 ! 14:KK self-cal status=correctable weight=1 ! 14:MK self-cal status=correctable weight=1 ! 14:WZ 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 18 telescope corrections were flagged in sub-array 1. ! A total of 28 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 15 telescope corrections were flagged in sub-array 4. ! A total of 12 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 8 telescope corrections were flagged in sub-array 7. ! A total of 17 telescope corrections were flagged in sub-array 9. ! A total of 171 telescope corrections were flagged in sub-array 10. ! A total of 35 telescope corrections were flagged in sub-array 11. ! A total of 6 telescope corrections were flagged in sub-array 12. ! A total of 38 telescope corrections were flagged in sub-array 13. ! A total of 8 telescope corrections were flagged in sub-array 14. ! ! Correcting IF 2. ! A total of 18 telescope corrections were flagged in sub-array 1. ! A total of 28 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 15 telescope corrections were flagged in sub-array 4. ! A total of 12 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 8 telescope corrections were flagged in sub-array 7. ! A total of 17 telescope corrections were flagged in sub-array 9. ! A total of 171 telescope corrections were flagged in sub-array 10. ! A total of 35 telescope corrections were flagged in sub-array 11. ! A total of 6 telescope corrections were flagged in sub-array 12. ! A total of 38 telescope corrections were flagged in sub-array 13. ! A total of 8 telescope corrections were flagged in sub-array 14. ! ! Correcting IF 3. ! A total of 18 telescope corrections were flagged in sub-array 1. ! A total of 28 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 15 telescope corrections were flagged in sub-array 4. ! A total of 12 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 8 telescope corrections were flagged in sub-array 7. ! A total of 17 telescope corrections were flagged in sub-array 9. ! A total of 171 telescope corrections were flagged in sub-array 10. ! A total of 35 telescope corrections were flagged in sub-array 11. ! A total of 6 telescope corrections were flagged in sub-array 12. ! A total of 38 telescope corrections were flagged in sub-array 13. ! A total of 8 telescope corrections were flagged in sub-array 14. ! ! Correcting IF 4. ! A total of 18 telescope corrections were flagged in sub-array 1. ! A total of 28 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 15 telescope corrections were flagged in sub-array 4. ! A total of 12 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 8 telescope corrections were flagged in sub-array 7. ! A total of 20 telescope corrections were flagged in sub-array 9. ! A total of 177 telescope corrections were flagged in sub-array 10. ! A total of 35 telescope corrections were flagged in sub-array 11. ! A total of 6 telescope corrections were flagged in sub-array 12. ! A total of 38 telescope corrections were flagged in sub-array 13. ! A total of 8 telescope corrections were flagged in sub-array 14. ! ! Fit before self-cal, rms=0.336946Jy sigma=0.961171 ! Fit after self-cal, rms=0.317882Jy sigma=0.522653 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 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 milli-arcsec ! Estimated beam: bmin=0.9064 mas, bmaj=0.9582 mas, bpa=62.25 degrees ! Estimated noise=2.41779 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 (-25, 10). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00395893 Jy ! Component: 100 - total flux cleaned = 0.00395498 Jy ! Component: 150 - total flux cleaned = 0.00395522 Jy ! Component: 200 - total flux cleaned = 0.00378393 Jy ! Total flux subtracted in 200 components = 0.00378393 Jy ! Clean residual min=-0.005597 max=0.005230 Jy/beam ! Clean residual mean=0.000000 rms=0.001206 Jy/beam ! Combined flux in latest and established models = 0.919161 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 174 components and 0.919161 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.317637Jy sigma=0.521629 ! Fit after self-cal, rms=0.317534Jy sigma=0.520075 ! Inverting map ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 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 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 milli-arcsec ! Estimated beam: bmin=1.156 mas, bmaj=1.246 mas, bpa=66.14 degrees ! Estimated noise=1.33934 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 (5.6, -5.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00171053 Jy ! Component: 100 - total flux cleaned = 0.00149907 Jy ! Component: 150 - total flux cleaned = 0.00140925 Jy ! Component: 200 - total flux cleaned = 0.0014106 Jy ! Component: 250 - total flux cleaned = 0.00133716 Jy ! Component: 300 - total flux cleaned = 0.00133651 Jy ! Component: 350 - total flux cleaned = 0.00133707 Jy ! Component: 400 - total flux cleaned = 0.00133906 Jy ! Total flux subtracted in 400 components = 0.00133906 Jy ! Clean residual min=-0.002831 max=0.002999 Jy/beam ! Clean residual mean=-0.000000 rms=0.000725 Jy/beam ! Combined flux in latest and established models = 0.9205 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 200 components and 0.9205 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.317300Jy sigma=0.518456 ! Fit after self-cal, rms=0.317297Jy sigma=0.518272 ! Inverting map ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 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.000409402 Jy ! Component: 100 - total flux cleaned = 0.000833996 Jy ! Component: 150 - total flux cleaned = 0.00089747 Jy ! Component: 200 - total flux cleaned = 0.000959418 Jy ! Component: 250 - total flux cleaned = 0.000960252 Jy ! Component: 300 - total flux cleaned = 0.00102072 Jy ! Component: 350 - total flux cleaned = 0.000961203 Jy ! Component: 400 - total flux cleaned = 0.000903887 Jy ! Total flux subtracted in 400 components = 0.000903887 Jy ! Clean residual min=-0.002801 max=0.003059 Jy/beam ! Clean residual mean=-0.000000 rms=0.000715 Jy/beam ! Combined flux in latest and established models = 0.921404 Jy ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 211 components and 0.921404 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.317108Jy sigma=0.517644 ! Fit after self-cal, rms=0.317120Jy sigma=0.517551 ! Inverting map ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 milli-arcsec ! Added new window around map position (14.4, -11.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0032649 Jy ! Component: 100 - total flux cleaned = 0.00462558 Jy ! Component: 150 - total flux cleaned = 0.00562793 Jy ! Component: 200 - total flux cleaned = 0.00674699 Jy ! Component: 250 - total flux cleaned = 0.00758962 Jy ! Component: 300 - total flux cleaned = 0.00807785 Jy ! Component: 350 - total flux cleaned = 0.00848957 Jy ! Component: 400 - total flux cleaned = 0.00860221 Jy ! Total flux subtracted in 400 components = 0.00860221 Jy ! Clean residual min=-0.002812 max=0.002661 Jy/beam ! Clean residual mean=-0.000002 rms=0.000664 Jy/beam ! Combined flux in latest and established models = 0.930006 Jy ! Performing phase self-cal ! Adding 53 model components to the UV plane model. ! The established model now contains 238 components and 0.930006 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.316954Jy sigma=0.515938 ! Fit after self-cal, rms=0.316946Jy sigma=0.515790 ! Inverting map ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 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. ! 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 18 telescope corrections were flagged in sub-array 1. ! A total of 42 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 20 telescope corrections were flagged in sub-array 4. ! A total of 12 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 8 telescope corrections were flagged in sub-array 7. ! A total of 18 telescope corrections were flagged in sub-array 8. ! A total of 18 telescope corrections were flagged in sub-array 9. ! A total of 225 telescope corrections were flagged in sub-array 10. ! A total of 49 telescope corrections were flagged in sub-array 11. ! A total of 12 telescope corrections were flagged in sub-array 12. ! A total of 46 telescope corrections were flagged in sub-array 13. ! A total of 24 telescope corrections were flagged in sub-array 14. ! ! Correcting IF 2. ! A total of 18 telescope corrections were flagged in sub-array 1. ! A total of 42 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 20 telescope corrections were flagged in sub-array 4. ! A total of 12 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 8 telescope corrections were flagged in sub-array 7. ! A total of 18 telescope corrections were flagged in sub-array 8. ! A total of 18 telescope corrections were flagged in sub-array 9. ! A total of 225 telescope corrections were flagged in sub-array 10. ! A total of 49 telescope corrections were flagged in sub-array 11. ! A total of 12 telescope corrections were flagged in sub-array 12. ! A total of 46 telescope corrections were flagged in sub-array 13. ! A total of 24 telescope corrections were flagged in sub-array 14. ! ! Correcting IF 3. ! A total of 18 telescope corrections were flagged in sub-array 1. ! A total of 42 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 20 telescope corrections were flagged in sub-array 4. ! A total of 12 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 8 telescope corrections were flagged in sub-array 7. ! A total of 18 telescope corrections were flagged in sub-array 8. ! A total of 18 telescope corrections were flagged in sub-array 9. ! A total of 225 telescope corrections were flagged in sub-array 10. ! A total of 49 telescope corrections were flagged in sub-array 11. ! A total of 12 telescope corrections were flagged in sub-array 12. ! A total of 46 telescope corrections were flagged in sub-array 13. ! A total of 24 telescope corrections were flagged in sub-array 14. ! ! Correcting IF 4. ! A total of 18 telescope corrections were flagged in sub-array 1. ! A total of 42 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 20 telescope corrections were flagged in sub-array 4. ! A total of 12 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 8 telescope corrections were flagged in sub-array 7. ! A total of 18 telescope corrections were flagged in sub-array 8. ! A total of 30 telescope corrections were flagged in sub-array 9. ! A total of 239 telescope corrections were flagged in sub-array 10. ! A total of 49 telescope corrections were flagged in sub-array 11. ! A total of 12 telescope corrections were flagged in sub-array 12. ! A total of 46 telescope corrections were flagged in sub-array 13. ! A total of 24 telescope corrections were flagged in sub-array 14. ! ! Fit before self-cal, rms=0.370527Jy sigma=0.953859 ! Fit after self-cal, rms=0.310457Jy sigma=0.515510 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 milli-arcsec ! Estimated beam: bmin=1.155 mas, bmaj=1.246 mas, bpa=65.15 degrees ! Estimated noise=1.33823 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000987262 Jy ! Component: 100 - total flux cleaned = 0.00154826 Jy ! Component: 150 - total flux cleaned = 0.00212402 Jy ! Component: 200 - total flux cleaned = 0.00254817 Jy ! Component: 250 - total flux cleaned = 0.00295432 Jy ! Component: 300 - total flux cleaned = 0.00317865 Jy ! Component: 350 - total flux cleaned = 0.00350356 Jy ! Component: 400 - total flux cleaned = 0.00376557 Jy ! Total flux subtracted in 400 components = 0.00376557 Jy ! Clean residual min=-0.002766 max=0.002530 Jy/beam ! Clean residual mean=0.000001 rms=0.000639 Jy/beam ! Combined flux in latest and established models = 0.933772 Jy selfcal ! Performing phase self-cal ! Adding 51 model components to the UV plane model. ! The established model now contains 248 components and 0.933772 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.310293Jy sigma=0.514805 ! Fit after self-cal, rms=0.310313Jy sigma=0.514739 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 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.000372337 Jy ! Component: 100 - total flux cleaned = 0.000861198 Jy ! Component: 150 - total flux cleaned = 0.00117332 Jy ! Component: 200 - total flux cleaned = 0.001375 Jy ! Component: 250 - total flux cleaned = 0.00152367 Jy ! Component: 300 - total flux cleaned = 0.00171762 Jy ! Component: 350 - total flux cleaned = 0.00186069 Jy ! Component: 400 - total flux cleaned = 0.00195451 Jy ! Total flux subtracted in 400 components = 0.00195451 Jy ! Clean residual min=-0.002735 max=0.002571 Jy/beam ! Clean residual mean=0.000001 rms=0.000630 Jy/beam ! Combined flux in latest and established models = 0.935726 Jy ! Performing phase self-cal ! Adding 68 model components to the UV plane model. ! The established model now contains 260 components and 0.935726 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.310190Jy sigma=0.514298 ! Fit after self-cal, rms=0.310220Jy sigma=0.514252 ! Inverting map ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 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.310220Jy sigma=0.514252 ! Fit after self-cal, rms=0.330563Jy sigma=0.514177 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 milli-arcsec ! Estimated beam: bmin=1.154 mas, bmaj=1.247 mas, bpa=64.38 degrees ! Estimated noise=1.33771 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000207008 Jy ! Component: 100 - total flux cleaned = 0.00048098 Jy ! Component: 150 - total flux cleaned = 0.000791432 Jy ! Component: 200 - total flux cleaned = 0.00113938 Jy ! Component: 250 - total flux cleaned = 0.0013298 Jy ! Component: 300 - total flux cleaned = 0.00146814 Jy ! Component: 350 - total flux cleaned = 0.00155826 Jy ! Component: 400 - total flux cleaned = 0.00169068 Jy ! Total flux subtracted in 400 components = 0.00169068 Jy ! Clean residual min=-0.002704 max=0.002647 Jy/beam ! Clean residual mean=0.000002 rms=0.000620 Jy/beam ! Combined flux in latest and established models = 0.937417 Jy selfcal ! Performing phase self-cal ! Adding 65 model components to the UV plane model. ! The established model now contains 269 components and 0.937417 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.330438Jy sigma=0.513775 ! Fit after self-cal, rms=0.330493Jy sigma=0.513736 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 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.000259566 Jy ! Component: 100 - total flux cleaned = 0.000530842 Jy ! Component: 150 - total flux cleaned = 0.000662874 Jy ! Component: 200 - total flux cleaned = 0.000792053 Jy ! Component: 250 - total flux cleaned = 0.000792757 Jy ! Component: 300 - total flux cleaned = 0.000918422 Jy ! Component: 350 - total flux cleaned = 0.00100057 Jy ! Component: 400 - total flux cleaned = 0.00108283 Jy ! Total flux subtracted in 400 components = 0.00108283 Jy ! Clean residual min=-0.002681 max=0.002660 Jy/beam ! Clean residual mean=0.000002 rms=0.000615 Jy/beam ! Combined flux in latest and established models = 0.9385 Jy ! Performing phase self-cal ! Adding 78 model components to the UV plane model. ! The established model now contains 283 components and 0.9385 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.330400Jy sigma=0.513442 ! Fit after self-cal, rms=0.330438Jy sigma=0.513412 ! Inverting map ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 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.330438Jy sigma=0.513412 ! Fit after self-cal, rms=0.422090Jy sigma=0.494862 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 milli-arcsec ! Estimated beam: bmin=2.534 mas, bmaj=2.936 mas, bpa=60.26 degrees ! Estimated noise=1.37712 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0011202 Jy ! Component: 100 - total flux cleaned = 0.00162015 Jy ! Component: 150 - total flux cleaned = 0.00195959 Jy ! Component: 200 - total flux cleaned = 0.00231291 Jy ! Component: 250 - total flux cleaned = 0.00248505 Jy ! Component: 300 - total flux cleaned = 0.00262357 Jy ! Component: 350 - total flux cleaned = 0.0027327 Jy ! Component: 400 - total flux cleaned = 0.00279448 Jy ! Total flux subtracted in 400 components = 0.00279448 Jy ! Clean residual min=-0.002808 max=0.002845 Jy/beam ! Clean residual mean=0.000005 rms=0.000693 Jy/beam ! Combined flux in latest and established models = 0.941294 Jy selfcal ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 324 components and 0.941294 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.422084Jy sigma=0.494596 ! Fit after self-cal, rms=0.422020Jy sigma=0.494570 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 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 (14.8, -5.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00266236 Jy ! Component: 100 - total flux cleaned = 0.00360515 Jy ! Component: 150 - total flux cleaned = 0.00437381 Jy ! Component: 200 - total flux cleaned = 0.00493031 Jy ! Component: 250 - total flux cleaned = 0.0053226 Jy ! Component: 300 - total flux cleaned = 0.00565213 Jy ! Component: 350 - total flux cleaned = 0.00589455 Jy ! Component: 400 - total flux cleaned = 0.00612253 Jy ! Total flux subtracted in 400 components = 0.00612253 Jy ! Clean residual min=-0.002449 max=0.002275 Jy/beam ! Clean residual mean=0.000009 rms=0.000624 Jy/beam ! Combined flux in latest and established models = 0.947417 Jy ! Performing phase self-cal ! Adding 52 model components to the UV plane model. ! The established model now contains 373 components and 0.947417 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.422055Jy sigma=0.493718 ! Fit after self-cal, rms=0.422022Jy sigma=0.493668 ! Inverting map ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 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 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 milli-arcsec ! Estimated beam: bmin=1.154 mas, bmaj=1.248 mas, bpa=64.25 degrees ! Estimated noise=1.33784 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.154 x 1.248 at 64.25 degrees (North through East) ! Clean map min=-0.0022813 max=0.61305 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.422022Jy sigma=0.493668 ! Fit after self-cal, rms=1.417334Jy sigma=0.447501 ! Inverting map and beam ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 milli-arcsec ! Estimated beam: bmin=1.154 mas, bmaj=1.247 mas, bpa=64.33 degrees ! Estimated noise=1.33703 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000263994 Jy ! Component: 100 - total flux cleaned = 0.000632507 Jy ! Component: 150 - total flux cleaned = 0.000905314 Jy ! Component: 200 - total flux cleaned = 0.000967312 Jy ! Component: 250 - total flux cleaned = 0.0010274 Jy ! Component: 300 - total flux cleaned = 0.00125944 Jy ! Component: 350 - total flux cleaned = 0.00137128 Jy ! Component: 400 - total flux cleaned = 0.00137024 Jy ! Total flux subtracted in 400 components = 0.00137024 Jy ! Clean residual min=-0.002316 max=0.002240 Jy/beam ! Clean residual mean=0.000003 rms=0.000543 Jy/beam ! Combined flux in latest and established models = 0.948787 Jy ! Performing phase self-cal ! Adding 56 model components to the UV plane model. ! The established model now contains 422 components and 0.948787 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.417273Jy sigma=0.446511 ! Fit after self-cal, rms=1.417296Jy sigma=0.446486 ! Inverting map ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 milli-arcsec ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 5.4581e-05 Jy ! Component: 100 - total flux cleaned = 0.000212987 Jy ! Component: 150 - total flux cleaned = 0.000212798 Jy ! Component: 200 - total flux cleaned = 0.000213513 Jy ! Component: 250 - total flux cleaned = 0.00021377 Jy ! Component: 300 - total flux cleaned = 0.000263381 Jy ! Component: 350 - total flux cleaned = 0.000310808 Jy ! Component: 400 - total flux cleaned = 0.000357558 Jy ! Total flux subtracted in 400 components = 0.000357558 Jy ! Clean residual min=-0.002284 max=0.002176 Jy/beam ! Clean residual mean=0.000003 rms=0.000535 Jy/beam ! Combined flux in latest and established models = 0.949145 Jy ! Performing phase self-cal ! Adding 87 model components to the UV plane model. ! The established model now contains 453 components and 0.949145 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.417246Jy sigma=0.445994 ! Fit after self-cal, rms=1.417328Jy sigma=0.445985 ! Inverting map ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 milli-arcsec ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.417328Jy sigma=0.445985 ! Fit after self-cal, rms=1.370178Jy sigma=0.445941 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 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 milli-arcsec ! Estimated beam: bmin=0.9054 mas, bmaj=0.9595 mas, bpa=59.41 degrees ! Estimated noise=2.41582 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.475086 Jy ! Component: 100 - total flux cleaned = 0.619035 Jy ! Component: 150 - total flux cleaned = 0.693232 Jy ! Component: 200 - total flux cleaned = 0.737153 Jy ! Component: 250 - total flux cleaned = 0.767145 Jy ! Component: 300 - total flux cleaned = 0.78954 Jy ! Component: 350 - total flux cleaned = 0.807143 Jy ! Component: 400 - total flux cleaned = 0.82142 Jy ! Total flux subtracted in 400 components = 0.82142 Jy ! Clean residual min=-0.006094 max=0.008541 Jy/beam ! Clean residual mean=0.000027 rms=0.001470 Jy/beam ! Combined flux in latest and established models = 0.82142 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 55 model components to the UV plane model. ! The established model now contains 55 components and 0.82142 Jy ! Inverting map and beam ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 milli-arcsec ! Estimated beam: bmin=1.154 mas, bmaj=1.247 mas, bpa=64.1 degrees ! Estimated noise=1.3371 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam keep new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0240207 Jy ! Component: 100 - total flux cleaned = 0.0404681 Jy ! Component: 150 - total flux cleaned = 0.0527028 Jy ! Component: 200 - total flux cleaned = 0.0621929 Jy ! Component: 250 - total flux cleaned = 0.0696588 Jy ! Component: 300 - total flux cleaned = 0.0756193 Jy ! Component: 350 - total flux cleaned = 0.0805787 Jy ! Component: 400 - total flux cleaned = 0.0848931 Jy ! Total flux subtracted in 400 components = 0.0848931 Jy ! Clean residual min=-0.003296 max=0.003629 Jy/beam ! Clean residual mean=0.000003 rms=0.000842 Jy/beam ! Combined flux in latest and established models = 0.906313 Jy ! Adding 122 model components to the UV plane model. ! The established model now contains 153 components and 0.906313 Jy ! Inverting map ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 milli-arcsec ! Added new window around map position (5.2, 11). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00407279 Jy ! Component: 100 - total flux cleaned = 0.00765408 Jy ! Component: 150 - total flux cleaned = 0.0109471 Jy ! Component: 200 - total flux cleaned = 0.0137485 Jy ! Component: 250 - total flux cleaned = 0.0160156 Jy ! Component: 300 - total flux cleaned = 0.0180301 Jy ! Component: 350 - total flux cleaned = 0.0199211 Jy ! Component: 400 - total flux cleaned = 0.0216988 Jy ! Total flux subtracted in 400 components = 0.0216988 Jy ! Clean residual min=-0.002646 max=0.002663 Jy/beam ! Clean residual mean=0.000004 rms=0.000650 Jy/beam ! Combined flux in latest and established models = 0.928012 Jy ! Adding 77 model components to the UV plane model. ! The established model now contains 193 components and 0.928012 Jy ! Inverting map ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 milli-arcsec ! Added new window around map position (-39.4, 35.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00231133 Jy ! Component: 100 - total flux cleaned = 0.00360334 Jy ! Component: 150 - total flux cleaned = 0.00482664 Jy ! Component: 200 - total flux cleaned = 0.00606668 Jy ! Component: 250 - total flux cleaned = 0.00725054 Jy ! Component: 300 - total flux cleaned = 0.00824288 Jy ! Component: 350 - total flux cleaned = 0.00919909 Jy ! Component: 400 - total flux cleaned = 0.009924 Jy ! Total flux subtracted in 400 components = 0.009924 Jy ! Clean residual min=-0.002358 max=0.002275 Jy/beam ! Clean residual mean=0.000001 rms=0.000584 Jy/beam ! Combined flux in latest and established models = 0.937936 Jy ! Adding 92 model components to the UV plane model. ! The established model now contains 243 components and 0.937936 Jy ! Inverting map ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 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 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 milli-arcsec ! Estimated beam: bmin=2.534 mas, bmaj=2.937 mas, bpa=60.05 degrees ! Estimated noise=1.37645 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00206936 Jy ! Component: 100 - total flux cleaned = 0.00359272 Jy ! Component: 150 - total flux cleaned = 0.00464287 Jy ! Component: 200 - total flux cleaned = 0.00555283 Jy ! Component: 250 - total flux cleaned = 0.00629867 Jy ! Component: 300 - total flux cleaned = 0.00704664 Jy ! Component: 350 - total flux cleaned = 0.00775128 Jy ! Component: 400 - total flux cleaned = 0.00841785 Jy ! Total flux subtracted in 400 components = 0.00841785 Jy ! Clean residual min=-0.002459 max=0.002277 Jy/beam ! Clean residual mean=0.000004 rms=0.000605 Jy/beam ! Combined flux in latest and established models = 0.946354 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.000554852 Jy ! Component: 100 - total flux cleaned = 0.00108341 Jy ! Component: 150 - total flux cleaned = 0.00151686 Jy ! Component: 200 - total flux cleaned = 0.0018642 Jy ! Component: 250 - total flux cleaned = 0.00213349 Jy ! Component: 300 - total flux cleaned = 0.00245827 Jy ! Component: 350 - total flux cleaned = 0.00270986 Jy ! Component: 400 - total flux cleaned = 0.00283222 Jy ! Total flux subtracted in 400 components = 0.00283222 Jy ! Clean residual min=-0.002265 max=0.002294 Jy/beam ! Clean residual mean=0.000006 rms=0.000579 Jy/beam ! Combined flux in latest and established models = 0.949186 Jy ! Adding 148 model components to the UV plane model. ! The established model now contains 380 components and 0.949186 Jy ! Inverting map ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 milli-arcsec ! Added new window around map position (21.2, -44.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00205764 Jy ! Component: 100 - total flux cleaned = 0.00235735 Jy ! Component: 150 - total flux cleaned = 0.00247267 Jy ! Component: 200 - total flux cleaned = 0.00265724 Jy ! Component: 250 - total flux cleaned = 0.00276363 Jy ! Component: 300 - total flux cleaned = 0.00290008 Jy ! Component: 350 - total flux cleaned = 0.00299937 Jy ! Component: 400 - total flux cleaned = 0.00318928 Jy ! Total flux subtracted in 400 components = 0.00318928 Jy ! Clean residual min=-0.002232 max=0.001915 Jy/beam ! Clean residual mean=0.000009 rms=0.000546 Jy/beam ! Combined flux in latest and established models = 0.952375 Jy ! Adding 73 model components to the UV plane model. ! The established model now contains 436 components and 0.952375 Jy ! Inverting map ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 milli-arcsec uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 milli-arcsec ! Estimated beam: bmin=1.154 mas, bmaj=1.247 mas, bpa=64.1 degrees ! Estimated noise=1.3371 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000390338 Jy ! Component: 100 - total flux cleaned = 0.000199403 Jy ! Component: 150 - total flux cleaned = -4.05666e-05 Jy ! Component: 200 - total flux cleaned = -0.000213378 Jy ! Component: 250 - total flux cleaned = -0.000325783 Jy ! Component: 300 - total flux cleaned = -0.000434101 Jy ! Component: 350 - total flux cleaned = -0.000537982 Jy ! Component: 400 - total flux cleaned = -0.000536691 Jy ! Total flux subtracted in 400 components = -0.000536691 Jy ! Clean residual min=-0.002287 max=0.002111 Jy/beam ! Clean residual mean=0.000004 rms=0.000512 Jy/beam ! Combined flux in latest and established models = 0.951838 Jy selfcal ! Performing phase self-cal ! Adding 76 model components to the UV plane model. ! The established model now contains 498 components and 0.951838 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.369372Jy sigma=0.444453 ! Fit after self-cal, rms=1.370421Jy sigma=0.444169 wmodel /scr/pima/rdv62_uvs/J1419+5423_X_map.mod ! Writing 498 model components to file: /scr/pima/rdv62_uvs/J1419+5423_X_map.mod wobs /scr/pima/rdv62_uvs/J1419+5423_X_uvs.fits ! Applying 418 buffered edits. ! Writing UV FITS file: /scr/pima/rdv62_uvs/J1419+5423_X_uvs.fits wwins /scr/pima/rdv62_uvs/J1419+5423_X_map.win ! wwins: Wrote 12 windows to /scr/pima/rdv62_uvs/J1419+5423_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 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 milli-arcsec ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.000172643 Jy ! Component: 100 - total flux cleaned = -0.000173584 Jy ! Component: 150 - total flux cleaned = -0.000174352 Jy ! Component: 200 - total flux cleaned = -2.85146e-05 Jy ! Component: 250 - total flux cleaned = 0.00015001 Jy ! Component: 300 - total flux cleaned = 0.000290155 Jy ! Component: 350 - total flux cleaned = 0.000427607 Jy ! Component: 400 - total flux cleaned = 0.00063063 Jy ! Component: 450 - total flux cleaned = 0.000797727 Jy ! Component: 500 - total flux cleaned = 0.000962883 Jy ! Component: 550 - total flux cleaned = 0.00109328 Jy ! Component: 600 - total flux cleaned = 0.00131921 Jy ! Component: 650 - total flux cleaned = 0.00151045 Jy ! Component: 700 - total flux cleaned = 0.00166816 Jy ! Component: 750 - total flux cleaned = 0.00169938 Jy ! Component: 800 - total flux cleaned = 0.00182406 Jy ! Component: 850 - total flux cleaned = 0.00200847 Jy ! Component: 900 - total flux cleaned = 0.00209948 Jy ! Component: 950 - total flux cleaned = 0.00219062 Jy ! Component: 1000 - total flux cleaned = 0.00237088 Jy ! Component: 1050 - total flux cleaned = 0.00246021 Jy ! Component: 1100 - total flux cleaned = 0.00266661 Jy ! Component: 1150 - total flux cleaned = 0.00269552 Jy ! Component: 1200 - total flux cleaned = 0.0026664 Jy ! Component: 1250 - total flux cleaned = 0.00295556 Jy ! Component: 1300 - total flux cleaned = 0.00292668 Jy ! Component: 1350 - total flux cleaned = 0.00312624 Jy ! Component: 1400 - total flux cleaned = 0.00318262 Jy ! Component: 1450 - total flux cleaned = 0.00321048 Jy ! Component: 1500 - total flux cleaned = 0.00323864 Jy ! Component: 1550 - total flux cleaned = 0.0032945 Jy ! Component: 1600 - total flux cleaned = 0.00334993 Jy ! Component: 1650 - total flux cleaned = 0.00340523 Jy ! Component: 1700 - total flux cleaned = 0.00340515 Jy ! Component: 1750 - total flux cleaned = 0.00356823 Jy ! Component: 1800 - total flux cleaned = 0.0036762 Jy ! Component: 1850 - total flux cleaned = 0.00375678 Jy ! Component: 1900 - total flux cleaned = 0.00375671 Jy ! Component: 1950 - total flux cleaned = 0.00383669 Jy ! Component: 2000 - total flux cleaned = 0.00396917 Jy ! Component: 2050 - total flux cleaned = 0.00415341 Jy ! Component: 2100 - total flux cleaned = 0.00420578 Jy ! Component: 2150 - total flux cleaned = 0.00428395 Jy ! Component: 2200 - total flux cleaned = 0.00443936 Jy ! Component: 2250 - total flux cleaned = 0.0044911 Jy ! Component: 2300 - total flux cleaned = 0.00459392 Jy ! Component: 2350 - total flux cleaned = 0.0047473 Jy ! Component: 2400 - total flux cleaned = 0.00484893 Jy ! Component: 2450 - total flux cleaned = 0.00489945 Jy ! Component: 2500 - total flux cleaned = 0.00507611 Jy ! Component: 2550 - total flux cleaned = 0.00505093 Jy ! Component: 2600 - total flux cleaned = 0.00517576 Jy ! Component: 2650 - total flux cleaned = 0.00522572 Jy ! Component: 2700 - total flux cleaned = 0.00527526 Jy ! Component: 2750 - total flux cleaned = 0.00542347 Jy ! Component: 2800 - total flux cleaned = 0.00544809 Jy ! Component: 2850 - total flux cleaned = 0.00552169 Jy ! Component: 2900 - total flux cleaned = 0.00559482 Jy ! Component: 2950 - total flux cleaned = 0.00547334 Jy ! Component: 3000 - total flux cleaned = 0.00554603 Jy ! Component: 3050 - total flux cleaned = 0.00571472 Jy ! Component: 3100 - total flux cleaned = 0.00569076 Jy ! Component: 3150 - total flux cleaned = 0.0057387 Jy ! Component: 3200 - total flux cleaned = 0.00576207 Jy ! Component: 3250 - total flux cleaned = 0.00573839 Jy ! Component: 3300 - total flux cleaned = 0.00576218 Jy ! Component: 3350 - total flux cleaned = 0.00571517 Jy ! Component: 3400 - total flux cleaned = 0.00573858 Jy ! Component: 3450 - total flux cleaned = 0.00580867 Jy ! Component: 3500 - total flux cleaned = 0.00587843 Jy ! Component: 3550 - total flux cleaned = 0.0058552 Jy ! Component: 3600 - total flux cleaned = 0.00592468 Jy ! Component: 3650 - total flux cleaned = 0.00590159 Jy ! Component: 3700 - total flux cleaned = 0.00594752 Jy ! Component: 3750 - total flux cleaned = 0.00606211 Jy ! Component: 3800 - total flux cleaned = 0.00601634 Jy ! Component: 3850 - total flux cleaned = 0.00597084 Jy ! Component: 3900 - total flux cleaned = 0.00610713 Jy ! Component: 3950 - total flux cleaned = 0.00606203 Jy ! Component: 4000 - total flux cleaned = 0.00603961 Jy ! Component: 4050 - total flux cleaned = 0.00615216 Jy ! Total flux subtracted in 4096 components = 0.00621929 Jy ! Clean residual min=-0.001251 max=0.001297 Jy/beam ! Clean residual mean=0.000001 rms=0.000378 Jy/beam ! Combined flux in latest and established models = 0.958058 Jy keep ! Adding 324 model components to the UV plane model. ! The established model now contains 822 components and 0.958058 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv62_uvs/J1419+5423_X\_map.fits" ! Inverting map ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 milli-arcsec ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.154 x 1.247 at 64.1 degrees (North through East) ! Clean map min=-0.001919 max=0.61303 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv62_uvs/J1419+5423_X_map.fits float peak_flux peak_flux = peak(flux,max) invert ! Inverting map ! Your choice of large map pixels excluded 3.42% of the data. ! The x-axis pixel size should ideally be below 0.168 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= 0.613 rms= 0.00037 Jy quit ! Quitting program ! Log file /scr/pima/rdv62_uvs/J1419+5423_X_uvs.log closed on Fri Nov 22 22:39:19 2013