! Started logfile: /scr/pima/rdv73_uvs/J1058+8114_X_uvs.log on Sat Jan 9 17:33:46 2016 obs /scr/pima/rdv73_uvs/J1058+8114_X_uva.fits ! Reading UV FITS file: /scr/pima/rdv73_uvs/J1058+8114_X_uva.fits ! AN table 1: 14 integrations on 45 of 45 possible baselines. ! AN table 2: 24 integrations on 15 of 15 possible baselines. ! AN table 3: 3 integrations on 3 of 3 possible baselines. ! AN table 4: 3 integrations on 3 of 3 possible baselines. ! AN table 5: 3 integrations on 1 of 1 possible baselines. ! AN table 6: 25 integrations on 91 of 91 possible baselines. ! AN table 7: 19 integrations on 3 of 3 possible baselines. ! AN table 8: 31 integrations on 78 of 78 possible baselines. ! Apparent sampling: 0.371116 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1058+8114 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 8.40599e+09 1.6e+07 1 1.6e+07 ! 02 2 8.47599e+09 1.6e+07 1 1.6e+07 ! 03 3 8.79099e+09 1.6e+07 1 1.6e+07 ! 04 4 8.89599e+09 1.6e+07 1 1.6e+07 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 8552 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/rdv73_uvs/J1058+8114_X] integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 800000 float soltime2 soltime2 = 15 float soltime3 soltime3 = 1 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float old_peak float new_peak float flux_cutoff selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:KP self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:MK self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:OV self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 1:ZC self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:HN self-cal status=correctable weight=1 ! 2:KK self-cal status=correctable weight=1 ! 2:NY 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:KK self-cal status=correctable weight=1 ! 3:MK self-cal status=correctable weight=1 ! 4:BR self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 5:WZ self-cal status=correctable weight=1 ! 5:ZC self-cal status=correctable weight=1 ! 6:BR self-cal status=correctable weight=1 ! 6:FD self-cal status=correctable weight=1 ! 6:HN 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:NL self-cal status=correctable weight=1 ! 6:NY self-cal status=correctable weight=1 ! 6:OV self-cal status=correctable weight=1 ! 6:PT self-cal status=correctable weight=1 ! 6:SC self-cal status=correctable weight=1 ! 6:WZ self-cal status=correctable weight=1 ! 6:ZC self-cal status=correctable weight=1 ! 7:KK self-cal status=correctable weight=1 ! 7:MK self-cal status=correctable weight=1 ! 7:ZC self-cal status=correctable weight=1 ! 8:BR self-cal status=correctable weight=1 ! 8:FD self-cal status=correctable weight=1 ! 8:HN self-cal status=correctable weight=1 ! 8:KK self-cal status=correctable weight=1 ! 8:KP self-cal status=correctable weight=1 ! 8:LA self-cal status=correctable weight=1 ! 8:MK self-cal status=correctable weight=1 ! 8:NY self-cal status=correctable weight=1 ! 8:OV self-cal status=correctable weight=1 ! 8:PT self-cal status=correctable weight=1 ! 8:SC self-cal status=correctable weight=1 ! 8:WZ self-cal status=correctable weight=1 ! 8:ZC self-cal status=correctable weight=1 system "rm /tmp/ff.ff" mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.200x0.200 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 65 telescope corrections were flagged in sub-array 1. ! A total of 70 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 6 telescope corrections were flagged in sub-array 5. ! A total of 193 telescope corrections were flagged in sub-array 6. ! A total of 24 telescope corrections were flagged in sub-array 7. ! A total of 229 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 2. ! A total of 65 telescope corrections were flagged in sub-array 1. ! A total of 70 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 6 telescope corrections were flagged in sub-array 5. ! A total of 193 telescope corrections were flagged in sub-array 6. ! A total of 24 telescope corrections were flagged in sub-array 7. ! A total of 229 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 3. ! A total of 65 telescope corrections were flagged in sub-array 1. ! A total of 70 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 6 telescope corrections were flagged in sub-array 5. ! A total of 193 telescope corrections were flagged in sub-array 6. ! A total of 24 telescope corrections were flagged in sub-array 7. ! A total of 229 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 4. ! A total of 65 telescope corrections were flagged in sub-array 1. ! A total of 70 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 6 telescope corrections were flagged in sub-array 5. ! A total of 193 telescope corrections were flagged in sub-array 6. ! A total of 24 telescope corrections were flagged in sub-array 7. ! A total of 229 telescope corrections were flagged in sub-array 8. ! ! Fit before self-cal, rms=1.293819Jy sigma=6.416723 ! Fit after self-cal, rms=0.326956Jy sigma=1.644389 ! 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:ZC ! Flagging baselines of antenna 1:ZC in the currently selected channels. flag 2:KK ! Flagging baselines of antenna 2:KK 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:KK ! Flagging baselines of antenna 3:KK in the currently selected channels. flag 5:WZ ! Flagging baselines of antenna 5:WZ in the currently selected channels. flag 5:ZC ! Flagging baselines of antenna 5:ZC in the currently selected channels. flag 6:KK ! Flagging baselines of antenna 6:KK in the currently selected channels. flag 6:NY ! Flagging baselines of antenna 6:NY in the currently selected channels. flag 6:WZ ! Flagging baselines of antenna 6:WZ in the currently selected channels. flag 6:ZC ! Flagging baselines of antenna 6:ZC in the currently selected channels. flag 7:KK ! Flagging baselines of antenna 7:KK in the currently selected channels. flag 7:ZC ! Flagging baselines of antenna 7:ZC in the currently selected channels. flag 8:KK ! Flagging baselines of antenna 8:KK 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. ![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.775 mas, bmaj=1.12 mas, bpa=-43.96 degrees ! Estimated noise=4.50839 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.580497 Jy ! Component: 100 - total flux cleaned = 0.710275 Jy ! Total flux subtracted in 100 components = 0.710275 Jy ! Clean residual min=-0.047270 max=0.121212 Jy/beam ! Clean residual mean=0.000263 rms=0.012090 Jy/beam ! Combined flux in latest and established models = 0.710275 Jy ! Performing phase self-cal ! Adding 4 model components to the UV plane model. ! The established model now contains 4 components and 0.710275 Jy ! ! Correcting IF 1. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 74 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 80 telescope corrections were flagged in sub-array 6. ! A total of 33 telescope corrections were flagged in sub-array 7. ! A total of 95 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 2. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 74 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 80 telescope corrections were flagged in sub-array 6. ! A total of 33 telescope corrections were flagged in sub-array 7. ! A total of 95 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 3. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 74 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 80 telescope corrections were flagged in sub-array 6. ! A total of 33 telescope corrections were flagged in sub-array 7. ! A total of 95 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 4. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 74 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 80 telescope corrections were flagged in sub-array 6. ! A total of 33 telescope corrections were flagged in sub-array 7. ! A total of 95 telescope corrections were flagged in sub-array 8. ! ! Fit before self-cal, rms=0.279239Jy sigma=1.122478 ! Fit after self-cal, rms=0.278909Jy sigma=1.119142 ! Inverting map and beam ! Estimated beam: bmin=0.7879 mas, bmaj=1.117 mas, bpa=-44.16 degrees ! Estimated noise=4.52921 mJy/beam. ! Added new window around map position (-1.4, -1). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.100949 Jy ! Component: 100 - total flux cleaned = 0.13626 Jy ! Total flux subtracted in 100 components = 0.13626 Jy ! Clean residual min=-0.042294 max=0.111987 Jy/beam ! Clean residual mean=0.000117 rms=0.010412 Jy/beam ! Combined flux in latest and established models = 0.846535 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 12 components and 0.846535 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.234434Jy sigma=0.939992 ! Fit after self-cal, rms=0.227720Jy sigma=0.913452 ! Inverting map ! Added new window around map position (1.4, 1). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.074629 Jy ! Component: 100 - total flux cleaned = 0.0987733 Jy ! Total flux subtracted in 100 components = 0.0987733 Jy ! Clean residual min=-0.039185 max=0.043802 Jy/beam ! Clean residual mean=0.000005 rms=0.008809 Jy/beam ! Combined flux in latest and established models = 0.945308 Jy ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 16 components and 0.945308 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.204273Jy sigma=0.820200 ! Fit after self-cal, rms=0.203767Jy sigma=0.818702 ! Inverting map print "=== Finished uniform weighting clean ===" ! === Finished uniform weighting clean === uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.9571 mas, bmaj=1.351 mas, bpa=-42.7 degrees ! Estimated noise=3.22472 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, -1.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0507016 Jy ! Component: 100 - total flux cleaned = 0.0673713 Jy ! Component: 150 - total flux cleaned = 0.0683989 Jy ! Component: 200 - total flux cleaned = 0.0683967 Jy ! Total flux subtracted in 200 components = 0.0683967 Jy ! Clean residual min=-0.034605 max=0.041390 Jy/beam ! Clean residual mean=-0.000164 rms=0.007704 Jy/beam ! Combined flux in latest and established models = 1.01371 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 36 components and 1.01371 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.189387Jy sigma=0.752674 ! Fit after self-cal, rms=0.188681Jy sigma=0.749183 ! 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: ! BR 0.93 FD 0.95 HN 1.33 KP 0.94 ! LA 1.05 MK 0.98 NL 1.23 OV 0.96 ! PT 0.99 ZC 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.00* HN 0.00* KK 0.00* NY 0.00* ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* KK 0.00* MK 0.00* ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! BR 0.83 FD 0.94 HN 1.88 KK 1.00* ! KP 0.97 LA 1.04 MK 0.94 NL 1.11 ! NY 1.00* OV 0.91 PT 0.97 SC 0.92 ! WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 7: ! KK 0.00* MK 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 8: ! BR 0.98 FD 0.93 HN 1.97 KK 1.00* ! KP 1.10 LA 1.08 MK 0.89 NY 1.00* ! OV 0.91 PT 0.98 SC 0.95 WZ 1.00* ! ZC 1.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.90 FD 0.96 HN 1.43 KP 0.89 ! LA 1.03 MK 1.00 NL 1.22 OV 0.92 ! PT 1.04 ZC 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.00* HN 0.00* KK 0.00* NY 0.00* ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* KK 0.00* MK 0.00* ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! BR 0.85 FD 0.93 HN 1.72 KK 1.00* ! KP 0.93 LA 1.03 MK 0.95 NL 1.07 ! NY 1.00* OV 0.88 PT 1.02 SC 0.95 ! WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 7: ! KK 0.00* MK 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 8: ! BR 0.97 FD 0.96 HN 1.95 KK 1.00* ! KP 1.05 LA 1.07 MK 0.88 NY 1.00* ! OV 0.88 PT 1.04 SC 0.97 WZ 1.00* ! ZC 1.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.93 FD 0.97 HN 1.41 KP 0.83 ! LA 1.07 MK 0.98 NL 1.22 OV 0.96 ! PT 1.03 ZC 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.00* HN 0.00* KK 0.00* NY 0.00* ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* KK 0.00* MK 0.00* ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! BR 0.88 FD 0.94 HN 1.73 KK 1.00* ! KP 0.84 LA 1.06 MK 0.94 NL 1.16 ! NY 1.00* OV 0.91 PT 1.02 SC 1.01 ! WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 7: ! KK 0.00* MK 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 8: ! BR 1.06 FD 0.93 HN 2.00 KK 1.00* ! KP 0.93 LA 1.12 MK 0.90 NY 1.00* ! OV 0.91 PT 1.05 SC 1.00 WZ 1.00* ! ZC 1.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 0.91 FD 1.00 HN 1.44 KP 0.88 ! LA 1.06 MK 1.01 NL 1.20 OV 0.96 ! PT 1.01 ZC 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.00* HN 0.00* KK 0.00* NY 0.00* ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* KK 0.00* MK 0.00* ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! BR 0.83 FD 0.94 HN 1.95 KK 1.00* ! KP 0.90 LA 1.08 MK 0.94 NL 1.15 ! NY 1.00* OV 0.88 PT 1.02 SC 1.03 ! WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 7: ! KK 0.00* MK 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 8: ! BR 0.99 FD 0.97 HN 2.26 KK 1.00* ! KP 1.01 LA 1.13 MK 0.94 NY 1.00* ! OV 0.91 PT 1.02 SC 1.02 WZ 1.00* ! ZC 1.00* ! ! ! Fit before self-cal, rms=0.188681Jy sigma=0.749183 ! Fit after self-cal, rms=0.133412Jy sigma=0.457903 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.9572 mas, bmaj=1.374 mas, bpa=-42.31 degrees ! Estimated noise=3.31462 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.0211925 Jy ! Component: 100 - total flux cleaned = -0.0250209 Jy ! Component: 150 - total flux cleaned = -0.0275171 Jy ! Component: 200 - total flux cleaned = -0.0296511 Jy ! Total flux subtracted in 200 components = -0.0296511 Jy ! Clean residual min=-0.018899 max=0.011977 Jy/beam ! Clean residual mean=-0.000010 rms=0.003092 Jy/beam ! Combined flux in latest and established models = 0.984054 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 48 components and 0.984054 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.126037Jy sigma=0.416133 ! Fit after self-cal, rms=0.125967Jy sigma=0.415196 ! 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.125967Jy sigma=0.415196 ! Fit after self-cal, rms=0.126290Jy sigma=0.409433 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.9571 mas, bmaj=1.38 mas, bpa=-42.2 degrees ! Estimated noise=3.30828 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0144804 Jy ! Component: 100 - total flux cleaned = -0.0170598 Jy ! Component: 150 - total flux cleaned = -0.0192021 Jy ! Component: 200 - total flux cleaned = -0.0210571 Jy ! Total flux subtracted in 200 components = -0.0210571 Jy ! Clean residual min=-0.014701 max=0.010103 Jy/beam ! Clean residual mean=-0.000005 rms=0.002529 Jy/beam ! Combined flux in latest and established models = 0.962997 Jy selfcal ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 48 components and 0.962997 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.122500Jy sigma=0.389298 ! Fit after self-cal, rms=0.122266Jy sigma=0.388912 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.00199539 Jy ! Component: 100 - total flux cleaned = -0.0031702 Jy ! Component: 150 - total flux cleaned = -0.00394478 Jy ! Component: 200 - total flux cleaned = -0.00462866 Jy ! Total flux subtracted in 200 components = -0.00462866 Jy ! Clean residual min=-0.013211 max=0.009705 Jy/beam ! Clean residual mean=0.000006 rms=0.002418 Jy/beam ! Combined flux in latest and established models = 0.958368 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 53 components and 0.958368 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.121140Jy sigma=0.383212 ! Fit after self-cal, rms=0.120962Jy sigma=0.382899 ! 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.120962Jy sigma=0.382899 ! Fit after self-cal, rms=0.121531Jy sigma=0.377483 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.9565 mas, bmaj=1.384 mas, bpa=-42.13 degrees ! Estimated noise=3.31018 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00992476 Jy ! Component: 100 - total flux cleaned = -0.0115278 Jy ! Component: 150 - total flux cleaned = -0.0131792 Jy ! Component: 200 - total flux cleaned = -0.0143962 Jy ! Total flux subtracted in 200 components = -0.0143962 Jy ! Clean residual min=-0.011012 max=0.008504 Jy/beam ! Clean residual mean=-0.000004 rms=0.002047 Jy/beam ! Combined flux in latest and established models = 0.943972 Jy selfcal ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 57 components and 0.943972 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.119718Jy sigma=0.369106 ! Fit after self-cal, rms=0.119597Jy sigma=0.368887 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.00161454 Jy ! Component: 100 - total flux cleaned = -0.00201268 Jy ! Component: 150 - total flux cleaned = -0.0023803 Jy ! Component: 200 - total flux cleaned = -0.00238226 Jy ! Total flux subtracted in 200 components = -0.00238226 Jy ! Clean residual min=-0.010221 max=0.008307 Jy/beam ! Clean residual mean=0.000002 rms=0.001990 Jy/beam ! Combined flux in latest and established models = 0.94159 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 63 components and 0.94159 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.118956Jy sigma=0.366235 ! Fit after self-cal, rms=0.118858Jy sigma=0.366068 ! 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.118858Jy sigma=0.366068 ! Fit after self-cal, rms=0.119441Jy sigma=0.363421 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.668 mas, bmaj=3.1 mas, bpa=-33.43 degrees ! Estimated noise=4.14189 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00741896 Jy ! Component: 100 - total flux cleaned = -0.0108242 Jy ! Component: 150 - total flux cleaned = -0.0114733 Jy ! Component: 200 - total flux cleaned = -0.0116675 Jy ! Total flux subtracted in 200 components = -0.0116675 Jy ! Clean residual min=-0.011170 max=0.007840 Jy/beam ! Clean residual mean=-0.000006 rms=0.002270 Jy/beam ! Combined flux in latest and established models = 0.929922 Jy selfcal ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 81 components and 0.929922 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.118994Jy sigma=0.360958 ! Fit after self-cal, rms=0.118931Jy sigma=0.360798 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.000620618 Jy ! Component: 100 - total flux cleaned = -0.000620365 Jy ! Component: 150 - total flux cleaned = -0.000697321 Jy ! Component: 200 - total flux cleaned = -0.000695984 Jy ! Total flux subtracted in 200 components = -0.000695984 Jy ! Clean residual min=-0.011008 max=0.007605 Jy/beam ! Clean residual mean=-0.000002 rms=0.002259 Jy/beam ! Combined flux in latest and established models = 0.929226 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 82 components and 0.929226 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.118778Jy sigma=0.360059 ! Fit after self-cal, rms=0.118746Jy sigma=0.360029 ! 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. ! Applying 2560 buffered edits. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. integer lun lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_selfant_process.csh /tmp/ff.ff" ![@/tmp/ff.ff] selfant 1:BR,true ! 1:BR self-cal status=fixed weight=1 selfant 1:FD,true ! 1:FD self-cal status=fixed weight=1 selfant 1:HN,true ! 1:HN self-cal status=fixed weight=1 selfant 1:KP,true ! 1:KP self-cal status=fixed weight=1 selfant 1:LA,true ! 1:LA self-cal status=fixed weight=1 selfant 1:MK,true ! 1:MK self-cal status=fixed weight=1 selfant 1:NL,true ! 1:NL self-cal status=fixed weight=1 selfant 1:OV,true ! 1:OV self-cal status=fixed weight=1 selfant 1:PT,true ! 1:PT self-cal status=fixed weight=1 selfant 2:BR,true ! 2:BR self-cal status=fixed weight=1 selfant 2:HN,true ! 2:HN self-cal status=fixed weight=1 selfant 3:BR,true ! 3:BR self-cal status=fixed weight=1 selfant 3:MK,true ! 3:MK self-cal status=fixed weight=1 selfant 4:BR,true ! 4:BR self-cal status=fixed weight=1 selfant 4:MK,true ! 4:MK self-cal status=fixed weight=1 selfant 4:OV,true ! 4:OV self-cal status=fixed weight=1 selfant 6:BR,true ! 6:BR self-cal status=fixed weight=1 selfant 6:FD,true ! 6:FD self-cal status=fixed weight=1 selfant 6:HN,true ! 6:HN 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:NL,true ! 6:NL 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 selfant 6:SC,true ! 6:SC self-cal status=fixed weight=1 selfant 7:MK,true ! 7:MK self-cal status=fixed weight=1 selfant 8:BR,true ! 8:BR self-cal status=fixed weight=1 selfant 8:FD,true ! 8:FD self-cal status=fixed weight=1 selfant 8:HN,true ! 8:HN self-cal status=fixed weight=1 selfant 8:KP,true ! 8:KP self-cal status=fixed weight=1 selfant 8:LA,true ! 8:LA self-cal status=fixed weight=1 selfant 8:MK,true ! 8:MK self-cal status=fixed weight=1 selfant 8:OV,true ! 8:OV self-cal status=fixed weight=1 selfant 8:PT,true ! 8:PT self-cal status=fixed weight=1 selfant 8:SC,true ! 8:SC self-cal status=fixed weight=1 ![Exited script file: /tmp/ff.ff] selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 32 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 5. ! A total of 44 telescope corrections were flagged in sub-array 6. ! A total of 38 telescope corrections were flagged in sub-array 7. ! A total of 44 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 2. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 32 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 5. ! A total of 44 telescope corrections were flagged in sub-array 6. ! A total of 38 telescope corrections were flagged in sub-array 7. ! A total of 44 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 3. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 32 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 5. ! A total of 44 telescope corrections were flagged in sub-array 6. ! A total of 38 telescope corrections were flagged in sub-array 7. ! A total of 44 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 4. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 32 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 5. ! A total of 44 telescope corrections were flagged in sub-array 6. ! A total of 38 telescope corrections were flagged in sub-array 7. ! A total of 44 telescope corrections were flagged in sub-array 8. ! ! Fit before self-cal, rms=0.216038Jy sigma=1.224936 ! Fit after self-cal, rms=0.162807Jy sigma=0.665478 selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:KP self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:MK self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:OV self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 1:ZC self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:HN self-cal status=correctable weight=1 ! 2:KK self-cal status=correctable weight=1 ! 2:NY 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:KK self-cal status=correctable weight=1 ! 3:MK self-cal status=correctable weight=1 ! 4:BR self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 5:WZ self-cal status=correctable weight=1 ! 5:ZC self-cal status=correctable weight=1 ! 6:BR self-cal status=correctable weight=1 ! 6:FD self-cal status=correctable weight=1 ! 6:HN 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:NL self-cal status=correctable weight=1 ! 6:NY self-cal status=correctable weight=1 ! 6:OV self-cal status=correctable weight=1 ! 6:PT self-cal status=correctable weight=1 ! 6:SC self-cal status=correctable weight=1 ! 6:WZ self-cal status=correctable weight=1 ! 6:ZC self-cal status=correctable weight=1 ! 7:KK self-cal status=correctable weight=1 ! 7:MK self-cal status=correctable weight=1 ! 7:ZC self-cal status=correctable weight=1 ! 8:BR self-cal status=correctable weight=1 ! 8:FD self-cal status=correctable weight=1 ! 8:HN self-cal status=correctable weight=1 ! 8:KK self-cal status=correctable weight=1 ! 8:KP self-cal status=correctable weight=1 ! 8:LA self-cal status=correctable weight=1 ! 8:MK self-cal status=correctable weight=1 ! 8:NY self-cal status=correctable weight=1 ! 8:OV self-cal status=correctable weight=1 ! 8:PT self-cal status=correctable weight=1 ! 8:SC self-cal status=correctable weight=1 ! 8:WZ self-cal status=correctable weight=1 ! 8:ZC self-cal status=correctable weight=1 selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 58 telescope corrections were flagged in sub-array 1. ! A total of 38 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 9 telescope corrections were flagged in sub-array 4. ! A total of 161 telescope corrections were flagged in sub-array 6. ! A total of 19 telescope corrections were flagged in sub-array 7. ! A total of 185 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 2. ! A total of 58 telescope corrections were flagged in sub-array 1. ! A total of 38 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 9 telescope corrections were flagged in sub-array 4. ! A total of 161 telescope corrections were flagged in sub-array 6. ! A total of 19 telescope corrections were flagged in sub-array 7. ! A total of 185 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 3. ! A total of 58 telescope corrections were flagged in sub-array 1. ! A total of 38 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 9 telescope corrections were flagged in sub-array 4. ! A total of 161 telescope corrections were flagged in sub-array 6. ! A total of 19 telescope corrections were flagged in sub-array 7. ! A total of 185 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 4. ! A total of 58 telescope corrections were flagged in sub-array 1. ! A total of 38 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 9 telescope corrections were flagged in sub-array 4. ! A total of 161 telescope corrections were flagged in sub-array 6. ! A total of 19 telescope corrections were flagged in sub-array 7. ! A total of 185 telescope corrections were flagged in sub-array 8. ! ! Fit before self-cal, rms=0.162807Jy sigma=0.665478 ! Fit after self-cal, rms=0.139630Jy sigma=0.565816 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 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 milli-arcsec ! Estimated beam: bmin=0.4861 mas, bmaj=0.6161 mas, bpa=-71.89 degrees ! Estimated noise=3.07929 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.0232691 Jy ! Component: 100 - total flux cleaned = -0.0265881 Jy ! Component: 150 - total flux cleaned = -0.029155 Jy ! Component: 200 - total flux cleaned = -0.0313703 Jy ! Total flux subtracted in 200 components = -0.0313703 Jy ! Clean residual min=-0.013012 max=0.012471 Jy/beam ! Clean residual mean=0.000007 rms=0.002213 Jy/beam ! Combined flux in latest and established models = 0.897856 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 94 components and 0.897856 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.130258Jy sigma=0.505344 ! Fit after self-cal, rms=0.129026Jy sigma=0.498466 ! Inverting map ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 milli-arcsec ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00196133 Jy ! Component: 100 - total flux cleaned = 0.00370324 Jy ! Component: 150 - total flux cleaned = 0.00521162 Jy ! Component: 200 - total flux cleaned = 0.00748591 Jy ! Total flux subtracted in 200 components = 0.00748591 Jy ! Clean residual min=-0.011129 max=0.010513 Jy/beam ! Clean residual mean=0.000003 rms=0.001952 Jy/beam ! Combined flux in latest and established models = 0.905342 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 104 components and 0.905342 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.126341Jy sigma=0.481572 ! Fit after self-cal, rms=0.126233Jy sigma=0.480546 ! Inverting map ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 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 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 milli-arcsec ! Estimated beam: bmin=0.5859 mas, bmaj=0.7565 mas, bpa=-69.82 degrees ! Estimated noise=2.08859 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 (-2, 4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0012675 Jy ! Component: 100 - total flux cleaned = 0.00088097 Jy ! Component: 150 - total flux cleaned = 0.000886624 Jy ! Component: 200 - total flux cleaned = 0.000879238 Jy ! Component: 250 - total flux cleaned = 0.00117556 Jy ! Component: 300 - total flux cleaned = 0.00144675 Jy ! Component: 350 - total flux cleaned = 0.00144337 Jy ! Component: 400 - total flux cleaned = 0.0016829 Jy ! Total flux subtracted in 400 components = 0.0016829 Jy ! Clean residual min=-0.009037 max=0.008402 Jy/beam ! Clean residual mean=0.000008 rms=0.001775 Jy/beam ! Combined flux in latest and established models = 0.907025 Jy ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 124 components and 0.907025 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.123881Jy sigma=0.463456 ! Fit after self-cal, rms=0.123936Jy sigma=0.462972 ! Inverting map ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 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.000304346 Jy ! Component: 100 - total flux cleaned = 0.000768113 Jy ! Component: 150 - total flux cleaned = 0.000982947 Jy ! Component: 200 - total flux cleaned = 0.000983567 Jy ! Component: 250 - total flux cleaned = 0.00117887 Jy ! Component: 300 - total flux cleaned = 0.00118203 Jy ! Component: 350 - total flux cleaned = 0.00118238 Jy ! Component: 400 - total flux cleaned = 0.00101261 Jy ! Total flux subtracted in 400 components = 0.00101261 Jy ! Clean residual min=-0.008460 max=0.007964 Jy/beam ! Clean residual mean=0.000006 rms=0.001673 Jy/beam ! Combined flux in latest and established models = 0.908038 Jy ! Performing phase self-cal ! Adding 53 model components to the UV plane model. ! The established model now contains 137 components and 0.908037 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.123001Jy sigma=0.457703 ! Fit after self-cal, rms=0.123071Jy sigma=0.457372 ! Inverting map ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 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 65 telescope corrections were flagged in sub-array 1. ! A total of 70 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 9 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! A total of 205 telescope corrections were flagged in sub-array 6. ! A total of 57 telescope corrections were flagged in sub-array 7. ! A total of 229 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 2. ! A total of 65 telescope corrections were flagged in sub-array 1. ! A total of 70 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 9 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! A total of 205 telescope corrections were flagged in sub-array 6. ! A total of 57 telescope corrections were flagged in sub-array 7. ! A total of 229 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 3. ! A total of 65 telescope corrections were flagged in sub-array 1. ! A total of 70 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 9 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! A total of 205 telescope corrections were flagged in sub-array 6. ! A total of 57 telescope corrections were flagged in sub-array 7. ! A total of 229 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 4. ! A total of 65 telescope corrections were flagged in sub-array 1. ! A total of 70 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 9 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! A total of 205 telescope corrections were flagged in sub-array 6. ! A total of 57 telescope corrections were flagged in sub-array 7. ! A total of 229 telescope corrections were flagged in sub-array 8. ! ! Fit before self-cal, rms=0.199964Jy sigma=1.029465 ! Fit after self-cal, rms=0.122668Jy sigma=0.450588 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 milli-arcsec ! Estimated beam: bmin=0.5884 mas, bmaj=0.7567 mas, bpa=-70.23 degrees ! Estimated noise=2.09558 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0089884 Jy ! Component: 100 - total flux cleaned = -0.010667 Jy ! Component: 150 - total flux cleaned = -0.0118078 Jy ! Component: 200 - total flux cleaned = -0.0122125 Jy ! Component: 250 - total flux cleaned = -0.0122123 Jy ! Component: 300 - total flux cleaned = -0.0123763 Jy ! Component: 350 - total flux cleaned = -0.0123738 Jy ! Component: 400 - total flux cleaned = -0.0122258 Jy ! Total flux subtracted in 400 components = -0.0122258 Jy ! Clean residual min=-0.006065 max=0.006408 Jy/beam ! Clean residual mean=0.000011 rms=0.001363 Jy/beam ! Combined flux in latest and established models = 0.895812 Jy selfcal ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 142 components and 0.895812 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.121310Jy sigma=0.442302 ! Fit after self-cal, rms=0.121349Jy sigma=0.441932 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 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 (0.8, -7.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00530386 Jy ! Component: 100 - total flux cleaned = 0.00583076 Jy ! Component: 150 - total flux cleaned = 0.00612755 Jy ! Component: 200 - total flux cleaned = 0.006413 Jy ! Component: 250 - total flux cleaned = 0.00654791 Jy ! Component: 300 - total flux cleaned = 0.00642341 Jy ! Component: 350 - total flux cleaned = 0.00642769 Jy ! Component: 400 - total flux cleaned = 0.00642963 Jy ! Total flux subtracted in 400 components = 0.00642963 Jy ! Clean residual min=-0.005745 max=0.005648 Jy/beam ! Clean residual mean=0.000008 rms=0.001279 Jy/beam ! Combined flux in latest and established models = 0.902241 Jy ! Performing phase self-cal ! Adding 52 model components to the UV plane model. ! The established model now contains 163 components and 0.902241 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.120652Jy sigma=0.438876 ! Fit after self-cal, rms=0.120643Jy sigma=0.438687 ! Inverting map ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 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.120643Jy sigma=0.438687 ! Fit after self-cal, rms=0.120932Jy sigma=0.437453 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 milli-arcsec ! Estimated beam: bmin=0.5892 mas, bmaj=0.7566 mas, bpa=-70.32 degrees ! Estimated noise=2.0969 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00529373 Jy ! Component: 100 - total flux cleaned = -0.00627081 Jy ! Component: 150 - total flux cleaned = -0.00627268 Jy ! Component: 200 - total flux cleaned = -0.00602059 Jy ! Component: 250 - total flux cleaned = -0.00590637 Jy ! Component: 300 - total flux cleaned = -0.00590743 Jy ! Component: 350 - total flux cleaned = -0.00580063 Jy ! Component: 400 - total flux cleaned = -0.00590107 Jy ! Total flux subtracted in 400 components = -0.00590107 Jy ! Clean residual min=-0.004991 max=0.005280 Jy/beam ! Clean residual mean=0.000007 rms=0.001187 Jy/beam ! Combined flux in latest and established models = 0.89634 Jy selfcal ! Performing phase self-cal ! Adding 54 model components to the UV plane model. ! The established model now contains 180 components and 0.89634 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.120384Jy sigma=0.434818 ! Fit after self-cal, rms=0.120385Jy sigma=0.434716 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 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 (19.8, -30). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00434041 Jy ! Component: 100 - total flux cleaned = 0.00530681 Jy ! Component: 150 - total flux cleaned = 0.00574427 Jy ! Component: 200 - total flux cleaned = 0.0060568 Jy ! Component: 250 - total flux cleaned = 0.00635387 Jy ! Component: 300 - total flux cleaned = 0.0064497 Jy ! Component: 350 - total flux cleaned = 0.00645057 Jy ! Component: 400 - total flux cleaned = 0.00654268 Jy ! Total flux subtracted in 400 components = 0.00654268 Jy ! Clean residual min=-0.004929 max=0.004714 Jy/beam ! Clean residual mean=0.000006 rms=0.001131 Jy/beam ! Combined flux in latest and established models = 0.902883 Jy ! Performing phase self-cal ! Adding 60 model components to the UV plane model. ! The established model now contains 200 components and 0.902883 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.119997Jy sigma=0.433109 ! Fit after self-cal, rms=0.119973Jy sigma=0.433015 ! Inverting map ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 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.119973Jy sigma=0.433015 ! Fit after self-cal, rms=0.119161Jy sigma=0.424402 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 milli-arcsec ! Estimated beam: bmin=2.409 mas, bmaj=2.578 mas, bpa=-27.79 degrees ! Estimated noise=3.40475 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00122497 Jy ! Component: 100 - total flux cleaned = 0.00231807 Jy ! Component: 150 - total flux cleaned = 0.00320704 Jy ! Component: 200 - total flux cleaned = 0.00383822 Jy ! Component: 250 - total flux cleaned = 0.00447307 Jy ! Component: 300 - total flux cleaned = 0.00488242 Jy ! Component: 350 - total flux cleaned = 0.00521018 Jy ! Component: 400 - total flux cleaned = 0.00546688 Jy ! Total flux subtracted in 400 components = 0.00546688 Jy ! Clean residual min=-0.007000 max=0.006433 Jy/beam ! Clean residual mean=0.000008 rms=0.002187 Jy/beam ! Combined flux in latest and established models = 0.90835 Jy selfcal ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 220 components and 0.90835 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.119039Jy sigma=0.423694 ! Fit after self-cal, rms=0.119024Jy sigma=0.423640 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 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.000455608 Jy ! Component: 100 - total flux cleaned = 0.000761635 Jy ! Component: 150 - total flux cleaned = 0.000953255 Jy ! Component: 200 - total flux cleaned = 0.00117874 Jy ! Component: 250 - total flux cleaned = 0.00130804 Jy ! Component: 300 - total flux cleaned = 0.00143183 Jy ! Component: 350 - total flux cleaned = 0.00154942 Jy ! Component: 400 - total flux cleaned = 0.00166366 Jy ! Total flux subtracted in 400 components = 0.00166366 Jy ! Clean residual min=-0.006780 max=0.006547 Jy/beam ! Clean residual mean=0.000004 rms=0.002152 Jy/beam ! Combined flux in latest and established models = 0.910013 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 232 components and 0.910013 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.118996Jy sigma=0.423546 ! Fit after self-cal, rms=0.118989Jy sigma=0.423532 ! Inverting map ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 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 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 milli-arcsec ! Estimated beam: bmin=0.5898 mas, bmaj=0.7566 mas, bpa=-70.43 degrees ! Estimated noise=2.09894 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.5898 x 0.7566 at -70.43 degrees (North through East) ! Clean map min=-0.011224 max=0.75411 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.118989Jy sigma=0.423532 ! Fit after self-cal, rms=0.114009Jy sigma=0.389085 ! Inverting map and beam ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 milli-arcsec ! Estimated beam: bmin=0.5898 mas, bmaj=0.757 mas, bpa=-70.29 degrees ! Estimated noise=2.0953 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00276524 Jy ! Component: 100 - total flux cleaned = -0.00398469 Jy ! Component: 150 - total flux cleaned = -0.00438789 Jy ! Component: 200 - total flux cleaned = -0.00476685 Jy ! Component: 250 - total flux cleaned = -0.00494759 Jy ! Component: 300 - total flux cleaned = -0.0051217 Jy ! Component: 350 - total flux cleaned = -0.00512272 Jy ! Component: 400 - total flux cleaned = -0.00512347 Jy ! Total flux subtracted in 400 components = -0.00512347 Jy ! Clean residual min=-0.004403 max=0.004638 Jy/beam ! Clean residual mean=0.000001 rms=0.001060 Jy/beam ! Combined flux in latest and established models = 0.90489 Jy ! Performing phase self-cal ! Adding 53 model components to the UV plane model. ! The established model now contains 246 components and 0.90489 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.113696Jy sigma=0.387549 ! Fit after self-cal, rms=0.113680Jy sigma=0.387492 ! Inverting map ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 milli-arcsec ! Added new window around map position (-1.8, -9.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00257997 Jy ! Component: 100 - total flux cleaned = 0.00288614 Jy ! Component: 150 - total flux cleaned = 0.00307104 Jy ! Component: 200 - total flux cleaned = 0.00298421 Jy ! Component: 250 - total flux cleaned = 0.00290126 Jy ! Component: 300 - total flux cleaned = 0.0029004 Jy ! Component: 350 - total flux cleaned = 0.00290029 Jy ! Component: 400 - total flux cleaned = 0.0029762 Jy ! Total flux subtracted in 400 components = 0.0029762 Jy ! Clean residual min=-0.004265 max=0.004322 Jy/beam ! Clean residual mean=0.000002 rms=0.001013 Jy/beam ! Combined flux in latest and established models = 0.907866 Jy ! Performing phase self-cal ! Adding 55 model components to the UV plane model. ! The established model now contains 271 components and 0.907866 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.113377Jy sigma=0.386192 ! Fit after self-cal, rms=0.113363Jy sigma=0.386111 ! Inverting map ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 milli-arcsec ! Added new window around map position (1.2, 7). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00312021 Jy ! Component: 100 - total flux cleaned = 0.00321909 Jy ! Component: 150 - total flux cleaned = 0.00358706 Jy ! Component: 200 - total flux cleaned = 0.00410217 Jy ! Component: 250 - total flux cleaned = 0.00459078 Jy ! Component: 300 - total flux cleaned = 0.00497947 Jy ! Component: 350 - total flux cleaned = 0.00528183 Jy ! Component: 400 - total flux cleaned = 0.00557495 Jy ! Total flux subtracted in 400 components = 0.00557495 Jy ! Clean residual min=-0.004231 max=0.004273 Jy/beam ! Clean residual mean=0.000000 rms=0.000978 Jy/beam ! Combined flux in latest and established models = 0.913441 Jy ! Performing phase self-cal ! Adding 54 model components to the UV plane model. ! The established model now contains 297 components and 0.913441 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.113117Jy sigma=0.385018 ! Fit after self-cal, rms=0.113109Jy sigma=0.384904 ! Inverting map ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 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.113109Jy sigma=0.384904 ! Fit after self-cal, rms=0.113075Jy sigma=0.384509 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 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 milli-arcsec ! Estimated beam: bmin=0.489 mas, bmaj=0.6159 mas, bpa=-72.76 degrees ! Estimated noise=3.10049 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.59272 Jy ! Component: 100 - total flux cleaned = 0.724566 Jy ! Component: 150 - total flux cleaned = 0.771367 Jy ! Component: 200 - total flux cleaned = 0.798403 Jy ! Component: 250 - total flux cleaned = 0.818394 Jy ! Component: 300 - total flux cleaned = 0.830684 Jy ! Component: 350 - total flux cleaned = 0.841093 Jy ! Component: 400 - total flux cleaned = 0.84959 Jy ! Total flux subtracted in 400 components = 0.84959 Jy ! Clean residual min=-0.009833 max=0.010161 Jy/beam ! Clean residual mean=0.000023 rms=0.001690 Jy/beam ! Combined flux in latest and established models = 0.849591 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 40 model components to the UV plane model. ! The established model now contains 40 components and 0.849591 Jy ! Inverting map and beam ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 milli-arcsec ! Estimated beam: bmin=0.5903 mas, bmaj=0.7567 mas, bpa=-70.37 degrees ! Estimated noise=2.09445 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.0127004 Jy ! Component: 100 - total flux cleaned = 0.0202709 Jy ! Component: 150 - total flux cleaned = 0.0257907 Jy ! Component: 200 - total flux cleaned = 0.0299011 Jy ! Component: 250 - total flux cleaned = 0.0328728 Jy ! Component: 300 - total flux cleaned = 0.0353741 Jy ! Component: 350 - total flux cleaned = 0.0379424 Jy ! Component: 400 - total flux cleaned = 0.0403474 Jy ! Total flux subtracted in 400 components = 0.0403474 Jy ! Clean residual min=-0.006055 max=0.005922 Jy/beam ! Clean residual mean=0.000012 rms=0.001297 Jy/beam ! Combined flux in latest and established models = 0.889938 Jy ! Adding 56 model components to the UV plane model. ! The established model now contains 74 components and 0.889938 Jy ! Inverting map ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 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 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 milli-arcsec ! Estimated beam: bmin=2.414 mas, bmaj=2.584 mas, bpa=-27.7 degrees ! Estimated noise=3.38677 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00781891 Jy ! Component: 100 - total flux cleaned = 0.0135642 Jy ! Component: 150 - total flux cleaned = 0.0176374 Jy ! Component: 200 - total flux cleaned = 0.0204442 Jy ! Component: 250 - total flux cleaned = 0.0223416 Jy ! Component: 300 - total flux cleaned = 0.0235453 Jy ! Component: 350 - total flux cleaned = 0.0247321 Jy ! Component: 400 - total flux cleaned = 0.0256113 Jy ! Total flux subtracted in 400 components = 0.0256113 Jy ! Clean residual min=-0.006871 max=0.006980 Jy/beam ! Clean residual mean=-0.000015 rms=0.002118 Jy/beam ! Combined flux in latest and established models = 0.915549 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) ! Added new window around map position (9.4, 40.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00216303 Jy ! Component: 100 - total flux cleaned = 0.00215742 Jy ! Component: 150 - total flux cleaned = 0.00255186 Jy ! Component: 200 - total flux cleaned = 0.00289576 Jy ! Component: 250 - total flux cleaned = 0.00309795 Jy ! Component: 300 - total flux cleaned = 0.00301046 Jy ! Component: 350 - total flux cleaned = 0.00318205 Jy ! Component: 400 - total flux cleaned = 0.00310393 Jy ! Total flux subtracted in 400 components = 0.00310393 Jy ! Clean residual min=-0.005635 max=0.005096 Jy/beam ! Clean residual mean=-0.000035 rms=0.001622 Jy/beam ! Combined flux in latest and established models = 0.918653 Jy ! Adding 105 model components to the UV plane model. ! The established model now contains 174 components and 0.918653 Jy ! Inverting map ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 milli-arcsec uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 milli-arcsec ! Estimated beam: bmin=0.5903 mas, bmaj=0.7567 mas, bpa=-70.37 degrees ! Estimated noise=2.09445 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00149615 Jy ! Component: 100 - total flux cleaned = -0.00205268 Jy ! Component: 150 - total flux cleaned = -0.00205547 Jy ! Component: 200 - total flux cleaned = -0.00189422 Jy ! Component: 250 - total flux cleaned = -0.00205296 Jy ! Component: 300 - total flux cleaned = -0.00219943 Jy ! Component: 350 - total flux cleaned = -0.00219683 Jy ! Component: 400 - total flux cleaned = -0.00205653 Jy ! Total flux subtracted in 400 components = -0.00205653 Jy ! Clean residual min=-0.004725 max=0.004543 Jy/beam ! Clean residual mean=-0.000005 rms=0.001005 Jy/beam ! Combined flux in latest and established models = 0.916597 Jy selfcal ! Performing phase self-cal ! Adding 67 model components to the UV plane model. ! The established model now contains 221 components and 0.916597 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.114817Jy sigma=0.405038 ! Fit after self-cal, rms=0.114458Jy sigma=0.401209 wmodel /scr/pima/rdv73_uvs/J1058+8114_X_map.mod ! Writing 221 model components to file: /scr/pima/rdv73_uvs/J1058+8114_X_map.mod wobs /scr/pima/rdv73_uvs/J1058+8114_X_uvs.fits ! Applying 1663 buffered edits. ! Writing UV FITS file: /scr/pima/rdv73_uvs/J1058+8114_X_uvs.fits wwins /scr/pima/rdv73_uvs/J1058+8114_X_map.win ! wwins: Wrote 10 windows to /scr/pima/rdv73_uvs/J1058+8114_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 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 milli-arcsec ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.000723466 Jy ! Component: 100 - total flux cleaned = -0.00128295 Jy ! Component: 150 - total flux cleaned = -0.00159453 Jy ! Component: 200 - total flux cleaned = -0.00181917 Jy ! Component: 250 - total flux cleaned = -0.0018913 Jy ! Component: 300 - total flux cleaned = -0.00182114 Jy ! Component: 350 - total flux cleaned = -0.00230993 Jy ! Component: 400 - total flux cleaned = -0.0025146 Jy ! Component: 450 - total flux cleaned = -0.00284946 Jy ! Component: 500 - total flux cleaned = -0.00291475 Jy ! Component: 550 - total flux cleaned = -0.00304455 Jy ! Component: 600 - total flux cleaned = -0.00336539 Jy ! Component: 650 - total flux cleaned = -0.00342846 Jy ! Component: 700 - total flux cleaned = -0.0036772 Jy ! Component: 750 - total flux cleaned = -0.00361583 Jy ! Component: 800 - total flux cleaned = -0.00385922 Jy ! Component: 850 - total flux cleaned = -0.00403713 Jy ! Component: 900 - total flux cleaned = -0.00385972 Jy ! Component: 950 - total flux cleaned = -0.00415151 Jy ! Component: 1000 - total flux cleaned = -0.0042677 Jy ! Component: 1050 - total flux cleaned = -0.00426732 Jy ! Component: 1100 - total flux cleaned = -0.0042676 Jy ! Component: 1150 - total flux cleaned = -0.00438022 Jy ! Component: 1200 - total flux cleaned = -0.00454764 Jy ! Component: 1250 - total flux cleaned = -0.00449294 Jy ! Component: 1300 - total flux cleaned = -0.0045467 Jy ! Component: 1350 - total flux cleaned = -0.00449238 Jy ! Component: 1400 - total flux cleaned = -0.00460015 Jy ! Component: 1450 - total flux cleaned = -0.00476011 Jy ! Component: 1500 - total flux cleaned = -0.00486565 Jy ! Component: 1550 - total flux cleaned = -0.00486543 Jy ! Component: 1600 - total flux cleaned = -0.00496962 Jy ! Component: 1650 - total flux cleaned = -0.00491738 Jy ! Component: 1700 - total flux cleaned = -0.00517498 Jy ! Component: 1750 - total flux cleaned = -0.00532777 Jy ! Component: 1800 - total flux cleaned = -0.00547957 Jy ! Component: 1850 - total flux cleaned = -0.00537842 Jy ! Component: 1900 - total flux cleaned = -0.00557852 Jy ! Component: 1950 - total flux cleaned = -0.00562793 Jy ! Component: 2000 - total flux cleaned = -0.00567707 Jy ! Component: 2050 - total flux cleaned = -0.00587346 Jy ! Component: 2100 - total flux cleaned = -0.00592242 Jy ! Component: 2150 - total flux cleaned = -0.00601945 Jy ! Component: 2200 - total flux cleaned = -0.00606782 Jy ! Component: 2250 - total flux cleaned = -0.00611559 Jy ! Component: 2300 - total flux cleaned = -0.00621007 Jy ! Component: 2350 - total flux cleaned = -0.0061629 Jy ! Component: 2400 - total flux cleaned = -0.00635066 Jy ! Component: 2450 - total flux cleaned = -0.00658426 Jy ! Component: 2500 - total flux cleaned = -0.00644499 Jy ! Component: 2550 - total flux cleaned = -0.00667592 Jy ! Component: 2600 - total flux cleaned = -0.00676769 Jy ! Component: 2650 - total flux cleaned = -0.00685895 Jy ! Component: 2700 - total flux cleaned = -0.0067677 Jy ! Component: 2750 - total flux cleaned = -0.00685834 Jy ! Component: 2800 - total flux cleaned = -0.00712845 Jy ! Component: 2850 - total flux cleaned = -0.00717341 Jy ! Component: 2900 - total flux cleaned = -0.00726249 Jy ! Component: 2950 - total flux cleaned = -0.00730691 Jy ! Component: 3000 - total flux cleaned = -0.00739504 Jy ! Component: 3050 - total flux cleaned = -0.00770248 Jy ! Component: 3100 - total flux cleaned = -0.00761449 Jy ! Component: 3150 - total flux cleaned = -0.00765818 Jy ! Component: 3200 - total flux cleaned = -0.00774461 Jy ! Component: 3250 - total flux cleaned = -0.00778796 Jy ! Component: 3300 - total flux cleaned = -0.0077018 Jy ! Component: 3350 - total flux cleaned = -0.0078303 Jy ! Component: 3400 - total flux cleaned = -0.00778733 Jy ! Component: 3450 - total flux cleaned = -0.0078723 Jy ! Component: 3500 - total flux cleaned = -0.00812604 Jy ! Component: 3550 - total flux cleaned = -0.00820979 Jy ! Component: 3600 - total flux cleaned = -0.00795843 Jy ! Component: 3650 - total flux cleaned = -0.00795857 Jy ! Component: 3700 - total flux cleaned = -0.00816655 Jy ! Component: 3750 - total flux cleaned = -0.00812477 Jy ! Component: 3800 - total flux cleaned = -0.00820695 Jy ! Component: 3850 - total flux cleaned = -0.00828889 Jy ! Component: 3900 - total flux cleaned = -0.00841193 Jy ! Component: 3950 - total flux cleaned = -0.00837103 Jy ! Component: 4000 - total flux cleaned = -0.00845236 Jy ! Component: 4050 - total flux cleaned = -0.00853361 Jy ! Total flux subtracted in 4096 components = -0.00845271 Jy ! Clean residual min=-0.002130 max=0.002284 Jy/beam ! Clean residual mean=-0.000001 rms=0.000650 Jy/beam ! Combined flux in latest and established models = 0.908144 Jy keep ! Adding 363 model components to the UV plane model. ! The established model now contains 584 components and 0.908145 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv73_uvs/J1058+8114_X\_map.fits" ! Inverting map ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 milli-arcsec ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.5903 x 0.7567 at -70.37 degrees (North through East) ! Clean map min=-0.0063459 max=0.75624 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv73_uvs/J1058+8114_X_map.fits float peak_flux peak_flux = peak(flux,max) invert ! Inverting map ! Your choice of large map pixels excluded 8.31% of the data. ! The x-axis pixel size should ideally be below 0.1682 milli-arcsec ! The y-axis pixel size should ideally be below 0.1746 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.7562 rms= 0.00064 Jy quit ! Quitting program ! Log file /scr/pima/rdv73_uvs/J1058+8114_X_uvs.log closed on Sat Jan 9 17:33:53 2016