from astropy.table import Table, vstack

t=Table.read('rc3_Tlt0_CZ_lt10k.fits')
tout=Table()

for i in range(len(t)):
	ra=t['_RAJ2000'][i]
	dec=t['_DEJ2000'][i]
	altname=t['altname'][i]	
	print  i,altname.replace(' ',''), ra, dec
	url="http://hla.stsci.edu/cgi-bin/hlaSIAP.cgi?POS=%f,%f&SIZE=0" % (ra,dec)
	thla=Table.read(url)
	if len(thla)!=0:
		print 'Find ',len(thla),' rows in HLA...'
		col_altname=Table.Column(name='altname',data=[altname.replace(' ','')]*len(thla))
		col_pgc=Table.Column(name='PGC',data=[t['PGC'][i].replace(' ','')]*len(thla))
		col_ra=Table.Column(name='_RAJ2000',data=[ra]*len(thla))
		col_dec=Table.Column(name='_DECJ2000',data=[dec]*len(thla))
		thla.add_columns([col_altname,col_pgc,col_ra,col_dec])
		tout=vstack([tout,thla])

tout.write('HLA_rc3_Tlt0_CZ_lt10k.csv',format='csv')