temp_outfile = open(survey_name + 'mask.pickle', 'wb') pickle.dump((mask_resolution, maskfile), temp_outfile) temp_outfile.close() ################################################################################ # # FILTER GALAXIES # ################################################################################ temp_infile = open(survey_name + 'mask.pickle', 'rb') mask_resolution, maskfile = pickle.load(temp_infile) temp_infile.close() coord_min_table, mask, ngrid = filter_galaxies(infile, maskfile, mask_resolution, min_dist, max_dist, survey_name, mag_cut, rm_isolated) temp_outfile = open(survey_name + "filter_galaxies_output.pickle", 'wb') pickle.dump((coord_min_table, mask, ngrid), temp_outfile) temp_outfile.close() ################################################################################ # # FIND VOIDS # ################################################################################ temp_infile = open(survey_name + "filter_galaxies_output.pickle", 'rb') coord_min_table, mask, ngrid = pickle.load(temp_infile) temp_infile.close()
''' print(infile.columns) print('Please rename columns') exit() ''' infile['magnitude'].name = 'rabsmag' # Calculate comoving distance if 'Rgal' not in infile.columns: infile['Rgal'] = Distance(infile['z'], Omega_M, h) ################################################################################ # # FILTER GALAXIES # ################################################################################ coord_min_table, mask, ngrid = filter_galaxies(infile, maskfile, min_dist, max_dist, survey_name) ################################################################################ # # FIND VOIDS # ################################################################################ find_voids(ngrid, min_dist, max_dist, coord_min_table, mask, out1_filename, out2_filename, survey_name)