! Started logfile: /d2/scr/pima/bp177b_uvs/J0019+7327_X_uvs.log on Sat Mar 22 08:06:26 2014 obs /d2/scr/pima/bp177b_uvs/J0019+7327_X_uva.fits ! Reading UV FITS file: /d2/scr/pima/bp177b_uvs/J0019+7327_X_uva.fits ! AN table 1: 3 integrations on 28 of 28 possible baselines. ! AN table 2: 3 integrations on 21 of 21 possible baselines. ! AN table 3: 8 integrations on 28 of 28 possible baselines. ! Apparent sampling: 1 visibilities/baseline/integration-bin. ! Found source: J0019+7327 ! ! There are 8 IFs, and a total of 8 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 7.392e+09 3.2e+07 1 3.2e+07 ! 02 2 7.424e+09 3.2e+07 1 3.2e+07 ! 03 3 7.456e+09 3.2e+07 1 3.2e+07 ! 04 4 7.552e+09 3.2e+07 1 3.2e+07 ! 05 5 7.744e+09 3.2e+07 1 3.2e+07 ! 06 6 7.776e+09 3.2e+07 1 3.2e+07 ! 07 7 7.808e+09 3.2e+07 1 3.2e+07 ! 08 8 7.84e+09 3.2e+07 1 3.2e+07 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 2968 visibilities. select RR ! Selecting polarization: RR, channels: 1..8 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 ! Reading IF 5 channels: 5..5 ! Reading IF 6 channels: 6..6 ! Reading IF 7 channels: 7..7 ! Reading IF 8 channels: 8..8 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/01_muppet /d2/scr/pima/bp177b_uvs/J0019+7327_X] float lim_self lim_self = 120.0 float lim_val lim_va = 1.2 integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 120 float soltime2 soltime2 = 20 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float time_av time_av = 16 float old_peak float new_peak float flux_cutoff 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. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=1.640308Jy sigma=57.975644 ! Fit after self-cal, rms=0.372009Jy sigma=12.342829 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. 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=1.111 mas, bmaj=1.39 mas, bpa=-4.289 degrees ! Estimated noise=0.611876 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.973172 Jy ! Component: 100 - total flux cleaned = 1.22103 Jy ! Total flux subtracted in 100 components = 1.22103 Jy ! Clean residual min=-0.025991 max=0.084094 Jy/beam ! Clean residual mean=-0.000003 rms=0.007814 Jy/beam ! Combined flux in latest and established models = 1.22104 Jy ! Performing phase self-cal ! Adding 3 model components to the UV plane model. ! The established model now contains 3 components and 1.22104 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.903492Jy sigma=30.374142 ! Fit after self-cal, rms=0.903420Jy sigma=30.373713 ! 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.334 mas, bmaj=1.906 mas, bpa=-18.38 degrees ! Estimated noise=0.387959 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.444782 Jy ! Component: 100 - total flux cleaned = -0.542127 Jy ! Component: 150 - total flux cleaned = -0.574632 Jy ! Component: 200 - total flux cleaned = -0.583749 Jy ! Total flux subtracted in 200 components = -0.583749 Jy ! Clean residual min=-0.099657 max=0.111860 Jy/beam ! Clean residual mean=0.000035 rms=0.024237 Jy/beam ! Combined flux in latest and established models = 0.637286 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 13 components and 0.637286 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.874988Jy sigma=31.419552 ! Fit after self-cal, rms=0.945477Jy sigma=33.978078 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.224231 Jy ! Component: 100 - total flux cleaned = 0.29333 Jy ! Component: 150 - total flux cleaned = 0.30589 Jy ! Component: 200 - total flux cleaned = 0.304878 Jy ! Total flux subtracted in 200 components = 0.304878 Jy ! Clean residual min=-0.138430 max=0.160681 Jy/beam ! Clean residual mean=0.000076 rms=0.034355 Jy/beam ! Combined flux in latest and established models = 0.942164 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 20 components and 0.942164 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=1.008322Jy sigma=35.415890 ! Fit after self-cal, rms=1.007929Jy sigma=35.400209 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! restore ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.334 x 1.906 at -18.38 degrees (North through East) ! Clean map min=-0.16365 max=0.77119 Jy/beam if ( peak(flux,max)/imstat(rms) < lim_self ) selflim lim_val ! Selfcal amplitude corrections will be limited to 0.833333 -> 1.2. ! Selfcal phase corrections will not be limited. gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! MK 1.00* NL 1.00* OV 1.00* PT 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! NL 1.00* OV 1.00* PT 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! NL 1.00* OV 1.00* PT 1.00* SC 1.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! MK 1.00* NL 1.00* OV 1.00* PT 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! NL 1.00* OV 1.00* PT 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! NL 1.00* OV 1.00* PT 1.00* SC 1.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! MK 1.00* NL 1.00* OV 1.00* PT 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! NL 1.00* OV 1.00* PT 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! NL 1.00* OV 1.00* PT 1.00* SC 1.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! MK 1.00* NL 1.00* OV 1.00* PT 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! NL 1.00* OV 1.00* PT 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! NL 1.00* OV 1.00* PT 1.00* SC 1.00* ! ! ! Correcting IF 5. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! MK 1.00* NL 1.00* OV 1.00* PT 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.11 FD 1.13 HN 1.07 LA 1.11 ! NL 1.06 OV 1.02 PT 0.99 ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.99 FD 0.99 HN 0.96 LA 0.96 ! NL 0.94 OV 0.93 PT 0.89 SC 1.04 ! ! ! Correcting IF 6. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! MK 1.00* NL 1.00* OV 1.00* PT 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! NL 1.00* OV 1.00* PT 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! NL 1.00* OV 1.00* PT 1.00* SC 1.00* ! ! ! Correcting IF 7. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! MK 1.00* NL 1.00* OV 1.00* PT 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! NL 1.00* OV 1.00* PT 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! NL 1.00* OV 1.00* PT 1.00* SC 1.00* ! ! ! Correcting IF 8. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! MK 1.00* NL 1.00* OV 1.00* PT 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! NL 1.00* OV 1.00* PT 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! NL 1.00* OV 1.00* PT 1.00* SC 1.00* ! ! ! Fit before self-cal, rms=1.007929Jy sigma=35.400209 ! Fit after self-cal, rms=1.003385Jy sigma=35.368200 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.335 mas, bmaj=1.909 mas, bpa=-18.56 degrees ! Estimated noise=0.386628 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.0728021 Jy ! Component: 100 - total flux cleaned = -0.101708 Jy ! Component: 150 - total flux cleaned = -0.106844 Jy ! Component: 200 - total flux cleaned = -0.108712 Jy ! Total flux subtracted in 200 components = -0.108712 Jy ! Clean residual min=-0.156512 max=0.176197 Jy/beam ! Clean residual mean=0.000129 rms=0.038047 Jy/beam ! Combined flux in latest and established models = 0.833452 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 27 components and 0.833452 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.995819Jy sigma=35.295735 ! Fit after self-cal, rms=0.995044Jy sigma=35.261828 ! Inverting map selfcal true, true, soltime1 ! Performing amp+phase self-cal over 120 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.995044Jy sigma=35.261828 ! Fit after self-cal, rms=0.995007Jy sigma=35.260022 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.335 mas, bmaj=1.91 mas, bpa=-18.55 degrees ! Estimated noise=0.386254 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0063333 Jy ! Component: 100 - total flux cleaned = -0.00863054 Jy ! Component: 150 - total flux cleaned = -0.00954191 Jy ! Component: 200 - total flux cleaned = -0.0102465 Jy ! Total flux subtracted in 200 components = -0.0102465 Jy ! Clean residual min=-0.157793 max=0.179498 Jy/beam ! Clean residual mean=0.000141 rms=0.037868 Jy/beam ! Combined flux in latest and established models = 0.823206 Jy selfcal ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 32 components and 0.823206 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.992133Jy sigma=35.174311 ! Fit after self-cal, rms=0.989507Jy sigma=35.061555 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (-42.8, -26.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0562939 Jy ! Component: 100 - total flux cleaned = 0.0195436 Jy ! Component: 150 - total flux cleaned = -0.00121444 Jy ! Component: 200 - total flux cleaned = -0.00869381 Jy ! Total flux subtracted in 200 components = -0.00869381 Jy ! Clean residual min=-0.162676 max=0.168751 Jy/beam ! Clean residual mean=0.000171 rms=0.035590 Jy/beam ! Combined flux in latest and established models = 0.814512 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 46 components and 0.814512 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.954457Jy sigma=33.756110 ! Fit after self-cal, rms=0.934963Jy sigma=32.992360 ! Inverting map ! Added new window around map position (-16.8, 5.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0775029 Jy ! Component: 100 - total flux cleaned = 0.0866972 Jy ! Component: 150 - total flux cleaned = 0.0915404 Jy ! Component: 200 - total flux cleaned = 0.091566 Jy ! Total flux subtracted in 200 components = 0.091566 Jy ! Clean residual min=-0.144237 max=0.157398 Jy/beam ! Clean residual mean=0.000269 rms=0.032375 Jy/beam ! Combined flux in latest and established models = 0.906078 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 62 components and 0.906078 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.871875Jy sigma=30.841740 ! Fit after self-cal, rms=0.826255Jy sigma=29.143935 ! Inverting map ! Added new window around map position (31, -43.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0106275 Jy ! Component: 100 - total flux cleaned = 0.0615479 Jy ! Component: 150 - total flux cleaned = 0.0811189 Jy ! Component: 200 - total flux cleaned = 0.0950352 Jy ! Total flux subtracted in 200 components = 0.0950352 Jy ! Clean residual min=-0.121307 max=0.130339 Jy/beam ! Clean residual mean=0.000048 rms=0.028250 Jy/beam ! Combined flux in latest and established models = 1.00111 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 80 components and 1.00111 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.776527Jy sigma=27.648103 ! Fit after self-cal, rms=0.748695Jy sigma=26.495461 ! Inverting map ! Added new window around map position (26.6, 19.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.103244 Jy ! Component: 100 - total flux cleaned = 0.132938 Jy ! Component: 150 - total flux cleaned = 0.14122 Jy ! Component: 200 - total flux cleaned = 0.151174 Jy ! Total flux subtracted in 200 components = 0.151174 Jy ! Clean residual min=-0.097735 max=0.113819 Jy/beam ! Clean residual mean=-0.000183 rms=0.025425 Jy/beam ! Combined flux in latest and established models = 1.15229 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 107 components and 1.15229 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.715012Jy sigma=25.333293 ! Fit after self-cal, rms=0.702602Jy sigma=24.756556 ! Inverting map ! Added new window around map position (36.2, 2.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0694561 Jy ! Component: 100 - total flux cleaned = 0.0687806 Jy ! Component: 150 - total flux cleaned = 0.0662488 Jy ! Component: 200 - total flux cleaned = 0.0662061 Jy ! Total flux subtracted in 200 components = 0.0662061 Jy ! Clean residual min=-0.090646 max=0.101063 Jy/beam ! Clean residual mean=-0.000017 rms=0.024060 Jy/beam ! Combined flux in latest and established models = 1.21849 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 133 components and 1.21849 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.680800Jy sigma=23.932437 ! Fit after self-cal, rms=0.678828Jy sigma=23.881559 ! Inverting map ! Added new window around map position (49, -45). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0754221 Jy ! Component: 100 - total flux cleaned = 0.073072 Jy ! Component: 150 - total flux cleaned = 0.0665627 Jy ! Component: 200 - total flux cleaned = 0.0665412 Jy ! Total flux subtracted in 200 components = 0.0665412 Jy ! Clean residual min=-0.089179 max=0.097064 Jy/beam ! Clean residual mean=-0.000070 rms=0.023054 Jy/beam ! Combined flux in latest and established models = 1.28504 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 158 components and 1.28504 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.660819Jy sigma=23.124252 ! Fit after self-cal, rms=0.655845Jy sigma=22.968804 ! Inverting map selfcal true, true, soltime2 ! Performing amp+phase self-cal over 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.655845Jy sigma=22.968804 ! Fit after self-cal, rms=0.655845Jy sigma=22.968804 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.335 mas, bmaj=1.91 mas, bpa=-18.55 degrees ! Estimated noise=0.386254 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0343768 Jy ! Component: 100 - total flux cleaned = 0.0425872 Jy ! Component: 150 - total flux cleaned = 0.0444208 Jy ! Component: 200 - total flux cleaned = 0.0496182 Jy ! Total flux subtracted in 200 components = 0.0496182 Jy ! Clean residual min=-0.085019 max=0.091991 Jy/beam ! Clean residual mean=-0.000152 rms=0.022014 Jy/beam ! Combined flux in latest and established models = 1.33465 Jy selfcal ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 172 components and 1.33465 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.643803Jy sigma=22.501259 ! Fit after self-cal, rms=0.640488Jy sigma=22.380118 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 (41.4, 16.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0790123 Jy ! Component: 100 - total flux cleaned = 0.106981 Jy ! Component: 150 - total flux cleaned = 0.112715 Jy ! Component: 200 - total flux cleaned = 0.109414 Jy ! Total flux subtracted in 200 components = 0.109414 Jy ! Clean residual min=-0.082937 max=0.086786 Jy/beam ! Clean residual mean=-0.000226 rms=0.021006 Jy/beam ! Combined flux in latest and established models = 1.44407 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 196 components and 1.44407 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.626091Jy sigma=21.799755 ! Fit after self-cal, rms=0.624695Jy sigma=21.741964 ! 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 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.619980Jy sigma=21.608524 ! Fit after self-cal, rms=0.617404Jy sigma=21.524783 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.219 mas, bmaj=3.675 mas, bpa=-15.09 degrees ! Estimated noise=0.440397 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0163345 Jy ! Component: 100 - total flux cleaned = 0.0200713 Jy ! Component: 150 - total flux cleaned = 0.0245853 Jy ! Component: 200 - total flux cleaned = 0.0274342 Jy ! Total flux subtracted in 200 components = 0.0274342 Jy ! Clean residual min=-0.097456 max=0.097981 Jy/beam ! Clean residual mean=-0.000341 rms=0.024171 Jy/beam ! Combined flux in latest and established models = 1.4715 Jy selfcal ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 223 components and 1.4715 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.608838Jy sigma=21.208912 ! Fit after self-cal, rms=0.607392Jy sigma=21.145750 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.749679 Jy ! Component: 100 - total flux cleaned = 1.08421 Jy ! Component: 150 - total flux cleaned = 1.29283 Jy ! Component: 200 - total flux cleaned = 1.42661 Jy ! Total flux subtracted in 200 components = 1.42661 Jy ! Clean residual min=-0.121328 max=0.135027 Jy/beam ! Clean residual mean=-0.000391 rms=0.034689 Jy/beam ! Combined flux in latest and established models = 2.89811 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 241 components and 2.89811 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.829070Jy sigma=29.982608 ! Fit after self-cal, rms=0.613083Jy sigma=21.768841 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.346156 Jy ! Component: 100 - total flux cleaned = 0.360717 Jy ! Component: 150 - total flux cleaned = 0.357229 Jy ! Component: 200 - total flux cleaned = 0.369943 Jy ! Total flux subtracted in 200 components = 0.369943 Jy ! Clean residual min=-0.165776 max=0.082885 Jy/beam ! Clean residual mean=-0.000084 rms=0.025686 Jy/beam ! Combined flux in latest and established models = 3.26805 Jy ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 269 components and 3.26805 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.950813Jy sigma=34.667564 ! Fit after self-cal, rms=0.906480Jy sigma=33.172225 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "************** FINISHED AMPLITUDE SELF-CAL **************" ! ************** FINISHED AMPLITUDE SELF-CAL ************** restore ! Inverting map and beam ! Estimated beam: bmin=1.334 mas, bmaj=1.913 mas, bpa=-19 degrees ! Estimated noise=0.383109 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.334 x 1.913 at -19 degrees (North through East) ! Clean map min=-0.36088 max=0.89469 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. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.435975Jy sigma=15.848608 ! Fit after self-cal, rms=0.432148Jy sigma=15.634124 ! Inverting map and beam ! Estimated beam: bmin=1.331 mas, bmaj=1.902 mas, bpa=-18.67 degrees ! Estimated noise=0.387188 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.219021 Jy ! Component: 100 - total flux cleaned = -0.33152 Jy ! Component: 150 - total flux cleaned = -0.39775 Jy ! Component: 200 - total flux cleaned = -0.422935 Jy ! Total flux subtracted in 200 components = -0.422935 Jy ! Clean residual min=-0.092456 max=0.046598 Jy/beam ! Clean residual mean=0.000064 rms=0.012628 Jy/beam ! Combined flux in latest and established models = 2.84512 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 283 components and 2.84512 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.312173Jy sigma=10.772118 ! Fit after self-cal, rms=0.286703Jy sigma=9.934282 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0687515 Jy ! Component: 100 - total flux cleaned = -0.0959934 Jy ! Component: 150 - total flux cleaned = -0.118101 Jy ! Component: 200 - total flux cleaned = -0.133006 Jy ! Total flux subtracted in 200 components = -0.133006 Jy ! Clean residual min=-0.042102 max=0.034888 Jy/beam ! Clean residual mean=0.000112 rms=0.008927 Jy/beam ! Combined flux in latest and established models = 2.71211 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 297 components and 2.71211 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.527184Jy sigma=18.156570 ! Fit after self-cal, rms=0.520854Jy sigma=17.919167 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.339836 Jy ! Component: 100 - total flux cleaned = 0.420414 Jy ! Component: 150 - total flux cleaned = 0.429012 Jy ! Component: 200 - total flux cleaned = 0.433342 Jy ! Total flux subtracted in 200 components = 0.433342 Jy ! Clean residual min=-0.050682 max=0.038324 Jy/beam ! Clean residual mean=0.000137 rms=0.009905 Jy/beam ! Combined flux in latest and established models = 3.14545 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 313 components and 3.14545 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.360537Jy sigma=12.819112 ! Fit after self-cal, rms=0.354172Jy sigma=12.578028 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.354172Jy sigma=12.578028 ! Fit after self-cal, rms=0.342895Jy sigma=12.159481 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=1.108 mas, bmaj=1.379 mas, bpa=-3.532 degrees ! Estimated noise=0.612483 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.964014 Jy ! Component: 100 - total flux cleaned = 1.2163 Jy ! Component: 150 - total flux cleaned = 1.32851 Jy ! Component: 200 - total flux cleaned = 1.39743 Jy ! Total flux subtracted in 200 components = 1.39743 Jy ! Clean residual min=-0.064924 max=0.048953 Jy/beam ! Clean residual mean=0.000011 rms=0.011434 Jy/beam ! Combined flux in latest and established models = 1.39743 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 29 model components to the UV plane model. ! The established model now contains 29 components and 1.39743 Jy ! Inverting map and beam ! Estimated beam: bmin=1.332 mas, bmaj=1.894 mas, bpa=-18.66 degrees ! Estimated noise=0.388513 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam keep new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0386511 Jy ! Component: 100 - total flux cleaned = 0.0519187 Jy ! Component: 150 - total flux cleaned = 0.0577958 Jy ! Component: 200 - total flux cleaned = 0.0618474 Jy ! Total flux subtracted in 200 components = 0.0618474 Jy ! Clean residual min=-0.075576 max=0.038766 Jy/beam ! Clean residual mean=0.000064 rms=0.010900 Jy/beam ! Combined flux in latest and established models = 1.45928 Jy ! Adding 32 model components to the UV plane model. ! The established model now contains 57 components and 1.45928 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.212 mas, bmaj=3.658 mas, bpa=-14.87 degrees ! Estimated noise=0.448597 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0139051 Jy ! Component: 100 - total flux cleaned = 0.0163968 Jy ! Component: 150 - total flux cleaned = 0.0163642 Jy ! Component: 200 - total flux cleaned = 0.0183932 Jy ! Total flux subtracted in 200 components = 0.0183932 Jy ! Clean residual min=-0.086313 max=0.040166 Jy/beam ! Clean residual mean=0.000038 rms=0.013455 Jy/beam ! Combined flux in latest and established models = 1.47767 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=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00382674 Jy ! Component: 100 - total flux cleaned = 0.00836349 Jy ! Component: 150 - total flux cleaned = 0.0126753 Jy ! Component: 200 - total flux cleaned = 0.0167582 Jy ! Total flux subtracted in 200 components = 0.0167582 Jy ! Clean residual min=-0.082022 max=0.037869 Jy/beam ! Clean residual mean=0.000020 rms=0.012896 Jy/beam ! Combined flux in latest and established models = 1.49443 Jy ! Adding 70 model components to the UV plane model. ! The established model now contains 121 components and 1.49443 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=1.332 mas, bmaj=1.894 mas, bpa=-18.66 degrees ! Estimated noise=0.388513 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0113922 Jy ! Component: 100 - total flux cleaned = 0.0183899 Jy ! Component: 150 - total flux cleaned = 0.022433 Jy ! Component: 200 - total flux cleaned = 0.0232031 Jy ! Total flux subtracted in 200 components = 0.0232031 Jy ! Clean residual min=-0.063472 max=0.035169 Jy/beam ! Clean residual mean=-0.000003 rms=0.009538 Jy/beam ! Combined flux in latest and established models = 1.51763 Jy selfcal ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 151 components and 1.51763 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.260318Jy sigma=8.899809 ! Fit after self-cal, rms=0.226427Jy sigma=7.638771 wmodel /d2/scr/pima/bp177b_uvs/J0019+7327_X_map.mod ! Writing 151 model components to file: /d2/scr/pima/bp177b_uvs/J0019+7327_X_map.mod wobs /d2/scr/pima/bp177b_uvs/J0019+7327_X_uvs.fits ! Writing UV FITS file: /d2/scr/pima/bp177b_uvs/J0019+7327_X_uvs.fits wwins /d2/scr/pima/bp177b_uvs/J0019+7327_X_map.win ! wwins: Wrote 8 windows to /d2/scr/pima/bp177b_uvs/J0019+7327_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.020749 Jy ! Component: 100 - total flux cleaned = -0.021916 Jy ! Component: 150 - total flux cleaned = -0.0264538 Jy ! Component: 200 - total flux cleaned = -0.0315351 Jy ! Component: 250 - total flux cleaned = -0.0359153 Jy ! Component: 300 - total flux cleaned = -0.0392483 Jy ! Component: 350 - total flux cleaned = -0.0416404 Jy ! Component: 400 - total flux cleaned = -0.0435516 Jy ! Component: 450 - total flux cleaned = -0.0450328 Jy ! Component: 500 - total flux cleaned = -0.0461149 Jy ! Component: 550 - total flux cleaned = -0.0475188 Jy ! Component: 600 - total flux cleaned = -0.0481953 Jy ! Component: 650 - total flux cleaned = -0.048867 Jy ! Component: 700 - total flux cleaned = -0.0498536 Jy ! Component: 750 - total flux cleaned = -0.0508184 Jy ! Component: 800 - total flux cleaned = -0.0511294 Jy ! Component: 850 - total flux cleaned = -0.0511297 Jy ! Component: 900 - total flux cleaned = -0.051131 Jy ! Component: 950 - total flux cleaned = -0.0517244 Jy ! Component: 1000 - total flux cleaned = -0.0517223 Jy ! Component: 1050 - total flux cleaned = -0.0508519 Jy ! Component: 1100 - total flux cleaned = -0.0497204 Jy ! Component: 1150 - total flux cleaned = -0.0497261 Jy ! Component: 1200 - total flux cleaned = -0.0494499 Jy ! Component: 1250 - total flux cleaned = -0.0483622 Jy ! Component: 1300 - total flux cleaned = -0.0483612 Jy ! Component: 1350 - total flux cleaned = -0.0473003 Jy ! Component: 1400 - total flux cleaned = -0.0465181 Jy ! Component: 1450 - total flux cleaned = -0.0457428 Jy ! Component: 1500 - total flux cleaned = -0.0449741 Jy ! Component: 1550 - total flux cleaned = -0.0434662 Jy ! Component: 1600 - total flux cleaned = -0.0422286 Jy ! Component: 1650 - total flux cleaned = -0.0417396 Jy ! Component: 1700 - total flux cleaned = -0.0405308 Jy ! Component: 1750 - total flux cleaned = -0.0400536 Jy ! Component: 1800 - total flux cleaned = -0.0400526 Jy ! Component: 1850 - total flux cleaned = -0.0391159 Jy ! Component: 1900 - total flux cleaned = -0.0391142 Jy ! Component: 1950 - total flux cleaned = -0.0382023 Jy ! Component: 2000 - total flux cleaned = -0.0373011 Jy ! Component: 2050 - total flux cleaned = -0.0379673 Jy ! Component: 2100 - total flux cleaned = -0.0368659 Jy ! Component: 2150 - total flux cleaned = -0.0364312 Jy ! Component: 2200 - total flux cleaned = -0.0368614 Jy ! Component: 2250 - total flux cleaned = -0.0366483 Jy ! Component: 2300 - total flux cleaned = -0.0364403 Jy ! Component: 2350 - total flux cleaned = -0.0368582 Jy ! Component: 2400 - total flux cleaned = -0.0358219 Jy ! Component: 2450 - total flux cleaned = -0.0356159 Jy ! Component: 2500 - total flux cleaned = -0.034601 Jy ! Component: 2550 - total flux cleaned = -0.0348047 Jy ! Component: 2600 - total flux cleaned = -0.0340083 Jy ! Component: 2650 - total flux cleaned = -0.0336168 Jy ! Component: 2700 - total flux cleaned = -0.0340074 Jy ! Component: 2750 - total flux cleaned = -0.0328427 Jy ! Component: 2800 - total flux cleaned = -0.0322679 Jy ! Component: 2850 - total flux cleaned = -0.0322678 Jy ! Component: 2900 - total flux cleaned = -0.0315137 Jy ! Component: 2950 - total flux cleaned = -0.0313294 Jy ! Component: 3000 - total flux cleaned = -0.0307724 Jy ! Component: 3050 - total flux cleaned = -0.0311401 Jy ! Component: 3100 - total flux cleaned = -0.0307732 Jy ! Component: 3150 - total flux cleaned = -0.0300481 Jy ! Component: 3200 - total flux cleaned = -0.0298664 Jy ! Component: 3250 - total flux cleaned = -0.0296874 Jy ! Component: 3300 - total flux cleaned = -0.0293334 Jy ! Component: 3350 - total flux cleaned = -0.0284561 Jy ! Component: 3400 - total flux cleaned = -0.0289786 Jy ! Component: 3450 - total flux cleaned = -0.0279417 Jy ! Component: 3500 - total flux cleaned = -0.0274278 Jy ! Component: 3550 - total flux cleaned = -0.0279376 Jy ! Component: 3600 - total flux cleaned = -0.0269252 Jy ! Component: 3650 - total flux cleaned = -0.0265901 Jy ! Component: 3700 - total flux cleaned = -0.0272543 Jy ! Component: 3750 - total flux cleaned = -0.0269245 Jy ! Component: 3800 - total flux cleaned = -0.0269239 Jy ! Component: 3850 - total flux cleaned = -0.026924 Jy ! Component: 3900 - total flux cleaned = -0.0270838 Jy ! Component: 3950 - total flux cleaned = -0.0280441 Jy ! Component: 4000 - total flux cleaned = -0.0272478 Jy ! Component: 4050 - total flux cleaned = -0.0274062 Jy ! Total flux subtracted in 4096 components = -0.0277179 Jy ! Clean residual min=-0.008808 max=0.012730 Jy/beam ! Clean residual mean=-0.000009 rms=0.002887 Jy/beam ! Combined flux in latest and established models = 1.48991 Jy keep ! Adding 335 model components to the UV plane model. ! The established model now contains 483 components and 1.48991 Jy ![Exited script file: /home/lpetrov/bin/01_muppet] wmap "/d2/scr/pima/bp177b_uvs/J0019+7327_X\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.332 x 1.894 at -18.66 degrees (North through East) ! Clean map min=-0.82809 max=2.1925 Jy/beam ! Writing clean map to FITS file: /d2/scr/pima/bp177b_uvs/J0019+7327_X_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 2.19249 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.0806924 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -11.0412 11.0412 22.0824 44.1648 88.3296 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.242 print "peak= ", int(1000*peak_flux)/1000 ! peak= 2.192 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 27.17 quit ! Quitting program ! Log file /d2/scr/pima/bp177b_uvs/J0019+7327_X_uvs.log closed on Sat Mar 22 08:06:31 2014