! Started logfile: /scr/pima/rdv73_uvs/J1415+1320_S_uvs.log on Sat Jan 9 17:50:20 2016 obs /scr/pima/rdv73_uvs/J1415+1320_S_uva.fits ! Reading UV FITS file: /scr/pima/rdv73_uvs/J1415+1320_S_uva.fits ! AN table 1: 30 integrations on 3 of 3 possible baselines. ! AN table 2: 9 integrations on 1 of 1 possible baselines. ! AN table 3: 32 integrations on 28 of 28 possible baselines. ! AN table 4: 31 integrations on 55 of 55 possible baselines. ! AN table 5: 31 integrations on 45 of 45 possible baselines. ! AN table 6: 9 integrations on 45 of 45 possible baselines. ! AN table 7: 7 integrations on 10 of 10 possible baselines. ! Apparent sampling: 0.854048 visibilities/baseline/integration-bin. ! Found source: J1415+1320 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 2.23299e+09 1.6e+07 1 1.6e+07 ! 02 2 2.26299e+09 1.6e+07 1 1.6e+07 ! 03 3 2.35299e+09 1.6e+07 1 1.6e+07 ! 04 4 2.37299e+09 1.6e+07 1 1.6e+07 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 15612 visibilities. select RR ! Selecting polarization: RR, channels: 1..4 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 print "field_size =",field_size ! field_size = 1024 print "field_cell =",field_cell ! field_cell = 0.4 print "taper_size =",taper_size ! taper_size = 20 ![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv73_uvs/J1415+1320_S] integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 800000 float soltime2 soltime2 = 15 float soltime3 soltime3 = 1 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float old_peak float new_peak float flux_cutoff selfant "",false ! 1:NY self-cal status=correctable weight=1 ! 1:WZ self-cal status=correctable weight=1 ! 1:ZC self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 2:ZC self-cal status=correctable weight=1 ! 3:FD self-cal status=correctable weight=1 ! 3:HN self-cal status=correctable weight=1 ! 3:LA self-cal status=correctable weight=1 ! 3:NY self-cal status=correctable weight=1 ! 3:PT self-cal status=correctable weight=1 ! 3:SC self-cal status=correctable weight=1 ! 3:WZ self-cal status=correctable weight=1 ! 3:ZC self-cal status=correctable weight=1 ! 4:BR self-cal status=correctable weight=1 ! 4:FD self-cal status=correctable weight=1 ! 4:HN self-cal status=correctable weight=1 ! 4:KP self-cal status=correctable weight=1 ! 4:LA self-cal status=correctable weight=1 ! 4:NY self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 4:PT self-cal status=correctable weight=1 ! 4:SC self-cal status=correctable weight=1 ! 4:WZ self-cal status=correctable weight=1 ! 4:ZC self-cal status=correctable weight=1 ! 5:BR self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:HN self-cal status=correctable weight=1 ! 5:KP self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:NY self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 5:SC self-cal status=correctable weight=1 ! 5:WZ 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:OV self-cal status=correctable weight=1 ! 6:PT self-cal status=correctable weight=1 ! 6:SC self-cal status=correctable weight=1 ! 7:BR self-cal status=correctable weight=1 ! 7:KK self-cal status=correctable weight=1 ! 7:KP self-cal status=correctable weight=1 ! 7:MK self-cal status=correctable weight=1 ! 7:OV self-cal status=correctable weight=1 system "rm /tmp/ff.ff" mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.400x0.400 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 18 telescope corrections were flagged in sub-array 2. ! A total of 24 telescope corrections were flagged in sub-array 3. ! A total of 26 telescope corrections were flagged in sub-array 4. ! A total of 39 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! A total of 7 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 2. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 18 telescope corrections were flagged in sub-array 2. ! A total of 24 telescope corrections were flagged in sub-array 3. ! A total of 26 telescope corrections were flagged in sub-array 4. ! A total of 39 telescope corrections were flagged in sub-array 5. ! A total of 25 telescope corrections were flagged in sub-array 6. ! A total of 13 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 3. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 18 telescope corrections were flagged in sub-array 2. ! A total of 24 telescope corrections were flagged in sub-array 3. ! A total of 26 telescope corrections were flagged in sub-array 4. ! A total of 39 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! A total of 7 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 4. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 18 telescope corrections were flagged in sub-array 2. ! A total of 24 telescope corrections were flagged in sub-array 3. ! A total of 56 telescope corrections were flagged in sub-array 4. ! A total of 39 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! A total of 7 telescope corrections were flagged in sub-array 7. ! ! Fit before self-cal, rms=1.070643Jy sigma=4.672295 ! Fit after self-cal, rms=0.759467Jy sigma=3.329419 ! 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:NY ! Flagging baselines of antenna 1:NY in the currently selected channels. flag 1:WZ ! Flagging baselines of antenna 1:WZ in the currently selected channels. flag 1:ZC ! Flagging baselines of antenna 1:ZC in the currently selected channels. flag 2:WZ ! Flagging baselines of antenna 2:WZ in the currently selected channels. flag 2:ZC ! Flagging baselines of antenna 2:ZC in the currently selected channels. flag 3:NY ! Flagging baselines of antenna 3:NY in the currently selected channels. flag 3:WZ ! Flagging baselines of antenna 3:WZ in the currently selected channels. flag 3:ZC ! Flagging baselines of antenna 3:ZC in the currently selected channels. flag 4:NY ! Flagging baselines of antenna 4:NY in the currently selected channels. flag 4:WZ ! Flagging baselines of antenna 4:WZ in the currently selected channels. flag 4:ZC ! Flagging baselines of antenna 4:ZC in the currently selected channels. flag 5:NY ! Flagging baselines of antenna 5:NY in the currently selected channels. flag 5:WZ ! Flagging baselines of antenna 5:WZ in the currently selected channels. flag 6:KK ! Flagging baselines of antenna 6:KK in the currently selected channels. flag 7:KK ! Flagging baselines of antenna 7:KK in the currently selected channels. ![Exited script file: /tmp/ff.ff] uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.477 mas, bmaj=5.52 mas, bpa=2.654 degrees ! Estimated noise=3.32204 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.210623 Jy ! Component: 100 - total flux cleaned = 0.277702 Jy ! Total flux subtracted in 100 components = 0.277702 Jy ! Clean residual min=-0.032010 max=0.060207 Jy/beam ! Clean residual mean=-0.000168 rms=0.009177 Jy/beam ! Combined flux in latest and established models = 0.277702 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 15 components and 0.277702 Jy ! ! Correcting IF 1. ! A total of 78 telescope corrections were flagged in sub-array 1. ! A total of 91 telescope corrections were flagged in sub-array 3. ! A total of 85 telescope corrections were flagged in sub-array 4. ! A total of 63 telescope corrections were flagged in sub-array 5. ! A total of 8 telescope corrections were flagged in sub-array 6. ! A total of 6 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 2. ! A total of 78 telescope corrections were flagged in sub-array 1. ! A total of 91 telescope corrections were flagged in sub-array 3. ! A total of 85 telescope corrections were flagged in sub-array 4. ! A total of 63 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 78 telescope corrections were flagged in sub-array 1. ! A total of 91 telescope corrections were flagged in sub-array 3. ! A total of 85 telescope corrections were flagged in sub-array 4. ! A total of 63 telescope corrections were flagged in sub-array 5. ! A total of 8 telescope corrections were flagged in sub-array 6. ! A total of 6 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 4. ! A total of 78 telescope corrections were flagged in sub-array 1. ! A total of 91 telescope corrections were flagged in sub-array 3. ! A total of 85 telescope corrections were flagged in sub-array 4. ! A total of 63 telescope corrections were flagged in sub-array 5. ! A total of 8 telescope corrections were flagged in sub-array 6. ! A total of 6 telescope corrections were flagged in sub-array 7. ! ! Fit before self-cal, rms=0.304322Jy sigma=1.278109 ! Fit after self-cal, rms=0.303134Jy sigma=1.269900 ! Inverting map and beam ! Estimated beam: bmin=3.477 mas, bmaj=5.52 mas, bpa=2.655 degrees ! Estimated noise=3.32227 mJy/beam. ! Added new window around map position (-6, -2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0544218 Jy ! Component: 100 - total flux cleaned = 0.0858239 Jy ! Total flux subtracted in 100 components = 0.0858239 Jy ! Clean residual min=-0.032494 max=0.058188 Jy/beam ! Clean residual mean=-0.000131 rms=0.008240 Jy/beam ! Combined flux in latest and established models = 0.363526 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 31 components and 0.363526 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.275267Jy sigma=1.139639 ! Fit after self-cal, rms=0.274090Jy sigma=1.129069 ! Inverting map ! Added new window around map position (17.6, 11.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0523942 Jy ! Component: 100 - total flux cleaned = 0.0864355 Jy ! Total flux subtracted in 100 components = 0.0864355 Jy ! Clean residual min=-0.025139 max=0.040337 Jy/beam ! Clean residual mean=-0.000111 rms=0.007290 Jy/beam ! Combined flux in latest and established models = 0.449962 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 51 components and 0.449962 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.246877Jy sigma=1.000626 ! Fit after self-cal, rms=0.244076Jy sigma=0.993892 ! 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=4.954 mas, bmaj=7.147 mas, bpa=16.9 degrees ! Estimated noise=1.86277 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.0921539 Jy ! Component: 100 - total flux cleaned = 0.156301 Jy ! Component: 150 - total flux cleaned = 0.203206 Jy ! Component: 200 - total flux cleaned = 0.238397 Jy ! Total flux subtracted in 200 components = 0.238397 Jy ! Clean residual min=-0.036054 max=0.059464 Jy/beam ! Clean residual mean=-0.000161 rms=0.012494 Jy/beam ! Combined flux in latest and established models = 0.688359 Jy ! Performing phase self-cal ! Adding 57 model components to the UV plane model. ! The established model now contains 98 components and 0.688359 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.190447Jy sigma=0.766734 ! Fit after self-cal, rms=0.184499Jy sigma=0.745326 ! Inverting map ! Added new window around map position (27.2, 8.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0668639 Jy ! Component: 100 - total flux cleaned = 0.103636 Jy ! Component: 150 - total flux cleaned = 0.13009 Jy ! Component: 200 - total flux cleaned = 0.151246 Jy ! Total flux subtracted in 200 components = 0.151246 Jy ! Clean residual min=-0.023637 max=0.041704 Jy/beam ! Clean residual mean=-0.000024 rms=0.007969 Jy/beam ! Combined flux in latest and established models = 0.839605 Jy ! Performing phase self-cal ! Adding 72 model components to the UV plane model. ! The established model now contains 164 components and 0.839606 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.152791Jy sigma=0.613720 ! Fit after self-cal, rms=0.145927Jy sigma=0.584605 ! 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: ! NY 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! WZ 0.00* ZC 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! FD 1.03 HN 1.19 LA 1.05 NY 1.00* ! PT 0.94 SC 0.91 WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 1.22 FD 0.85 HN 1.24 KP 1.00 ! LA 0.95 NY 1.00* OV 0.98 PT 0.98 ! SC 1.02 WZ 1.00* ZC 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.20 FD 0.85 HN 1.16 KP 0.95 ! LA 1.06 NY 1.00* OV 0.97 PT 1.03 ! SC 1.03 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 6: ! BR 0.92 FD 0.94 HN 1.16 KK 1.00* ! KP 1.20 LA 0.96 MK 0.82 OV 1.25 ! PT 1.04 SC 0.99 ! ! A total of 6 telescope corrections were flagged in sub-array 7. ! Telescope amplitude corrections in sub-array 7: ! BR 0.93 KK 1.00* KP 0.97 MK 0.82 ! OV 1.40 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! NY 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! WZ 0.00* ZC 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! FD 1.13 HN 1.22 LA 1.00 NY 1.00* ! PT 0.92 SC 0.90 WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 1.23 FD 0.92 HN 1.22 KP 1.04 ! LA 0.90 NY 1.00* OV 1.05 PT 0.99 ! SC 1.11 WZ 1.00* ZC 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.20 FD 0.92 HN 1.13 KP 0.94 ! LA 0.97 NY 1.00* OV 1.03 PT 1.07 ! SC 1.14 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 6: ! BR 0.95 FD 1.01 HN 1.08 KK 1.00* ! KP 1.25 LA 0.90 MK 0.89 OV 1.32 ! PT 1.03 SC 1.14 ! ! A total of 6 telescope corrections were flagged in sub-array 7. ! Telescope amplitude corrections in sub-array 7: ! BR 1.03 KK 1.00* KP 1.16 MK 0.84 ! OV 1.24 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! NY 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! WZ 0.00* ZC 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! FD 1.20 HN 1.21 LA 0.97 NY 1.00* ! PT 0.90 SC 0.81 WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 0.64 FD 0.94 HN 1.29 KP 1.06 ! LA 0.83 NY 1.00* OV 1.08 PT 1.02 ! SC 0.99 WZ 1.00* ZC 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.23 FD 0.92 HN 1.19 KP 0.97 ! LA 0.88 NY 1.00* OV 1.07 PT 1.09 ! SC 0.94 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 6: ! BR 0.94 FD 1.03 HN 1.07 KK 1.00* ! KP 1.21 LA 0.84 MK 0.89 OV 1.37 ! PT 1.08 SC 0.98 ! ! A total of 6 telescope corrections were flagged in sub-array 7. ! Telescope amplitude corrections in sub-array 7: ! BR 1.14 KK 1.00* KP 0.86 MK 0.89 ! OV 1.39 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! NY 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! WZ 0.00* ZC 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! FD 1.12 HN 1.07 LA 0.80 NY 1.00* ! PT 0.79 SC 0.86 WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 1.00* FD 0.93 HN 1.03 KP 0.94 ! LA 0.72 NY 1.00* OV 0.88 PT 0.89 ! SC 1.02 WZ 1.00* ZC 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.13 FD 0.90 HN 1.00 KP 0.87 ! LA 0.74 NY 1.00* OV 0.89 PT 0.95 ! SC 0.93 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 6: ! BR 0.84 FD 0.97 HN 0.99 KK 1.00* ! KP 1.05 LA 0.74 MK 1.01 OV 1.13 ! PT 0.94 SC 1.09 ! ! A total of 6 telescope corrections were flagged in sub-array 7. ! Telescope amplitude corrections in sub-array 7: ! BR 1.17 KK 1.00* KP 0.69 MK 0.93 ! OV 1.33 ! ! ! Fit before self-cal, rms=0.145927Jy sigma=0.584605 ! Fit after self-cal, rms=0.112422Jy sigma=0.480294 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=5.001 mas, bmaj=7.145 mas, bpa=17.78 degrees ! Estimated noise=1.86837 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 (36, 8.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0230507 Jy ! Component: 100 - total flux cleaned = 0.0270689 Jy ! Component: 150 - total flux cleaned = 0.027874 Jy ! Component: 200 - total flux cleaned = 0.0281771 Jy ! Total flux subtracted in 200 components = 0.0281771 Jy ! Clean residual min=-0.015419 max=0.023584 Jy/beam ! Clean residual mean=0.000022 rms=0.003891 Jy/beam ! Combined flux in latest and established models = 0.867783 Jy ! Performing phase self-cal ! Adding 71 model components to the UV plane model. ! The established model now contains 227 components and 0.867782 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.106440Jy sigma=0.452891 ! Fit after self-cal, rms=0.104914Jy sigma=0.446445 ! 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.104914Jy sigma=0.446445 ! Fit after self-cal, rms=0.104060Jy sigma=0.444740 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=4.997 mas, bmaj=7.142 mas, bpa=17.94 degrees ! Estimated noise=1.86259 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00793599 Jy ! Component: 100 - total flux cleaned = 0.0122349 Jy ! Component: 150 - total flux cleaned = 0.0122715 Jy ! Component: 200 - total flux cleaned = 0.0120275 Jy ! Total flux subtracted in 200 components = 0.0120275 Jy ! Clean residual min=-0.012887 max=0.019210 Jy/beam ! Clean residual mean=-0.000001 rms=0.003258 Jy/beam ! Combined flux in latest and established models = 0.87981 Jy selfcal ! Performing phase self-cal ! Adding 67 model components to the UV plane model. ! The established model now contains 281 components and 0.87981 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.101349Jy sigma=0.432317 ! Fit after self-cal, rms=0.100573Jy sigma=0.429172 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 (-22, -11.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00942962 Jy ! Component: 100 - total flux cleaned = 0.0153016 Jy ! Component: 150 - total flux cleaned = 0.0191128 Jy ! Component: 200 - total flux cleaned = 0.0220949 Jy ! Total flux subtracted in 200 components = 0.0220949 Jy ! Clean residual min=-0.009856 max=0.011426 Jy/beam ! Clean residual mean=0.000001 rms=0.002433 Jy/beam ! Combined flux in latest and established models = 0.901905 Jy ! Performing phase self-cal ! Adding 61 model components to the UV plane model. ! The established model now contains 329 components and 0.901904 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.097261Jy sigma=0.414889 ! Fit after self-cal, rms=0.097065Jy sigma=0.414176 ! Inverting map ! Added new window around map position (11.2, 10). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0055334 Jy ! Component: 100 - total flux cleaned = 0.00707893 Jy ! Component: 150 - total flux cleaned = 0.00848409 Jy ! Component: 200 - total flux cleaned = 0.00951606 Jy ! Total flux subtracted in 200 components = 0.00951606 Jy ! Clean residual min=-0.008635 max=0.008662 Jy/beam ! Clean residual mean=0.000004 rms=0.002100 Jy/beam ! Combined flux in latest and established models = 0.91142 Jy ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 351 components and 0.91142 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.095025Jy sigma=0.405541 ! Fit after self-cal, rms=0.094869Jy sigma=0.405031 ! 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.094869Jy sigma=0.405031 ! Fit after self-cal, rms=0.093385Jy sigma=0.401903 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=5.002 mas, bmaj=7.156 mas, bpa=17.22 degrees ! Estimated noise=1.85595 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00396992 Jy ! Component: 100 - total flux cleaned = 0.00554694 Jy ! Component: 150 - total flux cleaned = 0.00697385 Jy ! Component: 200 - total flux cleaned = 0.00764351 Jy ! Total flux subtracted in 200 components = 0.00764351 Jy ! Clean residual min=-0.007496 max=0.007003 Jy/beam ! Clean residual mean=0.000007 rms=0.001780 Jy/beam ! Combined flux in latest and established models = 0.919064 Jy selfcal ! Performing phase self-cal ! Adding 44 model components to the UV plane model. ! The established model now contains 375 components and 0.919064 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.092084Jy sigma=0.395955 ! Fit after self-cal, rms=0.091992Jy sigma=0.395682 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.00145262 Jy ! Component: 100 - total flux cleaned = 0.00185448 Jy ! Component: 150 - total flux cleaned = 0.00259915 Jy ! Component: 200 - total flux cleaned = 0.00312223 Jy ! Total flux subtracted in 200 components = 0.00312223 Jy ! Clean residual min=-0.006881 max=0.006599 Jy/beam ! Clean residual mean=0.000009 rms=0.001633 Jy/beam ! Combined flux in latest and established models = 0.922186 Jy ! Performing phase self-cal ! Adding 56 model components to the UV plane model. ! The established model now contains 405 components and 0.922186 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.091235Jy sigma=0.392044 ! Fit after self-cal, rms=0.091150Jy sigma=0.391818 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime2 ! Performing amp+phase self-cal over 15 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.091150Jy sigma=0.391818 ! Fit after self-cal, rms=0.091152Jy sigma=0.391250 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=9.786 mas, bmaj=10.58 mas, bpa=26.59 degrees ! Estimated noise=2.09055 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00219482 Jy ! Component: 100 - total flux cleaned = 0.00377193 Jy ! Component: 150 - total flux cleaned = 0.00513839 Jy ! Component: 200 - total flux cleaned = 0.00590522 Jy ! Total flux subtracted in 200 components = 0.00590522 Jy ! Clean residual min=-0.005534 max=0.005385 Jy/beam ! Clean residual mean=0.000014 rms=0.001679 Jy/beam ! Combined flux in latest and established models = 0.928091 Jy selfcal ! Performing phase self-cal ! Adding 74 model components to the UV plane model. ! The established model now contains 468 components and 0.928091 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.090622Jy sigma=0.388527 ! Fit after self-cal, rms=0.090529Jy sigma=0.388202 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.000987179 Jy ! Component: 100 - total flux cleaned = 0.00114496 Jy ! Component: 150 - total flux cleaned = 0.0015242 Jy ! Component: 200 - total flux cleaned = 0.00175917 Jy ! Total flux subtracted in 200 components = 0.00175917 Jy ! Clean residual min=-0.004995 max=0.005075 Jy/beam ! Clean residual mean=0.000016 rms=0.001565 Jy/beam ! Combined flux in latest and established models = 0.929851 Jy ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 493 components and 0.929851 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.090236Jy sigma=0.386463 ! Fit after self-cal, rms=0.090176Jy sigma=0.386255 ! 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 3:FD,true ! 3:FD self-cal status=fixed weight=1 selfant 3:HN,true ! 3:HN self-cal status=fixed weight=1 selfant 3:LA,true ! 3:LA self-cal status=fixed weight=1 selfant 3:PT,true ! 3:PT self-cal status=fixed weight=1 selfant 3:SC,true ! 3:SC self-cal status=fixed weight=1 selfant 4:BR,true ! 4:BR self-cal status=fixed weight=1 selfant 4:FD,true ! 4:FD self-cal status=fixed weight=1 selfant 4:HN,true ! 4:HN self-cal status=fixed weight=1 selfant 4:KP,true ! 4:KP self-cal status=fixed weight=1 selfant 4:LA,true ! 4:LA self-cal status=fixed weight=1 selfant 4:OV,true ! 4:OV self-cal status=fixed weight=1 selfant 4:PT,true ! 4:PT self-cal status=fixed weight=1 selfant 4:SC,true ! 4:SC self-cal status=fixed weight=1 selfant 5:BR,true ! 5:BR self-cal status=fixed weight=1 selfant 5:FD,true ! 5:FD self-cal status=fixed weight=1 selfant 5:HN,true ! 5:HN self-cal status=fixed weight=1 selfant 5:KP,true ! 5:KP self-cal status=fixed weight=1 selfant 5:LA,true ! 5:LA self-cal status=fixed weight=1 selfant 5:OV,true ! 5:OV self-cal status=fixed weight=1 selfant 5:PT,true ! 5:PT self-cal status=fixed weight=1 selfant 5:SC,true ! 5:SC self-cal status=fixed weight=1 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: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:BR,true ! 7:BR self-cal status=fixed weight=1 selfant 7:KP,true ! 7:KP self-cal status=fixed weight=1 selfant 7:MK,true ! 7:MK self-cal status=fixed weight=1 selfant 7:OV,true ! 7:OV 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 90 telescope corrections were flagged in sub-array 1. ! A total of 18 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 8 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 1 telescope corrections were flagged in sub-array 6. ! A total of 1 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 2. ! A total of 90 telescope corrections were flagged in sub-array 1. ! A total of 18 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 8 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 9 telescope corrections were flagged in sub-array 6. ! A total of 7 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 3. ! A total of 90 telescope corrections were flagged in sub-array 1. ! A total of 18 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 8 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 1 telescope corrections were flagged in sub-array 6. ! A total of 1 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 4. ! A total of 90 telescope corrections were flagged in sub-array 1. ! A total of 18 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 8 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 1 telescope corrections were flagged in sub-array 6. ! A total of 1 telescope corrections were flagged in sub-array 7. ! ! Fit before self-cal, rms=0.133522Jy sigma=0.512523 ! Fit after self-cal, rms=0.138343Jy sigma=0.445069 selfant "",false ! 1:NY self-cal status=correctable weight=1 ! 1:WZ self-cal status=correctable weight=1 ! 1:ZC self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 2:ZC self-cal status=correctable weight=1 ! 3:FD self-cal status=correctable weight=1 ! 3:HN self-cal status=correctable weight=1 ! 3:LA self-cal status=correctable weight=1 ! 3:NY self-cal status=correctable weight=1 ! 3:PT self-cal status=correctable weight=1 ! 3:SC self-cal status=correctable weight=1 ! 3:WZ self-cal status=correctable weight=1 ! 3:ZC self-cal status=correctable weight=1 ! 4:BR self-cal status=correctable weight=1 ! 4:FD self-cal status=correctable weight=1 ! 4:HN self-cal status=correctable weight=1 ! 4:KP self-cal status=correctable weight=1 ! 4:LA self-cal status=correctable weight=1 ! 4:NY self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 4:PT self-cal status=correctable weight=1 ! 4:SC self-cal status=correctable weight=1 ! 4:WZ self-cal status=correctable weight=1 ! 4:ZC self-cal status=correctable weight=1 ! 5:BR self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:HN self-cal status=correctable weight=1 ! 5:KP self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:NY self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 5:SC self-cal status=correctable weight=1 ! 5:WZ 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:OV self-cal status=correctable weight=1 ! 6:PT self-cal status=correctable weight=1 ! 6:SC self-cal status=correctable weight=1 ! 7:BR self-cal status=correctable weight=1 ! 7:KK self-cal status=correctable weight=1 ! 7:KP self-cal status=correctable weight=1 ! 7:MK self-cal status=correctable weight=1 ! 7:OV 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 16 telescope corrections were flagged in sub-array 3. ! A total of 18 telescope corrections were flagged in sub-array 4. ! A total of 38 telescope corrections were flagged in sub-array 5. ! A total of 16 telescope corrections were flagged in sub-array 6. ! A total of 6 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 2. ! A total of 16 telescope corrections were flagged in sub-array 3. ! A total of 18 telescope corrections were flagged in sub-array 4. ! A total of 38 telescope corrections were flagged in sub-array 5. ! A total of 16 telescope corrections were flagged in sub-array 6. ! A total of 12 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 3. ! A total of 16 telescope corrections were flagged in sub-array 3. ! A total of 18 telescope corrections were flagged in sub-array 4. ! A total of 38 telescope corrections were flagged in sub-array 5. ! A total of 16 telescope corrections were flagged in sub-array 6. ! A total of 6 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 4. ! A total of 16 telescope corrections were flagged in sub-array 3. ! A total of 48 telescope corrections were flagged in sub-array 4. ! A total of 38 telescope corrections were flagged in sub-array 5. ! A total of 16 telescope corrections were flagged in sub-array 6. ! A total of 6 telescope corrections were flagged in sub-array 7. ! ! Fit before self-cal, rms=0.138343Jy sigma=0.445069 ! Fit after self-cal, rms=0.135836Jy sigma=0.437457 uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.853 mas, bmaj=5.506 mas, bpa=-10.42 degrees ! Estimated noise=2.69837 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.00823227 Jy ! Component: 100 - total flux cleaned = -0.00676313 Jy ! Component: 150 - total flux cleaned = -0.00560356 Jy ! Component: 200 - total flux cleaned = -0.00453706 Jy ! Total flux subtracted in 200 components = -0.00453706 Jy ! Clean residual min=-0.009494 max=0.010640 Jy/beam ! Clean residual mean=-0.000002 rms=0.002196 Jy/beam ! Combined flux in latest and established models = 0.925313 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 523 components and 0.925314 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.129416Jy sigma=0.417923 ! Fit after self-cal, rms=0.127561Jy sigma=0.412935 ! Inverting map ! Added new window around map position (-10, -8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00262223 Jy ! Component: 100 - total flux cleaned = 0.00352457 Jy ! Component: 150 - total flux cleaned = 0.00354907 Jy ! Component: 200 - total flux cleaned = 0.0039279 Jy ! Total flux subtracted in 200 components = 0.0039279 Jy ! Clean residual min=-0.008320 max=0.008362 Jy/beam ! Clean residual mean=-0.000002 rms=0.001911 Jy/beam ! Combined flux in latest and established models = 0.929241 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 541 components and 0.929241 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.125045Jy sigma=0.404880 ! Fit after self-cal, rms=0.124859Jy sigma=0.404620 ! Inverting map ! Added new window around map position (-3.6, -10). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000737186 Jy ! Component: 100 - total flux cleaned = 0.000395897 Jy ! Component: 150 - total flux cleaned = -0.000566217 Jy ! Component: 200 - total flux cleaned = -0.0011591 Jy ! Total flux subtracted in 200 components = -0.0011591 Jy ! Clean residual min=-0.007647 max=0.007825 Jy/beam ! Clean residual mean=-0.000002 rms=0.001778 Jy/beam ! Combined flux in latest and established models = 0.928082 Jy ! Performing phase self-cal ! Adding 51 model components to the UV plane model. ! The established model now contains 556 components and 0.928082 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.123380Jy sigma=0.400108 ! Fit after self-cal, rms=0.123269Jy sigma=0.399966 ! Inverting map ! Added new window around map position (7.6, -35.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0039132 Jy ! Component: 100 - total flux cleaned = 0.00363325 Jy ! Component: 150 - total flux cleaned = 0.00362267 Jy ! Component: 200 - total flux cleaned = 0.00286875 Jy ! Total flux subtracted in 200 components = 0.00286875 Jy ! Clean residual min=-0.007003 max=0.007068 Jy/beam ! Clean residual mean=-0.000002 rms=0.001675 Jy/beam ! Combined flux in latest and established models = 0.930951 Jy ! Performing phase self-cal ! Adding 61 model components to the UV plane model. ! The established model now contains 578 components and 0.930951 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.122385Jy sigma=0.397354 ! Fit after self-cal, rms=0.122350Jy sigma=0.397268 ! Inverting map ! Added new window around map position (22.8, 74.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00299992 Jy ! Component: 100 - total flux cleaned = 0.00251924 Jy ! Component: 150 - total flux cleaned = 0.00183448 Jy ! Component: 200 - total flux cleaned = 0.00205136 Jy ! Total flux subtracted in 200 components = 0.00205136 Jy ! Clean residual min=-0.006422 max=0.006588 Jy/beam ! Clean residual mean=-0.000001 rms=0.001588 Jy/beam ! Combined flux in latest and established models = 0.933003 Jy ! Performing phase self-cal ! Adding 61 model components to the UV plane model. ! The established model now contains 599 components and 0.933003 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.121752Jy sigma=0.395197 ! Fit after self-cal, rms=0.121704Jy sigma=0.395109 ! Inverting map uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.09 mas, bmaj=6.267 mas, bpa=-10.23 degrees ! Estimated noise=1.5917 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00387718 Jy ! Component: 100 - total flux cleaned = 0.00470958 Jy ! Component: 150 - total flux cleaned = 0.00509547 Jy ! Component: 200 - total flux cleaned = 0.00545161 Jy ! Component: 250 - total flux cleaned = 0.00545589 Jy ! Component: 300 - total flux cleaned = 0.00577797 Jy ! Component: 350 - total flux cleaned = 0.00607571 Jy ! Component: 400 - total flux cleaned = 0.00607459 Jy ! Total flux subtracted in 400 components = 0.00607459 Jy ! Clean residual min=-0.005013 max=0.005026 Jy/beam ! Clean residual mean=0.000019 rms=0.001239 Jy/beam ! Combined flux in latest and established models = 0.939077 Jy ! Performing phase self-cal ! Adding 78 model components to the UV plane model. ! The established model now contains 654 components and 0.939077 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.120446Jy sigma=0.391352 ! Fit after self-cal, rms=0.120312Jy sigma=0.391064 ! Inverting map dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000417603 Jy ! Component: 100 - total flux cleaned = -0.000118553 Jy ! Component: 150 - total flux cleaned = 1.71534e-05 Jy ! Component: 200 - total flux cleaned = 0.000147728 Jy ! Component: 250 - total flux cleaned = 0.000148698 Jy ! Component: 300 - total flux cleaned = 0.000269859 Jy ! Component: 350 - total flux cleaned = 0.000385641 Jy ! Component: 400 - total flux cleaned = 0.000383411 Jy ! Total flux subtracted in 400 components = 0.000383411 Jy ! Clean residual min=-0.004682 max=0.004703 Jy/beam ! Clean residual mean=0.000020 rms=0.001165 Jy/beam ! Combined flux in latest and established models = 0.939461 Jy ! Performing phase self-cal ! Adding 90 model components to the UV plane model. ! The established model now contains 691 components and 0.939461 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.119652Jy sigma=0.389116 ! Fit after self-cal, rms=0.119570Jy sigma=0.388941 ! Inverting map print "=== Finished natural weighting clean for al sites ===" ! === Finished natural weighting clean for al sites === unflag * ! Unflagging all baselines in the currently selected channels. ! Applying 2560 buffered edits. ! Applying 2560 buffered edits. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 90 telescope corrections were flagged in sub-array 1. ! A total of 18 telescope corrections were flagged in sub-array 2. ! A total of 27 telescope corrections were flagged in sub-array 3. ! A total of 26 telescope corrections were flagged in sub-array 4. ! A total of 42 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! A total of 7 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 2. ! A total of 90 telescope corrections were flagged in sub-array 1. ! A total of 18 telescope corrections were flagged in sub-array 2. ! A total of 27 telescope corrections were flagged in sub-array 3. ! A total of 26 telescope corrections were flagged in sub-array 4. ! A total of 42 telescope corrections were flagged in sub-array 5. ! A total of 25 telescope corrections were flagged in sub-array 6. ! A total of 19 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 3. ! A total of 90 telescope corrections were flagged in sub-array 1. ! A total of 18 telescope corrections were flagged in sub-array 2. ! A total of 27 telescope corrections were flagged in sub-array 3. ! A total of 26 telescope corrections were flagged in sub-array 4. ! A total of 42 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! A total of 7 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 4. ! A total of 90 telescope corrections were flagged in sub-array 1. ! A total of 18 telescope corrections were flagged in sub-array 2. ! A total of 27 telescope corrections were flagged in sub-array 3. ! A total of 56 telescope corrections were flagged in sub-array 4. ! A total of 42 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! A total of 7 telescope corrections were flagged in sub-array 7. ! ! Fit before self-cal, rms=0.124059Jy sigma=0.402158 ! Fit after self-cal, rms=0.114830Jy sigma=0.385895 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.067 mas, bmaj=6.298 mas, bpa=-10.21 degrees ! Estimated noise=1.58372 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000496978 Jy ! Component: 100 - total flux cleaned = -0.000491615 Jy ! Component: 150 - total flux cleaned = -0.000201231 Jy ! Component: 200 - total flux cleaned = -6.70577e-05 Jy ! Component: 250 - total flux cleaned = -0.000195614 Jy ! Component: 300 - total flux cleaned = -0.00042898 Jy ! Component: 350 - total flux cleaned = -0.000316841 Jy ! Component: 400 - total flux cleaned = -0.000321426 Jy ! Total flux subtracted in 400 components = -0.000321426 Jy ! Clean residual min=-0.003622 max=0.004227 Jy/beam ! Clean residual mean=0.000012 rms=0.000969 Jy/beam ! Combined flux in latest and established models = 0.93914 Jy selfcal ! Performing phase self-cal ! Adding 68 model components to the UV plane model. ! The established model now contains 738 components and 0.93914 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.114104Jy sigma=0.383623 ! Fit after self-cal, rms=0.113987Jy sigma=0.383402 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 (62, -29.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00262606 Jy ! Component: 100 - total flux cleaned = 0.00263058 Jy ! Component: 150 - total flux cleaned = 0.00262965 Jy ! Component: 200 - total flux cleaned = 0.00285426 Jy ! Component: 250 - total flux cleaned = 0.0030674 Jy ! Component: 300 - total flux cleaned = 0.00316679 Jy ! Component: 350 - total flux cleaned = 0.00316781 Jy ! Component: 400 - total flux cleaned = 0.00345633 Jy ! Total flux subtracted in 400 components = 0.00345633 Jy ! Clean residual min=-0.003375 max=0.003753 Jy/beam ! Clean residual mean=0.000004 rms=0.000890 Jy/beam ! Combined flux in latest and established models = 0.942596 Jy ! Performing phase self-cal ! Adding 77 model components to the UV plane model. ! The established model now contains 779 components and 0.942596 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.113436Jy sigma=0.381615 ! Fit after self-cal, rms=0.113390Jy sigma=0.381480 ! Inverting map ! Added new window around map position (-23.6, 73.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00119163 Jy ! Component: 100 - total flux cleaned = 0.00131979 Jy ! Component: 150 - total flux cleaned = 0.00132782 Jy ! Component: 200 - total flux cleaned = 0.00122941 Jy ! Component: 250 - total flux cleaned = 0.00122847 Jy ! Component: 300 - total flux cleaned = 0.0012285 Jy ! Component: 350 - total flux cleaned = 0.00131539 Jy ! Component: 400 - total flux cleaned = 0.00131556 Jy ! Total flux subtracted in 400 components = 0.00131556 Jy ! Clean residual min=-0.003267 max=0.003590 Jy/beam ! Clean residual mean=0.000003 rms=0.000831 Jy/beam ! Combined flux in latest and established models = 0.943911 Jy ! Performing phase self-cal ! Adding 83 model components to the UV plane model. ! The established model now contains 817 components and 0.943911 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.112984Jy sigma=0.380145 ! Fit after self-cal, rms=0.112960Jy sigma=0.380042 ! Inverting map ! Added new window around map position (-22.8, -32.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00180033 Jy ! Component: 100 - total flux cleaned = 0.00199561 Jy ! Component: 150 - total flux cleaned = 0.00190148 Jy ! Component: 200 - total flux cleaned = 0.00206904 Jy ! Component: 250 - total flux cleaned = 0.00206926 Jy ! Component: 300 - total flux cleaned = 0.00214915 Jy ! Component: 350 - total flux cleaned = 0.00207272 Jy ! Component: 400 - total flux cleaned = 0.00207295 Jy ! Total flux subtracted in 400 components = 0.00207295 Jy ! Clean residual min=-0.003296 max=0.003431 Jy/beam ! Clean residual mean=0.000000 rms=0.000784 Jy/beam ! Combined flux in latest and established models = 0.945984 Jy ! Performing phase self-cal ! Adding 100 model components to the UV plane model. ! The established model now contains 853 components and 0.945984 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.112663Jy sigma=0.378995 ! Fit after self-cal, rms=0.112645Jy sigma=0.378936 ! 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.112645Jy sigma=0.378936 ! Fit after self-cal, rms=0.111251Jy sigma=0.378516 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.06 mas, bmaj=6.301 mas, bpa=-10.21 degrees ! Estimated noise=1.5785 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000528118 Jy ! Component: 100 - total flux cleaned = -0.000445204 Jy ! Component: 150 - total flux cleaned = -0.000189035 Jy ! Component: 200 - total flux cleaned = -0.000271179 Jy ! Component: 250 - total flux cleaned = -0.000112285 Jy ! Component: 300 - total flux cleaned = -0.000110888 Jy ! Component: 350 - total flux cleaned = 3.73073e-05 Jy ! Component: 400 - total flux cleaned = 0.000181234 Jy ! Total flux subtracted in 400 components = 0.000181234 Jy ! Clean residual min=-0.002786 max=0.003109 Jy/beam ! Clean residual mean=0.000004 rms=0.000724 Jy/beam ! Combined flux in latest and established models = 0.946166 Jy selfcal ! Performing phase self-cal ! Adding 93 model components to the UV plane model. ! The established model now contains 884 components and 0.946166 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.110981Jy sigma=0.377640 ! Fit after self-cal, rms=0.110941Jy sigma=0.377585 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 (-63.6, -22.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00234593 Jy ! Component: 100 - total flux cleaned = 0.00330151 Jy ! Component: 150 - total flux cleaned = 0.0037017 Jy ! Component: 200 - total flux cleaned = 0.00407932 Jy ! Component: 250 - total flux cleaned = 0.00436752 Jy ! Component: 300 - total flux cleaned = 0.00464715 Jy ! Component: 350 - total flux cleaned = 0.00485215 Jy ! Component: 400 - total flux cleaned = 0.00505418 Jy ! Total flux subtracted in 400 components = 0.00505418 Jy ! Clean residual min=-0.002642 max=0.002924 Jy/beam ! Clean residual mean=0.000003 rms=0.000681 Jy/beam ! Combined flux in latest and established models = 0.95122 Jy ! Performing phase self-cal ! Adding 99 model components to the UV plane model. ! The established model now contains 925 components and 0.95122 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.110702Jy sigma=0.376767 ! Fit after self-cal, rms=0.110684Jy sigma=0.376722 ! Inverting map ! Added new window around map position (2, -26). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0011896 Jy ! Component: 100 - total flux cleaned = 0.00149024 Jy ! Component: 150 - total flux cleaned = 0.00231452 Jy ! Component: 200 - total flux cleaned = 0.00287914 Jy ! Component: 250 - total flux cleaned = 0.00328645 Jy ! Component: 300 - total flux cleaned = 0.00368122 Jy ! Component: 350 - total flux cleaned = 0.0039368 Jy ! Component: 400 - total flux cleaned = 0.00425147 Jy ! Total flux subtracted in 400 components = 0.00425147 Jy ! Clean residual min=-0.002611 max=0.002688 Jy/beam ! Clean residual mean=0.000001 rms=0.000640 Jy/beam ! Combined flux in latest and established models = 0.955471 Jy ! Performing phase self-cal ! Adding 90 model components to the UV plane model. ! The established model now contains 959 components and 0.955471 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.110470Jy sigma=0.375996 ! Fit after self-cal, rms=0.110462Jy sigma=0.375962 ! Inverting map ! Added new window around map position (75.6, -18.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00168342 Jy ! Component: 100 - total flux cleaned = 0.00230661 Jy ! Component: 150 - total flux cleaned = 0.00275708 Jy ! Component: 200 - total flux cleaned = 0.00288057 Jy ! Component: 250 - total flux cleaned = 0.00306057 Jy ! Component: 300 - total flux cleaned = 0.00311822 Jy ! Component: 350 - total flux cleaned = 0.00311798 Jy ! Component: 400 - total flux cleaned = 0.00340341 Jy ! Total flux subtracted in 400 components = 0.00340341 Jy ! Clean residual min=-0.002472 max=0.002543 Jy/beam ! Clean residual mean=0.000001 rms=0.000603 Jy/beam ! Combined flux in latest and established models = 0.958875 Jy ! Performing phase self-cal ! Adding 108 model components to the UV plane model. ! The established model now contains 1002 components and 0.958875 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.110294Jy sigma=0.375419 ! Fit after self-cal, rms=0.110283Jy sigma=0.375398 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime3 ! Performing amp+phase self-cal over 1 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.110283Jy sigma=0.375398 ! Fit after self-cal, rms=0.109033Jy sigma=0.367030 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=9.375 mas, bmaj=10.52 mas, bpa=19.57 degrees ! Estimated noise=2.03566 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00199746 Jy ! Component: 100 - total flux cleaned = 0.00326257 Jy ! Component: 150 - total flux cleaned = 0.00370511 Jy ! Component: 200 - total flux cleaned = 0.0042139 Jy ! Component: 250 - total flux cleaned = 0.0046236 Jy ! Component: 300 - total flux cleaned = 0.00495796 Jy ! Component: 350 - total flux cleaned = 0.00518573 Jy ! Component: 400 - total flux cleaned = 0.00540383 Jy ! Total flux subtracted in 400 components = 0.00540383 Jy ! Clean residual min=-0.001945 max=0.001929 Jy/beam ! Clean residual mean=0.000012 rms=0.000584 Jy/beam ! Combined flux in latest and established models = 0.964279 Jy selfcal ! Performing phase self-cal ! Adding 63 model components to the UV plane model. ! The established model now contains 1057 components and 0.964279 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.109000Jy sigma=0.366702 ! Fit after self-cal, rms=0.109003Jy sigma=0.366647 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 (-88, 2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00146738 Jy ! Component: 100 - total flux cleaned = 0.00211122 Jy ! Component: 150 - total flux cleaned = 0.00263425 Jy ! Component: 200 - total flux cleaned = 0.00329685 Jy ! Component: 250 - total flux cleaned = 0.00364423 Jy ! Component: 300 - total flux cleaned = 0.0036847 Jy ! Component: 350 - total flux cleaned = 0.00372321 Jy ! Component: 400 - total flux cleaned = 0.00368479 Jy ! Total flux subtracted in 400 components = 0.00368479 Jy ! Clean residual min=-0.001837 max=0.001741 Jy/beam ! Clean residual mean=0.000014 rms=0.000518 Jy/beam ! Combined flux in latest and established models = 0.967964 Jy ! Performing phase self-cal ! Adding 76 model components to the UV plane model. ! The established model now contains 1117 components and 0.967963 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.108962Jy sigma=0.366402 ! Fit after self-cal, rms=0.108953Jy sigma=0.366382 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "=== Finished amplitude self-cal for all stations ===" ! === Finished amplitude self-cal for all stations === restore ! Inverting map and beam ! Estimated beam: bmin=2.059 mas, bmaj=6.302 mas, bpa=-10.21 degrees ! Estimated noise=1.57996 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.059 x 6.302 at -10.21 degrees (North through East) ! Clean map min=-0.005368 max=0.23135 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 print "=== CLearing model and starting over ===" ! === CLearing model and starting over === clrmod true ! clrmod: Cleared the established and tentative models. uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.829 mas, bmaj=5.527 mas, bpa=-10.34 degrees ! Estimated noise=2.68857 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.193619 Jy ! Component: 100 - total flux cleaned = 0.253713 Jy ! Component: 150 - total flux cleaned = 0.298513 Jy ! Component: 200 - total flux cleaned = 0.336713 Jy ! Component: 250 - total flux cleaned = 0.370712 Jy ! Component: 300 - total flux cleaned = 0.401544 Jy ! Component: 350 - total flux cleaned = 0.429809 Jy ! Component: 400 - total flux cleaned = 0.455869 Jy ! Total flux subtracted in 400 components = 0.455869 Jy ! Clean residual min=-0.012886 max=0.016697 Jy/beam ! Clean residual mean=-0.000060 rms=0.003692 Jy/beam ! Combined flux in latest and established models = 0.455868 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 108 model components to the UV plane model. ! The established model now contains 108 components and 0.455868 Jy ! Inverting map and beam ! Estimated beam: bmin=2.059 mas, bmaj=6.302 mas, bpa=-10.21 degrees ! Estimated noise=1.57996 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.0763832 Jy ! Component: 100 - total flux cleaned = 0.131147 Jy ! Component: 150 - total flux cleaned = 0.175819 Jy ! Component: 200 - total flux cleaned = 0.213091 Jy ! Component: 250 - total flux cleaned = 0.244515 Jy ! Component: 300 - total flux cleaned = 0.271233 Jy ! Component: 350 - total flux cleaned = 0.294079 Jy ! Component: 400 - total flux cleaned = 0.313838 Jy ! Total flux subtracted in 400 components = 0.313838 Jy ! Clean residual min=-0.009261 max=0.012304 Jy/beam ! Clean residual mean=-0.000045 rms=0.003796 Jy/beam ! Combined flux in latest and established models = 0.769706 Jy ! Adding 273 model components to the UV plane model. ! The established model now contains 358 components and 0.769706 Jy ! Inverting map ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0172098 Jy ! Component: 100 - total flux cleaned = 0.0322907 Jy ! Component: 150 - total flux cleaned = 0.045667 Jy ! Component: 200 - total flux cleaned = 0.0576275 Jy ! Component: 250 - total flux cleaned = 0.0683487 Jy ! Component: 300 - total flux cleaned = 0.077607 Jy ! Component: 350 - total flux cleaned = 0.0859682 Jy ! Component: 400 - total flux cleaned = 0.0935481 Jy ! Total flux subtracted in 400 components = 0.0935481 Jy ! Clean residual min=-0.004992 max=0.005601 Jy/beam ! Clean residual mean=0.000000 rms=0.001754 Jy/beam ! Combined flux in latest and established models = 0.863254 Jy ! Adding 197 model components to the UV plane model. ! The established model now contains 530 components and 0.863254 Jy ! Inverting map ! Added new window around map position (60, -62.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00668966 Jy ! Component: 100 - total flux cleaned = 0.0130023 Jy ! Component: 150 - total flux cleaned = 0.0180681 Jy ! Component: 200 - total flux cleaned = 0.0231841 Jy ! Component: 250 - total flux cleaned = 0.0274847 Jy ! Component: 300 - total flux cleaned = 0.0314544 Jy ! Component: 350 - total flux cleaned = 0.0353163 Jy ! Component: 400 - total flux cleaned = 0.0385824 Jy ! Total flux subtracted in 400 components = 0.0385824 Jy ! Clean residual min=-0.003505 max=0.003597 Jy/beam ! Clean residual mean=0.000011 rms=0.001052 Jy/beam ! Combined flux in latest and established models = 0.901837 Jy ! Adding 157 model components to the UV plane model. ! The established model now contains 651 components and 0.901837 Jy ! Inverting map ! Added new window around map position (11.6, -50). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00365299 Jy ! Component: 100 - total flux cleaned = 0.00659636 Jy ! Component: 150 - total flux cleaned = 0.00936424 Jy ! Component: 200 - total flux cleaned = 0.0114712 Jy ! Component: 250 - total flux cleaned = 0.0134891 Jy ! Component: 300 - total flux cleaned = 0.0154171 Jy ! Component: 350 - total flux cleaned = 0.0172691 Jy ! Component: 400 - total flux cleaned = 0.018497 Jy ! Total flux subtracted in 400 components = 0.018497 Jy ! Clean residual min=-0.002950 max=0.002873 Jy/beam ! Clean residual mean=0.000014 rms=0.000820 Jy/beam ! Combined flux in latest and established models = 0.920334 Jy ! Adding 147 model components to the UV plane model. ! The established model now contains 737 components and 0.920334 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=9.375 mas, bmaj=10.52 mas, bpa=19.57 degrees ! Estimated noise=2.03566 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00527463 Jy ! Component: 100 - total flux cleaned = 0.00869167 Jy ! Component: 150 - total flux cleaned = 0.0116803 Jy ! Component: 200 - total flux cleaned = 0.0147143 Jy ! Component: 250 - total flux cleaned = 0.0172613 Jy ! Component: 300 - total flux cleaned = 0.0193091 Jy ! Component: 350 - total flux cleaned = 0.0206472 Jy ! Component: 400 - total flux cleaned = 0.022174 Jy ! Total flux subtracted in 400 components = 0.022174 Jy ! Clean residual min=-0.002843 max=0.002848 Jy/beam ! Clean residual mean=0.000045 rms=0.000817 Jy/beam ! Combined flux in latest and established models = 0.942508 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 (10, -83.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000620685 Jy ! Component: 100 - total flux cleaned = 0.0018951 Jy ! Component: 150 - total flux cleaned = 0.00315693 Jy ! Component: 200 - total flux cleaned = 0.00403673 Jy ! Component: 250 - total flux cleaned = 0.00502818 Jy ! Component: 300 - total flux cleaned = 0.00598014 Jy ! Component: 350 - total flux cleaned = 0.00675398 Jy ! Component: 400 - total flux cleaned = 0.00750448 Jy ! Total flux subtracted in 400 components = 0.00750448 Jy ! Clean residual min=-0.002175 max=0.002173 Jy/beam ! Clean residual mean=0.000032 rms=0.000646 Jy/beam ! Combined flux in latest and established models = 0.950013 Jy ! Adding 204 model components to the UV plane model. ! The established model now contains 911 components and 0.950013 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.059 mas, bmaj=6.302 mas, bpa=-10.21 degrees ! Estimated noise=1.57996 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000123562 Jy ! Component: 100 - total flux cleaned = 2.92335e-05 Jy ! Component: 150 - total flux cleaned = 0.000394357 Jy ! Component: 200 - total flux cleaned = 0.000569936 Jy ! Component: 250 - total flux cleaned = 0.000570592 Jy ! Component: 300 - total flux cleaned = 0.000979832 Jy ! Component: 350 - total flux cleaned = 0.00129948 Jy ! Component: 400 - total flux cleaned = 0.00145449 Jy ! Total flux subtracted in 400 components = 0.00145449 Jy ! Clean residual min=-0.002361 max=0.002302 Jy/beam ! Clean residual mean=0.000010 rms=0.000576 Jy/beam ! Combined flux in latest and established models = 0.951467 Jy selfcal ! Performing phase self-cal ! Adding 99 model components to the UV plane model. ! The established model now contains 994 components and 0.951468 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.108857Jy sigma=0.367043 ! Fit after self-cal, rms=0.108708Jy sigma=0.366473 wmodel /scr/pima/rdv73_uvs/J1415+1320_S_map.mod ! Writing 994 model components to file: /scr/pima/rdv73_uvs/J1415+1320_S_map.mod wobs /scr/pima/rdv73_uvs/J1415+1320_S_uvs.fits ! Applying 1835 buffered edits. ! Writing UV FITS file: /scr/pima/rdv73_uvs/J1415+1320_S_uvs.fits wwins /scr/pima/rdv73_uvs/J1415+1320_S_map.win ! wwins: Wrote 21 windows to /scr/pima/rdv73_uvs/J1415+1320_S_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.000298872 Jy ! Component: 100 - total flux cleaned = -0.000577639 Jy ! Component: 150 - total flux cleaned = -0.000962756 Jy ! Component: 200 - total flux cleaned = -0.00133601 Jy ! Component: 250 - total flux cleaned = -0.0015549 Jy ! Component: 300 - total flux cleaned = -0.00180351 Jy ! Component: 350 - total flux cleaned = -0.0019079 Jy ! Component: 400 - total flux cleaned = -0.002113 Jy ! Component: 450 - total flux cleaned = -0.00221394 Jy ! Component: 500 - total flux cleaned = -0.00218117 Jy ! Component: 550 - total flux cleaned = -0.00221379 Jy ! Component: 600 - total flux cleaned = -0.00205354 Jy ! Component: 650 - total flux cleaned = -0.00208522 Jy ! Component: 700 - total flux cleaned = -0.00211638 Jy ! Component: 750 - total flux cleaned = -0.00217713 Jy ! Component: 800 - total flux cleaned = -0.00211539 Jy ! Component: 850 - total flux cleaned = -0.00202404 Jy ! Component: 900 - total flux cleaned = -0.00202409 Jy ! Component: 950 - total flux cleaned = -0.00205391 Jy ! Component: 1000 - total flux cleaned = -0.00205396 Jy ! Component: 1050 - total flux cleaned = -0.00205406 Jy ! Component: 1100 - total flux cleaned = -0.00214072 Jy ! Component: 1150 - total flux cleaned = -0.002141 Jy ! Component: 1200 - total flux cleaned = -0.00222609 Jy ! Component: 1250 - total flux cleaned = -0.00219811 Jy ! Component: 1300 - total flux cleaned = -0.00228186 Jy ! Component: 1350 - total flux cleaned = -0.00219823 Jy ! Component: 1400 - total flux cleaned = -0.002143 Jy ! Component: 1450 - total flux cleaned = -0.00211562 Jy ! Component: 1500 - total flux cleaned = -0.00211582 Jy ! Component: 1550 - total flux cleaned = -0.00208889 Jy ! Component: 1600 - total flux cleaned = -0.00198209 Jy ! Component: 1650 - total flux cleaned = -0.00200863 Jy ! Component: 1700 - total flux cleaned = -0.00198261 Jy ! Component: 1750 - total flux cleaned = -0.00200889 Jy ! Component: 1800 - total flux cleaned = -0.00198271 Jy ! Component: 1850 - total flux cleaned = -0.00198255 Jy ! Component: 1900 - total flux cleaned = -0.00203383 Jy ! Component: 1950 - total flux cleaned = -0.00195694 Jy ! Component: 2000 - total flux cleaned = -0.0018555 Jy ! Component: 2050 - total flux cleaned = -0.00203233 Jy ! Component: 2100 - total flux cleaned = -0.00195716 Jy ! Component: 2150 - total flux cleaned = -0.0019321 Jy ! Component: 2200 - total flux cleaned = -0.00195703 Jy ! Component: 2250 - total flux cleaned = -0.00183328 Jy ! Component: 2300 - total flux cleaned = -0.00185805 Jy ! Component: 2350 - total flux cleaned = -0.00198018 Jy ! Component: 2400 - total flux cleaned = -0.00188286 Jy ! Component: 2450 - total flux cleaned = -0.00188264 Jy ! Component: 2500 - total flux cleaned = -0.00173804 Jy ! Component: 2550 - total flux cleaned = -0.00180993 Jy ! Component: 2600 - total flux cleaned = -0.0017623 Jy ! Component: 2650 - total flux cleaned = -0.00169097 Jy ! Component: 2700 - total flux cleaned = -0.00169083 Jy ! Component: 2750 - total flux cleaned = -0.00171434 Jy ! Component: 2800 - total flux cleaned = -0.00169105 Jy ! Component: 2850 - total flux cleaned = -0.00164439 Jy ! Component: 2900 - total flux cleaned = -0.00171383 Jy ! Component: 2950 - total flux cleaned = -0.00169062 Jy ! Component: 3000 - total flux cleaned = -0.00162179 Jy ! Component: 3050 - total flux cleaned = -0.00150768 Jy ! Component: 3100 - total flux cleaned = -0.0014621 Jy ! Component: 3150 - total flux cleaned = -0.00141682 Jy ! Component: 3200 - total flux cleaned = -0.00139426 Jy ! Component: 3250 - total flux cleaned = -0.00134933 Jy ! Component: 3300 - total flux cleaned = -0.00123797 Jy ! Component: 3350 - total flux cleaned = -0.00128242 Jy ! Component: 3400 - total flux cleaned = -0.00108334 Jy ! Component: 3450 - total flux cleaned = -0.00117141 Jy ! Component: 3500 - total flux cleaned = -0.00108371 Jy ! Component: 3550 - total flux cleaned = -0.000996277 Jy ! Component: 3600 - total flux cleaned = -0.000952847 Jy ! Component: 3650 - total flux cleaned = -0.000909452 Jy ! Component: 3700 - total flux cleaned = -0.000866274 Jy ! Component: 3750 - total flux cleaned = -0.000801558 Jy ! Component: 3800 - total flux cleaned = -0.000844314 Jy ! Component: 3850 - total flux cleaned = -0.000673995 Jy ! Component: 3900 - total flux cleaned = -0.000567775 Jy ! Component: 3950 - total flux cleaned = -0.000652327 Jy ! Component: 4000 - total flux cleaned = -0.000610283 Jy ! Component: 4050 - total flux cleaned = -0.000463388 Jy ! Total flux subtracted in 4096 components = -0.000463337 Jy ! Clean residual min=-0.001398 max=0.001482 Jy/beam ! Clean residual mean=0.000009 rms=0.000393 Jy/beam ! Combined flux in latest and established models = 0.951004 Jy keep ! Adding 316 model components to the UV plane model. ! The established model now contains 1307 components and 0.951004 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv73_uvs/J1415+1320_S\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.059 x 6.302 at -10.21 degrees (North through East) ! Clean map min=-0.0023523 max=0.24171 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv73_uvs/J1415+1320_S_map.fits float peak_flux peak_flux = peak(flux,max) invert ! Inverting map float image_peak image_peak = int(10000*peak_flux)/10000 float image_rms image_rms = int(100000*imstat(rms))/100000 print "=========================================" ! ========================================= print "The image peak= ", image_peak, " rms= ", image_rms, " Jy" ! The image peak= 0.2417 rms= 0.00039 Jy quit ! Quitting program ! Log file /scr/pima/rdv73_uvs/J1415+1320_S_uvs.log closed on Sat Jan 9 17:50:30 2016