! Started logfile: /scr/pima/rdv62_uvs/J0927+3902_S_uvs.log on Fri Nov 22 22:35:11 2013 obs /scr/pima/rdv62_uvs/J0927+3902_S_uva.fits ! Reading UV FITS file: /scr/pima/rdv62_uvs/J0927+3902_S_uva.fits ! AN table 1: 21 integrations on 15 of 15 possible baselines. ! AN table 2: 5 integrations on 15 of 15 possible baselines. ! AN table 3: 8 integrations on 36 of 36 possible baselines. ! AN table 4: 5 integrations on 6 of 6 possible baselines. ! AN table 5: 6 integrations on 55 of 55 possible baselines. ! AN table 6: 3 integrations on 1 of 1 possible baselines. ! Apparent sampling: 0.591739 visibilities/baseline/integration-bin. ! Found source: J0927+3902 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 2.23299e+09 8e+06 1 8e+06 ! 02 2 2.26299e+09 8e+06 1 8e+06 ! 03 3 2.35299e+09 8e+06 1 8e+06 ! 04 4 2.37299e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 2464 visibilities. select RR ! Selecting polarization: RR, channels: 1..4 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 print "field_size =",field_size ! field_size = 1024 print "field_cell =",field_cell ! field_cell = 0.4 print "taper_size =",taper_size ! taper_size = 20 ![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv62_uvs/J0927+3902_S] integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 800000 float soltime2 soltime2 = 15 float soltime3 soltime3 = 1 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float old_peak float new_peak float flux_cutoff selfant "",false ! 1:HH self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:NY self-cal status=correctable weight=1 ! 1:SC self-cal status=correctable weight=1 ! 1:ZC self-cal status=correctable weight=1 ! 2:BR 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 ! 3:BR 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:MK self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 3:OV 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:NL self-cal status=correctable weight=1 ! 4:NY 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:NY self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 6:NY self-cal status=correctable weight=1 ! 6:WZ self-cal status=correctable weight=1 system "rm /tmp/ff.ff" mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.400x0.400 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 69 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 8 telescope corrections were flagged in sub-array 4. ! A total of 18 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 69 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 8 telescope corrections were flagged in sub-array 4. ! A total of 18 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 69 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 8 telescope corrections were flagged in sub-array 4. ! A total of 18 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 69 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 8 telescope corrections were flagged in sub-array 4. ! A total of 18 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=4.695375Jy sigma=228.909123 ! Fit after self-cal, rms=3.696993Jy sigma=183.813743 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. print "=== Flag out IVS stations" ! === Flag out IVS stations integer lun string cfg_file lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_flag_process.csh /tmp/ff.ff" ![@/tmp/ff.ff \] flag 1:HH ! Flagging baselines of antenna 1:HH in the currently selected channels. flag 1:NY ! Flagging baselines of antenna 1:NY in the currently selected channels. flag 1:ZC ! Flagging baselines of antenna 1:ZC in the currently selected channels. flag 4:NY ! Flagging baselines of antenna 4:NY in the currently selected channels. flag 5:KK ! Flagging baselines of antenna 5:KK in the currently selected channels. flag 5:NY ! Flagging baselines of antenna 5:NY in the currently selected channels. flag 6:NY ! Flagging baselines of antenna 6:NY in the currently selected channels. flag 6:WZ ! Flagging baselines of antenna 6:WZ in the currently selected channels. ![Exited script file: /tmp/ff.ff] uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.706 mas, bmaj=6.114 mas, bpa=-42.31 degrees ! Estimated noise=0.531767 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.76948 Jy ! Component: 100 - total flux cleaned = 4.64824 Jy ! Total flux subtracted in 100 components = 4.64824 Jy ! Clean residual min=-0.121151 max=0.328841 Jy/beam ! Clean residual mean=0.000583 rms=0.045964 Jy/beam ! Combined flux in latest and established models = 4.64824 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 11 components and 4.64824 Jy ! ! Correcting IF 1. ! A total of 30 telescope corrections were flagged in sub-array 1. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 30 telescope corrections were flagged in sub-array 1. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 30 telescope corrections were flagged in sub-array 1. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 30 telescope corrections were flagged in sub-array 1. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.963248Jy sigma=49.683498 ! Fit after self-cal, rms=0.965800Jy sigma=49.868533 ! Inverting map and beam ! Estimated beam: bmin=3.706 mas, bmaj=6.114 mas, bpa=-42.31 degrees ! Estimated noise=0.531831 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.337776 Jy ! Component: 100 - total flux cleaned = 0.51147 Jy ! Total flux subtracted in 100 components = 0.51147 Jy ! Clean residual min=-0.115505 max=0.116178 Jy/beam ! Clean residual mean=0.000285 rms=0.032274 Jy/beam ! Combined flux in latest and established models = 5.15971 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 25 components and 5.15971 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.835212Jy sigma=41.645649 ! Fit after self-cal, rms=0.835133Jy sigma=41.641528 ! Inverting map print "=== Finished uniform weighting clean ===" ! === Finished uniform weighting clean === uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=5.098 mas, bmaj=7.728 mas, bpa=-30.92 degrees ! Estimated noise=0.35075 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.134608 Jy ! Component: 100 - total flux cleaned = 0.210223 Jy ! Component: 150 - total flux cleaned = 0.254709 Jy ! Component: 200 - total flux cleaned = 0.2813 Jy ! Total flux subtracted in 200 components = 0.2813 Jy ! Clean residual min=-0.114252 max=0.110029 Jy/beam ! Clean residual mean=-0.000210 rms=0.033449 Jy/beam ! Combined flux in latest and established models = 5.44101 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 47 components and 5.44101 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.820241Jy sigma=40.296225 ! Fit after self-cal, rms=0.820220Jy sigma=40.295347 ! 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. ! A total of 27 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* HN 0.00* NL 0.00* NY 0.00* ! SC 0.00* ZC 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.03 FD 1.08 HN 0.93 LA 1.19 ! NL 1.04 PT 1.12 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.03 FD 0.98 HN 0.98 KP 0.97 ! LA 1.15 MK 0.92 NL 1.01 OV 0.99 ! PT 1.06 ! ! A total of 9 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! FD 0.00* HN 0.00* NL 0.00* NY 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.05 FD 1.01 HN 0.94 KK 1.00* ! KP 1.07 LA 1.20 MK 0.94 NL 1.01 ! NY 1.00* OV 0.98 PT 1.07 ! ! Telescope amplitude corrections in sub-array 6: ! NY 0.00* WZ 0.00* ! ! ! Correcting IF 2. ! A total of 27 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* HN 0.00* NL 0.00* NY 0.00* ! SC 0.00* ZC 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.05 FD 1.14 HN 1.06 LA 1.11 ! NL 1.07 PT 1.11 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.05 FD 1.05 HN 1.11 KP 1.02 ! LA 1.04 MK 1.06 NL 1.03 OV 1.03 ! PT 1.05 ! ! A total of 9 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! FD 0.00* HN 0.00* NL 0.00* NY 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.06 FD 1.07 HN 1.06 KK 1.00* ! KP 1.09 LA 1.09 MK 0.96 NL 1.26 ! NY 1.00* OV 1.04 PT 1.05 ! ! Telescope amplitude corrections in sub-array 6: ! NY 0.00* WZ 0.00* ! ! ! Correcting IF 3. ! A total of 27 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* HN 0.00* NL 0.00* NY 0.00* ! SC 0.00* ZC 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.08 FD 1.13 HN 0.97 LA 0.85 ! NL 1.11 PT 1.10 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.09 FD 1.03 HN 1.02 KP 0.99 ! LA 0.82 MK 1.05 NL 1.06 OV 1.04 ! PT 1.03 ! ! A total of 9 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! FD 0.00* HN 0.00* NL 0.00* NY 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.11 FD 1.06 HN 0.97 KK 1.00* ! KP 1.07 LA 0.85 MK 0.96 NL 1.06 ! NY 1.00* OV 1.04 PT 1.05 ! ! Telescope amplitude corrections in sub-array 6: ! NY 0.00* WZ 0.00* ! ! ! Correcting IF 4. ! A total of 27 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* HN 0.00* NL 0.00* NY 0.00* ! SC 0.00* ZC 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.96 FD 1.09 HN 0.92 LA 0.74 ! NL 1.00 PT 1.06 ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.99 FD 1.00 HN 0.96 KP 0.92 ! LA 0.71 MK 1.02 NL 0.96 OV 0.87 ! PT 0.97 ! ! A total of 9 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! FD 0.00* HN 0.00* NL 0.00* NY 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.03 FD 1.02 HN 0.92 KK 1.00* ! KP 0.99 LA 0.69 MK 0.95 NL 0.96 ! NY 1.00* OV 0.87 PT 0.98 ! ! Telescope amplitude corrections in sub-array 6: ! NY 0.00* WZ 0.00* ! ! ! Fit before self-cal, rms=0.820220Jy sigma=40.295347 ! Fit after self-cal, rms=0.167951Jy sigma=6.530927 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=5.099 mas, bmaj=7.791 mas, bpa=-27.22 degrees ! Estimated noise=0.355505 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.0132062 Jy ! Component: 100 - total flux cleaned = -0.0149956 Jy ! Component: 150 - total flux cleaned = -0.0165345 Jy ! Component: 200 - total flux cleaned = -0.0178718 Jy ! Total flux subtracted in 200 components = -0.0178718 Jy ! Clean residual min=-0.014432 max=0.018039 Jy/beam ! Clean residual mean=-0.000015 rms=0.003871 Jy/beam ! Combined flux in latest and established models = 5.42314 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 59 components and 5.42314 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.165496Jy sigma=6.407762 ! Fit after self-cal, rms=0.165524Jy sigma=6.405301 ! 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.165524Jy sigma=6.405301 ! Fit after self-cal, rms=0.164578Jy sigma=6.375229 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=5.102 mas, bmaj=7.801 mas, bpa=-27.15 degrees ! Estimated noise=0.355376 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00999441 Jy ! Component: 100 - total flux cleaned = -0.0121642 Jy ! Component: 150 - total flux cleaned = -0.013423 Jy ! Component: 200 - total flux cleaned = -0.014263 Jy ! Total flux subtracted in 200 components = -0.014263 Jy ! Clean residual min=-0.013183 max=0.016387 Jy/beam ! Clean residual mean=-0.000026 rms=0.003550 Jy/beam ! Combined flux in latest and established models = 5.40888 Jy selfcal ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 62 components and 5.40888 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.162979Jy sigma=6.292972 ! Fit after self-cal, rms=0.162962Jy sigma=6.292413 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (18.4, -7.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0114169 Jy ! Component: 100 - total flux cleaned = 0.0120871 Jy ! Component: 150 - total flux cleaned = 0.0114936 Jy ! Component: 200 - total flux cleaned = 0.0112141 Jy ! Total flux subtracted in 200 components = 0.0112141 Jy ! Clean residual min=-0.012559 max=0.014343 Jy/beam ! Clean residual mean=-0.000091 rms=0.003287 Jy/beam ! Combined flux in latest and established models = 5.42009 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 74 components and 5.42009 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.161005Jy sigma=6.200564 ! Fit after self-cal, rms=0.160958Jy sigma=6.197346 ! 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.160958Jy sigma=6.197346 ! Fit after self-cal, rms=0.160489Jy sigma=6.182136 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=5.104 mas, bmaj=7.808 mas, bpa=-27.09 degrees ! Estimated noise=0.355279 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00721211 Jy ! Component: 100 - total flux cleaned = -0.00749953 Jy ! Component: 150 - total flux cleaned = -0.00775238 Jy ! Component: 200 - total flux cleaned = -0.00822758 Jy ! Total flux subtracted in 200 components = -0.00822758 Jy ! Clean residual min=-0.012368 max=0.013048 Jy/beam ! Clean residual mean=-0.000113 rms=0.003073 Jy/beam ! Combined flux in latest and established models = 5.41186 Jy selfcal ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 81 components and 5.41186 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.159448Jy sigma=6.134358 ! Fit after self-cal, rms=0.159438Jy sigma=6.134100 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.000202369 Jy ! Component: 100 - total flux cleaned = -0.000615471 Jy ! Component: 150 - total flux cleaned = -0.000811567 Jy ! Component: 200 - total flux cleaned = -0.00118791 Jy ! Total flux subtracted in 200 components = -0.00118791 Jy ! Clean residual min=-0.012138 max=0.012755 Jy/beam ! Clean residual mean=-0.000105 rms=0.003019 Jy/beam ! Combined flux in latest and established models = 5.41068 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 88 components and 5.41068 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.158913Jy sigma=6.109983 ! Fit after self-cal, rms=0.158899Jy sigma=6.109720 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime2 ! Performing amp+phase self-cal over 15 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.158899Jy sigma=6.109720 ! Fit after self-cal, rms=0.158582Jy sigma=6.097413 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=10.12 mas, bmaj=10.34 mas, bpa=-53.39 degrees ! Estimated noise=0.387122 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00642728 Jy ! Component: 100 - total flux cleaned = -0.00772556 Jy ! Component: 150 - total flux cleaned = -0.00851803 Jy ! Component: 200 - total flux cleaned = -0.00908614 Jy ! Total flux subtracted in 200 components = -0.00908614 Jy ! Clean residual min=-0.011437 max=0.011873 Jy/beam ! Clean residual mean=-0.000181 rms=0.003077 Jy/beam ! Combined flux in latest and established models = 5.40159 Jy selfcal ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 104 components and 5.40159 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.158152Jy sigma=6.077281 ! Fit after self-cal, rms=0.158126Jy sigma=6.075567 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.000263125 Jy ! Component: 100 - total flux cleaned = -0.000998338 Jy ! Component: 150 - total flux cleaned = -0.00156327 Jy ! Component: 200 - total flux cleaned = -0.00173843 Jy ! Total flux subtracted in 200 components = -0.00173843 Jy ! Clean residual min=-0.011431 max=0.012135 Jy/beam ! Clean residual mean=-0.000162 rms=0.003068 Jy/beam ! Combined flux in latest and established models = 5.39985 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 111 components and 5.39985 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.157994Jy sigma=6.068804 ! Fit after self-cal, rms=0.157983Jy sigma=6.067930 ! 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:HN,true ! 1:HN self-cal status=fixed weight=1 selfant 1:NL,true ! 1:NL self-cal status=fixed weight=1 selfant 1:SC,true ! 1:SC self-cal status=fixed weight=1 selfant 2:BR,true ! 2:BR self-cal status=fixed weight=1 selfant 2:FD,true ! 2:FD self-cal status=fixed weight=1 selfant 2: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 3:BR,true ! 3:BR 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:MK,true ! 3:MK self-cal status=fixed weight=1 selfant 3:NL,true ! 3:NL self-cal status=fixed weight=1 selfant 3:OV,true ! 3:OV 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:NL,true ! 4:NL 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 ![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 44 telescope corrections were flagged in sub-array 1. ! A total of 2 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 44 telescope corrections were flagged in sub-array 1. ! A total of 2 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 44 telescope corrections were flagged in sub-array 1. ! A total of 2 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 44 telescope corrections were flagged in sub-array 1. ! A total of 2 telescope corrections were flagged in sub-array 4. ! A total of 6 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=1.521445Jy sigma=30.594340 ! Fit after self-cal, rms=0.991609Jy sigma=23.739763 selfant "",false ! 1:HH self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:NY self-cal status=correctable weight=1 ! 1:SC self-cal status=correctable weight=1 ! 1:ZC self-cal status=correctable weight=1 ! 2:BR 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 ! 3:BR 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:MK self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 3:OV 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:NL self-cal status=correctable weight=1 ! 4:NY 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:NY self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 6:NY self-cal status=correctable weight=1 ! 6:WZ 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 34 telescope corrections were flagged in sub-array 1. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 12 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 34 telescope corrections were flagged in sub-array 1. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 12 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 34 telescope corrections were flagged in sub-array 1. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 12 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 34 telescope corrections were flagged in sub-array 1. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 12 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.991609Jy sigma=23.739763 ! Fit after self-cal, rms=0.754492Jy sigma=6.909593 uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.411 mas, bmaj=4.275 mas, bpa=-50.68 degrees ! Estimated noise=0.612419 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.00138978 Jy ! Component: 100 - total flux cleaned = -0.00639525 Jy ! Component: 150 - total flux cleaned = -0.0142716 Jy ! Component: 200 - total flux cleaned = -0.0183529 Jy ! Total flux subtracted in 200 components = -0.0183529 Jy ! Clean residual min=-0.026479 max=0.028130 Jy/beam ! Clean residual mean=-0.000041 rms=0.006246 Jy/beam ! Combined flux in latest and established models = 5.3815 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 123 components and 5.3815 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.736908Jy sigma=6.630690 ! Fit after self-cal, rms=0.717854Jy sigma=6.609368 ! Inverting map ! Added new window around map position (-15.2, 10.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00866627 Jy ! Component: 100 - total flux cleaned = 0.00874049 Jy ! Component: 150 - total flux cleaned = 0.0077279 Jy ! Component: 200 - total flux cleaned = 0.00673766 Jy ! Total flux subtracted in 200 components = 0.00673766 Jy ! Clean residual min=-0.023655 max=0.024076 Jy/beam ! Clean residual mean=-0.000064 rms=0.005819 Jy/beam ! Combined flux in latest and established models = 5.38823 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 138 components and 5.38823 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.708962Jy sigma=6.563647 ! Fit after self-cal, rms=0.708908Jy sigma=6.557484 ! 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=4.612 mas, bmaj=5.582 mas, bpa=-28.98 degrees ! Estimated noise=0.36775 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.0087757 Jy ! Component: 100 - total flux cleaned = 0.0112333 Jy ! Component: 150 - total flux cleaned = 0.0134627 Jy ! Component: 200 - total flux cleaned = 0.0134628 Jy ! Component: 250 - total flux cleaned = 0.0138693 Jy ! Component: 300 - total flux cleaned = 0.013865 Jy ! Component: 350 - total flux cleaned = 0.0142742 Jy ! Component: 400 - total flux cleaned = 0.0138771 Jy ! Total flux subtracted in 400 components = 0.0138771 Jy ! Clean residual min=-0.015017 max=0.014686 Jy/beam ! Clean residual mean=-0.000123 rms=0.003727 Jy/beam ! Combined flux in latest and established models = 5.40211 Jy ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 161 components and 5.40211 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.703643Jy sigma=6.497016 ! Fit after self-cal, rms=0.703607Jy sigma=6.494637 ! 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.000448161 Jy ! Component: 100 - total flux cleaned = 0.00084612 Jy ! Component: 150 - total flux cleaned = 0.000471718 Jy ! Component: 200 - total flux cleaned = 0.000847876 Jy ! Component: 250 - total flux cleaned = 0.000841823 Jy ! Component: 300 - total flux cleaned = 0.000476743 Jy ! Component: 350 - total flux cleaned = 0.000842124 Jy ! Component: 400 - total flux cleaned = 0.000494687 Jy ! Total flux subtracted in 400 components = 0.000494687 Jy ! Clean residual min=-0.014410 max=0.014341 Jy/beam ! Clean residual mean=-0.000124 rms=0.003681 Jy/beam ! Combined flux in latest and established models = 5.40261 Jy ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 173 components and 5.40261 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.699051Jy sigma=6.469881 ! Fit after self-cal, rms=0.699034Jy sigma=6.468555 ! 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 78 telescope corrections were flagged in sub-array 1. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 8 telescope corrections were flagged in sub-array 4. ! A total of 18 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 78 telescope corrections were flagged in sub-array 1. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 8 telescope corrections were flagged in sub-array 4. ! A total of 18 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 78 telescope corrections were flagged in sub-array 1. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 8 telescope corrections were flagged in sub-array 4. ! A total of 18 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 78 telescope corrections were flagged in sub-array 1. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 8 telescope corrections were flagged in sub-array 4. ! A total of 18 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=1.233735Jy sigma=24.012387 ! Fit after self-cal, rms=0.686548Jy sigma=6.278378 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=4.614 mas, bmaj=5.556 mas, bpa=-27.73 degrees ! Estimated noise=0.367768 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00140349 Jy ! Component: 100 - total flux cleaned = -0.00271948 Jy ! Component: 150 - total flux cleaned = -0.00391515 Jy ! Component: 200 - total flux cleaned = -0.0050194 Jy ! Component: 250 - total flux cleaned = -0.00607825 Jy ! Component: 300 - total flux cleaned = -0.00706656 Jy ! Component: 350 - total flux cleaned = -0.00801894 Jy ! Component: 400 - total flux cleaned = -0.00890745 Jy ! Total flux subtracted in 400 components = -0.00890745 Jy ! Clean residual min=-0.012633 max=0.011825 Jy/beam ! Clean residual mean=-0.000168 rms=0.002991 Jy/beam ! Combined flux in latest and established models = 5.3937 Jy selfcal ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 179 components and 5.3937 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.675930Jy sigma=6.173285 ! Fit after self-cal, rms=0.675887Jy sigma=6.171819 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.000842381 Jy ! Component: 100 - total flux cleaned = -0.0016747 Jy ! Component: 150 - total flux cleaned = -0.0024905 Jy ! Component: 200 - total flux cleaned = -0.00326576 Jy ! Component: 250 - total flux cleaned = -0.0044198 Jy ! Component: 300 - total flux cleaned = -0.00479519 Jy ! Component: 350 - total flux cleaned = -0.00515092 Jy ! Component: 400 - total flux cleaned = -0.00586158 Jy ! Total flux subtracted in 400 components = -0.00586158 Jy ! Clean residual min=-0.012748 max=0.011486 Jy/beam ! Clean residual mean=-0.000140 rms=0.002972 Jy/beam ! Combined flux in latest and established models = 5.38784 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 183 components and 5.38784 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.669246Jy sigma=6.153450 ! Fit after self-cal, rms=0.669229Jy sigma=6.152779 ! 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.669229Jy sigma=6.152779 ! Fit after self-cal, rms=0.662208Jy sigma=6.082017 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=4.615 mas, bmaj=5.541 mas, bpa=-26.85 degrees ! Estimated noise=0.367835 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000480133 Jy ! Component: 100 - total flux cleaned = -0.00157242 Jy ! Component: 150 - total flux cleaned = -0.00255643 Jy ! Component: 200 - total flux cleaned = -0.00348823 Jy ! Component: 250 - total flux cleaned = -0.00439209 Jy ! Component: 300 - total flux cleaned = -0.00563204 Jy ! Component: 350 - total flux cleaned = -0.00645855 Jy ! Component: 400 - total flux cleaned = -0.00724176 Jy ! Total flux subtracted in 400 components = -0.00724176 Jy ! Clean residual min=-0.011480 max=0.010537 Jy/beam ! Clean residual mean=-0.000188 rms=0.002658 Jy/beam ! Combined flux in latest and established models = 5.3806 Jy selfcal ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 186 components and 5.3806 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.653533Jy sigma=6.024066 ! Fit after self-cal, rms=0.653528Jy sigma=6.023189 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.000795115 Jy ! Component: 100 - total flux cleaned = -0.00189186 Jy ! Component: 150 - total flux cleaned = -0.00224927 Jy ! Component: 200 - total flux cleaned = -0.00260248 Jy ! Component: 250 - total flux cleaned = -0.00327603 Jy ! Component: 300 - total flux cleaned = -0.00392649 Jy ! Component: 350 - total flux cleaned = -0.00393125 Jy ! Component: 400 - total flux cleaned = -0.0045537 Jy ! Total flux subtracted in 400 components = -0.0045537 Jy ! Clean residual min=-0.011656 max=0.010184 Jy/beam ! Clean residual mean=-0.000166 rms=0.002634 Jy/beam ! Combined flux in latest and established models = 5.37604 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 191 components and 5.37604 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.648157Jy sigma=6.009246 ! Fit after self-cal, rms=0.648164Jy sigma=6.008719 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime3 ! Performing amp+phase self-cal over 1 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.648164Jy sigma=6.008719 ! Fit after self-cal, rms=0.635123Jy sigma=5.362539 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=9.616 mas, bmaj=10.5 mas, bpa=-62.75 degrees ! Estimated noise=0.378313 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00503385 Jy ! Component: 100 - total flux cleaned = -0.00832697 Jy ! Component: 150 - total flux cleaned = -0.0104406 Jy ! Component: 200 - total flux cleaned = -0.0111221 Jy ! Component: 250 - total flux cleaned = -0.0112771 Jy ! Component: 300 - total flux cleaned = -0.0114223 Jy ! Component: 350 - total flux cleaned = -0.011487 Jy ! Component: 400 - total flux cleaned = -0.0115502 Jy ! Total flux subtracted in 400 components = -0.0115502 Jy ! Clean residual min=-0.008483 max=0.006592 Jy/beam ! Clean residual mean=-0.000288 rms=0.002303 Jy/beam ! Combined flux in latest and established models = 5.36449 Jy selfcal ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 229 components and 5.36449 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.634947Jy sigma=5.346635 ! Fit after self-cal, rms=0.634975Jy sigma=5.345451 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.00100602 Jy ! Component: 100 - total flux cleaned = -0.00108178 Jy ! Component: 150 - total flux cleaned = -0.00114514 Jy ! Component: 200 - total flux cleaned = -0.00114558 Jy ! Component: 250 - total flux cleaned = -0.0011466 Jy ! Component: 300 - total flux cleaned = -0.00120047 Jy ! Component: 350 - total flux cleaned = -0.0012003 Jy ! Component: 400 - total flux cleaned = -0.0012032 Jy ! Total flux subtracted in 400 components = -0.0012032 Jy ! Clean residual min=-0.008281 max=0.006484 Jy/beam ! Clean residual mean=-0.000279 rms=0.002282 Jy/beam ! Combined flux in latest and established models = 5.36329 Jy ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 248 components and 5.36329 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.634902Jy sigma=5.338074 ! Fit after self-cal, rms=0.634909Jy sigma=5.337658 ! 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=4.613 mas, bmaj=5.528 mas, bpa=-26.2 degrees ! Estimated noise=0.367787 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 4.613 x 5.528 at -26.2 degrees (North through East) ! Clean map min=-0.012237 max=4.8818 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.634909Jy sigma=5.337658 ! Fit after self-cal, rms=3.853110Jy sigma=3.740584 ! Inverting map and beam ! Estimated beam: bmin=4.616 mas, bmaj=5.564 mas, bpa=-26.55 degrees ! Estimated noise=0.367981 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00150721 Jy ! Component: 100 - total flux cleaned = -0.00251002 Jy ! Component: 150 - total flux cleaned = -0.00343443 Jy ! Component: 200 - total flux cleaned = -0.00402354 Jy ! Component: 250 - total flux cleaned = -0.00429293 Jy ! Component: 300 - total flux cleaned = -0.00456948 Jy ! Component: 350 - total flux cleaned = -0.00483103 Jy ! Component: 400 - total flux cleaned = -0.00534827 Jy ! Total flux subtracted in 400 components = -0.00534827 Jy ! Clean residual min=-0.009500 max=0.009397 Jy/beam ! Clean residual mean=-0.000205 rms=0.002173 Jy/beam ! Combined flux in latest and established models = 5.35794 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 260 components and 5.35794 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=3.852955Jy sigma=3.695090 ! Fit after self-cal, rms=3.852990Jy sigma=3.694471 ! Inverting map ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -1.21072e-08 Jy ! Component: 100 - total flux cleaned = 6.27578e-06 Jy ! Component: 150 - total flux cleaned = -0.000479817 Jy ! Component: 200 - total flux cleaned = -0.000721243 Jy ! Component: 250 - total flux cleaned = -0.000958449 Jy ! Component: 300 - total flux cleaned = -0.00119308 Jy ! Component: 350 - total flux cleaned = -0.00118917 Jy ! Component: 400 - total flux cleaned = -0.00140783 Jy ! Total flux subtracted in 400 components = -0.00140783 Jy ! Clean residual min=-0.009361 max=0.009428 Jy/beam ! Clean residual mean=-0.000200 rms=0.002155 Jy/beam ! Combined flux in latest and established models = 5.35653 Jy ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 275 components and 5.35653 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=3.852931Jy sigma=3.673948 ! Fit after self-cal, rms=3.852973Jy sigma=3.673366 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=3.852973Jy sigma=3.673366 ! Fit after self-cal, rms=3.857183Jy sigma=3.652956 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=3.435 mas, bmaj=4.217 mas, bpa=-49.24 degrees ! Estimated noise=0.617139 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 3.75105 Jy ! Component: 100 - total flux cleaned = 4.61577 Jy ! Component: 150 - total flux cleaned = 4.94254 Jy ! Component: 200 - total flux cleaned = 5.10149 Jy ! Component: 250 - total flux cleaned = 5.18691 Jy ! Component: 300 - total flux cleaned = 5.23831 Jy ! Component: 350 - total flux cleaned = 5.27285 Jy ! Component: 400 - total flux cleaned = 5.28787 Jy ! Total flux subtracted in 400 components = 5.28787 Jy ! Clean residual min=-0.018567 max=0.024996 Jy/beam ! Clean residual mean=-0.000029 rms=0.004308 Jy/beam ! Combined flux in latest and established models = 5.28787 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 35 model components to the UV plane model. ! The established model now contains 35 components and 5.28787 Jy ! Inverting map and beam ! Estimated beam: bmin=4.617 mas, bmaj=5.561 mas, bpa=-26.29 degrees ! Estimated noise=0.368093 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 (6.4, -0.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0253971 Jy ! Component: 100 - total flux cleaned = 0.0416986 Jy ! Component: 150 - total flux cleaned = 0.0504076 Jy ! Component: 200 - total flux cleaned = 0.0547113 Jy ! Component: 250 - total flux cleaned = 0.0550921 Jy ! Component: 300 - total flux cleaned = 0.055834 Jy ! Component: 350 - total flux cleaned = 0.0561891 Jy ! Component: 400 - total flux cleaned = 0.0565298 Jy ! Total flux subtracted in 400 components = 0.0565298 Jy ! Clean residual min=-0.008934 max=0.009656 Jy/beam ! Clean residual mean=-0.000177 rms=0.002335 Jy/beam ! Combined flux in latest and established models = 5.3444 Jy ! Adding 41 model components to the UV plane model. ! The established model now contains 73 components and 5.3444 Jy ! Inverting map ! Added new window around map position (-7.2, 0.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00467191 Jy ! Component: 100 - total flux cleaned = 0.00363993 Jy ! Component: 150 - total flux cleaned = 0.00266242 Jy ! Component: 200 - total flux cleaned = 0.00267175 Jy ! Component: 250 - total flux cleaned = 0.00237366 Jy ! Component: 300 - total flux cleaned = 0.00297141 Jy ! Component: 350 - total flux cleaned = 0.00268043 Jy ! Component: 400 - total flux cleaned = 0.00296105 Jy ! Total flux subtracted in 400 components = 0.00296105 Jy ! Clean residual min=-0.008212 max=0.008460 Jy/beam ! Clean residual mean=-0.000188 rms=0.002111 Jy/beam ! Combined flux in latest and established models = 5.34736 Jy ! Adding 51 model components to the UV plane model. ! The established model now contains 108 components and 5.34737 Jy ! Inverting map ! Added new window around map position (95.2, -24.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0028858 Jy ! Component: 100 - total flux cleaned = 0.00203145 Jy ! Component: 150 - total flux cleaned = 0.00285352 Jy ! Component: 200 - total flux cleaned = 0.00284206 Jy ! Component: 250 - total flux cleaned = 0.00310408 Jy ! Component: 300 - total flux cleaned = 0.00284675 Jy ! Component: 350 - total flux cleaned = 0.00335867 Jy ! Component: 400 - total flux cleaned = 0.00311025 Jy ! Total flux subtracted in 400 components = 0.00311025 Jy ! Clean residual min=-0.007926 max=0.007540 Jy/beam ! Clean residual mean=-0.000202 rms=0.001972 Jy/beam ! Combined flux in latest and established models = 5.35048 Jy ! Adding 58 model components to the UV plane model. ! The established model now contains 131 components and 5.35048 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=9.624 mas, bmaj=10.51 mas, bpa=-62.49 degrees ! Estimated noise=0.378074 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000549636 Jy ! Component: 100 - total flux cleaned = 0.001107 Jy ! Component: 150 - total flux cleaned = 0.00111175 Jy ! Component: 200 - total flux cleaned = 0.00103309 Jy ! Component: 250 - total flux cleaned = 0.00111409 Jy ! Component: 300 - total flux cleaned = 0.00118965 Jy ! Component: 350 - total flux cleaned = 0.00133143 Jy ! Component: 400 - total flux cleaned = 0.00146569 Jy ! Total flux subtracted in 400 components = 0.00146569 Jy ! Clean residual min=-0.007483 max=0.006060 Jy/beam ! Clean residual mean=-0.000349 rms=0.001924 Jy/beam ! Combined flux in latest and established models = 5.35194 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 = -7.68072e-07 Jy ! Component: 100 - total flux cleaned = 6.16303e-05 Jy ! Component: 150 - total flux cleaned = 0.00012073 Jy ! Component: 200 - total flux cleaned = 0.000176693 Jy ! Component: 250 - total flux cleaned = 6.84196e-05 Jy ! Component: 300 - total flux cleaned = 1.62e-05 Jy ! Component: 350 - total flux cleaned = 1.5631e-05 Jy ! Component: 400 - total flux cleaned = -3.34673e-05 Jy ! Total flux subtracted in 400 components = -3.34673e-05 Jy ! Clean residual min=-0.007319 max=0.005903 Jy/beam ! Clean residual mean=-0.000347 rms=0.001894 Jy/beam ! Combined flux in latest and established models = 5.35191 Jy ! Adding 77 model components to the UV plane model. ! The established model now contains 207 components and 5.35191 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=4.617 mas, bmaj=5.561 mas, bpa=-26.29 degrees ! Estimated noise=0.368093 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00344681 Jy ! Component: 100 - total flux cleaned = -0.0041721 Jy ! Component: 150 - total flux cleaned = -0.00394172 Jy ! Component: 200 - total flux cleaned = -0.00371265 Jy ! Component: 250 - total flux cleaned = -0.00348888 Jy ! Component: 300 - total flux cleaned = -0.00326772 Jy ! Component: 350 - total flux cleaned = -0.0030476 Jy ! Component: 400 - total flux cleaned = -0.00305369 Jy ! Total flux subtracted in 400 components = -0.00305369 Jy ! Clean residual min=-0.007152 max=0.006608 Jy/beam ! Clean residual mean=-0.000186 rms=0.001843 Jy/beam ! Combined flux in latest and established models = 5.34886 Jy selfcal ! Performing phase self-cal ! Adding 58 model components to the UV plane model. ! The established model now contains 231 components and 5.34886 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=3.858078Jy sigma=3.501996 ! Fit after self-cal, rms=3.858045Jy sigma=3.496519 wmodel /scr/pima/rdv62_uvs/J0927+3902_S_map.mod ! Writing 231 model components to file: /scr/pima/rdv62_uvs/J0927+3902_S_map.mod wobs /scr/pima/rdv62_uvs/J0927+3902_S_uvs.fits ! Applying 2168 buffered edits. ! Writing UV FITS file: /scr/pima/rdv62_uvs/J0927+3902_S_uvs.fits wwins /scr/pima/rdv62_uvs/J0927+3902_S_map.win ! wwins: Wrote 6 windows to /scr/pima/rdv62_uvs/J0927+3902_S_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.00331417 Jy ! Component: 100 - total flux cleaned = -0.0054972 Jy ! Component: 150 - total flux cleaned = -0.00734694 Jy ! Component: 200 - total flux cleaned = -0.00866642 Jy ! Component: 250 - total flux cleaned = -0.00995036 Jy ! Component: 300 - total flux cleaned = -0.0107471 Jy ! Component: 350 - total flux cleaned = -0.0117526 Jy ! Component: 400 - total flux cleaned = -0.0125191 Jy ! Component: 450 - total flux cleaned = -0.0133821 Jy ! Component: 500 - total flux cleaned = -0.0140183 Jy ! Component: 550 - total flux cleaned = -0.0148557 Jy ! Component: 600 - total flux cleaned = -0.0157799 Jy ! Component: 650 - total flux cleaned = -0.0165888 Jy ! Component: 700 - total flux cleaned = -0.0171871 Jy ! Component: 750 - total flux cleaned = -0.0179744 Jy ! Component: 800 - total flux cleaned = -0.0186553 Jy ! Component: 850 - total flux cleaned = -0.0192313 Jy ! Component: 900 - total flux cleaned = -0.0198019 Jy ! Component: 950 - total flux cleaned = -0.0202718 Jy ! Component: 1000 - total flux cleaned = -0.0210165 Jy ! Component: 1050 - total flux cleaned = -0.0215676 Jy ! Component: 1100 - total flux cleaned = -0.0221129 Jy ! Component: 1150 - total flux cleaned = -0.0226527 Jy ! Component: 1200 - total flux cleaned = -0.023009 Jy ! Component: 1250 - total flux cleaned = -0.0236287 Jy ! Component: 1300 - total flux cleaned = -0.0242426 Jy ! Component: 1350 - total flux cleaned = -0.0249369 Jy ! Component: 1400 - total flux cleaned = -0.0254535 Jy ! Component: 1450 - total flux cleaned = -0.025795 Jy ! Component: 1500 - total flux cleaned = -0.0265575 Jy ! Component: 1550 - total flux cleaned = -0.0271462 Jy ! Component: 1600 - total flux cleaned = -0.0276475 Jy ! Component: 1650 - total flux cleaned = -0.0282281 Jy ! Component: 1700 - total flux cleaned = -0.028805 Jy ! Component: 1750 - total flux cleaned = -0.0296225 Jy ! Component: 1800 - total flux cleaned = -0.0301095 Jy ! Component: 1850 - total flux cleaned = -0.0306741 Jy ! Component: 1900 - total flux cleaned = -0.031236 Jy ! Component: 1950 - total flux cleaned = -0.0317143 Jy ! Component: 2000 - total flux cleaned = -0.0323477 Jy ! Component: 2050 - total flux cleaned = -0.032821 Jy ! Component: 2100 - total flux cleaned = -0.0335262 Jy ! Component: 2150 - total flux cleaned = -0.0339158 Jy ! Component: 2200 - total flux cleaned = -0.0346129 Jy ! Component: 2250 - total flux cleaned = -0.0349205 Jy ! Component: 2300 - total flux cleaned = -0.0353038 Jy ! Component: 2350 - total flux cleaned = -0.0357612 Jy ! Component: 2400 - total flux cleaned = -0.0360643 Jy ! Component: 2450 - total flux cleaned = -0.0364413 Jy ! Component: 2500 - total flux cleaned = -0.0367418 Jy ! Component: 2550 - total flux cleaned = -0.0368913 Jy ! Component: 2600 - total flux cleaned = -0.0373374 Jy ! Component: 2650 - total flux cleaned = -0.0376335 Jy ! Component: 2700 - total flux cleaned = -0.0380012 Jy ! Component: 2750 - total flux cleaned = -0.0380745 Jy ! Component: 2800 - total flux cleaned = -0.0382934 Jy ! Component: 2850 - total flux cleaned = -0.0385841 Jy ! Component: 2900 - total flux cleaned = -0.0386565 Jy ! Component: 2950 - total flux cleaned = -0.0388732 Jy ! Component: 3000 - total flux cleaned = -0.0388735 Jy ! Component: 3050 - total flux cleaned = -0.0393739 Jy ! Component: 3100 - total flux cleaned = -0.0395162 Jy ! Component: 3150 - total flux cleaned = -0.0395868 Jy ! Component: 3200 - total flux cleaned = -0.0398693 Jy ! Component: 3250 - total flux cleaned = -0.0400096 Jy ! Component: 3300 - total flux cleaned = -0.0400096 Jy ! Component: 3350 - total flux cleaned = -0.0402888 Jy ! Component: 3400 - total flux cleaned = -0.0404278 Jy ! Component: 3450 - total flux cleaned = -0.0404974 Jy ! Component: 3500 - total flux cleaned = -0.0406357 Jy ! Component: 3550 - total flux cleaned = -0.0407734 Jy ! Component: 3600 - total flux cleaned = -0.0409103 Jy ! Component: 3650 - total flux cleaned = -0.0411153 Jy ! Component: 3700 - total flux cleaned = -0.0412513 Jy ! Component: 3750 - total flux cleaned = -0.0413192 Jy ! Component: 3800 - total flux cleaned = -0.0415892 Jy ! Component: 3850 - total flux cleaned = -0.041589 Jy ! Component: 3900 - total flux cleaned = -0.0416558 Jy ! Component: 3950 - total flux cleaned = -0.0418564 Jy ! Component: 4000 - total flux cleaned = -0.0418564 Jy ! Component: 4050 - total flux cleaned = -0.0421218 Jy ! Total flux subtracted in 4096 components = -0.0422541 Jy ! Clean residual min=-0.003686 max=0.003552 Jy/beam ! Clean residual mean=-0.000006 rms=0.001152 Jy/beam ! Combined flux in latest and established models = 5.3066 Jy keep ! Adding 243 model components to the UV plane model. ! The established model now contains 465 components and 5.30661 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv62_uvs/J0927+3902_S\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 4.617 x 5.561 at -26.29 degrees (North through East) ! Clean map min=-0.012111 max=4.9003 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv62_uvs/J0927+3902_S_map.fits float peak_flux peak_flux = peak(flux,max) invert ! Inverting map float image_peak image_peak = int(10000*peak_flux)/10000 float image_rms image_rms = int(100000*imstat(rms))/100000 print "=========================================" ! ========================================= print "The image peak= ", image_peak, " rms= ", image_rms, " Jy" ! The image peak= 4.9002 rms= 0.00115 Jy quit ! Quitting program ! Log file /scr/pima/rdv62_uvs/J0927+3902_S_uvs.log closed on Fri Nov 22 22:35:17 2013