! Started logfile: /scr/pima/rv119_uvs/J1800+7828_S_dfm.log on Sun Dec 17 11:38:19 2017 float field_size field_size = 1024 float field_cell field_cell = 0.5000 float taper_size taper_size = 20.0 float freq freq = 2.23275 obs /scr/pima/rv119_uvs/J1800+7828_S_uva.fits ! Reading UV FITS file: /scr/pima/rv119_uvs/J1800+7828_S_uva.fits ! AN table 1: 2 integrations on 3 of 3 possible baselines. ! AN table 2: 2 integrations on 6 of 6 possible baselines. ! AN table 3: 2 integrations on 3 of 3 possible baselines. ! AN table 4: 4 integrations on 3 of 3 possible baselines. ! AN table 5: 6 integrations on 45 of 45 possible baselines. ! AN table 6: 48 integrations on 55 of 55 possible baselines. ! Apparent sampling: 0.388323 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1800+7828 ! ! 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.23275e+09 8e+06 1 8e+06 ! 02 2 2.26275e+09 8e+06 1 8e+06 ! 03 3 2.35275e+09 8e+06 1 8e+06 ! 04 4 2.37275e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 79 lines of history. ! ! Reading 4576 visibilities. select rr, 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16 ! Restricting channel ranges to the available 4 channels. ! Selecting polarization: RR, channels: 1..4 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 print "field_size =", field_size ! field_size = 1024 print "field_cell =", field_cell ! field_cell = 0.5 print "taper_size =", taper_size ! taper_size = 20 ![@/opt64/bin/pima_mupet_01.dfm J1800+7828_S] 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.500x0.500 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 34 telescope corrections were flagged in sub-array 5. ! A total of 238 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 34 telescope corrections were flagged in sub-array 5. ! A total of 238 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 2 telescope corrections were flagged in sub-array 2. ! A total of 36 telescope corrections were flagged in sub-array 5. ! A total of 270 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 34 telescope corrections were flagged in sub-array 5. ! A total of 238 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=1.692699Jy sigma=20.739267 ! Fit after self-cal, rms=0.668217Jy sigma=7.991724 ! 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=3.034 mas, bmaj=3.325 mas, bpa=69.66 degrees ! Estimated noise=1.32693 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 = 1.0092 Jy ! Component: 100 - total flux cleaned = 1.2564 Jy ! Total flux subtracted in 100 components = 1.2564 Jy ! Clean residual min=-0.036128 max=0.093420 Jy/beam ! Clean residual mean=-0.000322 rms=0.009547 Jy/beam ! Combined flux in latest and established models = 1.2564 Jy ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 8 components and 1.2564 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.605836Jy sigma=8.048382 ! Fit after self-cal, rms=0.605420Jy sigma=8.038028 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0914898 Jy ! Component: 100 - total flux cleaned = 0.129541 Jy ! Total flux subtracted in 100 components = 0.129541 Jy ! Clean residual min=-0.032352 max=0.051654 Jy/beam ! Clean residual mean=-0.000260 rms=0.009021 Jy/beam ! Combined flux in latest and established models = 1.38594 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 17 components and 1.38594 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.602451Jy sigma=8.432673 ! Fit after self-cal, rms=0.602417Jy sigma=8.432365 ! 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=3.609 mas, bmaj=4.594 mas, bpa=-72.22 degrees ! Estimated noise=0.880551 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.12065 Jy ! Component: 100 - total flux cleaned = -0.148085 Jy ! Component: 150 - total flux cleaned = -0.148117 Jy ! Component: 200 - total flux cleaned = -0.149449 Jy ! Total flux subtracted in 200 components = -0.149449 Jy ! Clean residual min=-0.100132 max=0.223448 Jy/beam ! Clean residual mean=-0.001963 rms=0.038799 Jy/beam ! Combined flux in latest and established models = 1.23649 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 36 components and 1.23649 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.600479Jy sigma=7.860346 ! Fit after self-cal, rms=0.600356Jy sigma=7.858798 ! Inverting map ! Added new window around map position (-15.5, -2.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.198622 Jy ! Component: 100 - total flux cleaned = 0.287276 Jy ! Component: 150 - total flux cleaned = 0.320245 Jy ! Component: 200 - total flux cleaned = 0.333623 Jy ! Total flux subtracted in 200 components = 0.333623 Jy ! Clean residual min=-0.107696 max=0.168271 Jy/beam ! Clean residual mean=-0.000799 rms=0.029785 Jy/beam ! Combined flux in latest and established models = 1.57011 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 69 components and 1.57011 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.556327Jy sigma=6.657251 ! Fit after self-cal, rms=0.545770Jy sigma=6.513342 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! restore ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.609 x 4.594 at -72.22 degrees (North through East) ! Clean map min=-0.10027 max=1.1758 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 6 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! MK 0.00* NY 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! HN 0.93 NL 1.13 ON 0.94 SC 0.94 ! ! A total of 6 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! HN 0.00* KP 0.00* ON 0.00* ! ! A total of 12 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! NY 0.00* ON 0.00* WZ 0.00* ! ! A total of 6 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 0.77 FD 0.95 KP 0.99 LA 0.93 ! MK 1.08 NL 0.99 NY 1.11 ON 1.22 ! OV 0.94 PT 0.93 ! ! A total of 6 telescope corrections were flagged in sub-array 6. ! Telescope amplitude corrections in sub-array 6: ! BR 0.98 FD 0.96 HN 1.01 KP 0.89 ! LA 0.94 MK 1.04 NL 0.99 NY 1.02 ! OV 1.88 PT 0.96 SC 1.01 ! ! ! Correcting IF 2. ! A total of 6 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! MK 0.00* NY 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! HN 0.90 NL 1.13 ON 0.89 SC 0.96 ! ! A total of 6 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! HN 0.00* KP 0.00* ON 0.00* ! ! A total of 12 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! NY 0.00* ON 0.00* WZ 0.00* ! ! A total of 6 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 0.77 FD 0.97 KP 0.96 LA 0.91 ! MK 1.12 NL 1.01 NY 1.16 ON 1.22 ! OV 0.93 PT 0.91 ! ! A total of 6 telescope corrections were flagged in sub-array 6. ! Telescope amplitude corrections in sub-array 6: ! BR 1.00 FD 0.98 HN 1.01 KP 0.90 ! LA 0.94 MK 1.05 NL 0.97 NY 1.05 ! OV 1.75 PT 0.95 SC 1.02 ! ! ! Correcting IF 3. ! A total of 6 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! MK 0.00* NY 0.00* WZ 0.00* ! ! A total of 6 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! HN 0.00* NL 0.00* ON 0.00* SC 0.00* ! ! A total of 6 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! HN 0.00* KP 0.00* ON 0.00* ! ! A total of 12 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! NY 0.00* ON 0.00* WZ 0.00* ! ! A total of 6 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 0.79 FD 0.94 KP 0.83 LA 0.93 ! MK 0.76 NL 1.00* NY 1.24 ON 1.14 ! OV 0.95 PT 0.93 ! ! A total of 12 telescope corrections were flagged in sub-array 6. ! Telescope amplitude corrections in sub-array 6: ! BR 1.04 FD 0.94 HN 1.05 KP 0.73 ! LA 0.95 MK 0.78 NL 1.00* NY 1.11 ! OV 1.72 PT 0.95 SC 0.91 ! ! ! Correcting IF 4. ! A total of 6 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! MK 0.00* NY 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! HN 0.92 NL 1.08 ON 0.82 SC 1.03 ! ! A total of 6 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! HN 0.00* KP 0.00* ON 0.00* ! ! A total of 12 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! NY 0.00* ON 0.00* WZ 0.00* ! ! A total of 6 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 0.76 FD 0.97 KP 0.99 LA 0.96 ! MK 1.09 NL 1.01 NY 1.22 ON 1.19 ! OV 0.97 PT 0.87 ! ! A total of 6 telescope corrections were flagged in sub-array 6. ! Telescope amplitude corrections in sub-array 6: ! BR 1.00 FD 0.95 HN 1.01 KP 0.93 ! LA 0.94 MK 1.05 NL 0.97 NY 1.01 ! OV 1.89 PT 0.96 SC 1.01 ! ! ! Fit before self-cal, rms=0.545770Jy sigma=6.513342 ! Fit after self-cal, rms=0.635159Jy sigma=5.405312 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.602 mas, bmaj=4.649 mas, bpa=-73.81 degrees ! Estimated noise=0.948351 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.000660625 Jy ! Component: 100 - total flux cleaned = 0.00313311 Jy ! Component: 150 - total flux cleaned = 0.00846175 Jy ! Component: 200 - total flux cleaned = 0.0130148 Jy ! Total flux subtracted in 200 components = 0.0130148 Jy ! Clean residual min=-0.073540 max=0.068978 Jy/beam ! Clean residual mean=-0.000199 rms=0.018346 Jy/beam ! Combined flux in latest and established models = 1.58313 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 84 components and 1.58313 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.627741Jy sigma=5.367647 ! Fit after self-cal, rms=0.627748Jy sigma=5.362499 ! 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. ! ! Fit before self-cal, rms=0.627748Jy sigma=5.362499 ! Fit after self-cal, rms=1.338206Jy sigma=2.604862 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.6 mas, bmaj=4.583 mas, bpa=-74.53 degrees ! Estimated noise=0.969859 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0416723 Jy ! Component: 100 - total flux cleaned = -0.0447904 Jy ! Component: 150 - total flux cleaned = -0.0450706 Jy ! Component: 200 - total flux cleaned = -0.045152 Jy ! Total flux subtracted in 200 components = -0.045152 Jy ! Clean residual min=-0.020447 max=0.029204 Jy/beam ! Clean residual mean=-0.000435 rms=0.004480 Jy/beam ! Combined flux in latest and established models = 1.53797 Jy selfcal ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 97 components and 1.53797 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.330090Jy sigma=2.327591 ! Fit after self-cal, rms=1.329503Jy sigma=2.284491 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.0111193 Jy ! Component: 100 - total flux cleaned = -0.0111318 Jy ! Component: 150 - total flux cleaned = -0.0136175 Jy ! Component: 200 - total flux cleaned = -0.0158392 Jy ! Total flux subtracted in 200 components = -0.0158392 Jy ! Clean residual min=-0.017924 max=0.023503 Jy/beam ! Clean residual mean=-0.000506 rms=0.004020 Jy/beam ! Combined flux in latest and established models = 1.52213 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 103 components and 1.52214 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.326890Jy sigma=2.222083 ! Fit after self-cal, rms=1.326588Jy sigma=2.215334 ! Inverting map ! Added new window around map position (-45.5, -3). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0121418 Jy ! Component: 100 - total flux cleaned = 0.0126814 Jy ! Component: 150 - total flux cleaned = 0.0126758 Jy ! Component: 200 - total flux cleaned = 0.0130976 Jy ! Total flux subtracted in 200 components = 0.0130976 Jy ! Clean residual min=-0.015771 max=0.020987 Jy/beam ! Clean residual mean=-0.000457 rms=0.003694 Jy/beam ! Combined flux in latest and established models = 1.53523 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 121 components and 1.53523 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.324857Jy sigma=2.175548 ! Fit after self-cal, rms=1.324697Jy sigma=2.170919 ! Inverting map ! Added new window around map position (13.5, 3). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0129767 Jy ! Component: 100 - total flux cleaned = 0.0153542 Jy ! Component: 150 - total flux cleaned = 0.0153498 Jy ! Component: 200 - total flux cleaned = 0.015348 Jy ! Total flux subtracted in 200 components = 0.015348 Jy ! Clean residual min=-0.014957 max=0.019833 Jy/beam ! Clean residual mean=-0.000389 rms=0.003454 Jy/beam ! Combined flux in latest and established models = 1.55058 Jy ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 137 components and 1.55058 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.324450Jy sigma=2.135020 ! Fit after self-cal, rms=1.324337Jy sigma=2.130217 ! Inverting map ! Added new window around map position (-29.5, -7.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0159016 Jy ! Component: 100 - total flux cleaned = 0.0229505 Jy ! Component: 150 - total flux cleaned = 0.0270918 Jy ! Component: 200 - total flux cleaned = 0.0295234 Jy ! Total flux subtracted in 200 components = 0.0295234 Jy ! Clean residual min=-0.012729 max=0.014132 Jy/beam ! Clean residual mean=-0.000268 rms=0.003007 Jy/beam ! Combined flux in latest and established models = 1.5801 Jy ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 160 components and 1.5801 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.323461Jy sigma=2.095281 ! Fit after self-cal, rms=1.323384Jy sigma=2.093014 ! Inverting map ! Added new window around map position (-5.5, 0.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00702258 Jy ! Component: 100 - total flux cleaned = 0.00820838 Jy ! Component: 150 - total flux cleaned = 0.00857492 Jy ! Component: 200 - total flux cleaned = 0.0092484 Jy ! Total flux subtracted in 200 components = 0.0092484 Jy ! Clean residual min=-0.011033 max=0.012625 Jy/beam ! Clean residual mean=-0.000230 rms=0.002861 Jy/beam ! Combined flux in latest and established models = 1.58935 Jy ! Performing phase self-cal ! Adding 44 model components to the UV plane model. ! The established model now contains 186 components and 1.58935 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.322749Jy sigma=2.074036 ! Fit after self-cal, rms=1.322733Jy sigma=2.072488 ! Inverting map ! Added new window around map position (-36.5, -20). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0124474 Jy ! Component: 100 - total flux cleaned = 0.0138964 Jy ! Component: 150 - total flux cleaned = 0.0145241 Jy ! Component: 200 - total flux cleaned = 0.0154026 Jy ! Total flux subtracted in 200 components = 0.0154026 Jy ! Clean residual min=-0.010818 max=0.011383 Jy/beam ! Clean residual mean=-0.000179 rms=0.002705 Jy/beam ! Combined flux in latest and established models = 1.60476 Jy ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 204 components and 1.60475 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.322338Jy sigma=2.054719 ! Fit after self-cal, rms=1.322295Jy sigma=2.053361 ! 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. ! ! Fit before self-cal, rms=1.322295Jy sigma=2.053361 ! Fit after self-cal, rms=1.576247Jy sigma=1.949267 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.589 mas, bmaj=4.545 mas, bpa=-74.99 degrees ! Estimated noise=0.971431 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0057845 Jy ! Component: 100 - total flux cleaned = -0.00529366 Jy ! Component: 150 - total flux cleaned = -0.00387654 Jy ! Component: 200 - total flux cleaned = -0.00136422 Jy ! Total flux subtracted in 200 components = -0.00136422 Jy ! Clean residual min=-0.009683 max=0.010065 Jy/beam ! Clean residual mean=-0.000143 rms=0.002256 Jy/beam ! Combined flux in latest and established models = 1.60339 Jy selfcal ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 223 components and 1.60339 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.574077Jy sigma=1.914783 ! Fit after self-cal, rms=1.574050Jy sigma=1.912496 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 = -2.03998e-05 Jy ! Component: 100 - total flux cleaned = 0.00316348 Jy ! Component: 150 - total flux cleaned = 0.00575494 Jy ! Component: 200 - total flux cleaned = 0.00755786 Jy ! Total flux subtracted in 200 components = 0.00755786 Jy ! Clean residual min=-0.009163 max=0.009358 Jy/beam ! Clean residual mean=-0.000115 rms=0.002124 Jy/beam ! Combined flux in latest and established models = 1.61095 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 243 components and 1.61095 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.573031Jy sigma=1.901237 ! Fit after self-cal, rms=1.573011Jy sigma=1.900495 ! 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 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.573011Jy sigma=1.900495 ! Fit after self-cal, rms=1.645381Jy sigma=1.887255 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.889 mas, bmaj=9.983 mas, bpa=-78.26 degrees ! Estimated noise=1.28971 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00456096 Jy ! Component: 100 - total flux cleaned = 0.00647567 Jy ! Component: 150 - total flux cleaned = 0.0080376 Jy ! Component: 200 - total flux cleaned = 0.010036 Jy ! Total flux subtracted in 200 components = 0.010036 Jy ! Clean residual min=-0.009469 max=0.010542 Jy/beam ! Clean residual mean=-0.000216 rms=0.002745 Jy/beam ! Combined flux in latest and established models = 1.62098 Jy selfcal ! Performing phase self-cal ! Adding 52 model components to the UV plane model. ! The established model now contains 284 components and 1.62098 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.644026Jy sigma=1.876386 ! Fit after self-cal, rms=1.644015Jy sigma=1.874939 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.00387978 Jy ! Component: 100 - total flux cleaned = 0.00498842 Jy ! Component: 150 - total flux cleaned = 0.00574042 Jy ! Component: 200 - total flux cleaned = 0.00619145 Jy ! Total flux subtracted in 200 components = 0.00619145 Jy ! Clean residual min=-0.009361 max=0.009832 Jy/beam ! Clean residual mean=-0.000164 rms=0.002612 Jy/beam ! Combined flux in latest and established models = 1.62718 Jy ! Performing phase self-cal ! Adding 49 model components to the UV plane model. ! The established model now contains 314 components and 1.62718 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.643426Jy sigma=1.872955 ! Fit after self-cal, rms=1.643413Jy sigma=1.872580 ! 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=3.582 mas, bmaj=4.524 mas, bpa=-75.19 degrees ! Estimated noise=0.972156 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.582 x 4.524 at -75.19 degrees (North through East) ! Clean map min=-0.01864 max=1.2538 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=1.643413Jy sigma=1.872580 ! Fit after self-cal, rms=1.874325Jy sigma=0.927078 ! Inverting map and beam ! Estimated beam: bmin=3.581 mas, bmaj=4.521 mas, bpa=-75.27 degrees ! Estimated noise=0.970446 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000265237 Jy ! Component: 100 - total flux cleaned = 0.000574365 Jy ! Component: 150 - total flux cleaned = 0.00129419 Jy ! Component: 200 - total flux cleaned = 0.00239636 Jy ! Total flux subtracted in 200 components = 0.00239636 Jy ! Clean residual min=-0.006394 max=0.006990 Jy/beam ! Clean residual mean=-0.000075 rms=0.001480 Jy/beam ! Combined flux in latest and established models = 1.62957 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 325 components and 1.62957 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.874154Jy sigma=0.906369 ! Fit after self-cal, rms=1.874175Jy sigma=0.905372 ! Inverting map ! Added new window around map position (-35.5, 10). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00469389 Jy ! Component: 100 - total flux cleaned = 0.00597182 Jy ! Component: 150 - total flux cleaned = 0.00695582 Jy ! Component: 200 - total flux cleaned = 0.00787588 Jy ! Total flux subtracted in 200 components = 0.00787588 Jy ! Clean residual min=-0.005893 max=0.005439 Jy/beam ! Clean residual mean=-0.000049 rms=0.001369 Jy/beam ! Combined flux in latest and established models = 1.63745 Jy ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 337 components and 1.63745 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.873932Jy sigma=0.892584 ! Fit after self-cal, rms=1.873951Jy sigma=0.891947 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.873951Jy sigma=0.891947 ! Fit after self-cal, rms=1.851785Jy sigma=0.887614 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.093 mas, bmaj=3.511 mas, bpa=69.2 degrees ! Estimated noise=1.29131 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.981863 Jy ! Component: 100 - total flux cleaned = 1.21451 Jy ! Component: 150 - total flux cleaned = 1.30225 Jy ! Component: 200 - total flux cleaned = 1.35141 Jy ! Total flux subtracted in 200 components = 1.35141 Jy ! Clean residual min=-0.011921 max=0.025823 Jy/beam ! Clean residual mean=-0.000171 rms=0.002902 Jy/beam ! Combined flux in latest and established models = 1.35141 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 20 model components to the UV plane model. ! The established model now contains 20 components and 1.35141 Jy ! Inverting map and beam ! Estimated beam: bmin=3.579 mas, bmaj=4.514 mas, bpa=-75.42 degrees ! Estimated noise=0.970716 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.0512809 Jy ! Component: 100 - total flux cleaned = 0.0866995 Jy ! Component: 150 - total flux cleaned = 0.114427 Jy ! Component: 200 - total flux cleaned = 0.137369 Jy ! Total flux subtracted in 200 components = 0.137369 Jy ! Clean residual min=-0.010515 max=0.015238 Jy/beam ! Clean residual mean=-0.000653 rms=0.003070 Jy/beam ! Combined flux in latest and established models = 1.48878 Jy ! Adding 51 model components to the UV plane model. ! The established model now contains 65 components and 1.48878 Jy ! Inverting map ! Added new window around map position (-28.5, -1.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0194968 Jy ! Component: 100 - total flux cleaned = 0.0361225 Jy ! Component: 150 - total flux cleaned = 0.0505686 Jy ! Component: 200 - total flux cleaned = 0.0632831 Jy ! Total flux subtracted in 200 components = 0.0632831 Jy ! Clean residual min=-0.008252 max=0.009384 Jy/beam ! Clean residual mean=-0.000400 rms=0.002161 Jy/beam ! Combined flux in latest and established models = 1.55206 Jy ! Adding 72 model components to the UV plane model. ! The established model now contains 115 components and 1.55206 Jy ! Inverting map ! Added new window around map position (-44, 13). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0116094 Jy ! Component: 100 - total flux cleaned = 0.0198838 Jy ! Component: 150 - total flux cleaned = 0.0274375 Jy ! Component: 200 - total flux cleaned = 0.0340154 Jy ! Total flux subtracted in 200 components = 0.0340154 Jy ! Clean residual min=-0.007191 max=0.007798 Jy/beam ! Clean residual mean=-0.000269 rms=0.001731 Jy/beam ! Combined flux in latest and established models = 1.58607 Jy ! Adding 69 model components to the UV plane model. ! The established model now contains 153 components and 1.58607 Jy ! Inverting map ! Added new window around map position (-23.5, -18). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00745083 Jy ! Component: 100 - total flux cleaned = 0.0128598 Jy ! Component: 150 - total flux cleaned = 0.0178551 Jy ! Component: 200 - total flux cleaned = 0.0224853 Jy ! Total flux subtracted in 200 components = 0.0224853 Jy ! Clean residual min=-0.006149 max=0.006305 Jy/beam ! Clean residual mean=-0.000183 rms=0.001492 Jy/beam ! Combined flux in latest and established models = 1.60856 Jy ! Adding 61 model components to the UV plane model. ! The established model now contains 181 components and 1.60856 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=7.874 mas, bmaj=9.96 mas, bpa=-78.24 degrees ! Estimated noise=1.29159 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0102562 Jy ! Component: 100 - total flux cleaned = 0.0175713 Jy ! Component: 150 - total flux cleaned = 0.0218383 Jy ! Component: 200 - total flux cleaned = 0.0248533 Jy ! Total flux subtracted in 200 components = 0.0248533 Jy ! Clean residual min=-0.006508 max=0.005473 Jy/beam ! Clean residual mean=-0.000205 rms=0.001768 Jy/beam ! Combined flux in latest and established models = 1.63341 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.00226919 Jy ! Component: 100 - total flux cleaned = 0.00396506 Jy ! Component: 150 - total flux cleaned = 0.00520997 Jy ! Component: 200 - total flux cleaned = 0.00664663 Jy ! Total flux subtracted in 200 components = 0.00664663 Jy ! Clean residual min=-0.006387 max=0.005318 Jy/beam ! Clean residual mean=-0.000142 rms=0.001633 Jy/beam ! Combined flux in latest and established models = 1.64006 Jy ! Adding 109 model components to the UV plane model. ! The established model now contains 278 components and 1.64006 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.579 mas, bmaj=4.514 mas, bpa=-75.42 degrees ! Estimated noise=0.970716 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0010453 Jy ! Component: 100 - total flux cleaned = -0.00141149 Jy ! Component: 150 - total flux cleaned = -0.00141833 Jy ! Component: 200 - total flux cleaned = -0.00111923 Jy ! Total flux subtracted in 200 components = -0.00111923 Jy ! Clean residual min=-0.005439 max=0.005398 Jy/beam ! Clean residual mean=-0.000063 rms=0.001239 Jy/beam ! Combined flux in latest and established models = 1.63894 Jy selfcal ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 301 components and 1.63894 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.850375Jy sigma=0.874510 ! Fit after self-cal, rms=1.849779Jy sigma=0.860042 wmodel J1800+7828_S_map.mod ! Writing 301 model components to file: J1800+7828_S_map.mod wobs J1800+7828_S_uvs.fits ! Writing UV FITS file: J1800+7828_S_uvs.fits wwins J1800+7828_S_map.win ! wwins: Wrote 11 windows to J1800+7828_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.00066758 Jy ! Component: 100 - total flux cleaned = -0.000841149 Jy ! Component: 150 - total flux cleaned = -0.00084107 Jy ! Component: 200 - total flux cleaned = -0.000922108 Jy ! Component: 250 - total flux cleaned = -0.00115273 Jy ! Component: 300 - total flux cleaned = -0.00130211 Jy ! Component: 350 - total flux cleaned = -0.00159495 Jy ! Component: 400 - total flux cleaned = -0.00152228 Jy ! Component: 450 - total flux cleaned = -0.00166504 Jy ! Component: 500 - total flux cleaned = -0.00173453 Jy ! Component: 550 - total flux cleaned = -0.00166504 Jy ! Component: 600 - total flux cleaned = -0.00166523 Jy ! Component: 650 - total flux cleaned = -0.00159823 Jy ! Component: 700 - total flux cleaned = -0.00166475 Jy ! Component: 750 - total flux cleaned = -0.0015981 Jy ! Component: 800 - total flux cleaned = -0.00159757 Jy ! Component: 850 - total flux cleaned = -0.0015982 Jy ! Component: 900 - total flux cleaned = -0.00204541 Jy ! Component: 950 - total flux cleaned = -0.00210857 Jy ! Component: 1000 - total flux cleaned = -0.00261138 Jy ! Component: 1050 - total flux cleaned = -0.00273479 Jy ! Component: 1100 - total flux cleaned = -0.00285817 Jy ! Component: 1150 - total flux cleaned = -0.00310411 Jy ! Component: 1200 - total flux cleaned = -0.00316526 Jy ! Component: 1250 - total flux cleaned = -0.00334695 Jy ! Component: 1300 - total flux cleaned = -0.00346689 Jy ! Component: 1350 - total flux cleaned = -0.00340803 Jy ! Component: 1400 - total flux cleaned = -0.00376356 Jy ! Component: 1450 - total flux cleaned = -0.00370526 Jy ! Component: 1500 - total flux cleaned = -0.00393836 Jy ! Component: 1550 - total flux cleaned = -0.00376444 Jy ! Component: 1600 - total flux cleaned = -0.0041106 Jy ! Component: 1650 - total flux cleaned = -0.00422552 Jy ! Component: 1700 - total flux cleaned = -0.00428253 Jy ! Component: 1750 - total flux cleaned = -0.00428233 Jy ! Component: 1800 - total flux cleaned = -0.0046208 Jy ! Component: 1850 - total flux cleaned = -0.00473246 Jy ! Component: 1900 - total flux cleaned = -0.0047326 Jy ! Component: 1950 - total flux cleaned = -0.00500995 Jy ! Component: 2000 - total flux cleaned = -0.00528576 Jy ! Component: 2050 - total flux cleaned = -0.00528556 Jy ! Component: 2100 - total flux cleaned = -0.00555895 Jy ! Component: 2150 - total flux cleaned = -0.00550481 Jy ! Component: 2200 - total flux cleaned = -0.0056674 Jy ! Component: 2250 - total flux cleaned = -0.00588271 Jy ! Component: 2300 - total flux cleaned = -0.0059895 Jy ! Component: 2350 - total flux cleaned = -0.00609602 Jy ! Component: 2400 - total flux cleaned = -0.0061491 Jy ! Component: 2450 - total flux cleaned = -0.00614897 Jy ! Component: 2500 - total flux cleaned = -0.00630709 Jy ! Component: 2550 - total flux cleaned = -0.00651662 Jy ! Component: 2600 - total flux cleaned = -0.00646454 Jy ! Component: 2650 - total flux cleaned = -0.00651647 Jy ! Component: 2700 - total flux cleaned = -0.00672381 Jy ! Component: 2750 - total flux cleaned = -0.00656947 Jy ! Component: 2800 - total flux cleaned = -0.00656939 Jy ! Component: 2850 - total flux cleaned = -0.00667136 Jy ! Component: 2900 - total flux cleaned = -0.00656933 Jy ! Component: 2950 - total flux cleaned = -0.00692431 Jy ! Component: 3000 - total flux cleaned = -0.00657139 Jy ! Component: 3050 - total flux cleaned = -0.0066718 Jy ! Component: 3100 - total flux cleaned = -0.00662202 Jy ! Component: 3150 - total flux cleaned = -0.00677197 Jy ! Component: 3200 - total flux cleaned = -0.00662311 Jy ! Component: 3250 - total flux cleaned = -0.00672244 Jy ! Component: 3300 - total flux cleaned = -0.0069692 Jy ! Component: 3350 - total flux cleaned = -0.00677241 Jy ! Component: 3400 - total flux cleaned = -0.00701732 Jy ! Component: 3450 - total flux cleaned = -0.00687121 Jy ! Component: 3500 - total flux cleaned = -0.00711421 Jy ! Component: 3550 - total flux cleaned = -0.00701763 Jy ! Component: 3600 - total flux cleaned = -0.00701754 Jy ! Component: 3650 - total flux cleaned = -0.00735412 Jy ! Component: 3700 - total flux cleaned = -0.00725843 Jy ! Component: 3750 - total flux cleaned = -0.00735382 Jy ! Component: 3800 - total flux cleaned = -0.00754383 Jy ! Component: 3850 - total flux cleaned = -0.0074016 Jy ! Component: 3900 - total flux cleaned = -0.00735459 Jy ! Component: 3950 - total flux cleaned = -0.00763695 Jy ! Component: 4000 - total flux cleaned = -0.00773073 Jy ! Component: 4050 - total flux cleaned = -0.00787096 Jy ! Total flux subtracted in 4096 components = -0.00763821 Jy ! Clean residual min=-0.002611 max=0.002782 Jy/beam ! Clean residual mean=-0.000030 rms=0.000858 Jy/beam ! Combined flux in latest and established models = 1.6313 Jy keep ! Adding 315 model components to the UV plane model. ! The established model now contains 613 components and 1.6313 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /scr/pima/rv119_uvs/J1800+7828_S_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.579 x 4.514 at -75.42 degrees (North through East) ! Clean map min=-0.005691 max=1.2582 Jy/beam ! Writing clean map to FITS file: /scr/pima/rv119_uvs/J1800+7828_S_map.fits wmodel /scr/pima/rv119_uvs/J1800+7828_S_map.mod ! Writing 613 model components to file: /scr/pima/rv119_uvs/J1800+7828_S_map.mod wobs /scr/pima/rv119_uvs/J1800+7828_S_uvs.fits ! Writing UV FITS file: /scr/pima/rv119_uvs/J1800+7828_S_uvs.fits wwins /scr/pima/rv119_uvs/J1800+7828_S_map.win ! wwins: Wrote 12 windows to /scr/pima/rv119_uvs/J1800+7828_S_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 1.25819 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000857921 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.204561 0.204561 0.409122 0.818245 1.63649 3.27298 6.54596 13.0919 26.1838 52.3677 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0025 print "peak= ", int(1000*peak_flux)/1000 ! peak= 1.258 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 1466.55 quit ! Quitting program ! Log file /scr/pima/rv119_uvs/J1800+7828_S_dfm.log closed on Sun Dec 17 11:38:24 2017