-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathipython-20141029-152045.log
85 lines (85 loc) · 2.62 KB
/
ipython-20141029-152045.log
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
#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-152045.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")
_ip.system("ls -F *.txt")
_ip.system("ls -F -d *.image")
execfile('imAnalysis.py')
ims
execfile('imAnalysis.py')
viewer()
_ip.system("ls -F *.py")
_ip.system("mkdir testresults")
_ip.system("ls -F -d *.image")
cp timesteps*tp?.image timesteps*tp?.residual timesteps*tp?.psf timesteps*tp?.model timesteps*tp?.model timesteps*tp?_subHD80606.* testresults/
cp timesteps*tp?.image timesteps*tp?.resid* testresults/
man cp
clean(vis='L85562_all_flagged.MS/
execfile('imcalcs_script.py')
execfile('imcalcs_script.py')
execfile('imcalcs_script.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()
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
dbpixsz = synthbeamasec/3
AprimeB = (FOVasec/2)**2*np.pi
Apix = (dbpixsz/2)**2*np.pi
maxbaselines = baselinesm1.max()*np.ones((np.shape(centerlambda)[0], 1))/centerlambda
maxbaselines = baselinesm1.max()*np.ones((np.shape(centerlambda)[0], 1))/centerlambda
maxbaselines
maxbaselines.shape()
shape(maxbaselines)
len(maxbaselines)
maxbaselines
maxbaselines[1]
maxbaselines[15]
maxbaselines[0][15]
nwplanes = (0.31/dA)*baselinesm1.max()*centerlambda/diamsavg
nwplanes
FOVamin = FOVasec/60
FOVamin
FOVdeg
execfile('timeClean_testing.py')
execfile('timeClean_testing.py')
execfile('timeClean_testing.py')