! Started logfile: /scr/pima/rdv93_uvs/J0339-0146_X_uvs.log on Tue Oct 8 15:43:27 2013 obs /scr/pima/rdv93_uvs/J0339-0146_X_uva.fits ! Reading UV FITS file: /scr/pima/rdv93_uvs/J0339-0146_X_uva.fits ! AN table 1: 2 integrations on 45 of 45 possible baselines. ! AN table 2: 11 integrations on 1 of 1 possible baselines. ! AN table 3: 9 integrations on 15 of 15 possible baselines. ! AN table 4: 4 integrations on 36 of 36 possible baselines. ! AN table 5: 4 integrations on 21 of 21 possible baselines. ! AN table 6: 8 integrations on 55 of 55 possible baselines. ! Apparent sampling: 0.747788 visibilities/baseline/integration-bin. ! Found source: J0339-0146 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 8.40599e+09 1.6e+07 1 1.6e+07 ! 02 2 8.47599e+09 1.6e+07 1 1.6e+07 ! 03 3 8.79099e+09 1.6e+07 1 1.6e+07 ! 04 4 8.89599e+09 1.6e+07 1 1.6e+07 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 2704 visibilities. select RR ! Selecting polarization: RR, channels: 1..4 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 print "field_size =",field_size ! field_size = 1024 print "field_cell =",field_cell ! field_cell = 0.2 print "taper_size =",taper_size ! taper_size = 70 ![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv93_uvs/J0339-0146_X] integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 800000 float soltime2 soltime2 = 15 float soltime3 soltime3 = 1 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float old_peak float new_peak float flux_cutoff selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:KP self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:MK self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:OV self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 1:SC self-cal status=correctable weight=1 ! 2:FT self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 3:FT self-cal status=correctable weight=1 ! 3:HN self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 3:NY self-cal status=correctable weight=1 ! 3:SC self-cal status=correctable weight=1 ! 3:WZ self-cal status=correctable weight=1 ! 4:FD self-cal status=correctable weight=1 ! 4:FT 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:NL 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 ! 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:PT 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:FT self-cal status=correctable weight=1 ! 6:HN 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.200x0.200 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 22 telescope corrections were flagged in sub-array 2. ! A total of 14 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 48 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 22 telescope corrections were flagged in sub-array 2. ! A total of 14 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 48 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 22 telescope corrections were flagged in sub-array 2. ! A total of 14 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 48 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 22 telescope corrections were flagged in sub-array 2. ! A total of 14 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 48 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=2.498394Jy sigma=117.494440 ! Fit after self-cal, rms=1.371857Jy sigma=70.085641 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. print "=== Flag out IVS stations" ! === Flag out IVS stations integer lun string cfg_file lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_flag_process.csh /tmp/ff.ff" ![@/tmp/ff.ff \] flag 2:FT ! Flagging baselines of antenna 2:FT in the currently selected channels. flag 2:WZ ! Flagging baselines of antenna 2:WZ in the currently selected channels. flag 3:FT ! Flagging baselines of antenna 3:FT 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 4:FT ! Flagging baselines of antenna 4:FT in the currently selected channels. flag 4:WZ ! Flagging baselines of antenna 4:WZ in the currently selected channels. flag 5:WZ ! Flagging baselines of antenna 5:WZ in the currently selected channels. flag 6:FT ! Flagging baselines of antenna 6:FT in the currently selected channels. ![Exited script file: /tmp/ff.ff] uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.9223 mas, bmaj=1.912 mas, bpa=-1.181 degrees ! Estimated noise=0.48796 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.90193 Jy ! Component: 100 - total flux cleaned = 2.34175 Jy ! Total flux subtracted in 100 components = 2.34175 Jy ! Clean residual min=-0.087689 max=0.208829 Jy/beam ! Clean residual mean=0.000343 rms=0.021980 Jy/beam ! Combined flux in latest and established models = 2.34176 Jy ! Performing phase self-cal ! Adding 5 model components to the UV plane model. ! The established model now contains 5 components and 2.34176 Jy ! ! Correcting IF 1. ! A total of 16 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 16 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 16 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 16 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.475577Jy sigma=25.089021 ! Fit after self-cal, rms=0.475319Jy sigma=25.084029 ! Inverting map and beam ! Estimated beam: bmin=0.9223 mas, bmaj=1.912 mas, bpa=-1.181 degrees ! Estimated noise=0.48796 mJy/beam. ! Added new window around map position (3.4, 2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.215648 Jy ! Component: 100 - total flux cleaned = 0.351602 Jy ! Total flux subtracted in 100 components = 0.351602 Jy ! Clean residual min=-0.076358 max=0.089820 Jy/beam ! Clean residual mean=0.000209 rms=0.015314 Jy/beam ! Combined flux in latest and established models = 2.69336 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 25 components and 2.69336 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.354459Jy sigma=17.225725 ! Fit after self-cal, rms=0.338462Jy sigma=16.740924 ! Inverting map ! Added new window around map position (2.2, 1.8). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.113003 Jy ! Component: 100 - total flux cleaned = 0.193391 Jy ! Total flux subtracted in 100 components = 0.193391 Jy ! Clean residual min=-0.059936 max=0.059949 Jy/beam ! Clean residual mean=0.000136 rms=0.011809 Jy/beam ! Combined flux in latest and established models = 2.88675 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 49 components and 2.88675 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.289136Jy sigma=13.293946 ! Fit after self-cal, rms=0.282034Jy sigma=13.060110 ! Inverting map print "=== Finished uniform weighting clean ===" ! === Finished uniform weighting clean === uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.117 mas, bmaj=2.226 mas, bpa=-0.5144 degrees ! Estimated noise=0.379107 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.101738 Jy ! Component: 100 - total flux cleaned = 0.165877 Jy ! Component: 150 - total flux cleaned = 0.209813 Jy ! Component: 200 - total flux cleaned = 0.240767 Jy ! Total flux subtracted in 200 components = 0.240767 Jy ! Clean residual min=-0.046044 max=0.052092 Jy/beam ! Clean residual mean=0.000081 rms=0.009760 Jy/beam ! Combined flux in latest and established models = 3.12752 Jy ! Performing phase self-cal ! Adding 52 model components to the UV plane model. ! The established model now contains 90 components and 3.12752 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.247510Jy sigma=10.297101 ! Fit after self-cal, rms=0.240031Jy sigma=9.964941 ! Inverting map print "=== Finished natural weighting clean for VLBI sites ===" ! === Finished natural weighting clean for VLBI sites === gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.96 FD 0.97 HN 0.97 KP 0.98 ! LA 0.99 MK 0.92 NL 0.96 OV 0.99 ! PT 0.95 SC 1.05 ! ! Telescope amplitude corrections in sub-array 2: ! FT 0.00* WZ 0.00* ! ! A total of 24 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! FT 0.00* HN 0.00* NL 0.00* NY 0.00* ! SC 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! FD 0.99 FT 1.00* HN 0.98 KP 1.06 ! LA 1.06 NL 1.00 PT 1.02 SC 0.97 ! WZ 1.00* ! ! Telescope amplitude corrections in sub-array 5: ! FD 0.97 HN 1.03 KP 1.02 LA 1.03 ! NL 1.01 PT 0.99 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 6: ! BR 0.99 FD 0.98 FT 1.00* HN 0.98 ! KP 1.02 LA 0.99 MK 1.15 NL 0.96 ! OV 0.99 PT 0.97 SC 0.99 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.97 FD 0.97 HN 0.99 KP 0.93 ! LA 0.95 MK 0.94 NL 0.97 OV 0.97 ! PT 0.95 SC 1.07 ! ! Telescope amplitude corrections in sub-array 2: ! FT 0.00* WZ 0.00* ! ! A total of 24 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! FT 0.00* HN 0.00* NL 0.00* NY 0.00* ! SC 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! FD 1.00 FT 1.00* HN 0.99 KP 1.00 ! LA 1.04 NL 1.00 PT 1.03 SC 1.00 ! WZ 1.00* ! ! Telescope amplitude corrections in sub-array 5: ! FD 0.98 HN 1.03 KP 0.96 LA 1.00 ! NL 1.02 PT 0.99 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 6: ! BR 1.00 FD 0.98 FT 1.00* HN 0.99 ! KP 0.96 LA 0.96 MK 1.18 NL 0.97 ! OV 0.98 PT 0.98 SC 1.01 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.02 FD 0.97 HN 1.01 KP 0.84 ! LA 1.02 MK 0.94 NL 1.01 OV 1.03 ! PT 0.96 SC 1.06 ! ! Telescope amplitude corrections in sub-array 2: ! FT 0.00* WZ 0.00* ! ! A total of 24 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! FT 0.00* HN 0.00* NL 0.00* NY 0.00* ! SC 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! FD 1.00 FT 1.00* HN 1.01 KP 0.89 ! LA 1.10 NL 1.07 PT 1.02 SC 0.97 ! WZ 1.00* ! ! Telescope amplitude corrections in sub-array 5: ! FD 0.98 HN 1.05 KP 0.87 LA 1.08 ! NL 1.07 PT 0.98 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 6: ! BR 1.05 FD 0.97 FT 1.00* HN 1.01 ! KP 0.87 LA 1.02 MK 1.13 NL 1.02 ! OV 1.04 PT 0.98 SC 1.00 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00 FD 1.03 HN 1.00 KP 0.89 ! LA 1.01 MK 0.98 NL 1.02 OV 1.01 ! PT 1.00 SC 1.07 ! ! Telescope amplitude corrections in sub-array 2: ! FT 0.00* WZ 0.00* ! ! A total of 24 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! FT 0.00* HN 0.00* NL 0.00* NY 0.00* ! SC 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! FD 1.06 FT 1.00* HN 1.01 KP 0.96 ! LA 1.08 NL 1.08 PT 1.07 SC 1.01 ! WZ 1.00* ! ! Telescope amplitude corrections in sub-array 5: ! FD 1.03 HN 1.05 KP 0.93 LA 1.05 ! NL 1.09 PT 1.03 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 6: ! BR 1.03 FD 1.04 FT 1.00* HN 1.02 ! KP 0.91 LA 1.01 MK 1.20 NL 1.05 ! OV 1.01 PT 1.02 SC 1.03 ! ! ! Fit before self-cal, rms=0.240031Jy sigma=9.964941 ! Fit after self-cal, rms=0.120147Jy sigma=5.074556 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.111 mas, bmaj=2.336 mas, bpa=-0.08765 degrees ! Estimated noise=0.37943 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 (4.4, 2.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0323056 Jy ! Component: 100 - total flux cleaned = 0.0437717 Jy ! Component: 150 - total flux cleaned = 0.0471096 Jy ! Component: 200 - total flux cleaned = 0.0493397 Jy ! Total flux subtracted in 200 components = 0.0493397 Jy ! Clean residual min=-0.011559 max=0.017502 Jy/beam ! Clean residual mean=-0.000005 rms=0.003164 Jy/beam ! Combined flux in latest and established models = 3.17686 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 107 components and 3.17686 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.110323Jy sigma=4.477760 ! Fit after self-cal, rms=0.107832Jy sigma=4.374275 ! Inverting map ! Added new window around map position (-2, 10.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0181035 Jy ! Component: 100 - total flux cleaned = 0.0211344 Jy ! Component: 150 - total flux cleaned = 0.021974 Jy ! Component: 200 - total flux cleaned = 0.0219802 Jy ! Total flux subtracted in 200 components = 0.0219802 Jy ! Clean residual min=-0.010217 max=0.011188 Jy/beam ! Clean residual mean=-0.000010 rms=0.002462 Jy/beam ! Combined flux in latest and established models = 3.19884 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 122 components and 3.19884 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.102722Jy sigma=4.124426 ! Fit after self-cal, rms=0.102231Jy sigma=4.105618 ! 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.102231Jy sigma=4.105618 ! Fit after self-cal, rms=0.101334Jy sigma=4.053443 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.112 mas, bmaj=2.337 mas, bpa=-0.07648 degrees ! Estimated noise=0.37968 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0118075 Jy ! Component: 100 - total flux cleaned = 0.0118467 Jy ! Component: 150 - total flux cleaned = 0.0121362 Jy ! Component: 200 - total flux cleaned = 0.0116028 Jy ! Total flux subtracted in 200 components = 0.0116028 Jy ! Clean residual min=-0.007876 max=0.009615 Jy/beam ! Clean residual mean=-0.000001 rms=0.002041 Jy/beam ! Combined flux in latest and established models = 3.21044 Jy selfcal ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 131 components and 3.21044 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.098859Jy sigma=3.934826 ! Fit after self-cal, rms=0.098442Jy sigma=3.919045 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 (-6.8, -19.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00811133 Jy ! Component: 100 - total flux cleaned = 0.0102553 Jy ! Component: 150 - total flux cleaned = 0.0107148 Jy ! Component: 200 - total flux cleaned = 0.0111437 Jy ! Total flux subtracted in 200 components = 0.0111437 Jy ! Clean residual min=-0.008116 max=0.008578 Jy/beam ! Clean residual mean=-0.000006 rms=0.001871 Jy/beam ! Combined flux in latest and established models = 3.22158 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 141 components and 3.22158 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.096841Jy sigma=3.849587 ! Fit after self-cal, rms=0.096687Jy sigma=3.844831 ! 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.096687Jy sigma=3.844831 ! Fit after self-cal, rms=0.096213Jy sigma=3.829692 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.113 mas, bmaj=2.337 mas, bpa=-0.05479 degrees ! Estimated noise=0.379696 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00612953 Jy ! Component: 100 - total flux cleaned = 0.00748711 Jy ! Component: 150 - total flux cleaned = 0.0078897 Jy ! Component: 200 - total flux cleaned = 0.00770681 Jy ! Total flux subtracted in 200 components = 0.00770681 Jy ! Clean residual min=-0.006952 max=0.008009 Jy/beam ! Clean residual mean=-0.000001 rms=0.001709 Jy/beam ! Combined flux in latest and established models = 3.22929 Jy selfcal ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 151 components and 3.22929 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.095080Jy sigma=3.783041 ! Fit after self-cal, rms=0.094947Jy sigma=3.778902 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 (-6.8, -36.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00706991 Jy ! Component: 100 - total flux cleaned = 0.00889374 Jy ! Component: 150 - total flux cleaned = 0.00994784 Jy ! Component: 200 - total flux cleaned = 0.0102767 Jy ! Total flux subtracted in 200 components = 0.0102767 Jy ! Clean residual min=-0.006677 max=0.006805 Jy/beam ! Clean residual mean=0.000001 rms=0.001571 Jy/beam ! Combined flux in latest and established models = 3.23957 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 161 components and 3.23957 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.093936Jy sigma=3.738381 ! Fit after self-cal, rms=0.093904Jy sigma=3.736474 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime2 ! Performing amp+phase self-cal over 15 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.093904Jy sigma=3.736474 ! Fit after self-cal, rms=0.093661Jy sigma=3.730278 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.368 mas, bmaj=3.585 mas, bpa=-6.151 degrees ! Estimated noise=0.416224 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00508348 Jy ! Component: 100 - total flux cleaned = 0.00598292 Jy ! Component: 150 - total flux cleaned = 0.00624599 Jy ! Component: 200 - total flux cleaned = 0.00659697 Jy ! Total flux subtracted in 200 components = 0.00659697 Jy ! Clean residual min=-0.007242 max=0.007029 Jy/beam ! Clean residual mean=0.000003 rms=0.001743 Jy/beam ! Combined flux in latest and established models = 3.24616 Jy selfcal ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 180 components and 3.24616 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.093363Jy sigma=3.710582 ! Fit after self-cal, rms=0.093315Jy sigma=3.708921 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.000900833 Jy ! Component: 100 - total flux cleaned = 0.000792714 Jy ! Component: 150 - total flux cleaned = 0.000695802 Jy ! Component: 200 - total flux cleaned = 0.000792488 Jy ! Total flux subtracted in 200 components = 0.000792488 Jy ! Clean residual min=-0.006895 max=0.007228 Jy/beam ! Clean residual mean=0.000002 rms=0.001691 Jy/beam ! Combined flux in latest and established models = 3.24696 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 193 components and 3.24696 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.093118Jy sigma=3.699164 ! Fit after self-cal, rms=0.093081Jy sigma=3.698214 ! Inverting map ! Added new window around map position (-22, 37.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00652104 Jy ! Component: 100 - total flux cleaned = 0.00846795 Jy ! Component: 150 - total flux cleaned = 0.00835756 Jy ! Component: 200 - total flux cleaned = 0.00790679 Jy ! Total flux subtracted in 200 components = 0.00790679 Jy ! Clean residual min=-0.005867 max=0.005812 Jy/beam ! Clean residual mean=0.000003 rms=0.001560 Jy/beam ! Combined flux in latest and established models = 3.25486 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 222 components and 3.25486 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.092520Jy sigma=3.672369 ! Fit after self-cal, rms=0.092489Jy sigma=3.670945 ! 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 1:BR,true ! 1:BR self-cal status=fixed weight=1 selfant 1:FD,true ! 1:FD self-cal status=fixed weight=1 selfant 1:HN,true ! 1:HN self-cal status=fixed weight=1 selfant 1:KP,true ! 1:KP self-cal status=fixed weight=1 selfant 1:LA,true ! 1:LA self-cal status=fixed weight=1 selfant 1:MK,true ! 1:MK self-cal status=fixed weight=1 selfant 1:NL,true ! 1:NL self-cal status=fixed weight=1 selfant 1:OV,true ! 1:OV self-cal status=fixed weight=1 selfant 1:PT,true ! 1:PT self-cal status=fixed weight=1 selfant 1:SC,true ! 1:SC self-cal status=fixed weight=1 selfant 3:HN,true ! 3:HN self-cal status=fixed weight=1 selfant 3:NL,true ! 3:NL self-cal status=fixed weight=1 selfant 3:SC,true ! 3:SC 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:NL,true ! 4:NL self-cal status=fixed weight=1 selfant 4:PT,true ! 4:PT self-cal status=fixed weight=1 selfant 4:SC,true ! 4:SC self-cal status=fixed weight=1 selfant 5:FD,true ! 5:FD self-cal status=fixed weight=1 selfant 5:HN,true ! 5:HN self-cal status=fixed weight=1 selfant 5:KP,true ! 5:KP self-cal status=fixed weight=1 selfant 5:LA,true ! 5:LA self-cal status=fixed weight=1 selfant 5:NL,true ! 5:NL self-cal status=fixed weight=1 selfant 5:PT,true ! 5:PT self-cal status=fixed weight=1 selfant 6:BR,true ! 6:BR self-cal status=fixed weight=1 selfant 6: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 22 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 22 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 22 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 22 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.277378Jy sigma=7.299666 ! Fit after self-cal, rms=0.183466Jy sigma=4.315677 selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:KP self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:MK self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:OV self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 1:SC self-cal status=correctable weight=1 ! 2:FT self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 3:FT self-cal status=correctable weight=1 ! 3:HN self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 3:NY self-cal status=correctable weight=1 ! 3:SC self-cal status=correctable weight=1 ! 3:WZ self-cal status=correctable weight=1 ! 4:FD self-cal status=correctable weight=1 ! 4:FT 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:NL 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 ! 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:PT 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:FT self-cal status=correctable weight=1 ! 6:HN 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 3 telescope corrections were flagged in sub-array 3. ! A total of 40 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 40 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 40 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 3 telescope corrections were flagged in sub-array 3. ! A total of 40 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.183466Jy sigma=4.315677 ! Fit after self-cal, rms=0.153207Jy sigma=3.933982 uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.7158 mas, bmaj=1.877 mas, bpa=-7.608 degrees ! Estimated noise=0.48383 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.00523734 Jy ! Component: 100 - total flux cleaned = 0.00579929 Jy ! Component: 150 - total flux cleaned = 0.00743636 Jy ! Component: 200 - total flux cleaned = 0.00843538 Jy ! Total flux subtracted in 200 components = 0.00843538 Jy ! Clean residual min=-0.012583 max=0.012982 Jy/beam ! Clean residual mean=-0.000008 rms=0.002640 Jy/beam ! Combined flux in latest and established models = 3.2633 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 244 components and 3.2633 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.147944Jy sigma=3.859200 ! Fit after self-cal, rms=0.146545Jy sigma=3.836392 ! 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=0.8351 mas, bmaj=2.149 mas, bpa=-7.518 degrees ! Estimated noise=0.377874 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.000347477 Jy ! Component: 100 - total flux cleaned = -0.00102982 Jy ! Component: 150 - total flux cleaned = -0.00135805 Jy ! Component: 200 - total flux cleaned = -0.0019614 Jy ! Component: 250 - total flux cleaned = -0.00224898 Jy ! Component: 300 - total flux cleaned = -0.00279638 Jy ! Component: 350 - total flux cleaned = -0.00332831 Jy ! Component: 400 - total flux cleaned = -0.0038356 Jy ! Total flux subtracted in 400 components = -0.0038356 Jy ! Clean residual min=-0.008587 max=0.007796 Jy/beam ! Clean residual mean=-0.000002 rms=0.001893 Jy/beam ! Combined flux in latest and established models = 3.25946 Jy ! Performing phase self-cal ! Adding 44 model components to the UV plane model. ! The established model now contains 271 components and 3.25946 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.142332Jy sigma=3.758634 ! Fit after self-cal, rms=0.142254Jy sigma=3.756297 ! 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.000564374 Jy ! Component: 100 - total flux cleaned = 8.48584e-05 Jy ! Component: 150 - total flux cleaned = -0.0001483 Jy ! Component: 200 - total flux cleaned = -0.000379409 Jy ! Component: 250 - total flux cleaned = -0.00103776 Jy ! Component: 300 - total flux cleaned = -0.0008168 Jy ! Component: 350 - total flux cleaned = -0.0010284 Jy ! Component: 400 - total flux cleaned = -0.00123059 Jy ! Total flux subtracted in 400 components = -0.00123059 Jy ! Clean residual min=-0.008118 max=0.007453 Jy/beam ! Clean residual mean=-0.000003 rms=0.001827 Jy/beam ! Combined flux in latest and established models = 3.25823 Jy ! Performing phase self-cal ! Adding 60 model components to the UV plane model. ! The established model now contains 293 components and 3.25823 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.139776Jy sigma=3.716319 ! Fit after self-cal, rms=0.139755Jy sigma=3.715387 ! 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. 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 22 telescope corrections were flagged in sub-array 2. ! A total of 14 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 48 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 22 telescope corrections were flagged in sub-array 2. ! A total of 14 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 48 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 22 telescope corrections were flagged in sub-array 2. ! A total of 14 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 48 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 22 telescope corrections were flagged in sub-array 2. ! A total of 14 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 48 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.204114Jy sigma=5.045654 ! Fit after self-cal, rms=0.139265Jy sigma=3.704749 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.8351 mas, bmaj=2.149 mas, bpa=-7.505 degrees ! Estimated noise=0.37781 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00221881 Jy ! Component: 100 - total flux cleaned = 0.00221788 Jy ! Component: 150 - total flux cleaned = 0.00200811 Jy ! Component: 200 - total flux cleaned = 0.00160942 Jy ! Component: 250 - total flux cleaned = 0.00141469 Jy ! Component: 300 - total flux cleaned = 0.00103847 Jy ! Component: 350 - total flux cleaned = 0.000855001 Jy ! Component: 400 - total flux cleaned = 0.000494636 Jy ! Total flux subtracted in 400 components = 0.000494636 Jy ! Clean residual min=-0.007851 max=0.006997 Jy/beam ! Clean residual mean=-0.000002 rms=0.001750 Jy/beam ! Combined flux in latest and established models = 3.25873 Jy selfcal ! Performing phase self-cal ! Adding 65 model components to the UV plane model. ! The established model now contains 315 components and 3.25872 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.137118Jy sigma=3.673238 ! Fit after self-cal, rms=0.137093Jy sigma=3.672097 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.00097587 Jy ! Component: 100 - total flux cleaned = 0.000622153 Jy ! Component: 150 - total flux cleaned = 0.000107792 Jy ! Component: 200 - total flux cleaned = 0.000111254 Jy ! Component: 250 - total flux cleaned = -5.18595e-05 Jy ! Component: 300 - total flux cleaned = -0.000212616 Jy ! Component: 350 - total flux cleaned = -0.000532877 Jy ! Component: 400 - total flux cleaned = -0.000686034 Jy ! Total flux subtracted in 400 components = -0.000686034 Jy ! Clean residual min=-0.007777 max=0.006916 Jy/beam ! Clean residual mean=-0.000002 rms=0.001717 Jy/beam ! Combined flux in latest and established models = 3.25804 Jy ! Performing phase self-cal ! Adding 79 model components to the UV plane model. ! The established model now contains 332 components and 3.25804 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.135544Jy sigma=3.652301 ! Fit after self-cal, rms=0.135543Jy sigma=3.651677 ! 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.135543Jy sigma=3.651677 ! Fit after self-cal, rms=0.112199Jy sigma=3.471128 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.8352 mas, bmaj=2.15 mas, bpa=-7.482 degrees ! Estimated noise=0.377974 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0009895 Jy ! Component: 100 - total flux cleaned = 0.00138669 Jy ! Component: 150 - total flux cleaned = 0.0020878 Jy ! Component: 200 - total flux cleaned = 0.00257155 Jy ! Component: 250 - total flux cleaned = 0.00303639 Jy ! Component: 300 - total flux cleaned = 0.00333064 Jy ! Component: 350 - total flux cleaned = 0.00347237 Jy ! Component: 400 - total flux cleaned = 0.00375188 Jy ! Total flux subtracted in 400 components = 0.00375188 Jy ! Clean residual min=-0.006490 max=0.005641 Jy/beam ! Clean residual mean=0.000004 rms=0.001196 Jy/beam ! Combined flux in latest and established models = 3.26179 Jy selfcal ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 366 components and 3.26179 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.110780Jy sigma=3.437502 ! Fit after self-cal, rms=0.110666Jy sigma=3.436066 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.000815558 Jy ! Component: 100 - total flux cleaned = 0.00122513 Jy ! Component: 150 - total flux cleaned = 0.00135668 Jy ! Component: 200 - total flux cleaned = 0.00161169 Jy ! Component: 250 - total flux cleaned = 0.00198096 Jy ! Component: 300 - total flux cleaned = 0.00210232 Jy ! Component: 350 - total flux cleaned = 0.00245168 Jy ! Component: 400 - total flux cleaned = 0.00256548 Jy ! Total flux subtracted in 400 components = 0.00256548 Jy ! Clean residual min=-0.005986 max=0.005507 Jy/beam ! Clean residual mean=0.000004 rms=0.001150 Jy/beam ! Combined flux in latest and established models = 3.26436 Jy ! Performing phase self-cal ! Adding 64 model components to the UV plane model. ! The established model now contains 386 components and 3.26436 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.109859Jy sigma=3.418241 ! Fit after self-cal, rms=0.109812Jy sigma=3.417777 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, 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.109812Jy sigma=3.417777 ! Fit after self-cal, rms=0.106988Jy sigma=3.292400 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.369 mas, bmaj=3.345 mas, bpa=-6.518 degrees ! Estimated noise=0.415344 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00136687 Jy ! Component: 100 - total flux cleaned = 0.00194575 Jy ! Component: 150 - total flux cleaned = 0.00251753 Jy ! Component: 200 - total flux cleaned = 0.0030363 Jy ! Component: 250 - total flux cleaned = 0.00327765 Jy ! Component: 300 - total flux cleaned = 0.00356386 Jy ! Component: 350 - total flux cleaned = 0.00388854 Jy ! Component: 400 - total flux cleaned = 0.00409489 Jy ! Total flux subtracted in 400 components = 0.00409489 Jy ! Clean residual min=-0.004182 max=0.004402 Jy/beam ! Clean residual mean=0.000004 rms=0.001128 Jy/beam ! Combined flux in latest and established models = 3.26845 Jy selfcal ! Performing phase self-cal ! Adding 46 model components to the UV plane model. ! The established model now contains 424 components and 3.26845 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.106934Jy sigma=3.286084 ! Fit after self-cal, rms=0.106909Jy sigma=3.285669 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 (-42.4, -48.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000218835 Jy ! Component: 100 - total flux cleaned = -2.44684e-05 Jy ! Component: 150 - total flux cleaned = -0.000202369 Jy ! Component: 200 - total flux cleaned = -0.000583659 Jy ! Component: 250 - total flux cleaned = -0.000931015 Jy ! Component: 300 - total flux cleaned = -0.00112655 Jy ! Component: 350 - total flux cleaned = -0.00137372 Jy ! Component: 400 - total flux cleaned = -0.0014888 Jy ! Total flux subtracted in 400 components = -0.0014888 Jy ! Clean residual min=-0.003420 max=0.004024 Jy/beam ! Clean residual mean=0.000009 rms=0.001009 Jy/beam ! Combined flux in latest and established models = 3.26696 Jy ! Performing phase self-cal ! Adding 53 model components to the UV plane model. ! The established model now contains 472 components and 3.26696 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.106664Jy sigma=3.271292 ! Fit after self-cal, rms=0.106663Jy sigma=3.270580 ! Inverting map ! Added new window around map position (-22.6, -37.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00370655 Jy ! Component: 100 - total flux cleaned = 0.00615477 Jy ! Component: 150 - total flux cleaned = 0.0082173 Jy ! Component: 200 - total flux cleaned = 0.00955447 Jy ! Component: 250 - total flux cleaned = 0.0105156 Jy ! Component: 300 - total flux cleaned = 0.0112256 Jy ! Component: 350 - total flux cleaned = 0.0119569 Jy ! Component: 400 - total flux cleaned = 0.0124234 Jy ! Total flux subtracted in 400 components = 0.0124234 Jy ! Clean residual min=-0.003346 max=0.003656 Jy/beam ! Clean residual mean=0.000007 rms=0.000947 Jy/beam ! Combined flux in latest and established models = 3.27939 Jy ! Performing phase self-cal ! Adding 61 model components to the UV plane model. ! The established model now contains 523 components and 3.27939 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.106437Jy sigma=3.257473 ! Fit after self-cal, rms=0.106403Jy sigma=3.256529 ! Inverting map ! Added new window around map position (49.4, 38.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00229666 Jy ! Component: 100 - total flux cleaned = 0.00279042 Jy ! Component: 150 - total flux cleaned = 0.00304343 Jy ! Component: 200 - total flux cleaned = 0.00342477 Jy ! Component: 250 - total flux cleaned = 0.00370609 Jy ! Component: 300 - total flux cleaned = 0.00390736 Jy ! Component: 350 - total flux cleaned = 0.0041613 Jy ! Component: 400 - total flux cleaned = 0.00440246 Jy ! Total flux subtracted in 400 components = 0.00440246 Jy ! Clean residual min=-0.003499 max=0.003525 Jy/beam ! Clean residual mean=0.000002 rms=0.000885 Jy/beam ! Combined flux in latest and established models = 3.28379 Jy ! Performing phase self-cal ! Adding 64 model components to the UV plane model. ! The established model now contains 578 components and 3.28379 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.106136Jy sigma=3.243358 ! Fit after self-cal, rms=0.106113Jy sigma=3.242678 ! Inverting map ! Added new window around map position (-8.2, 30.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00320039 Jy ! Component: 100 - total flux cleaned = 0.00413381 Jy ! Component: 150 - total flux cleaned = 0.00529702 Jy ! Component: 200 - total flux cleaned = 0.00644645 Jy ! Component: 250 - total flux cleaned = 0.00753283 Jy ! Component: 300 - total flux cleaned = 0.00878954 Jy ! Component: 350 - total flux cleaned = 0.00998601 Jy ! Component: 400 - total flux cleaned = 0.0110663 Jy ! Total flux subtracted in 400 components = 0.0110663 Jy ! Clean residual min=-0.003123 max=0.002928 Jy/beam ! Clean residual mean=0.000007 rms=0.000792 Jy/beam ! Combined flux in latest and established models = 3.29485 Jy ! Performing phase self-cal ! Adding 68 model components to the UV plane model. ! The established model now contains 629 components and 3.29485 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.105874Jy sigma=3.230713 ! Fit after self-cal, rms=0.105876Jy sigma=3.230320 ! 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=0.8353 mas, bmaj=2.15 mas, bpa=-7.478 degrees ! Estimated noise=0.377985 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.8353 x 2.15 at -7.478 degrees (North through East) ! Clean map min=-0.0091393 max=2.3476 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.105876Jy sigma=3.230320 ! Fit after self-cal, rms=0.094499Jy sigma=2.740320 ! Inverting map and beam ! Estimated beam: bmin=0.8352 mas, bmaj=2.151 mas, bpa=-7.487 degrees ! Estimated noise=0.377999 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000132285 Jy ! Component: 100 - total flux cleaned = 0.000733063 Jy ! Component: 150 - total flux cleaned = 0.0020909 Jy ! Component: 200 - total flux cleaned = 0.0030067 Jy ! Component: 250 - total flux cleaned = 0.00414854 Jy ! Component: 300 - total flux cleaned = 0.00464472 Jy ! Component: 350 - total flux cleaned = 0.00536745 Jy ! Component: 400 - total flux cleaned = 0.00595514 Jy ! Total flux subtracted in 400 components = 0.00595514 Jy ! Clean residual min=-0.004385 max=0.004209 Jy/beam ! Clean residual mean=0.000005 rms=0.000901 Jy/beam ! Combined flux in latest and established models = 3.30081 Jy ! Performing phase self-cal ! Adding 78 model components to the UV plane model. ! The established model now contains 690 components and 3.30081 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.093456Jy sigma=2.716731 ! Fit after self-cal, rms=0.093445Jy sigma=2.716226 ! Inverting map ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000601747 Jy ! Component: 100 - total flux cleaned = 0.00117737 Jy ! Component: 150 - total flux cleaned = 0.00208703 Jy ! Component: 200 - total flux cleaned = 0.00241912 Jy ! Component: 250 - total flux cleaned = 0.00318673 Jy ! Component: 300 - total flux cleaned = 0.0036188 Jy ! Component: 350 - total flux cleaned = 0.00404706 Jy ! Component: 400 - total flux cleaned = 0.00478276 Jy ! Total flux subtracted in 400 components = 0.00478276 Jy ! Clean residual min=-0.004187 max=0.004132 Jy/beam ! Clean residual mean=0.000005 rms=0.000867 Jy/beam ! Combined flux in latest and established models = 3.30559 Jy ! Performing phase self-cal ! Adding 119 model components to the UV plane model. ! The established model now contains 736 components and 3.30559 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.092688Jy sigma=2.699873 ! Fit after self-cal, rms=0.092683Jy sigma=2.699626 ! Inverting map ! Added new window around map position (13.6, 1.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00165558 Jy ! Component: 100 - total flux cleaned = 0.00220702 Jy ! Component: 150 - total flux cleaned = 0.00316581 Jy ! Component: 200 - total flux cleaned = 0.00368117 Jy ! Component: 250 - total flux cleaned = 0.00388516 Jy ! Component: 300 - total flux cleaned = 0.00438671 Jy ! Component: 350 - total flux cleaned = 0.00438791 Jy ! Component: 400 - total flux cleaned = 0.00438845 Jy ! Total flux subtracted in 400 components = 0.00438845 Jy ! Clean residual min=-0.003783 max=0.003717 Jy/beam ! Clean residual mean=0.000003 rms=0.000833 Jy/beam ! Combined flux in latest and established models = 3.30998 Jy ! Performing phase self-cal ! Adding 110 model components to the UV plane model. ! The established model now contains 776 components and 3.30998 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.091959Jy sigma=2.682815 ! Fit after self-cal, rms=0.091941Jy sigma=2.682518 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.091941Jy sigma=2.682518 ! Fit after self-cal, rms=0.091816Jy sigma=2.680583 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=0.716 mas, bmaj=1.88 mas, bpa=-7.579 degrees ! Estimated noise=0.48415 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 1.838 Jy ! Component: 100 - total flux cleaned = 2.29069 Jy ! Component: 150 - total flux cleaned = 2.52142 Jy ! Component: 200 - total flux cleaned = 2.67314 Jy ! Component: 250 - total flux cleaned = 2.78486 Jy ! Component: 300 - total flux cleaned = 2.87048 Jy ! Component: 350 - total flux cleaned = 2.938 Jy ! Component: 400 - total flux cleaned = 2.99192 Jy ! Total flux subtracted in 400 components = 2.99192 Jy ! Clean residual min=-0.016361 max=0.032080 Jy/beam ! Clean residual mean=0.000082 rms=0.004572 Jy/beam ! Combined flux in latest and established models = 2.99192 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 68 model components to the UV plane model. ! The established model now contains 68 components and 2.99192 Jy ! Inverting map and beam ! Estimated beam: bmin=0.8352 mas, bmaj=2.151 mas, bpa=-7.488 degrees ! Estimated noise=0.378057 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.0702573 Jy ! Component: 100 - total flux cleaned = 0.11809 Jy ! Component: 150 - total flux cleaned = 0.15169 Jy ! Component: 200 - total flux cleaned = 0.175785 Jy ! Component: 250 - total flux cleaned = 0.193522 Jy ! Component: 300 - total flux cleaned = 0.207014 Jy ! Component: 350 - total flux cleaned = 0.217604 Jy ! Component: 400 - total flux cleaned = 0.224423 Jy ! Total flux subtracted in 400 components = 0.224423 Jy ! Clean residual min=-0.006894 max=0.008316 Jy/beam ! Clean residual mean=0.000016 rms=0.001779 Jy/beam ! Combined flux in latest and established models = 3.21634 Jy ! Adding 109 model components to the UV plane model. ! The established model now contains 151 components and 3.21634 Jy ! Inverting map ! Added new window around map position (0.6, 3.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00963451 Jy ! Component: 100 - total flux cleaned = 0.0156228 Jy ! Component: 150 - total flux cleaned = 0.0199302 Jy ! Component: 200 - total flux cleaned = 0.0233822 Jy ! Component: 250 - total flux cleaned = 0.0268524 Jy ! Component: 300 - total flux cleaned = 0.0301149 Jy ! Component: 350 - total flux cleaned = 0.0332337 Jy ! Component: 400 - total flux cleaned = 0.0364331 Jy ! Total flux subtracted in 400 components = 0.0364331 Jy ! Clean residual min=-0.006416 max=0.005851 Jy/beam ! Clean residual mean=0.000008 rms=0.001325 Jy/beam ! Combined flux in latest and established models = 3.25278 Jy ! Adding 104 model components to the UV plane model. ! The established model now contains 231 components and 3.25278 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.369 mas, bmaj=3.345 mas, bpa=-6.518 degrees ! Estimated noise=0.415397 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00662378 Jy ! Component: 100 - total flux cleaned = 0.0111124 Jy ! Component: 150 - total flux cleaned = 0.0153832 Jy ! Component: 200 - total flux cleaned = 0.0191382 Jy ! Component: 250 - total flux cleaned = 0.0224556 Jy ! Component: 300 - total flux cleaned = 0.0255578 Jy ! Component: 350 - total flux cleaned = 0.0280317 Jy ! Component: 400 - total flux cleaned = 0.0303714 Jy ! Total flux subtracted in 400 components = 0.0303714 Jy ! Clean residual min=-0.004819 max=0.004640 Jy/beam ! Clean residual mean=0.000009 rms=0.001226 Jy/beam ! Combined flux in latest and established models = 3.28315 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.0020943 Jy ! Component: 100 - total flux cleaned = 0.00409004 Jy ! Component: 150 - total flux cleaned = 0.00623522 Jy ! Component: 200 - total flux cleaned = 0.00806393 Jy ! Component: 250 - total flux cleaned = 0.00982404 Jy ! Component: 300 - total flux cleaned = 0.0115255 Jy ! Component: 350 - total flux cleaned = 0.0130687 Jy ! Component: 400 - total flux cleaned = 0.0143669 Jy ! Total flux subtracted in 400 components = 0.0143669 Jy ! Clean residual min=-0.004057 max=0.004086 Jy/beam ! Clean residual mean=0.000007 rms=0.001058 Jy/beam ! Combined flux in latest and established models = 3.29751 Jy ! Adding 163 model components to the UV plane model. ! The established model now contains 378 components and 3.29752 Jy ! Inverting map ! Added new window around map position (15.8, -51.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000691614 Jy ! Component: 100 - total flux cleaned = 0.00230694 Jy ! Component: 150 - total flux cleaned = 0.00409442 Jy ! Component: 200 - total flux cleaned = 0.00578564 Jy ! Component: 250 - total flux cleaned = 0.00702516 Jy ! Component: 300 - total flux cleaned = 0.00858408 Jy ! Component: 350 - total flux cleaned = 0.0097366 Jy ! Component: 400 - total flux cleaned = 0.0108559 Jy ! Total flux subtracted in 400 components = 0.0108559 Jy ! Clean residual min=-0.003583 max=0.003318 Jy/beam ! Clean residual mean=0.000004 rms=0.000885 Jy/beam ! Combined flux in latest and established models = 3.30837 Jy ! Adding 65 model components to the UV plane model. ! The established model now contains 425 components and 3.30837 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=0.8352 mas, bmaj=2.151 mas, bpa=-7.488 degrees ! Estimated noise=0.378057 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000435424 Jy ! Component: 100 - total flux cleaned = -0.000249549 Jy ! Component: 150 - total flux cleaned = 8.29501e-05 Jy ! Component: 200 - total flux cleaned = 0.00103345 Jy ! Component: 250 - total flux cleaned = 0.0013369 Jy ! Component: 300 - total flux cleaned = 0.00163745 Jy ! Component: 350 - total flux cleaned = 0.0019294 Jy ! Component: 400 - total flux cleaned = 0.002901 Jy ! Total flux subtracted in 400 components = 0.002901 Jy ! Clean residual min=-0.005348 max=0.004850 Jy/beam ! Clean residual mean=0.000003 rms=0.000931 Jy/beam ! Combined flux in latest and established models = 3.31127 Jy selfcal ! Performing phase self-cal ! Adding 76 model components to the UV plane model. ! The established model now contains 486 components and 3.31127 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.092607Jy sigma=2.716495 ! Fit after self-cal, rms=0.092206Jy sigma=2.705291 wmodel /scr/pima/rdv93_uvs/J0339-0146_X_map.mod ! Writing 486 model components to file: /scr/pima/rdv93_uvs/J0339-0146_X_map.mod wobs /scr/pima/rdv93_uvs/J0339-0146_X_uvs.fits ! Applying 1877 buffered edits. ! Writing UV FITS file: /scr/pima/rdv93_uvs/J0339-0146_X_uvs.fits wwins /scr/pima/rdv93_uvs/J0339-0146_X_map.win ! wwins: Wrote 15 windows to /scr/pima/rdv93_uvs/J0339-0146_X_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.000863237 Jy ! Component: 100 - total flux cleaned = -0.00142123 Jy ! Component: 150 - total flux cleaned = -0.00214494 Jy ! Component: 200 - total flux cleaned = -0.00284448 Jy ! Component: 250 - total flux cleaned = -0.00334073 Jy ! Component: 300 - total flux cleaned = -0.0039459 Jy ! Component: 350 - total flux cleaned = -0.00442032 Jy ! Component: 400 - total flux cleaned = -0.0048287 Jy ! Component: 450 - total flux cleaned = -0.00500151 Jy ! Component: 500 - total flux cleaned = -0.00545528 Jy ! Component: 550 - total flux cleaned = -0.00579036 Jy ! Component: 600 - total flux cleaned = -0.00606587 Jy ! Component: 650 - total flux cleaned = -0.00644663 Jy ! Component: 700 - total flux cleaned = -0.00666098 Jy ! Component: 750 - total flux cleaned = -0.00682004 Jy ! Component: 800 - total flux cleaned = -0.00703034 Jy ! Component: 850 - total flux cleaned = -0.00708269 Jy ! Component: 900 - total flux cleaned = -0.00723748 Jy ! Component: 950 - total flux cleaned = -0.00733983 Jy ! Component: 1000 - total flux cleaned = -0.00754197 Jy ! Component: 1050 - total flux cleaned = -0.00764259 Jy ! Component: 1100 - total flux cleaned = -0.00779237 Jy ! Component: 1150 - total flux cleaned = -0.00808859 Jy ! Component: 1200 - total flux cleaned = -0.00823525 Jy ! Component: 1250 - total flux cleaned = -0.00833176 Jy ! Component: 1300 - total flux cleaned = -0.00842819 Jy ! Component: 1350 - total flux cleaned = -0.008619 Jy ! Component: 1400 - total flux cleaned = -0.00852475 Jy ! Component: 1450 - total flux cleaned = -0.00880712 Jy ! Component: 1500 - total flux cleaned = -0.00890074 Jy ! Component: 1550 - total flux cleaned = -0.0088538 Jy ! Component: 1600 - total flux cleaned = -0.00899245 Jy ! Component: 1650 - total flux cleaned = -0.00922158 Jy ! Component: 1700 - total flux cleaned = -0.00903886 Jy ! Component: 1750 - total flux cleaned = -0.00908435 Jy ! Component: 1800 - total flux cleaned = -0.00926426 Jy ! Component: 1850 - total flux cleaned = -0.00921975 Jy ! Component: 1900 - total flux cleaned = -0.00917562 Jy ! Component: 1950 - total flux cleaned = -0.00922009 Jy ! Component: 2000 - total flux cleaned = -0.00926407 Jy ! Component: 2050 - total flux cleaned = -0.00943853 Jy ! Component: 2100 - total flux cleaned = -0.00943791 Jy ! Component: 2150 - total flux cleaned = -0.00956714 Jy ! Component: 2200 - total flux cleaned = -0.00935249 Jy ! Component: 2250 - total flux cleaned = -0.00939521 Jy ! Component: 2300 - total flux cleaned = -0.00926771 Jy ! Component: 2350 - total flux cleaned = -0.00926771 Jy ! Component: 2400 - total flux cleaned = -0.00951965 Jy ! Component: 2450 - total flux cleaned = -0.00939388 Jy ! Component: 2500 - total flux cleaned = -0.00939391 Jy ! Component: 2550 - total flux cleaned = -0.00947647 Jy ! Component: 2600 - total flux cleaned = -0.00931176 Jy ! Component: 2650 - total flux cleaned = -0.00951692 Jy ! Component: 2700 - total flux cleaned = -0.00947627 Jy ! Component: 2750 - total flux cleaned = -0.00939501 Jy ! Component: 2800 - total flux cleaned = -0.00943531 Jy ! Component: 2850 - total flux cleaned = -0.00939536 Jy ! Component: 2900 - total flux cleaned = -0.00947527 Jy ! Component: 2950 - total flux cleaned = -0.00947483 Jy ! Component: 3000 - total flux cleaned = -0.00943517 Jy ! Component: 3050 - total flux cleaned = -0.00943516 Jy ! Component: 3100 - total flux cleaned = -0.00947447 Jy ! Component: 3150 - total flux cleaned = -0.00951381 Jy ! Component: 3200 - total flux cleaned = -0.00963039 Jy ! Component: 3250 - total flux cleaned = -0.00963035 Jy ! Component: 3300 - total flux cleaned = -0.00970706 Jy ! Component: 3350 - total flux cleaned = -0.0096687 Jy ! Component: 3400 - total flux cleaned = -0.00970693 Jy ! Component: 3450 - total flux cleaned = -0.00974495 Jy ! Component: 3500 - total flux cleaned = -0.00989647 Jy ! Component: 3550 - total flux cleaned = -0.0099341 Jy ! Component: 3600 - total flux cleaned = -0.00993414 Jy ! Component: 3650 - total flux cleaned = -0.00993398 Jy ! Component: 3700 - total flux cleaned = -0.0101576 Jy ! Component: 3750 - total flux cleaned = -0.0101205 Jy ! Component: 3800 - total flux cleaned = -0.0100467 Jy ! Component: 3850 - total flux cleaned = -0.0101203 Jy ! Component: 3900 - total flux cleaned = -0.0103397 Jy ! Component: 3950 - total flux cleaned = -0.0101573 Jy ! Component: 4000 - total flux cleaned = -0.010339 Jy ! Component: 4050 - total flux cleaned = -0.0101944 Jy ! Total flux subtracted in 4096 components = -0.0102666 Jy ! Clean residual min=-0.002126 max=0.002223 Jy/beam ! Clean residual mean=0.000002 rms=0.000620 Jy/beam ! Combined flux in latest and established models = 3.30101 Jy keep ! Adding 282 model components to the UV plane model. ! The established model now contains 759 components and 3.301 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv93_uvs/J0339-0146_X\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.8352 x 2.151 at -7.488 degrees (North through East) ! Clean map min=-0.0044334 max=2.3449 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv93_uvs/J0339-0146_X_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= 2.3448 rms= 0.00061 Jy quit ! Quitting program ! Log file /scr/pima/rdv93_uvs/J0339-0146_X_uvs.log closed on Tue Oct 8 15:43:34 2013