! Started logfile: /g0/scr/pima/npcsc_uvs/J1000+8127_S_uvs.log on Sun Oct  6 14:27:57 2013
obs /g0/scr/pima/npcsc_uvs/J1000+8127_S_uva.fits
! Reading UV FITS file: /g0/scr/pima/npcsc_uvs/J1000+8127_S_uva.fits
! AN table 1: 45 integrations on 15 of 15 possible baselines.
! Apparent sampling: 0.8 visibilities/baseline/integration-bin.
! Found source: J1000+8127
! 
! 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 2160 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.5 
print "taper_size =",taper_size
! taper_size = 20 
![@/home/lpetrov/bin/01_muppet /g0/scr/pima/npcsc_uvs/J1000+8127_S]
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.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=1.021680Jy  sigma=1.833647
! Fit after  self-cal, rms=0.821838Jy  sigma=1.275354
! 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=10.1 mas, bmaj=25.65 mas, bpa=45.34 degrees
! Estimated noise=26.9791 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=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.162646 Jy
! Component: 100  -  total flux cleaned = 0.227049 Jy
! Total flux subtracted in 100 components = 0.227049 Jy
! Clean residual min=-0.027286 max=0.030870 Jy/beam
! Clean residual mean=-0.000356 rms=0.011486 Jy/beam
! Combined flux in latest and established models = 0.227049 Jy
! Performing phase self-cal
! Adding 57 model components to the UV plane model.
! The established model now contains 57 components and 0.227049 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.102401Jy  sigma=0.231800
! Fit after  self-cal, rms=0.102448Jy  sigma=0.231547
! 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=12.99 mas, bmaj=27.99 mas, bpa=42.4 degrees
! Estimated noise=9.77364 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.0446088 Jy
! Component: 100  -  total flux cleaned = 0.0551009 Jy
! Component: 150  -  total flux cleaned = 0.0524247 Jy
! Component: 200  -  total flux cleaned = 0.0511389 Jy
! Total flux subtracted in 200 components = 0.0511389 Jy
! Clean residual min=-0.009011 max=0.009039 Jy/beam
! Clean residual mean=-0.000174 rms=0.004988 Jy/beam
! Combined flux in latest and established models = 0.278188 Jy
! Performing phase self-cal
! Adding 128 model components to the UV plane model.
! The established model now contains 181 components and 0.278188 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.085243Jy  sigma=0.152120
! Fit after  self-cal, rms=0.085111Jy  sigma=0.151501
! Inverting map 
print "************** FINISHED NATURAL WEIGHTING CLEAN **************"
! ************** FINISHED NATURAL WEIGHTING CLEAN ************** 
! 
gscale true
! Performing overall amplitude self-cal
! 
! Correcting IF 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.76     FD        0.96     KP        1.00     LA        1.06     
!   OV        0.93     PT        0.97     
! 
! 
! Correcting IF 2.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.82     FD        0.98     KP        1.01     LA        1.01     
!   OV        0.96     PT        0.90     
! 
! 
! Correcting IF 3.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.89     FD        1.08     KP        1.02     LA        0.88     
!   OV        1.09     PT        1.01     
! 
! 
! Correcting IF 4.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.80     FD        0.99     KP        0.96     LA        0.77     
!   OV        0.83     PT        0.89     
! 
! 
! Fit before self-cal, rms=0.085111Jy  sigma=0.151501
! Fit after  self-cal, rms=0.077326Jy  sigma=0.127088
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=12.76 mas, bmaj=27.45 mas, bpa=43.67 degrees
! Estimated noise=8.89396 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.00665383 Jy
! Component: 100  -  total flux cleaned = -0.00820832 Jy
! Component: 150  -  total flux cleaned = -0.00857709 Jy
! Component: 200  -  total flux cleaned = -0.0093818 Jy
! Total flux subtracted in 200 components = -0.0093818 Jy
! Clean residual min=-0.001523 max=0.001518 Jy/beam
! Clean residual mean=-0.000002 rms=0.000599 Jy/beam
! Combined flux in latest and established models = 0.268807 Jy
! Performing phase self-cal
! Adding 57 model components to the UV plane model.
! The established model now contains 238 components and 0.268806 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.073168Jy  sigma=0.128318
! Fit after  self-cal, rms=0.073219Jy  sigma=0.128206
! 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.073219Jy  sigma=0.128206
! Fit after  self-cal, rms=0.071744Jy  sigma=0.126353
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=12.4 mas, bmaj=26.93 mas, bpa=44.99 degrees
! Estimated noise=8.3504 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.00929618 Jy
! Component: 100  -  total flux cleaned = -0.0109271 Jy
! Component: 150  -  total flux cleaned = -0.0112229 Jy
! Component: 200  -  total flux cleaned = -0.0109524 Jy
! Total flux subtracted in 200 components = -0.0109524 Jy
! Clean residual min=-0.001054 max=0.001058 Jy/beam
! Clean residual mean=-0.000002 rms=0.000430 Jy/beam
! Combined flux in latest and established models = 0.257854 Jy
selfcal
! Performing phase self-cal
! Adding 47 model components to the UV plane model.
! The established model now contains 281 components and 0.257854 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.067972Jy  sigma=0.128228
! Fit after  self-cal, rms=0.067970Jy  sigma=0.128217
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.54348e-05 Jy
! Component: 100  -  total flux cleaned = 3.20899e-05 Jy
! Component: 150  -  total flux cleaned = 0.000267169 Jy
! Component: 200  -  total flux cleaned = 0.000834069 Jy
! Total flux subtracted in 200 components = 0.000834069 Jy
! Clean residual min=-0.000707 max=0.000709 Jy/beam
! Clean residual mean=-0.000002 rms=0.000288 Jy/beam
! Combined flux in latest and established models = 0.258688 Jy
! Performing phase self-cal
! Adding 68 model components to the UV plane model.
! The established model now contains 340 components and 0.258688 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.067774Jy  sigma=0.128234
! Fit after  self-cal, rms=0.067778Jy  sigma=0.128231
! 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=0.067778Jy  sigma=0.128231
! Fit after  self-cal, rms=0.066305Jy  sigma=0.126069
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=12.03 mas, bmaj=26.44 mas, bpa=46.2 degrees
! Estimated noise=7.79386 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.00845876 Jy
! Component: 100  -  total flux cleaned = -0.00975301 Jy
! Component: 150  -  total flux cleaned = -0.0105491 Jy
! Component: 200  -  total flux cleaned = -0.010786 Jy
! Total flux subtracted in 200 components = -0.010786 Jy
! Clean residual min=-0.000843 max=0.000839 Jy/beam
! Clean residual mean=0.000002 rms=0.000325 Jy/beam
! Combined flux in latest and established models = 0.247902 Jy
selfcal
! Performing phase self-cal
! Adding 62 model components to the UV plane model.
! The established model now contains 398 components and 0.247902 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.062945Jy  sigma=0.127997
! Fit after  self-cal, rms=0.062941Jy  sigma=0.127990
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.00015748 Jy
! Component: 100  -  total flux cleaned = 0.000490486 Jy
! Component: 150  -  total flux cleaned = 0.00101857 Jy
! Component: 200  -  total flux cleaned = 0.00166124 Jy
! Total flux subtracted in 200 components = 0.00166124 Jy
! Clean residual min=-0.000576 max=0.000582 Jy/beam
! Clean residual mean=-0.000000 rms=0.000225 Jy/beam
! Combined flux in latest and established models = 0.249563 Jy
! Performing phase self-cal
! Adding 71 model components to the UV plane model.
! The established model now contains 466 components and 0.249563 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.062774Jy  sigma=0.128083
! Fit after  self-cal, rms=0.062778Jy  sigma=0.128080
! 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=0.062778Jy  sigma=0.128080
! Fit after  self-cal, rms=0.061482Jy  sigma=0.125920
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=13.12 mas, bmaj=27.56 mas, bpa=43.2 degrees
! Estimated noise=7.06583 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.00786915 Jy
! Component: 100  -  total flux cleaned = -0.00923599 Jy
! Component: 150  -  total flux cleaned = -0.0101091 Jy
! Component: 200  -  total flux cleaned = -0.010271 Jy
! Total flux subtracted in 200 components = -0.010271 Jy
! Clean residual min=-0.000767 max=0.000766 Jy/beam
! Clean residual mean=0.000002 rms=0.000303 Jy/beam
! Combined flux in latest and established models = 0.239292 Jy
selfcal
! Performing phase self-cal
! Adding 62 model components to the UV plane model.
! The established model now contains 514 components and 0.239292 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.058450Jy  sigma=0.127843
! Fit after  self-cal, rms=0.058447Jy  sigma=0.127836
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 = -4.48077e-05 Jy
! Component: 100  -  total flux cleaned = 0.000228857 Jy
! Component: 150  -  total flux cleaned = 0.000703907 Jy
! Component: 200  -  total flux cleaned = 0.00124443 Jy
! Total flux subtracted in 200 components = 0.00124443 Jy
! Clean residual min=-0.000482 max=0.000492 Jy/beam
! Clean residual mean=0.000000 rms=0.000198 Jy/beam
! Combined flux in latest and established models = 0.240537 Jy
! Performing phase self-cal
! Adding 77 model components to the UV plane model.
! The established model now contains 583 components and 0.240537 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.058291Jy  sigma=0.127968
! Fit after  self-cal, rms=0.058295Jy  sigma=0.127964
! 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=11.65 mas, bmaj=25.97 mas, bpa=47.46 degrees
! Estimated noise=7.25066 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 11.65 x 25.97 at 47.46 degrees (North through East)
! Clean map  min=-0.013806  max=0.16168 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
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=9.216 mas, bmaj=23.93 mas, bpa=50.02 degrees
! Estimated noise=18.9418 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.11886 Jy
! Component: 100  -  total flux cleaned = 0.172177 Jy
! Component: 150  -  total flux cleaned = 0.199506 Jy
! Component: 200  -  total flux cleaned = 0.212561 Jy
! Total flux subtracted in 200 components = 0.212561 Jy
! Clean residual min=-0.012018 max=0.012165 Jy/beam
! Clean residual mean=-0.000195 rms=0.005299 Jy/beam
! Combined flux in latest and established models = 0.212561 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 113 model components to the UV plane model.
! The established model now contains 113 components and 0.212561 Jy
! Inverting map and beam 
! Estimated beam: bmin=11.65 mas, bmaj=25.97 mas, bpa=47.46 degrees
! Estimated noise=7.25066 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.0129915 Jy
! Component: 100  -  total flux cleaned = 0.0144762 Jy
! Component: 150  -  total flux cleaned = 0.0124965 Jy
! Component: 200  -  total flux cleaned = 0.00877563 Jy
! Total flux subtracted in 200 components = 0.00877563 Jy
! Clean residual min=-0.005435 max=0.005450 Jy/beam
! Clean residual mean=-0.000151 rms=0.002932 Jy/beam
! Combined flux in latest and established models = 0.221336 Jy
! Adding 144 model components to the UV plane model.
! The established model now contains 253 components and 0.221336 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=13.12 mas, bmaj=27.56 mas, bpa=43.2 degrees
! Estimated noise=7.06583 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0044141 Jy
! Component: 100  -  total flux cleaned = 0.00461281 Jy
! Component: 150  -  total flux cleaned = 0.00357362 Jy
! Component: 200  -  total flux cleaned = 0.00342635 Jy
! Total flux subtracted in 200 components = 0.00342635 Jy
! Clean residual min=-0.002065 max=0.002060 Jy/beam
! Clean residual mean=-0.000048 rms=0.001137 Jy/beam
! Combined flux in latest and established models = 0.224763 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.000708842 Jy
! Component: 100  -  total flux cleaned = 0.00126732 Jy
! Component: 150  -  total flux cleaned = 0.00126588 Jy
! Component: 200  -  total flux cleaned = 0.00137699 Jy
! Total flux subtracted in 200 components = 0.00137699 Jy
! Clean residual min=-0.001834 max=0.001846 Jy/beam
! Clean residual mean=-0.000010 rms=0.000698 Jy/beam
! Combined flux in latest and established models = 0.22614 Jy
! Adding 226 model components to the UV plane model.
! The established model now contains 476 components and 0.22614 Jy
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
clean
! Inverting map and beam 
! Estimated beam: bmin=11.65 mas, bmaj=25.97 mas, bpa=47.46 degrees
! Estimated noise=7.25066 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00130103 Jy
! Component: 100  -  total flux cleaned = 0.00179537 Jy
! Component: 150  -  total flux cleaned = 0.00212186 Jy
! Component: 200  -  total flux cleaned = 0.00255255 Jy
! Total flux subtracted in 200 components = 0.00255255 Jy
! Clean residual min=-0.000684 max=0.000688 Jy/beam
! Clean residual mean=-0.000010 rms=0.000274 Jy/beam
! Combined flux in latest and established models = 0.228692 Jy
selfcal
! Performing phase self-cal
! Adding 82 model components to the UV plane model.
! The established model now contains 556 components and 0.228692 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.058223Jy  sigma=0.128490
! Fit after  self-cal, rms=0.058195Jy  sigma=0.128418
wmodel /g0/scr/pima/npcsc_uvs/J1000+8127_S_map.mod
! Writing 556 model components to file: /g0/scr/pima/npcsc_uvs/J1000+8127_S_map.mod
wobs /g0/scr/pima/npcsc_uvs/J1000+8127_S_uvs.fits
! Writing UV FITS file: /g0/scr/pima/npcsc_uvs/J1000+8127_S_uvs.fits
wwins /g0/scr/pima/npcsc_uvs/J1000+8127_S_map.win
! wwins: No CLEAN windows to write.
! Error occured in command: wwins
![Exited script file: /home/lpetrov/bin/01_muppet]
! 
wmap "/g0/scr/pima/npcsc_uvs/J1000+8127_S\_map.fits"
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 11.65 x 25.97 at 47.46 degrees (North through East)
! Clean map  min=-0.0092083  max=0.1496 Jy/beam
! Writing clean map to FITS file: /g0/scr/pima/npcsc_uvs/J1000+8127_S_map.fits
float peak_flux
peak_flux = peak(flux,max)
print "The peak flux is",peak_flux
! The peak flux is 0.149597 
invert
! Inverting map 
float image_rms
image_rms = imstat(rms)
print "The final image rms is",image_rms
! The final image rms is 0.000260827 
loglevs 100*3*image_rms/peak_flux
! The new contour levels are:
!  -0.52306 0.52306 1.04612 2.09224 4.18448 8.36896 16.7379 33.4758 66.9516
print "clev=",int(3*image_rms*10000)/10000
! clev= 0.0007 
print "peak=",int(1000*peak_flux)/1000
! peak= 0.149 
quit
! Quitting program
! Log file /g0/scr/pima/npcsc_uvs/J1000+8127_S_uvs.log closed on Sun Oct  6 14:28:00 2013