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])