! Started logfile: /scr/pima/rdv73_uvs/J0730-1141_X_uvs.log on Sat Jan 9 17:32:03 2016 obs /scr/pima/rdv73_uvs/J0730-1141_X_uva.fits ! Reading UV FITS file: /scr/pima/rdv73_uvs/J0730-1141_X_uva.fits ! AN table 1: 2 integrations on 3 of 3 possible baselines. ! AN table 2: 3 integrations on 21 of 21 possible baselines. ! AN table 3: 2 integrations on 15 of 15 possible baselines. ! AN table 4: 3 integrations on 36 of 36 possible baselines. ! AN table 5: 3 integrations on 55 of 55 possible baselines. ! AN table 6: 3 integrations on 45 of 45 possible baselines. ! Apparent sampling: 0.895464 visibilities/baseline/integration-bin. ! Found source: J0730-1141 ! ! 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 1816 visibilities. select RR ! Selecting polarization: RR, channels: 1..4 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 print "field_size =",field_size ! field_size = 1024 print "field_cell =",field_cell ! field_cell = 0.2 print "taper_size =",taper_size ! taper_size = 70 ![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv73_uvs/J0730-1141_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:SC self-cal status=correctable weight=1 ! 1:WZ self-cal status=correctable weight=1 ! 1:ZC self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:HN self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 3:FD self-cal status=correctable weight=1 ! 3:HN self-cal status=correctable weight=1 ! 3:KP self-cal status=correctable weight=1 ! 3:LA self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 3:PT self-cal status=correctable weight=1 ! 4:FD self-cal status=correctable weight=1 ! 4:HN self-cal status=correctable weight=1 ! 4:KP self-cal status=correctable weight=1 ! 4:LA self-cal status=correctable weight=1 ! 4:NL self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 4:PT self-cal status=correctable weight=1 ! 4:SC self-cal status=correctable weight=1 ! 4:WZ self-cal status=correctable weight=1 ! 5:BR self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:HN self-cal status=correctable weight=1 ! 5:KK self-cal status=correctable weight=1 ! 5:KP self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:MK self-cal status=correctable weight=1 ! 5:NL self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 5:SC self-cal status=correctable weight=1 ! 6:BR self-cal status=correctable weight=1 ! 6:FD self-cal status=correctable weight=1 ! 6:HN self-cal status=correctable weight=1 ! 6:KK self-cal status=correctable weight=1 ! 6:KP self-cal status=correctable weight=1 ! 6:LA self-cal status=correctable weight=1 ! 6:MK self-cal status=correctable weight=1 ! 6: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 6 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 5. ! A total of 1 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 6 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 5. ! A total of 1 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 6 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 5. ! A total of 1 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 6 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 5. ! A total of 1 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=4.805268Jy sigma=10.570790 ! Fit after self-cal, rms=3.751685Jy sigma=8.069608 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. print "=== Flag out IVS stations" ! === Flag out IVS stations integer lun string cfg_file lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_flag_process.csh /tmp/ff.ff" ![@/tmp/ff.ff \] flag 1:WZ ! Flagging baselines of antenna 1:WZ in the currently selected channels. flag 1:ZC ! Flagging baselines of antenna 1:ZC in the currently selected channels. flag 2:WZ ! Flagging baselines of antenna 2:WZ in the currently selected channels. flag 4:WZ ! Flagging baselines of antenna 4:WZ in the currently selected channels. flag 5:KK ! Flagging baselines of antenna 5:KK in the currently selected channels. flag 6:KK ! Flagging baselines of antenna 6:KK in the currently selected channels. ![Exited script file: /tmp/ff.ff] uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.8463 mas, bmaj=1.885 mas, bpa=-4.177 degrees ! Estimated noise=11.0648 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 = 3.19685 Jy ! Component: 100 - total flux cleaned = 4.10544 Jy ! Total flux subtracted in 100 components = 4.10544 Jy ! Clean residual min=-0.255939 max=0.727029 Jy/beam ! Clean residual mean=-0.000748 rms=0.076378 Jy/beam ! Combined flux in latest and established models = 4.10544 Jy ! Performing phase self-cal ! Adding 3 model components to the UV plane model. ! The established model now contains 3 components and 4.10544 Jy ! ! Correcting IF 1. ! A total of 3 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 2 telescope corrections were flagged in sub-array 5. ! A total of 2 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 3 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 2 telescope corrections were flagged in sub-array 5. ! A total of 2 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 3 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 2 telescope corrections were flagged in sub-array 5. ! A total of 2 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 3 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 2 telescope corrections were flagged in sub-array 5. ! A total of 2 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=1.674309Jy sigma=3.758741 ! Fit after self-cal, rms=1.674082Jy sigma=3.758482 ! Inverting map and beam ! Estimated beam: bmin=0.8463 mas, bmaj=1.885 mas, bpa=-4.177 degrees ! Estimated noise=11.0648 mJy/beam. ! Added new window around map position (-1, -1.6). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.655649 Jy ! Component: 100 - total flux cleaned = 1.03631 Jy ! Total flux subtracted in 100 components = 1.03631 Jy ! Clean residual min=-0.262013 max=0.550055 Jy/beam ! Clean residual mean=-0.000408 rms=0.062189 Jy/beam ! Combined flux in latest and established models = 5.14175 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 11 components and 5.14175 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.277722Jy sigma=3.029977 ! Fit after self-cal, rms=1.264933Jy sigma=2.997698 ! Inverting map ! Added new window around map position (-2.4, -3.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.456722 Jy ! Component: 100 - total flux cleaned = 0.726008 Jy ! Total flux subtracted in 100 components = 0.726008 Jy ! Clean residual min=-0.267516 max=0.384186 Jy/beam ! Clean residual mean=-0.000223 rms=0.052144 Jy/beam ! Combined flux in latest and established models = 5.86776 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 18 components and 5.86776 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.101487Jy sigma=2.582925 ! Fit after self-cal, rms=1.073573Jy sigma=2.521818 ! Inverting map ! Added new window around map position (1.2, 1.6). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.292881 Jy ! Component: 100 - total flux cleaned = 0.405235 Jy ! Total flux subtracted in 100 components = 0.405235 Jy ! Clean residual min=-0.253734 max=0.344943 Jy/beam ! Clean residual mean=-0.000101 rms=0.046220 Jy/beam ! Combined flux in latest and established models = 6.27299 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 26 components and 6.27299 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.974277Jy sigma=2.313301 ! Fit after self-cal, rms=0.969164Jy sigma=2.297036 ! Inverting map ! Added new window around map position (1.6, 5.6). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.27538 Jy ! Component: 100 - total flux cleaned = 0.387884 Jy ! Total flux subtracted in 100 components = 0.387884 Jy ! Clean residual min=-0.244236 max=0.234353 Jy/beam ! Clean residual mean=-0.000052 rms=0.040687 Jy/beam ! Combined flux in latest and established models = 6.66088 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 35 components and 6.66088 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.898287Jy sigma=2.095928 ! Fit after self-cal, rms=0.891549Jy sigma=2.074345 ! Inverting map print "=== Finished uniform weighting clean ===" ! === Finished uniform weighting clean === uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.9828 mas, bmaj=2.219 mas, bpa=-6.719 degrees ! Estimated noise=9.78659 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.0747992 Jy ! Component: 100 - total flux cleaned = 0.113253 Jy ! Component: 150 - total flux cleaned = 0.141823 Jy ! Component: 200 - total flux cleaned = 0.167058 Jy ! Total flux subtracted in 200 components = 0.167058 Jy ! Clean residual min=-0.262041 max=0.216794 Jy/beam ! Clean residual mean=-0.000096 rms=0.043762 Jy/beam ! Combined flux in latest and established models = 6.82794 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 44 components and 6.82794 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.845252Jy sigma=1.928666 ! Fit after self-cal, rms=0.845186Jy sigma=1.925154 ! 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: ! SC 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.82 HN 1.25 LA 1.23 NL 1.19 ! PT 1.10 SC 0.88 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 0.94 HN 1.31 KP 0.89 LA 1.03 ! NL 0.96 PT 1.00 ! ! Telescope amplitude corrections in sub-array 4: ! FD 0.90 HN 1.36 KP 0.84 LA 0.98 ! NL 0.96 OV 1.03 PT 1.00 SC 0.88 ! WZ 1.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 0.89 FD 0.98 HN 1.49 KK 1.00* ! KP 0.87 LA 0.99 MK 0.96 NL 0.97 ! OV 0.95 PT 0.96 SC 0.97 ! ! Telescope amplitude corrections in sub-array 6: ! BR 0.96 FD 0.98 HN 1.37 KK 1.00* ! KP 0.95 LA 1.11 MK 0.87 OV 0.94 ! PT 1.00 SC 0.93 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! SC 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.80 HN 1.29 LA 1.22 NL 1.30 ! PT 1.12 SC 0.94 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 0.95 HN 1.32 KP 0.82 LA 1.03 ! NL 0.94 PT 1.09 ! ! Telescope amplitude corrections in sub-array 4: ! FD 0.91 HN 1.32 KP 0.79 LA 1.04 ! NL 0.98 OV 0.96 PT 1.06 SC 0.90 ! WZ 1.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 0.87 FD 1.01 HN 1.41 KK 1.00* ! KP 0.82 LA 0.99 MK 0.97 NL 1.03 ! OV 0.88 PT 1.06 SC 1.00 ! ! Telescope amplitude corrections in sub-array 6: ! BR 0.96 FD 0.94 HN 1.41 KK 1.00* ! KP 0.94 LA 1.09 MK 0.97 OV 0.96 ! PT 1.02 SC 0.92 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! SC 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.82 HN 1.30 LA 1.26 NL 1.18 ! PT 1.18 SC 0.94 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 1.02 HN 1.47 KP 0.83 LA 1.03 ! NL 0.91 PT 1.04 ! ! Telescope amplitude corrections in sub-array 4: ! FD 0.94 HN 1.42 KP 0.74 LA 1.02 ! NL 0.98 OV 1.01 PT 1.04 SC 0.90 ! WZ 1.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 0.97 FD 0.98 HN 1.39 KK 1.00* ! KP 0.81 LA 1.07 MK 0.97 NL 0.92 ! OV 0.94 PT 0.97 SC 1.03 ! ! Telescope amplitude corrections in sub-array 6: ! BR 1.01 FD 0.96 HN 1.36 KK 1.00* ! KP 0.87 LA 1.11 MK 0.93 OV 0.99 ! PT 1.05 SC 0.98 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! SC 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.85 HN 1.33 LA 1.31 NL 1.14 ! PT 1.14 SC 0.95 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 1.01 HN 1.43 KP 0.81 LA 1.12 ! NL 0.95 PT 1.08 ! ! Telescope amplitude corrections in sub-array 4: ! FD 0.99 HN 1.41 KP 0.76 LA 1.01 ! NL 1.00 OV 0.99 PT 1.05 SC 0.92 ! WZ 1.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 0.92 FD 1.01 HN 1.37 KK 1.00* ! KP 0.81 LA 1.02 MK 0.99 NL 0.96 ! OV 1.01 PT 0.95 SC 1.01 ! ! Telescope amplitude corrections in sub-array 6: ! BR 0.96 FD 1.00 HN 1.47 KK 1.00* ! KP 0.91 LA 1.13 MK 0.99 OV 0.93 ! PT 1.01 SC 1.01 ! ! ! Fit before self-cal, rms=0.845186Jy sigma=1.925154 ! Fit after self-cal, rms=0.597348Jy sigma=1.060034 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.997 mas, bmaj=2.197 mas, bpa=-5.41 degrees ! Estimated noise=10.1042 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.172979 Jy ! Component: 100 - total flux cleaned = -0.227144 Jy ! Component: 150 - total flux cleaned = -0.266291 Jy ! Component: 200 - total flux cleaned = -0.300588 Jy ! Total flux subtracted in 200 components = -0.300588 Jy ! Clean residual min=-0.149174 max=0.077062 Jy/beam ! Clean residual mean=-0.000441 rms=0.019361 Jy/beam ! Combined flux in latest and established models = 6.52735 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 56 components and 6.52735 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.511489Jy sigma=0.929917 ! Fit after self-cal, rms=0.507434Jy sigma=0.923372 ! 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.507434Jy sigma=0.923372 ! Fit after self-cal, rms=0.502566Jy sigma=0.869158 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.001 mas, bmaj=2.197 mas, bpa=-4.973 degrees ! Estimated noise=10.1587 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.125684 Jy ! Component: 100 - total flux cleaned = -0.175937 Jy ! Component: 150 - total flux cleaned = -0.209881 Jy ! Component: 200 - total flux cleaned = -0.235194 Jy ! Total flux subtracted in 200 components = -0.235194 Jy ! Clean residual min=-0.129921 max=0.068296 Jy/beam ! Clean residual mean=-0.000373 rms=0.016479 Jy/beam ! Combined flux in latest and established models = 6.29215 Jy selfcal ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 63 components and 6.29215 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.446705Jy sigma=0.791430 ! Fit after self-cal, rms=0.444498Jy sigma=0.787253 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.0265268 Jy ! Component: 100 - total flux cleaned = -0.0329491 Jy ! Component: 150 - total flux cleaned = -0.0329694 Jy ! Component: 200 - total flux cleaned = -0.0355013 Jy ! Total flux subtracted in 200 components = -0.0355013 Jy ! Clean residual min=-0.123268 max=0.067318 Jy/beam ! Clean residual mean=-0.000375 rms=0.016028 Jy/beam ! Combined flux in latest and established models = 6.25665 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 74 components and 6.25665 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.432334Jy sigma=0.770502 ! Fit after self-cal, rms=0.431576Jy sigma=0.769022 ! 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.431576Jy sigma=0.769022 ! Fit after self-cal, rms=0.423698Jy sigma=0.713959 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.005 mas, bmaj=2.198 mas, bpa=-4.577 degrees ! Estimated noise=10.2262 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0926523 Jy ! Component: 100 - total flux cleaned = -0.134664 Jy ! Component: 150 - total flux cleaned = -0.157588 Jy ! Component: 200 - total flux cleaned = -0.176422 Jy ! Total flux subtracted in 200 components = -0.176422 Jy ! Clean residual min=-0.109019 max=0.058012 Jy/beam ! Clean residual mean=-0.000308 rms=0.013832 Jy/beam ! Combined flux in latest and established models = 6.08023 Jy selfcal ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 77 components and 6.08023 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.387253Jy sigma=0.665575 ! Fit after self-cal, rms=0.386134Jy sigma=0.663477 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.0180952 Jy ! Component: 100 - total flux cleaned = -0.0241151 Jy ! Component: 150 - total flux cleaned = -0.0283721 Jy ! Component: 200 - total flux cleaned = -0.032295 Jy ! Total flux subtracted in 200 components = -0.032295 Jy ! Clean residual min=-0.103359 max=0.057329 Jy/beam ! Clean residual mean=-0.000310 rms=0.013528 Jy/beam ! Combined flux in latest and established models = 6.04794 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 82 components and 6.04794 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.378341Jy sigma=0.652172 ! Fit after self-cal, rms=0.377719Jy sigma=0.651232 ! 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.377719Jy sigma=0.651232 ! Fit after self-cal, rms=0.370689Jy sigma=0.613933 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.138 mas, bmaj=3.567 mas, bpa=-10.49 degrees ! Estimated noise=12.4609 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.080301 Jy ! Component: 100 - total flux cleaned = -0.128971 Jy ! Component: 150 - total flux cleaned = -0.162187 Jy ! Component: 200 - total flux cleaned = -0.178727 Jy ! Total flux subtracted in 200 components = -0.178727 Jy ! Clean residual min=-0.094423 max=0.055857 Jy/beam ! Clean residual mean=-0.000493 rms=0.014926 Jy/beam ! Combined flux in latest and established models = 5.86921 Jy selfcal ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 101 components and 5.86921 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.354031Jy sigma=0.598445 ! Fit after self-cal, rms=0.353066Jy sigma=0.596271 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.0128807 Jy ! Component: 100 - total flux cleaned = -0.0173784 Jy ! Component: 150 - total flux cleaned = -0.0213086 Jy ! Component: 200 - total flux cleaned = -0.0236524 Jy ! Total flux subtracted in 200 components = -0.0236524 Jy ! Clean residual min=-0.090935 max=0.055034 Jy/beam ! Clean residual mean=-0.000492 rms=0.014458 Jy/beam ! Combined flux in latest and established models = 5.84556 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 105 components and 5.84556 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.350687Jy sigma=0.595086 ! Fit after self-cal, rms=0.350275Jy sigma=0.594398 ! 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:SC,true ! 1:SC self-cal status=fixed weight=1 selfant 2:FD,true ! 2:FD self-cal status=fixed weight=1 selfant 2:HN,true ! 2:HN self-cal status=fixed weight=1 selfant 2:LA,true ! 2:LA self-cal status=fixed weight=1 selfant 2:NL,true ! 2:NL self-cal status=fixed weight=1 selfant 2:PT,true ! 2:PT self-cal status=fixed weight=1 selfant 2:SC,true ! 2:SC self-cal status=fixed weight=1 selfant 3:FD,true ! 3:FD self-cal status=fixed weight=1 selfant 3:HN,true ! 3:HN self-cal status=fixed weight=1 selfant 3:KP,true ! 3:KP self-cal status=fixed weight=1 selfant 3:LA,true ! 3:LA self-cal status=fixed weight=1 selfant 3:NL,true ! 3:NL self-cal status=fixed weight=1 selfant 3:PT,true ! 3:PT self-cal status=fixed weight=1 selfant 4: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:OV,true ! 4:OV self-cal status=fixed weight=1 selfant 4:PT,true ! 4:PT self-cal status=fixed weight=1 selfant 4:SC,true ! 4:SC self-cal status=fixed weight=1 selfant 5:BR,true ! 5:BR self-cal status=fixed weight=1 selfant 5:FD,true ! 5:FD self-cal status=fixed weight=1 selfant 5:HN,true ! 5:HN self-cal status=fixed weight=1 selfant 5:KP,true ! 5:KP self-cal status=fixed weight=1 selfant 5:LA,true ! 5:LA self-cal status=fixed weight=1 selfant 5:MK,true ! 5:MK self-cal status=fixed weight=1 selfant 5:NL,true ! 5:NL self-cal status=fixed weight=1 selfant 5:OV,true ! 5:OV self-cal status=fixed weight=1 selfant 5:PT,true ! 5:PT self-cal status=fixed weight=1 selfant 5:SC,true ! 5:SC self-cal status=fixed weight=1 selfant 6:BR,true ! 6:BR self-cal status=fixed weight=1 selfant 6:FD,true ! 6:FD self-cal status=fixed weight=1 selfant 6:HN,true ! 6:HN self-cal status=fixed weight=1 selfant 6:KP,true ! 6:KP self-cal status=fixed weight=1 selfant 6:LA,true ! 6:LA self-cal status=fixed weight=1 selfant 6:MK,true ! 6:MK self-cal status=fixed weight=1 selfant 6: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 4 telescope corrections were flagged in sub-array 1. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 1 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 4 telescope corrections were flagged in sub-array 1. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 1 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 4 telescope corrections were flagged in sub-array 1. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 1 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 4 telescope corrections were flagged in sub-array 1. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 1 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.820039Jy sigma=2.050005 ! Fit after self-cal, rms=0.467429Jy sigma=0.869998 selfant "",false ! 1:SC self-cal status=correctable weight=1 ! 1:WZ self-cal status=correctable weight=1 ! 1:ZC self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:HN self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 3:FD self-cal status=correctable weight=1 ! 3:HN self-cal status=correctable weight=1 ! 3:KP self-cal status=correctable weight=1 ! 3:LA self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 3:PT self-cal status=correctable weight=1 ! 4:FD self-cal status=correctable weight=1 ! 4:HN self-cal status=correctable weight=1 ! 4:KP self-cal status=correctable weight=1 ! 4:LA self-cal status=correctable weight=1 ! 4:NL self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 4:PT self-cal status=correctable weight=1 ! 4:SC self-cal status=correctable weight=1 ! 4:WZ self-cal status=correctable weight=1 ! 5:BR self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:HN self-cal status=correctable weight=1 ! 5:KK self-cal status=correctable weight=1 ! 5:KP self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:MK self-cal status=correctable weight=1 ! 5:NL self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 5:SC self-cal status=correctable weight=1 ! 6:BR self-cal status=correctable weight=1 ! 6:FD self-cal status=correctable weight=1 ! 6:HN self-cal status=correctable weight=1 ! 6:KK self-cal status=correctable weight=1 ! 6:KP self-cal status=correctable weight=1 ! 6:LA self-cal status=correctable weight=1 ! 6:MK self-cal status=correctable weight=1 ! 6: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 2 telescope corrections were flagged in sub-array 1. ! A total of 5 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 2 telescope corrections were flagged in sub-array 1. ! A total of 5 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 2 telescope corrections were flagged in sub-array 1. ! A total of 5 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 2 telescope corrections were flagged in sub-array 1. ! A total of 5 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.467429Jy sigma=0.869998 ! Fit after self-cal, rms=0.388143Jy sigma=0.595863 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.6208 mas, bmaj=1.955 mas, bpa=-6.485 degrees ! Estimated noise=10.8588 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.0593555 Jy ! Component: 100 - total flux cleaned = -0.0776472 Jy ! Component: 150 - total flux cleaned = -0.0888124 Jy ! Component: 200 - total flux cleaned = -0.0938618 Jy ! Total flux subtracted in 200 components = -0.0938618 Jy ! Clean residual min=-0.073268 max=0.051300 Jy/beam ! Clean residual mean=-0.000069 rms=0.010828 Jy/beam ! Combined flux in latest and established models = 5.75169 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 116 components and 5.75169 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.357253Jy sigma=0.560776 ! Fit after self-cal, rms=0.355882Jy sigma=0.556195 ! 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.7245 mas, bmaj=2.248 mas, bpa=-7.547 degrees ! Estimated noise=9.48163 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.0131975 Jy ! Component: 100 - total flux cleaned = -0.0255382 Jy ! Component: 150 - total flux cleaned = -0.0350532 Jy ! Component: 200 - total flux cleaned = -0.0393299 Jy ! Component: 250 - total flux cleaned = -0.0422942 Jy ! Component: 300 - total flux cleaned = -0.044148 Jy ! Component: 350 - total flux cleaned = -0.0450191 Jy ! Component: 400 - total flux cleaned = -0.0458536 Jy ! Total flux subtracted in 400 components = -0.0458536 Jy ! Clean residual min=-0.066829 max=0.042161 Jy/beam ! Clean residual mean=-0.000137 rms=0.009710 Jy/beam ! Combined flux in latest and established models = 5.70584 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 126 components and 5.70584 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.337972Jy sigma=0.532778 ! Fit after self-cal, rms=0.337388Jy sigma=0.531711 ! 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.00549456 Jy ! Component: 100 - total flux cleaned = -0.00389578 Jy ! Component: 150 - total flux cleaned = -0.00315673 Jy ! Component: 200 - total flux cleaned = -0.00317233 Jy ! Component: 250 - total flux cleaned = -0.00318735 Jy ! Component: 300 - total flux cleaned = -0.00385188 Jy ! Component: 350 - total flux cleaned = -0.003847 Jy ! Component: 400 - total flux cleaned = -0.00323316 Jy ! Total flux subtracted in 400 components = -0.00323316 Jy ! Clean residual min=-0.064056 max=0.040532 Jy/beam ! Clean residual mean=-0.000137 rms=0.009514 Jy/beam ! Combined flux in latest and established models = 5.70261 Jy ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 139 components and 5.70261 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.330205Jy sigma=0.522125 ! Fit after self-cal, rms=0.329757Jy sigma=0.521599 ! 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 6 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 5. ! A total of 1 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 6 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 5. ! A total of 1 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 6 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 5. ! A total of 1 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 6 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 5. ! A total of 1 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.783866Jy sigma=2.000813 ! Fit after self-cal, rms=0.317248Jy sigma=0.489880 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.7253 mas, bmaj=2.244 mas, bpa=-7.423 degrees ! Estimated noise=9.50965 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0479899 Jy ! Component: 100 - total flux cleaned = -0.0639937 Jy ! Component: 150 - total flux cleaned = -0.0748929 Jy ! Component: 200 - total flux cleaned = -0.0823765 Jy ! Component: 250 - total flux cleaned = -0.0840379 Jy ! Component: 300 - total flux cleaned = -0.0848119 Jy ! Component: 350 - total flux cleaned = -0.0862419 Jy ! Component: 400 - total flux cleaned = -0.0875722 Jy ! Total flux subtracted in 400 components = -0.0875722 Jy ! Clean residual min=-0.057278 max=0.033849 Jy/beam ! Clean residual mean=-0.000096 rms=0.008383 Jy/beam ! Combined flux in latest and established models = 5.61503 Jy selfcal ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 142 components and 5.61504 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.297943Jy sigma=0.465392 ! Fit after self-cal, rms=0.297262Jy sigma=0.464580 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.00303684 Jy ! Component: 100 - total flux cleaned = -0.0036489 Jy ! Component: 150 - total flux cleaned = -0.00364373 Jy ! Component: 200 - total flux cleaned = -0.00421756 Jy ! Component: 250 - total flux cleaned = -0.00529178 Jy ! Component: 300 - total flux cleaned = -0.00475651 Jy ! Component: 350 - total flux cleaned = -0.00525736 Jy ! Component: 400 - total flux cleaned = -0.00574772 Jy ! Total flux subtracted in 400 components = -0.00574772 Jy ! Clean residual min=-0.055629 max=0.032941 Jy/beam ! Clean residual mean=-0.000096 rms=0.008234 Jy/beam ! Combined flux in latest and established models = 5.60929 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 146 components and 5.60929 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.292443Jy sigma=0.457216 ! Fit after self-cal, rms=0.292017Jy sigma=0.456864 ! 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.292017Jy sigma=0.456864 ! Fit after self-cal, rms=0.285902Jy sigma=0.440803 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.7251 mas, bmaj=2.239 mas, bpa=-7.346 degrees ! Estimated noise=9.5248 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0366175 Jy ! Component: 100 - total flux cleaned = -0.0479312 Jy ! Component: 150 - total flux cleaned = -0.0564363 Jy ! Component: 200 - total flux cleaned = -0.0622522 Jy ! Component: 250 - total flux cleaned = -0.0622469 Jy ! Component: 300 - total flux cleaned = -0.062862 Jy ! Component: 350 - total flux cleaned = -0.0639851 Jy ! Component: 400 - total flux cleaned = -0.0639918 Jy ! Total flux subtracted in 400 components = -0.0639918 Jy ! Clean residual min=-0.050257 max=0.030639 Jy/beam ! Clean residual mean=-0.000073 rms=0.007550 Jy/beam ! Combined flux in latest and established models = 5.5453 Jy selfcal ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 148 components and 5.54529 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.273611Jy sigma=0.424041 ! Fit after self-cal, rms=0.273041Jy sigma=0.423454 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.00248776 Jy ! Component: 100 - total flux cleaned = -0.00248274 Jy ! Component: 150 - total flux cleaned = -0.0034046 Jy ! Component: 200 - total flux cleaned = -0.00384378 Jy ! Component: 250 - total flux cleaned = -0.00426435 Jy ! Component: 300 - total flux cleaned = -0.00468515 Jy ! Component: 350 - total flux cleaned = -0.00548876 Jy ! Component: 400 - total flux cleaned = -0.00548062 Jy ! Total flux subtracted in 400 components = -0.00548062 Jy ! Clean residual min=-0.048888 max=0.029815 Jy/beam ! Clean residual mean=-0.000073 rms=0.007446 Jy/beam ! Combined flux in latest and established models = 5.53981 Jy ! Performing phase self-cal ! Adding 39 model components to the UV plane model. ! The established model now contains 155 components and 5.53981 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.269922Jy sigma=0.418460 ! Fit after self-cal, rms=0.269572Jy sigma=0.418211 ! 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.269572Jy sigma=0.418211 ! Fit after self-cal, rms=0.263816Jy sigma=0.405727 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.157 mas, bmaj=3.599 mas, bpa=-10.2 degrees ! Estimated noise=12.4127 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0316192 Jy ! Component: 100 - total flux cleaned = -0.0418504 Jy ! Component: 150 - total flux cleaned = -0.0492364 Jy ! Component: 200 - total flux cleaned = -0.0541671 Jy ! Component: 250 - total flux cleaned = -0.057413 Jy ! Component: 300 - total flux cleaned = -0.0590432 Jy ! Component: 350 - total flux cleaned = -0.0594039 Jy ! Component: 400 - total flux cleaned = -0.0597538 Jy ! Total flux subtracted in 400 components = -0.0597538 Jy ! Clean residual min=-0.047509 max=0.034507 Jy/beam ! Clean residual mean=-0.000118 rms=0.009543 Jy/beam ! Combined flux in latest and established models = 5.48006 Jy selfcal ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 163 components and 5.48006 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.260884Jy sigma=0.404450 ! Fit after self-cal, rms=0.260629Jy sigma=0.403183 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.00262086 Jy ! Component: 100 - total flux cleaned = -0.00501305 Jy ! Component: 150 - total flux cleaned = -0.00560334 Jy ! Component: 200 - total flux cleaned = -0.00613805 Jy ! Component: 250 - total flux cleaned = -0.00663128 Jy ! Component: 300 - total flux cleaned = -0.00663587 Jy ! Component: 350 - total flux cleaned = -0.00663858 Jy ! Component: 400 - total flux cleaned = -0.00663897 Jy ! Total flux subtracted in 400 components = -0.00663897 Jy ! Clean residual min=-0.046722 max=0.034047 Jy/beam ! Clean residual mean=-0.000118 rms=0.009348 Jy/beam ! Combined flux in latest and established models = 5.47342 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 168 components and 5.47342 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.261788Jy sigma=0.407326 ! Fit after self-cal, rms=0.261774Jy sigma=0.406895 ! 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.7245 mas, bmaj=2.235 mas, bpa=-7.288 degrees ! Estimated noise=9.53527 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.7245 x 2.235 at -7.288 degrees (North through East) ! Clean map min=-0.10865 max=4.1461 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.261774Jy sigma=0.406895 ! Fit after self-cal, rms=0.251649Jy sigma=0.382352 ! Inverting map and beam ! Estimated beam: bmin=0.7241 mas, bmaj=2.234 mas, bpa=-7.26 degrees ! Estimated noise=9.5407 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0264966 Jy ! Component: 100 - total flux cleaned = -0.0288829 Jy ! Component: 150 - total flux cleaned = -0.0289136 Jy ! Component: 200 - total flux cleaned = -0.0289266 Jy ! Component: 250 - total flux cleaned = -0.0283665 Jy ! Component: 300 - total flux cleaned = -0.02837 Jy ! Component: 350 - total flux cleaned = -0.0278694 Jy ! Component: 400 - total flux cleaned = -0.0274094 Jy ! Total flux subtracted in 400 components = -0.0274094 Jy ! Clean residual min=-0.042438 max=0.028773 Jy/beam ! Clean residual mean=-0.000029 rms=0.006605 Jy/beam ! Combined flux in latest and established models = 5.44601 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 170 components and 5.44601 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.243548Jy sigma=0.367836 ! Fit after self-cal, rms=0.243032Jy sigma=0.367245 ! Inverting map ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000358124 Jy ! Component: 100 - total flux cleaned = 0.000559088 Jy ! Component: 150 - total flux cleaned = 0.000557014 Jy ! Component: 200 - total flux cleaned = 0.000558798 Jy ! Component: 250 - total flux cleaned = 0.000165341 Jy ! Component: 300 - total flux cleaned = 0.000169322 Jy ! Component: 350 - total flux cleaned = 0.000164673 Jy ! Component: 400 - total flux cleaned = -0.000194797 Jy ! Total flux subtracted in 400 components = -0.000194797 Jy ! Clean residual min=-0.041261 max=0.027752 Jy/beam ! Clean residual mean=-0.000029 rms=0.006512 Jy/beam ! Combined flux in latest and established models = 5.44582 Jy ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 172 components and 5.44582 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.240167Jy sigma=0.362275 ! Fit after self-cal, rms=0.239757Jy sigma=0.362030 ! 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.239757Jy sigma=0.362030 ! Fit after self-cal, rms=0.235328Jy sigma=0.355264 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.6185 mas, bmaj=1.936 mas, bpa=-6.145 degrees ! Estimated noise=10.9307 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 3.38632 Jy ! Component: 100 - total flux cleaned = 4.15094 Jy ! Component: 150 - total flux cleaned = 4.48627 Jy ! Component: 200 - total flux cleaned = 4.70207 Jy ! Component: 250 - total flux cleaned = 4.85581 Jy ! Component: 300 - total flux cleaned = 4.97082 Jy ! Component: 350 - total flux cleaned = 5.06137 Jy ! Component: 400 - total flux cleaned = 5.12645 Jy ! Total flux subtracted in 400 components = 5.12645 Jy ! Clean residual min=-0.048765 max=0.069916 Jy/beam ! Clean residual mean=-0.000072 rms=0.010304 Jy/beam ! Combined flux in latest and established models = 5.12644 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 41 model components to the UV plane model. ! The established model now contains 41 components and 5.12644 Jy ! Inverting map and beam ! Estimated beam: bmin=0.7231 mas, bmaj=2.231 mas, bpa=-7.235 degrees ! Estimated noise=9.5447 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 (-1.2, 1). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0933685 Jy ! Component: 100 - total flux cleaned = 0.137366 Jy ! Component: 150 - total flux cleaned = 0.16024 Jy ! Component: 200 - total flux cleaned = 0.18206 Jy ! Component: 250 - total flux cleaned = 0.196727 Jy ! Component: 300 - total flux cleaned = 0.210065 Jy ! Component: 350 - total flux cleaned = 0.220947 Jy ! Component: 400 - total flux cleaned = 0.23097 Jy ! Total flux subtracted in 400 components = 0.23097 Jy ! Clean residual min=-0.035252 max=0.036090 Jy/beam ! Clean residual mean=-0.000032 rms=0.007114 Jy/beam ! Combined flux in latest and established models = 5.35741 Jy ! Adding 48 model components to the UV plane model. ! The established model now contains 82 components and 5.35741 Jy ! Inverting map ! Added new window around map position (2.4, 3). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00209451 Jy ! Component: 100 - total flux cleaned = 0.00327862 Jy ! Component: 150 - total flux cleaned = 0.00321564 Jy ! Component: 200 - total flux cleaned = 0.00428719 Jy ! Component: 250 - total flux cleaned = 0.00636119 Jy ! Component: 300 - total flux cleaned = 0.00828569 Jy ! Component: 350 - total flux cleaned = 0.00828124 Jy ! Component: 400 - total flux cleaned = 0.00829117 Jy ! Total flux subtracted in 400 components = 0.00829117 Jy ! Clean residual min=-0.028846 max=0.024541 Jy/beam ! Clean residual mean=-0.000027 rms=0.006036 Jy/beam ! Combined flux in latest and established models = 5.3657 Jy ! Adding 41 model components to the UV plane model. ! The established model now contains 113 components and 5.3657 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.166 mas, bmaj=3.602 mas, bpa=-9.942 degrees ! Estimated noise=12.4242 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0109991 Jy ! Component: 100 - total flux cleaned = -0.0111178 Jy ! Component: 150 - total flux cleaned = -0.0103166 Jy ! Component: 200 - total flux cleaned = -0.00818231 Jy ! Component: 250 - total flux cleaned = -0.00693221 Jy ! Component: 300 - total flux cleaned = -0.0063553 Jy ! Component: 350 - total flux cleaned = -0.00634886 Jy ! Component: 400 - total flux cleaned = -0.00632515 Jy ! Total flux subtracted in 400 components = -0.00632515 Jy ! Clean residual min=-0.032322 max=0.031597 Jy/beam ! Clean residual mean=-0.000050 rms=0.007880 Jy/beam ! Combined flux in latest and established models = 5.35938 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 = 6.08638e-06 Jy ! Component: 100 - total flux cleaned = 0.000449087 Jy ! Component: 150 - total flux cleaned = 0.000460808 Jy ! Component: 200 - total flux cleaned = 0.000479214 Jy ! Component: 250 - total flux cleaned = 0.000488288 Jy ! Component: 300 - total flux cleaned = 0.000494286 Jy ! Component: 350 - total flux cleaned = 0.000498115 Jy ! Component: 400 - total flux cleaned = 0.000498795 Jy ! Total flux subtracted in 400 components = 0.000498795 Jy ! Clean residual min=-0.030224 max=0.031250 Jy/beam ! Clean residual mean=-0.000050 rms=0.007729 Jy/beam ! Combined flux in latest and established models = 5.35988 Jy ! Adding 67 model components to the UV plane model. ! The established model now contains 169 components and 5.35988 Jy ! Inverting map ! Added new window around map position (-1.6, -7.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000363377 Jy ! Component: 100 - total flux cleaned = 0.000438727 Jy ! Component: 150 - total flux cleaned = 0.000344125 Jy ! Component: 200 - total flux cleaned = 0.000927755 Jy ! Component: 250 - total flux cleaned = 0.00147633 Jy ! Component: 300 - total flux cleaned = 0.000433535 Jy ! Component: 350 - total flux cleaned = -0.000524288 Jy ! Component: 400 - total flux cleaned = -0.000527498 Jy ! Total flux subtracted in 400 components = -0.000527498 Jy ! Clean residual min=-0.024676 max=0.024075 Jy/beam ! Clean residual mean=-0.000054 rms=0.006949 Jy/beam ! Combined flux in latest and established models = 5.35935 Jy ! Adding 51 model components to the UV plane model. ! The established model now contains 217 components and 5.35935 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=0.7231 mas, bmaj=2.231 mas, bpa=-7.235 degrees ! Estimated noise=9.5447 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00073202 Jy ! Component: 100 - total flux cleaned = -0.000714358 Jy ! Component: 150 - total flux cleaned = -0.00225816 Jy ! Component: 200 - total flux cleaned = -0.00439883 Jy ! Component: 250 - total flux cleaned = -0.00640681 Jy ! Component: 300 - total flux cleaned = -0.00703594 Jy ! Component: 350 - total flux cleaned = -0.00764921 Jy ! Component: 400 - total flux cleaned = -0.00824752 Jy ! Total flux subtracted in 400 components = -0.00824752 Jy ! Clean residual min=-0.020637 max=0.018783 Jy/beam ! Clean residual mean=-0.000023 rms=0.004817 Jy/beam ! Combined flux in latest and established models = 5.3511 Jy selfcal ! Performing phase self-cal ! Adding 59 model components to the UV plane model. ! The established model now contains 268 components and 5.3511 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.201895Jy sigma=0.301511 ! Fit after self-cal, rms=0.195274Jy sigma=0.292710 wmodel /scr/pima/rdv73_uvs/J0730-1141_X_map.mod ! Writing 268 model components to file: /scr/pima/rdv73_uvs/J0730-1141_X_map.mod wobs /scr/pima/rdv73_uvs/J0730-1141_X_uvs.fits ! Applying 1030 buffered edits. ! Writing UV FITS file: /scr/pima/rdv73_uvs/J0730-1141_X_uvs.fits wwins /scr/pima/rdv73_uvs/J0730-1141_X_map.win ! wwins: Wrote 8 windows to /scr/pima/rdv73_uvs/J0730-1141_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.00640779 Jy ! Component: 100 - total flux cleaned = -0.0071055 Jy ! Component: 150 - total flux cleaned = -0.00810684 Jy ! Component: 200 - total flux cleaned = -0.00875247 Jy ! Component: 250 - total flux cleaned = -0.00719386 Jy ! Component: 300 - total flux cleaned = -0.00659365 Jy ! Component: 350 - total flux cleaned = -0.00659434 Jy ! Component: 400 - total flux cleaned = -0.00630042 Jy ! Component: 450 - total flux cleaned = -0.00572553 Jy ! Component: 500 - total flux cleaned = -0.00600798 Jy ! Component: 550 - total flux cleaned = -0.00573175 Jy ! Component: 600 - total flux cleaned = -0.00710026 Jy ! Component: 650 - total flux cleaned = -0.00818111 Jy ! Component: 700 - total flux cleaned = -0.00951233 Jy ! Component: 750 - total flux cleaned = -0.0100349 Jy ! Component: 800 - total flux cleaned = -0.0110702 Jy ! Component: 850 - total flux cleaned = -0.0118407 Jy ! Component: 900 - total flux cleaned = -0.0118405 Jy ! Component: 950 - total flux cleaned = -0.0125921 Jy ! Component: 1000 - total flux cleaned = -0.012592 Jy ! Component: 1050 - total flux cleaned = -0.0128383 Jy ! Component: 1100 - total flux cleaned = -0.0121081 Jy ! Component: 1150 - total flux cleaned = -0.0128305 Jy ! Component: 1200 - total flux cleaned = -0.0121191 Jy ! Component: 1250 - total flux cleaned = -0.012824 Jy ! Component: 1300 - total flux cleaned = -0.0114171 Jy ! Component: 1350 - total flux cleaned = -0.0121143 Jy ! Component: 1400 - total flux cleaned = -0.0116551 Jy ! Component: 1450 - total flux cleaned = -0.0116543 Jy ! Component: 1500 - total flux cleaned = -0.0118808 Jy ! Component: 1550 - total flux cleaned = -0.0125501 Jy ! Component: 1600 - total flux cleaned = -0.0121069 Jy ! Component: 1650 - total flux cleaned = -0.0114497 Jy ! Component: 1700 - total flux cleaned = -0.0110128 Jy ! Component: 1750 - total flux cleaned = -0.0125308 Jy ! Component: 1800 - total flux cleaned = -0.0129608 Jy ! Component: 1850 - total flux cleaned = -0.0133859 Jy ! Component: 1900 - total flux cleaned = -0.0133817 Jy ! Component: 1950 - total flux cleaned = -0.0135919 Jy ! Component: 2000 - total flux cleaned = -0.0140081 Jy ! Component: 2050 - total flux cleaned = -0.0135935 Jy ! Component: 2100 - total flux cleaned = -0.014619 Jy ! Component: 2150 - total flux cleaned = -0.0148249 Jy ! Component: 2200 - total flux cleaned = -0.0148276 Jy ! Component: 2250 - total flux cleaned = -0.015834 Jy ! Component: 2300 - total flux cleaned = -0.0162336 Jy ! Component: 2350 - total flux cleaned = -0.0166319 Jy ! Component: 2400 - total flux cleaned = -0.0174211 Jy ! Component: 2450 - total flux cleaned = -0.0183975 Jy ! Component: 2500 - total flux cleaned = -0.0180091 Jy ! Component: 2550 - total flux cleaned = -0.0185856 Jy ! Component: 2600 - total flux cleaned = -0.0183929 Jy ! Component: 2650 - total flux cleaned = -0.0195363 Jy ! Component: 2700 - total flux cleaned = -0.0200999 Jy ! Component: 2750 - total flux cleaned = -0.0195383 Jy ! Component: 2800 - total flux cleaned = -0.0197251 Jy ! Component: 2850 - total flux cleaned = -0.0199106 Jy ! Component: 2900 - total flux cleaned = -0.0197275 Jy ! Component: 2950 - total flux cleaned = -0.0197271 Jy ! Component: 3000 - total flux cleaned = -0.0202722 Jy ! Component: 3050 - total flux cleaned = -0.0193701 Jy ! Component: 3100 - total flux cleaned = -0.0195493 Jy ! Component: 3150 - total flux cleaned = -0.0195514 Jy ! Component: 3200 - total flux cleaned = -0.019906 Jy ! Component: 3250 - total flux cleaned = -0.0204337 Jy ! Component: 3300 - total flux cleaned = -0.0200847 Jy ! Component: 3350 - total flux cleaned = -0.0206077 Jy ! Component: 3400 - total flux cleaned = -0.0206072 Jy ! Component: 3450 - total flux cleaned = -0.0209524 Jy ! Component: 3500 - total flux cleaned = -0.0206121 Jy ! Component: 3550 - total flux cleaned = -0.0221422 Jy ! Component: 3600 - total flux cleaned = -0.0212953 Jy ! Component: 3650 - total flux cleaned = -0.021462 Jy ! Component: 3700 - total flux cleaned = -0.0222971 Jy ! Component: 3750 - total flux cleaned = -0.0219637 Jy ! Component: 3800 - total flux cleaned = -0.0219648 Jy ! Component: 3850 - total flux cleaned = -0.0221285 Jy ! Component: 3900 - total flux cleaned = -0.0222921 Jy ! Component: 3950 - total flux cleaned = -0.0237543 Jy ! Component: 4000 - total flux cleaned = -0.0244004 Jy ! Component: 4050 - total flux cleaned = -0.0242393 Jy ! Total flux subtracted in 4096 components = -0.0245594 Jy ! Clean residual min=-0.008371 max=0.009221 Jy/beam ! Clean residual mean=-0.000033 rms=0.002844 Jy/beam ! Combined flux in latest and established models = 5.32654 Jy keep ! Adding 346 model components to the UV plane model. ! The established model now contains 609 components and 5.32654 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv73_uvs/J0730-1141_X\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.7231 x 2.231 at -7.235 degrees (North through East) ! Clean map min=-0.055674 max=4.2356 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv73_uvs/J0730-1141_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= 4.2355 rms= 0.00283 Jy quit ! Quitting program ! Log file /scr/pima/rdv73_uvs/J0730-1141_X_uvs.log closed on Sat Jan 9 17:32:10 2016