CANDELS_Combining_Fields_SFR.py

from astropy.table import vstack, Table
from astropy.io import ascii

fields = ['COSMOS', 'EGS', 'GOODSN', 'GOODSS', 'UDS']
SFR_Combined_Fields = Table()

for each_field in fields:
    catalog = ascii.read('/Users/laurenhiggins/Dropbox/GEG_Research/Key_Parameters_Catalogs_Per_Field/Star_Formation_Rate_Catalogs_old/CANDELS.%s.F160W.v1.IRUV_SFRs.cat'%each_field)
    SFR_subset = catalog[['id','flag_prox', 'flag_lum', 'id_mips', 'sfr_best', 'sfr_flag']]
    SFR_subset['Field_flag'] = each_field
    SFR_Combined_Fields = vstack([SFR_Combined_Fields, SFR_subset])


ascii.write(SFR_Combined_Fields, '/Users/laurenhiggins/Dropbox/GEG_Research/Key_Parameters_Catalogs_Per_Field/Star_Formation_Rate_Catalogs_old/CANDELS.ALL.FIELDS.F160W.v1.IRUV_SFRs.cat', 
            names=('id', 'flag_prox', 'flag_lum', 'id_mips', 'sfr_best', 'sfr_flag', 'Field_Flag'), format='csv', overwrite=True)

'''
2. id
3. flag_prox
4. flag_lum
6. id_mips
8. sfr_best
7. sfr_flag
'''

'''
for each_field in fields:
    catalog = ascii.read('/Users/laurenhiggins/Dropbox/GEG_Research/Key_Parameters_Catalogs_Per_Field/Star_Formation_Rate_Catalogs_old/CANDELS.%s.F160W.v1.IRUV_SFRs.cat'%each_field)
    SFR_subset = Table(catalog[['id','flag_prox', 'flag_lum', 'id_mips', 'sfr_best', 'sfr_flag']], names=('id','flag_prox', 'flag_lum', 'id_mips', 'sfr_best', 'sfr_flag'))
    SFR_subset['Field_flag'] = each_field
    SFR_Combined_Fields = vstack([SFR_subset])
    print(SFR_Combined_Fields.dtype)
print(SFR_Combined_Fields)
'''