args = read_command_line(sys.argv)
print args
outfiles = []
try:
seqno, restFreqHz, outcube = make_cube(args) # make the image cube
outfiles.append(outcube)
except ValueError, msg:
print 'ERROR: ', msg
print 'Please run this command on: ',
print 'arcturus.gb.nrao.edu'
print 'If you are on arcturus, please report this error.'
sys.exit(-1)
except:
print(
'ERROR: There was a problem making the image cube. Please '
'see the log files for details')
sys.exit(-3)
# make the moment 0 map
if not args.noave:
outcont = make_average_map(restFreqHz, args.uniqueid)
outfiles.append(outcont)
if not args.noave and not args.noline:
outline = make_cube_minus_continuum(restFreqHz, seqno, args.uniqueid)
outfiles.append(outline)
cat.show()
fixAipsImages(outfiles)
imlin.box[1][4]=round(nChan*0.89)
imlin.order=0 # beware using imlin.order=1; use imlin.order=0
print imlin.box
imlin.go()
#Run trans task on sdgrd file
trans.inname=AIPSCat()[mydisk][-1].name
trans.inclass=AIPSCat()[mydisk][-1].klass
trans.inseq=AIPSCat()[mydisk][-1].seq
trans.transc= '231'
trans.outdi=mydisk
trans.outcl='baseli'
trans.go()
## and write the last thing now in the catalog to disk
fittp.indisk=mydisk
fittp.inname=AIPSCat()[mydisk][-1].name
fittp.inclass=AIPSCat()[mydisk][-1].klass
fittp.inseq=AIPSCat()[mydisk][-1].seq
outline = outName+'_line.fits'
if os.path.exists(outline):
os.remove(outline)
print 'Removed existing file to make room for new one :',outline
fittp.dataout='PWD:'+outline
fittp.go()
# clean up output images: reset FREQ axis to appropriate reference frame
fixAipsImages([outcube,outline,outcont])
if __name__ == "__main__":
args = read_command_line(sys.argv)
print args
outfiles = []
try:
seqno, restFreqHz, outcube = make_cube(args) # make the image cube
outfiles.append(outcube)
except ValueError, msg:
print "ERROR: ", msg
print "Please run this command on: ",
print "arcturus.gb.nrao.edu"
print "If you are on arcturus, please report this error."
sys.exit(-1)
except:
print ("ERROR: There was a problem making the image cube. Please " "see the log files for details")
sys.exit(-3)
# make the moment 0 map
if not args.noave:
outcont = make_average_map(restFreqHz, args.uniqueid)
outfiles.append(outcont)
if not args.noave and not args.noline:
outline = make_cube_minus_continuum(restFreqHz, seqno, args.uniqueid)
outfiles.append(outline)
cat.show()
fixAipsImages(outfiles)
imlin.box[1][4] = round(nChan * 0.89)
imlin.order = 0 # beware using imlin.order=1; use imlin.order=0
print imlin.box
imlin.go()
#Run trans task on sdgrd file
trans.inname = AIPSCat()[mydisk][-1].name
trans.inclass = AIPSCat()[mydisk][-1].klass
trans.inseq = AIPSCat()[mydisk][-1].seq
trans.transc = '231'
trans.outdi = mydisk
trans.outcl = 'baseli'
trans.go()
## and write the last thing now in the catalog to disk
fittp.indisk = mydisk
fittp.inname = AIPSCat()[mydisk][-1].name
fittp.inclass = AIPSCat()[mydisk][-1].klass
fittp.inseq = AIPSCat()[mydisk][-1].seq
outline = outName + '_line.fits'
if os.path.exists(outline):
os.remove(outline)
print 'Removed existing file to make room for new one :', outline
fittp.dataout = 'PWD:' + outline
fittp.go()
# clean up output images: reset FREQ axis to appropriate reference frame
fixAipsImages([outcube, outline, outcont])
