! Started logfile: /scr/pima/rdv62_uvs/J1608+1029_S_uvs.log on Fri Nov 22 22:38:37 2013 obs /scr/pima/rdv62_uvs/J1608+1029_S_uva.fits ! Reading UV FITS file: /scr/pima/rdv62_uvs/J1608+1029_S_uva.fits ! AN table 1: 11 integrations on 6 of 6 possible baselines. ! AN table 2: 11 integrations on 15 of 15 possible baselines. ! AN table 3: 3 integrations on 10 of 10 possible baselines. ! AN table 4: 15 integrations on 36 of 36 possible baselines. ! AN table 5: 22 integrations on 45 of 45 possible baselines. ! AN table 6: 4 integrations on 55 of 55 possible baselines. ! Apparent sampling: 0.415714 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1608+1029 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 2.23299e+09 8e+06 1 8e+06 ! 02 2 2.26299e+09 8e+06 1 8e+06 ! 03 3 2.35299e+09 8e+06 1 8e+06 ! 04 4 2.37299e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 3344 visibilities. select RR ! Selecting polarization: RR, channels: 1..4 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 print "field_size =",field_size ! field_size = 1024 print "field_cell =",field_cell ! field_cell = 0.4 print "taper_size =",taper_size ! taper_size = 20 ![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv62_uvs/J1608+1029_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:HH self-cal status=correctable weight=1 ! 1:NY self-cal status=correctable weight=1 ! 1:WZ self-cal status=correctable weight=1 ! 1:ZC self-cal status=correctable weight=1 ! 2:HH self-cal status=correctable weight=1 ! 2:HN self-cal status=correctable weight=1 ! 2:NY self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 2:ZC self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:KP self-cal status=correctable weight=1 ! 3:LA self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 3:PT 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:OV self-cal status=correctable weight=1 ! 4:PT self-cal status=correctable weight=1 ! 4:SC 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:NL 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 ! 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:OV self-cal status=correctable weight=1 ! 6:PT self-cal status=correctable weight=1 ! 6:SC self-cal status=correctable weight=1 system "rm /tmp/ff.ff" mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.400x0.400 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 25 telescope corrections were flagged in sub-array 2. ! A total of 5 telescope corrections were flagged in sub-array 3. ! A total of 29 telescope corrections were flagged in sub-array 4. ! A total of 169 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 25 telescope corrections were flagged in sub-array 2. ! A total of 5 telescope corrections were flagged in sub-array 3. ! A total of 29 telescope corrections were flagged in sub-array 4. ! A total of 169 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 25 telescope corrections were flagged in sub-array 2. ! A total of 5 telescope corrections were flagged in sub-array 3. ! A total of 29 telescope corrections were flagged in sub-array 4. ! A total of 169 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 25 telescope corrections were flagged in sub-array 2. ! A total of 5 telescope corrections were flagged in sub-array 3. ! A total of 29 telescope corrections were flagged in sub-array 4. ! A total of 169 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=1.812027Jy sigma=26.197317 ! Fit after self-cal, rms=0.695367Jy sigma=9.788939 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. print "=== Flag out IVS stations" ! === Flag out IVS stations integer lun string cfg_file lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_flag_process.csh /tmp/ff.ff" ![@/tmp/ff.ff \] flag 1:HH ! Flagging baselines of antenna 1:HH in the currently selected channels. flag 1:NY ! Flagging baselines of antenna 1:NY in the currently selected channels. flag 1:WZ ! Flagging baselines of antenna 1:WZ in the currently selected channels. flag 1:ZC ! Flagging baselines of antenna 1:ZC in the currently selected channels. flag 2:HH ! Flagging baselines of antenna 2:HH in the currently selected channels. flag 2:NY ! Flagging baselines of antenna 2:NY in the currently selected channels. flag 2:WZ ! Flagging baselines of antenna 2:WZ in the currently selected channels. flag 2:ZC ! Flagging baselines of antenna 2:ZC in the currently selected channels. flag 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 6:KK ! Flagging baselines of antenna 6:KK in the currently selected channels. ![Exited script file: /tmp/ff.ff] uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.619 mas, bmaj=6.39 mas, bpa=-0.8771 degrees ! Estimated noise=1.66154 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 = 1.17964 Jy ! Component: 100 - total flux cleaned = 1.49721 Jy ! Total flux subtracted in 100 components = 1.49721 Jy ! Clean residual min=-0.042190 max=0.127612 Jy/beam ! Clean residual mean=-0.000178 rms=0.013759 Jy/beam ! Combined flux in latest and established models = 1.49721 Jy ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 8 components and 1.49721 Jy ! ! Correcting IF 1. ! A total of 32 telescope corrections were flagged in sub-array 1. ! A total of 41 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 4. ! A total of 3 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 32 telescope corrections were flagged in sub-array 1. ! A total of 41 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 4. ! A total of 3 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 32 telescope corrections were flagged in sub-array 1. ! A total of 41 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 4. ! A total of 3 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 32 telescope corrections were flagged in sub-array 1. ! A total of 41 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 4. ! A total of 3 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.366595Jy sigma=5.557552 ! Fit after self-cal, rms=0.368670Jy sigma=5.594075 ! Inverting map and beam ! Estimated beam: bmin=3.61 mas, bmaj=6.43 mas, bpa=-0.6372 degrees ! Estimated noise=1.67004 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.131981 Jy ! Component: 100 - total flux cleaned = 0.197881 Jy ! Total flux subtracted in 100 components = 0.197881 Jy ! Clean residual min=-0.032807 max=0.034333 Jy/beam ! Clean residual mean=-0.000138 rms=0.008550 Jy/beam ! Combined flux in latest and established models = 1.69509 Jy ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 34 components and 1.69509 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.279435Jy sigma=4.221469 ! Fit after self-cal, rms=0.279412Jy sigma=4.220941 ! 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=5.552 mas, bmaj=9.238 mas, bpa=-11.65 degrees ! Estimated noise=0.995375 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.0771788 Jy ! Component: 100 - total flux cleaned = 0.114154 Jy ! Component: 150 - total flux cleaned = 0.134312 Jy ! Component: 200 - total flux cleaned = 0.145681 Jy ! Total flux subtracted in 200 components = 0.145681 Jy ! Clean residual min=-0.038077 max=0.028394 Jy/beam ! Clean residual mean=-0.000090 rms=0.009978 Jy/beam ! Combined flux in latest and established models = 1.84077 Jy ! Performing phase self-cal ! Adding 57 model components to the UV plane model. ! The established model now contains 83 components and 1.84077 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.258049Jy sigma=3.924596 ! Fit after self-cal, rms=0.257987Jy sigma=3.923915 ! Inverting map print "=== Finished natural weighting clean for VLBI sites ===" ! === Finished natural weighting clean for VLBI sites === gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* NY 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! HH 0.00* HN 0.00* NY 0.00* SC 0.00* ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.04 KP 1.03 LA 1.15 NL 1.15 ! PT 1.11 ! ! A total of 3 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.99 FD 0.97 HN 0.95 KP 0.99 ! LA 1.14 OV 0.98 PT 1.04 SC 0.91 ! ZC 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.00 FD 1.01 HN 0.96 KP 1.00 ! LA 1.19 NL 1.04 NY 1.00* OV 0.99 ! PT 1.06 SC 0.90 ! ! Telescope amplitude corrections in sub-array 6: ! BR 0.99 FD 0.97 HN 0.93 KK 1.00* ! KP 1.08 LA 1.14 MK 1.01 NL 1.01 ! OV 1.00 PT 1.05 SC 0.93 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* NY 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! HH 0.00* HN 0.00* NY 0.00* SC 0.00* ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.08 KP 1.04 LA 1.06 NL 1.15 ! PT 1.09 ! ! A total of 3 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.01 FD 1.03 HN 1.08 KP 1.02 ! LA 1.04 OV 1.04 PT 1.04 SC 1.03 ! ZC 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.03 FD 1.09 HN 1.05 KP 1.01 ! LA 1.10 NL 1.07 NY 1.00* OV 1.04 ! PT 1.05 SC 0.99 ! ! Telescope amplitude corrections in sub-array 6: ! BR 1.02 FD 1.05 HN 1.08 KK 1.00* ! KP 1.04 LA 1.05 MK 1.02 NL 1.03 ! OV 1.05 PT 1.04 SC 1.01 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* NY 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! HH 0.00* HN 0.00* NY 0.00* SC 0.00* ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.12 KP 1.04 LA 0.87 NL 1.21 ! PT 1.10 ! ! A total of 3 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.09 FD 1.05 HN 1.02 KP 0.99 ! LA 0.83 OV 1.07 PT 1.05 SC 1.02 ! ZC 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.02 FD 1.09 HN 1.02 KP 0.99 ! LA 0.86 NL 1.14 NY 1.00* OV 1.08 ! PT 1.07 SC 0.97 ! ! Telescope amplitude corrections in sub-array 6: ! BR 1.09 FD 1.06 HN 1.00 KK 1.00* ! KP 1.02 LA 0.84 MK 1.05 NL 1.09 ! OV 1.09 PT 1.06 SC 1.03 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* NY 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! HH 0.00* HN 0.00* NY 0.00* SC 0.00* ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.82 KP 0.98 LA 0.75 NL 1.04 ! PT 1.03 ! ! A total of 3 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.98 FD 1.02 HN 0.98 KP 0.94 ! LA 0.72 OV 0.90 PT 1.00 SC 1.10 ! ZC 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 0.98 FD 1.07 HN 0.97 KP 0.94 ! LA 0.75 NL 1.04 NY 1.00* OV 0.90 ! PT 1.01 SC 1.08 ! ! Telescope amplitude corrections in sub-array 6: ! BR 0.99 FD 1.04 HN 0.96 KK 1.00* ! KP 0.96 LA 0.74 MK 1.04 NL 0.99 ! OV 0.91 PT 1.01 SC 1.13 ! ! ! Fit before self-cal, rms=0.257987Jy sigma=3.923915 ! Fit after self-cal, rms=0.089151Jy sigma=1.298884 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=5.577 mas, bmaj=9.274 mas, bpa=-11.4 degrees ! Estimated noise=1.00335 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.00561316 Jy ! Component: 100 - total flux cleaned = -0.00632669 Jy ! Component: 150 - total flux cleaned = -0.00668667 Jy ! Component: 200 - total flux cleaned = -0.00699658 Jy ! Total flux subtracted in 200 components = -0.00699658 Jy ! Clean residual min=-0.007920 max=0.009183 Jy/beam ! Clean residual mean=-0.000006 rms=0.002117 Jy/beam ! Combined flux in latest and established models = 1.83377 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 91 components and 1.83377 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.088655Jy sigma=1.291539 ! Fit after self-cal, rms=0.088722Jy sigma=1.290617 ! 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.088722Jy sigma=1.290617 ! Fit after self-cal, rms=0.088607Jy sigma=1.288777 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=5.582 mas, bmaj=9.286 mas, bpa=-11.38 degrees ! Estimated noise=1.00313 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00424859 Jy ! Component: 100 - total flux cleaned = -0.00485225 Jy ! Component: 150 - total flux cleaned = -0.00525976 Jy ! Component: 200 - total flux cleaned = -0.00552939 Jy ! Total flux subtracted in 200 components = -0.00552939 Jy ! Clean residual min=-0.007723 max=0.009146 Jy/beam ! Clean residual mean=-0.000005 rms=0.002029 Jy/beam ! Combined flux in latest and established models = 1.82824 Jy selfcal ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 94 components and 1.82824 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.088249Jy sigma=1.283755 ! Fit after self-cal, rms=0.088220Jy sigma=1.283457 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 (10.8, -38.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00807903 Jy ! Component: 100 - total flux cleaned = 0.00972069 Jy ! Component: 150 - total flux cleaned = 0.0105865 Jy ! Component: 200 - total flux cleaned = 0.0112159 Jy ! Total flux subtracted in 200 components = 0.0112159 Jy ! Clean residual min=-0.006810 max=0.006818 Jy/beam ! Clean residual mean=-0.000013 rms=0.001712 Jy/beam ! Combined flux in latest and established models = 1.83946 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 119 components and 1.83946 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.087277Jy sigma=1.270009 ! Fit after self-cal, rms=0.087240Jy sigma=1.269484 ! Inverting map ! Added new window around map position (-5.2, 9.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00395409 Jy ! Component: 100 - total flux cleaned = 0.00361366 Jy ! Component: 150 - total flux cleaned = 0.00360297 Jy ! Component: 200 - total flux cleaned = 0.0035971 Jy ! Total flux subtracted in 200 components = 0.0035971 Jy ! Clean residual min=-0.006819 max=0.005940 Jy/beam ! Clean residual mean=-0.000000 rms=0.001535 Jy/beam ! Combined flux in latest and established models = 1.84306 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 154 components and 1.84306 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.086505Jy sigma=1.258366 ! Fit after self-cal, rms=0.086447Jy sigma=1.257579 ! 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.086447Jy sigma=1.257579 ! Fit after self-cal, rms=0.086378Jy sigma=1.256137 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=5.585 mas, bmaj=9.296 mas, bpa=-11.38 degrees ! Estimated noise=1.00308 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000513716 Jy ! Component: 100 - total flux cleaned = -8.79486e-05 Jy ! Component: 150 - total flux cleaned = -8.35952e-05 Jy ! Component: 200 - total flux cleaned = -0.000424255 Jy ! Total flux subtracted in 200 components = -0.000424255 Jy ! Clean residual min=-0.006449 max=0.005556 Jy/beam ! Clean residual mean=0.000001 rms=0.001459 Jy/beam ! Combined flux in latest and established models = 1.84263 Jy selfcal ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 170 components and 1.84263 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.085946Jy sigma=1.250246 ! Fit after self-cal, rms=0.085918Jy sigma=1.249924 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 = 3.22194e-05 Jy ! Component: 100 - total flux cleaned = -0.000280019 Jy ! Component: 150 - total flux cleaned = -0.00067278 Jy ! Component: 200 - total flux cleaned = -0.000949843 Jy ! Total flux subtracted in 200 components = -0.000949843 Jy ! Clean residual min=-0.006380 max=0.005284 Jy/beam ! Clean residual mean=0.000002 rms=0.001416 Jy/beam ! Combined flux in latest and established models = 1.84168 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 179 components and 1.84168 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.085648Jy sigma=1.246685 ! Fit after self-cal, rms=0.085627Jy sigma=1.246491 ! 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.085627Jy sigma=1.246491 ! Fit after self-cal, rms=0.085570Jy sigma=1.245593 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=9.881 mas, bmaj=13.09 mas, bpa=-33.4 degrees ! Estimated noise=1.08005 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000605654 Jy ! Component: 100 - total flux cleaned = -0.000410914 Jy ! Component: 150 - total flux cleaned = -0.000414869 Jy ! Component: 200 - total flux cleaned = -0.000619345 Jy ! Total flux subtracted in 200 components = -0.000619345 Jy ! Clean residual min=-0.005401 max=0.005304 Jy/beam ! Clean residual mean=0.000002 rms=0.001466 Jy/beam ! Combined flux in latest and established models = 1.84106 Jy selfcal ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 207 components and 1.84106 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.085426Jy sigma=1.243085 ! Fit after self-cal, rms=0.085400Jy sigma=1.242774 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.000149511 Jy ! Component: 100 - total flux cleaned = 2.55102e-05 Jy ! Component: 150 - total flux cleaned = -0.000209901 Jy ! Component: 200 - total flux cleaned = -0.000426749 Jy ! Total flux subtracted in 200 components = -0.000426749 Jy ! Clean residual min=-0.005305 max=0.005129 Jy/beam ! Clean residual mean=0.000003 rms=0.001439 Jy/beam ! Combined flux in latest and established models = 1.84064 Jy ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 231 components and 1.84064 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.085334Jy sigma=1.241548 ! Fit after self-cal, rms=0.085318Jy sigma=1.241377 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "=== Finished amplitude self-cal for VLBA stations ===" ! === Finished amplitude self-cal for VLBA stations === unflag * ! Unflagging all baselines in the currently selected channels. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. integer lun lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_selfant_process.csh /tmp/ff.ff" ![@/tmp/ff.ff] selfant 2:HN,true ! 2:HN self-cal status=fixed weight=1 selfant 2:SC,true ! 2:SC self-cal status=fixed weight=1 selfant 3:BR,true ! 3:BR self-cal status=fixed weight=1 selfant 3:KP,true ! 3:KP self-cal status=fixed weight=1 selfant 3:LA,true ! 3:LA self-cal status=fixed weight=1 selfant 3:NL,true ! 3:NL self-cal status=fixed weight=1 selfant 3:PT,true ! 3:PT 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:NL,true ! 5:NL 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: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 ![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 12 telescope corrections were flagged in sub-array 1. ! A total of 24 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 19 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 24 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 19 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 24 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 19 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 24 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 19 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.511882Jy sigma=5.650025 ! Fit after self-cal, rms=0.484790Jy sigma=5.501381 selfant "",false ! 1:HH self-cal status=correctable weight=1 ! 1:NY self-cal status=correctable weight=1 ! 1:WZ self-cal status=correctable weight=1 ! 1:ZC self-cal status=correctable weight=1 ! 2:HH self-cal status=correctable weight=1 ! 2:HN self-cal status=correctable weight=1 ! 2:NY self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 2:ZC self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:KP self-cal status=correctable weight=1 ! 3:LA self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 3:PT 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:OV self-cal status=correctable weight=1 ! 4:PT self-cal status=correctable weight=1 ! 4:SC 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:NL 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 ! 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:OV self-cal status=correctable weight=1 ! 6:PT self-cal status=correctable weight=1 ! 6:SC self-cal status=correctable weight=1 selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 13 telescope corrections were flagged in sub-array 2. ! A total of 5 telescope corrections were flagged in sub-array 3. ! A total of 29 telescope corrections were flagged in sub-array 4. ! A total of 153 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 13 telescope corrections were flagged in sub-array 2. ! A total of 5 telescope corrections were flagged in sub-array 3. ! A total of 29 telescope corrections were flagged in sub-array 4. ! A total of 153 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 13 telescope corrections were flagged in sub-array 2. ! A total of 5 telescope corrections were flagged in sub-array 3. ! A total of 29 telescope corrections were flagged in sub-array 4. ! A total of 153 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 13 telescope corrections were flagged in sub-array 2. ! A total of 5 telescope corrections were flagged in sub-array 3. ! A total of 29 telescope corrections were flagged in sub-array 4. ! A total of 153 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.484790Jy sigma=5.501381 ! Fit after self-cal, rms=0.292162Jy sigma=1.420043 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=2.737 mas, bmaj=5.059 mas, bpa=-9.933 degrees ! Estimated noise=1.88435 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.00182091 Jy ! Component: 100 - total flux cleaned = -0.00345638 Jy ! Component: 150 - total flux cleaned = -0.00503287 Jy ! Component: 200 - total flux cleaned = -0.00652459 Jy ! Total flux subtracted in 200 components = -0.00652459 Jy ! Clean residual min=-0.012216 max=0.015362 Jy/beam ! Clean residual mean=0.000004 rms=0.002660 Jy/beam ! Combined flux in latest and established models = 1.83411 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 258 components and 1.83411 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.287620Jy sigma=1.411833 ! Fit after self-cal, rms=0.284548Jy sigma=1.403073 ! Inverting map ! Added new window around map position (4.8, -4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00016666 Jy ! Component: 100 - total flux cleaned = -0.000730934 Jy ! Component: 150 - total flux cleaned = -0.000236034 Jy ! Component: 200 - total flux cleaned = 0.000708667 Jy ! Total flux subtracted in 200 components = 0.000708667 Jy ! Clean residual min=-0.010171 max=0.010409 Jy/beam ! Clean residual mean=0.000004 rms=0.002428 Jy/beam ! Combined flux in latest and established models = 1.83482 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 273 components and 1.83482 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.280904Jy sigma=1.396663 ! Fit after self-cal, rms=0.280868Jy sigma=1.396332 ! Inverting map ! Added new window around map position (-2, 100.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00222227 Jy ! Component: 100 - total flux cleaned = 0.0026449 Jy ! Component: 150 - total flux cleaned = 0.00306626 Jy ! Component: 200 - total flux cleaned = 0.00426142 Jy ! Total flux subtracted in 200 components = 0.00426142 Jy ! Clean residual min=-0.009581 max=0.009784 Jy/beam ! Clean residual mean=0.000002 rms=0.002317 Jy/beam ! Combined flux in latest and established models = 1.83908 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 284 components and 1.83908 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.278396Jy sigma=1.392278 ! Fit after self-cal, rms=0.278399Jy sigma=1.392024 ! 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=3.762 mas, bmaj=7.12 mas, bpa=-12.89 degrees ! Estimated noise=1.06548 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 (16.8, -64.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00404148 Jy ! Component: 100 - total flux cleaned = 0.0047524 Jy ! Component: 150 - total flux cleaned = 0.00519812 Jy ! Component: 200 - total flux cleaned = 0.00584717 Jy ! Component: 250 - total flux cleaned = 0.00604743 Jy ! Component: 300 - total flux cleaned = 0.00624341 Jy ! Component: 350 - total flux cleaned = 0.00644025 Jy ! Component: 400 - total flux cleaned = 0.00662671 Jy ! Total flux subtracted in 400 components = 0.00662671 Jy ! Clean residual min=-0.005472 max=0.005410 Jy/beam ! Clean residual mean=0.000000 rms=0.001359 Jy/beam ! Combined flux in latest and established models = 1.84571 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 308 components and 1.84571 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.275670Jy sigma=1.385922 ! Fit after self-cal, rms=0.275565Jy sigma=1.385572 ! 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.00101535 Jy ! Component: 100 - total flux cleaned = 0.00156898 Jy ! Component: 150 - total flux cleaned = 0.00158286 Jy ! Component: 200 - total flux cleaned = 0.00193617 Jy ! Component: 250 - total flux cleaned = 0.00210341 Jy ! Component: 300 - total flux cleaned = 0.00210898 Jy ! Component: 350 - total flux cleaned = 0.00243715 Jy ! Component: 400 - total flux cleaned = 0.00259872 Jy ! Total flux subtracted in 400 components = 0.00259872 Jy ! Clean residual min=-0.005156 max=0.005170 Jy/beam ! Clean residual mean=0.000005 rms=0.001292 Jy/beam ! Combined flux in latest and established models = 1.84831 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 322 components and 1.84831 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.273635Jy sigma=1.382450 ! Fit after self-cal, rms=0.273620Jy sigma=1.382253 ! Inverting map ! Added new window around map position (49.2, 26.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00214526 Jy ! Component: 100 - total flux cleaned = 0.0036868 Jy ! Component: 150 - total flux cleaned = 0.00433734 Jy ! Component: 200 - total flux cleaned = 0.00465246 Jy ! Component: 250 - total flux cleaned = 0.00465144 Jy ! Component: 300 - total flux cleaned = 0.00480286 Jy ! Component: 350 - total flux cleaned = 0.00510612 Jy ! Component: 400 - total flux cleaned = 0.00525103 Jy ! Total flux subtracted in 400 components = 0.00525103 Jy ! Clean residual min=-0.004819 max=0.004582 Jy/beam ! Clean residual mean=0.000004 rms=0.001192 Jy/beam ! Combined flux in latest and established models = 1.85356 Jy ! Performing phase self-cal ! Adding 63 model components to the UV plane model. ! The established model now contains 347 components and 1.85356 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.272081Jy sigma=1.378568 ! Fit after self-cal, rms=0.272038Jy sigma=1.378367 ! 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. 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 12 telescope corrections were flagged in sub-array 1. ! A total of 37 telescope corrections were flagged in sub-array 2. ! A total of 5 telescope corrections were flagged in sub-array 3. ! A total of 32 telescope corrections were flagged in sub-array 4. ! A total of 172 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 37 telescope corrections were flagged in sub-array 2. ! A total of 5 telescope corrections were flagged in sub-array 3. ! A total of 32 telescope corrections were flagged in sub-array 4. ! A total of 172 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 37 telescope corrections were flagged in sub-array 2. ! A total of 5 telescope corrections were flagged in sub-array 3. ! A total of 32 telescope corrections were flagged in sub-array 4. ! A total of 172 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 37 telescope corrections were flagged in sub-array 2. ! A total of 5 telescope corrections were flagged in sub-array 3. ! A total of 32 telescope corrections were flagged in sub-array 4. ! A total of 172 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.482226Jy sigma=5.446634 ! Fit after self-cal, rms=0.265039Jy sigma=1.375318 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.721 mas, bmaj=7.105 mas, bpa=-12.68 degrees ! Estimated noise=1.06312 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000487568 Jy ! Component: 100 - total flux cleaned = 0.000494952 Jy ! Component: 150 - total flux cleaned = 0.000500495 Jy ! Component: 200 - total flux cleaned = 0.000504072 Jy ! Component: 250 - total flux cleaned = 0.000669054 Jy ! Component: 300 - total flux cleaned = 0.000661512 Jy ! Component: 350 - total flux cleaned = 0.000832322 Jy ! Component: 400 - total flux cleaned = 0.00115403 Jy ! Total flux subtracted in 400 components = 0.00115403 Jy ! Clean residual min=-0.004624 max=0.004429 Jy/beam ! Clean residual mean=0.000006 rms=0.001144 Jy/beam ! Combined flux in latest and established models = 1.85471 Jy selfcal ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 361 components and 1.85471 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.262858Jy sigma=1.372713 ! Fit after self-cal, rms=0.262864Jy sigma=1.372634 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000308554 Jy ! Component: 100 - total flux cleaned = 0.000463878 Jy ! Component: 150 - total flux cleaned = 0.000760646 Jy ! Component: 200 - total flux cleaned = 0.00104204 Jy ! Component: 250 - total flux cleaned = 0.00118526 Jy ! Component: 300 - total flux cleaned = 0.00146092 Jy ! Component: 350 - total flux cleaned = 0.00132077 Jy ! Component: 400 - total flux cleaned = 0.0014574 Jy ! Total flux subtracted in 400 components = 0.0014574 Jy ! Clean residual min=-0.004465 max=0.004259 Jy/beam ! Clean residual mean=0.000004 rms=0.001102 Jy/beam ! Combined flux in latest and established models = 1.85617 Jy ! Performing phase self-cal ! Adding 49 model components to the UV plane model. ! The established model now contains 371 components and 1.85617 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.261373Jy sigma=1.370842 ! Fit after self-cal, rms=0.261360Jy sigma=1.370766 ! 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.261360Jy sigma=1.370766 ! Fit after self-cal, rms=0.260204Jy sigma=1.368555 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.706 mas, bmaj=7.108 mas, bpa=-12.6 degrees ! Estimated noise=1.06214 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -3.76717e-06 Jy ! Component: 100 - total flux cleaned = 4.7386e-06 Jy ! Component: 150 - total flux cleaned = 2.7603e-06 Jy ! Component: 200 - total flux cleaned = 0.000151624 Jy ! Component: 250 - total flux cleaned = 0.000425595 Jy ! Component: 300 - total flux cleaned = 0.000692994 Jy ! Component: 350 - total flux cleaned = 0.000823532 Jy ! Component: 400 - total flux cleaned = 0.000817029 Jy ! Total flux subtracted in 400 components = 0.000817029 Jy ! Clean residual min=-0.004325 max=0.003966 Jy/beam ! Clean residual mean=0.000003 rms=0.001055 Jy/beam ! Combined flux in latest and established models = 1.85699 Jy selfcal ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 383 components and 1.85699 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.258751Jy sigma=1.366859 ! Fit after self-cal, rms=0.258811Jy sigma=1.366808 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000113568 Jy ! Component: 100 - total flux cleaned = 0.000115191 Jy ! Component: 150 - total flux cleaned = 0.000112356 Jy ! Component: 200 - total flux cleaned = 0.000232516 Jy ! Component: 250 - total flux cleaned = 0.000231015 Jy ! Component: 300 - total flux cleaned = 0.000349209 Jy ! Component: 350 - total flux cleaned = 0.000575408 Jy ! Component: 400 - total flux cleaned = 0.000685815 Jy ! Total flux subtracted in 400 components = 0.000685815 Jy ! Clean residual min=-0.004153 max=0.003839 Jy/beam ! Clean residual mean=0.000000 rms=0.001022 Jy/beam ! Combined flux in latest and established models = 1.85767 Jy ! Performing phase self-cal ! Adding 62 model components to the UV plane model. ! The established model now contains 394 components and 1.85767 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.257864Jy sigma=1.365253 ! Fit after self-cal, rms=0.257859Jy sigma=1.365195 ! 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.257859Jy sigma=1.365195 ! Fit after self-cal, rms=0.380145Jy sigma=1.272911 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=9.67 mas, bmaj=12.8 mas, bpa=-31.45 degrees ! Estimated noise=1.06886 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000554096 Jy ! Component: 100 - total flux cleaned = 0.00102105 Jy ! Component: 150 - total flux cleaned = 0.00144463 Jy ! Component: 200 - total flux cleaned = 0.00177059 Jy ! Component: 250 - total flux cleaned = 0.00189539 Jy ! Component: 300 - total flux cleaned = 0.00207482 Jy ! Component: 350 - total flux cleaned = 0.00224825 Jy ! Component: 400 - total flux cleaned = 0.00258338 Jy ! Total flux subtracted in 400 components = 0.00258338 Jy ! Clean residual min=-0.003474 max=0.003258 Jy/beam ! Clean residual mean=-0.000004 rms=0.000905 Jy/beam ! Combined flux in latest and established models = 1.86026 Jy selfcal ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 436 components and 1.86026 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.380106Jy sigma=1.271098 ! Fit after self-cal, rms=0.383248Jy sigma=1.270329 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000114686 Jy ! Component: 100 - total flux cleaned = 0.000278298 Jy ! Component: 150 - total flux cleaned = 0.000328612 Jy ! Component: 200 - total flux cleaned = 0.000329321 Jy ! Component: 250 - total flux cleaned = 0.000379551 Jy ! Component: 300 - total flux cleaned = 0.000475783 Jy ! Component: 350 - total flux cleaned = 0.000568353 Jy ! Component: 400 - total flux cleaned = 0.000612757 Jy ! Total flux subtracted in 400 components = 0.000612757 Jy ! Clean residual min=-0.003447 max=0.003304 Jy/beam ! Clean residual mean=-0.000005 rms=0.000885 Jy/beam ! Combined flux in latest and established models = 1.86087 Jy ! Performing phase self-cal ! Adding 48 model components to the UV plane model. ! The established model now contains 461 components and 1.86087 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.383223Jy sigma=1.269426 ! Fit after self-cal, rms=0.383235Jy sigma=1.269395 ! 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=3.827 mas, bmaj=7.164 mas, bpa=-13.34 degrees ! Estimated noise=1.06938 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.827 x 7.164 at -13.34 degrees (North through East) ! Clean map min=-0.0085114 max=1.5635 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.383235Jy sigma=1.269395 ! Fit after self-cal, rms=1.137445Jy sigma=0.993262 ! Inverting map and beam ! Estimated beam: bmin=3.788 mas, bmaj=7.151 mas, bpa=-13.1 degrees ! Estimated noise=1.06639 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00150984 Jy ! Component: 100 - total flux cleaned = -0.00182347 Jy ! Component: 150 - total flux cleaned = -0.00210842 Jy ! Component: 200 - total flux cleaned = -0.00224771 Jy ! Component: 250 - total flux cleaned = -0.00237581 Jy ! Component: 300 - total flux cleaned = -0.00249712 Jy ! Component: 350 - total flux cleaned = -0.00273758 Jy ! Component: 400 - total flux cleaned = -0.00285564 Jy ! Total flux subtracted in 400 components = -0.00285564 Jy ! Clean residual min=-0.004180 max=0.003989 Jy/beam ! Clean residual mean=-0.000009 rms=0.000971 Jy/beam ! Combined flux in latest and established models = 1.85801 Jy ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 492 components and 1.85801 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.137147Jy sigma=0.989325 ! Fit after self-cal, rms=1.137091Jy sigma=0.989249 ! Inverting map ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000134487 Jy ! Component: 100 - total flux cleaned = -0.000134869 Jy ! Component: 150 - total flux cleaned = -0.000131661 Jy ! Component: 200 - total flux cleaned = -0.000130243 Jy ! Component: 250 - total flux cleaned = -0.000235677 Jy ! Component: 300 - total flux cleaned = -0.000441978 Jy ! Component: 350 - total flux cleaned = -0.000640023 Jy ! Component: 400 - total flux cleaned = -0.000738031 Jy ! Total flux subtracted in 400 components = -0.000738031 Jy ! Clean residual min=-0.004042 max=0.003987 Jy/beam ! Clean residual mean=-0.000011 rms=0.000949 Jy/beam ! Combined flux in latest and established models = 1.85728 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 509 components and 1.85728 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.136874Jy sigma=0.987313 ! Fit after self-cal, rms=1.136840Jy sigma=0.987275 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.136840Jy sigma=0.987275 ! Fit after self-cal, rms=1.155115Jy sigma=0.986501 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=2.748 mas, bmaj=5.093 mas, bpa=-10.13 degrees ! Estimated noise=1.89111 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 1.15609 Jy ! Component: 100 - total flux cleaned = 1.47994 Jy ! Component: 150 - total flux cleaned = 1.61864 Jy ! Component: 200 - total flux cleaned = 1.68694 Jy ! Component: 250 - total flux cleaned = 1.72904 Jy ! Component: 300 - total flux cleaned = 1.75701 Jy ! Component: 350 - total flux cleaned = 1.77385 Jy ! Component: 400 - total flux cleaned = 1.78306 Jy ! Total flux subtracted in 400 components = 1.78306 Jy ! Clean residual min=-0.009542 max=0.010099 Jy/beam ! Clean residual mean=-0.000009 rms=0.002310 Jy/beam ! Combined flux in latest and established models = 1.78306 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 54 model components to the UV plane model. ! The established model now contains 54 components and 1.78306 Jy ! Inverting map and beam ! Estimated beam: bmin=3.797 mas, bmaj=7.154 mas, bpa=-13.14 degrees ! Estimated noise=1.06672 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.0185367 Jy ! Component: 100 - total flux cleaned = 0.0294451 Jy ! Component: 150 - total flux cleaned = 0.037755 Jy ! Component: 200 - total flux cleaned = 0.0444492 Jy ! Component: 250 - total flux cleaned = 0.0498195 Jy ! Component: 300 - total flux cleaned = 0.0537748 Jy ! Component: 350 - total flux cleaned = 0.0563963 Jy ! Component: 400 - total flux cleaned = 0.0580947 Jy ! Total flux subtracted in 400 components = 0.0580947 Jy ! Clean residual min=-0.005359 max=0.005313 Jy/beam ! Clean residual mean=-0.000012 rms=0.001417 Jy/beam ! Combined flux in latest and established models = 1.84115 Jy ! Adding 63 model components to the UV plane model. ! The established model now contains 113 components and 1.84115 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.674 mas, bmaj=12.79 mas, bpa=-31.32 degrees ! Estimated noise=1.06716 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00259944 Jy ! Component: 100 - total flux cleaned = 0.00487881 Jy ! Component: 150 - total flux cleaned = 0.00675605 Jy ! Component: 200 - total flux cleaned = 0.00842468 Jy ! Component: 250 - total flux cleaned = 0.0097319 Jy ! Component: 300 - total flux cleaned = 0.0104732 Jy ! Component: 350 - total flux cleaned = 0.0109951 Jy ! Component: 400 - total flux cleaned = 0.0114866 Jy ! Total flux subtracted in 400 components = 0.0114866 Jy ! Clean residual min=-0.004079 max=0.003803 Jy/beam ! Clean residual mean=-0.000022 rms=0.001165 Jy/beam ! Combined flux in latest and established models = 1.85264 Jy flux_cutoff = imstat(rms) * dynam repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam keep new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000388541 Jy ! Component: 100 - total flux cleaned = 0.000610029 Jy ! Component: 150 - total flux cleaned = 0.000893356 Jy ! Component: 200 - total flux cleaned = 0.00109548 Jy ! Component: 250 - total flux cleaned = 0.0012225 Jy ! Component: 300 - total flux cleaned = 0.00134579 Jy ! Component: 350 - total flux cleaned = 0.00140569 Jy ! Component: 400 - total flux cleaned = 0.00157601 Jy ! Total flux subtracted in 400 components = 0.00157601 Jy ! Clean residual min=-0.003984 max=0.003573 Jy/beam ! Clean residual mean=-0.000020 rms=0.001041 Jy/beam ! Combined flux in latest and established models = 1.85421 Jy ! Adding 109 model components to the UV plane model. ! The established model now contains 222 components and 1.85421 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.797 mas, bmaj=7.154 mas, bpa=-13.14 degrees ! Estimated noise=1.06672 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000273137 Jy ! Component: 100 - total flux cleaned = -0.000274869 Jy ! Component: 150 - total flux cleaned = -0.000543942 Jy ! Component: 200 - total flux cleaned = -0.000411612 Jy ! Component: 250 - total flux cleaned = -0.000536151 Jy ! Component: 300 - total flux cleaned = -0.00053572 Jy ! Component: 350 - total flux cleaned = -0.000418638 Jy ! Component: 400 - total flux cleaned = -0.00041926 Jy ! Total flux subtracted in 400 components = -0.00041926 Jy ! Clean residual min=-0.004355 max=0.004093 Jy/beam ! Clean residual mean=-0.000012 rms=0.001022 Jy/beam ! Combined flux in latest and established models = 1.85379 Jy selfcal ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 247 components and 1.85379 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.153806Jy sigma=0.990263 ! Fit after self-cal, rms=1.153577Jy sigma=0.989052 wmodel /scr/pima/rdv62_uvs/J1608+1029_S_map.mod ! Writing 247 model components to file: /scr/pima/rdv62_uvs/J1608+1029_S_map.mod wobs /scr/pima/rdv62_uvs/J1608+1029_S_uvs.fits ! Applying 1591 buffered edits. ! Writing UV FITS file: /scr/pima/rdv62_uvs/J1608+1029_S_uvs.fits wwins /scr/pima/rdv62_uvs/J1608+1029_S_map.win ! wwins: Wrote 7 windows to /scr/pima/rdv62_uvs/J1608+1029_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.00116616 Jy ! Component: 100 - total flux cleaned = -0.00159242 Jy ! Component: 150 - total flux cleaned = -0.00192981 Jy ! Component: 200 - total flux cleaned = -0.0024518 Jy ! Component: 250 - total flux cleaned = -0.00283194 Jy ! Component: 300 - total flux cleaned = -0.00307794 Jy ! Component: 350 - total flux cleaned = -0.00331927 Jy ! Component: 400 - total flux cleaned = -0.00355615 Jy ! Component: 450 - total flux cleaned = -0.003614 Jy ! Component: 500 - total flux cleaned = -0.00378374 Jy ! Component: 550 - total flux cleaned = -0.00394928 Jy ! Component: 600 - total flux cleaned = -0.00411321 Jy ! Component: 650 - total flux cleaned = -0.00433042 Jy ! Component: 700 - total flux cleaned = -0.00454521 Jy ! Component: 750 - total flux cleaned = -0.00496813 Jy ! Component: 800 - total flux cleaned = -0.00522919 Jy ! Component: 850 - total flux cleaned = -0.00564218 Jy ! Component: 900 - total flux cleaned = -0.0059999 Jy ! Component: 950 - total flux cleaned = -0.00645471 Jy ! Component: 1000 - total flux cleaned = -0.00680504 Jy ! Component: 1050 - total flux cleaned = -0.00710304 Jy ! Component: 1100 - total flux cleaned = -0.00754545 Jy ! Component: 1150 - total flux cleaned = -0.00793541 Jy ! Component: 1200 - total flux cleaned = -0.0085142 Jy ! Component: 1250 - total flux cleaned = -0.00884878 Jy ! Component: 1300 - total flux cleaned = -0.00918124 Jy ! Component: 1350 - total flux cleaned = -0.00974598 Jy ! Component: 1400 - total flux cleaned = -0.0100732 Jy ! Component: 1450 - total flux cleaned = -0.010536 Jy ! Component: 1500 - total flux cleaned = -0.0108575 Jy ! Component: 1550 - total flux cleaned = -0.0112677 Jy ! Component: 1600 - total flux cleaned = -0.0116297 Jy ! Component: 1650 - total flux cleaned = -0.0120337 Jy ! Component: 1700 - total flux cleaned = -0.0123899 Jy ! Component: 1750 - total flux cleaned = -0.0126553 Jy ! Component: 1800 - total flux cleaned = -0.0128749 Jy ! Component: 1850 - total flux cleaned = -0.0132671 Jy ! Component: 1900 - total flux cleaned = -0.0134839 Jy ! Component: 1950 - total flux cleaned = -0.0138276 Jy ! Component: 2000 - total flux cleaned = -0.0143834 Jy ! Component: 2050 - total flux cleaned = -0.014595 Jy ! Component: 2100 - total flux cleaned = -0.014975 Jy ! Component: 2150 - total flux cleaned = -0.0153525 Jy ! Component: 2200 - total flux cleaned = -0.0156446 Jy ! Component: 2250 - total flux cleaned = -0.0159759 Jy ! Component: 2300 - total flux cleaned = -0.0162232 Jy ! Component: 2350 - total flux cleaned = -0.0164686 Jy ! Component: 2400 - total flux cleaned = -0.0167538 Jy ! Component: 2450 - total flux cleaned = -0.0171991 Jy ! Component: 2500 - total flux cleaned = -0.0173602 Jy ! Component: 2550 - total flux cleaned = -0.0175201 Jy ! Component: 2600 - total flux cleaned = -0.0179588 Jy ! Component: 2650 - total flux cleaned = -0.0180779 Jy ! Component: 2700 - total flux cleaned = -0.0183538 Jy ! Component: 2750 - total flux cleaned = -0.0186676 Jy ! Component: 2800 - total flux cleaned = -0.0188628 Jy ! Component: 2850 - total flux cleaned = -0.0189404 Jy ! Component: 2900 - total flux cleaned = -0.0192497 Jy ! Component: 2950 - total flux cleaned = -0.0195957 Jy ! Component: 3000 - total flux cleaned = -0.0198255 Jy ! Component: 3050 - total flux cleaned = -0.0201303 Jy ! Component: 3100 - total flux cleaned = -0.0202062 Jy ! Component: 3150 - total flux cleaned = -0.0204706 Jy ! Component: 3200 - total flux cleaned = -0.0206962 Jy ! Component: 3250 - total flux cleaned = -0.0208834 Jy ! Component: 3300 - total flux cleaned = -0.0209204 Jy ! Component: 3350 - total flux cleaned = -0.0209576 Jy ! Component: 3400 - total flux cleaned = -0.021327 Jy ! Component: 3450 - total flux cleaned = -0.0215476 Jy ! Component: 3500 - total flux cleaned = -0.0215839 Jy ! Component: 3550 - total flux cleaned = -0.0219122 Jy ! Component: 3600 - total flux cleaned = -0.0220938 Jy ! Component: 3650 - total flux cleaned = -0.0222746 Jy ! Component: 3700 - total flux cleaned = -0.0223108 Jy ! Component: 3750 - total flux cleaned = -0.0225619 Jy ! Component: 3800 - total flux cleaned = -0.0227405 Jy ! Component: 3850 - total flux cleaned = -0.022883 Jy ! Component: 3900 - total flux cleaned = -0.0230957 Jy ! Component: 3950 - total flux cleaned = -0.0231664 Jy ! Component: 4000 - total flux cleaned = -0.0233425 Jy ! Component: 4050 - total flux cleaned = -0.0235528 Jy ! Total flux subtracted in 4096 components = -0.0236574 Jy ! Clean residual min=-0.001888 max=0.001894 Jy/beam ! Clean residual mean=-0.000002 rms=0.000607 Jy/beam ! Combined flux in latest and established models = 1.83014 Jy keep ! Adding 227 model components to the UV plane model. ! The established model now contains 474 components and 1.83014 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv62_uvs/J1608+1029_S\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.797 x 7.154 at -13.14 degrees (North through East) ! Clean map min=-0.0054594 max=1.5642 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv62_uvs/J1608+1029_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= 1.5642 rms= 0.0006 Jy quit ! Quitting program ! Log file /scr/pima/rdv62_uvs/J1608+1029_S_uvs.log closed on Fri Nov 22 22:38:43 2013