! Started logfile: /scr/pima/rdv94_uvs/J1231-3533_S_uvs.log on Tue Oct 8 12:19:00 2013 obs /scr/pima/rdv94_uvs/J1231-3533_S_uva.fits ! Reading UV FITS file: /scr/pima/rdv94_uvs/J1231-3533_S_uva.fits ! AN table 1: 22 integrations on 15 of 15 possible baselines. ! AN table 2: 22 integrations on 10 of 10 possible baselines. ! Apparent sampling: 0.88 visibilities/baseline/integration-bin. ! Found source: J1231-3533 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 2.23299e+09 1.6e+07 1 1.6e+07 ! 02 2 2.26299e+09 1.6e+07 1 1.6e+07 ! 03 3 2.35299e+09 1.6e+07 1 1.6e+07 ! 04 4 2.37299e+09 1.6e+07 1 1.6e+07 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 1936 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.5 print "taper_size =",taper_size ! taper_size = 20 ![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv94_uvs/J1231-3533_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:FD self-cal status=correctable weight=1 ! 1:KP self-cal status=correctable weight=1 ! 1:LA 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 ! 2:BR self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:KP self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:OV self-cal status=correctable weight=1 system "rm /tmp/ff.ff" mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.500x0.500 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. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.998472Jy sigma=0.722135 ! Fit after self-cal, rms=0.958733Jy sigma=0.693783 ! 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 \] ![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=10.4 mas, bmaj=31.53 mas, bpa=-5.577 degrees ! Estimated noise=42.7862 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=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0341102 Jy ! Component: 100 - total flux cleaned = 0.040868 Jy ! Total flux subtracted in 100 components = 0.040868 Jy ! Clean residual min=-0.004295 max=0.004291 Jy/beam ! Clean residual mean=0.000103 rms=0.001512 Jy/beam ! Combined flux in latest and established models = 0.040868 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 16 components and 0.040868 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.047750Jy sigma=0.033937 ! Fit after self-cal, rms=0.046514Jy sigma=0.033090 ! 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=11.51 mas, bmaj=34.93 mas, bpa=-5.428 degrees ! Estimated noise=22.3957 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.00491878 Jy ! Component: 100 - total flux cleaned = 0.00686659 Jy ! Component: 150 - total flux cleaned = 0.00739273 Jy ! Component: 200 - total flux cleaned = 0.00821667 Jy ! Total flux subtracted in 200 components = 0.00821667 Jy ! Clean residual min=-0.001823 max=0.001819 Jy/beam ! Clean residual mean=0.000060 rms=0.000739 Jy/beam ! Combined flux in latest and established models = 0.0490847 Jy ! Performing phase self-cal ! Adding 62 model components to the UV plane model. ! The established model now contains 77 components and 0.0490847 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.045542Jy sigma=0.032419 ! Fit after self-cal, rms=0.045460Jy sigma=0.032379 ! 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: ! FD 1.18 KP 1.07 LA 1.38 NL 0.75 ! OV 1.01 PT 0.81 ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.98 FD 1.01 KP 1.14 LA 1.24 ! OV 0.92 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! FD 1.16 KP 1.26 LA 1.23 NL 0.81 ! OV 0.86 PT 0.85 ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.84 FD 1.13 KP 1.24 LA 1.06 ! OV 1.07 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! FD 0.99 KP 1.02 LA 0.89 NL 1.00 ! OV 1.01 PT 1.04 ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.91 FD 1.09 KP 0.91 LA 1.02 ! OV 0.89 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! FD 1.18 KP 0.82 LA 0.67 NL 0.80 ! OV 0.91 PT 1.03 ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.84 FD 1.12 KP 0.93 LA 0.72 ! OV 0.96 ! ! ! Fit before self-cal, rms=0.045460Jy sigma=0.032379 ! Fit after self-cal, rms=0.042916Jy sigma=0.031506 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=11.32 mas, bmaj=34.24 mas, bpa=-6.134 degrees ! Estimated noise=21.6298 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.00094917 Jy ! Component: 100 - total flux cleaned = 0.00141047 Jy ! Component: 150 - total flux cleaned = 0.00166114 Jy ! Component: 200 - total flux cleaned = 0.00198347 Jy ! Total flux subtracted in 200 components = 0.00198347 Jy ! Clean residual min=-0.000558 max=0.000560 Jy/beam ! Clean residual mean=-0.000003 rms=0.000231 Jy/beam ! Combined flux in latest and established models = 0.0510681 Jy ! Performing phase self-cal ! Adding 80 model components to the UV plane model. ! The established model now contains 157 components and 0.0510681 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.042805Jy sigma=0.031443 ! Fit after self-cal, rms=0.043112Jy sigma=0.031340 ! 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.043112Jy sigma=0.031340 ! Fit after self-cal, rms=0.043917Jy sigma=0.031317 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=11.28 mas, bmaj=34.14 mas, bpa=-6.341 degrees ! Estimated noise=21.6633 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000364343 Jy ! Component: 100 - total flux cleaned = 0.000538269 Jy ! Component: 150 - total flux cleaned = 0.000759994 Jy ! Component: 200 - total flux cleaned = 0.00095562 Jy ! Total flux subtracted in 200 components = 0.00095562 Jy ! Clean residual min=-0.000451 max=0.000453 Jy/beam ! Clean residual mean=-0.000004 rms=0.000190 Jy/beam ! Combined flux in latest and established models = 0.0520237 Jy selfcal ! Performing phase self-cal ! Adding 66 model components to the UV plane model. ! The established model now contains 217 components and 0.0520237 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.043849Jy sigma=0.031289 ! Fit after self-cal, rms=0.043960Jy sigma=0.031280 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 = -5.98607e-05 Jy ! Component: 100 - total flux cleaned = 0.000104713 Jy ! Component: 150 - total flux cleaned = 0.000308455 Jy ! Component: 200 - total flux cleaned = 0.000544588 Jy ! Total flux subtracted in 200 components = 0.000544588 Jy ! Clean residual min=-0.000381 max=0.000380 Jy/beam ! Clean residual mean=-0.000004 rms=0.000161 Jy/beam ! Combined flux in latest and established models = 0.0525683 Jy ! Performing phase self-cal ! Adding 61 model components to the UV plane model. ! The established model now contains 261 components and 0.0525683 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.043921Jy sigma=0.031265 ! Fit after self-cal, rms=0.043931Jy sigma=0.031264 ! 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.043931Jy sigma=0.031264 ! Fit after self-cal, rms=0.044748Jy sigma=0.031256 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=11.26 mas, bmaj=34.09 mas, bpa=-6.347 degrees ! Estimated noise=21.7992 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000290357 Jy ! Component: 100 - total flux cleaned = -5.92171e-05 Jy ! Component: 150 - total flux cleaned = 0.00021605 Jy ! Component: 200 - total flux cleaned = 0.000493323 Jy ! Total flux subtracted in 200 components = 0.000493323 Jy ! Clean residual min=-0.000343 max=0.000344 Jy/beam ! Clean residual mean=-0.000001 rms=0.000142 Jy/beam ! Combined flux in latest and established models = 0.0530616 Jy selfcal ! Performing phase self-cal ! Adding 66 model components to the UV plane model. ! The established model now contains 313 components and 0.0530616 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.044713Jy sigma=0.031242 ! Fit after self-cal, rms=0.044766Jy sigma=0.031240 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 = 1.87498e-05 Jy ! Component: 100 - total flux cleaned = 0.000226954 Jy ! Component: 150 - total flux cleaned = 0.000421826 Jy ! Component: 200 - total flux cleaned = 0.000625126 Jy ! Total flux subtracted in 200 components = 0.000625126 Jy ! Clean residual min=-0.000300 max=0.000302 Jy/beam ! Clean residual mean=-0.000002 rms=0.000125 Jy/beam ! Combined flux in latest and established models = 0.0536867 Jy ! Performing phase self-cal ! Adding 66 model components to the UV plane model. ! The established model now contains 357 components and 0.0536867 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.044744Jy sigma=0.031231 ! Fit after self-cal, rms=0.044749Jy sigma=0.031230 ! 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.044749Jy sigma=0.031230 ! Fit after self-cal, rms=0.045108Jy sigma=0.031224 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=12.36 mas, bmaj=36.63 mas, bpa=-4.453 degrees ! Estimated noise=22.0883 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 7.42093e-05 Jy ! Component: 100 - total flux cleaned = 0.000296255 Jy ! Component: 150 - total flux cleaned = 0.000461954 Jy ! Component: 200 - total flux cleaned = 0.000636196 Jy ! Total flux subtracted in 200 components = 0.000636196 Jy ! Clean residual min=-0.000282 max=0.000283 Jy/beam ! Clean residual mean=0.000000 rms=0.000120 Jy/beam ! Combined flux in latest and established models = 0.0543229 Jy selfcal ! Performing phase self-cal ! Adding 87 model components to the UV plane model. ! The established model now contains 419 components and 0.0543229 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.045086Jy sigma=0.031216 ! Fit after self-cal, rms=0.045107Jy sigma=0.031214 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 = 8.59776e-05 Jy ! Component: 100 - total flux cleaned = 0.000218401 Jy ! Component: 150 - total flux cleaned = 0.000345657 Jy ! Component: 200 - total flux cleaned = 0.000483318 Jy ! Total flux subtracted in 200 components = 0.000483318 Jy ! Clean residual min=-0.000250 max=0.000250 Jy/beam ! Clean residual mean=-0.000001 rms=0.000108 Jy/beam ! Combined flux in latest and established models = 0.0548063 Jy ! Performing phase self-cal ! Adding 66 model components to the UV plane model. ! The established model now contains 458 components and 0.0548063 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.045092Jy sigma=0.031208 ! Fit after self-cal, rms=0.045095Jy sigma=0.031208 ! 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:FD,true ! 1:FD 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:NL,true ! 1:NL self-cal status=fixed weight=1 selfant 1:OV,true ! 1:OV self-cal status=fixed weight=1 selfant 1:PT,true ! 1:PT self-cal status=fixed weight=1 selfant 2:BR,true ! 2:BR self-cal status=fixed weight=1 selfant 2:FD,true ! 2:FD self-cal status=fixed weight=1 selfant 2:KP,true ! 2:KP self-cal status=fixed weight=1 selfant 2:LA,true ! 2:LA self-cal status=fixed weight=1 selfant 2:OV,true ! 2:OV self-cal status=fixed weight=1 ![Exited script file: /tmp/ff.ff] selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.045095Jy sigma=0.031208 ! Fit after self-cal, rms=0.045095Jy sigma=0.031208 selfant "",false ! 1:FD self-cal status=correctable weight=1 ! 1:KP self-cal status=correctable weight=1 ! 1:LA 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 ! 2:BR self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:KP self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:OV self-cal status=correctable weight=1 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.045095Jy sigma=0.031208 ! Fit after self-cal, rms=0.044744Jy sigma=0.031204 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=10.11 mas, bmaj=30.82 mas, bpa=-7.23 degrees ! Estimated noise=40.8382 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.000210425 Jy ! Component: 100 - total flux cleaned = 0.000193554 Jy ! Component: 150 - total flux cleaned = 0.000846169 Jy ! Component: 200 - total flux cleaned = 0.00125937 Jy ! Total flux subtracted in 200 components = 0.00125937 Jy ! Clean residual min=-0.000799 max=0.000801 Jy/beam ! Clean residual mean=-0.000004 rms=0.000354 Jy/beam ! Combined flux in latest and established models = 0.0560656 Jy ! Performing phase self-cal ! Adding 85 model components to the UV plane model. ! The established model now contains 543 components and 0.0560656 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.044775Jy sigma=0.031247 ! Fit after self-cal, rms=0.044776Jy sigma=0.031229 ! 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=11.26 mas, bmaj=34.1 mas, bpa=-6.449 degrees ! Estimated noise=21.7384 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.000517885 Jy ! Component: 100 - total flux cleaned = 0.000487819 Jy ! Component: 150 - total flux cleaned = 0.000687404 Jy ! Component: 200 - total flux cleaned = 0.000953622 Jy ! Component: 250 - total flux cleaned = 0.00123853 Jy ! Component: 300 - total flux cleaned = 0.00160322 Jy ! Component: 350 - total flux cleaned = 0.00196561 Jy ! Component: 400 - total flux cleaned = 0.00225839 Jy ! Total flux subtracted in 400 components = 0.00225839 Jy ! Clean residual min=-0.000280 max=0.000282 Jy/beam ! Clean residual mean=0.000003 rms=0.000113 Jy/beam ! Combined flux in latest and established models = 0.058324 Jy ! Performing phase self-cal ! Adding 118 model components to the UV plane model. ! The established model now contains 655 components and 0.058324 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.044728Jy sigma=0.031192 ! Fit after self-cal, rms=0.044715Jy sigma=0.031188 ! Inverting map ! Added new window around map position (-92.5, -73.5). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -8.84056e-05 Jy ! Component: 100 - total flux cleaned = -0.000127612 Jy ! Component: 150 - total flux cleaned = -9.46718e-05 Jy ! Component: 200 - total flux cleaned = -8.46184e-05 Jy ! Component: 250 - total flux cleaned = -0.000101992 Jy ! Component: 300 - total flux cleaned = -0.000133597 Jy ! Component: 350 - total flux cleaned = -0.000178124 Jy ! Component: 400 - total flux cleaned = -0.00021285 Jy ! Total flux subtracted in 400 components = -0.00021285 Jy ! Clean residual min=-0.000363 max=0.000374 Jy/beam ! Clean residual mean=0.000003 rms=0.000105 Jy/beam ! Combined flux in latest and established models = 0.0581111 Jy ! Performing phase self-cal ! Adding 94 model components to the UV plane model. ! The established model now contains 749 components and 0.0581111 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.044702Jy sigma=0.031182 ! Fit after self-cal, rms=0.044702Jy sigma=0.031181 ! 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 = 4.85756e-05 Jy ! Component: 100 - total flux cleaned = 4.62313e-05 Jy ! Component: 150 - total flux cleaned = 4.61104e-05 Jy ! Component: 200 - total flux cleaned = 3.21628e-05 Jy ! Component: 250 - total flux cleaned = 2.56929e-05 Jy ! Component: 300 - total flux cleaned = -4.85447e-06 Jy ! Component: 350 - total flux cleaned = -1.07206e-05 Jy ! Component: 400 - total flux cleaned = -3.91877e-05 Jy ! Total flux subtracted in 400 components = -3.91877e-05 Jy ! Clean residual min=-0.000393 max=0.000378 Jy/beam ! Clean residual mean=0.000003 rms=0.000103 Jy/beam ! Combined flux in latest and established models = 0.0580719 Jy ! Performing phase self-cal ! Adding 78 model components to the UV plane model. ! The established model now contains 810 components and 0.0580719 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.044695Jy sigma=0.031179 ! Fit after self-cal, rms=0.044696Jy sigma=0.031178 ! 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. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.044696Jy sigma=0.031178 ! Fit after self-cal, rms=0.044457Jy sigma=0.031176 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=11.25 mas, bmaj=34.09 mas, bpa=-6.558 degrees ! Estimated noise=21.646 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000159738 Jy ! Component: 100 - total flux cleaned = 0.000228287 Jy ! Component: 150 - total flux cleaned = 0.000266655 Jy ! Component: 200 - total flux cleaned = 0.00028064 Jy ! Component: 250 - total flux cleaned = 0.000280683 Jy ! Component: 300 - total flux cleaned = 0.000274505 Jy ! Component: 350 - total flux cleaned = 0.000257374 Jy ! Component: 400 - total flux cleaned = 0.000246565 Jy ! Total flux subtracted in 400 components = 0.000246565 Jy ! Clean residual min=-0.000424 max=0.000402 Jy/beam ! Clean residual mean=0.000004 rms=0.000120 Jy/beam ! Combined flux in latest and established models = 0.0583185 Jy selfcal ! Performing phase self-cal ! Adding 104 model components to the UV plane model. ! The established model now contains 904 components and 0.0583185 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.044450Jy sigma=0.031174 ! Fit after self-cal, rms=0.044448Jy sigma=0.031173 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.000151633 Jy ! Component: 100 - total flux cleaned = 0.000164627 Jy ! Component: 150 - total flux cleaned = 0.000158552 Jy ! Component: 200 - total flux cleaned = 0.000152654 Jy ! Component: 250 - total flux cleaned = 0.000142199 Jy ! Component: 300 - total flux cleaned = 0.000142415 Jy ! Component: 350 - total flux cleaned = 0.000127739 Jy ! Component: 400 - total flux cleaned = 9.97316e-05 Jy ! Total flux subtracted in 400 components = 9.97316e-05 Jy ! Clean residual min=-0.000422 max=0.000391 Jy/beam ! Clean residual mean=0.000004 rms=0.000115 Jy/beam ! Combined flux in latest and established models = 0.0584182 Jy ! Performing phase self-cal ! Adding 80 model components to the UV plane model. ! The established model now contains 956 components and 0.0584182 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.044444Jy sigma=0.031172 ! Fit after self-cal, rms=0.044445Jy sigma=0.031172 ! 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.044445Jy sigma=0.031172 ! Fit after self-cal, rms=0.044847Jy sigma=0.031169 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=11.24 mas, bmaj=34.05 mas, bpa=-6.53 degrees ! Estimated noise=21.7481 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000136859 Jy ! Component: 100 - total flux cleaned = 0.000160089 Jy ! Component: 150 - total flux cleaned = 0.000187157 Jy ! Component: 200 - total flux cleaned = 0.000205622 Jy ! Component: 250 - total flux cleaned = 0.00022847 Jy ! Component: 300 - total flux cleaned = 0.000239169 Jy ! Component: 350 - total flux cleaned = 0.000244321 Jy ! Component: 400 - total flux cleaned = 0.000244382 Jy ! Total flux subtracted in 400 components = 0.000244382 Jy ! Clean residual min=-0.000512 max=0.000426 Jy/beam ! Clean residual mean=0.000003 rms=0.000121 Jy/beam ! Combined flux in latest and established models = 0.0586626 Jy selfcal ! Performing phase self-cal ! Adding 87 model components to the UV plane model. ! The established model now contains 1032 components and 0.0586626 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.044841Jy sigma=0.031168 ! Fit after self-cal, rms=0.044854Jy sigma=0.031167 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.000113976 Jy ! Component: 100 - total flux cleaned = 0.000113721 Jy ! Component: 150 - total flux cleaned = 0.000134243 Jy ! Component: 200 - total flux cleaned = 0.000143994 Jy ! Component: 250 - total flux cleaned = 0.000162537 Jy ! Component: 300 - total flux cleaned = 0.000166945 Jy ! Component: 350 - total flux cleaned = 0.000167016 Jy ! Component: 400 - total flux cleaned = 0.000175308 Jy ! Total flux subtracted in 400 components = 0.000175308 Jy ! Clean residual min=-0.000496 max=0.000410 Jy/beam ! Clean residual mean=0.000004 rms=0.000117 Jy/beam ! Combined flux in latest and established models = 0.0588379 Jy ! Performing phase self-cal ! Adding 67 model components to the UV plane model. ! The established model now contains 1075 components and 0.0588379 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.044850Jy sigma=0.031166 ! Fit after self-cal, rms=0.044853Jy sigma=0.031166 ! 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.044853Jy sigma=0.031166 ! Fit after self-cal, rms=0.049986Jy sigma=0.030154 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=12.34 mas, bmaj=36.74 mas, bpa=-4.861 degrees ! Estimated noise=23.2419 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00126624 Jy ! Component: 100 - total flux cleaned = 0.00208632 Jy ! Component: 150 - total flux cleaned = 0.00260507 Jy ! Component: 200 - total flux cleaned = 0.00288885 Jy ! Component: 250 - total flux cleaned = 0.00299782 Jy ! Component: 300 - total flux cleaned = 0.0030778 Jy ! Component: 350 - total flux cleaned = 0.00316241 Jy ! Component: 400 - total flux cleaned = 0.00317488 Jy ! Total flux subtracted in 400 components = 0.00317488 Jy ! Clean residual min=-0.000630 max=0.001024 Jy/beam ! Clean residual mean=0.000003 rms=0.000234 Jy/beam ! Combined flux in latest and established models = 0.0620128 Jy selfcal ! Performing phase self-cal ! Adding 125 model components to the UV plane model. ! The established model now contains 1192 components and 0.0620128 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.049941Jy sigma=0.030109 ! Fit after self-cal, rms=0.050319Jy sigma=0.029984 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 (101.5, -128). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000601705 Jy ! Component: 100 - total flux cleaned = 0.000831068 Jy ! Component: 150 - total flux cleaned = 0.000830353 Jy ! Component: 200 - total flux cleaned = 0.000830014 Jy ! Component: 250 - total flux cleaned = 0.000860882 Jy ! Component: 300 - total flux cleaned = 0.000917348 Jy ! Component: 350 - total flux cleaned = 0.00101981 Jy ! Component: 400 - total flux cleaned = 0.00105585 Jy ! Total flux subtracted in 400 components = 0.00105585 Jy ! Clean residual min=-0.000533 max=0.000666 Jy/beam ! Clean residual mean=0.000002 rms=0.000151 Jy/beam ! Combined flux in latest and established models = 0.0630687 Jy ! Performing phase self-cal ! Adding 106 model components to the UV plane model. ! The established model now contains 1295 components and 0.0630687 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.050286Jy sigma=0.029962 ! Fit after self-cal, rms=0.050285Jy sigma=0.029959 ! Inverting map ! Added new window around map position (127.5, 84). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000258248 Jy ! Component: 100 - total flux cleaned = 0.000347316 Jy ! Component: 150 - total flux cleaned = 0.000422373 Jy ! Component: 200 - total flux cleaned = 0.000504774 Jy ! Component: 250 - total flux cleaned = 0.000580865 Jy ! Component: 300 - total flux cleaned = 0.000581302 Jy ! Component: 350 - total flux cleaned = 0.000536386 Jy ! Component: 400 - total flux cleaned = 0.000515089 Jy ! Total flux subtracted in 400 components = 0.000515089 Jy ! Clean residual min=-0.000437 max=0.000386 Jy/beam ! Clean residual mean=-0.000003 rms=0.000115 Jy/beam ! Combined flux in latest and established models = 0.0635838 Jy ! Performing phase self-cal ! Adding 114 model components to the UV plane model. ! The established model now contains 1397 components and 0.0635838 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.050265Jy sigma=0.029948 ! Fit after self-cal, rms=0.050260Jy sigma=0.029946 ! 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=11.23 mas, bmaj=34.21 mas, bpa=-6.787 degrees ! Estimated noise=22.9919 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 11.23 x 34.21 at -6.787 degrees (North through East) ! Clean map min=-0.0024391 max=0.040441 Jy/beam if(peak(flux) > thresh) selfcal true, true clean clean_niter,clean_gain selfcal flux_cutoff = imstat(rms) * dynam repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) selfcal true, true end if print "=== CLearing model and starting over ===" ! === CLearing model and starting over === clrmod true ! clrmod: Cleared the established and tentative models. uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=10.1 mas, bmaj=31.17 mas, bpa=-7.672 degrees ! Estimated noise=43.4627 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0246746 Jy ! Component: 100 - total flux cleaned = 0.0332558 Jy ! Component: 150 - total flux cleaned = 0.0399384 Jy ! Component: 200 - total flux cleaned = 0.0449951 Jy ! Component: 250 - total flux cleaned = 0.0477063 Jy ! Component: 300 - total flux cleaned = 0.0465452 Jy ! Component: 350 - total flux cleaned = 0.0473939 Jy ! Component: 400 - total flux cleaned = 0.0473979 Jy ! Total flux subtracted in 400 components = 0.0473979 Jy ! Clean residual min=-0.008094 max=0.015184 Jy/beam ! Clean residual mean=0.000061 rms=0.002930 Jy/beam ! Combined flux in latest and established models = 0.0473979 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 95 model components to the UV plane model. ! The established model now contains 95 components and 0.0473979 Jy ! Inverting map and beam ! Estimated beam: bmin=11.23 mas, bmaj=34.21 mas, bpa=-6.787 degrees ! Estimated noise=22.9919 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) ! Added new window around map position (-0.5, 2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0104143 Jy ! Component: 100 - total flux cleaned = 0.00991922 Jy ! Component: 150 - total flux cleaned = 0.0101303 Jy ! Component: 200 - total flux cleaned = 0.00963987 Jy ! Component: 250 - total flux cleaned = 0.00878371 Jy ! Component: 300 - total flux cleaned = 0.00751502 Jy ! Component: 350 - total flux cleaned = 0.00635312 Jy ! Component: 400 - total flux cleaned = 0.00527997 Jy ! Total flux subtracted in 400 components = 0.00527997 Jy ! Clean residual min=-0.002752 max=0.003347 Jy/beam ! Clean residual mean=-0.000051 rms=0.000949 Jy/beam ! Combined flux in latest and established models = 0.0526778 Jy ! Adding 133 model components to the UV plane model. ! The established model now contains 227 components and 0.0526778 Jy ! Inverting map ! Added new window around map position (-116.5, 57.5). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00289619 Jy ! Component: 100 - total flux cleaned = 0.00375757 Jy ! Component: 150 - total flux cleaned = 0.00410082 Jy ! Component: 200 - total flux cleaned = 0.00426357 Jy ! Component: 250 - total flux cleaned = 0.00388333 Jy ! Component: 300 - total flux cleaned = 0.00366415 Jy ! Component: 350 - total flux cleaned = 0.00337867 Jy ! Component: 400 - total flux cleaned = 0.00310399 Jy ! Total flux subtracted in 400 components = 0.00310399 Jy ! Clean residual min=-0.001901 max=0.001828 Jy/beam ! Clean residual mean=-0.000042 rms=0.000598 Jy/beam ! Combined flux in latest and established models = 0.0557818 Jy ! Adding 143 model components to the UV plane model. ! The established model now contains 353 components and 0.0557818 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=12.34 mas, bmaj=36.74 mas, bpa=-4.861 degrees ! Estimated noise=23.2419 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000571927 Jy ! Component: 100 - total flux cleaned = -0.00103989 Jy ! Component: 150 - total flux cleaned = -0.00131139 Jy ! Component: 200 - total flux cleaned = -0.00141378 Jy ! Component: 250 - total flux cleaned = -0.00180616 Jy ! Component: 300 - total flux cleaned = -0.00199598 Jy ! Component: 350 - total flux cleaned = -0.00231958 Jy ! Component: 400 - total flux cleaned = -0.00259139 Jy ! Total flux subtracted in 400 components = -0.00259139 Jy ! Clean residual min=-0.001427 max=0.001269 Jy/beam ! Clean residual mean=-0.000031 rms=0.000422 Jy/beam ! Combined flux in latest and established models = 0.0531904 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.000444126 Jy ! Component: 100 - total flux cleaned = -0.000791041 Jy ! Component: 150 - total flux cleaned = -0.00113296 Jy ! Component: 200 - total flux cleaned = -0.00155283 Jy ! Component: 250 - total flux cleaned = -0.00184138 Jy ! Component: 300 - total flux cleaned = -0.00204423 Jy ! Component: 350 - total flux cleaned = -0.0022856 Jy ! Component: 400 - total flux cleaned = -0.00240615 Jy ! Total flux subtracted in 400 components = -0.00240615 Jy ! Clean residual min=-0.001332 max=0.001007 Jy/beam ! Clean residual mean=-0.000031 rms=0.000366 Jy/beam ! Combined flux in latest and established models = 0.0507843 Jy ! Adding 205 model components to the UV plane model. ! The established model now contains 532 components and 0.0507843 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=11.23 mas, bmaj=34.21 mas, bpa=-6.787 degrees ! Estimated noise=22.9919 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000627013 Jy ! Component: 100 - total flux cleaned = -0.00059285 Jy ! Component: 150 - total flux cleaned = -0.000526763 Jy ! Component: 200 - total flux cleaned = -0.000401988 Jy ! Component: 250 - total flux cleaned = -0.000373174 Jy ! Component: 300 - total flux cleaned = -0.000286499 Jy ! Component: 350 - total flux cleaned = -0.000230554 Jy ! Component: 400 - total flux cleaned = -0.000148176 Jy ! Total flux subtracted in 400 components = -0.000148176 Jy ! Clean residual min=-0.001050 max=0.000876 Jy/beam ! Clean residual mean=-0.000020 rms=0.000262 Jy/beam ! Combined flux in latest and established models = 0.0506361 Jy selfcal ! Performing phase self-cal ! Adding 122 model components to the UV plane model. ! The established model now contains 645 components and 0.0506361 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.050486Jy sigma=0.030168 ! Fit after self-cal, rms=0.050384Jy sigma=0.030055 wmodel /scr/pima/rdv94_uvs/J1231-3533_S_map.mod ! Writing 645 model components to file: /scr/pima/rdv94_uvs/J1231-3533_S_map.mod wobs /scr/pima/rdv94_uvs/J1231-3533_S_uvs.fits ! Applying 1100 buffered edits. ! Writing UV FITS file: /scr/pima/rdv94_uvs/J1231-3533_S_uvs.fits wwins /scr/pima/rdv94_uvs/J1231-3533_S_map.win ! wwins: Wrote 5 windows to /scr/pima/rdv94_uvs/J1231-3533_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.000242242 Jy ! Component: 100 - total flux cleaned = -0.000479375 Jy ! Component: 150 - total flux cleaned = -0.000647572 Jy ! Component: 200 - total flux cleaned = -0.000785711 Jy ! Component: 250 - total flux cleaned = -0.000943902 Jy ! Component: 300 - total flux cleaned = -0.00109611 Jy ! Component: 350 - total flux cleaned = -0.00121096 Jy ! Component: 400 - total flux cleaned = -0.0013067 Jy ! Component: 450 - total flux cleaned = -0.00140789 Jy ! Component: 500 - total flux cleaned = -0.00150654 Jy ! Component: 550 - total flux cleaned = -0.00161772 Jy ! Component: 600 - total flux cleaned = -0.00169763 Jy ! Component: 650 - total flux cleaned = -0.00176197 Jy ! Component: 700 - total flux cleaned = -0.00184617 Jy ! Component: 750 - total flux cleaned = -0.00192208 Jy ! Component: 800 - total flux cleaned = -0.00199692 Jy ! Component: 850 - total flux cleaned = -0.00206395 Jy ! Component: 900 - total flux cleaned = -0.00213007 Jy ! Component: 950 - total flux cleaned = -0.00219535 Jy ! Component: 1000 - total flux cleaned = -0.00226631 Jy ! Component: 1050 - total flux cleaned = -0.0023365 Jy ! Component: 1100 - total flux cleaned = -0.00240595 Jy ! Component: 1150 - total flux cleaned = -0.00247466 Jy ! Component: 1200 - total flux cleaned = -0.00254888 Jy ! Component: 1250 - total flux cleaned = -0.00261008 Jy ! Component: 1300 - total flux cleaned = -0.00268279 Jy ! Component: 1350 - total flux cleaned = -0.00276086 Jy ! Component: 1400 - total flux cleaned = -0.0028383 Jy ! Component: 1450 - total flux cleaned = -0.00292113 Jy ! Component: 1500 - total flux cleaned = -0.00300925 Jy ! Component: 1550 - total flux cleaned = -0.00309088 Jy ! Component: 1600 - total flux cleaned = -0.00316618 Jy ! Component: 1650 - total flux cleaned = -0.00324675 Jy ! Component: 1700 - total flux cleaned = -0.00330967 Jy ! Component: 1750 - total flux cleaned = -0.00336656 Jy ! Component: 1800 - total flux cleaned = -0.00341186 Jy ! Component: 1850 - total flux cleaned = -0.00346264 Jy ! Component: 1900 - total flux cleaned = -0.00350759 Jy ! Component: 1950 - total flux cleaned = -0.00355237 Jy ! Component: 2000 - total flux cleaned = -0.00358585 Jy ! Component: 2050 - total flux cleaned = -0.00364699 Jy ! Component: 2100 - total flux cleaned = -0.00366915 Jy ! Component: 2150 - total flux cleaned = -0.00370229 Jy ! Component: 2200 - total flux cleaned = -0.00373531 Jy ! Component: 2250 - total flux cleaned = -0.00375724 Jy ! Component: 2300 - total flux cleaned = -0.0037901 Jy ! Component: 2350 - total flux cleaned = -0.00381739 Jy ! Component: 2400 - total flux cleaned = -0.00383917 Jy ! Component: 2450 - total flux cleaned = -0.00385557 Jy ! Component: 2500 - total flux cleaned = -0.00388819 Jy ! Component: 2550 - total flux cleaned = -0.00392069 Jy ! Component: 2600 - total flux cleaned = -0.00395314 Jy ! Component: 2650 - total flux cleaned = -0.00396395 Jy ! Component: 2700 - total flux cleaned = -0.00398554 Jy ! Component: 2750 - total flux cleaned = -0.00400714 Jy ! Component: 2800 - total flux cleaned = -0.00403945 Jy ! Component: 2850 - total flux cleaned = -0.00404482 Jy ! Component: 2900 - total flux cleaned = -0.00407709 Jy ! Component: 2950 - total flux cleaned = -0.00409861 Jy ! Component: 3000 - total flux cleaned = -0.00408785 Jy ! Component: 3050 - total flux cleaned = -0.00409861 Jy ! Component: 3100 - total flux cleaned = -0.004104 Jy ! Component: 3150 - total flux cleaned = -0.00412021 Jy ! Component: 3200 - total flux cleaned = -0.00412563 Jy ! Component: 3250 - total flux cleaned = -0.00413648 Jy ! Component: 3300 - total flux cleaned = -0.00413648 Jy ! Component: 3350 - total flux cleaned = -0.00414197 Jy ! Component: 3400 - total flux cleaned = -0.00414198 Jy ! Component: 3450 - total flux cleaned = -0.00415289 Jy ! Component: 3500 - total flux cleaned = -0.00414192 Jy ! Component: 3550 - total flux cleaned = -0.00413091 Jy ! Component: 3600 - total flux cleaned = -0.00414196 Jy ! Component: 3650 - total flux cleaned = -0.0041419 Jy ! Component: 3700 - total flux cleaned = -0.00413631 Jy ! Component: 3750 - total flux cleaned = -0.00413072 Jy ! Component: 3800 - total flux cleaned = -0.00414196 Jy ! Component: 3850 - total flux cleaned = -0.00413631 Jy ! Component: 3900 - total flux cleaned = -0.00413059 Jy ! Component: 3950 - total flux cleaned = -0.00413056 Jy ! Component: 4000 - total flux cleaned = -0.00413625 Jy ! Component: 4050 - total flux cleaned = -0.0041305 Jy ! Total flux subtracted in 4096 components = -0.00413052 Jy ! Clean residual min=-0.000309 max=0.000304 Jy/beam ! Clean residual mean=-0.000010 rms=0.000135 Jy/beam ! Combined flux in latest and established models = 0.0465056 Jy keep ! Adding 402 model components to the UV plane model. ! The established model now contains 1046 components and 0.0465056 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv94_uvs/J1231-3533_S\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 11.23 x 34.21 at -6.787 degrees (North through East) ! Clean map min=-0.0081177 max=0.021881 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv94_uvs/J1231-3533_S_map.fits float peak_flux peak_flux = peak(flux,max) invert ! Inverting map float image_peak image_peak = int(10000*peak_flux)/10000 float image_rms image_rms = int(100000*imstat(rms))/100000 print "=========================================" ! ========================================= print "The image peak= ", image_peak, " rms= ", image_rms, " Jy" ! The image peak= 0.0218 rms= 0.00013 Jy quit ! Quitting program ! Log file /scr/pima/rdv94_uvs/J1231-3533_S_uvs.log closed on Tue Oct 8 12:19:09 2013