#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
#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
# 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?