diff --git a/Metallicity_Stack_Commons/analysis/composite_indv_detect.py b/Metallicity_Stack_Commons/analysis/composite_indv_detect.py
index 8098777..82c38d4 100644
--- a/Metallicity_Stack_Commons/analysis/composite_indv_detect.py
+++ b/Metallicity_Stack_Commons/analysis/composite_indv_detect.py
@@ -53,85 +53,3 @@ def main(fitspath, dataset, composite_file, outfile):
 
     # Write Astropy ASCII table containing composite T_e and derived metallicity
     det3_table.write(outfile, format='ascii.fixed_width_two_line')
-
-
-def run_ind_detection(fitspath, dataset, average_value_ascii):
-
-    N_gal_tab = asc.read(average_value_ascii)
-    ID = N_gal_tab['ID']
-    for aa in range(len(ID)):
-        ind_detection(fitspath, dataset, ID[aa])
-    new_name = fitspath + 'Individual_ratio_temperature.tbl'
-    vertical_stacking(fitspath, dataset, new_name)
-    print('run complete')
-
-
-def ind_detection(fitspath, dataset, bin_id):
-    get_det3_tab = asc.read(fitspath + 'get_det3_table2.tbl')
-    bin_tab = asc.read(fitspath + dataset + '_2d_binning_datadet3.tbl')
-    N_gal_tab = asc.read(fitspath + dataset + '_Average_R23_O32_Values.tbl')
-    stackmeas_tab = asc.read(fitspath + dataset + '_temperatures_metalicity.tbl')
-
-    # From tables
-    Source_id = get_det3_tab['Individual_IDs']
-    O4959 = get_det3_tab['O4959']
-    O5007 = get_det3_tab['O5007']
-    Bin_number = bin_tab['Bin_number']
-    O2 = get_det3_tab['O2']
-    O3 = get_det3_tab['O3']
-    Hb = get_det3_tab['Hb']
-    N_Galaxies = N_gal_tab['N_Galaxies']
-    temp_bin = stackmeas_tab['Temperature']
-
-    R23 = get_det3_tab['R23']
-    O32 = get_det3_tab['O32']
-
-    # Initializing Arrays
-
-    Source_IDs = []
-    Bin_ID = []
-    two_beta = []
-    three_beta = []
-    OIII4959 = []
-    OIII5007 = []
-    HBeta = []
-    average_temp = []
-    R23_ind = []
-    O32_ind = []
-
-    for ii in range(len(O2)):
-        if Bin_number[ii] == bin_id:
-            # print 'Bin_number:', Bin_number[ii], 'O2:', O2[ii], 'O3:', O3[ii], 'Hb:', Hb[ii]
-            Source_ID.append(Source_id[ii])
-            Bin_ID.append(bin_id)
-            R23_ind.append(R23[ii])
-            O32_ind.append(O32[ii])
-            two_beta.append(O2[ii] / Hb[ii])
-            three_beta.append(O3[ii] / Hb[ii])
-            OIII4959.append(O4959[ii])
-            OIII5007.append(O5007[ii])
-            HBeta.append(Hb[ii])
-            average_temp.append(temp_bin[bin_id])
-
-    individual_ascii = '/Users/reagenleimbach/Desktop/Zcalbase_gal/individual_detection/' + str(
-        bin_id) + '_individual_ratios_temp.tbl'
-    n = ('Source_ID', 'Bin_ID', 'Individual_R23', 'Individual_O32', 'two_beta', 'three_beta', 'OIII4959', 'OIII5007',
-         'HBeta', 'Temperature')  # 'ID', 'R23_Average', 'O32_Average'
-    ind_tab = Table(
-        [Source_ID, Bin_ID, R23_ind, O32_ind, two_beta, three_beta, OIII4959, OIII5007, HBeta, average_temp],
-        names=n)  # ID, R23, O32,
-    asc.write(ind_tab, individual_ascii, format='fixed_width_two_line')
-
-
-def individual_galaxy_table_stacking(fitspath, dataset, new_name):
-    individual_ascii = '/Users/reagenleimbach/Desktop/Zcalbase_gal/individual_detection/*_individual_ratios_temp.tbl'
-    table_files = glob.glob(individual_ascii)
-    table_files.sort()
-
-    for ii in range(len(table_files)):
-        asc_tab = asc.read(table_files[ii])
-        if ii == 0:
-            vstacking = asc_tab
-        else:
-            vstacking = vstack([vstacking, asc_tab])
-    asc.write(vstacking, new_name, format='fixed_width_two_line', overwrite=True)