! Started logfile: /scr/pima/rdv95_uvs/J1257-3155_X_uvs.log on Wed Nov 6 20:46:50 2013 obs /scr/pima/rdv95_uvs/J1257-3155_X_uva.fits ! Reading UV FITS file: /scr/pima/rdv95_uvs/J1257-3155_X_uva.fits ! AN table 1: 7 integrations on 45 of 45 possible baselines. ! AN table 2: 24 integrations on 3 of 3 possible baselines. ! AN table 3: 18 integrations on 6 of 6 possible baselines. ! AN table 4: 13 integrations on 3 of 3 possible baselines. ! AN table 5: 4 integrations on 1 of 1 possible baselines. ! AN table 6: 17 integrations on 15 of 15 possible baselines. ! AN table 7: 17 integrations on 28 of 28 possible baselines. ! AN table 8: 26 integrations on 21 of 21 possible baselines. ! AN table 9: 2 integrations on 28 of 28 possible baselines. ! Apparent sampling: 0.57349 visibilities/baseline/integration-bin. ! Found source: J1257-3155 ! ! 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 8e+06 1 8e+06 ! 02 2 8.47599e+09 8e+06 1 8e+06 ! 03 3 8.79099e+09 8e+06 1 8e+06 ! 04 4 8.89599e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 4292 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/rdv95_uvs/J1257-3155_X] integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 800000 float soltime2 soltime2 = 15 float soltime3 soltime3 = 1 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float old_peak float new_peak float flux_cutoff selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:KP self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:MK self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:OV self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 1:SC self-cal status=correctable weight=1 ! 2:FT self-cal status=correctable weight=1 ! 2:MA self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 3:FT self-cal status=correctable weight=1 ! 3:MA self-cal status=correctable weight=1 ! 3:SC self-cal status=correctable weight=1 ! 3:WZ self-cal status=correctable weight=1 ! 4:FT self-cal status=correctable weight=1 ! 4:MA self-cal status=correctable weight=1 ! 4:SC self-cal status=correctable weight=1 ! 5:FT self-cal status=correctable weight=1 ! 5:SC self-cal status=correctable weight=1 ! 6:FD self-cal status=correctable weight=1 ! 6:FT self-cal status=correctable weight=1 ! 6:HN self-cal status=correctable weight=1 ! 6:LA self-cal status=correctable weight=1 ! 6:PT self-cal status=correctable weight=1 ! 6:SC self-cal status=correctable weight=1 ! 7:FD self-cal status=correctable weight=1 ! 7:FT self-cal status=correctable weight=1 ! 7:HN self-cal status=correctable weight=1 ! 7:KP self-cal status=correctable weight=1 ! 7:LA self-cal status=correctable weight=1 ! 7:NL self-cal status=correctable weight=1 ! 7:PT self-cal status=correctable weight=1 ! 7:SC self-cal status=correctable weight=1 ! 8:FD self-cal status=correctable weight=1 ! 8:FT self-cal status=correctable weight=1 ! 8:HN self-cal status=correctable weight=1 ! 8:KP self-cal status=correctable weight=1 ! 8:LA self-cal status=correctable weight=1 ! 8:NL self-cal status=correctable weight=1 ! 8:PT self-cal status=correctable weight=1 ! 9:FD self-cal status=correctable weight=1 ! 9:HN self-cal status=correctable weight=1 ! 9:KP self-cal status=correctable weight=1 ! 9:LA self-cal status=correctable weight=1 ! 9:NL self-cal status=correctable weight=1 ! 9:OV self-cal status=correctable weight=1 ! 9:PT self-cal status=correctable weight=1 ! 9: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 40 telescope corrections were flagged in sub-array 1. ! A total of 72 telescope corrections were flagged in sub-array 2. ! A total of 39 telescope corrections were flagged in sub-array 3. ! A total of 39 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! A total of 54 telescope corrections were flagged in sub-array 6. ! A total of 50 telescope corrections were flagged in sub-array 7. ! A total of 82 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 2. ! A total of 40 telescope corrections were flagged in sub-array 1. ! A total of 72 telescope corrections were flagged in sub-array 2. ! A total of 39 telescope corrections were flagged in sub-array 3. ! A total of 39 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! A total of 54 telescope corrections were flagged in sub-array 6. ! A total of 50 telescope corrections were flagged in sub-array 7. ! A total of 82 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 3. ! A total of 40 telescope corrections were flagged in sub-array 1. ! A total of 72 telescope corrections were flagged in sub-array 2. ! A total of 39 telescope corrections were flagged in sub-array 3. ! A total of 39 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! A total of 54 telescope corrections were flagged in sub-array 6. ! A total of 50 telescope corrections were flagged in sub-array 7. ! A total of 82 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 4. ! A total of 40 telescope corrections were flagged in sub-array 1. ! A total of 72 telescope corrections were flagged in sub-array 2. ! A total of 39 telescope corrections were flagged in sub-array 3. ! A total of 39 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! A total of 54 telescope corrections were flagged in sub-array 6. ! A total of 50 telescope corrections were flagged in sub-array 7. ! A total of 82 telescope corrections were flagged in sub-array 8. ! ! Fit before self-cal, rms=1.381687Jy sigma=15.550057 ! Fit after self-cal, rms=0.361428Jy sigma=2.335605 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. print "=== Flag out IVS stations" ! === Flag out IVS stations integer lun string cfg_file lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_flag_process.csh /tmp/ff.ff" ![@/tmp/ff.ff \] flag 2:FT ! Flagging baselines of antenna 2:FT in the currently selected channels. flag 2:MA ! Flagging baselines of antenna 2:MA in the currently selected channels. flag 2:WZ ! Flagging baselines of antenna 2:WZ in the currently selected channels. flag 3:FT ! Flagging baselines of antenna 3:FT in the currently selected channels. flag 3:MA ! Flagging baselines of antenna 3:MA in the currently selected channels. flag 3:WZ ! Flagging baselines of antenna 3:WZ in the currently selected channels. flag 4:FT ! Flagging baselines of antenna 4:FT in the currently selected channels. flag 4:MA ! Flagging baselines of antenna 4:MA in the currently selected channels. flag 5:FT ! Flagging baselines of antenna 5:FT in the currently selected channels. flag 6:FT ! Flagging baselines of antenna 6:FT in the currently selected channels. flag 7:FT ! Flagging baselines of antenna 7:FT in the currently selected channels. flag 8:FT ! Flagging baselines of antenna 8:FT in the currently selected channels. ![Exited script file: /tmp/ff.ff] uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.9345 mas, bmaj=2.733 mas, bpa=0.1503 degrees ! Estimated noise=2.76101 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (0, 0). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.682074 Jy ! Component: 100 - total flux cleaned = 0.884023 Jy ! Total flux subtracted in 100 components = 0.884023 Jy ! Clean residual min=-0.030133 max=0.085816 Jy/beam ! Clean residual mean=0.000110 rms=0.011045 Jy/beam ! Combined flux in latest and established models = 0.884023 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 13 components and 0.884023 Jy ! ! Correcting IF 1. ! A total of 33 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 2. ! A total of 33 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 3. ! A total of 33 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 4. ! A total of 33 telescope corrections were flagged in sub-array 3. ! ! Fit before self-cal, rms=0.252007Jy sigma=3.118557 ! Fit after self-cal, rms=0.252006Jy sigma=3.118497 ! Inverting map and beam ! Estimated beam: bmin=0.9345 mas, bmaj=2.733 mas, bpa=0.1503 degrees ! Estimated noise=2.76101 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0990414 Jy ! Component: 100 - total flux cleaned = 0.15927 Jy ! Total flux subtracted in 100 components = 0.15927 Jy ! Clean residual min=-0.023891 max=0.040627 Jy/beam ! Clean residual mean=0.000005 rms=0.006783 Jy/beam ! Combined flux in latest and established models = 1.04329 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 37 components and 1.04329 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.182581Jy sigma=1.914073 ! Fit after self-cal, rms=0.182552Jy sigma=1.913986 ! Inverting map ! Added new window around map position (-1.4, 1.8). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0407803 Jy ! Component: 100 - total flux cleaned = 0.0674234 Jy ! Total flux subtracted in 100 components = 0.0674234 Jy ! Clean residual min=-0.022313 max=0.033283 Jy/beam ! Clean residual mean=-0.000039 rms=0.005773 Jy/beam ! Combined flux in latest and established models = 1.11072 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 56 components and 1.11072 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.165981Jy sigma=1.546749 ! Fit after self-cal, rms=0.165815Jy sigma=1.544152 ! Inverting map print "=== Finished uniform weighting clean ===" ! === Finished uniform weighting clean === uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.571 mas, bmaj=4.244 mas, bpa=-4.813 degrees ! Estimated noise=1.42263 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.0484676 Jy ! Component: 100 - total flux cleaned = 0.0697952 Jy ! Component: 150 - total flux cleaned = 0.0795395 Jy ! Component: 200 - total flux cleaned = 0.0840241 Jy ! Total flux subtracted in 200 components = 0.0840241 Jy ! Clean residual min=-0.032312 max=0.040904 Jy/beam ! Clean residual mean=-0.000478 rms=0.010565 Jy/beam ! Combined flux in latest and established models = 1.19474 Jy ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 83 components and 1.19474 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.158714Jy sigma=1.320873 ! Fit after self-cal, rms=0.158734Jy sigma=1.319199 ! Inverting map print "=== Finished natural weighting clean for VLBI sites ===" ! === Finished natural weighting clean for VLBI sites === gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.96 FD 0.97 HN 1.43 KP 1.01 ! LA 1.07 MK 0.93 NL 0.98 OV 1.01 ! PT 0.99 SC 0.98 ! ! Telescope amplitude corrections in sub-array 2: ! FT 0.00* MA 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FT 0.00* MA 0.00* SC 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! FT 0.00* MA 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! FT 0.00* SC 0.00* ! ! A total of 21 telescope corrections were flagged in sub-array 6. ! Telescope amplitude corrections in sub-array 6: ! FD 0.98 FT 1.00* HN 1.19 LA 1.25 ! PT 1.17 SC 0.81 ! ! A total of 3 telescope corrections were flagged in sub-array 7. ! Telescope amplitude corrections in sub-array 7: ! FD 0.95 FT 1.00* HN 1.21 KP 1.05 ! LA 1.06 NL 1.22 PT 0.97 SC 0.85 ! ! Telescope amplitude corrections in sub-array 8: ! FD 0.94 FT 1.00* HN 1.17 KP 1.04 ! LA 1.04 NL 0.95 PT 0.94 ! ! Telescope amplitude corrections in sub-array 9: ! FD 0.99 HN 1.30 KP 1.01 LA 1.06 ! NL 0.92 OV 1.02 PT 0.99 SC 0.94 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 FD 0.96 HN 1.41 KP 0.95 ! LA 1.03 MK 0.98 NL 0.96 OV 0.97 ! PT 1.00 SC 1.04 ! ! Telescope amplitude corrections in sub-array 2: ! FT 0.00* MA 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FT 0.00* MA 0.00* SC 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! FT 0.00* MA 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! FT 0.00* SC 0.00* ! ! A total of 21 telescope corrections were flagged in sub-array 6. ! Telescope amplitude corrections in sub-array 6: ! FD 1.01 FT 1.00* HN 1.32 LA 1.19 ! PT 1.19 SC 0.85 ! ! A total of 3 telescope corrections were flagged in sub-array 7. ! Telescope amplitude corrections in sub-array 7: ! FD 0.96 FT 1.00* HN 1.22 KP 0.98 ! LA 1.03 NL 1.23 PT 0.97 SC 0.89 ! ! Telescope amplitude corrections in sub-array 8: ! FD 0.96 FT 1.00* HN 1.26 KP 0.98 ! LA 0.99 NL 0.94 PT 0.95 ! ! Telescope amplitude corrections in sub-array 9: ! FD 1.00 HN 1.28 KP 0.96 LA 1.03 ! NL 0.91 OV 1.00 PT 0.99 SC 1.01 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.03 FD 0.99 HN 1.39 KP 0.86 ! LA 1.09 MK 1.01 NL 0.99 OV 1.02 ! PT 1.00 SC 1.01 ! ! Telescope amplitude corrections in sub-array 2: ! FT 0.00* MA 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FT 0.00* MA 0.00* SC 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! FT 0.00* MA 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! FT 0.00* SC 0.00* ! ! A total of 21 telescope corrections were flagged in sub-array 6. ! Telescope amplitude corrections in sub-array 6: ! FD 0.95 FT 1.00* HN 1.27 LA 1.29 ! PT 1.28 SC 0.83 ! ! A total of 3 telescope corrections were flagged in sub-array 7. ! Telescope amplitude corrections in sub-array 7: ! FD 0.97 FT 1.00* HN 1.28 KP 0.91 ! LA 1.11 NL 1.28 PT 0.96 SC 0.87 ! ! Telescope amplitude corrections in sub-array 8: ! FD 0.95 FT 1.00* HN 1.31 KP 0.90 ! LA 1.07 NL 0.98 PT 0.95 ! ! Telescope amplitude corrections in sub-array 9: ! FD 0.99 HN 1.35 KP 0.87 LA 1.10 ! NL 0.94 OV 1.05 PT 1.01 SC 0.97 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 1.01 FD 1.04 HN 1.51 KP 0.92 ! LA 1.09 MK 1.08 NL 1.04 OV 1.01 ! PT 1.01 SC 1.04 ! ! Telescope amplitude corrections in sub-array 2: ! FT 0.00* MA 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FT 0.00* MA 0.00* SC 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! FT 0.00* MA 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! FT 0.00* SC 0.00* ! ! A total of 21 telescope corrections were flagged in sub-array 6. ! Telescope amplitude corrections in sub-array 6: ! FD 1.04 FT 1.00* HN 1.38 LA 1.28 ! PT 1.26 SC 0.87 ! ! A total of 3 telescope corrections were flagged in sub-array 7. ! Telescope amplitude corrections in sub-array 7: ! FD 1.03 FT 1.00* HN 1.28 KP 0.95 ! LA 1.08 NL 1.32 PT 1.01 SC 0.90 ! ! Telescope amplitude corrections in sub-array 8: ! FD 1.01 FT 1.00* HN 1.22 KP 0.93 ! LA 1.06 NL 1.01 PT 0.98 ! ! Telescope amplitude corrections in sub-array 9: ! FD 1.07 HN 1.42 KP 0.93 LA 1.10 ! NL 0.96 OV 1.03 PT 1.05 SC 1.03 ! ! ! Fit before self-cal, rms=0.158734Jy sigma=1.319199 ! Fit after self-cal, rms=0.122806Jy sigma=0.861153 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.569 mas, bmaj=4.258 mas, bpa=-4.29 degrees ! Estimated noise=1.45183 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.00461989 Jy ! Component: 100 - total flux cleaned = -0.0051249 Jy ! Component: 150 - total flux cleaned = -0.005347 Jy ! Component: 200 - total flux cleaned = -0.00534508 Jy ! Total flux subtracted in 200 components = -0.00534508 Jy ! Clean residual min=-0.016098 max=0.027212 Jy/beam ! Clean residual mean=-0.000013 rms=0.005390 Jy/beam ! Combined flux in latest and established models = 1.1894 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 90 components and 1.1894 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.120042Jy sigma=0.857448 ! Fit after self-cal, rms=0.120387Jy sigma=0.856934 ! Inverting map ! Added new window around map position (3.2, 8.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0251005 Jy ! Component: 100 - total flux cleaned = 0.0365177 Jy ! Component: 150 - total flux cleaned = 0.0401366 Jy ! Component: 200 - total flux cleaned = 0.040634 Jy ! Total flux subtracted in 200 components = 0.040634 Jy ! Clean residual min=-0.012995 max=0.013382 Jy/beam ! Clean residual mean=-0.000009 rms=0.003550 Jy/beam ! Combined flux in latest and established models = 1.23003 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 122 components and 1.23003 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.115467Jy sigma=0.762317 ! Fit after self-cal, rms=0.115190Jy sigma=0.753505 ! 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.115190Jy sigma=0.753505 ! Fit after self-cal, rms=0.115561Jy sigma=0.738595 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.57 mas, bmaj=4.26 mas, bpa=-4.241 degrees ! Estimated noise=1.44884 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00983834 Jy ! Component: 100 - total flux cleaned = 0.0116962 Jy ! Component: 150 - total flux cleaned = 0.0129396 Jy ! Component: 200 - total flux cleaned = 0.0131302 Jy ! Total flux subtracted in 200 components = 0.0131302 Jy ! Clean residual min=-0.009046 max=0.007479 Jy/beam ! Clean residual mean=0.000008 rms=0.002246 Jy/beam ! Combined flux in latest and established models = 1.24316 Jy selfcal ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 130 components and 1.24316 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.112908Jy sigma=0.716886 ! Fit after self-cal, rms=0.112856Jy sigma=0.715042 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.00312963 Jy ! Component: 100 - total flux cleaned = 0.00313898 Jy ! Component: 150 - total flux cleaned = 0.00298337 Jy ! Component: 200 - total flux cleaned = 0.00298899 Jy ! Total flux subtracted in 200 components = 0.00298899 Jy ! Clean residual min=-0.009391 max=0.006918 Jy/beam ! Clean residual mean=0.000008 rms=0.002184 Jy/beam ! Combined flux in latest and established models = 1.24615 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 137 components and 1.24615 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.111555Jy sigma=0.710676 ! Fit after self-cal, rms=0.111564Jy sigma=0.710401 ! 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.111564Jy sigma=0.710401 ! Fit after self-cal, rms=0.112086Jy sigma=0.703960 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.57 mas, bmaj=4.256 mas, bpa=-4.214 degrees ! Estimated noise=1.44724 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00273964 Jy ! Component: 100 - total flux cleaned = 0.00342878 Jy ! Component: 150 - total flux cleaned = 0.00389292 Jy ! Component: 200 - total flux cleaned = 0.00418646 Jy ! Total flux subtracted in 200 components = 0.00418646 Jy ! Clean residual min=-0.006920 max=0.005761 Jy/beam ! Clean residual mean=0.000011 rms=0.001672 Jy/beam ! Combined flux in latest and established models = 1.25033 Jy selfcal ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 142 components and 1.25033 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.110594Jy sigma=0.697643 ! Fit after self-cal, rms=0.110538Jy sigma=0.697164 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.00104284 Jy ! Component: 100 - total flux cleaned = 0.00117692 Jy ! Component: 150 - total flux cleaned = 0.00129601 Jy ! Component: 200 - total flux cleaned = 0.00129369 Jy ! Total flux subtracted in 200 components = 0.00129369 Jy ! Clean residual min=-0.007000 max=0.005839 Jy/beam ! Clean residual mean=0.000010 rms=0.001661 Jy/beam ! Combined flux in latest and established models = 1.25163 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 146 components and 1.25163 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.109675Jy sigma=0.695230 ! Fit after self-cal, rms=0.109668Jy sigma=0.695099 ! 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.109668Jy sigma=0.695099 ! Fit after self-cal, rms=0.110406Jy sigma=0.692510 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.191 mas, bmaj=8.863 mas, bpa=-17.28 degrees ! Estimated noise=1.37693 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00209867 Jy ! Component: 100 - total flux cleaned = 0.00221091 Jy ! Component: 150 - total flux cleaned = 0.00196869 Jy ! Component: 200 - total flux cleaned = 0.00182706 Jy ! Total flux subtracted in 200 components = 0.00182706 Jy ! Clean residual min=-0.006191 max=0.005294 Jy/beam ! Clean residual mean=0.000014 rms=0.001540 Jy/beam ! Combined flux in latest and established models = 1.25346 Jy selfcal ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 156 components and 1.25346 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.110142Jy sigma=0.689944 ! Fit after self-cal, rms=0.110076Jy sigma=0.689615 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.000686823 Jy ! Component: 100 - total flux cleaned = 0.000621628 Jy ! Component: 150 - total flux cleaned = 0.000498372 Jy ! Component: 200 - total flux cleaned = 0.000442143 Jy ! Total flux subtracted in 200 components = 0.000442143 Jy ! Clean residual min=-0.006234 max=0.005245 Jy/beam ! Clean residual mean=0.000014 rms=0.001543 Jy/beam ! Combined flux in latest and established models = 1.2539 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 167 components and 1.2539 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.109917Jy sigma=0.689096 ! Fit after self-cal, rms=0.109914Jy sigma=0.689024 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "=== Finished amplitude self-cal for VLBA stations ===" ! === Finished amplitude self-cal for VLBA stations === unflag * ! Unflagging all baselines in the currently selected channels. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. integer lun lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_selfant_process.csh /tmp/ff.ff" ![@/tmp/ff.ff] selfant 1:BR,true ! 1:BR self-cal status=fixed weight=1 selfant 1:FD,true ! 1:FD self-cal status=fixed weight=1 selfant 1:HN,true ! 1:HN self-cal status=fixed weight=1 selfant 1:KP,true ! 1:KP self-cal status=fixed weight=1 selfant 1:LA,true ! 1:LA self-cal status=fixed weight=1 selfant 1:MK,true ! 1:MK self-cal status=fixed weight=1 selfant 1:NL,true ! 1:NL self-cal status=fixed weight=1 selfant 1:OV,true ! 1:OV self-cal status=fixed weight=1 selfant 1:PT,true ! 1:PT self-cal status=fixed weight=1 selfant 1:SC,true ! 1:SC self-cal status=fixed weight=1 selfant 3:SC,true ! 3:SC self-cal status=fixed weight=1 selfant 4:SC,true ! 4:SC self-cal status=fixed weight=1 selfant 5:SC,true ! 5:SC 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:LA,true ! 6:LA self-cal status=fixed weight=1 selfant 6:PT,true ! 6:PT self-cal status=fixed weight=1 selfant 6:SC,true ! 6:SC self-cal status=fixed weight=1 selfant 7:FD,true ! 7:FD self-cal status=fixed weight=1 selfant 7:HN,true ! 7:HN self-cal status=fixed weight=1 selfant 7:KP,true ! 7:KP self-cal status=fixed weight=1 selfant 7:LA,true ! 7:LA self-cal status=fixed weight=1 selfant 7:NL,true ! 7:NL self-cal status=fixed weight=1 selfant 7:PT,true ! 7:PT self-cal status=fixed weight=1 selfant 7:SC,true ! 7:SC self-cal status=fixed weight=1 selfant 8:FD,true ! 8:FD self-cal status=fixed weight=1 selfant 8:HN,true ! 8:HN self-cal status=fixed weight=1 selfant 8:KP,true ! 8:KP self-cal status=fixed weight=1 selfant 8:LA,true ! 8:LA self-cal status=fixed weight=1 selfant 8:NL,true ! 8:NL self-cal status=fixed weight=1 selfant 8:PT,true ! 8:PT self-cal status=fixed weight=1 selfant 9:FD,true ! 9:FD self-cal status=fixed weight=1 selfant 9:HN,true ! 9:HN self-cal status=fixed weight=1 selfant 9:KP,true ! 9:KP self-cal status=fixed weight=1 selfant 9:LA,true ! 9:LA self-cal status=fixed weight=1 selfant 9:NL,true ! 9:NL self-cal status=fixed weight=1 selfant 9:OV,true ! 9:OV self-cal status=fixed weight=1 selfant 9:PT,true ! 9:PT self-cal status=fixed weight=1 selfant 9:SC,true ! 9: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 72 telescope corrections were flagged in sub-array 2. ! A total of 54 telescope corrections were flagged in sub-array 3. ! A total of 26 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! A total of 17 telescope corrections were flagged in sub-array 7. ! A total of 26 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 2. ! A total of 72 telescope corrections were flagged in sub-array 2. ! A total of 54 telescope corrections were flagged in sub-array 3. ! A total of 26 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! A total of 17 telescope corrections were flagged in sub-array 7. ! A total of 26 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 3. ! A total of 72 telescope corrections were flagged in sub-array 2. ! A total of 54 telescope corrections were flagged in sub-array 3. ! A total of 26 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! A total of 17 telescope corrections were flagged in sub-array 7. ! A total of 26 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 4. ! A total of 72 telescope corrections were flagged in sub-array 2. ! A total of 54 telescope corrections were flagged in sub-array 3. ! A total of 26 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! A total of 17 telescope corrections were flagged in sub-array 7. ! A total of 26 telescope corrections were flagged in sub-array 8. ! ! Fit before self-cal, rms=0.369736Jy sigma=1.871383 ! Fit after self-cal, rms=0.162605Jy sigma=1.366271 selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:KP self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:MK self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:OV self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 1:SC self-cal status=correctable weight=1 ! 2:FT self-cal status=correctable weight=1 ! 2:MA self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 3:FT self-cal status=correctable weight=1 ! 3:MA self-cal status=correctable weight=1 ! 3:SC self-cal status=correctable weight=1 ! 3:WZ self-cal status=correctable weight=1 ! 4:FT self-cal status=correctable weight=1 ! 4:MA self-cal status=correctable weight=1 ! 4:SC self-cal status=correctable weight=1 ! 5:FT self-cal status=correctable weight=1 ! 5:SC self-cal status=correctable weight=1 ! 6:FD self-cal status=correctable weight=1 ! 6:FT self-cal status=correctable weight=1 ! 6:HN self-cal status=correctable weight=1 ! 6:LA self-cal status=correctable weight=1 ! 6:PT self-cal status=correctable weight=1 ! 6:SC self-cal status=correctable weight=1 ! 7:FD self-cal status=correctable weight=1 ! 7:FT self-cal status=correctable weight=1 ! 7:HN self-cal status=correctable weight=1 ! 7:KP self-cal status=correctable weight=1 ! 7:LA self-cal status=correctable weight=1 ! 7:NL self-cal status=correctable weight=1 ! 7:PT self-cal status=correctable weight=1 ! 7:SC self-cal status=correctable weight=1 ! 8:FD self-cal status=correctable weight=1 ! 8:FT self-cal status=correctable weight=1 ! 8:HN self-cal status=correctable weight=1 ! 8:KP self-cal status=correctable weight=1 ! 8:LA self-cal status=correctable weight=1 ! 8:NL self-cal status=correctable weight=1 ! 8:PT self-cal status=correctable weight=1 ! 9:FD self-cal status=correctable weight=1 ! 9:HN self-cal status=correctable weight=1 ! 9:KP self-cal status=correctable weight=1 ! 9:LA self-cal status=correctable weight=1 ! 9:NL self-cal status=correctable weight=1 ! 9:OV self-cal status=correctable weight=1 ! 9:PT self-cal status=correctable weight=1 ! 9: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 40 telescope corrections were flagged in sub-array 1. ! A total of 18 telescope corrections were flagged in sub-array 3. ! A total of 13 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 58 telescope corrections were flagged in sub-array 6. ! A total of 36 telescope corrections were flagged in sub-array 7. ! A total of 56 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 2. ! A total of 40 telescope corrections were flagged in sub-array 1. ! A total of 18 telescope corrections were flagged in sub-array 3. ! A total of 13 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 58 telescope corrections were flagged in sub-array 6. ! A total of 36 telescope corrections were flagged in sub-array 7. ! A total of 56 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 3. ! A total of 40 telescope corrections were flagged in sub-array 1. ! A total of 18 telescope corrections were flagged in sub-array 3. ! A total of 13 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 58 telescope corrections were flagged in sub-array 6. ! A total of 36 telescope corrections were flagged in sub-array 7. ! A total of 56 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 4. ! A total of 40 telescope corrections were flagged in sub-array 1. ! A total of 18 telescope corrections were flagged in sub-array 3. ! A total of 13 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 58 telescope corrections were flagged in sub-array 6. ! A total of 36 telescope corrections were flagged in sub-array 7. ! A total of 56 telescope corrections were flagged in sub-array 8. ! ! Fit before self-cal, rms=0.162605Jy sigma=1.366271 ! Fit after self-cal, rms=0.109922Jy sigma=0.687866 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.9325 mas, bmaj=2.748 mas, bpa=0.6059 degrees ! Estimated noise=2.8201 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.00669791 Jy ! Component: 100 - total flux cleaned = -0.00593222 Jy ! Component: 150 - total flux cleaned = -0.00479041 Jy ! Component: 200 - total flux cleaned = -0.00396096 Jy ! Total flux subtracted in 200 components = -0.00396096 Jy ! Clean residual min=-0.007751 max=0.007625 Jy/beam ! Clean residual mean=0.000001 rms=0.001937 Jy/beam ! Combined flux in latest and established models = 1.24994 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 186 components and 1.24994 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.108464Jy sigma=0.686982 ! Fit after self-cal, rms=0.108378Jy sigma=0.686793 ! 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=1.573 mas, bmaj=4.265 mas, bpa=-4.179 degrees ! Estimated noise=1.4459 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.00385654 Jy ! Component: 100 - total flux cleaned = 0.004659 Jy ! Component: 150 - total flux cleaned = 0.00465801 Jy ! Component: 200 - total flux cleaned = 0.00465896 Jy ! Component: 250 - total flux cleaned = 0.00458761 Jy ! Component: 300 - total flux cleaned = 0.00459053 Jy ! Component: 350 - total flux cleaned = 0.00459006 Jy ! Component: 400 - total flux cleaned = 0.00458819 Jy ! Total flux subtracted in 400 components = 0.00458819 Jy ! Clean residual min=-0.005709 max=0.005198 Jy/beam ! Clean residual mean=0.000011 rms=0.001307 Jy/beam ! Combined flux in latest and established models = 1.25452 Jy ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 199 components and 1.25452 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.107590Jy sigma=0.684009 ! Fit after self-cal, rms=0.107583Jy sigma=0.683931 ! 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.000122837 Jy ! Component: 100 - total flux cleaned = -6.21365e-06 Jy ! Component: 150 - total flux cleaned = -7.01362e-05 Jy ! Component: 200 - total flux cleaned = -0.00013124 Jy ! Component: 250 - total flux cleaned = -0.000188052 Jy ! Component: 300 - total flux cleaned = -0.000187146 Jy ! Component: 350 - total flux cleaned = -0.000128295 Jy ! Component: 400 - total flux cleaned = -0.000184135 Jy ! Total flux subtracted in 400 components = -0.000184135 Jy ! Clean residual min=-0.005735 max=0.005212 Jy/beam ! Clean residual mean=0.000011 rms=0.001299 Jy/beam ! Combined flux in latest and established models = 1.25434 Jy ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 214 components and 1.25434 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.107113Jy sigma=0.683235 ! Fit after self-cal, rms=0.107106Jy sigma=0.683193 ! Inverting map print "=== Finished natural weighting clean for al sites ===" ! === Finished natural weighting clean for al sites === unflag * ! Unflagging all baselines in the currently selected channels. ! Applying 2560 buffered edits. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 40 telescope corrections were flagged in sub-array 1. ! A total of 72 telescope corrections were flagged in sub-array 2. ! A total of 72 telescope corrections were flagged in sub-array 3. ! A total of 39 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! A total of 75 telescope corrections were flagged in sub-array 6. ! A total of 53 telescope corrections were flagged in sub-array 7. ! A total of 82 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 2. ! A total of 40 telescope corrections were flagged in sub-array 1. ! A total of 72 telescope corrections were flagged in sub-array 2. ! A total of 72 telescope corrections were flagged in sub-array 3. ! A total of 39 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! A total of 75 telescope corrections were flagged in sub-array 6. ! A total of 53 telescope corrections were flagged in sub-array 7. ! A total of 82 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 3. ! A total of 40 telescope corrections were flagged in sub-array 1. ! A total of 72 telescope corrections were flagged in sub-array 2. ! A total of 72 telescope corrections were flagged in sub-array 3. ! A total of 39 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! A total of 75 telescope corrections were flagged in sub-array 6. ! A total of 53 telescope corrections were flagged in sub-array 7. ! A total of 82 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 4. ! A total of 40 telescope corrections were flagged in sub-array 1. ! A total of 72 telescope corrections were flagged in sub-array 2. ! A total of 72 telescope corrections were flagged in sub-array 3. ! A total of 39 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! A total of 75 telescope corrections were flagged in sub-array 6. ! A total of 53 telescope corrections were flagged in sub-array 7. ! A total of 82 telescope corrections were flagged in sub-array 8. ! ! Fit before self-cal, rms=0.369776Jy sigma=1.872012 ! Fit after self-cal, rms=0.108158Jy sigma=0.682180 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.576 mas, bmaj=4.274 mas, bpa=-4.071 degrees ! Estimated noise=1.44673 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000370479 Jy ! Component: 100 - total flux cleaned = -0.000371714 Jy ! Component: 150 - total flux cleaned = -0.000217672 Jy ! Component: 200 - total flux cleaned = -7.60833e-05 Jy ! Component: 250 - total flux cleaned = -8.39009e-06 Jy ! Component: 300 - total flux cleaned = -6.75441e-06 Jy ! Component: 350 - total flux cleaned = -6.66279e-06 Jy ! Component: 400 - total flux cleaned = -6.50582e-06 Jy ! Total flux subtracted in 400 components = -6.50582e-06 Jy ! Clean residual min=-0.005426 max=0.004981 Jy/beam ! Clean residual mean=0.000011 rms=0.001233 Jy/beam ! Combined flux in latest and established models = 1.25433 Jy selfcal ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 216 components and 1.25433 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.107471Jy sigma=0.681107 ! Fit after self-cal, rms=0.107446Jy sigma=0.681042 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 = 6.12209e-05 Jy ! Component: 100 - total flux cleaned = 0.000121941 Jy ! Component: 150 - total flux cleaned = 0.000123346 Jy ! Component: 200 - total flux cleaned = 0.000180195 Jy ! Component: 250 - total flux cleaned = 0.000181073 Jy ! Component: 300 - total flux cleaned = 0.000181727 Jy ! Component: 350 - total flux cleaned = 0.000182961 Jy ! Component: 400 - total flux cleaned = 0.000180767 Jy ! Total flux subtracted in 400 components = 0.000180767 Jy ! Clean residual min=-0.005403 max=0.004982 Jy/beam ! Clean residual mean=0.000011 rms=0.001229 Jy/beam ! Combined flux in latest and established models = 1.25451 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 218 components and 1.25451 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.107072Jy sigma=0.680338 ! Fit after self-cal, rms=0.107069Jy sigma=0.680303 ! 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.107069Jy sigma=0.680303 ! Fit after self-cal, rms=0.107709Jy sigma=0.679710 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.578 mas, bmaj=4.281 mas, bpa=-3.992 degrees ! Estimated noise=1.44726 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000441508 Jy ! Component: 100 - total flux cleaned = -0.000441494 Jy ! Component: 150 - total flux cleaned = -0.000379942 Jy ! Component: 200 - total flux cleaned = -0.000380575 Jy ! Component: 250 - total flux cleaned = -0.000324029 Jy ! Component: 300 - total flux cleaned = -0.000267989 Jy ! Component: 350 - total flux cleaned = -0.000322205 Jy ! Component: 400 - total flux cleaned = -0.000269288 Jy ! Total flux subtracted in 400 components = -0.000269288 Jy ! Clean residual min=-0.005152 max=0.004799 Jy/beam ! Clean residual mean=0.000012 rms=0.001192 Jy/beam ! Combined flux in latest and established models = 1.25425 Jy selfcal ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 221 components and 1.25425 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.107246Jy sigma=0.678892 ! Fit after self-cal, rms=0.107225Jy sigma=0.678843 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 = -1.05012e-05 Jy ! Component: 100 - total flux cleaned = -1.01527e-05 Jy ! Component: 150 - total flux cleaned = -1.02289e-05 Jy ! Component: 200 - total flux cleaned = -1.2293e-05 Jy ! Component: 250 - total flux cleaned = -5.92033e-05 Jy ! Component: 300 - total flux cleaned = -5.80573e-05 Jy ! Component: 350 - total flux cleaned = -5.86151e-05 Jy ! Component: 400 - total flux cleaned = -1.35525e-05 Jy ! Total flux subtracted in 400 components = -1.35525e-05 Jy ! Clean residual min=-0.005135 max=0.004777 Jy/beam ! Clean residual mean=0.000012 rms=0.001186 Jy/beam ! Combined flux in latest and established models = 1.25423 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 224 components and 1.25423 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.106949Jy sigma=0.678266 ! Fit after self-cal, rms=0.106947Jy sigma=0.678240 ! 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.106947Jy sigma=0.678240 ! Fit after self-cal, rms=0.106400Jy sigma=0.657793 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.207 mas, bmaj=8.913 mas, bpa=-17.28 degrees ! Estimated noise=1.37394 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000351619 Jy ! Component: 100 - total flux cleaned = 0.000277633 Jy ! Component: 150 - total flux cleaned = 0.000276493 Jy ! Component: 200 - total flux cleaned = 0.000334607 Jy ! Component: 250 - total flux cleaned = 0.000391051 Jy ! Component: 300 - total flux cleaned = 0.000416853 Jy ! Component: 350 - total flux cleaned = 0.000470524 Jy ! Component: 400 - total flux cleaned = 0.000496984 Jy ! Total flux subtracted in 400 components = 0.000496984 Jy ! Clean residual min=-0.004970 max=0.004412 Jy/beam ! Clean residual mean=0.000029 rms=0.001250 Jy/beam ! Combined flux in latest and established models = 1.25473 Jy selfcal ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 234 components and 1.25473 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.106348Jy sigma=0.657421 ! Fit after self-cal, rms=0.106327Jy sigma=0.657373 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.000141219 Jy ! Component: 100 - total flux cleaned = 0.000141748 Jy ! Component: 150 - total flux cleaned = 0.00016604 Jy ! Component: 200 - total flux cleaned = 0.000189628 Jy ! Component: 250 - total flux cleaned = 0.000213198 Jy ! Component: 300 - total flux cleaned = 0.000235915 Jy ! Component: 350 - total flux cleaned = 0.000258226 Jy ! Component: 400 - total flux cleaned = 0.000258119 Jy ! Total flux subtracted in 400 components = 0.000258119 Jy ! Clean residual min=-0.005002 max=0.004351 Jy/beam ! Clean residual mean=0.000029 rms=0.001248 Jy/beam ! Combined flux in latest and established models = 1.25499 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 241 components and 1.25499 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.106298Jy sigma=0.657174 ! Fit after self-cal, rms=0.106276Jy sigma=0.657161 ! 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=1.579 mas, bmaj=4.286 mas, bpa=-3.917 degrees ! Estimated noise=1.44767 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.579 x 4.286 at -3.917 degrees (North through East) ! Clean map min=-0.0066852 max=0.95563 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.106276Jy sigma=0.657161 ! Fit after self-cal, rms=0.102273Jy sigma=0.567470 ! Inverting map and beam ! Estimated beam: bmin=1.581 mas, bmaj=4.291 mas, bpa=-3.948 degrees ! Estimated noise=1.44745 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000724389 Jy ! Component: 100 - total flux cleaned = -0.000724957 Jy ! Component: 150 - total flux cleaned = -0.000724243 Jy ! Component: 200 - total flux cleaned = -0.000723133 Jy ! Component: 250 - total flux cleaned = -0.000722793 Jy ! Component: 300 - total flux cleaned = -0.000724153 Jy ! Component: 350 - total flux cleaned = -0.000767569 Jy ! Component: 400 - total flux cleaned = -0.000767003 Jy ! Total flux subtracted in 400 components = -0.000767003 Jy ! Clean residual min=-0.004913 max=0.004456 Jy/beam ! Clean residual mean=0.000013 rms=0.001123 Jy/beam ! Combined flux in latest and established models = 1.25422 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 250 components and 1.25422 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.101960Jy sigma=0.566659 ! Fit after self-cal, rms=0.101935Jy sigma=0.566605 ! Inverting map ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -5.66099e-05 Jy ! Component: 100 - total flux cleaned = -5.89646e-05 Jy ! Component: 150 - total flux cleaned = -5.98886e-05 Jy ! Component: 200 - total flux cleaned = -9.90246e-05 Jy ! Component: 250 - total flux cleaned = -0.000138105 Jy ! Component: 300 - total flux cleaned = -0.000137767 Jy ! Component: 350 - total flux cleaned = -0.000136598 Jy ! Component: 400 - total flux cleaned = -0.00017399 Jy ! Total flux subtracted in 400 components = -0.00017399 Jy ! Clean residual min=-0.004865 max=0.004443 Jy/beam ! Clean residual mean=0.000013 rms=0.001116 Jy/beam ! Combined flux in latest and established models = 1.25405 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 254 components and 1.25405 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.101735Jy sigma=0.566053 ! Fit after self-cal, rms=0.101724Jy sigma=0.566029 ! 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.101724Jy sigma=0.566029 ! Fit after self-cal, rms=0.102219Jy sigma=0.565810 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.9337 mas, bmaj=2.769 mas, bpa=0.8989 degrees ! Estimated noise=2.8408 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.698505 Jy ! Component: 100 - total flux cleaned = 0.898506 Jy ! Component: 150 - total flux cleaned = 0.992894 Jy ! Component: 200 - total flux cleaned = 1.04789 Jy ! Component: 250 - total flux cleaned = 1.08734 Jy ! Component: 300 - total flux cleaned = 1.11683 Jy ! Component: 350 - total flux cleaned = 1.13997 Jy ! Component: 400 - total flux cleaned = 1.15883 Jy ! Total flux subtracted in 400 components = 1.15883 Jy ! Clean residual min=-0.010592 max=0.013591 Jy/beam ! Clean residual mean=0.000040 rms=0.002603 Jy/beam ! Combined flux in latest and established models = 1.15883 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 64 model components to the UV plane model. ! The established model now contains 64 components and 1.15883 Jy ! Inverting map and beam ! Estimated beam: bmin=1.582 mas, bmaj=4.296 mas, bpa=-3.889 degrees ! Estimated noise=1.44771 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, 1.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0308171 Jy ! Component: 100 - total flux cleaned = 0.0504604 Jy ! Component: 150 - total flux cleaned = 0.0640743 Jy ! Component: 200 - total flux cleaned = 0.0736984 Jy ! Component: 250 - total flux cleaned = 0.080689 Jy ! Component: 300 - total flux cleaned = 0.0859152 Jy ! Component: 350 - total flux cleaned = 0.088859 Jy ! Component: 400 - total flux cleaned = 0.0902437 Jy ! Total flux subtracted in 400 components = 0.0902437 Jy ! Clean residual min=-0.005382 max=0.004075 Jy/beam ! Clean residual mean=0.000016 rms=0.001247 Jy/beam ! Combined flux in latest and established models = 1.24908 Jy ! Adding 122 model components to the UV plane model. ! The established model now contains 173 components and 1.24908 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=3.211 mas, bmaj=8.931 mas, bpa=-17.29 degrees ! Estimated noise=1.37285 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00294026 Jy ! Component: 100 - total flux cleaned = 0.00382262 Jy ! Component: 150 - total flux cleaned = 0.00398062 Jy ! Component: 200 - total flux cleaned = 0.0041256 Jy ! Component: 250 - total flux cleaned = 0.00425993 Jy ! Component: 300 - total flux cleaned = 0.00432235 Jy ! Component: 350 - total flux cleaned = 0.00438366 Jy ! Component: 400 - total flux cleaned = 0.00444245 Jy ! Total flux subtracted in 400 components = 0.00444245 Jy ! Clean residual min=-0.005033 max=0.003870 Jy/beam ! Clean residual mean=0.000021 rms=0.001239 Jy/beam ! Combined flux in latest and established models = 1.25352 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 = 5.40532e-05 Jy ! Component: 100 - total flux cleaned = 0.000107045 Jy ! Component: 150 - total flux cleaned = 0.000159864 Jy ! Component: 200 - total flux cleaned = 0.000207769 Jy ! Component: 250 - total flux cleaned = 0.000255367 Jy ! Component: 300 - total flux cleaned = 0.000303284 Jy ! Component: 350 - total flux cleaned = 0.000305286 Jy ! Component: 400 - total flux cleaned = 0.000351351 Jy ! Total flux subtracted in 400 components = 0.000351351 Jy ! Clean residual min=-0.004975 max=0.003879 Jy/beam ! Clean residual mean=0.000021 rms=0.001230 Jy/beam ! Combined flux in latest and established models = 1.25387 Jy ! Adding 66 model components to the UV plane model. ! The established model now contains 234 components and 1.25387 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=1.582 mas, bmaj=4.296 mas, bpa=-3.889 degrees ! Estimated noise=1.44771 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000532757 Jy ! Component: 100 - total flux cleaned = -0.0006074 Jy ! Component: 150 - total flux cleaned = -0.000537859 Jy ! Component: 200 - total flux cleaned = -0.000534891 Jy ! Component: 250 - total flux cleaned = -0.000469075 Jy ! Component: 300 - total flux cleaned = -0.00047098 Jy ! Component: 350 - total flux cleaned = -0.000469987 Jy ! Component: 400 - total flux cleaned = -0.000409273 Jy ! Total flux subtracted in 400 components = -0.000409273 Jy ! Clean residual min=-0.004743 max=0.004013 Jy/beam ! Clean residual mean=0.000016 rms=0.001107 Jy/beam ! Combined flux in latest and established models = 1.25346 Jy selfcal ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 261 components and 1.25346 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.103003Jy sigma=0.562354 ! Fit after self-cal, rms=0.102922Jy sigma=0.561803 wmodel /scr/pima/rdv95_uvs/J1257-3155_X_map.mod ! Writing 261 model components to file: /scr/pima/rdv95_uvs/J1257-3155_X_map.mod wobs /scr/pima/rdv95_uvs/J1257-3155_X_uvs.fits ! Applying 1398 buffered edits. ! Writing UV FITS file: /scr/pima/rdv95_uvs/J1257-3155_X_uvs.fits wwins /scr/pima/rdv95_uvs/J1257-3155_X_map.win ! wwins: Wrote 4 windows to /scr/pima/rdv95_uvs/J1257-3155_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.00143206 Jy ! Component: 100 - total flux cleaned = -0.001925 Jy ! Component: 150 - total flux cleaned = -0.00243863 Jy ! Component: 200 - total flux cleaned = -0.0026795 Jy ! Component: 250 - total flux cleaned = -0.00268098 Jy ! Component: 300 - total flux cleaned = -0.00268106 Jy ! Component: 350 - total flux cleaned = -0.00257521 Jy ! Component: 400 - total flux cleaned = -0.00237182 Jy ! Component: 450 - total flux cleaned = -0.00217612 Jy ! Component: 500 - total flux cleaned = -0.00198446 Jy ! Component: 550 - total flux cleaned = -0.00189206 Jy ! Component: 600 - total flux cleaned = -0.00189219 Jy ! Component: 650 - total flux cleaned = -0.00184768 Jy ! Component: 700 - total flux cleaned = -0.00184756 Jy ! Component: 750 - total flux cleaned = -0.00184691 Jy ! Component: 800 - total flux cleaned = -0.00176119 Jy ! Component: 850 - total flux cleaned = -0.00180291 Jy ! Component: 900 - total flux cleaned = -0.00180187 Jy ! Component: 950 - total flux cleaned = -0.00163849 Jy ! Component: 1000 - total flux cleaned = -0.00163843 Jy ! Component: 1050 - total flux cleaned = -0.00155889 Jy ! Component: 1100 - total flux cleaned = -0.00151979 Jy ! Component: 1150 - total flux cleaned = -0.00148114 Jy ! Component: 1200 - total flux cleaned = -0.00144261 Jy ! Component: 1250 - total flux cleaned = -0.00136716 Jy ! Component: 1300 - total flux cleaned = -0.00140466 Jy ! Component: 1350 - total flux cleaned = -0.00129294 Jy ! Component: 1400 - total flux cleaned = -0.00121961 Jy ! Component: 1450 - total flux cleaned = -0.00121983 Jy ! Component: 1500 - total flux cleaned = -0.00114754 Jy ! Component: 1550 - total flux cleaned = -0.00111212 Jy ! Component: 1600 - total flux cleaned = -0.00097125 Jy ! Component: 1650 - total flux cleaned = -0.00107638 Jy ! Component: 1700 - total flux cleaned = -0.00100688 Jy ! Component: 1750 - total flux cleaned = -0.00104126 Jy ! Component: 1800 - total flux cleaned = -0.000972932 Jy ! Component: 1850 - total flux cleaned = -0.000938852 Jy ! Component: 1900 - total flux cleaned = -0.000837932 Jy ! Component: 1950 - total flux cleaned = -0.000771153 Jy ! Component: 2000 - total flux cleaned = -0.000770997 Jy ! Component: 2050 - total flux cleaned = -0.000738056 Jy ! Component: 2100 - total flux cleaned = -0.000639829 Jy ! Component: 2150 - total flux cleaned = -0.000510148 Jy ! Component: 2200 - total flux cleaned = -0.000478042 Jy ! Component: 2250 - total flux cleaned = -0.000477844 Jy ! Component: 2300 - total flux cleaned = -0.000382323 Jy ! Component: 2350 - total flux cleaned = -0.000382328 Jy ! Component: 2400 - total flux cleaned = -0.00031937 Jy ! Component: 2450 - total flux cleaned = -0.000350718 Jy ! Component: 2500 - total flux cleaned = -0.000257438 Jy ! Component: 2550 - total flux cleaned = -0.000319407 Jy ! Component: 2600 - total flux cleaned = -0.000350045 Jy ! Component: 2650 - total flux cleaned = -0.000197474 Jy ! Component: 2700 - total flux cleaned = -0.000167335 Jy ! Component: 2750 - total flux cleaned = -4.6763e-05 Jy ! Component: 2800 - total flux cleaned = -0.000136792 Jy ! Component: 2850 - total flux cleaned = 1.2368e-05 Jy ! Component: 2900 - total flux cleaned = 0.000101364 Jy ! Component: 2950 - total flux cleaned = 7.18065e-05 Jy ! Component: 3000 - total flux cleaned = 0.00018926 Jy ! Component: 3050 - total flux cleaned = 0.000218167 Jy ! Component: 3100 - total flux cleaned = 0.00021809 Jy ! Component: 3150 - total flux cleaned = 0.000276099 Jy ! Component: 3200 - total flux cleaned = 0.000304936 Jy ! Component: 3250 - total flux cleaned = 0.000419833 Jy ! Component: 3300 - total flux cleaned = 0.00044857 Jy ! Component: 3350 - total flux cleaned = 0.000505532 Jy ! Component: 3400 - total flux cleaned = 0.000505521 Jy ! Component: 3450 - total flux cleaned = 0.00050525 Jy ! Component: 3500 - total flux cleaned = 0.00050525 Jy ! Component: 3550 - total flux cleaned = 0.000617662 Jy ! Component: 3600 - total flux cleaned = 0.000589654 Jy ! Component: 3650 - total flux cleaned = 0.000561662 Jy ! Component: 3700 - total flux cleaned = 0.000617304 Jy ! Component: 3750 - total flux cleaned = 0.000617013 Jy ! Component: 3800 - total flux cleaned = 0.000700145 Jy ! Component: 3850 - total flux cleaned = 0.000644668 Jy ! Component: 3900 - total flux cleaned = 0.00064483 Jy ! Component: 3950 - total flux cleaned = 0.000644825 Jy ! Component: 4000 - total flux cleaned = 0.000672286 Jy ! Component: 4050 - total flux cleaned = 0.00058992 Jy ! Total flux subtracted in 4096 components = 0.000671916 Jy ! Clean residual min=-0.001604 max=0.001560 Jy/beam ! Clean residual mean=0.000018 rms=0.000488 Jy/beam ! Combined flux in latest and established models = 1.25413 Jy keep ! Adding 225 model components to the UV plane model. ! The established model now contains 486 components and 1.25413 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv95_uvs/J1257-3155_X\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.582 x 4.296 at -3.889 degrees (North through East) ! Clean map min=-0.008052 max=0.9699 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv95_uvs/J1257-3155_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= 0.9698 rms= 0.00048 Jy quit ! Quitting program ! Log file /scr/pima/rdv95_uvs/J1257-3155_X_uvs.log closed on Wed Nov 6 20:46:57 2013