예제 #1
0
#changes made to get to OrionS_rawACSmod
#modifications to label sig/ref positions

os.environ['CASAPATH'] = casapath

s = sd.scantable(
    'OrionS_rawACSmod',
    False)  #load the data without averaging		# filein,'Orion-S.raw.fits'

#s.summary()				#summary info					# summary
# fileout,'Orion-S-reduced.fits'
s.set_fluxunit('K')  # make 'K' default unit

#scal = sd.calps(s, [20,21,22,23])		# Calibrate CH3OH scans				# for i=21,24,2 do begin getps,i,ifnum=2,plnum=0,units='Ta*',
scal = sd.calps(
    s, [21, 22, 23, 24]
)  # Calibrate CH3OH scans				# for i=21,24,2 do begin getps,i,ifnum=2,plnum=0,units='Ta*',
del s  # remove s from memory
# recalculate az/el (NOT needed for GBT data)
antennaname = scal.get_antennaname()
if (antennaname != 'GBT'):
    scal.recalc_azel(
    )  # recalculate az/el to 		# tau=0.09 & accum & getps, i, ifnum=2,plnum=1,units='Ta*',
scal.opacity(0.09)  # do opacity correction				# tau=0.09 & accum & end & ave
sel = sd.selector()  # Prepare a selection
sel.set_ifs(2)  # select CH3OH IF
scal.set_selection(sel)  # get this IF
stave = sd.average_time(scal, weight='tintsys')  # average in time
spave = stave.average_pol(
    weight='tsys')  # average polarizations;Tsys-weighted (1/Tsys**2) average
sd.plotter.plot(spave)  # plot
예제 #2
0
					#scan numbers (zero-based) as compared to GBTIDL

					#changes made to get to OrionS_rawACSmod
					#modifications to label sig/ref positions

os.environ['CASAPATH'] = casapath


s = sd.scantable('OrionS_rawACSmod', False)#load the data without averaging		# filein,'Orion-S.raw.fits'

#s.summary()				#summary info					# summary
											# fileout,'Orion-S-reduced.fits'
s.set_fluxunit('K')         		# make 'K' default unit

#scal = sd.calps(s, [20,21,22,23])		# Calibrate CH3OH scans				# for i=21,24,2 do begin getps,i,ifnum=2,plnum=0,units='Ta*',
scal = sd.calps(s, [21,22,23,24])		# Calibrate CH3OH scans				# for i=21,24,2 do begin getps,i,ifnum=2,plnum=0,units='Ta*',
del s                                   # remove s from memory
# recalculate az/el (NOT needed for GBT data)
antennaname = scal.get_antennaname()
if ( antennaname != 'GBT'): scal.recalc_azel()      # recalculate az/el to 		# tau=0.09 & accum & getps, i, ifnum=2,plnum=1,units='Ta*',
scal.opacity(0.09)			# do opacity correction				# tau=0.09 & accum & end & ave
sel = sd.selector()			# Prepare a selection
sel.set_ifs(2)				# select CH3OH IF
scal.set_selection(sel)			# get this IF
stave = sd.average_time(scal, weight='tintsys')	# average in time
spave = stave.average_pol(weight='tsys')	# average polarizations;Tsys-weighted (1/Tsys**2) average
sd.plotter.plot(spave)			# plot

# what is going on with autoscaling?

					# do some smoothing
예제 #3
0
# changes made to get to OrionS_rawACSmod
# modifications to label sig/ref positions

os.environ["CASAPATH"] = casapath


s = sd.scantable("OrionS_rawACSmod", False)  # load the data without averaging		# filein,'Orion-S.raw.fits'

# s.summary()				#summary info					# summary
# fileout,'Orion-S-reduced.fits'
s.set_fluxunit("K")  # make 'K' default unit

# scal = sd.calps(s, [24,25,26,27])		# Calibrate SiO scans				# for i=25,28,2 do begin getps,i,ifnum=3,plnum=0,units='Ta*',
scal = sd.calps(
    s, [25, 26, 27, 28]
)  # Calibrate SiO scans				# for i=25,28,2 do begin getps,i,ifnum=3,plnum=0,units='Ta*',
del s  # remove s from memory
# recalculate az/el (NOT needed for GBT data)
antennaname = scal.get_antennaname()
if antennaname != "GBT":
    scal.recalc_azel()  # recalculate az/el to 		# tau=0.09 & accum & getps, i, ifnum=3,plnum=1,units='Ta*',
scal.opacity(0.09)  # do opacity correction				# tau=0.09 & accum & end & ave
sel = sd.selector()  # Prepare a selection
sel.set_ifs(15)  # select SiO IF
scal.set_selection(sel)  # get this IF
stave = sd.average_time(scal, weight="tintsys")  # average in time
spave = stave.average_pol(weight="tsys")  # average polarizations;Tsys-weighted (1/Tsys**2) average
sd.plotter.plot(spave)  # plot

# what is going on with autoscaling?