os.system("rm smooth_dir/smooth_*.fits") ### Go through the folders, identify and reduce each smooth image smooth_files_list = [] stars_to_try = ["ltt4364","eg131"] for star in stars_to_try: count = 1 for i in range(len(folders_to_search)): file_path_current = file_path + folders_to_search[i] ### Find any smooth star exposures ccdlist_info = iraf.ccdlist(file_path_current + "*.fits", Stdout = 1) ccdlist_info = functions.ccdlist_extract(ccdlist_info) ccdlist_match = functions.ccdlist_identify(ccdlist_info,star) ### Reduce any such exposures if len(ccdlist_match) > 0: for j in range(len(ccdlist_match)): file_name = ccdlist_info[ccdlist_match[j]][2] + ".fits" file_path_temp = file_path_current + "temp/" file_path_reduced = file_path_current + "reduced/" print "currently reducing image " + file_path_current + file_name smooth_files_list,count = reduce_image(file_path_current,file_name,file_path_reduced,smooth_files_list,count) ### If exposures were found for starX, then break the for loop if count > 1: break
######################## ### Start of program ### ######################## ### Set file_path file_path = sys.argv[1] #file_path = "/priv/miner3/hat-south/george/Honours/data/wifes/2010/" file_path_temp = file_path + "temp/" file_path_reduced = file_path + "reduced/" ############################################ ### Open file_path and find all SpecPhot ### ############################################ ### Get a list of fits files with objects begining in "HD" object_list = functions.ccdlist_extract(iraf.ccdlist(file_path +"*.fits",Stdout = 1)) HD_match = functions.ccdlist_identify_HD(object_list) HD_match = [] for i in range(len(object_list)): HD_match.append(i) ### Open those images and check if NOTES says "RV Standard" SP_match = check_head(HD_match,"NOTES","SpecPhot") ### Write the file_names to temp textfile SP_list = open(file_path_temp + "SpecPhot_list","w") for i in range(len(SP_match)): file_name = object_list[SP_match[i]][2] + "\n" SP_list.write(file_name)