def coord(refmap,params): #takes input reference map (.sdf) and parameter file as produced by FellWalker) (.fit) data = pyfits.getdata(params) array = [] i = 0 for i in range(len(data)): ra = RaDecDegs.timeRa(data[i][1],'TRUE') dec = RaDecDegs.timeDec(data[i][2],'TRUE') coords = "%10.7f %s"%(ra, dec) cmd = stardir+"kappa/wcstran %s posin=\"\'%s\'\" framein=SKY frameout=PIXEL > /dev/null"%(refmap,coords) os.system(cmd) posout = (os.popen(stardir+"kappa/parget posout wcstran").readlines()) #print posout[0] (pixx,pixy,junk) = re.split("\s+",posout[0],2) #print pixx, pixy #print 'Ra '+str(ra)+' and Dec '+ str(dec) +' are coordinates '+ str(pixx)+' and ' + str(pixy) array.append([pixx,pixy]) i = i + 1 return array
err_Ji = Ji * (0.05) #frac. error in temp. if float(Ji) == 0: equ = 0 else: equ = float(m850i) / float(Ji) #Error in mass ratio: frac = ((dMi_850/float(m850i))**2.)+((err_Ji/float(Ji))**2.) err_R = equ * ((frac)**0.5) #Convert Ra/Dec from degrees to Time format. Ra = float(data[i-1][1]) Dec = float(data[i-1][2]) Ra_T = RaDecDegs.timeRa(Ra,'FALSE') Dec_T = RaDecDegs.timeDec(Dec,'FALSE') #Write to files #open files for appending SMM = open("SMM.txt","a") SMM_LTX = open("SMM_LTX.txt","a") #RAW data file SMM.write(str(i)+'\t') SMM.write(str(round(float(data[i-1][1]),2))+'\t') SMM.write(str(round(float(data[i-1][2]),4))+'\t') SMM.write(str(c450i)+'\t') SMM.write(str(c850i)+'\t') SMM.write(str(m450i)+'\t') SMM.write(str(m850i)+'\t') SMM.write(str(dMi_850)+'\t')
# Extract core mass as number cmd = "%s/stats ndf=%s %s" % (kapdir, m450, "> /dev/null") os.system(cmd) cmd = "%s/parget parname=%s applic=%s" % (kapdir, "total", "stats") status, output = commands.getstatusoutput(cmd) m450i = output cmd = "%s/stats ndf=%s %s" % (kapdir, m850, "> /dev/null") os.system(cmd) cmd = "%s/parget parname=%s applic=%s" % (kapdir, "total", "stats") status, output = commands.getstatusoutput(cmd) m850i = output print "mass of clump 850 =", m850i # Convert Ra/Dec from degrees to Time format. Ra = float(data[i - 1][1]) Dec = float(data[i - 1][2]) Ra_T = RaDecDegs.timeRa(Ra, "FALSE") Dec_T = RaDecDegs.timeDec(Dec, "FALSE") # Write to files #LaTex data file core_LTX.write(str(i) + "\t&\t") # core_LTX.write(Ra_T+'\t&\t') # core_LTX.write(Dec_T+'\t&\t') core_LTX.write(str(round(float(c450i), 1)) + "\t&\t") core_LTX.write(str(round(float(c850i), 2)) + "\t&\t") # core_LTX.write(str(round(float(m450i),1))+'('+str(round(float(dMi_450),1))+')\t&\t') core_LTX.write(str(round(float(m850i), 2)) + "(" + str(round(float(dMi_850), 2)) + ")\t&\t") core_LTX.write(str(round(float(temp_mi), 1)) + "(" + str(round(float(0.05 * temp_mi), 1)) + ") \n")