ipython-20141029-191114.log

#log# Automatic Logger file. *** THIS MUST BE THE FIRST LINE ***
#log# DO NOT CHANGE THIS LINE OR THE TWO BELOW -* coding: UTF-8 *-
#log# opts = Struct({'__allownew': True,
 'autocall': 2,
 'banner': 0,
 'colors': 'LightBG',
 'ipythondir': '/home/mknapp/.casa/ipython',
 'logfile': 'ipython-20141029-191114.log',
 'messages': 0,
 'prompt_in1': 'CASA <\\#>: '})
#log# args = []
#log# It is safe to make manual edits below here.
#log#-----------------------------------------------------------------------
execfile("/usr/local/bin/lib64/python2.7/assignmentFilter.py")
import numpy as np
import scipy as sp
import scipy.constants as const
import glob
import os
import matplotlib.pyplot as plt
import time
scriptmode=True
vis = 'L85562_all_flagged.MS'
r2asec = 206264.806
# Get list of subband center frequencies
tb.open(vis+'/SPECTRAL_WINDOW')
centerfqs = tb.getcol('REF_FREQUENCY')
tb.close()
# Convert to wavelength
centerlambda = const.c/centerfqs
# Get station diameter
tb.open(vis+'/ANTENNA')
diams = tb.getcol('DISH_DIAMETER')
tb.close()
diamsavg = np.mean(diams)
baselinesm = au.getBaselineLengths(vis, sort=True)
baselinesm1 = np.array(tuple(x[1] for x in baselinesm))
# Convert to wavelength
maxbaselines = baselinesm1.max()*np.ones((np.shape(centerlambda)[0], 1))/centerlambda
# Number of w-planes
dA = 0.5
nwplanes = (0.31/dA)*baselinesm1.max()*centerlambda/diamsavg
# Expected beamsize
FOV = centerlambda/diamsavg
FOVasec = FOV*r2asec
FOVamin = FOVasec/60
FOVdeg = FOVamin/60
synthbeam = centerlambda/baselinesm1.max()
synthbeamasec = synthbeam*r2asec
# How many pixels across primary beam?
dbpixsz = synthbeamasec/3
AprimeB = (FOVasec/2)**2*np.pi
Apix = (dbpixsz/2)**2*np.pi
execfile('timeClean_testing.py')
execfile('timeClean_testing.py')
_ip.system("ls -F ")
casaviewer
viewer90
viewer()
execfile('timeClean_testing.py')
execfile('timeClean_testing.py')